Arabidopsis hybrid speciation processes

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.

Rapid large-scale genomic introgression in Arabidopsis suecica via an autoallohexaploid bridge

Gene flow between species in the genus Arabidopsis occurs in significant amounts, but how exactly gene flow is achieved is not well understood. Polyploidization may be one avenue to explain gene flow between species. One problem, however, with polyploidization as a satisfying explanation is the occurrence of lethal genomic instabilities in neopolyploids as a result of genomic exchange, erratic meiotic behavior, and genomic shock. We have created an autoallohexaploid by pollinating naturally co-occurring diploid Arabidopsis thaliana with allotetraploid Arabidopsis suecica (an allotetraploid composed of A. thaliana and Arabidopsis arenosa). Its triploid offspring underwent spontaneous genome duplication and was used to generate a multigenerational pedigree. Using genome resequencing, we show that 2 major mechanisms promote stable genomic exchange in this population. Legitimate meiotic recombination and chromosome segregation between the autopolyploid chromosomes of the 2 A. thaliana genomes occur without any obvious bias for the parental origin and combine the A. thaliana haplotypes from the A. thaliana parent with the A. thaliana haplotypes from A. suecica similar to purely autopolyploid plants. In addition, we repeatedly observed that occasional exchanges between regions of the homoeologous chromosomes are tolerated. The combination of these mechanisms may result in gene flow leading to stable introgression in natural populations. Unlike the previously reported resynthesized neoallotetraploid A. suecica, this population of autoallohexaploids contains mostly vigorous, and genetically, cytotypically, and phenotypically variable individuals. We propose that naturally formed autoallohexaploid populations might serve as an intermediate bridge between diploid and polyploid species, which can facilitate gene flow rapidly and efficiently.

Inference of Polyploid Origin and Inheritance Mode from Population Genomic Data

Whole-genome duplications yield varied chromosomal pairing patterns, ranging from strictly bivalent to multivalent, resulting in disomic and polysomic inheritance modes. In the bivalent case, homeologous chromosomes form pairs, where in a multivalent pattern all copies are homologous and are therefore free to pair and recombine. As sufficient sequencing data is more readily available than high-quality cytological assessments of meiotic behavior or population genetic assessment of allelic segregation, especially for non-model organisms, bioinformatics approaches to infer origins and inheritance modes of polyploids using short-read sequencing data are attractive. Here we describe two such approaches, where the first is based on distributions of allelic read depth at heterozygous sites within an individual, as the expectations of such distributions are different for disomic and polysomic inheritance modes. The second approach is more laborious and based on a phylogenetic assessment of partially phased haplotypes of a polyploid in comparison to the closest diploid relatives. We discuss the sources of deviations from expected inheritance patterns, advantages and pitfalls of both methods, effects of mating types on the performance of the methods, and possible future developments.

Adaptive introgression: how polyploidy reshapes gene flow landscapes

Rare yet accumulating evidence in both plants and animals shows that whole genome duplication (WGD, leading to polyploidy) can break down reproductive barriers, facilitating gene flow between otherwise isolated species. Recent population genomic studies in wild, outcrossing Arabidopsis arenosa and Arabidopsis lyrata indicate that this WGD-potentiated gene flow can be adaptive and highly specific in response to particular environmental and intracellular challenges. The mechanistic basis of WGD-mediated easing of species barrier strength seems to primarily lie in the relative dosage of each parental genome in the endosperm. While generalisations about polyploids can be fraught, this evidence indicates that the breakdown of these barriers, combined with diploid to polyploid gene flow and gene flow between polyploids, allows some polyploids to act as adaptable 'allelic sponges', enjoying increased potential to respond to challenging environments.

Novelty and convergence in adaptation to whole genome duplication

