Interspecific and interploidal gene flow in Central European Arabidopsis (Brassicaceae)

Unambiguous chromosome identification reveals the factors impacting irregular chromosome behaviors in allotriploid AAC Brassica

KEY MESSAGE

The major irregular chromosome pairing and mis-segregation were detected during meiosis through unambiguous chromosome identification and found that allotriploid Brassica can undergo meiosis successfully and produce mostly viable aneuploid gametes. Triploids have played a crucial role in the evolution of species by forming polyploids and facilitating interploidy gene transfer. It is widely accepted that triploids cannot undergo meiosis normally and predominantly produce nonfunctional aneuploid gametes, which restricts their role in species evolution. In this study, we demonstrated that natural and synthetic allotriploid Brassica (AAC), produced by crossing natural and synthetic Brassica napus (AACC) with Brassica rapa (AA), exhibits basically normal chromosome pairing and segregation during meiosis. Homologous A chromosomes paired faithfully and generally segregated equally. Monosomic C chromosomes were largely retained as univalents and randomly entered daughter cells. The primary irregular meiotic behaviors included associations of homoeologs and 45S rDNA loci at diakinesis, as well as homoeologous chromosome replacement and premature sister chromatid separation at anaphase I. Preexisting homoeologous arrangements altered meiotic behaviors in both chromosome irregular pairing and mis-segregation by increasing the formation of A-genomic univalents and A-C bivalents, as well as premature sister chromatid separation and homologous chromosome nondisjunction. Meiotic behaviors depended significantly on the genetic background and heterozygous homoeologous rearrangement. AAC triploids mainly generated aneuploid gametes, most of which were viable. These results demonstrate that allotriploid Brassica containing an intact karyotype can proceed through meiosis successfully, broadening our current understanding of the inheritance and role in species evolution of allotriploid.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Interspecific transfer of genetic information through polyploid bridges

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Genetic and environmental manipulation of Arabidopsis hybridization barriers uncovers antagonistic functions in endosperm cellularization

Speciation involves reproductive isolation, which can occur by hybridization barriers acting in the endosperm of the developing seed. The nuclear endosperm is a nutrient sink, accumulating sugars from surrounding tissues, and undergoes coordinated cellularization, switching to serve as a nutrient source for the developing embryo. Tight regulation of cellularization is therefore vital for seed and embryonic development. Here we show that hybrid seeds from crosses between Arabidopsis thaliana as maternal contributor and A. arenosa or A. lyrata as pollen donors result in an endosperm based post-zygotic hybridization barrier that gives rise to a reduced seed germination rate. Hybrid seeds display opposite endosperm cellularization phenotypes, with late cellularization in crosses with A. arenosa and early cellularization in crosses with A. lyrata. Stage specific endosperm reporters display temporally ectopic expression in developing hybrid endosperm, in accordance with the early and late cellularization phenotypes, confirming a disturbance of the source-sink endosperm phase change. We demonstrate that the hybrid barrier is under the influence of abiotic factors, and show that a temperature gradient leads to diametrically opposed cellularization phenotype responses in hybrid endosperm with A. arenosa or A. lyrata as pollen donors. Furthermore, different A. thaliana accession genotypes also enhance or diminish seed viability in the two hybrid cross-types, emphasizing that both genetic and environmental cues control the hybridization barrier. We have identified an A. thaliana MADS-BOX type I family single locus that is required for diametrically opposed cellularization phenotype responses in hybrid endosperm. Loss of AGAMOUS-LIKE 35 significantly affects the germination rate of hybrid seeds in opposite directions when transmitted through the A. thaliana endosperm, and is suggested to be a locus that promotes cellularization as part of an endosperm based mechanism involved in post-zygotic hybrid barriers. The role of temperature in hybrid speciation and the identification of distinct loci in control of hybrid failure have great potential to aid the introduction of advantageous traits in breeding research and to support models to predict hybrid admixture in a changing global climate.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

New insights into polyploid evolution and dynamic nature of Ludwigia section Isnardia (Onagraceae)

BACKGROUND

While polyploids are common in plants, the evolutionary history and natural dynamics of most polyploid groups are still unclear. Owing to plentiful earlier systematic studies, Ludwigia sect. Isnardia (comprising 22 wetland taxa) is an ideal allopolyploid complex to investigate polyploid evolution and natural dynamics within and among taxa. With a considerable sampling, we concentrated on revisiting earlier phylogenies of Isnardia, reevaluating the earlier estimated age of the most recent common ancestor (TMRCA), exploring the correlation between infraspecific genetic diversity and ploidy levels, and inspecting interspecific gene flows among taxa.

RESULTS

Phylogenetic trees and network concurred with earlier phylogenies and hypothesized genomes by incorporating 192 atpB-rbcL and ITS sequences representing 91% of Isnardia taxa. Moreover, we detected three multi-origin taxa. Our findings on L. repens and L. sphaerocarpa were consistent with earlier studies; L. arcuata was reported as a multi-origin taxon here, and an additional evolutionary scenario of L. sphaerocarpa was uncovered, both for the first time. Furthermore, estimated Isnardia TMRCA ages based on our data (5.9 or 8.9 million years ago) are in accordance with earlier estimates, although younger than fossil dates (Middle Miocene). Surprisingly, infraspecific genetic variations of Isnardia taxa did not increase with ploidy levels as anticipated from many other polyploid groups. In addition, the exuberant, low, and asymmetrical gene flows among Isnardia taxa indicated that the reproductive barriers may be weakened owing to allopolyploidization, which has rarely been reported.

CONCLUSIONS

The present research gives new perceptions of the reticulate evolution and dynamic nature of Isnardia and points to gaps in current knowledge about allopolyploid evolution.

Allopolyploid Speciation Accompanied by Gene Flow in a Tree Fern

