Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites

Abstract

The characteristic secondary metabolites in tea (theanine, caffeine, and catechins) are important factors contributing to unique tea flavors. However, there has been relatively little research on molecular markers related to these metabolites. Thus, we conducted a genome-wide association analysis of the levels of these tea flavor-related metabolites in three seasons. The theanine, caffeine, and catechin levels in Population 1 comprising 191 tea plant germplasms were examined, which revealed that their heritability exceeded 0.5 in the analyzed seasons, with the following rank order (highest to lowest heritabilities): (+)-catechin > (−)-gallocatechin gallate > caffeine = (−)-epicatechin > (−)-epigallocatechin-3-gallate > theanine > (−)-epigallocatechin > (−)-epicatechin-3-gallate > catechin gallate > (+)-gallocatechin. The SNPs detected by amplified-fragment SNP and methylation sequencing divided Population 1 into three groups and seven subgroups. An association analysis yielded 307 SNP markers related to theanine, caffeine, and catechins that were common to all three seasons. Some of the markers were pleiotropic. The functional annotation of 180 key genes at the SNP loci revealed that FLS, UGT, MYB, and WD40 domain-containing proteins, as well as ATP-binding cassette transporters, may be important for catechin synthesis. KEGG and GO analyses indicated that these genes are associated with metabolic pathways and secondary metabolite biosynthesis. Moreover, in Population 2 (98 tea plant germplasm resources), 30 candidate SNPs were verified, including 17 SNPs that were significantly or extremely significantly associated with specific metabolite levels. These results will provide a foundation for future research on important flavor-related metabolites and may help accelerate the breeding of new tea varieties.

References

[1] Yamanishi, T. Flavor of tea. Food Rev. Int 11, 371-546 (1995).

[2] Chen, G. H., Yang, C. Y., Lee, S. J., Wu, C. C. & Tzen, J. T. C. Catechin content and the degree of its galloylation in oolong tea are inversely correlated with cultivation altitude. J. Food Drug Anal. 22, 303-309 (2014).

[3] Han, Z. et al. Data on green tea flavor determinants as affected by cultivars and manufacturing processes. Data Brief. 10, 492-498 (2017).

[4] Zhang, W. et al. The effects and mechanisms of epigallocatechin-3-gallate on reversing multidrug resistance in cancer. Trends Food Sci. Tech. 93, 221-233 (2019).

[5] Liang, Y. & Shi, M. Advances in tea plant genetics and breeding. J. Tea Sci. 35, 103-109 (2015).

[6] Ma, J. et al. Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis) . PLoS ONE 9, e93131 (2014).

[7] Xu, L. Y. et al. High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing . BMC Genomics 19, 955 (2018).

[8] Ma, J. et al. Quantitative trait loci mapping for theobromine and caffeine contents in tea plant (Camellia sinensis) . J. Agric Food Chem. 66, 13321-13327 (2018).

[9] Li, X., Ma, J., Yao, M. & Chen, L. SNP detection and mapping of theanine synthetase gene in tea plant. J. Tea Sci. 37, 251-257 (2017).

[10] Hazra, A., Kumar, R., Sengupta, C. & Das, S. Genome-wide SNP detection in Darjeeling tea: unravelling genetic structure, functional impact and trait associations. Preprint at https://doi.org/10.1101/20200902.280461 (2020).

[11] Wang, R., Gao, X., Yang, J. & Kong X. Genome-wide association study to identify favorable SNP allelic variations and candidate genes that control the timing of Spring bud flush of Tea (Camellia sinensis) using SLAF-seq . J Agr Food Chem. 67, 10380-10391 (2019).

[12] Zhang, W., Zhang, Y., Qiu, H., Guo, Y. & Wen, W. Genome assembly of wild tea tree dasz reveals pedigree and selection history of tea varieties. Nat Commun. 11, 3719 (2020).

[13] Jiang, C. et al. Identification and distribution of a single nucleotide polymorphism responsible for the catechin content in tea plants. Hortic. Res 7, 24 (2020).

[14] Gaut, B. S. & Long, A. D. The lowdown on linkage disequilibrium. Plant Cell 15, 1502-1506 (2003).

[15] Xia, Z., Zou, M., Zhang, S., Feng, B. & Wang, W. AFSM sequencing approach: a simple and rapid method for genome-wide SNP and methylation site discovery and genetic mapping. Sci. Rep. 4, 7300 (2014).

[16] Espigolan, R. et al. Study of whole genome linkage disequilibrium in Nellore cattle. BMC Genomics 14, 305 (2013).

[17] Huang, H., Yao, Q., Xia, E. & Gao, L. Metabolomics and transcriptomics analyses reveal nitrogen influences on the accumulation of flavonoids and amino acids in young shoots of tea plant (Camellia sinensis L.) associated with tea flavor . J. Agric Food Chem. 66, 9828-9838 (2018).

[18] Zhu, J. et al. Transcriptome analysis reveals the mechanism of fluoride treatment affecting biochemical components in Camellia sinensis . Int J. Mol. Sci. 20, 237 (2019).

[19] Wei, K. et al. A coupled role for CsMYB75 and CsGSTF1 in anthocyanin hyperaccumulation in purple tea. Plant J. 97, 825-840 (2019).

[20] Zhang, Y. et al. Identification of key genes involved in catechin metabolism in tea seedlings based on transcriptomic and HPLC analysis. Plant Physiol. Biochem 133, 107-115 (2018).

[21] Wu, L. Y. et al. Complementary iTRAQ proteomic and transcriptomic analyses of leaves in tea plant (Camellia sinensis L.) with different maturity and regulatory network of flavonoid biosynthesis . Front Plant Sci. 7, 332 (2019).

