R-gene variation across Arabidopsis lyrata subspecies: effects of population structure, selection and mating system

Transposable elements maintain genome-wide heterozygosity in inbred populations

Elevated levels of inbreeding increase the risk of inbreeding depression and extinction, yet many inbred species are widespread, suggesting that inbreeding has little impact on evolutionary potential. Here, we explore the potential for transposable elements (TEs) to maintain genetic variation in functional genomic regions under extreme inbreeding. Capitalizing on the mixed mating system of Arabidopsis lyrata, we assess genome-wide heterozygosity and signatures of selection at single nucleotide polymorphisms near transposable elements across an inbreeding gradient. Under intense inbreeding, we find systematically elevated heterozygosity downstream of several TE superfamilies, associated with signatures of balancing selection. In addition, we demonstrate increased heterozygosity in stress-responsive genes that consistently occur downstream of TEs. We finally reveal that TE superfamilies are associated with specific signatures of selection that are reproducible across independent evolutionary lineages of A. lyrata. Together, our study provides an important hypothesis for the success of self-fertilizing species.

Early evolution of reproductive isolation: A case of weak inbreeder/strong outbreeder leads to an intraspecific hybridization barrier in Arabidopsis lyrata

Reproductive strategies play a major role in plant speciation. Notably, transitions from outcrossing to selfing may lead to relaxed sexual selection and parental conflict. Shifts in mating systems can affect maternal and paternal interests, and thus parent-specific influence on endosperm development, leading to reproductive isolation: if selfing and outcrossing species hybridize, the resulting seeds may not be viable due to endosperm failure. Nevertheless, it remains unclear how the switch in mating systems can impact reproductive isolation between recently diverged lineages, that is, during the process of speciation. We investigated this question using Arabidopsis lyrata, which recently transitioned to selfing (10,000 years ago) in certain North American populations, where European populations remain outcrossing. We performed reciprocal crosses between selfers and outcrossers, and measured seed viability and endosperm development. We show that parental genomes in the hybrid seed negatively interact, as predicted by parental conflict. This leads to extensive hybrid seed lethality associated with endosperm cellularization disturbance. Our results suggest that this is primarily driven by divergent evolution of the paternal genome between selfers and outcrossers. In addition, we observed other hybrid seed defects, suggesting that sex-specific interests are not the only processes contributing to postzygotic reproductive isolation.

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.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

Changing environments and genetic variation: natural variation in inbreeding does not compromise short-term physiological responses

Selfing plant lineages are surprisingly widespread and successful in a broad range of environments, despite showing reduced genetic diversity, which is predicted to reduce their long-term evolutionary potential. However, appropriate short-term plastic responses to new environmental conditions might not require high levels of standing genetic variation. In this study, we tested whether mating system variation among populations, and associated changes in genetic variability, affected short-term responses to environmental challenges. We compared relative fitness and metabolome profiles of naturally outbreeding (genetically diverse) and inbreeding (genetically depauperate) populations of a perennial plant, Arabidopsis lyrata, under constant growth chamber conditions and an outdoor common garden environment outside its native range. We found no effect of inbreeding on survival, flowering phenology or short-term physiological responses. Specifically, naturally occurring inbreeding had no significant effects on the plasticity of metabolome profiles, using either multivariate approaches or analysis of variation in individual metabolites, with inbreeding populations showing similar physiological responses to outbreeding populations over time in both growing environments. We conclude that low genetic diversity in naturally inbred populations may not always compromise fitness or short-term physiological capacity to respond to environmental change, which could help to explain the global success of selfing mating strategies.

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.

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.

Non-host Resistance Induced by the Xanthomonas Effector XopQ Is Widespread within the Genus Nicotiana and Functionally Depends on EDS1

Most Gram-negative plant pathogenic bacteria translocate effector proteins (T3Es) directly into plant cells via a conserved type III secretion system, which is essential for pathogenicity in susceptible plants. In resistant plants, recognition of some T3Es is mediated by corresponding resistance (R) genes or R proteins and induces effector triggered immunity (ETI) that often results in programmed cell death reactions. The identification of R genes and understanding their evolution/distribution bears great potential for the generation of resistant crop plants. We focus on T3Es from Xanthomonas campestris pv. vesicatoria (Xcv), the causal agent of bacterial spot disease on pepper and tomato plants. Here, 86 Solanaceae lines mainly of the genus Nicotiana were screened for phenotypical reactions after Agrobacterium tumefaciens-mediated transient expression of 21 different Xcv effectors to (i) identify new plant lines for T3E characterization, (ii) analyze conservation/evolution of putative R genes and (iii) identify promising plant lines as repertoire for R gene isolation. The effectors provoked different reactions on closely related plant lines indicative of a high variability and evolution rate of potential R genes. In some cases, putative R genes were conserved within a plant species but not within superordinate phylogenetical units. Interestingly, the effector XopQ was recognized by several Nicotiana spp. lines, and Xcv infection assays revealed that XopQ is a host range determinant in many Nicotiana species. Non-host resistance against Xcv and XopQ recognition in N. benthamiana required EDS1, strongly suggesting the presence of a TIR domain-containing XopQ-specific R protein in these plant lines. XopQ is a conserved effector among most xanthomonads, pointing out the XopQ-recognizing R_xopQ as candidate for targeted crop improvement.