Functional variation in a disease resistance gene in populations ofArabidopsis thaliana

A Genome-Wide Association study in Arabidopsis thaliana to decipher the adaptive genetics of quantitative disease resistance in a native heterogeneous environment

Pathogens are often the main selective agents acting in plant communities, thereby influencing the distribution of polymorphism at loci affecting resistance within and among natural plant populations. In addition, the outcome of plant-pathogen interactions can be drastically affected by abiotic and biotic factors at different spatial and temporal grains. The characterization of the adaptive genetic architecture of disease resistance in native heterogeneous environments is however still missing. In this study, we conducted an in situ Genome-Wide Association study in the spatially heterogeneous native habitat of a highly genetically polymorphic local mapping population of Arabidopsis thaliana, to unravel the adaptive genetic architecture of quantitative disease resistance. Disease resistance largely differed among three native soils and was affected by the presence of the grass Poa annua. The observation of strong crossing reactions norms among the 195 A. thaliana genotypes for disease resistance among micro-habitats, combined with a negative fecundity-disease resistance relationship in each micro-habitat, suggest that alternative local genotypes of A. thaliana are favored under contrasting environmental conditions at the scale of few meters. A complex genetic architecture was detected for disease resistance and fecundity. However, only few QTLs were common between these two traits. Heterogeneous selection in this local population should therefore promote the maintenance of polymorphism at only few candidate resistance genes.

Toward Unifying Evolutionary Ecology and Genomics to Understand Positive Plant-Plant Interactions Within Wild Species

In a local environment, plant networks include interactions among individuals of different species and among genotypes of the same species. While interspecific interactions are recognized as main drivers of plant community patterns, intraspecific interactions have recently gained attention in explaining plant community dynamics. However, an overview of intraspecific genotype-by-genotype interaction patterns within wild plant species is still missing. From the literature, we identified 91 experiments that were mainly designed to investigate the presence of positive interactions based on two contrasting hypotheses. Kin selection theory predicts partisan help given to a genealogical relative. The rationale behind this hypothesis relies on kin/non-kin recognition, with the positive outcome of kin cooperation substantiating it. On the other hand, the elbow-room hypothesis supports intraspecific niche partitioning leading to positive outcome when genetically distant genotypes interact. Positive diversity-productivity relationship rationalizes this hypothesis, notably with the outcome of overyielding. We found that both these hypotheses have been highly supported in experimental studies despite their opposite predictions between the extent of genetic relatedness among neighbors and the level of positive interactions. Interestingly, we identified a highly significant effect of breeding system, with a high proportion of selfing species associated with the presence of kin cooperation. Nonetheless, we identified several shortcomings regardless of the species considered, such as the lack of a reliable estimate of genetic relatedness among genotypes and ecological characterization of the natural habitats from which genotypes were collected, thereby impeding the identification of selective drivers of positive interactions. We therefore propose a framework combining evolutionary ecology and genomics to establish the eco-genomic landscape of positive GxG interactions in wild plant species.

Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population

Background

Pathogens are key components in natural and agricultural plant systems. There is evidence of evolutionary changes in disease susceptibility as a consequence of climate change, but we know little about the underlying genetic basis of this evolution. To address this, we took advantage of a historical seed collection of a Brassica rapa population, which we previously demonstrated evolved an increase in disease susceptibility to a necrotrophic fungal pathogen following a drought.

Results

Previously, we combined a resurrection experiment with genome-wide sequencing of 124 pooled ancestral and descendant plants. Here, using these previously generated sequence data (Franks et al. in Mol Ecol 25(15):3622–3631, 2016), we show that well-characterized necrotrophic fungal pathogen response (NFPR) genes have evolved, as indicated by changes in allele frequency, between ancestors and descendants, with several of them identified as extreme F_ST outliers. The jasmonic acid (JA) signaling pathway in particular seems to underlie the evolution of disease susceptibility, in addition to its well characterized role in plastic disease response. We identify a list of 260 genes that are both NFPR genes and are differentially expressed in response to drought, based on publicly available data. We present evidence that five of these genes evolved between ancestors and descendants, suggesting that the drought acted as the evolutionary driver, and that the accompanying increase in disease susceptibility may have been a consequence of genetic pleiotropy.

Conclusions

