Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection.

Gene duplication is increasingly recognized as an important mechanism for the origination of new genes, as revealed by comparative genomic analysis. However, how new duplicate genes contribute to phenotypic evolution remains largely unknown, especially in plants. Here, we identified the new gene EXOV, derived from a partial gene duplication of its parental gene EXOVL in Arabidopsis thaliana. EXOV is a species-specific gene that originated within the last 3.5 million years and shows strong signals of positive selection. Unexpectedly, RNA-sequencing analyses revealed that, despite its young age, EXOV has acquired many novel direct and indirect interactions in which the parental gene does not engage. This observation is consistent with the high, selection-driven substitution rate of its encoded protein, in contrast to the slowly evolving EXOVL, suggesting an important role for EXOV in phenotypic evolution. We observed significant differentiation of morphological changes for all phenotypes assessed in genome-edited and T-DNA insertional single mutants and in double T-DNA insertion mutants in EXOV and EXOVL. We discovered a substantial divergence of phenotypic effects by principal component analyses, suggesting neofunctionalization of the new gene. These results reveal a young gene that plays critical roles in biological processes that underlie morphological evolution in A. thaliana.

Plant genetic effects on microbial hubs impact host fitness in repeated field trials.

Although complex interactions between hosts and microbial associates are increasingly well documented, we still know little about how and why hosts shape microbial communities in nature. In addition, host genetic effects on microbial communities vary widely depending on the environment, obscuring conclusions about which microbes are impacted and which plant functions are important. We characterized the leaf microbiota of 200 Arabidopsis thaliana genotypes in eight field experiments and detected consistent host effects on specific, broadly distributed microbial species (operational taxonomic unit [OTUs]). Host genetic effects disproportionately influenced central ecological hubs, with heritability of particular OTUs declining with their distance from the nearest hub within the microbial network. These host effects could reflect either OTUs preferentially associating with specific genotypes or differential microbial success within them. Host genetics associated with microbial hubs explained over 10% of the variation in lifetime seed production among host genotypes across sites and years. We successfully cultured one of these microbial hubs and demonstrated its growth-promoting effects on plants in sterile conditions. Finally, genome-wide association mapping identified many putatively causal genes with small effects on the relative abundance of microbial hubs across sites and years, and these genes were enriched for those involved in the synthesis of specialized metabolites, auxins, and the immune system. Using untargeted metabolomics, we corroborate the consistent association between variation in specialized metabolites and microbial hubs across field sites. Together, our results reveal that host genetic variation impacts the microbial communities in consistent ways across environments and that these effects contribute to fitness variation among host genotypes.

Multiple FLC haplotypes defined by independent cis-regulatory variation underpin life history diversity in Arabidopsis thaliana.

Relating molecular variation to phenotypic diversity is a central goal in evolutionary biology. In Arabidopsis thaliana, FLOWERING LOCUS C (FLC) is a major determinant of variation in vernalization--the acceleration of flowering by prolonged cold. Here, through analysis of 1307 A. thaliana accessions, we identify five predominant FLC haplotypes defined by noncoding sequence variation. Genetic and transgenic experiments show that they are functionally distinct, varying in FLC expression level and rate of epigenetic silencing. Allelic heterogeneity at this single locus accounts for a large proportion of natural variation in vernalization that contributes to adaptation of A. thaliana.

Two-way mixed-effects methods for joint association analysis using both host and pathogen genomes.

Infectious diseases are often affected by specific pairings of hosts and pathogens and therefore by both of their genomes. The integration of a pair of genomes into genome-wide association mapping can provide an exquisitely detailed view of the genetic landscape of complex traits. We present a statistical method, ATOMM (Analysis with a Two-Organism Mixed Model), that maps a trait of interest to a pair of genomes simultaneously; this method makes use of whole-genome sequence data for both host and pathogen organisms. ATOMM uses a two-way mixed-effect model to test for genetic associations and cross-species genetic interactions while accounting for sample structure including interactions between the genetic backgrounds of the two organisms. We demonstrate the applicability of ATOMM to a joint association study of quantitative disease resistance (QDR) in the Arabidopsis thaliana-Xanthomonas arboricola pathosystem. Our method uncovers a clear host-strain specificity in QDR and provides a powerful approach to identify genetic variants on both genomes that contribute to phenotypic variation.

