Interactions of Tomato and Botrytis cinerea Genetic Diversity: Parsing the Contributions of Host Differentiation, Domestication, and Pathogen Variation.

Although the impacts of crop domestication on specialist pathogens are well known, less is known about the interaction of crop variation and generalist pathogens. To study how genetic variation within a crop affects plant resistance to generalist pathogens, we infected a collection of wild and domesticated tomato accessions with a genetically diverse population of the generalist pathogen Botrytis cinerea We quantified variation in lesion size of 97 B. cinerea genotypes (isolates) on six domesticated tomato genotypes (Solanum lycopersicum) and six wild tomato genotypes (Solanum pimpinellifolium). Lesion size was significantly affected by large effects of the host and pathogen's genotype, with a much smaller contribution of domestication. This pathogen collection also enables genome-wide association mapping of B. cinerea Genome-wide association mapping of the pathogen showed that virulence is highly polygenic and involves a diversity of mechanisms. Breeding against this pathogen would likely require the use of diverse isolates to capture all possible mechanisms. Critically, we identified a subset of B. cinerea genes where allelic variation was linked to altered virulence against wild versus domesticated tomato, as well as loci that could handle both groups. This generalist pathogen already has a large collection of allelic variation that must be considered when designing a breeding program.

"Missing" G x E Variation Controls Flowering Time in Arabidopsis thaliana.

Understanding how genetic variation interacts with the environment is essential for understanding adaptation. In particular, the life cycle of plants is tightly coordinated with local environmental signals through complex interactions with the genetic variation (G x E). The mechanistic basis for G x E is almost completely unknown. We collected flowering time data for 173 natural inbred lines of Arabidopsis thaliana from Sweden under two growth temperatures (10°C and 16°C), and observed massive G x E variation. To identify the genetic polymorphisms underlying this variation, we conducted genome-wide scans using both SNPs and local variance components. The SNP-based scan identified several variants that had common effects in both environments, but found no trace of G x E effects, whereas the scan using local variance components found both. Furthermore, the G x E effects appears to be concentrated in a small fraction of the genome (0.5%). Our conclusion is that G x E effects in this study are mostly due to large numbers of allele or haplotypes at a small number of loci, many of which correspond to previously identified flowering time genes.

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.

Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis.

Transcriptional reprogramming forms a major part of a plant's response to pathogen infection. Many individual components and pathways operating during plant defense have been identified, but our knowledge of how these different components interact is still rudimentary. We generated a high-resolution time series of gene expression profiles from a single Arabidopsis thaliana leaf during infection by the necrotrophic fungal pathogen Botrytis cinerea. Approximately one-third of the Arabidopsis genome is differentially expressed during the first 48 h after infection, with the majority of changes in gene expression occurring before significant lesion development. We used computational tools to obtain a detailed chronology of the defense response against B. cinerea, highlighting the times at which signaling and metabolic processes change, and identify transcription factor families operating at different times after infection. Motif enrichment and network inference predicted regulatory interactions, and testing of one such prediction identified a role for TGA3 in defense against necrotrophic pathogens. These data provide an unprecedented level of detail about transcriptional changes during a defense response and are suited to systems biology analyses to generate predictive models of the gene regulatory networks mediating the Arabidopsis response to B. cinerea.

Genome-wide survey of Arabidopsis natural variation in downy mildew resistance using combined association and linkage mapping.

The model plant Arabidopsis thaliana exhibits extensive natural variation in resistance to parasites. Immunity is often conferred by resistance (R) genes that permit recognition of specific races of a disease. The number of such R genes and their distribution are poorly understood. In this study, we investigated the basis for resistance to the downy mildew agent Hyaloperonospora arabidopsidis ex parasitica (Hpa) in a global sample of A. thaliana. We implemented a combined genome-wide mapping of resistance using populations of recombinant inbred lines and a collection of wild A. thaliana accessions. We tested the interaction between 96 host genotypes collected worldwide and five strains of Hpa. Then, a fraction of the species-wide resistance was genetically dissected using six recently constructed populations of recombinant inbred lines. We found that resistance is usually governed by single dominant R genes that are concentrated in four genomic regions only. We show that association genetics of resistance to diseases such as downy mildew enables increased mapping resolution from quantitative trait loci interval to candidate gene level. Association patterns in quantitative trait loci intervals indicate that the pool of A. thaliana resistance sources against the tested Hpa isolates may be predominantly confined to six RPP (Resistance to Hpa) loci isolated in previous studies. Our results suggest that combining association and linkage mapping could accelerate resistance gene discovery in plants.