Our study provides evidence that for this population, standing variation in NFPR genes is affected by natural selection related to climate change. Our results reveal potentially important candidates that may underlie trait evolution in both crops and natural systems. Additionally, this trade-off between adaptation to biotic and abiotic stresses is an example of how climate change can have diverse and unexpected consequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01789-7.

Arabidopsis thaliana Genes Associated with Cucumber mosaic virus Virulence and Their Link to Virus Seed Transmission

Virulence, the effect of pathogen infection on progeny production, is a major determinant of host and pathogen fitness as it affects host fecundity and pathogen transmission. In plant-virus interactions, ample evidence indicates that virulence is genetically controlled by both partners. However, the host genetic determinants are poorly understood. Through a genome-wide association study (GWAS) of 154 Arabidopsis thaliana genotypes infected by Cucumber mosaic virus (CMV), we identified eight host genes associated with virulence, most of them involved in response to biotic stresses and in cell wall biogenesis in plant reproductive structures. Given that virulence is a main determinant of the efficiency of plant virus seed transmission, we explored the link between this trait and the genetic regulation of virulence. Our results suggest that the same functions that control virulence are also important for CMV transmission through seeds. In sum, this work provides evidence of a novel role for some previously known plant defense genes and for the cell wall metabolism in plant virus interactions.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

Atypical Resistance Protein RPW8/HR Triggers Oligomerization of the NLR Immune Receptor RPP7 and Autoimmunity

In certain plant hybrids, immunity signaling is initiated when immune components interact in the absence of a pathogen trigger. In Arabidopsis thaliana, such autoimmunity and cell death are linked to variants of the NLR RPP7 and the RPW8 proteins involved in broad-spectrum resistance. We uncover the molecular basis for this autoimmunity and demonstrate that a homolog of RPW8, HR4Fei-0, can trigger the assembly of a higher-order RPP7 complex, with autoimmunity signaling as a consequence. HR4Fei-0-mediated RPP7 oligomerization occurs via the RPP7 C-terminal leucine-rich repeat (LRR) domain and ATP-binding P-loop. RPP7 forms a higher-order complex only in the presence of HR4Fei-0 and not with the standard HR4 variant, which is distinguished from HR4Fei-0 by length variation in C-terminal repeats. Additionally, HR4Fei-0 can independently form self-oligomers, which directly kill cells in an RPP7-independent manner. Our work provides evidence for a plant resistosome complex and the mechanisms by which RPW8/HR proteins trigger cell death.

RPW8/HR repeats control NLR activation in Arabidopsis thaliana

In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant.

Cucumber mosaic virus infection as a potential selective pressure on Arabidopsis thaliana populations

It has been proposed that in wild ecosystems viruses are often plant mutualists, whereas agroecosystems favour pathogenicity. We seek evidence for virus pathogenicity in wild ecosystems through the analysis of plant-virus coevolution, which requires a negative effect of infection on the host fitness. We focus on the interaction between Arabidopsis thaliana and Cucumber mosaic virus (CMV), which is significant in nature. We studied the genetic diversity of A. thaliana for two defence traits, resistance and tolerance, to CMV. A set of 185 individuals collected in 76 A. thaliana Iberian wild populations were inoculated with different CMV strains. Resistance was estimated from the level of virus multiplication in infected plants, and tolerance from the effect of infection on host progeny production. Resistance and tolerance to CMV showed substantial genetic variation within and between host populations, and depended on the virus x host genotype interaction, two conditions for coevolution. Resistance and tolerance were co-occurring independent traits that have evolved independently from related life-history traits involved in adaptation to climate. The comparison of the genetic structure for resistance and tolerance with that for neutral traits (QST/FST analyses) indicated that both defence traits are likely under uniform selection. These results strongly suggest that CMV infection selects for defence on A. thaliana populations, and support plant-virus coevolution. Thus, we propose that CMV infection reduces host fitness under the field conditions of the wild A. thaliana populations studied.

The Genetics Underlying Natural Variation in the Biotic Interactions of Arabidopsis thaliana: The Challenges of Linking Evolutionary Genetics and Community Ecology