The rate and potential relevance of new mutations in a colonizing plant lineage.

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by "natural evolution experiments" that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions-the subset of mutations that survived natural selection and drift-, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

Genome-wide association studies on the phyllosphere microbiome: Embracing complexity in host-microbe interactions.

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe-related traits were developed in the relatively simple context of binary host-microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome-wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant-microbe interactions will influence the choice of traits to be mapped.

Mechanisms to Mitigate the Trade-Off between Growth and Defense.

Plants have evolved an array of defenses against pathogens. However, mounting a defense response frequently comes with the cost of a reduction in growth and reproduction, carrying critical implications for natural and agricultural populations. This review focuses on how costs are generated and whether and how they can be mitigated. Most well-characterized growth-defense trade-offs stem from antagonistic crosstalk among hormones rather than an identified metabolic expenditure. A primary way plants mitigate such costs is through restricted expression of resistance; this can be achieved through inducible expression of defense genes or by the concentration of defense to particular times or tissues. Defense pathways can be primed for more effective induction, and primed states can be transmitted to offspring. We examine the resistance (R) genes as a case study of how the toll of defense can be generated and ameliorated. The fine-scale regulation of R genes is critical to alleviate the burden of their expression, and the genomic organization of R genes into coregulatory modules reduces costs. Plants can also recruit protection from other species. Exciting new evidence indicates that a plant's genotype influences the microbiome composition, lending credence to the hypothesis that plants shape their microbiome to enhance defense.

The long-term maintenance of a resistance polymorphism through diffuse interactions.

Plant resistance (R) genes are a crucial component in plant defence against pathogens. Although R genes often fail to provide durable resistance in an agricultural context, they frequently persist as long-lived balanced polymorphisms in nature. Standard theory explains the maintenance of such polymorphisms through a balance of the costs and benefits of resistance and virulence in a tightly coevolving host-pathogen pair. However, many plant-pathogen interactions lack such specificity. Whether, and how, balanced polymorphisms are maintained in diffusely interacting species is unknown. Here we identify a naturally interacting R gene and effector pair in Arabidopsis thaliana and its facultative plant pathogen, Pseudomonas syringae. The protein encoded by the R gene RPS5 recognizes an AvrPphB homologue (AvrPphB2) and exhibits a balanced polymorphism that has been maintained for over 2 million years (ref. 3). Consistent with the presence of an ancient balanced polymorphism, the R gene confers a benefit when plants are infected with P. syringae carrying avrPphB2 but also incurs a large cost in the absence of infection. RPS5 alleles are maintained at intermediate frequencies in populations globally, suggesting ubiquitous selection for resistance. However, the presence of P. syringae carrying avrPphB is probably insufficient to explain the RPS5 polymorphism. First, avrPphB homologues occur at very low frequencies in P. syringae populations on A. thaliana. Second, AvrPphB only rarely confers a virulence benefit to P. syringae on A. thaliana. Instead, we find evidence that selection for RPS5 involves multiple non-homologous effectors and multiple pathogen species. These results and an associated model suggest that the R gene polymorphism in A. thaliana may not be maintained through a tightly coupled interaction involving a single coevolved R gene and effector pair. More likely, the stable polymorphism is maintained through complex and diffuse community-wide interactions.

Genome-wide association study reveals new loci involved in Arabidopsis thaliana and Turnip mosaic virus (TuMV) interactions in the field.

The genetic architecture of plant response to viruses has often been studied in model nonnatural pathosystems under controlled conditions. There is an urgent need to elucidate the genetic architecture of the response to viruses in a natural setting. A field experiment was performed in each of two years. In total, 317 Arabidopsis thaliana accessions were inoculated with its natural Turnip mosaic virus (TuMV). The accessions were phenotyped for viral accumulation, frequency of infected plants, stem length and symptoms. Genome-wide association mapping was performed. Arabidopsis thaliana exhibits extensive natural variation in its response to TuMV in the field. The underlying genetic architecture reveals a more quantitative picture than in controlled conditions. Ten genomic regions were consistently identified across the two years. RTM3 (Restricted TEV Movement 3) is a major candidate for the response to TuMV in the field. New candidate genes include Dead box helicase 1, a Tim Barrel domain protein and the eukaryotic translation initiation factor eIF3b. To our knowledge, this study is the first to report the genetic architecture of quantitative response of A. thaliana to a naturally occurring virus in a field environment, thereby highlighting relevant candidate genes involved in plant virus interactions in nature.

