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

The self-compatibility is acquired after polyploidization: a case study of Brassica napus self-incompatible trilinear hybrid breeding system

Self-incompatibility plays a vital role in angiosperms, by preventing inbreeding depression and maintaining genetic diversity within populations. Following polyploidization, many angiosperm species transition from self-incompatibility to self-compatibility. Here, we investigated the S-locus in Brassicaceae and identified distinct origins for the sRNA loci, SMI and SMI2 (SCR Methylation Inducer 1 and 2), within the S-locus. The SMI loci were found to be widespread in Cruciferae, whereas the SMI2 loci were exclusive to Brassica species. Additionally, we discovered four major S-haplotypes (BnS-1, BnS-6, BnS-7, and BnS-1300) in rapeseed. Overexpression of BnSMI-1 in self-incompatible Brassica napus ('S-70S1300S6 ') resulted in a significant increase in DNA methylation in the promoter regions of BnSCR-6 and BnSCR-1300, leading to self-compatibility. Conversely, by overexpressing a point mutation of BnSmi-1 in the 'S-70S1300S6 ' line, we observed lower levels of DNA methylation in BnSCR-6 and BnSCR-1300 promoters. Furthermore, the overexpression of BnSMI2-1300 in the 'SI-326S7S6 ' line inhibited the expression of BnSCR-7 through transcriptional repression of the Smi2 sRNA from the BnS-1300 haplotype. Our study demonstrates that the self-compatibility of rapeseed is determined by S-locus sRNA-mediated silencing of SCR after polyploidization, which helps to further breed self-incompatible or self-compatible rapeseed lines, thereby facilitating the utilization of heterosis.

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.].

Breakdown of self-incompatibility due to genetic interaction between a specific S-allele and an unlinked modifier

Breakdown of self-incompatibility has frequently been attributed to loss-of-function mutations of alleles at the locus responsible for recognition of self-pollen (i.e. the S-locus). However, other potential causes have rarely been tested. Here, we show that self-compatibility of S1S1-homozygotes in selfing populations of the otherwise self-incompatible Arabidopsis lyrata is not due to S-locus mutation. Between-breeding-system cross-progeny are self-compatible if they combine S1 from the self-compatible cross-partner with recessive S1 from the self-incompatible cross-partner, but self-incompatible with dominant S-alleles. Because S1S1 homozygotes in outcrossing populations are self-incompatible, mutation of S1 cannot explain self-compatibility in S1S1 cross-progeny. This supports the hypothesis that an S1-specific modifier unlinked to the S-locus causes self-compatibility by functionally disrupting S1. Self-compatibility in S19S19 homozygotes may also be caused by an S19-specific modifier, but we cannot rule out a loss-of-function mutation of S19. Taken together, our findings indicate that breakdown of self-incompatibility is possible without disruptive mutations at the S-locus.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

Uncovering the first complete plastome genomics, comparative analyses, and phylogenetic dispositions of endemic medicinal plant Ziziphus hajarensis (Rhamnaceae)

BACKGROUND

Ziziphus hajarensis is an endemic plant species well-distributed in the Western Hajar mountains of Oman. Despite its potential medicinal uses, little is known regarding its genomic architecture, phylogenetic position, or evolution. Here we sequenced and analyzed the entire chloroplast (cp) genome of Z. hajarensis to understand its genetic organization, structure, and phylogenomic disposition among Rhamnaceae species.

RESULTS

The results revealed the genome of Z. hajarensis cp comprised 162,162 bp and exhibited a typical quadripartite structure, with a large single copy (LSC) region of 895,67 bp, a small single copy (SSC) region of 19,597 bp and an inverted repeat (IR) regions of 26,499 bp. In addition, the cp genome of Z. hajarensis comprises 126 genes, including 82 protein-coding genes, eight rRNA genes, and 36 tRNA genes. Furthermore, the analysis revealed 208 microsatellites, 96.6% of which were mononucleotides. Similarly, a total of 140 repeats were identified, including 11 palindromic, 24 forward, 14 reverse, and 104 tandem repeats. The whole cp genome comparison of Z. hajarensis and nine other species from family Rhamnaceae showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. Comparative phylogenetic analysis based on the complete cp genome, 66 shared genes and matK gene revealed that Z. hajarensis shares a clade with Z. jujuba and that the family Rhamnaceae is the closest family to Barbeyaceae and Elaeagnaceae.

CONCLUSION

All the genome features such as genome size, GC content, genome organization and gene order were highly conserved compared to the other related genomes. The whole cp genome of Z. hajarensis gives fascinating insights and valuable data that may be used to identify related species and reconstruct the phylogeny of the species.