In the context of global change, predicting the responses of plant communities in an ever-changing biotic environment calls for a multipronged approach at the interface of evolutionary genetics and community ecology. However, our understanding of the genetic basis of natural variation involved in mediating biotic interactions, and associated adaptive dynamics of focal plants in their natural communities, is still in its infancy. Here, we review the genetic and molecular bases of natural variation in the response to biotic interactions (viruses, bacteria, fungi, oomycetes, herbivores, and plants) in the model plant Arabidopsis thaliana as well as the adaptive value of these bases. Among the 60 identified genes are a number that encode nucleotide-binding site leucine-rich repeat (NBS-LRR)-type proteins, consistent with early examples of plant defense genes. However, recent studies have revealed an extensive diversity in the molecular mechanisms of defense. Many types of genetic variants associate with phenotypic variation in biotic interactions, even among the genes of large effect that tend to be identified. In general, we found that (i) balancing selection rather than directional selection explains the observed patterns of genetic diversity within A. thaliana and (ii) the cost/benefit tradeoffs of adaptive alleles can be strongly dependent on both genomic and environmental contexts. Finally, because A. thaliana rarely interacts with only one biotic partner in nature, we highlight the benefit of exploring diffuse biotic interactions rather than tightly associated host-enemy pairs. This challenge would help to improve our understanding of coevolutionary quantitative genetics within the context of realistic community complexity.

Lagging adaptation to warming climate in Arabidopsis thaliana

If climate change outpaces the rate of adaptive evolution within a site, populations previously well adapted to local conditions may decline or disappear, and banked seeds from those populations will be unsuitable for restoring them. However, if such adaptational lag has occurred, immigrants from historically warmer climates will outperform natives and may provide genetic potential for evolutionary rescue. We tested for lagging adaptation to warming climate using banked seeds of the annual weed Arabidopsis thaliana in common garden experiments in four sites across the species' native European range: Valencia, Spain; Norwich, United Kingdom; Halle, Germany; and Oulu, Finland. Genotypes originating from geographic regions near the planting site had high relative fitness in each site, direct evidence for broad-scale geographic adaptation in this model species. However, genotypes originating in sites historically warmer than the planting site had higher average relative fitness than local genotypes in every site, especially at the northern range limit in Finland. This result suggests that local adaptive optima have shifted rapidly with recent warming across the species' native range. Climatic optima also differed among seasonal germination cohorts within the Norwich site, suggesting that populations occurring where summer germination is common may have greater evolutionary potential to persist under future warming. If adaptational lag has occurred over just a few decades in banked seeds of an annual species, it may be an important consideration for managing longer-lived species, as well as for attempts to conserve threatened populations through ex situ preservation.

Loss of compatibility might explain resistance of the Arabidopsis thaliana accession Te-0 to Golovinomyces cichoracearum

BACKGROUND

The establishment of compatibility between plants and pathogens requires compliance with various conditions, such as recognition of the right host, suppression of defence mechanisms, and maintenance of an environment allowing pathogen reproduction. To date, most of the plant factors required to sustain compatibility remain unknown, with the few best characterized being those interfering with defence responses. A suitable system to study host compatibility factors is the interaction between Arabidopsis thaliana and the powdery mildew (PM) Golovinomyces cichoracearum. As an obligate biotrophic pathogen, this fungus must establish compatibility in order to perpetuate. In turn, A. thaliana displays natural variation for susceptibility to this invader, with some accessions showing full susceptibility (Col-0), and others monogenic dominant resistance (Kas-1). Interestingly, Te-0, among other accessions, displays recessive partial resistance to this PM.

RESULTS

In this study, we characterized the interaction of G. cichoracearum with Te-0 plants to investigate the basis of this plant resistance. We found that Te-0's incompatibility was not associated with hyper-activation of host inducible defences. Te-0 plants allowed germination of conidia and development of functional haustoria, but could not support the formation of mature conidiophores. Using a suppressive subtractive hybridization technique, we identified plant genes showing differential expression between resistant Te-0 and susceptible Col-0 plants at the fungal pre-conidiation stage.

CONCLUSIONS

Te-0 resistance is likely caused by loss of host compatibility and not by stimulation of inducible defences. Conidiophores formation is the main constraint for completion of fungal life cycle in Te-0 plants. The system here described allowed the identification of genes proposed as markers for susceptibility to this PM.

The effect of environmental heterogeneity on RPW8-mediated resistance to powdery mildews in Arabidopsis thaliana

BACKGROUND AND AIMS