Hybridization in plants may result in hybrid speciation or introgression and, thus, is now widely understood to be an important mechanism of species diversity on an evolutionary timescale. Hybridization is particularly common in ferns, as is polyploidy, which often results from hybrid crosses. Nevertheless, hybrid speciation as an evolutionary process in fern lineages remains poorly understood. Here, we employ flow cytometry, phylogeny, genomewide single nucleotide polymorphism data sets, and admixture and coalescent modeling to show that the scaly tree fern, Gymnosphaera metteniana is a naturally occurring allotetraploid species derived from hybridization between the diploids, G. denticulata and G. gigantea. Moreover, we detected ongoing gene flow between the hybrid species and its progenitors, and we found that G. gigantea and G. metteniana inhabit distinct niches, whereas climatic niches of G. denticulata and G. metteniana largely overlap. Taken together, these results suggest that either some degree of intrinsic genetic isolation between the hybrid species and its parental progenitors or ecological isolation over short distances may be playing an important role in the evolution of reproductive barriers. Historical climate change may have facilitated the origin of G. metteniana, with the timing of hybridization coinciding with a period of intensification of the East Asian monsoon during the Pliocene and Pleistocene periods in southern China. Our study of allotetraploid G. metteniana represents the first genomic-level documentation of hybrid speciation in scaly tree ferns and, thus, provides a new perspective on evolution in the lineage.

Comparative Phylogeography of Veronica spicata and V. longifolia (Plantaginaceae) Across Europe: Integrating Hybridization and Polyploidy in Phylogeography

Climatic fluctuations in the Pleistocene caused glacial expansion-contraction cycles in Eurasia and other parts of the world. Consequences of these cycles, such as population expansion and subsequent subdivision, have been studied in many taxa at intraspecific population level across much of the Northern Hemisphere. However, the consequences for the potential of hybridization and polyploidization are poorly understood. Here, we investigated the phylogeographic structure of two widespread, closely related species, Veronica spicata and Veronica longifolia, across their European distribution ranges. We assessed the extent and the geographic pattern of polyploidization in both species and hybridization between them. We used genome-scale SNP data to clarify phylogenetic relationships and detect possible hybridization/introgression events. In addition, crossing experiments were performed in different combination between V. spicata and V. longifolia individuals of two ploidy levels and of different geographic origins. Finally, we employed ecological niche modeling to infer macroclimatic differences between both species and both ploidy levels. We found a clear genetic structure reflecting the geographical distribution patterns in both species, with V. spicata showing higher genetic differentiation than V. longifolia. We retrieved significant signals of hybridization and introgression in natural populations from the genetic data and corroborated this with crossing experiments. However, there were no clear phylogeographic patterns and unequivocal macroclimatic niche differences between diploid and tetraploid lineages. This favors the hypothesis, that autopolyploidization has happened frequently and in different regions. The crossing experiments produced viable hybrids when the crosses were made between plants of the same ploidy levels but not in the interploidy crosses. The results suggest that hybridization occurs across the overlapping areas of natural distribution ranges of both species, with apparently directional introgression from V. spicata to V. longifolia. Nevertheless, the two species maintain their species-level separation due to their adaptation to different habitats and spatial isolation rather than reproductive isolation.

Population Genetics of Odontarrhena (Brassicaceae) from Albania: The Effects of Anthropic Habitat Disturbance, Soil, and Altitude on a Ni-Hyperaccumulator Plant Group from a Major Serpentine Hotspot

Albanian taxa and populations of the genus Odontarrhena are most promising candidates for research on metal tolerance and Ni-agromining, but their genetic structure remains unknown. We investigated phylogenetic relationships and genetic differentiation in relation to distribution and ploidy of the taxa, anthropic site disturbance, elevation, soil type, and trace metals at each population site. After performing DNA sequencing of selected accessions, we applied DNA-fingerprinting to analyze the genetic structure of 32 populations from ultramafic and non-ultramafic outcrops across Albania. Low sequence divergence resulted in poorly resolved phylograms, but supported affinity between the two diploid serpentine endemics O. moravensis and O. rigida. Analysis of molecular variance (AMOVA) revealed significant population differentiation, but no isolation by distance. Among-population variation was higher in polyploids than in diploids, in which genetic distances were lower. Genetic admixing at population and individual level occurred especially in the polyploids O. chalcidica, O. decipiens, and O. smolikana. Admixing increased with site disturbance. Outlier loci were higher in serpentine populations but decreased along altitude with lower drought and heat stress. Genetic variability gained by gene flow and hybridization at contact zones with "resident" species of primary ultramafic habitats promoted expansion of the tetraploid O. chalcidica across anthropogenic sites.

Maintenance of Species Differences in Closely Related Tetraploid Parasitic Euphrasia (Orobanchaceae) on an Isolated Island

Polyploidy is pervasive in angiosperm evolution and plays important roles in adaptation and speciation. However, polyploid groups are understudied due to complex sequence homology, challenging genome assembly, and taxonomic complexity. Here, we study adaptive divergence in taxonomically complex eyebrights (Euphrasia), where recent divergence, phenotypic plasticity, and hybridization blur species boundaries. We focus on three closely related tetraploid species with contrasting ecological preferences that are sympatric on Fair Isle, a small isolated island in the British Isles. Using a common garden experiment, we show a genetic component to the morphological differences present between these species. Using whole-genome sequencing and a novel k-mer approach we call "Tetmer", we demonstrate that the species are of allopolyploid origin, with a sub-genome divergence of approximately 5%. Using ∼2 million SNPs, we show sub-genome homology across species, with a very low sequence divergence characteristic of recent speciation. This genetic variation is broadly structured by species, with clear divergence of Fair Isle heathland Euphrasia micrantha, while grassland Euphrasia arctica and coastal Euphrasia foulaensis are more closely related. Overall, we show that tetraploid Euphrasia is a system of allopolyploids of postglacial species divergence, where adaptation to novel environments may be conferred by old variants rearranged into new genetic lineages.

Genetic similarities versus morphological resemblance: Unraveling a polyploid complex in a Mediterranean biodiversity hotspot