Quantitative interactions: the disease outcome of Botrytis cinerea across the plant kingdom.

Botrytis cinerea is a fungal pathogen that causes necrotic disease on more than a thousand known hosts widely spread across the plant kingdom. How B. cinerea interacts with such extensive host diversity remains largely unknown. To address this question, we generated an infectivity matrix of 98 strains of B. cinerea on 90 genotypes representing eight host plants. This experimental infectivity matrix revealed that the disease outcome is largely explained by variations in either the host resistance or pathogen virulence. However, the specific interactions between host and pathogen account for 16% of the disease outcome. Furthermore, the disease outcomes cluster among genotypes of a species but are independent of the relatedness between hosts. When analyzing the host specificity and virulence of B. cinerea, generalist strains are predominant. In this fungal necrotroph, specialization may happen by a loss in virulence on most hosts rather than an increase of virulence on a specific host. To uncover the genetic architecture of Botrytis host specificity and virulence, a genome-wide association study (GWAS) was performed and revealed up to 1492 genes of interest. The genetic architecture of these traits is widespread across the B. cinerea genome. The complexity of the disease outcome might be explained by hundreds of functionally diverse genes putatively involved in adjusting the infection to diverse hosts.

Pathogen Genetic Control of Transcriptome Variation in the Arabidopsis thaliana - Botrytis cinerea Pathosystem.

In plant-pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host-pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms' transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity.

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

The Glucosinolate Biosynthetic Gene AOP2 Mediates Feed-back Regulation of Jasmonic Acid Signaling in Arabidopsis.

Survival in changing and challenging environments requires an organism to efficiently obtain and use its resources. Due to their sessile nature, it is particularly critical for plants to dynamically optimize their metabolism. In plant primary metabolism, metabolic fine-tuning involves feed-back mechanisms whereby the output of a pathway controls its input to generate a precise and robust response to environmental changes. By contrast, few studies have addressed the potential for feed-back regulation of secondary metabolism. In Arabidopsis, accumulation of the defense compounds glucosinolates has previously been linked to genetic variation in the glucosinolate biosynthetic gene AOP2. AOP2 expression can increase the transcript levels of two known regulators (MYB28 and MYB29) of the pathway, suggesting that AOP2 plays a role in positive feed-back regulation controlling glucosinolate biosynthesis. We generated mutants affecting AOP2, MYB28/29, or both. Transcriptome analysis of these mutants identified a so far unrecognized link between AOP2 and jasmonic acid (JA) signaling independent of MYB28 and MYB29. Thus, AOP2 is part of a regulatory feed-back loop linking glucosinolate biosynthesis and JA signaling and thereby allows the glucosinolate pathway to influence JA sensitivity. The discovery of this regulatory feed-back loop provides insight into how plants optimize the use of resources for defensive metabolites.

Acetylation of cell wall is required for structural integrity of the leaf surface and exerts a global impact on plant stress responses.

The epidermis on leaves protects plants from pathogen invasion and provides a waterproof barrier. It consists of a layer of cells that is surrounded by thick cell walls, which are partially impregnated by highly hydrophobic cuticular components. We show that the Arabidopsis T-DNA insertion mutants of REDUCED WALL ACETYLATION 2 (rwa2), previously identified as having reduced O-acetylation of both pectins and hemicelluloses, exhibit pleiotrophic phenotype on the leaf surface. The cuticle layer appeared diffused and was significantly thicker and underneath cell wall layer was interspersed with electron-dense deposits. A large number of trichomes were collapsed and surface permeability of the leaves was enhanced in rwa2 as compared to the wild type. A massive reprogramming of the transcriptome was observed in rwa2 as compared to the wild type, including a coordinated up-regulation of genes involved in responses to abiotic stress, particularly detoxification of reactive oxygen species and defense against microbial pathogens (e.g., lipid transfer proteins, peroxidases). In accordance, peroxidase activities were found to be elevated in rwa2 as compared to the wild type. These results indicate that cell wall acetylation is essential for maintaining the structural integrity of leaf epidermis, and that reduction of cell wall acetylation leads to global stress responses in Arabidopsis.

Meta-analysis of metabolome QTLs in Arabidopsis: trying to estimate the network size controlling genetic variation of the metabolome.