The biotic and abiotic environment of interacting hosts and parasites may vary considerably over small spatial and temporal scales. It is essential to understand how different environments affect host disease resistance because this determines frequency of disease and, importantly, heterogeneous environments can retard direct selection and potentially maintain genetic variation for resistance in natural populations.

METHODS

The effect of different temperatures and soil nutrient conditions on the outcome of infection by a pathogen was quantified in Arabidopsis thaliana. Expression levels of a gene conferring resistance to powdery mildews, RPW8, were compared with levels of disease to test a possible mechanism behind variation in resistance.

KEY RESULTS

Most host genotypes changed from susceptible to resistant across environments with the ranking of genotypes differing between treatments. Transcription levels of RPW8 increased after infection and varied between environments, but there was no tight association between transcription and resistance levels.

CONCLUSIONS

There is a strong potential for a heterogeneous environment to change the resistance capacity of A. thaliana genotypes and hence the direction and magnitude of selection in the presence of the pathogen. Possible causative links between resistance gene expression and disease resistance are discussed in light of the present results on RPW8.

Complex evolutionary events at a tandem cluster of Arabidopsis thaliana genes resulting in a single-locus genetic incompatibility

Non-additive interactions between genomes have important implications, not only for practical applications such as breeding, but also for understanding evolution. In extreme cases, genes from different genomic backgrounds may be incompatible and compromise normal development or physiology. Of particular interest are non-additive interactions of alleles at the same locus. For example, overdominant behavior of alleles, with respect to plant fitness, has been proposed as an important component of hybrid vigor, while underdominance may lead to reproductive isolation. Despite their importance, only a few cases of genetic over- or underdominance affecting plant growth or fitness are understood at the level of individual genes. Moreover, the relationship between biochemical and fitness effects may be complex: genetic overdominance, that is, increased or novel activity of a gene may lead to evolutionary underdominance expressed as hybrid weakness. Here, we describe a non-additive interaction between alleles at the Arabidopsis thaliana OAK (OUTGROWTH-ASSOCIATED PROTEIN KINASE) gene. OAK alleles from two different accessions interact in F₁ hybrids to cause a variety of aberrant growth phenotypes that depend on a recently acquired promoter with a novel expression pattern. The OAK gene, which is located in a highly variable tandem array encoding closely related receptor-like kinases, is found in one third of A. thaliana accessions, but not in the reference accession Col-0. Besides recruitment of exons from nearby genes as promoter sequences, key events in OAK evolution include gene duplication and divergence of a potential ligand-binding domain. OAK kinase activity is required for the aberrant phenotypes, indicating it is not recognition of an aberrant protein, but rather a true gain of function, or overdominance for gene activity, that leads to this underdominance for fitness. Our work provides insights into how tandem arrays, which are particularly prone to frequent, complex rearrangements, can produce genetic novelty.

While intraspecific hybrids are vitally important in modern agriculture because they often perform better than their inbred parents, certain hybrid combinations fail to develop normally and are inferior to their parents. We have identified an Arabidopsis thaliana hybrid with several aberrant growth phenotypes that are caused by divergence at a single locus encoding the receptor-like kinase OUTGROWTH-ASSOCIATED PROTEIN KINASE (OAK). OAK belongs to a group of similar genes arranged in a tandem cluster that varies substantially between A. thaliana strains. OAK originated through duplication within the cluster with concurrent recruitment of coding sequences from nearby genes to form a new promoter with a novel expression pattern. Kinase activity of OAK is required for its effects, indicating that it is not recognition of an aberrant protein but rather a true gain of function that leads to the incompatibility. Most of the incompatibility seems to come from divergence within the extracellular ligand-binding domain of the OAK protein, indicating that heterodimers of OAK may have higher affinity for a natural substrate compared to either homodimer. Finally, mis-expression of the incompatible OAK alleles from the promoter present in the reference strain of A. thaliana also leads to genetic incompatibility, but with different phenotypic outcomes.

Spatial variation in disease resistance: from molecules to metapopulations