The Balkan Peninsula is recognized as one of the hotspots of biodiversity in Europe. This area has shown since the Last Glacial Maximum appropriate conditions for species diversification and hybridization, which has led to the existence of numerous taxonomically unresolved entities. Here, we focus on the Western Balkans and explore the genetic structure and relationships among species belonging to the V. austriaca - V. orbiculata diploid-polyploid complex, including populations showing intermediate morphologies. A combination of nuclear markers (microsatellites), plastid DNA regions (trnH-psbA, ycf6-psbM) and ploidy level estimations using flow cytometry are employed to assess the genetic structure and evolutionary dynamics of this polyploid complex. To reconstruct the evolutionary history, an approximate Bayesian computation approach is combined with projections of the species distribution models onto the climatic scenarios of the Mid-Holocene (6 ka BP) and Last Glacial Maximum (22 ka BP). Four main groups were found: one well-established entity within the diploid level, V. dalmatica, a second diploid-tetraploid group which corresponds to V. orbiculata, a hexaploid cluster harboring V. austriaca subsp. jacquinii individuals, and an enigmatic tetraploid group. According to the molecular data obtained, this latter cluster represents an allopolyploid cryptic lineage −with V. orbiculata and V. dalmatica as putative parents− morphologically similar to V. orbiculata, but genetically more related to V. austriaca subsp. jacquinii. Veronica dalmatica and this “uncertain tetraploid” group are involved in the formation of the hexaploid taxon V. austriaca subsp. jacquinii, with the possibility of recent gene flow among different cytotypes. The present study supports a scenario of diversification from a diploid common ancestor leading to two different but interrelated lineages. The first one would correspond with the diploid V. orbiculata plus tetraploid individuals of this species arising through allo- and autopolyploidization, and the second one would involve all ploidy levels with allopolyploidization being prevalent.

A novel allele of ASY3 is associated with greater meiotic stability in autotetraploid Arabidopsis lyrata

In this study we performed a genotype-phenotype association analysis of meiotic stability in 10 autotetraploid Arabidopsis lyrata and A. lyrata/A. arenosa hybrid populations collected from the Wachau region and East Austrian Forealps. The aim was to determine the effect of eight meiosis genes under extreme selection upon adaptation to whole genome duplication. Individual plants were genotyped by high-throughput sequencing of the eight meiosis genes (ASY1, ASY3, PDS5b, PRD3, REC8, SMC3, ZYP1a/b) implicated in synaptonemal complex formation and phenotyped by assessing meiotic metaphase I chromosome configurations. Our results reveal that meiotic stability varied greatly (20-100%) between individual tetraploid plants and associated with segregation of a novel ASYNAPSIS3 (ASY3) allele derived from A. lyrata. The ASY3 allele that associates with meiotic stability possesses a putative in-frame tandem duplication (TD) of a serine-rich region upstream of the coiled-coil domain that appears to have arisen at sites of DNA microhomology. The frequency of multivalents observed in plants homozygous for the ASY3 TD haplotype was significantly lower than in plants heterozygous for ASY3 TD/ND (non-duplicated) haplotypes. The chiasma distribution was significantly altered in the stable plants compared to the unstable plants with a shift from proximal and interstitial to predominantly distal locations. The number of HEI10 foci at pachytene that mark class I crossovers was significantly reduced in a plant homozygous for ASY3 TD compared to a plant heterozygous for ASY3 ND/TD. Fifty-eight alleles of the 8 meiosis genes were identified from the 10 populations analysed, demonstrating dynamic population variability at these loci. Widespread chimerism between alleles originating from A. lyrata/A. arenosa and diploid/tetraploids indicates that this group of rapidly evolving genes may provide precise adaptive control over meiotic recombination in the tetraploids, the very process that gave rise to them.

Genetic variation and temperature affects hybrid barriers during interspecific hybridization

Genomic imprinting regulates parent-specific transcript dosage during seed development and is mainly confined to the endosperm. Elucidation of the function of many imprinted genes has been hampered by the lack of corresponding mutant phenotypes, and the role of imprinting is mainly associated with genome dosage regulation or allocation of resources. Disruption of imprinted genes has also been suggested to mediate endosperm-based post-zygotic hybrid barriers depending on genetic variation and gene dosage. Here, we have analyzed the conservation of a clade from the MADS-box type I class transcription factors in the closely related species Arabidopsis arenosa, A. lyrata, and A. thaliana, and show that AGL36-like genes are imprinted and maternally expressed in seeds of Arabidopsis species and in hybrid seeds between outbreeding species. In hybridizations between outbreeding and inbreeding species the paternally silenced allele of the AGL36-like gene is reactivated in the hybrid, demonstrating that also maternally expressed imprinted genes are perturbed during hybridization and that such effects on imprinted genes are specific to the species combination. Furthermore, we also demonstrate a quantitative effect of genetic diversity and temperature on the strength of the post-zygotic hybridization barrier. Markedly, a small decrease in temperature during seed development increases the survival of hybrid F1 seeds, suggesting that abiotic and genetic parameters play important roles in post-zygotic species barriers, pointing at evolutionary scenarios favoring such effects. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA562212. All sequences generated in this study have been deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with project number PRJNA562212.

Interspecific introgression mediates adaptation to whole genome duplication

Adaptive gene flow is a consequential phenomenon across all kingdoms. Although recognition is increasing, there is no study showing that bidirectional gene flow mediates adaptation at loci that manage core processes. We previously discovered concerted molecular changes among interacting members of the meiotic machinery controlling crossover number upon adaptation to whole-genome duplication (WGD) in Arabidopsis arenosa. Here we conduct a population genomic study to test the hypothesis that adaptation to WGD has been mediated by adaptive gene flow between A. arenosa and A. lyrata. We find that A. lyrata underwent WGD more recently than A. arenosa, suggesting that pre-adapted alleles have rescued nascent A. lyrata, but we also detect gene flow in the opposite direction at functionally interacting loci under the most extreme levels of selection. These data indicate that bidirectional gene flow allowed for survival after WGD, and that the merger of these species is greater than the sum of their parts.

Whole genome duplication (WGD) presents new challenges to the establishment of optimal allelic combinations and to the meiotic machinery. Here, the authors show that adaptive gene flow from Arabidopsis arenosa could rescue the nascent A. lyrata from extinction following WGD.

Autopolyploid lineage shows climatic niche expansion but not divergence in Arabidopsis arenosa

PREMISE OF THE STUDY

