Metabolomic and genome-wide association studies drive genetic dissection and gene mining in tea plant

Abstract

As a vital reproductive organ, flowers significantly influence the yield, sensory quality, and breeding efficacy of Camellia sinensis. Many biological characteristics of tea plants are influenced by metabolites; however, our knowledge of metabolites in tea flowers is limited. To investigate the physiological basis and molecular mechanisms underlying tea flower metabolism, we integrated metabolomics and genome-wide association studies (GWAS) to analyze the metabolites present in the flowers of 171 tea genotypes. Untargeted metabolomic analysis detected 581 and 295 metabolites in positive and negative ionization modes, respectively. Twenty-seven distinct metabolites were observed between C. sinensis var. assamica (CSA) and C. sinensis var. sinensis (CSS). GWAS identified 1238 quantitative trait loci (QTL) associated with 505 metabolites. Some structurally related metabolites tended to share common QTL. Integrating GWAS findings with secondary mass spectrometry (MS/MS) fragmentation and haplotype analysis for metabolites (-)-epigallocatechin-3-(3"-O-methyl) gallate (EGCG-3''-O-ME), (-)-Epicatechin-3-(3''-O-methyl) gallate (ECG-3''-O-ME), Pos_1118, and Neg_365 (p-coumaroylquinic acid) resulted in the identification of three candidate genes ( W07g015551, W08g018636, and W01g002625). Taken together, our findings provide a foundation for exploring comprehensive metabolic pathways in various tissues of C. sinensis.

Keywords

Camellia sinensis / Flower / Metabolite / Genome-wide association study / Candidate genes

References

[1]	Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Hum Genet. 2018;103:338–48. https://doi.org/10.1016/j.ajhg.2018.07.015.
[2]	Chen W, Wang W, Peng M, Gong L, Gao Y, Wan J, et al. Comparative and parallel genome-wide association studies for metabolic and agronomic traits in cereals. Nat Commun. 2016;7:12767. https://doi.org/10.1038/ncomms12767.
[3]	Chen Y, Zhou Y, Zeng L, Dong F, Tu Y, Yang Z. Occurrence of functional molecules in the flowers of tea (Camellia sinensis) plants: evidence for a second resource. Molecules. 2018;23:790. https://doi.org/10.3390/molecules23040790.
[4]	Cui J, Zhou J, Du W, Guo D, Tang X, Zhao W, et al. Distribution of and temporal variation in volatiles in tea (Camellia sinensis) flowers during the opening stages. J Agric Food Chem. 2023;71:19682–93. https://doi.org/10.1021/acs.jafc.3c02690.
[5]	Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and vcftools. Bioinformatics. 2011;27:2156–8. https://doi.org/10.1093/bioinformatics/btr330.
[6]	Dong NQ, Sun Y, Guo T, Shi CL, Zhang YM, Kan Y, et al. UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice. Nat Commun. 2020;11:2629. https://doi.org/10.1038/s41467-020-16403-5.
[7]	Fang K, Xia Z, Li H, Jiang X, Qin D, Wang Q, et al. Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites. Hortic Res. 2021;8:42. https://doi.org/10.1038/s41438-021-00477-3.
[8]	He D, Singh SK, Peng L, Kaushal R, Vílchez JI, Shao C, et al. Flavonoid-attracted Aeromonas sp. from the Arabidopsis root microbiome enhances plant dehydration resistance. Isme J. 2022;16:2622–32. https://doi.org/10.1038/s41396-022-01288-7.
[9]	Hu Y, Zhang M, Lu M, Wu Y, Jing T, Zhao M, et al. Salicylic acid carboxyl glucosyltransferase ugt87e7 regulates disease resistance in camellia sinensis. Plant Physiol. 2022;188:1507–20. https://doi.org/10.1093/plphys/kiab569.
[10]	Huang R, Wang JY, Yao MZ, Ma CL, Chen L. Quantitative trait loci mapping for free amino acid content using an albino population and SNP markers provides insight into the genetic improvement of tea plants. Hortic Res. 2022;9:uhab29. https://doi.org/10.1093/hr/uhab029.
[11]	Jia S, Wang Y, Hu J, Ding Z, Liang Q, Zhang Y, et al. Mineral and metabolic profiles in tea leaves and flowers during flower development. Plant Physiol Biochem. 2016;106:316–26. https://doi.org/10.1016/j.plaphy.2016.06.013.
[12]	Jia X, Zhang W, Fernie AR, Wen W. Camellia sinensis (tea). Trends Genet. 2020;37:201–2. https://doi.org/10.1016/j.tig.2020.10.002.
[13]	Jiang X, Zhang W, Fernie AR, Wen W. Combining novel technologies with interdisciplinary basic research to enhance horticultural crops. Plant J. 2021;109:35–46. https://doi.org/10.1111/tpj.15553.
[14]	Jin J, Zhao M, Jing T, Wang J, Lu M, Pan Y, et al. (Z)-3-Hexenol integrates drought and cold stress signaling by activating abscisic acid glucosylation in tea plants. Plant Physiol. 2023a;193:1491–507. https://doi.org/10.1093/plphys/kiad346.
[15]	Jin JQ, Qu FR, Huang H, Liu QS, Wei MY, Zhou Y, et al. Characterization of two O-methyltransferases involved in the biosynthesis of O-methylated catechins in tea plant. Nat Commun. 2023b;14:5075. https://doi.org/10.1038/s41467-023-40868-9.
[16]	Jing T, Du W, Gao T, Wu Y, Zhang N, Zhao M. Herbivore-induced dmnt catalyzed by cyp82d47 plays an important role in the induction of ja-dependent herbivore resistance of neighboring tea plants. Plant Cell Environ. 2021;44:1178–91. https://doi.org/10.1111/pce.13861.
[17]	Kang HM, Sul JH, Service SK, Zaitlen NA, Kong SY, Freimer NB, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42:348–54. https://doi.org/10.1038/ng.548.
[18]	Kong W, Jiang M, Wang Y, Chen S, Zhang S, Lei W, et al. Pan-transcriptome assembly combined with multiple association analysis provides new insights into the regulatory network of specialized metabolites in the tea plant camellia sinensis. Hortic Res. 2022;9:uhac100. https://doi.org/10.1093/hr/uhac100.
[19]	Li MX, Yeung JM, Cherny SS, Sham PC. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum Genet. 2012;131:747–56. https://doi.org/10.1007/s00439-011-1118-2.
[20]	Li CF, Ma JQ, Huang DJ, Ma CL, Jin JQ, Yao MZ, et al. Comprehensive dissection of metabolic changes in albino and green tea cultivars. J Agric Food Chem. 2018;66:2040–8. https://doi.org/10.1021/acs.jafc.7b05623.
[21]	Li Y, Chen Y, Zhou L, You S, Deng H, Chen Y, et al. Microtom metabolic network: rewiring tomato metabolic regulatory network throughout the growth cycle. Mol Plant. 2020;13:1203–18. https://doi.org/10.1016/j.molp.2020.06.005.
[22]	Li HL, Wu L, Dong Z, Jiang Y, Jiang S, Xing H, et al. Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway. Hortic Res. 2021;8:189. https://doi.org/10.1038/s41438-021-00627-7.
[23]	Li X, Shi Z, Gao J, Wang X, Guo K. Candihap: a haplotype analysis toolkit for natural variation study. Mol Breeding. 2023;43:21. https://doi.org/10.1007/s11032-023-01366-4.
[24]	Liang X, Liu S, Wang T, Li F, Cheng J, Lai J, et al. Metabolomics-driven gene mining and genetic improvement of tolerance to salt-induced osmotic stress in maize. New Phytol. 2021;230:2355–70. https://doi.org/10.1111/nph.17323.
[25]	Liu Y, Hao X, Lu Q, Zhang W, Zhang H, Wang L, et al. Genome-wide identification and expression analysis of flowering-related genes reveal putative floral induction and differentiation mechanisms in tea plant (camellia sinensis). Genomics. 2020;112:2318–26. https://doi.org/10.1016/j.ygeno.2020.01.003.
[26]	Liu G, Huang K, Ke J, Chen C, Bao GH, Wan X. Novel camellia sinensis O-methyltransferase regulated by CsMADSL1 specifically methylates EGCG in cultivar “GZMe4.” J Agric Food Chem. 2023;71:6706–16. https://doi.org/10.1021/acs.jafc.2c06031.
[27]	Qin D, Wang Q, Li H, Jiang X, Fang K, Wang Q, et al. Identification of key metabolites based on non-targeted metabolomics and chemometrics analyses provides insights into bitterness in Kucha [camellia kucha (Chang et Wang) Chang]. Food Res Int. 2020;138:109789. https://doi.org/10.1016/j.foodres.2020.109789.
[28]	Qiu H, Zhu X, Wan H, Xu L, Zhang Q, Hou P, et al. Parallel metabolomic and transcriptomic analysis reveals key factors for quality improvement of tea plants. J Agric Food Chem. 2020;68:5483–95. https://doi.org/10.1021/acs.jafc.0c00434.
[29]	Ren N, Chen L, Li B, Rankin GO, Chen YC, Tu Y. Purified tea (camellia sinensis (L.) Kuntze) flower saponins induce the p53-dependent intrinsic apoptosis of cisplatin-resistant ovarian cancer cells. Int J Mol Sci. 2020;21:4324. https://doi.org/10.3390/ijms21124324.
[30]	Visscher PM, Wray NR, Zhang Q, Sklar P, McCarthy MI, Brown MA, et al. 10 years of gwas discovery: Biology function and translation. Am J Hum Genet. 2017;101:5–22. https://doi.org/10.1016/j.ajhg.2017.06.005.
[31]	Wan H, Liu H, Zhang J, Lyu Y, Li Z, He Y, et al. Lipidomic and transcriptomic analysis reveals reallocation of carbon flux from cuticular wax into plastid membrane lipids in a glossy “Newhall” navel orange mutant. Hortic Res. 2020a;7:41. https://doi.org/10.1038/s41438-020-0262-z.
[32]	Wan X, Wu S, Li Z, An X, Tian Y. Lipid metabolism: critical roles in male fertility and other aspects of reproductive development in plants. Mol Plant. 2020b;13:955–83. https://doi.org/10.1016/j.molp.2020.05.009.
[33]	Wang Y, Ren N, Rankin GO, Li B, Rojanasakul Y, Tu Y, et al. Anti-proliferative effect and cell cycle arrest induced by saponins extracted from tea (camellia sinensis) flower in human ovarian cancer cells. Journal of Functional Foods. 2017;37:310–21. https://doi.org/10.1016/j.jff.2017.08.001.
[34]	Wang H, Li S, Li Y, Xu Y, Wang Y, Zhang R, et al. MED25 connects enhancer–promoter looping and MYC2-dependent activation of jasmonate signalling. Nature Plants. 2019;5:616–25. https://doi.org/10.1038/s41477-019-0441-9.
[35]	Wang Q, Jiang X, Qin D, Liu S, Li H, Fang K, et al. Metabolic profiling of flavor compounds in black teas with almond odor during processing. Eur Food Res Technol. 2020a;246:2039–53. https://doi.org/10.1007/s00217-020-03554-4.
[36]	Wang X, Feng H, Chang Y, Ma C, Wang L, Hao X, et al. Population sequencing enhances understanding of tea plant evolution. Nat Commun. 2020b;11:4447. https://doi.org/10.1038/s41467-020-18228-8.
[37]	Way TD, Lin HY, Hua KT, Lee JC, Li WH, Lee MR, et al. Beneficial effects of different tea flowers against human breast cancer MCF-7 cells. Food Chem. 2009;114:1231–6. https://doi.org/10.1016/j.foodchem.2008.10.084.
[38]	Wen W, Li D, Li X, Gao Y, Li W, Li H, et al. Metabolome-based genome-wide association study of maize kernel leads to novel biochemical insights. Nat Commun. 2014;5:3438. https://doi.org/10.1038/ncomms4438.
[39]	Xia E, Tong W, Hou Y, An Y, Chen L, Wu Q, et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into its genome evolution and adaptation. Mol Plant. 2020;13:1013–26. https://doi.org/10.1016/j.molp.2020.04.010.
[40]	Ying S, Su M, Wu Y, Zhou L, Fu R, Li Y, et al. Trichome regulator slmixta-like directly manipulates primary metabolism in tomato fruit. Plant Biotechnol J. 2020;18:354–63. https://doi.org/10.1111/pbi.13202.
[41]	Yu A, Zhao J, Wang Z, Cheng K, Zhang P, Tian G, et al. Transcriptome and metabolite analysis reveal the drought tolerance of foxtail millet significantly correlated with phenylpropanoids-related pathways during germination process under PEG stress. BMC Plant Biol. 2020a;20:274. https://doi.org/10.1186/s12870-020-02483-4.
[42]	Yu X, Xiao J, Chen S, Yu Y, Ma J, Lin Y, et al. Metabolite signatures of diverse camellia sinensis tea populations. Nat Commun. 2020b;11:5586. https://doi.org/10.1038/s41467-020-19441-1.
[43]	Zeng L, Zhou X, Liao Y, Yang Z. Roles of specialized metabolites in biological function and environmental adaptability of tea plant (camellia sinensis) as a metabolite studying model. J Adv Res. 2021;34:159–71. https://doi.org/10.1016/j.jare.2020.11.004.
[44]	Zhang W, Zhang Y, Qiu H, Guo Y, Wan H, Zhang X, et al. Genome assembly of wild tea tree dasz reveals pedigree and selection history of tea varieties. Nat Commun. 2020;11:3719. https://doi.org/10.1038/s41467-020-17498-6.
[45]	Zhang W, Luo C, Scossa F, Zhang Q, Usadel B, Fernie AR, et al. A phased genome based on single sperm sequencing reveals crossover pattern and complex relatedness in tea plants. Plant J. 2021a;105:197–208. https://doi.org/10.1111/tpj.15051.
[46]	Zhang X, Chen S, Shi L, Gong D, Zhang S, Zhao Q, et al. Haplotype-resolved genome assembly provides insights into evolutionary history of the tea plant camellia sinensis. Nat Gene. 2021b;53:1250–9. https://doi.org/10.1038/s41588-021-00895-y.
[47]	Zhao M, Jin J, Wang J, Gao T, Luo Y, Jing T, et al. Eugenol functions as a signal mediating cold and drought tolerance via UGT71A59-mediated glucosylation in tea plants. Plant J. 2022;109:1489–506. https://doi.org/10.1111/tpj.15647.
[48]	Zhou S, Kremling KA, Bandillo N, Richter A, Zhang YK, Ahern KR, et al. Metabolome-scale genome-wide association studies reveal chemical diversity and genetic control of maize specialized metabolites. Plant Cell. 2019;31:937–55. https://doi.org/10.1105/tpc.18.00772.
[49]	Zhou B, Wang Z, Yin P, Ma B, Ma C, Xu C, et al. Impact of prolonged withering on phenolic compounds and antioxidant capability in white tea using LC-MS-based metabolomics and HPLC analysis: comparison with green tea. Food Chem. 2022;368:130855. https://doi.org/10.1016/j.foodchem.2021.130855.
[50]	Zhuang J, Dai X, Zhu M, Zhang S, Dai Q, Jiang X, et al. Evaluation of astringent taste of green tea through mass spectrometry-based targeted metabolic profiling of polyphenols. Food Chem. 2020;305:125507. https://doi.org/10.1016/j.foodchem.2019.125507.

Funding

National Natural Science Foundation of China NSFC-DFG collaborative project(3211101118); Fundamental Research Funds for the Central Universities(2662023PY011); Huazhong Agricultural University(SZYJY2021004); Guangdong Science and Technology Project(2019B030316026)