Cucumber Mosaic Virus Infection in Arabidopsis: A Conditional Mutualistic Symbiont?

A cucumber mosaic virus isolate, named Ho [CMV(Ho)], was isolated from a symptomless Arabidopsis halleri field sample containing low virus titers. An analysis of CMV(Ho) RNA molecules indicated that the virus isolate, besides the usual cucumovirus tripartite RNA genome, additionally contained defective RNA3 molecules and a satellite RNA. To study the underlying mechanism of the persistent CMV(Ho) infection in perennial A. halleri, infectious cDNA clones were generated for all its genetic elements. CMV, which consists of synthetic transcripts from the infectious tripartite RNA genomes, and designated CMV(Ho)tr, multiplied in A. halleri and annual Arabidopsis thaliana Col-0 to a similar level as the virulent strain CMV(Y), but did not induce any symptoms in them. The response of Col-0 to a series of reassortant CMVs between CMV(Ho)tr and CMV(Y) suggested that the establishment of an asymptomatic phenotype of CMV(Ho) infection was due to the 2b gene of CMV RNA2, but not due to the presence of the defective RNA3 and satellite RNA. The accumulation of CMV(Ho) 2b protein tagged with the FLAG epitope (2b.Ho-FLAG) in 2b.Ho-FLAG-transformed Col-0 did not induce any symptoms, suggesting a 2b-dependent persistency of CMV(Ho)tr infection in Arabidopsis. The 2b protein interacted with Argonaute 4, which is known to regulate the cytosine methylation levels of host genomic DNA. Whole genomic bisulfite sequencing analysis of CMV(Ho)tr- and mock-inoculated Col-0 revealed that cytosine hypomethylation in the promoter regions of 82 genes, including two genes encoding transcriptional regulators (DOF1.7 and CBP1), was induced in response to CMV(Ho)tr infection. Moreover, the increased levels of hypomethylation in the promoter region of both genes, during CMV(Ho)tr infection, were correlated with the up- or down-regulation of their expression. Taken altogether, the results indicate that during persistent CMV(Ho) infection in Arabidopsis, host gene expression may be epigenetically modulated resulting from a 2b-mediated cytosine hypomethylation of host genomic DNA.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Flow Cytometry and Sorting in Arabidopsis

Flow cytometry and sorting represents a valuable and mature experimental platform for the analysis of cellular populations. Applications involving higher plants started to emerge around 40 years ago and are now widely employed both to provide unique information regarding basic and applied questions in the biosciences and to advance agricultural productivity in practical ways. Further development of this platform is being actively pursued, and this promises additional progress in our understanding of the interactions of cells within complex tissues and organs. Higher plants offer unique challenges in terms of flow cytometric analysis, first since their organs and tissues are, almost without exception, three-dimensional assemblies of different cell types held together by tough cell walls, and, second, because individual plant cells are generally larger than those of mammals.This chapter, which updates work last reviewed in 2014 [Galbraith DW (2014) Flow cytometry and sorting in Arabidopsis. In: Sanchez Serrano JJ, Salinas J (eds) Arabidopsis Protocols, 3rd ed. Methods in molecular biology, vol 1062. Humana Press, Totowa, pp 509-537], describes the application of techniques of flow cytometry and sorting to the model plant species Arabidopsis thaliana, in particular emphasizing (a) fluorescence labeling in vivo of specific cell types and of subcellular components, (b) analysis using both conventional cytometers and spectral analyzers, (c) fluorescence-activated sorting of protoplasts and nuclei, and (d) transcriptome analyses using sorted protoplasts and nuclei, focusing on population analyses at the level of single protoplasts and nuclei. Since this is an update, details of new experimental methods are emphasized.

Unraveling the Chloroplast Genomes of Two Prosopis Species to Identify Its Genomic Information, Comparative Analyses and Phylogenetic Relationship