Successful establishment of neopolyploids, and therefore polyploid speciation, is thought to be contingent on environmental niche shifts from their progenitors. We explore this niche shift hypothesis in the obligate outcrosser Arabidopsis arenosa complex, which includes diploid and recently formed autotetraploid populations.

METHODS

To characterize the climatic niches for both cytotypes in Arabidopsis arenosa, we first gathered climatic data from localities with known ploidy types. We then estimated the climatic niches for diploids and autotetraploids and calculated niche overlap. Using this niche overlap statistic, we tested for niche equivalency and similarity. We explored differences in niches by estimating and comparing niche optimum and breadth and then calculated indices of niche expansion and unfilling.

KEY RESULTS

Climatic niche overlap between diploids and autotetraploids is substantial. Although the two niche models are not significantly divergent, they are not identical as they differ in both optimum and breadth along two environmental gradients. Autotetraploids fill nearly the entire niche space of diploids and have expanded into novel environments.

CONCLUSIONS

We find climatic niche expansion but not divergence, together with a moderate change in the niche optimum, in the autotetraploid lineage of Arabidopsis arenosa. These results indicate that the climatic niche shift hypothesis alone cannot explain the coexistence of tetraploid and diploid cytotypes.

Genetic basis and evolution of rapid cycling in railway populations of tetraploid Arabidopsis arenosa

Spatially structured plant populations with diverse adaptations provide powerful models to investigate evolution. Human-generated ruderal habitats are abundant and low-competition, but are challenging for plants not adapted to them. Ruderal habitats also sometimes form networked corridors (e.g. roadsides and railways) that allow rapid long-distance spread of successfully adapted variants. Here we use transcriptomic and genomic analyses, coupled with genetic mapping and transgenic follow-up, to understand the evolution of rapid cycling during adaptation to railway sites in autotetraploid Arabidopsis arenosa. We focus mostly on a hybrid population that is likely a secondary colonist of a railway site. These mountain railway plants are phenotypically similar to their cosmopolitan cousins. We thus hypothesized that colonization primarily involved the flow of adaptive alleles from the cosmopolitan railway variant. But our data shows that it is not that simple: while there is evidence of selection having acted on introgressed alleles, selection also acted on rare standing variation, and new mutations may also contribute. Among the genes we show have allelic divergence with functional relevance to flowering time are known regulators of flowering, including FLC and CONSTANS. Prior implications of these genes in weediness and rapid cycling supports the idea that these are “evolutionary hotspots” for these traits. We also find that one of two alleles of CONSTANS under selection in the secondary colonist was selected from rare standing variation in mountain populations, while the other was introgressed from the cosmopolitan railway populations. The latter allele likely arose in diploid populations over 700km away, highlighting how ruderal populations could act as allele conduits and thus influence local adaptation.

Polyploid Arabidopsis species originated around recent glaciation maxima

Polyploidy may provide adaptive advantages and is considered to be important for evolution and speciation. Polyploidy events are found throughout the evolutionary history of plants, however they do not seem to be uniformly distributed along the time axis. For example, many of the detected ancient whole-genome duplications (WGDs) seem to cluster around the K/Pg boundary (∼66Mya), which corresponds to a drastic climate change event and a mass extinction. Here, we discuss more recent polyploidy events using Arabidopsis as the most developed plant model at the level of the entire genus. We review the history of the origin of allotetraploid species A. suecica and A. kamchatica, and tetraploid lineages of A. lyrata, A. arenosa and A. thaliana, and discuss potential adaptive advantages. Also, we highlight an association between recent glacial maxima and estimated times of origins of polyploidy in Arabidopsis. Such association might further support a link between polyploidy and environmental challenge, which has been observed now for different time-scales and for both ancient and recent polyploids.

An Arabidopsis introgression zone studied at high spatio-temporal resolution: interglacial and multiple genetic contact exemplified using whole nuclear and plastid genomes

Background

Gene flow between species, across ploidal levels, and even between evolutionary lineages is a common phenomenon in the genus Arabidopsis. However, apart from two genetically fully stabilized allotetraploid species that have been investigated in detail, the extent and temporal dynamics of hybridization are not well understood. An introgression zone, with tetraploid A. arenosa introgressing into A. lyrata subsp. petraea in the Eastern Austrian Forealps and subsequent expansion towards pannonical lowlands, was described previously based on morphological observations as well as molecular data using microsatellite and plastid DNA markers. Here we investigate the spatio-temporal context of this suture zone, making use of the potential of next-generation sequencing and whole-genome data. By utilizing a combination of nuclear and plastid genomic data, the extent, direction and temporal dynamics of gene flow are elucidated in detail and Late Pleistocene evolutionary processes are resolved.

Results

Analysis of nuclear genomic data significantly recognizes the clinal structure of the introgression zone, but also reveals that hybridization and introgression is more common and substantial than previously thought. Also tetraploid A. lyrata and A. arenosa subsp. borbasii from outside the previously defined suture zone show genomic signals of past introgression. A. lyrata is shown to serve usually as the maternal parent in these hybridizations, but one exception is identified from plastome-based phylogenetic reconstruction. Using plastid phylogenomics with secondary time calibration, the origin of A. lyrata and A. arenosa lineages is pre-dating the last three glaciation complexes (approx. 550,000 years ago). Hybridization and introgression followed during the last two glacial-interglacial periods (since approx. 300,000 years ago) with later secondary contact at the northern and southern border of the introgression zone during the Holocene.

Conclusions

Footprints of adaptive introgression in the Northeastern Forealps are older than expected and predate the Last Glaciation Maximum. This correlates well with high genetic diversity found within areas that served as refuge area multiple times. Our data also provide some first hints that early introgressed and presumably preadapted populations account for successful and rapid postglacial re-colonization and range expansion.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4220-6) contains supplementary material, which is available to authorized users.

Abnormal meiosis in an intersectional allotriploid of Populus L. and segregation of ploidy levels in 2x × 3x progeny