Whole genome duplication (WGD) can promote adaptation but is disruptive to conserved processes, especially meiosis. Studies in Arabidopsis arenosa revealed a coordinated evolutionary response to WGD involving interacting proteins controlling meiotic crossovers, which are minimised in an autotetraploid (within-species polyploid) to avoid mis-segregation. Here we test whether this surprising flexibility of a conserved essential process, meiosis, is recapitulated in an independent WGD system, Cardamine amara, 17 million years diverged from A. arenosa. We assess meiotic stability and perform population-based scans for positive selection, contrasting the genomic response to WGD in C. amara with that of A. arenosa. We found in C. amara the strongest selection signals at genes with predicted functions thought important to adaptation to WGD: meiosis, chromosome remodelling, cell cycle, and ion transport. However, genomic responses to WGD in the two species differ: minimal ortholog-level convergence emerged, with none of the meiosis genes found in A. arenosa exhibiting strong signal in C. amara. This is consistent with our observations of lower meiotic stability and occasional clonal spreading in diploid C. amara, suggesting that nascent C. amara autotetraploid lineages were preadapted by their diploid lifestyle to survive while enduring reduced meiotic fidelity. However, in contrast to a lack of ortholog convergence, we see process-level and network convergence in DNA management, chromosome organisation, stress signalling, and ion homeostasis processes. This gives the first insight into the salient adaptations required to meet the challenges of a WGD state and shows that autopolyploids can utilize multiple evolutionary trajectories to adapt to WGD.

Low dispersal and ploidy differences in a grass maintain photosynthetic diversity despite gene flow and habitat overlap

Geographical isolation facilitates the emergence of distinct phenotypes within a single species, but reproductive barriers or selection are needed to maintain the polymorphism after secondary contact. Here, we explore the processes that maintain intraspecific variation of C₄ photosynthesis, a complex trait that results from the combined action of multiple genes. The grass Alloteropsis semialata includes C₄ and non-C₄ populations, which have coexisted as a polyploid series for more than 1 million years in the miombo woodlands of Africa. Using population genomics, we show that there is genome-wide divergence for the photosynthetic types, but the current geographical distribution does not reflect a simple habitat displacement scenario as the genetic clusters overlap, being occasionally mixed within a given habitat. Despite evidence of recurrent introgression between non-C₄ and C₄ groups, in both diploids and polyploids, the distinct genetic lineages retain their identity, potentially because of selection against hybrids. Coupled with strong isolation by distance within each genetic group, this selection created a geographical mosaic of photosynthetic types. Diploid C₄ and non-C₄ types never grew together, and the C₄ type from mixed populations constantly belonged to the hexaploid lineage. By limiting reproductive interactions between photosynthetic types, the ploidy difference probably allows their co-occurrence, reinforcing the functional diversity within this species. Together, these factors enabled the persistence of divergent physiological traits of ecological importance within a single species despite gene flow and habitat overlap.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

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.

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.

Transgenerational effects of inter-ploidy cross direction on reproduction and F2 seed development of Arabidopsis thaliana F1 hybrid triploids

KEY MESSAGE

Reproduction in triploid plants is important for understanding polyploid population dynamics. We show that genetically identical reciprocal F1 hybrid triploids can display transgenerational epigenetic effects on viable F2 seed development. The success or failure of reproductive outcomes from intra-species crosses between plants of different ploidy levels is an important factor in flowering plant evolution and crop breeding. However, the effects of inter-ploidy cross directions on F1 hybrid offspring fitness are poorly understood. In Arabidopsis thaliana, hybridization between diploid and tetraploid plants can produce viable F1 triploid plants. When selfed, such F1 triploid plants act as aneuploid gamete production "machines" where the vast majority of gametes generated are aneuploid which, following sexual reproduction, can generate aneuploid swarms of F2 progeny (Henry et al. 2009). There is potential for some aneuploids to cause gametophyte abortion and/or F2 seed abortion (Henry et al. 2009). In this study, we analyse the reproductive success of 178 self-fertilized inter-accession F1 hybrid triploids and demonstrate that the proportions of aborted or normally developed F2 seeds from the selfed F1 triploids depend upon a combination of natural variation and cross direction, with strong interaction between these factors. Single-seed ploidy analysis indicates that the embryonic DNA content of phenotypically normal F2 seeds is highly variable and that these DNA content distributions are also affected by genotype and cross direction. Notably, genetically identical reciprocal F1 hybrid triploids display grandparent-of-origin effects on F2 seed set, and hence on the ability to tolerate aneuploidy in F2 seed. There are differences between reciprocal F1 hybrid triploids regarding the proportions of normal and aborted F2 seeds generated, and also for the DNA content averages and distributions of the F2 seeds. To identify genetic variation for tolerance of aneuploidy in F2 seeds, we carried out a GWAS which identified two SNPs, termed MOT and POT, which represent candidate loci for genetic control of the proportion of normal F2 seeds obtained from selfed F1 triploids. Parental and grandparental effects on F2 seeds obtained from selfed F1 triploids can have transgenerational consequences for asymmetric gene flow, emergence of novel genotypes in polyploid populations, and for control of F2 seed set in triploid crops.

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.