Differentiation between MAMP Triggered Defenses in Arabidopsis thaliana.

A first line of defense against pathogen attack for both plants and animals involves the detection of microbe-associated molecular patterns (MAMPs), followed by the induction of a complex immune response. Plants, like animals, encode several receptors that recognize different MAMPs. While these receptors are thought to function largely redundantly, the physiological responses to different MAMPs can differ in detail. Responses to MAMP exposure evolve quantitatively in natural populations of Arabidopsis thaliana, perhaps in response to environment specific differences in microbial threat. Here, we sought to determine the extent to which the detection of two canonical MAMPs were evolving redundantly or distinctly within natural populations. Our results reveal negligible correlation in plant growth responses between the bacterial MAMPs EF-Tu and flagellin. Further investigation of the genetic bases of differences in seedling growth inhibition and validation of 11 candidate genes reveal substantial differences in the genetic loci that underlie variation in response to these two MAMPs. Our results indicate that natural variation in MAMP recognition is largely MAMP-specific, indicating an ability to differentially tailor responses to EF-Tu and flagellin in A. thaliana populations.

Coselected genes determine adaptive variation in herbivore resistance throughout the native range of Arabidopsis thaliana.

The "mustard oil bomb" is a major defense mechanism in the Brassicaceae, which includes crops such as canola and the model plant Arabidopsis thaliana. These plants produce and store blends of amino acid-derived secondary metabolites called glucosinolates. Upon tissue rupture by natural enemies, the myrosinase enzyme hydrolyses glucosinolates, releasing defense molecules. Brassicaceae display extensive variation in the mixture of glucosinolates that they produce. To investigate the genetics underlying natural variation in glucosinolate profiles, we conducted a large genome-wide association study of 22 methionine-derived glucosinolates using A. thaliana accessions from across Europe. We found that 36% of among accession variation in overall glucosinolate profile was explained by genetic differentiation at only three known loci from the glucosinolate pathway. Glucosinolate-related SNPs were up to 490-fold enriched in the extreme tail of the genome-wide [Formula: see text] scan, indicating strong selection on loci controlling this pathway. Glucosinolate profiles displayed a striking longitudinal gradient with alkenyl and hydroxyalkenyl glucosinolates enriched in the West. We detected a significant contribution of glucosinolate loci toward general herbivore resistance and lifetime fitness in common garden experiments conducted in France, where accessions are enriched in hydroxyalkenyls. In addition to demonstrating the adaptive value of glucosinolate profile variation, we also detected long-distance linkage disequilibrium at two underlying loci, GS-OH and GS-ELONG. Locally cooccurring alleles at these loci display epistatic effects on herbivore resistance and fitness in ecologically realistic conditions. Together, our results suggest that natural selection has favored a locally adaptive configuration of physically unlinked loci in Western Europe.

Century-scale methylome stability in a recently diverged Arabidopsis thaliana lineage.

There has been much excitement about the possibility that exposure to specific environments can induce an ecological memory in the form of whole-sale, genome-wide epigenetic changes that are maintained over many generations. In the model plant Arabidopsis thaliana, numerous heritable DNA methylation differences have been identified in greenhouse-grown isogenic lines, but it remains unknown how natural, highly variable environments affect the rate and spectrum of such changes. Here we present detailed methylome analyses in a geographically dispersed A. thaliana population that constitutes a collection of near-isogenic lines, diverged for at least a century from a common ancestor. Methylome variation largely reflected genetic distance, and was in many aspects similar to that of lines raised in uniform conditions. Thus, even when plants are grown in varying and diverse natural sites, genome-wide epigenetic variation accumulates mostly in a clock-like manner, and epigenetic divergence thus parallels the pattern of genome-wide DNA sequence divergence.

A genome-wide survey reveals abundant rice blast R genes in resistant cultivars.