Triploid plants are usually highly aborted owing to unbalanced meiotic chromosome segregation, but limited viable gametes can participate in the transition to different ploidy levels. In this study, numerous meiotic abnormalities were found with high frequency in an intersectional allotriploid poplar (Populus alba × P. berolinensis 'Yinzhong'), including univalents, precocious chromosome migration, lagging chromosomes, chromosome bridges, micronuclei, and precocious cytokinesis, indicating high genetic imbalance in this allotriploid. Some micronuclei trigger mini-spindle formation in metaphase II and participate in cytokinesis to form polyads with microcytes. Unbalanced chromosome segregation and chromosome elimination resulted in the formation of microspores with aneuploid chromosome sets. Fusion of sister nuclei occurs in microsporocytes with precocious cytokinesis, which could form second meiotic division restitution (SDR)-type gametes. However, SDR-type gametes likely contain incomplete chromosome sets due to unbalanced segregation of homologous chromosomes during the first meiotic division in triploids. Misorientation of spindles during the second meiotic division, such as fused and tripolar spindles with low frequency, could result in the formation of first meiotic division restitution (FDR)-type unreduced gametes, which most likely contain three complete chromosome sets. Although 'Yinzhong' yields 88.7% stainable pollen grains with wide diameter variation from 23.9 to 61.3 μm, the pollen viability is poor (2.78% ± 0.38). A cross of 'Yinzhong' pollen with a diploid female clone produced progeny with extensive segregation of ploidy levels, including 29 diploids, 18 triploids, 4 tetraploids, and 48 aneuploids, suggesting the formation of viable aneuploidy and unreduced pollen in 'Yinzhong'. Individuals with different chromosome compositions are potential to analyze chromosomal function and to integrate the chromosomal dosage variation into breeding programs of Populus.

Endosperm-based hybridization barriers explain the pattern of gene flow between Arabidopsis lyrata and Arabidopsis arenosa in Central Europe

Based on the biological species concept, two species are considered distinct if reproductive barriers prevent gene flow between them. In Central Europe, the diploid species Arabidopsis lyrata and Arabidopsis arenosa are genetically isolated, thus fitting this concept as "good species." Nonetheless, interspecific gene flow involving their tetraploid forms has been described. The reasons for this ploidy-dependent reproductive isolation remain unknown. Here, we show that hybridization between diploid A. lyrata and A. arenosa causes mainly inviable seed formation, revealing a strong postzygotic reproductive barrier separating these two species. Although viability of hybrid seeds was impaired in both directions of hybridization, the cause for seed arrest differed. Hybridization of A. lyrata seed parents with A. arenosa pollen donors resulted in failure of endosperm cellularization, whereas the endosperm of reciprocal hybrids cellularized precociously. Endosperm cellularization failure in both hybridization directions is likely causal for the embryo arrest. Importantly, natural tetraploid A. lyrata was able to form viable hybrid seeds with diploid and tetraploid A. arenosa, associated with the reestablishment of normal endosperm cellularization. Conversely, the defects of hybrid seeds between tetraploid A. arenosa and diploid A. lyrata were aggravated. According to these results, we hypothesize that a tetraploidization event in A. lyrata allowed the production of viable hybrid seeds with A. arenosa, enabling gene flow between the two species.

Polyploidy can Confer Superiority to West African Acacia senegal (L.) Willd. Trees

Polyploidy is a common phenomenon in the evolution of angiosperms. It has been suggested that polyploids manage harsh environments better than their diploid relatives but empirical data supporting this hypothesis are scarce, especially for trees. Using microsatellite markers and flow cytometry, we examine the frequency of polyploids and diploids in a progeny trial testing four different populations of Acacia senegal, a species native to sub-Saharan regions of Africa. We compare growth between cytotypes and test whether polyploid seedlings grow better than diploids. Our results show that polyploids coexist with diploids in highly variable proportions among populations in Senegal. Acacia senegal genotypes were predominantly diploid and tetraploid, but triploid, pentaploid, hexaploid, and octaploid forms were also found. We find that polyploids show faster growth than diploids under our test conditions: in an 18 years old field trial, polyploid superiority was estimated to be 17% in trunk diameter and 9% in height while in a growth chamber experiment, polyploids grew 28% taller, but only after being exposed to drought stress. The results suggest that polyploid A. senegal can have an adaptive advantage in some regions of Africa.

Unidirectional diploid-tetraploid introgression among British birch trees with shifting ranges shown by restriction site-associated markers

Hybridization may lead to introgression of genes among species. Introgression may be bidirectional or unidirectional, depending on factors such as the demography of the hybridizing species, or the nature of reproductive barriers between them. Previous microsatellite studies suggested bidirectional introgression between diploid Betula nana (dwarf birch) and tetraploid B. pubescens (downy birch) and also between B. pubescens and diploid B. pendula (silver birch) in Britain. Here, we analyse introgression among these species using 51 237 variants in restriction site‐associated (RAD) markers in 194 individuals, called with allele dosages in the tetraploids. In contrast to the microsatellite study, we found unidirectional introgression into B. pubescens from both of the diploid species. This pattern fits better with the expected nature of the reproductive barrier between diploids and tetraploids. As in the microsatellite study, introgression into B. pubescens showed clear clines with increasing introgression from B. nana in the north and from B. pendula in the south. Unlike B. pendula alleles, introgression of B. nana alleles was found far from the current area of sympatry or allopatry between B. nana and B. pubescens. This pattern fits a shifting zone of hybridization due to Holocene reduction in the range of B. nana and expansion in the range of B. pubescens.

Habitat-Associated Life History and Stress-Tolerance Variation in Arabidopsis arenosa

Weediness in ephemeral plants is commonly characterized by rapid cycling, prolific all-in flowering, and loss of perenniality. Many species made transitions to weediness of this sort, which can be advantageous in high-disturbance or human-associated habitats. The molecular basis of this shift, however, remains mostly mysterious. Here, we use transcriptome sequencing, genome resequencing scans for selection, and stress tolerance assays to study a weedy population of the otherwise nonweedy Arabidopsis arenosa, an obligately outbreeding relative of Arabidopsis thaliana Although weedy A. arenosa is widespread, a single genetic lineage colonized railways throughout central and northern Europe. We show that railway plants, in contrast to plants from sheltered outcrops in hill/mountain regions, are rapid cycling, have lost the vernalization requirement, show prolific flowering, and do not return to vegetative growth. Comparing transcriptomes of railway and mountain plants across time courses with and without vernalization, we found that railway plants have sharply abrogated vernalization responsiveness and high constitutive expression of heat- and cold-responsive genes. Railway plants also have strong constitutive heat shock and freezing tolerance compared with mountain plants, where tolerance must be induced. We found 20 genes with good evidence of selection in the railway population. One of these, LATE ELONGATED HYPOCOTYL, is known in A. thaliana to regulate many stress-response genes that we found to be differentially regulated among the distinct habitats. Our data suggest that, beyond life history regulation, other traits like basal stress tolerance also are associated with the evolution of weediness in A. arenosa.