Diverse paths to hybrid incompatibility in Arabidopsis

One of the most essential questions of biology is to understand how different species have evolved. Hybrid incompatibility, a phenomenon in which hybrids show reduced fitness in comparison with their parents, can result in reproductive isolation and speciation. Therefore, studying hybrid incompatibility provides an entry point in understanding speciation. Hybrid incompatibilities are known throughout taxa, and the underlying mechanisms have mystified scientists since the theory of evolution by means of natural selection was introduced. In plants, it is only in recent years that the high-throughput genetic and molecular tools have become available for the Arabidopsis genus, thus helping to shed light on the different genes and molecular and evolutionary mechanisms that underlie hybrid incompatibilities. In this review, we highlight the current knowledge of diverse mechanisms that are known to contribute to hybrid incompatibility.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

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.

Transposon-derived small RNAs triggered by miR845 mediate genome dosage response in Arabidopsis

Chromosome dosage has substantial effects on reproductive isolation and speciation in both plants and animals, but the underlying mechanisms are largely obscure 1 . Transposable elements in animals can regulate hybridity through maternal small RNA 2 , whereas small RNAs in plants have been postulated to regulate dosage response via neighboring imprinted genes3,4. Here we show that a highly conserved microRNA in plants, miR845, targets the tRNAMet primer-binding site (PBS) of long terminal repeat (LTR) retrotransposons in Arabidopsis pollen, and triggers the accumulation of 21-22-nucleotide (nt) small RNAs in a dose-dependent fashion via RNA polymerase IV. We show that these epigenetically activated small interfering RNAs (easiRNAs) mediate hybridization barriers between diploid seed parents and tetraploid pollen parents (the 'triploid block'), and that natural variation for miR845 may account for 'endosperm balance' allowing the formation of triploid seeds. Targeting of the PBS with small RNA is a common mechanism for transposon control in mammals and plants, and provides a uniquely sensitive means to monitor chromosome dosage and imprinting in the developing seed.

Cryptic gene pools in the Hypericum perforatum-H. maculatum complex: diploid persistence versus trapped polyploid melting

BACKGROUND AND AIMS

In Central Europe Hypericum perforatum and Hypericum maculatum show significant hybridization and introgression as a consequence of Pleistocene range fluctuations, and their gene pools are merging on higher ploidy levels. This paper discusses whether polyploid hybrid gene pools are trapped in the ecological climatic niche space of their diploid ancestors, and tests the idea of geographical parthenogenesis.

METHODS

DNA sequence information of nuclear ribosomal DNA and plastid loci, ploidy level estimates and ecological niche modelling are used to characterize the various diploid and polyploid gene pools and unravel spatio-temporal patterns of gene flow among them.

KEY RESULTS

On the diploid level, the three gene pools are clearly distinct between and within species of H. perforatum (two gene pools) and H. maculatum, and their divergence dates back to the first half of the Pleistocene. All polyploids in Central Europe show high levels of past and contemporary gene flow between all three gene pools. The correlation of genetic and geographical distances breaks down if the latter is larger than 250 km, indicating recent and ongoing gene flow. The two species are ecologically differentiated, but in particular hybrids among all three gene pools do not show significant niche differences compared to their parental gene pools, except for some combinations with H. maculatum.

CONCLUSIONS

