Anne Genissel
Johann Confais
Marc-Henri Lebrun- Unité Mixte de Recherche Bioger, Institut National de la Recherche Agronomique, AgroParisTech Institut des Sciences et Industries du Vivant et de L'environnement, Université Paris-Saclay, Thiverval-Grignon, France
Introduction
One of the main goal in phytopathology is to better understand the molecular basis of plant–pathogen co-evolution through the identification of effectors and effector targets that play a role in natural phenotypic variation. Fortunately, next generation sequencing (NGS)—which can measure genetic variation at hundreds of thousands of markers across a genome, including for non-model organisms—is now helping to reach this goal. Among all possible strategies using NGS data, we expect that genome-wide association studies (GWAS) have the most potential to revolutionize the field of phytopathology. In contrast to QTL mapping, GWAS use outbred populations to capture the standing genetic variation, thus characterizing the raw material for evolution. By examining the natural phenotypic and genetic variation, association mapping can elucidate the genetic basis underlying complex traits. In the two decades since association mapping successfully detected common variants for human complex diseases (Risch and Merikangas, 1996) and with the publication of the first successful GWAS in humans in 2005 (Klein et al., 2005), the number of published GWAS keeps increasing. Researchers in the field of plant pathogens are now embarking on GWAS, with the promise to open new frontiers of research.
Association Mapping is Emerging in the Field of Phytopathology
The feasibility of association genetics (within candidate regions or genome-wide) for identifying traits related to plant–pathogen interactions, has been demonstrated for various species including both host plants and pathogens. Paving the way, a few recent studies detected significant associations targeting resistance or (a)virulence loci. The first GWAS targeting phenotypes related to plant–pathogen interactions were done in the model species Arabidopsis (A.) thaliana. Even for this selfing species with an extent of linkage disequilibrium (LD) over a distance of 20 kb, Atwell et al. (2010) demonstrated the power of GWAS. A total of 107 phenotypes were tested using a high density SNP map for 191 natural accessions. This study revealed significant associations for disease resistance genes against the bacteria Pseudomonas (P.) syringae (genes known to mediate pathogenic bacteria recognition such as RPM1, RPS2, and RPS5), and confirmed the contribution of these loci to the natural variation of resistance of A. thaliana against P. synrigae. Subsequent studies reported more bacterial resistance loci in A. thaliana (Nemri et al., 2010; Huard-Chauveau et al., 2013; Debieu et al., 2016). In addition, many other resistance genes against pathogens have been identified in agriculturally important plant species such as in maize (Olukolu et al., 2014; Li et al., 2016), soybean (Chang et al., 2016), rice (Yoshida et al., 2009), wheat (Gurung et al., 2014), and barley (Turuspekov et al., 2016). Despite the promise of GWAS, until now there has been few association studies on the pathogen side, and in addition the approaches used to identify candidate loci were suboptimal [Xanthomonas campestris (Guy et al., 2013); Ralstonia solanacearum (Pensec et al., 2015); Heterobasidion annosum (Dalman et al., 2013); Fusarium graminearum (Talas et al., 2016); Parastagonospora nodorum (Gao et al., 2016)].
Overall, these findings demonstrate that GWAS is a successful method in the field of plant–pathogen interactions. In addition, when conducted with a sufficiently high-density marker map single-gene resolution of GWAS reaches far beyond the resolution of QTL mapping.
Association Mapping in Plant Pathogens: Lessons from Human GWAS
Based upon the premise that the contribution of a genetic variant to the phenotype is expressed, association mapping measures the significance of the association between the phenotype and the genotype. A detailed discussion of the statistical aspects of GWAS is beyond the scope of this article, a thorough discussion can be found in Balding (2006). Apart from the complicated statistics related to GWAS, we advice to read a recent review from Power et al. (2017), which considers in details many key elements on the methodology for GWAS. We believe that these considerations have to be considered carefully for plant pathogen GWAS as well. We argue that among them, sample size and LD are the most crucial.
Sample size must be considered because greater statistical power is obtained when increasing the sample size (Long and Langley, 1999). However, a suitable sample size to detect common variants with large phenotypic effect is not necessarily the same between species. In fact, different sample sizes can achieve the same power due to differences in ploidy level, recombination rate, and genetic diversity between the species under study. In addition, when sample size is not very large (i.e., less than few hundreds of individuals), rare variants with a minor allele frequency <5% should be excluded from the study because significant rare alleles represented by very few individuals are often spurious.
When carefully examined, LD greatly helps to improve the power of the study. Pairwise LD, as measured in association studies, corresponds to the non-random association of alleles at two loci. It is usually estimated from the coefficient correlation R2 and Lewontin's D′ (Gaut and Long, 2003; also see Mangin et al. (2012) for corrected R2 by considering the sample structure). The measure of LD is fundamental to predict the resolution of association studies. Lastly, to clarify how many significant markers are dragged along with the putatively causal marker, LD also should be well described in genomic regions covering top significant markers. The worst-case scenario is when there is a group of significant markers in “complete” LD. In that situation, complementary sequencing approaches using larger sample size may help to identify among associated markers which one(s) is causal. Complementary sequencing approaches also are beneficial when performing GWAS using RAD (Restriction site associated DNA)-sequencing that partially surveys the genome. Indeed if the LD drops quickly between adjacent (but potentially distant) markers, additional sequencing should be performed to call more variants in the region and identify the ones with the lowest associated P-values.
GWAS are now emerging in phytopathology; however, the vast heterogeneity of species that are under study and the lack of a uniform design for all systems remains a challenge. Thus, thorough analyses of GWAS are necessary to positively impact the reliability of this approach in the field of phytopathology.
Plant Pathogen GWAS: the Need for Considering Structural Variants
Association genetics is commonly performed on SNPs called from sequencing data. SNPs represent the highest frequency variants in genomes; nonetheless, insertions, deletions, and other structural variants (e.g., inversions, translocations within and between chromosomes) are not rare in micro-organisms and have been largely neglected so far. In our opinion, structural variants likely contribute to the potential of plant pathogen evolution and thus should be carefully examined. Hartmann et al. (2017) recently found that segmental deletion polymorphism within an effector gene is associated with virulence in the fungus Zymoseptoria tritici against wheat. Since structural variants can be identified using current genomic tools, we suggest structural variants need to be considered for GWAS.
Interpreting GWAS and the Search for Causality
GWAS studies yield a number of significant SNPs that must be interpreted cautiously regarding the inference of causality. For example, a candidate SNP can be thought of as DNA variation at a given position, which contributes to a proportion of the natural phenotypic variation in the population. It is a good candidate (i.e., putatively causal) because we observe the most significant difference of disease symptoms between the two alleles at that position, within a given region. For instance, for a haploid species with two alleles “A” and “B” at the site, the phenotypic value varies the most significantly between the two groups “A” vs. “B.” Once we have this result, however, additional work to elucidate the molecular function associated with the SNP needs to be conducted to support causality (i.e., the allelic variation at the SNP truly causes the variation of the phenotype). This step, the so-called “functional validation” is commonly feasible in micro-organisms. This resource should be better exploited to examine fully the molecular function hidden behind a top significant marker.
The functional approach can be optimized depending on the nature of the putative causal variant. Different strategies may be undertaken if the variant is regulatory or coding. For a SNP in a coding region, mutant lines can be produced from the wild type by knockout, and complement mutants can be produced to rescue the phenotype. If a variant is regulatory and within or nearby a gene, we need to test whether it has an effect on the tempo and mode of expression of that gene, by looking at natural variation of the gene expression within the sample. When molecular tools are available to transform the species, a more tedious strategy can be attempted. Regardless of the nature of the variant, mutant lines using allele swapping (for bi- or polyallelic sites) can be constructed in several genetic backgrounds and phenotypic effects of each mutation can be quantified. We think this strategy should be prioritized for a better understanding on the molecular mechanisms associated with plant pathogen evolution. Using GWAS and reverse genetics one demonstration of the molecular co-evolutionary dynamics between host and pathogens was done (Huard-Chauveau et al., 2013; Wang et al., 2015). The identification of a pseudokinase (RKS1) for the resistance of A. thaliana against X. campestris pv. campestris by GWAS (Huard-Chauveau et al., 2013) was supported by complemented study from Wang et al. (2015) who elucidated the role of that pseudokinase in AvrAC effector recognition.
The use of different, interconnected approaches is clearly an advantage to validate the GWAS results. For example, QTL mapping with recombinant populations can be used to confirm a region (Nemri et al., 2010). This strategy is more applicable when recombinant populations with high-density marker map are available. Alternatively, independent association genetics validation is also a relevant strategy. Often neglected except for human GWA, replication of association mapping greatly contributes to elucidating causality, with the advantage of examining the natural genetic variation segregating within populations that have different recombination history. There is indeed a strong hint for causality if two independent association studies reveal the same top variant, even though the two replicate studies use different population samples.
Conclusion
Examining the genotype-phenotype link may lead to an attractive list of genes useful for molecular studies in plant pathology and resistance management strategies. Yet, association genetics that aims at finding causal variants also provides the opportunity to address fundamental questions in evolutionary quantitative genetics. These questions are the following: (1) What molecular variants contribute the most to the standing phenotypic variation? (2) What is their phenotypic effect? (3) Are they regulatory or coding? (4) Are there many loci in action with small effects? (5) What molecular functions are associated with the phenotypic differences among individuals? (6) What evolutionary forces are maintaining this variation?
Importantly, GWA has limitations that may vary depending on the species and the phenotype under study and it is not trivial to perform a successful GWAS with good statistical power to detect causal variants. Even if molecular variants with large phenotypic effect are often discovered, certainly many more variants that contribute to the phenotypic variation remain to be found. The ability to detect these variants will require alternative strategies. We suggest that using very large sample sizes and combining population genomics with population transcriptomics and epigenetics should greatly contribute to elucidating the complex genotype-phenotype map in plant pathogen interactions.
Author Contributions
AG conceived and wrote the article with significant intellectual input from all authors: JC, ML, and LG.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors are thankful to Jana U'Ren for helpful suggestions on this manuscript.
References
Atwell, S., Huang, Y. S., Vilhjálmsson, B. J., Willems, G., Horton, M., Li, Y., et al. (2010). Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines. Nature 465, 627–631. doi: 10.1038/nature08800
Balding, D. J. (2006). A tutorial on statistical methods for population association studies. Nat. Rev. Genet. 7, 781–791. doi: 10.1038/nrg1916
Chang, H.-X., Brown, P. J., Lipka, A. E., Domier, L. L., and Hartman, G. L. (2016). Genome-wide association and genomic prediction identifies associated loci and predicts the sensitivity of Tobacco ringspot virus in soybean plant introductions. BMC Genomics 17:153. doi: 10.1186/s12864-016-2487-7
Dalman, K., Himmelstrand, K., Olson, A., Lind, M., Brandstrom-Durling, M., and Stenlid, J. (2013). A genome-wide association study identifies genomic regions for virulence in the non-model organism Heterobasidion annosum s.s. PLoS ONE 8:e53525. doi: 10.1371/journal.pone.0053525
Debieu, M., Huard-Chauveau, C., Genissel, A., Roux, F., and Roby, D. (2016). Quantitative disease resistance to the bacterial pathogen Xanthomonas campestris involves an Arabidopsis immune receptor pair and a gene of unknown function. Mol. Plant Pathol. 17, 510–520. doi: 10.1111/mpp.12298
Gao, Y., Liu, Z., Faris, J. D., Richards, J., Brueggeman, R. S., Li, X., et al. (2016). Validation of genome-wide association studies as a tool to identify virulence factors in Parastagonospora nodorum. Phytopathology 106, 1177–1185. doi: 10.1094/PHYTO-02-16-0113-FI
Gaut, B. S., and Long, A. D. (2003). The lowdown on linkage disequilibrium. Plant Cell Online 15, 1502–1506. doi: 10.1105/tpc.150730
Gurung, S., Mamidi, S., Bonman, J. M., Xiong, M., Brown-Guedira, G., and Adhikari, T. B. (2014). Genome-wide association study reveals novel quantitative trait loci associated with resistance to multiple leaf spot diseases of spring wheat. PLoS ONE 9:e108179. doi: 10.1371/journal.pone.0108179
Guy, E., Genissel, A., Hajri, A., Chabannes, M., David, P., Carrere, S., et al. (2013). Natural genetic variation of Xanthomonas campestris pv. campestris pathogenicity on Arabidopsis revealed by association and reverse genetics. mBio 4, e00538–e00512. doi: 10.1128/mBio.00538-12
Hartmann, F. E., Sanchez-Vallet, A., McDonald, B., and Croll, D. (2017). A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements. ISME J. 11, 1189–1204. doi: 10.1038/ismej.2016.196
Huard-Chauveau, C., Perchepied, L., Debieu, M., Rivas, S., Kroj, T., Kars, I., et al. (2013). An atypical kinase under balancing selection confers broad-spectrum disease resistance in Arabidopsis. PLoS Genet. 9:e1003766. doi: 10.1371/journal.pgen.1003766
Klein, R. J., Zeiss, C., Chew, E. Y., Tsai, J.-Y., Sackler, R. S., Haynes, C., et al. (2005). Complement factor H polymorphism in age-related macular degeneration. Science 308, 385–389. doi: 10.1126/science.1109557
Li, X., Zhou, Z., Ding, J., Wu, Y., Zhou, B., Wang, R., et al. (2016). Combined linkage and association mapping reveals QTL and candidate genes for plant and ear height in maize. Front. Plant Sci. 7:833. doi: 10.3389/fpls.2016.00833
Long, A. D., and Langley, C. (1999). The power of association studies to detect the contribution of candidate genetic loci to variation in complex traits. Genome Res. 9, 720–731. doi: 10.1101/gr.9.8.720
Mangin, B., Siberchicot, A., Nicolas, S., Doligez, A., This, P., and Cierco-Ayrolles, C. (2012). Novel measures of linkage disequilibrium that correct the bias due to population structure and relatedness. Heredity 108, 285–291. doi: 10.1038/hdy.2011.73
Nemri, A., Atwell, S., Tarone, A. M., Huang, Y. S., Zhao, K., Studholme, D. J., et al. (2010). Genome-wide survey of Arabidopsis natural variation in downy mildew resistance using combined association and linkage mapping. Proc. Natl. Acad. Sci. U.S.A. 107, 10302–10307. doi: 10.1073/pnas.0913160107
Olukolu, B. A., Wang, G.-F., Vontimitta, V., Venkata, B. P., Marla, S., Ji, J., et al. (2014). A genome-wide association study of the maize hypersensitive defense response identifies genes that cluster in related pathways. PLoS Genet. 10:e1004562. doi: 10.1371/journal.pgen.1004562
Pensec, F., Lebeau, A., Daunay, M. C., Chiroleu, F., Guidot, A., and Wicker, E. (2015). Towards the identification of type III effectors associated with Ralstonia solanacearum virulence on tomato and eggplant. Phytopathology 105, 1529–1544. doi: 10.1094/PHYTO-06-15-0140-R
Power, R. A., Parkhill, J., and de Oliveira, T. (2017). Microbial genome-wide association studies: lessons from human GWAS. Nat. Rev. Genet. 18, 41–50. doi: 10.1038/nrg.2016.132
Risch, N., and Merikangas, K. (1996). The future of genetic studies of complex human diseases. Science 273, 1516–1517. doi: 10.1126/science.273.5281.1516
Talas, F., Kalih, R., Miedaner, T., and McDonald, B. A. (2016). Genome-wide association study identifies novel candidate genes for aggressiveness, deoxynivalenol production, and azole sensitivity in natural field populations of Fusarium graminearum. Mol. Plant Microbe Interact. 29, 417–430. doi: 10.1094/MPMI-09-15-0218-R
Turuspekov, Y., Ormanbekova, D., Rsaliev, A., and Abugalieva, S. (2016). Genome-wide association study on stem rust resistance in Kazakh sring barley lines. BMC Plant Biol. 16:6. doi: 10.1186/s12870-015-0686-z
Wang, G., Roux, B., Feng, F., Guy, E., Li, L., Li, N., et al. (2015). The decoy substrate of a pathogen effector and a pseudokinase specify pathogen-induced modified-self recognition and immunity in plants. Cell Host Microbe 18, 285–295. doi: 10.1016/j.chom.2015.08.004
Keywords: plant–pathogen interaction, natural variation, GWAS, phenotype-genotype, linkage disequilibrium, functional validation
Citation: Genissel A, Confais J, Lebrun M-H and Gout L (2017) Association Genetics in Plant Pathogens: Minding the Gap between the Natural Variation and the Molecular Function. Front. Plant Sci. 8:1301. doi: 10.3389/fpls.2017.01301
Received: 01 February 2017; Accepted: 11 July 2017;
Published: 25 July 2017.
Edited by:
Richard Berthomé, Institut National de la Recherche Agronomique de Toulouse, FranceReviewed by:
Daniel J. Kliebenstein, University of California, Davis, United StatesPierre Pétriacq, University of Sheffield, United Kingdom
Copyright © 2017 Genissel, Confais, Lebrun and Gout. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Anne Genissel, anne.genissel@inra.fr