Phenotypic and Genotypic Analysis of Newly Obtained Interspecific Hybrids in the Campanula Genus

Interspecific hybridisation creates new phenotypes within several ornamental plant species including the Campanula genus. We have employed phenotypic and genotypic methods to analyse and evaluate interspecific hybridisation among cultivars of four Campanula species, i.e. C. cochleariifolia, C. isophylla, C. medium and C. formanekiana. Hybrids were analysed using amplified fragment length polymorphism (AFLP), flow cytometry and biometrical measurements. Results of correlation matrices demonstrated heterogeneous phenotypes for the parental species, which confirmed our basic premise for new phenotypes of interspecific hybrids. AFLP assays confirmed the hybridity and identified self-pollinated plants. Limitation of flow cytometry analysis detection was observed while detecting the hybridity status of two closely related parents, e.g. C. cochleariiafolia × C. isophylla. Phenotypic characteristics such as shoot habitus and flower colour were strongly influenced by one of the parental species in most crosses. Rooting analysis revealed that inferior rooting quality occurred more often in interspecific hybrids than in the parental species. Only interspecific hybrid lines of C. formanekiana ‘White’ × C. medium ‘Pink’ showed a high rooting level. Phenotype analyses demonstrated a separation from the interspecific hybrid lines of C. formanekiana ‘White’ × C. medium ‘Pink’ to the other clustered hybrids of C. formanekiana and C. medium. In our study we demonstrated that the use of correlation matrices is a suitable tool for identifying suitable cross material. This study presents a comprehensive overview for analysing newly obtained interspecific hybrids. The chosen methods can be used as guidance for analyses for further interspecific hybrids in Campanula, as well as in other ornamental species.

Frequent introgressions from diploid species contribute to the adaptation of the tetraploid Shepherd's purse (Capsella bursa-pastoris)

Adaptation is the most important ability for organisms to survive in diverse habitats. Animals have the option to escape from stressful environments, but plants do not. In plants, polyploids consist of about 30%-70% angiosperms and 95% ferns, of which some are important crops such as cotton and wheat. How polyploid plants adapt to various habitats has been a fundamental question remained largely unanswered. The tetraploid Shepherd's purse (Capsella bursa-pastoris) is one of the most successful plants on earth and has been distributed across the world, thus being an ideal model system for studying the adaptation of polyploids. We found that there are frequent introgressions from congeneric diploids to Shepherd's purse. Ecological niche modeling suggests that ecological differentiation is evident between the introgressed and non-introgressed C. bursa-pastoris, and the introgressions are a source of adaptation. This result links an evolutionary process to the adaptation of polyploids, and sheds light on the breeding strategy of polyploids as well. We conclude that frequent introgressions from congeneric diploids contributed to the acquisition of adequate genetic variations, thereby allowing C. bursa-pastoris to adapt to various habitats across the world. Our results highlight how a polyploid could have successfully established after it originated.

Epigenetic mechanisms of postzygotic reproductive isolation in plants

Hybrid incompatibility is generally viewed as a consequence of negative epistatic interactions between alleles that do not cause negative fitness effects in their parents. Substantial evidence in support of the model has accumulated over recent years. Nevertheless, there is evidence that hybrid incompatibility can have an epigenetic basis and results from deregulated small RNAs (sRNAs), causing changes in DNA methylation and transposable element (TE) activation. Epigenetically regulated loci can impact on the expression of nearby located genes. Alteration of dosage-sensitive gene expression builds hybridization barriers in the endosperm; however, it may also offer an explanation for transgressive effects in plant hybrids. In this review we highlight recent advances that illuminate the role of epigenetic pathways in establishing hybrid incompatibility in plants.

Selection on meiosis genes in diploid and tetraploid Arabidopsis arenosa

Meiotic chromosome segregation is critical for fertility across eukaryotes, and core meiotic processes are well conserved even between kingdoms. Nevertheless, recent work in animals has shown that at least some meiosis genes are highly diverse or strongly differentiated among populations. What drives this remains largely unknown. We previously showed that autotetraploid Arabidopsis arenosa evolved stable meiosis, likely through reduced crossover rates, and that associated with this there is strong evidence for selection in a subset of meiosis genes known to affect axis formation, synapsis, and crossover frequency. Here, we use genome-wide data to study the molecular evolution of 70 meiosis genes in a much wider sample of A. arenosa. We sample the polyploid lineage, a diploid lineage from the Carpathian Mountains, and a more distantly related diploid lineage from the adjacent, but biogeographically distinct Pannonian Basin. We find that not only did selection act on meiosis genes in the polyploid lineage but also independently on a smaller subset of meiosis genes in Pannonian diploids. Functionally related genes are targeted by selection in these distinct contexts, and in two cases, independent sweeps occurred in the same loci. The tetraploid lineage has sustained selection on more genes, has more amino acid changes in each, and these more often affect conserved or potentially functional sites. We hypothesize that Pannonian diploid and tetraploid A. arenosa experienced selection on structural proteins that mediate sister chromatid cohesion, the formation of meiotic chromosome axes, and synapsis, likely for different underlying reasons.

Loss and retention of resistance genes in five species of the Brassicaceae family

BACKGROUND

Plants have evolved disease resistance (R) genes encoding for nucleotide-binding site (NB) and leucine-rich repeat (LRR) proteins with N-terminals represented by either Toll/Interleukin-1 receptor (TIR) or coiled-coil (CC) domains. Here, a genome-wide study of presence and diversification of CC-NB-LRR and TIR-NB-LRR encoding genes, and shorter domain combinations in 19 Arabidopsis thaliana accessions and Arabidopsis lyrata, Capsella rubella, Brassica rapa and Eutrema salsugineum are presented.