Genus Prosopis (family Fabaceae) are shrubby trees, native to arid and semi-arid regions of Asia, Africa, and America and known for nitrogen fixation. Here, we have sequenced the complete chloroplast (cp) genomes of two Prosopis species (P. juliflora and P. cineraria) and compared them with previously sequenced P. glandulosa, Adenanthera microsperma, and Parkia javanica belonging to the same family. The complete genome sequences of Prosopis species and related species ranged from 159,389 bp (A. microsperma) to 163,677 bp (P. cineraria). The overall GC contents of the genomes were almost the similar (35.9–36.6%). The P. juliflora and P. cineraria genomes encoded 132 and 131 genes, respectively, whereas both the species comprised of 85 protein-coding genes higher than other compared species. About 140, 134, and 129 repeats were identified in P. juliflora, P. cineraria and P. glandulosa cp genomes, respectively. Similarly, the maximum number of simple sequence repeats were determined in P. juliflora (88), P. cineraria (84), and P. glandulosa (78). Moreover, complete cp genome comparison determined a high degree of sequence similarity among P. juliflora, P. cineraria, and P. glandulosa, however some divergence in the intergenic spacers of A. microsperma and Parkia javanica were observed. The phylogenetic analysis showed that P. juliflora is closer to P. cineraria than P. glandulosa.

Expanded inverted repeat region with large scale inversion in the first complete plastid genome sequence of Plantago ovata

Plantago ovata (Plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from Plantaginaceae. The results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. The divergence time estimation as employed in BEAST revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other Plantaginaceae species.

Zahora, a new monotypic genus from tribe Brassiceae (Brassicaceae) endemic to the Moroccan Sahara

Zahora ait-atta Lemmel & M.Koch, a new species from the Moroccan Sahara, is described and documented here and constitutes a monotypic new genus. The new taxon belongs to the tribe Brassiceae (Brassicaceae), and cytogenetic and phylogenetic analyses reveal that this diploid species has a remote status of Miocene origin in the northwestern Sahara Desert. We examined the morphological differences between morphologically related genera and provide photographs of the new species. The new genus may play a key role in future Brassica-Raphanus crop research since it is placed phylogenetically at the base of a generically highly diverse clade including Raphanus sativus, and it shows affinities to various Brassica species.

Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation

Local adaptation is assumed to occur when populations differ in a phenotypic trait or a set of traits, and such variation has a genetic basis. Here, we introduce Arabidopsis halleri and its life history as a perennial model system to study population differentiation and local adaptation. Studies on altitudinal adaptation have been conducted in two regions: Mt. Ibuki in Japan and the European Alps. Several studies have demonstrated altitudinal adaptation in ultraviolet-B (UV-B) tolerance, leaf water repellency against spring frost and anti-herbivore defences. Studies on population differentiation in A. halleri have also focused on metal hyperaccumulation and tolerance to heavy metal contamination. In these study systems, genome scans to identify candidate genes under selection have been applied. Lastly, we briefly discuss how RNA-Seq can broaden phenotypic space and serve as a link to underlying mechanisms. In conclusion, A. halleri provides us with opportunities to study population differentiation and local adaptation, and relate these to the genetic systems underlying target functional traits.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.

Genetic architecture of a plant adaptive trait: QTL mapping of intraspecific variation for tolerance to metal pollution in Arabidopsis halleri

Anthropogenic activities are among the main drivers of global change and result in drastic habitat modifications, which represent strong evolutionary challenges for biological species that can either migrate, adapt, or disappear. In this context, understanding the genetics of adaptive traits is a prerequisite to enable long-term maintenance of populations under strong environmental constraints. To examine these processes, a QTL approach was developed here using the pseudometallophyte Arabidopsis halleri, which displays among-population adaptive divergence for tolerance to metallic pollution in soils. An F2 progeny was obtained by crossing individuals from metallicolous and non-metallicolous populations from Italian Alps, where intense metallurgic activities have created strong landscape heterogeneity. Then, we combined genome de novo assembly and genome resequencing of parental genotypes to obtain single-nucleotide polymorphism markers and achieve high-throughput genotyping of the progeny. QTL analysis was performed using growth parameters and photosynthetic yield to assess zinc tolerance levels. One major QTL was identified for photosynthetic yield. It explained about 27% of the phenotypic variance. Functional annotation of the QTL and gene expression analyses highlighted putative candidate genes. Our study represents a successful approach combining evolutionary genetics and advanced molecular tools, helping to better understand how a species can face new selective pressures of anthropogenic origin.

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.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4806-7) contains supplementary material, which is available to authorized users.

Evolution of interspecies unilateral incompatibility in the relatives of Arabidopsis thaliana