Inter- and intraspecific gene flow between inter- and intra-species gene pools is limited on the diploid level, and the geographical distribution of the diploids largely reflects Pleistocene evolutionary history. Secondary contact promoted hybridization and introgression on the polyploid level, enabling offspring to escape the diploid gene pools. However, the hybrid polyploids do not show significant niche differences compared to their diploid progenitors. It is concluded that the observed absence of niche divergence has precluded further differentiation and geographical partitioning of new polyploid lineages being effectively separated from the parental lines. The predominantly apomictic reproducing polyploids are trapped in the polyploid gene pool and the ecological climatic niche space of their diploid ancestors.

Hybrids and horizontal transfer: introgression allows adaptive allele discovery

Evolution has devised countless remarkable solutions to diverse challenges. Understanding the mechanistic basis of these solutions provides insights into how biological systems can be subtly tweaked without maladaptive consequences. The knowledge gained from illuminating these mechanisms is equally important to our understanding of fundamental evolutionary mechanisms as it is to our hopes of developing truly rational plant breeding and synthetic biology. In particular, modern population genomic approaches are proving very powerful in the detection of candidate alleles for mediating consequential adaptations that can be tested functionally. Especially striking are signals gained from contexts involving genetic transfers between populations, closely related species, or indeed between kingdoms. Here we discuss two major classes of these scenarios, adaptive introgression and horizontal gene flow, illustrating discoveries made across kingdoms.

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.

Altered chromatin compaction and histone methylation drive non-additive gene expression in an interspecific Arabidopsis hybrid

BACKGROUND

The merging of two diverged genomes can result in hybrid offspring that phenotypically differ greatly from both parents. In plants, interspecific hybridization plays important roles in evolution and speciation. In addition, many agricultural and horticultural species are derived from interspecific hybridization. However, the detailed mechanisms responsible for non-additive phenotypic novelty in hybrids remain elusive.

RESULTS

In an interspecific hybrid between Arabidopsis thaliana and A. lyrata, the vast majority of genes that become upregulated or downregulated relative to the parents originate from A. thaliana. Among all differentially expressed A. thaliana genes, the majority is downregulated in the hybrid. To understand why parental origin affects gene expression in this system, we compare chromatin packing patterns and epigenomic landscapes in the hybrid and parents. We find that the chromatin of A. thaliana, but not that of A. lyrata, becomes more compact in the hybrid. Parental patterns of DNA methylation and H3K27me3 deposition are mostly unaltered in the hybrid, with the exception of higher CHH DNA methylation in transposon-rich regions. However, A. thaliana genes enriched for the H3K27me3 mark are particularly likely to differ in expression between the hybrid and parent.

CONCLUSIONS

It has long been suspected that genome-scale properties cause the differential responses of genes from one or the other parent to hybridization. Our work links global chromatin compactness and H3K27me3 histone modification to global differences in gene expression in an interspecific Arabidopsis hybrid.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Chloroplast genomes of Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea: Structures and comparative analysis

We investigated the complete chloroplast (cp) genomes of non-model Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea using Illumina paired-end sequencing to understand their genetic organization and structure. Detailed bioinformatics analysis revealed genome sizes of both subspecies ranging between 154.4~154.5 kbp, with a large single-copy region (84,197~84,158 bp), a small single-copy region (17,738~17,813 bp) and pair of inverted repeats (IRa/IRb; 26,264~26,259 bp). Both cp genomes encode 130 genes, including 85 protein-coding genes, eight ribosomal RNA genes and 37 transfer RNA genes. Whole cp genome comparison of A. halleri ssp. gemmifera and A. lyrata ssp. petraea, along with ten other Arabidopsis species, showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. The location and distribution of repeat sequences were determined, and sequence divergences of shared genes were calculated among related species. Comparative phylogenetic analysis of the entire genomic data set and 70 shared genes between both cp genomes confirmed the previous phylogeny and generated phylogenetic trees with the same topologies. The sister species of A. halleri ssp. gemmifera is A. umezawana, whereas the closest relative of A. lyrata spp. petraea is A. arenicola.

Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis Species