[22] Jiang, X. et al. Tissue-specific, development-dependent phenolic compounds accumulation profile and gene expression pattern in tea plant [ Camellia sinensis] . PLoS One 8, E62315 (2013).

[23] Song, L. et al. Molecular link between leaf coloration and gene expression of flavonoid and carotenoid biosynthesis in Camellia sinensis cultivar ‘Huangjinya’. Front Plant Sci. 8, 803 (2017).

[24] Yan, W. et al. Genome-wide identification of genes probably relevant to the uniqueness of tea plant (Camellia sinensis) and its cultivars . Int J. Genomics 2015, 527054 (2015).

[25] Liu, L. A. et al. Metabolite profiling and transcriptomic analyses reveal an essential role of UVR8-mediated signal transduction pathway in regulating flavonoid biosynthesis in tea plants (Camellia sinensis) in response to shading . BMC Plant Biol. 18, 233 (2018).

[26] Tai, Y. et al. Gene co-expression network analysis reveals coordinated regulation of three characteristic secondary biosynthetic pathways in tea plant (Camellia sinensis) . BMC Genomics 19, 616 (2018).

[27] Di, T. et al. Transcriptomic and metabolic insights into the distinctive effects of exogenous melatonin and gibberellin on terpenoid synthesis and plant hormone signal transduction pathway in Camellia sinensis . J. Agric Food Chem. 67, 4689-4699 (2018).

[28] Wang, W. W. et al. Transcriptome and metabolome analysis reveal candidate genes and biochemicals involved in tea geometrid defense in Camellia sinensis . PLoS One 13, e0201670 (2018).

[29] Zhao, J., Li, P., Xia, T. & Wan X. Exploring plant metabolic genomics: chemical diversity, metabolic complexity in the biosynthesis and transport of specialized metabolites with the tea plant as a model. Crit. Rev. Biotechnol. https://doi.org/10.1080/07388551.2020.1752617 (2020).

[30] Zhu, J. et al. Comprehensive co-expression analysis provides novel insights into temporal variation of flavonoids in fresh leaves of the tea plant (Camellia sinensis) . Plant Sci. 290, 110306 (2019).

[31] Fang, K., Jiang, X. & Wu, H. Research progress on theanine metabolism and its content breeding in tea. Acta Horticulturae Sin. 43, 1791-1802 (2016).

[32] Flint-Garcia, S. et al. Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 44, 1054-1064 (2006).

[33] Yao, M., Ma, C., Qiao, T., Jin, J. & Chen, L. Diversity distribution and population structure of tea germplasms in China revealed by EST-SSR markers. Tree Genet genomes 8, 205-220 (2012).

[34] Klein, R. J. et al. Complement factor H polymorphism in age-related macular degeneration. Science 308, 385-389 (2005).

[35] Luo, J. Metabolite-based genome-wide association studies in plants. Curr. Opin. Plant Biol. 24, 31-38 (2015).

[36] Li, C., Wang, Y., Ai, N., Li, Y. & Song, J. A genome-wide association study of early-maturation traits in upland cotton based on the CottonSNP80K array. J. Integr. Plant Biol. 60, 970-985 (2018).

[37] Lander, E. & Kruglyak, L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat. Genet 11, 241-247 (1995).

[38] Manolio, T. A. et al. Finding the missing heritability of complex diseases. Nature 461, 747-753 (2009).

[39] Myles, S. et al. Association mapping: critical considerations shift from genotyping to experimental design. Plant Cell 21, 2194-2202 (2009).

[40] Olsen, K. M. & Wendel, J. F. A bountiful harvest: genomic insights into crop domestication phenotypes. Annu Rev. Plant Biol. 64, 47-70 (2013).

[41] Korte, A. & Farlow, A. The advantages and limitations of trait analysis with GWAS: a review. Plant Methods 9, 29 (2013).

[42] Kato, M., Mizuno, K., Crozier, A., Fujimura, T. & Ashihara, H. Caffeine synthase gene from tea leaves. Nature 406, 956-957 (2000).

[43] Xia, E. et al. Tea plant genomics: achievements, challenges and perspectives. Horticulture Res. 7, 7 (2020).

[44] Shen, B. R. et al. Engineering a new chloroplastic photorespiratory bypass to increase photosynthetic efficiency and productivity in rice. Mol. Plant 12, 199-214 (2019).

[45] Li, R. et al. Comparative analysis of nitrogen spatial heterogeneity responses in low nitrogen susceptible and tolerant tea plants (Camellia sinensis) . Sci. Hortic.-Amst. 246, 182-189 (2019).

[46] Sun, M. et al. Metabolic flux enhancement and transcriptomic analysis displayed the changes of catechins following long-term pruning in tea trees (Camellia sinensis) . J. Agric Food Chem. 66, 8566-8573 (2018).

[47] Chen, Y. et al. Proteolysis of chloroplast proteins is responsible for accumulation of free amino acids in dark-treated tea (Camellia sinensis) leaves . J. Proteom. 157, 10-17 (2017).

[48] Liu, Y. et al. A WD40 repeat protein from Camellia sinensis regulates anthocyanin and proanthocyanidin accumulation through the formation of MYB-bHLH-WD40 ternary complexes. Int J. Mol. Sci. 19, 1686 (2018).

[49] Xia, E. H. et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol. Plant 10, 866-877 (2017).

[50] Zou, B. et al. Atmyb44 positively modulates disease resistance to pseudomonas syringae through the salicylic acid signalling pathway in arabidopsis. Funct. Plant Biol. 40, 304 (2013).

[51] Mohot, P. R. C. Tea-determination of tea polyphenols content. All China Federation of Supply and Marketing Cooperatives, China. 149-151 (2002).

[52] Wei, C. et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality . Proc. Natl Acad. Sci. USA 115, E4151-E4158 (2018).