The evolutionary concurrence of intraspecies self-incompatibility (SI) and explosive angiosperm radiation in the Cretaceous have led to the hypothesis that SI was one of the predominant drivers of rapid speciation in angiosperms. Interspecies unilateral incompatibility (UI) usually occurs when pollen from a self-compatible (SC) species is rejected by the pistils of a SI species, while the reciprocal pollination is compatible (UC). Although this SI × SC type UI is most prevalent and viewed as a prezygotic isolation barrier to promote incipient speciation of angiosperms, comparative evidence to support such a role is lacking. We show that SI × SI type UI in SI species pairs is also common in the well-characterized accessions representing the four major lineages of the Arabidopsis genus and is developmentally regulated. This allowed us to reveal a strong correlation between UI strength and species divergence in these representative accessions. In addition, analyses of a SC accession and the pseudo-self-compatible (psc) spontaneous mutant of Arabidopsis lyrata indicate that UI shares, at least, common pollen rejection pathway with SI. Furthermore, genetic and genomic analyses of SI × SI type UI in A. lyrata × A. arenosa species pair showed that two major-effect quantitative trait loci are the stigma and pollen-side determinant of UI, respectively, which could be involved in heterospecies pollen discrimination. By revealing a close link between UI and SI pathway, particularly between UI and species divergence in these representative accessions, our findings establish a connection between SI and speciation. Thus, the pre-existence of SI system would have facilitated the evolution of UI and accordingly promote speciation.

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.

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

Divergent biology of facultative heavy metal plants

Among heavy metal plants (the metallophytes), facultative species can live both in soils contaminated by an excess of heavy metals and in non-affected sites. In contrast, obligate metallophytes are restricted to polluted areas. Metallophytes offer a fascinating biology, due to the fact that species have developed different strategies to cope with the adverse conditions of heavy metal soils. The literature distinguishes between hyperaccumulating, accumulating, tolerant and excluding metallophytes, but the borderline between these categories is blurred. Due to the fact that heavy metal soils are dry, nutrient limited and are not uniform but have a patchy distribution in many instances, drought-tolerant or low nutrient demanding species are often regarded as metallophytes in the literature. In only a few cases, the concentrations of heavy metals in soils are so toxic that only a few specifically adapted plants, the genuine metallophytes, can cope with these adverse soil conditions. Current molecular biological studies focus on the genetically amenable and hyperaccumulating Arabidopsis halleri and Noccaea (Thlaspi) caerulescens of the Brassicaceae. Armeria maritima ssp. halleri utilizes glands for the excretion of heavy metals and is, therefore, a heavy metal excluder. The two endemic zinc violets of Western Europe, Viola lutea ssp. calaminaria of the Aachen-Liège area and Viola lutea ssp. westfalica of the Pb-Cu-ditch of Blankenrode, Eastern Westphalia, as well as Viola tricolor ecotypes of Eastern Europe, keep their cells free of excess heavy metals by arbuscular mycorrhizal fungi which bind heavy metals. The Caryophyllaceae, Silene vulgaris f. humilis and Minuartia verna, apparently discard leaves when overloaded with heavy metals. All Central European metallophytes have close relatives that grow in areas outside of heavy metal soils, mainly in the Alps, and have, therefore, been considered as relicts of the glacial epoch in the past. However, the current literature favours the idea that hyperaccumulation of heavy metals serves plants as deterrent against attack by feeding animals (termed elemental defense hypothesis). The capability to hyperaccumulate heavy metals in A. halleri and N. caerulescens is achieved by duplications and alterations of the cis-regulatory properties of genes coding for heavy metal transporting/excreting proteins. Several metallophytes have developed ecotypes with a varying content of such heavy metal transporters as an adaption to the specific toxicity of a heavy metal site.

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.

Molecular Resources from Transcriptomes in the Brassicaceae Family

The rapidly falling costs and the increasing availability of large DNA sequence data sets facilitate the fast and affordable mining of large molecular markers data sets for comprehensive evolutionary studies. The Brassicaceae (mustards) are an important species-rich family in the plant kingdom with taxa distributed worldwide and a complex evolutionary history. We performed Simple Sequence Repeats (SSRs) mining using de novo assembled transcriptomes from 19 species across the Brassicaceae in order to study SSR evolution and provide comprehensive sets of molecular markers for genetic studies within the family. Moreover, we selected the genus Cochlearia to test the transferability and polymorphism of these markers among species. Additionally, we annotated Cochlearia pyrenaica transcriptome in order to identify the position of each of the mined SSRs. While we introduce a new set of tools that will further enable evolutionary studies across the Brassicaceae, we also discuss some broader aspects of SSR evolution. Overall, we developed 2012 ready-to-use SSR markers with their respective primers in 19 Brassicaceae species and a high quality annotated transcriptome for C. pyrenaica. As indicated by our transferability test with the genus Cochlearia these SSRs are transferable to species within the genus increasing exponentially the number of targeted species. Also, our polymorphism results showed substantial levels of variability for these markers. Finally, despite its complex evolutionary history, SSR evolution across the Brassicaceae family is highly conserved and we found no deviation from patterns reported in other Angiosperms.

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.