RESULTS

Out of 528 R genes analyzed, 12 CC-NB-LRR and 17 TIR-NB-LRR genes were conserved among the 19 A. thaliana genotypes, while only two CC-NB-LRRs, including ZAR1, and three TIR-NB-LRRs were conserved when comparing the five species. The RESISTANCE TO LEPTOSPHAERIA MACULANS 1 (RLM1) locus confers resistance to the Brassica pathogen L. maculans the causal agent of blackleg disease and has undergone conservation and diversification events particularly in B. rapa. On the contrary, the RLM3 locus important in the immune response towards Botrytis cinerea and Alternaria spp. has recently evolved in the Arabidopsis genus.

CONCLUSION

Our genome-wide analysis of the R gene repertoire revealed a large sequence variation in the 23 cruciferous genomes. The data provides further insights into evolutionary processes impacting this important gene family.

Taming the wild: resolving the gene pools of non-model Arabidopsis lineages

BACKGROUND

Wild relatives in the genus Arabidopsis are recognized as useful model systems to study traits and evolutionary processes in outcrossing species, which are often difficult or even impossible to investigate in the selfing and annual Arabidopsis thaliana. However, Arabidopsis as a genus is littered with sub-species and ecotypes which make realizing the potential of these non-model Arabidopsis lineages problematic. There are relatively few evolutionary studies which comprehensively characterize the gene pools across all of the Arabidopsis supra-groups and hypothesized evolutionary lineages and none include sampling at a world-wide scale. Here we explore the gene pools of these various taxa using various molecular markers and cytological analyses.

RESULTS

Based on ITS, microsatellite, chloroplast and nuclear DNA content data we demonstrate the presence of three major evolutionary groups broadly characterized as A. lyrata group, A. halleri group and A. arenosa group. All are composed of further species and sub-species forming larger aggregates. Depending on the resolution of the marker, a few closely related taxa such as A. pedemontana, A. cebennensis and A. croatica are also clearly distinct evolutionary lineages. ITS sequences and a population-based screen based on microsatellites were highly concordant. The major gene pools identified by ITS sequences were also significantly differentiated by their homoploid nuclear DNA content estimated by flow cytometry. The chloroplast genome provided less resolution than the nuclear data, and it remains unclear whether the extensive haplotype sharing apparent between taxa results from gene flow or incomplete lineage sorting in this relatively young group of species with Pleistocene origins.

CONCLUSIONS

Our study provides a comprehensive overview of the genetic variation within and among the various taxa of the genus Arabidopsis. The resolved gene pools and evolutionary lineages will set the framework for future comparative studies on genetic diversity. Extensive population-based phylogeographic studies will also be required, however, in particular for A. arenosa and their affiliated taxa and cytotypes.

Genome size variation and evolution in allotetraploid Arabidopsis kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri

Polyploidization and subsequent changes in genome size are fundamental processes in evolution and diversification. Little is currently known about the extent of genome size variation within taxa and the evolutionary forces acting on this variation. Arabidopsis kamchatica has been reported to contain both diploid and tetraploid individuals. The aim of this study was to determine the genome size of A. kamchatica, whether there is variation in ploidy and/or genome size in A. kamchatica and to study how genome size has evolved. We used propidium iodide flow cytometry to measure 2C DNA content of 73 plants from 25 geographically diverse populations of the putative allotetraploid A. kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri. All A. kamchatica plants appear to be tetraploids. The mean 2C DNA content of A. kamchatica was 1.034 pg (1011 Mbp), which is slightly smaller than the sum of its diploid parents (A. lyrata: 0.502 pg; A. halleri: 0.571 pg). Arabidopsis kamchatica appears to have lost ∼37.594 Mbp (3.6 %) of DNA from its 2C genome. Tetraploid A. lyrata from Germany and Austria appears to have lost ∼70.366 Mbp (7.2 %) of DNA from the 2C genome, possibly due to hybridization with A. arenosa, which has a smaller genome than A. lyrata. We did find genome size differences among A. kamchatica populations, which varied up to 7 %. Arabidopsis kamchatica ssp. kawasakiana from Japan appears to have a slightly larger genome than A. kamchatica ssp. kamchatica from North America, perhaps due to multiple allopolyploid origins or hybridization with A. halleri. However, the among-population coefficient of variation in 2C DNA content is lower in A. kamchatica than in other Arabidopsis taxa. Due to its close relationship to A. thaliana, A. kamchatica has the potential to be very useful in the study of polyploidy and genome evolution.

Phylogeny of a genomically diverse group of elymus (poaceae) allopolyploids reveals multiple levels of reticulation

The grass tribe Triticeae (=Hordeeae) comprises only about 300 species, but it is well known for the economically important crop plants wheat, barley, and rye. The group is also recognized as a fascinating example of evolutionary complexity, with a history shaped by numerous events of auto- and allopolyploidy and apparent introgression involving diploids and polyploids. The genus Elymus comprises a heterogeneous collection of allopolyploid genome combinations, all of which include at least one set of homoeologs, designated St, derived from Pseudoroegneria. The current analysis includes a geographically and genomically diverse collection of 21 tetraploid Elymus species, and a single hexaploid species. Diploid and polyploid relationships were estimated using four molecular data sets, including one that combines two regions of the chloroplast genome, and three from unlinked nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Four gene trees were generated using maximum likelihood, and the phylogenetic placement of the polyploid sequences reveals extensive reticulation beyond allopolyploidy alone. The trees were interpreted with reference to numerous phenomena known to complicate allopolyploid phylogenies, and introgression was identified as a major factor in their history. The work illustrates the interpretation of complicated phylogenetic results through the sequential consideration of numerous possible explanations, and the results highlight the value of careful inspection of multiple independent molecular phylogenetic estimates, with particular focus on the differences among them.

Meiotic adaptation to genome duplication in Arabidopsis arenosa