The generation and maintenance of functional variation in the pathogen defense system of plants is central to the constant evolutionary battle between hosts and parasites. If a species is susceptible to a given pathogen, hybridization and subsequent introgression of a resistance allele from a related species can potentially be an important source of new immunity and is therefore expected to be selected for in a process referred to as adaptive introgression. Here, we survey sequence variation in 10 resistance (R-) genes and compare them with 37 reference genes in natural populations of the two closely related and interfertile species: Arabidopsis lyrata and A. halleri. The R-genes are highly polymorphic in both species and show clear signs of trans-species polymorphisms. We show that A. lyrata and A. halleri have had a history of limited introgression for the reference genes. For the R-genes, the introgression rate has been significantly higher than for the reference genes, resulting in fewer fixed differences between species and a higher sharing of identical haplotypes. We conclude that R-genes likely cross the species boundaries at a higher rate than reference genes and therefore also that some of the increased diversity and trans-specific polymorphisms in R-genes is due to adaptive introgression.

Genomic studies of adaptive evolution in outcrossing Arabidopsis species

Large-scale population genomic approaches have very recently been fruitfully applied to the Arabidopsis relatives Arabidopsis halleri, A. lyrata and especially A. arenosa. In contrast to A. thaliana, these species are obligately outcrossing and thus the footprints of natural selection are more straightforward to detect. Furthermore, both theoretical and empirical studies indicate that outcrossers are better able to evolve in response to selection pressure. As a result, recent work in these species serves as a paradigm of population genomic studies of adaptation both to environmental as well as intracellular challenges.

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.

Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism

The notion of species as reproductively isolated units related through a bifurcating tree implies that gene trees should generally agree with the species tree and that sister taxa should not share polymorphisms unless they diverged recently and should be equally closely related to outgroups. It is now possible to evaluate this model systematically. We sequenced multiple individuals from 27 described taxa representing the entire Arabidopsis genus. Cluster analysis identified seven groups, corresponding to described species that capture the structure of the genus. However, at the level of gene trees, only the separation of Arabidopsis thaliana from the remaining species was universally supported, and, overall, the amount of shared polymorphism demonstrated that reproductive isolation was considerably more recent than the estimated divergence times. We uncovered multiple cases of past gene flow that contradict a bifurcating species tree. Finally, we showed that the pattern of divergence differs between gene ontologies, suggesting a role for selection.

Borrowed alleles and convergence in serpentine adaptation

Serpentine barrens represent extreme hazards for plant colonists. These sites are characterized by high porosity leading to drought, lack of essential mineral nutrients, and phytotoxic levels of metals. Nevertheless, nature forged populations adapted to these challenges. Here, we use a population-based evolutionary genomic approach coupled with elemental profiling to assess how autotetraploid Arabidopsis arenosa adapted to a multichallenge serpentine habitat in the Austrian Alps. We first demonstrate that serpentine-adapted plants exhibit dramatically altered elemental accumulation levels in common conditions, and then resequence 24 autotetraploid individuals from three populations to perform a genome scan. We find evidence for highly localized selective sweeps that point to a polygenic, multitrait basis for serpentine adaptation. Comparing our results to a previous study of independent serpentine colonizations in the closely related diploid Arabidopsis lyrata in the United Kingdom and United States, we find the highest levels of differentiation in 11 of the same loci, providing candidate alleles for mediating convergent evolution. This overlap between independent colonizations in different species suggests that a limited number of evolutionary strategies are suited to overcome the multiple challenges of serpentine adaptation. Interestingly, we detect footprints of selection in A. arenosa in the context of substantial gene flow from nearby off-serpentine populations of A. arenosa, as well as from A. lyrata In several cases, quantitative tests of introgression indicate that some alleles exhibiting strong selective sweep signatures appear to have been introgressed from A. lyrata This finding suggests that migrant alleles may have facilitated adaptation of A. arenosa to this multihazard environment.

Environmental Regulation of Heterosis in the Allopolyploid Arabidopsis suecica