Plant resistance genes (R genes) harbor tremendous allelic diversity, constituting a robust immune system effective against microbial pathogens. Nevertheless, few functional R genes have been identified for even the best-studied pathosystems. Does this limited repertoire reflect specificity, with most R genes having been defeated by former pests, or do plants harbor a rich diversity of functional R genes, the composite behavior of which is yet to be characterized? Here, we survey 332 NBS-LRR genes cloned from five resistant Oryza sativa (rice) cultivars for their ability to confer recognition of 12 rice blast isolates when transformed into susceptible cultivars. Our survey reveals that 48.5% of the 132 NBS-LRR loci tested contain functional rice blast R genes, with most R genes deriving from multi-copy clades containing especially diversified loci. Each R gene recognized, on average, 2.42 of the 12 isolates screened. The abundant R genes identified in resistant genomes provide extraordinary redundancy in the ability of host genotypes to recognize particular isolates. If the same is true for other pathogens, many extant NBS-LRR genes retain functionality. Our success at identifying rice blast R genes also validates a highly efficient cloning and screening strategy.

Modulation of R-gene expression across environments.

Some environments are more conducive to pathogen growth than others, and, as a consequence, plants might be expected to invest more in resistance when pathogen growth is favored. Resistance (R-) genes in Arabidopsis thaliana have unusually extensive variation in basal expression when comparing the same R-gene among accessions collected from different environments. R-gene expression variation was characterized to explore whether R-gene expression is up-regulated in environments favoring pathogen proliferation and down-regulated when risks of infection are low; down-regulation would follow if costs of R-gene expression negatively impact plant fitness in the absence of disease. Quantitative reverse transcription-PCR was used to quantify the expression of 13 R-gene loci in plants grown in eight environmental conditions for each of 12 A. thaliana accessions, and large effects of the environment on R-gene expression were found. Surprisingly, almost every change in the environment--be it a change in biotic or abiotic conditions--led to an increase in R-gene expression, a response that was distinct from the average transcriptome response and from that of other stress response genes. These changes in expression are functional in that environmental change prior to infection affected levels of specific disease resistance to isolates of Pseudomonas syringae. In addition, there are strong latitudinal clines in basal R-gene expression and clines in R-gene expression plasticity correlated with drought and high temperatures. These results suggest that variation in R-gene expression across environments may be shaped by natural selection to reduce fitness costs of R-gene expression in permissive or predictable environments.

Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana.

A major challenge in evolutionary biology and plant breeding is to identify the genetic basis of complex quantitative traits, including those that contribute to adaptive variation. Here we review the development of new methods and resources to fine-map intraspecific genetic variation that underlies natural phenotypic variation in plants. In particular, the analysis of 107 quantitative traits reported in the first genome-wide association mapping study in Arabidopsis thaliana sets the stage for an exciting time in our understanding of plant adaptation. We also argue for the need to place phenotype-genotype association studies in an ecological context if one is to predict the evolutionary trajectories of plant species.

Natural allelic variation underlying a major fitness trade-off in Arabidopsis thaliana.

Plants can defend themselves against a wide array of enemies, from microbes to large animals, yet there is great variability in the effectiveness of such defences, both within and between species. Some of this variation can be explained by conflicting pressures from pathogens with different modes of attack. A second explanation comes from an evolutionary 'tug of war', in which pathogens adapt to evade detection, until the plant has evolved new recognition capabilities for pathogen invasion. If selection is, however, sufficiently strong, susceptible hosts should remain rare. That this is not the case is best explained by costs incurred from constitutive defences in a pest-free environment. Using a combination of forward genetics and genome-wide association analyses, we demonstrate that allelic diversity at a single locus, ACCELERATED CELL DEATH 6 (ACD6), underpins marked pleiotropic differences in both vegetative growth and resistance to microbial infection and herbivory among natural Arabidopsis thaliana strains. A hyperactive ACD6 allele, compared to the reference allele, strongly enhances resistance to a broad range of pathogens from different phyla, but at the same time slows the production of new leaves and greatly reduces the biomass of mature leaves. This allele segregates at intermediate frequency both throughout the worldwide range of A. thaliana and within local populations, consistent with this allele providing substantial fitness benefits despite its marked impact on growth.

Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines.

Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases, genome-wide association (GWA) studies have, owing to advances in genotyping and sequencing technology, become an obvious general approach for studying the genetics of natural variation and traits of agricultural importance. They are particularly useful when inbred lines are available, because once these lines have been genotyped they can be phenotyped multiple times, making it possible (as well as extremely cost effective) to study many different traits in many different environments, while replicating the phenotypic measurements to reduce environmental noise. Here we demonstrate the power of this approach by carrying out a GWA study of 107 phenotypes in Arabidopsis thaliana, a widely distributed, predominantly self-fertilizing model plant known to harbour considerable genetic variation for many adaptively important traits. Our results are dramatically different from those of human GWA studies, in that we identify many common alleles of major effect, but they are also, in many cases, harder to interpret because confounding by complex genetics and population structure make it difficult to distinguish true associations from false. However, a-priori candidates are significantly over-represented among these associations as well, making many of them excellent candidates for follow-up experiments. Our study demonstrates the feasibility of GWA studies in A. thaliana and suggests that the approach will be appropriate for many other organisms.