Variation in disease resistance is a widespread phenomenon in wild plant-pathogen associations. Here, we review current literature on natural plant-pathogen associations to determine how diversity in disease resistance is distributed at different hierarchical levels - within host individuals, within host populations, among host populations at the metapopulation scale and at larger regional scales.We find diversity in resistance across all spatial scales examined. Furthermore, variability seems to be the best counter-defence of plants against their rapidly evolving pathogens. We find that higher diversity of resistance phenotypes also results in higher levels of resistance at the population level.Overall, we find that wild plant populations are more likely to be susceptible than resistant to their pathogens. However, the degree of resistance differs strikingly depending on the origin of the pathogen strains used in experimental inoculation studies. Plant populations are on average 16% more resistant to allopatric pathogen strains than they are to strains that occur within the same population (48 % vs. 32 % respectively).Pathogen dispersal mode affects levels of resistance in natural plant populations with lowest levels detected for hosts of airborne pathogens and highest for waterborne pathogens.Detailed analysis of two model systems, Linum marginale infected by Melampsora lini, and Plantago lanceolata infected by Podosphaera plantaginis, show that the amount of variation in disease resistance declines towards higher spatial scales as we move from individual hosts to metapopulations, but evaluation of multiple spatial scales is needed to fully capture the structure of disease resistance.Synthesis: Variation in disease resistance is ubiquitous in wild plant-pathogen associations. While the debate over whether the resistance structure of plant populations is determined by pathogen-imposed selection versus non-adaptive processes remains unresolved, we do report examples of pathogen-imposed selection on host resistance. Here we highlight the importance of measuring resistance across multiple spatial scales, and of using sympatric strains when looking for signs of coevolution in wild plant-pathogen interactions.

Local-scale patterns of genetic variability, outcrossing, and spatial structure in natural stands of Arabidopsis thaliana

As Arabidopsis thaliana is increasingly employed in evolutionary and ecological studies, it is essential to understand patterns of natural genetic variation and the forces that shape them. Previous work focusing mostly on global and regional scales has demonstrated the importance of historical events such as long-distance migration and colonization. Far less is known about the role of contemporary factors or environmental heterogeneity in generating diversity patterns at local scales. We sampled 1,005 individuals from 77 closely spaced stands in diverse settings around Tübingen, Germany. A set of 436 SNP markers was used to characterize genome-wide patterns of relatedness and recombination. Neighboring genotypes often shared mosaic blocks of alternating marker identity and divergence. We detected recent outcrossing as well as stretches of residual heterozygosity in largely homozygous recombinants. As has been observed for several other selfing species, there was considerable heterogeneity among sites in diversity and outcrossing, with rural stands exhibiting greater diversity and heterozygosity than urban stands. Fine-scale spatial structure was evident as well. Within stands, spatial structure correlated negatively with observed heterozygosity, suggesting that the high homozygosity of natural A. thaliana may be partially attributable to nearest-neighbor mating of related individuals. The large number of markers and extensive local sampling employed here afforded unusual power to characterize local genetic patterns. Contemporary processes such as ongoing outcrossing play an important role in determining distribution of genetic diversity at this scale. Local "outcrossing hotspots" appear to reshuffle genetic information at surprising rates, while other stands contribute comparatively little. Our findings have important implications for sampling and interpreting diversity among A. thaliana accessions.

RPW8 and resistance to powdery mildew pathogens in natural populations of Arabidopsis lyrata

It is not clear to what extent the orthologues of genes that are adaptively important in one species also contribute to adaptive variation in others. Here, we examine Arabidopsis lyrata to assess the functional and evolutionary significance of natural variation in an orthologue of the gene RPW8 known to be a major determinant of powdery mildew resistance in Arabidopsis thaliana. We assessed the sequence variation at RPW8 and the associated resistance reaction in populations of A. lyrata ssp. petraea. Neutrality tests were performed to understand the importance of local adaptation in maintaining variation at the locus. Highly truncated RPW8 proteins were frequent in all populations and were associated with an increased risk of susceptibility. Haplotypes encoding full-length proteins were highly significantly associated with resistance. There were no signatures of selection at the species-wide level, but some evidence for positive selection in two populations. RPW8 in A. lyrata appears to have a role in powdery mildew resistance, similar to its orthologue in A. thaliana. Unlike A. thaliana, A. lyrata contains a genetic component that can act independently of RPW8 to confer resistance to powdery mildew pathogens. Infrequent local selective sweeps may favour different alleles in different populations, and thereby contribute to the maintenance of species-wide variation at the locus.