737155V1.Full.Pdf
Total Page:16
File Type:pdf, Size:1020Kb
bioRxiv preprint doi: https://doi.org/10.1101/737155; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. Ancestral reconstruction of sunflower karyotypes reveals dramatic chromosomal evolution Kate L. Ostevik1,2, Kieran Samuk1, and Loren H. Rieseberg2 1. Department of Biology, Duke University, Durham, NC, 27701 2. Department of Botany, University of British Columbia, Vancouver, BC, Canada, V6T 1Z4 Correspondence to: [email protected] 1 bioRxiv preprint doi: https://doi.org/10.1101/737155; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. 1 Abstract 2 3 Mapping the chromosomal rearrangements between species can inform our understanding of genome 4 evolution, reproductive isolation, and speciation. Here we present a systematic survey of chromosomal 5 rearrangements in the annual sunflowers, which is a group known for extreme karyotypic diversity. We 6 build high-density genetic maps for two subspecies of the prairie sunflower, Helianthus petiolaris ssp. 7 petiolaris and H. petiolaris ssp. fallax. Using a novel algorithm implemented in the accompanying R 8 package syntR, we identify blocks of synteny between these two subspecies and previously published 9 high-density genetic maps. We reconstruct ancestral karyotypes for annual sunflowers using those 10 synteny blocks and conservatively estimate that there have been 9.7 chromosomal rearrangements 11 per million years – a high rate of chromosomal evolution. Although the rate of inversion is even higher 12 than the rate of translocation in this group, we further find that every extant karyotype is distinguished 13 by between 1 and 3 translocations involving only 8 of the 17 chromosomes. This non-random sampling 14 suggests that certain chromosomes are prone to translocation and may thus contribute 15 disproportionately to widespread hybrid sterility in sunflowers. These data deepen our understanding 16 of chromosome evolution and confirm that Helianthus has an exceptional rate of chromosomal 17 rearrangement that likely facilitates similarly rapid diversification. 18 Introduction 19 20 Organisms vary widely in the number and arrangement of their chromosomes – i.e. their karyotype. 21 Interestingly, karyotypic differences are often associated with species boundaries and, therefore, 22 suggest a link between chromosomal evolution and speciation (White 1978, King 1993). Indeed, it is 23 well established that chromosomal rearrangements can contribute to reproductive isolation. Although 24 not a universal phenomenon, individuals heterozygous for divergent karyotypes are often sterile or 25 inviable, illustrating a direct link to reproductive isolation (King 1987, Lai et al. 2005, Stathos and 26 Fishman 2014). Apart from being underdominant per se, chromosomal rearrangements can also 27 facilitate the evolution of reproductive barriers by extending genomic regions that are protected from 2 bioRxiv preprint doi: https://doi.org/10.1101/737155; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. 28 introgression (Noor et al. 2001, Rieseberg 2001), accumulating genetic incompatibilities (Navarro and 29 Barton 2003), and simplifying reinforcement (Trickett and Butlin 1994). 30 31 Despite their prevalence, the fixation of underdominant chromosomal rearrangements is still not well 32 understood. The simplest model is that in sufficiently small populations, underdominant 33 rearrangements can be stochastically fixed by genetic drift (Lewis 1966, Templeton 1981, Grant 1981). 34 However, several other mechanisms have been proposed to explain chromosomal evolution (Rieseberg 35 2001). These include individually neutral chromosome fissions and fusions (Baker and Bickham 1986, 36 Britton-Davidian 2000), meiotic drive (Novitski 1951, White 1978, Chmátal et al. 2014) and association 37 with locally adapted genes (Kirkpatrick and Barton 2006, Lowry and Willis 2010). The importance of 38 each of these mechanisms will largely depend on the types of rearrangements fixed (e.g., chromosome 39 fusions, pericentric inversions, or reciprocal translocations) and the cost/benefit of the genes 40 associated with them. Consequently, mapping and characterizing chromosomal rearrangements is a 41 critical step towards understanding their evolutionary dynamics. 42 43 The genus Helianthus is well known to have particularly labile genome structure. These sunflowers 44 have several paleopolyploidy events in their evolutionary history (Barker et al. 2008), have given rise to 45 three homoploid hybrid species (Rieseberg 1991), and are prone to transposable element activity 46 (Kawakami et al. 2011, Staton et al. 2012). Evidence in the form of hybrid pollen inviability, abnormal 47 chromosome pairings during meiosis, and genetic map comparisons suggests that Helianthus 48 karyotypes are extremely diverse (Heiser 1947, Heiser 1951, Heiser 1961, Whelan 1979, Chandler 49 1986, Rieseberg et al. 1995, Quillet et al. 1995, Burke et al. 2004, Heesacker et al. 2009, Barb et al. 50 2014). In fact, annual sunflowers have one of the highest described rates of chromosomal evolution 51 across all plants and animals (Burke et al. 2004). Furthermore, while Helianthus species generally have 52 large effective population sizes (but see Stebbins and Daly 1961, Carney et al. 2000), where drift should 53 be relatively weak, chromosomal rearrangements appear to be strongly underdominant in this group 54 (Chandler 1986, Lai et al. 2005). This means that chromosomal evolution is less likely to be driven by 55 drift alone and more likely to be an important driver of reproductive isolation in sunflower. 56 3 bioRxiv preprint doi: https://doi.org/10.1101/737155; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. 57 Studying chromosomal evolution within any group requires high-density genetic maps and a method 58 for systematic comparison of karyotypes. Recently, high-density genetic maps of H. niveus ssp. 59 tephrodes and H. argophyllus were built and compared to H. annuus (Barb et al. 2014). This analysis 60 clearly mapped previously inferred karyotypes (Heiser 1951, Chandler 1986, Quillet et al. 1995) and 61 showed that chromosomal rearrangements affect introgression between sunflowers. That said, 62 creating additional high-density maps for related species would yield a fuller picture of chromosomal 63 evolution. 64 65 Previously constructed genetic maps with limited marker density suggest that several chromosomal 66 rearrangements differentiate H. annuus and H. petiolaris (Rieseberg et al. 1995, Burke et al. 2004) and 67 evidence from cytological surveys suggests that subspecies within H. petiolaris subspecies carry 68 divergent karyotypes (Heiser 1961). Creating high-density genetic maps for H. petiolaris and adding 69 them to the above analysis will allow us to (1) precisely track additional rearrangements and (2) 70 reconstruct ancestral karyotypes for the group, untangling overlapping rearrangements that can be 71 obscured by directly comparing present-day karyotypes and that are likely to occur in lineages with 72 extensive chromosomal rearrangement like sunflower. At the same time, these H. petiolaris maps will 73 be useful to breeders looking to introgress beneficial alleles from crop wild relatives into domesticated 74 sunflower lineages (Chetelat and Meglic 2000, Foulongne et al. 2003, Dirlewanger et al. 2004). 75 76 A key part of a multi-species comparative study of chromosome evolution is a systematic and 77 repeatable method for identifying syntenic chromosomal regions (sensu Pevzner and Tesler 2003) 78 using genetic map data. This is especially important for cases with high marker density, because 79 breakpoints between synteny blocks can be blurred by mapping error, micro-rearrangements, and 80 paralogy (Hackett and Broadfoot 2003, Choi et al. 2007, Barb et al. 2014, Bilton et al. 2018). In previous 81 studies, synteny blocks have been found by counting all differences in marker order (Wu and Tanksley 82 2010), by visual inspection (Burke et al. 2004), or by manually applying simple rules like size thresholds 83 (Heesacker et al. 2009, Barb et al. 2014) and Spearman’s rank comparisons (Berdan et al 2014, 84 Schlautman et al 2017). However, these methods become intractable and prone to error when applied 85 to very dense genetic maps. Furthermore, to our knowledge, there is no software available that 4 bioRxiv preprint doi: https://doi.org/10.1101/737155; this version posted August 15, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. 86 identifies