An atypical kinase under balancing selection confers broad-spectrum disease resistance in Arabidopsis.

The failure of gene-for-gene resistance traits to provide durable and broad-spectrum resistance in an agricultural context has led to the search for genes underlying quantitative resistance in plants. Such genes have been identified in only a few cases, all for fungal or nematode resistance, and encode diverse molecular functions. However, an understanding of the molecular mechanisms of quantitative resistance variation to other enemies and the associated evolutionary forces shaping this variation remain largely unknown. We report the identification, map-based cloning and functional validation of QRX3 (RKS1, Resistance related KinaSe 1), conferring broad-spectrum resistance to Xanthomonas campestris (Xc), a devastating worldwide bacterial vascular pathogen of crucifers. RKS1 encodes an atypical kinase that mediates a quantitative resistance mechanism in plants by restricting bacterial spread from the infection site. Nested Genome-Wide Association mapping revealed a major locus corresponding to an allelic series at RKS1 at the species level. An association between variation in resistance and RKS1 transcription was found using various transgenic lines as well as in natural accessions, suggesting that regulation of RKS1 expression is a major component of quantitative resistance to Xc. The co-existence of long lived RKS1 haplotypes in A. thaliana is shared with a variety of genes involved in pathogen recognition, suggesting common selective pressures. The identification of RKS1 constitutes a starting point for deciphering the mechanisms underlying broad spectrum quantitative disease resistance that is effective against a devastating and vascular crop pathogen. Because putative RKS1 orthologous have been found in other Brassica species, RKS1 provides an exciting opportunity for plant breeders to improve resistance to black rot in crops.

The role of glucosinolates and the jasmonic acid pathway in resistance of Arabidopsis thaliana against molluscan herbivores.

Although slugs and snails play important roles in terrestrial ecosystems and cause considerable damage on a variety of crop plants, knowledge about the mechanisms of plant immunity to molluscs is limited. We found slugs to be natural herbivores of Arabidopsis thaliana and therefore investigated possible resistance mechanisms of this species against several molluscan herbivores. Treating wounded leaves with the mucus residue ('slime trail') of the Spanish slug Arion lusitanicus increased wound-induced jasmonate levels, suggesting the presence of defence elicitors in the mucus. Plants deficient in jasmonate biosynthesis and signalling suffered more damage by molluscan herbivores in the laboratory and in the field, demonstrating that JA-mediated defences protect A. thaliana against slugs and snails. Furthermore, experiments using A. thaliana mutants with altered levels of specific glucosinolate classes revealed the importance of aliphatic glucosinolates in defending leaves and reproductive structures against molluscs. The presence in mollusc faeces of known and novel metabolites arising from glutathione conjugation with glucosinolate hydrolysis products suggests that molluscan herbivores actively detoxify glucosinolates. Higher levels of aliphatic glucosinolates were found in plants during the night compared to the day, which correlated well with the nocturnal activity rhythms of slugs and snails. Our data highlight the function of well-known antiherbivore defence pathways in resistance against slugs and snails and suggest an important role for the diurnal regulation of defence metabolites against nocturnal molluscan herbivores.