Allopolyploids are organisms possessing more than two complete sets of chromosomes from two or more species and are frequently more vigorous than their progenitors. To address the question why allopolyploids display hybrid vigor, we compared the natural allopolyploid Arabidopsis suecica to its progenitor species Arabidopsis thaliana and Arabidopsis arenosa. We measured chlorophyll content, CO2 assimilation, and carbohydrate production under varying light conditions and found that the allopolyploid assimilates more CO2 per unit chlorophyll than either of the two progenitor species in high intensity light. The increased carbon assimilation corresponds with greater starch accumulation, but only in strong light, suggesting that the strength of hybrid vigor is dependent on environmental conditions. In weaker light A. suecica tends to produce as much primary metabolites as the better progenitor. We found that gene expression of LIMIT DEXTRINASE1, a debranching enzyme that cleaves branch points within starch molecules, is at the same level in the allopolyploid as in the maternal progenitor A. thaliana and significantly more expressed than in the paternal progenitor A. arenosa. However, expression differences of β-amylases and GLUCAN-WATER DIKINASE1 were not statistically significantly elevated in the allopolyploid over progenitor expression levels. In contrast to allopolyploids, autopolyploid A. thaliana showed the same photosynthetic rate as diploids, indicating that polyploidization alone is likely not the reason for enhanced vigor in the allopolyploid. Taken together, our data suggest that the magnitude of heterosis in A. suecica is environmentally regulated, arises from more efficient photosynthesis, and, under specific conditions, leads to greater starch accumulation than in its progenitor species.

A subunit of the oligosaccharyltransferase complex is required for interspecific gametophyte recognition in Arabidopsis

Species-specific gamete recognition is a key premise to ensure reproductive success and the maintenance of species boundaries. During plant pollen tube (PT) reception, gametophyte interactions likely allow the species-specific recognition of signals from the PT (male gametophyte) by the embryo sac (female gametophyte), resulting in PT rupture, sperm release, and double fertilization. This process is impaired in interspecific crosses between Arabidopsis thaliana and related species, leading to PT overgrowth and a failure to deliver the sperm cells. Here we show that ARTUMES (ARU) specifically regulates the recognition of interspecific PTs in A. thaliana. ARU, identified in a genome-wide association study (GWAS), exclusively influences interspecific--but not intraspecific--gametophyte interactions. ARU encodes the OST3/6 subunit of the oligosaccharyltransferase complex conferring protein N-glycosylation. Our results suggest that glycosylation patterns of cell surface proteins may represent an important mechanism of gametophyte recognition and thus speciation.

Exogenous selection rather than cytonuclear incompatibilities shapes asymmetrical fitness of reciprocal Arabidopsis hybrids

Reciprocal crosses between species often display an asymmetry in the fitness of F1 hybrids. This pattern, referred to as isolation asymmetry or Darwin's corollary to Haldane's rule, is a general feature of reproductive isolation in plants, yet factors determining its magnitude and direction remain unclear. We evaluated reciprocal species crosses between two naturally hybridizing diploid species of Arabidopsis to assess the degree of isolation asymmetry at different postmating life stages. We found that pollen from Arabidopsis arenosa will usually fertilize ovules from Arabidopsis lyrata; the reverse receptivity being less complete. Maternal A. lyrata parents set more F1 hybrid seed, but germinate at lower frequency, reversing the asymmetry. As predicted by theory, A. lyrata (the maternal parent with lower seed viability in crosses) exhibited accelerated chloroplast evolution, indicating that cytonuclear incompatibilities may play a role in reproductive isolation. However, this direction of asymmetrical reproductive isolation is not replicated in natural suture zones, where delayed hybrid breakdown of fertility at later developmental stages, or later-acting selection against A. arenosa maternal hybrids (unrelated to hybrid fertility, e.g., substrate adaptation) may be responsible for an excess of A. lyrata maternal hybrids. Exogenous selection rather than cytonuclear incompatibilities thus shapes the asymmetrical postmating isolation in nature.

Cold tolerance in Arabidopsis kamchatica

•

PREMISE OF THE STUDY

Cold tolerance is a critically important factor determining how plants will be influenced by climate change, including changes in snowcover and extreme weather events. Although a great deal is known about cold tolerance in Arabidopsis thaliana, it is not highly cold tolerant. This study examined cold tolerance and its genetic diversity in an herbaceous subarctic relative, Arabidopsis kamchatica, which generally occurs in much colder climates.•