A central goal of systems biology is to develop models that are both predictive and accurately describe the biological system. One complexity to this endeavor is that it is possible to develop models that appear predictive even if they use far fewer components than the biological system itself uses for the same process. This problem also occurs in quantitative genetics where it is often possible to describe the variation in a system using fewer genes than are actually variable due to the complications of linkage between causal polymorphisms and population structure. Thus, there is a crucial need to begin an empirical investigation into the true number of components that are used by biological systems to determine a phenotypic outcome. In this study, we use a meta-analysis of directly comparable metabolomics quantitative studies using quantitative trait locus mapping and genome wide association mapping to show that it is currently not possible to estimate how many genetic loci are truly polymorphic within Arabidopsis thaliana. Our analysis shows that it would require the analysis of at least a 1000 line bi-parental population to begin to estimate how many polymorphic loci control metabolic variation within Arabidopsis. Understanding the base number of loci that are actually involved in determining variation in metabolic systems is fundamental to developing systems models that are truly reflective of how metabolism is modulated within a living organism.

Exploiting natural variation of secondary metabolism identifies a gene controlling the glycosylation diversity of dihydroxybenzoic acids in Arabidopsis thaliana.

Plant secondary metabolism is an active research area because of the unique and important roles the specialized metabolites have in the interaction of plants with their biotic and abiotic environment, the diversity and complexity of the compounds and their importance to human medicine. Thousands of natural accessions of Arabidopsis thaliana characterized with increasing genomic precision are available, providing new opportunities to explore the biochemical and genetic mechanisms affecting variation in secondary metabolism within this model species. In this study, we focused on four aromatic metabolites that were differentially accumulated among 96 Arabidopsis natural accessions as revealed by leaf metabolic profiling. Using UV, mass spectrometry, and NMR data, we identified these four compounds as different dihydroxybenzoic acid (DHBA) glycosides, namely 2,5-dihydroxybenzoic acid (gentisic acid) 5-O-ß-D-glucoside, 2,3-dihydroxybenzoic acid 3-O-ß-D-glucoside, 2,5-dihydroxybenzoic acid 5-O-ß-D-xyloside, and 2,3-dihydroxybenzoic acid 3-O-ß-D-xyloside. Quantitative trait locus (QTL) mapping using recombinant inbred lines generated from C24 and Col-0 revealed a major-effect QTL controlling the relative proportion of xylosides vs. glucosides. Association mapping identified markers linked to a gene encoding a UDP glycosyltransferase gene. Analysis of Transfer DNA (T-DNA) knockout lines verified that this gene is required for DHBA xylosylation in planta and recombinant protein was able to xylosylate DHBA in vitro. This study demonstrates that exploiting natural variation of secondary metabolism is a powerful approach for gene function discovery.

Parallel Loss-of-Function at the RPM1 Bacterial Resistance Locus in Arabidopsis thaliana.

Dimorphism at the Resistance to Pseudomonas syringae pv. maculicola 1 (RPM1) locus is well documented in natural populations of Arabidopsis thaliana and has been portrayed as a long-term balanced polymorphism. The haplotype from resistant plants contains the RPM1 gene, which enables these plants to recognize at least two structurally unrelated bacterial effector proteins (AvrB and AvrRpm1) from bacterial crop pathogens. A complete deletion of the RPM1 coding sequence has been interpreted as a single event resulting in susceptibility in these individuals. Consequently, the ability to revert to resistance or for alternative R-gene specificities to evolve at this locus has also been lost in these individuals. Our survey of variation at the RPM1 locus in a large species-wide sample of A. thaliana has revealed four new loss-of-function alleles that contain most of the intervening sequence of the RPM1 open reading frame. Multiple loss-of-function alleles may have originated due to the reported intrinsic cost to plants expressing the RPM1 protein. The frequency and geographic distribution of rpm1 alleles observed in our survey indicate the parallel origin and maintenance of these loss-of-function mutations and reveal a more complex history of natural selection at this locus than previously thought.

Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel.

Arabidopsis thaliana is native to Eurasia and is naturalized across the world. Its ability to be easily propagated and its high phenotypic variability make it an ideal model system for functional, ecological and evolutionary genetics. To date, analyses of the natural genetic variation of A. thaliana have involved small numbers of individual plants or genetic markers. Here we genotype 1,307 worldwide accessions, including several regional samples, using a 250K SNP chip. This allowed us to produce a high-resolution description of the global pattern of genetic variation. We applied three complementary selection tests and identified new targets of selection. Further, we characterized the pattern of historical recombination in A. thaliana and observed an enrichment of hotspots in its intergenic regions and repetitive DNA, which is consistent with the pattern that is observed for humans but which is strikingly different from that observed in other plant species. We have made the seeds we used to produce this Regional Mapping (RegMap) panel publicly available. This panel comprises one of the largest genomic mapping resources currently available for global natural isolates of a non-human species.