Comparative analysis of complete plastid genomes from wild soybean (Glycine soja) and nine other Glycine species

The plastid genomes of different plant species exhibit significant variation, thereby providing valuable markers for exploring evolutionary relationships and population genetics. Glycine soja (wild soybean) is recognized as the wild ancestor of cultivated soybean (G. max), representing a valuable genetic resource for soybean breeding programmes. In the present study, the complete plastid genome of G. soja was sequenced using Illumina paired-end sequencing and then compared it for the first time with previously reported plastid genome sequences from nine other Glycine species. The G. soja plastid genome was 152,224 bp in length and possessed a typical quadripartite structure, consisting of a pair of inverted repeats (IRa/IRb; 25,574 bp) separated by small (178,963 bp) and large (83,181 bp) single-copy regions, with a 51-kb inversion in the large single-copy region. The genome encoded 134 genes, including 87 protein-coding genes, eight ribosomal RNA genes, and 39 transfer RNA genes, and possessed 204 randomly distributed microsatellites, including 15 forward, 25 tandem, and 34 palindromic repeats. Whole-plastid genome comparisons revealed an overall high degree of sequence similarity between G. max and G. gracilis and some divergence in the intergenic spacers of other species. Greater numbers of indels and SNP substitutions were observed compared with G. cyrtoloba. The sequence of the accD gene from G. soja was highly divergent from those of the other species except for G. max and G. gracilis. Phylogenomic analyses of the complete plastid genomes and 76 shared genes yielded an identical topology and indicated that G. soja is closely related to G. max and G. gracilis. The complete G. soja genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of Glycine species and can be used to identify related 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.

Phylogeography of Arabidopsis halleri (Brassicaceae) in mountain regions of Central Europe inferred from cpDNA variation and ecological niche modelling

The present study aimed to investigate phylogeographical patterns present within A. halleri in Central Europe. 1,281 accessions sampled from 52 populations within the investigated area were used in the study of genetic variation based on chloroplast DNA. Over 500 high-quality species occurrence records were used in ecological niche modelling experiments. We evidenced the presence of a clear phylogeographic structure within A. halleri in Central Europe. Our results showed that two genetically different groups of populations are present in western and eastern part of the Carpathians. The hypothesis of the existence of a glacial refugium in the Western Carpathians adn the Bohemian Forest cannot be rejected from our data. It seems, however, that the evidence collected during the present study is not conclusive. The area of Sudetes was colonised after LGM probably by migrants from the Bohemian Forest.

A Time-Calibrated Road Map of Brassicaceae Species Radiation and Evolutionary History

The Brassicaceae include several major crop plants and numerous important model species in comparative evolutionary research such as Arabidopsis, Brassica, Boechera, Thellungiella, and Arabis species. As any evolutionary hypothesis needs to be placed in a temporal context, reliably dated major splits within the evolution of Brassicaceae are essential. We present a comprehensive time-calibrated framework with important divergence time estimates based on whole-chloroplast sequence data for 29 Brassicaceae species. Diversification of the Brassicaceae crown group started at the Eocene-to-Oligocene transition. Subsequent major evolutionary splits are dated to ∼20 million years ago, coinciding with the Oligocene-to-Miocene transition, with increasing drought and aridity and transient glaciation events. The age of the Arabidopsis thaliana crown group is 6 million years ago, at the Miocene and Pliocene border. The overall species richness of the family is well explained by high levels of neopolyploidy (43% in total), but this trend is neither directly associated with an increase in genome size nor is there a general lineage-specific constraint. Our results highlight polyploidization as an important source for generating new evolutionary lineages adapted to changing environments. We conclude that species radiation, paralleled by high levels of neopolyploidization, follows genome size decrease, stabilization, and genetic diploidization.

Planting molecular functions in an ecological context with Arabidopsis thaliana

The vascular plant Arabidopsis thaliana is a central genetic model and universal reference organism in plant and crop science. The successful integration of different fields of research in the study of A. thaliana has made a large contribution to our molecular understanding of key concepts in biology. The availability and active development of experimental tools and resources, in combination with the accessibility of a wealth of cumulatively acquired knowledge about this plant, support the most advanced systems biology approaches among all land plants. Research in molecular ecology and evolution has also brought the natural history of A. thaliana into the limelight. This article showcases our current knowledge of the natural history of A. thaliana from the perspective of the most closely related plant species, providing an evolutionary framework for interpreting novel findings and for developing new hypotheses based on our knowledge of this plant.

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.