METHODS

Thermal analysis and electrolyte leakage were used to estimate supercooling points and lethal temperatures (LT50) in cold-acclimated and nonacclimated families from three populations of A. kamchatica.•

KEY RESULTS

Arabidopsis kamchatica was highly cold tolerant, with a mean LT50 of -10.8°C when actively growing, and -21.8°C when cold acclimated. It also was able to supercool to very low temperatures. Surprisingly, actively growing plants supercooled more than acclimated plants (-14.7 vs. -12.7°C). There was significant genetic variation for cold tolerance both within and among populations. However, both cold tolerance and genetic diversity were highest in the midlatitude population rather than in the far north, indicating that adaptations to climate change are most likely to arise in the center of the species range rather than at the edges.•

CONCLUSIONS

Arabidopsis kamchatica is highly cold tolerant throughout its range. It is far more freeze tolerant than A. thaliana, and supercooled to lower temperatures, suggesting that A. kamchatica provides a valuable complement to A. thaliana for cold tolerance research.

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.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae)

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

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.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

The BOY NAMED SUE quantitative trait locus confers increased meiotic stability to an adapted natural allopolyploid of Arabidopsis

Whole-genome duplication resulting from polyploidy is ubiquitous in the evolutionary history of plant species. Yet, polyploids must overcome the meiotic challenge of pairing, recombining, and segregating more than two sets of chromosomes. Using genomic sequencing of synthetic and natural allopolyploids of Arabidopsis thaliana and Arabidopsis arenosa, we determined that dosage variation and chromosomal translocations consistent with homoeologous pairing were more frequent in the synthetic allopolyploids. To test the role of structural chromosomal differentiation versus genetic regulation of meiotic pairing, we performed sequenced-based, high-density genetic mapping in F2 hybrids between synthetic and natural lines. This F2 population displayed frequent dosage variation and deleterious homoeologous recombination. The genetic map derived from this population provided no indication of structural evolution of the genome of the natural allopolyploid Arabidopsis suecica, compared with its predicted parents. The F2 population displayed variation in meiotic regularity and pollen viability that correlated with a single quantitative trait locus, which we named BOY NAMED SUE, and whose beneficial allele was contributed by A. suecica. This demonstrates that an additive, gain-of-function allele contributes to meiotic stability and fertility in a recently established allopolyploid and provides an Arabidopsis system to decipher evolutionary and molecular mechanisms of meiotic regularity in polyploids.

Taxonomy and systematics are key to biological information: Arabidopsis, Eutrema (Thellungiella), Noccaea and Schrenkiella (Brassicaceae) as examples

Taxonomy and systematics provide the names and evolutionary framework for any biological study. Without these names there is no access to a biological context of the evolutionary processes which gave rise to a given taxon: close relatives and sister species (hybridization), more distantly related taxa (ancestral states), for example. This is not only true for the single species a research project is focusing on, but also for its relatives, which might be selected for comparative approaches and future research. Nevertheless, taxonomical and systematic knowledge is rarely fully explored and considered across biological disciplines. One would expect the situation to be more developed with model organisms such as Noccaea, Arabidopsis, Schrenkiella and Eutrema (Thellungiella). However, we show the reverse. Using Arabidopsis halleri and Noccaea caerulescens, two model species among metal accumulating taxa, we summarize and reflect past taxonomy and systematics of Arabidopsis and Noccaea and provide a modern synthesis of taxonomic, systematic and evolutionary perspectives. The same is presented for several species of Eutrema s. l. and Schrenkiella recently appeared as models for studying stress tolerance in plants and widely known under the name Thellungiella.

A population genetic model to infer allotetraploid speciation and long-term evolution applied to two yarrow species