Whole genome duplication (WGD) is a major factor in the evolution of multicellular eukaryotes, yet by doubling the number of homologs, WGD severely challenges reliable chromosome segregation, a process conserved across kingdoms. Despite this, numerous genome-duplicated (polyploid) species persist in nature, indicating early problems can be overcome. Little is known about which genes are involved--only one has been molecularly characterized. To gain new insights into the molecular basis of adaptation to polyploidy, we investigated genome-wide patterns of differentiation between natural diploids and tetraploids of Arabidopsis arenosa, an outcrossing relative of A. thaliana. We first show that diploids are not preadapted to polyploid meiosis. We then use a genome scanning approach to show that although polymorphism is extensively shared across ploidy levels, there is strong ploidy-specific differentiation in 39 regions spanning 44 genes. These are discrete, mostly single-gene peaks of sharply elevated differentiation. Among these peaks are eight meiosis genes whose encoded proteins coordinate a specific subset of early meiotic functions, suggesting these genes comprise a polygenic solution to WGD-associated chromosome segregation challenges. Our findings indicate that even conserved meiotic processes can be capable of nimble evolutionary shifts when required.

Interploidy hybridization in sympatric zones: the formation of Epidendrum fulgens × E. puniceoluteum hybrids (Epidendroideae, Orchidaceae)

Interspecific hybridization is a primary cause of extensive morphological and chromosomal variation and plays an important role in plant species diversification. However, the role of interploidal hybridization in the formation of hybrid swarms is less clear. Epidendrum encompasses wide variation in chromosome number and lacks strong premating barriers, making the genus a good model for clarifying the role of chromosomes in postzygotic barriers in interploidal hybrids. In this sense, hybrids from the interploidal sympatric zone between E. fulgens (2n = 2x = 24) and E. puniceoluteum (2n = 4x = 56) were analyzed using cytogenetic techniques to elucidate the formation and establishment of interploidal hybrids. Hybrids were not a uniform group: two chromosome numbers were observed, with the variation being a consequence of severe hybrid meiotic abnormalities and backcrossing with E. puniceoluteum. The hybrids were triploids (2n = 3x = 38 and 40) and despite the occurrence of enormous meiotic problems associated with triploidy, the hybrids were able to backcross, producing successful hybrid individuals with broad ecological distributions. In spite of the nonpolyploidization of the hybrid, its formation is a long-term evolutionary process rather than a product of a recent disturbance, and considering other sympatric zones in Epidendrum, these events could be recurrent.

Genetic adaptation associated with genome-doubling in autotetraploid Arabidopsis arenosa

Genome duplication, which results in polyploidy, is disruptive to fundamental biological processes. Genome duplications occur spontaneously in a range of taxa and problems such as sterility, aneuploidy, and gene expression aberrations are common in newly formed polyploids. In mammals, genome duplication is associated with cancer and spontaneous abortion of embryos. Nevertheless, stable polyploid species occur in both plants and animals. Understanding how natural selection enabled these species to overcome early challenges can provide important insights into the mechanisms by which core cellular functions can adapt to perturbations of the genomic environment. Arabidopsis arenosa includes stable tetraploid populations and is related to well-characterized diploids A. lyrata and A. thaliana. It thus provides a rare opportunity to leverage genomic tools to investigate the genetic basis of polyploid stabilization. We sequenced the genomes of twelve A. arenosa individuals and found signatures suggestive of recent and ongoing selective sweeps throughout the genome. Many of these are at genes implicated in genome maintenance functions, including chromosome cohesion and segregation, DNA repair, homologous recombination, transcriptional regulation, and chromatin structure. Numerous encoded proteins are predicted to interact with one another. For a critical meiosis gene, ASYNAPSIS1, we identified a non-synonymous mutation that is highly differentiated by cytotype, but present as a rare variant in diploid A. arenosa, indicating selection may have acted on standing variation already present in the diploid. Several genes we identified that are implicated in sister chromatid cohesion and segregation are homologous to genes identified in a yeast mutant screen as necessary for survival of polyploid cells, and also implicated in genome instability in human diseases including cancer. This points to commonalities across kingdoms and supports the hypothesis that selection has acted on genes controlling genome integrity in A. arenosa as an adaptive response to genome doubling.

Extending coalescent theory to autotetraploids

We develop coalescent models for autotetraploid species with tetrasomic inheritance. We show that the ancestral genetic process in a large population without recombination may be approximated using Kingman's standard coalescent, with a coalescent effective population size 4N. Numerical results suggest that this approximation is accurate for population sizes on the order of hundreds of individuals. Therefore, existing coalescent simulation programs can be adapted to study population history in autotetraploids simply by interpreting the timescale in units of 4N generations. We also consider the possibility of double reduction, a phenomenon unique to polysomic inheritance, and show that its effects on gene genealogies are similar to partial self-fertilization.

The evolutionary history of the Arabidopsis arenosa complex: diverse tetraploids mask the Western Carpathian center of species and genetic diversity

The Arabidopsis arenosa complex is closely related to the model plant Arabidopsis thaliana. Species and subspecies in the complex are mainly biennial, predominantly outcrossing, herbaceous, and with a distribution range covering most parts of latitudes and the eastern reaches of Europe. In this study we present the first comprehensive evolutionary history of the A. arenosa species complex, covering its natural range, by using chromosome counts, nuclear AFLP data, and a maternally inherited marker from the chloroplast genome [trnL intron (trnL) and trnL/F intergenic spacer (trnL/F-IGS) of tRNA(Leu) and tRNA(Phe), respectively]. We unravel the broad-scale cytogeographic and phylogeographic patterns of diploids and tetraploids. Diploid cytotypes were exclusively found on the Balkan Peninsula and in the Carpathians while tetraploid cytotypes were found throughout the remaining distribution range of the A. arenosa complex. Three centers of genetic diversity were identified: the Balkan Peninsula, the Carpathians, and the unglaciated Eastern and Southeastern Alps. All three could have served as long-term refugia during Pleistocene climate oscillations. We hypothesize that the Western Carpathians were and still are the cradle of speciation within the A. arenosa complex due to the high species number and genetic diversity and the concurrence of both cytotypes there.