Allotetraploid speciation, that is, the generation of a hybrid tetraploid species from two diploid species, and the long-term evolution of tetraploid populations and species are important in plants. We developed a population genetic model to infer population genetic parameters of tetraploid populations from data of the progenitor and descendant species. Two yarrow species, Achillea alpina-4x and A. wilsoniana-4x, arose by allotetraploidization from the diploid progenitors, A. acuminata-2x and A. asiatica-2x. Yet, the population genetic process has not been studied in detail. We applied the model to sequences of three nuclear genes in populations of the four yarrow species and compared their pattern of variability with that in four plastid regions. The plastid data indicated that the two tetraploid species probably originated from multiple independent allopolyploidization events and have accumulated many mutations since. With the nuclear data, we found a low rate of homeologous recombination or gene conversion and a reduction in diversity relative to the level of both diploid species combined. The present analysis with a novel probabilistic model suggests a genetic bottleneck during tetraploid speciation, that the two tetraploid species have a long evolutionary history, and that they have a small amount of genetic exchange between the homeologous genomes.

Evolution of cryptic gene pools in Hypericum perforatum: the influence of reproductive system and gene flow

BACKGROUND AND AIMS

Hypericum perforatum (St. John's wort) is a widespread Eurasian perennial plant species with remarkable variation in its morphology, ploidy and breeding system, which ranges from sex to apomixis. Here, hypotheses on the evolutionary origin of St. John's wort are tested and contrasted with the subsequent history of interspecific gene flow.

METHODS

Extensive field collections were analysed for quantitative morphological variation, ploidy, chromosome numbers and genetic diversity using nuclear (amplified fragment length polymorphism) and plastid (trnL-trnF) markers. The mode of reproduction was analysed by FCSS (flow cytometric seed screen).

KEY RESULTS

It is demonstrated that H. perforatum is not of hybrid origin, and for the first time wild diploid populations are documented. Pseudogamous facultative apomictic reproduction is prevalent in the polyploids, whereas diploids are predominantly sexual, a phenomenon which also characterizes its sister species H. maculatum. Both molecular markers characterize identical major gene pools, distinguishing H. perforatum from H. maculatum and two genetic groups in H. perforatum. All three gene pools are in close geographical contact. Extensive gene flow and hybridization throughout Europe within and between gene pools and species is exemplified by the molecular data and confirmed by morphometric analyses.

CONCLUSIONS

Hypericum perforatum is of a single evolutionary origin and later split into two major gene pools. Subsequently, independent and recurrent polyploidization occurred in all lineages and was accompanied by substantial gene flow within and between H. perforatum and H. maculatum. These processes are highly influenced by the reproductive system in both species, with a switch to predominantly apomictic reproduction in polyploids, irrespective of their origin.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries

PREMISE OF THE STUDY

Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. •

METHODS

We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. •

KEY RESULTS

Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. •

CONCLUSIONS

This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Genetic changes in flowering and morphology in response to adaptation to a high-latitude environment in Arabidopsis lyrata

BACKGROUND AND AIMS

The adaptive plastic reactions of plant populations to changing climatic factors, such as winter temperatures and photoperiod, have changed during range shifts after the last glaciation. Timing of flowering is an adaptive trait regulated by environmental cues. Its genetics has been intensively studied in annual plants, but in perennials it is currently not well characterized. This study examined the genetic basis of differentiation in flowering time, morphology, and their plastic responses to vernalization in two locally adapted populations of the perennial Arabidopsis lyrata: (1) to determine whether the two populations differ in their vernalization responses for flowering phenology and morphology; and (2) to determine the genomic areas governing differentiation and vernalization responses.

METHODS

Two A. lyrata populations, from central Europe and Scandinavia, were grown in growth-chamber conditions with and without cold treatment. A QTL analysis was performed to find genomic regions that interact with vernalization.

KEY RESULTS

The population from central Europe flowered more rapidly and invested more in inflorescence growth than the population from alpine Scandinavia, especially after vernalization. The alpine population had consistently a low number of inflorescences and few flowers, suggesting strong constraints due to a short growing season, but instead had longer leaves and higher leaf rosettes. QTL mapping in the F2 population revealed genomic regions governing differentiation in flowering time and morphology and, in some cases, the allelic effects from the two populations on a trait were influenced by vernalization (QTL × vernalization interactions).

CONCLUSIONS

The results indicate that many potentially adaptive genetic changes have occurred during colonization; the two populations have diverged in their plastic responses to vernalization in traits closely connected to fitness through changes in many genomic areas.

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.

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.