Pan-transcriptome assembly combined with multiple association analysis provides new insights into the regulatory network of specialized metabolites in the tea plant Camellia sinensis

Abstract

Specialized metabolites not only play important roles in biotic and abiotic stress adaptation of tea plants (Camellia sinensis (L.) O. Kuntze) but also contribute to the unique flavor of tea, the most important nonalcoholic beverage. However, the molecular networks and major genes that regulate specialized metabolites in tea plants are not well understood. Here, we constructed a population-level pan-transcriptome of the tea plant leaf using second-leaf transcriptome data from 134 accessions to investigate global expression differences in the population, expression presence or absence variations (ePAVs), and differentially expressed genes (DEGs) between pure Camellia sinensis var. assamica (CSA) and pure Camellia sinensis var. sinensis (CSS) accessions. Next, we used a genome-wide association study, a quantitative trait transcript study, and a transcriptome-wide association study to integrate genotypes, accumulation levels of specialized metabolites, and expression levels of pan-transcriptome genes to identify candidate regulatory genes for flavor-related metabolites and to construct a regulatory network for specialized metabolites in tea plants. The pan-transcriptome contains 30 482 expressed genes, 4940 and 5506 of which were newly annotated from a de novo transcriptome assembly without a reference and a genome reference-based assembly, respectively. DEGs and ePAVs indicated that CSA and CSS were clearly differentiated at the population transcriptome level, and they were closely related to abiotic tolerance and secondary metabolite synthesis phenotypes of CSA and CSS based on gene annotations. The regulatory network contained 212 specialized metabolites, 3843 candidate genes, and 3407 eQTLs, highlighting many pleiotropic candidate genes, candidate gene-rich eQTLs, and potential regulators of specialized metabolites. These included important transcription factors in the AP2/ERF-ERF, MYB, WD40, and bHLH families. CsTGY14G0001296, an ortholog of AtANS, appeared to be directly related to variation in proanthocyanins in the tea plant population, and the CsTGY11G0002074 gene encoding F3^′5^′H was found to contribute to the biased distribution of catechins between pure CSAs and pure CSSs. Together, these results provide a new understanding of the metabolite diversity in tea plants and offer new insights for more effective breeding of better-flavored tea varieties.

References

[1] Yu XM, Xiao J, Chen S et al. Metabolite signatures of diverse Camellia sinensis tea populations. Nat Commun. 2020; 11: 5586.

[2] Huang H, Yao QY, Xia EH et al. 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. 2018; 66: 9828-38.

[3] Zhang SH, Zhang L, Tai Y et al. Gene discovery of characteristic metabolic pathways in the tea plant (Camellia sinensis) using ’Omics’-based network approaches: a future perspective . Front Plant Sci. 2018; 9: 480.

[4] Liao YY, Fu XM, Zeng LT et al. Strategies for studying in vivo biochemical formation pathways and multilevel distributions of quality or function-related specialized metabolites in tea (Camellia sinensis) . Crit Rev Food Sci Nutr. 2022; 62: 429-42.

[5] Fang KX, Xia Z, Li H et al. Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites. Hortic Res. 2021; 8: 42.

[6] Narukawa M, Toda Y, Nakagita T et al. L-Theanine elicits umami taste via the T1R1+T1R3 umami taste receptor. Amino Acids. 2014; 46: 1583-7.

[7] Scharbert S, Holzmann N, Hofmann T . Identification of the astringent taste compounds in black tea infusions by combining instrumental analysis and human bioresponse. J Agric Food Chem. 2004; 52: 3498-508.

[8] Yang ZY, Baldermann S, Watanabe N . Recent studies of the volatile compounds in tea. Food Res Int. 2013; 53: 585-99.

[9] Wei CL, Yang H, Wang S 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 U S A. 2018; 115: E4151-8.

[10] Xia EH, Zhang HB, Sheng J et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol Plant. 2017; 10: 866-77.

[11] Dixon RA, Pasinetti GM . Flavonoids and isoflavonoids: from plant biology to agriculture and neuroscience. Plant Physiol. 2010; 154: 453-7.

[12] Li CF, Zhu Y, Zhao QY et al. Global transcriptome and gene regulation network for secondary metabolite biosynthesis of tea plant (Camellia sinensis) . BMC Genomics. 2015; 16: 560.

[13] Suzuki T, Takahashi E . Biosynthesis of caffeine by tea-leaf extracts. Enzymic formation of theobromine from 7-methylxanthine and of caffeine from theobromine. The Biochemical Journal. 1975; 146: 87-96.

[14] Sasaoka K, Kito M, Onishi Y . Some properties of the theanine synthesizing enzyme in tea seedlings. Agric Biol Chem. 1965; 29: 984-8.

[15] Zhu B, Guo J, Dong C et al. CsAlaDC and CsTSI work coordinately to determine theanine biosynthesis in tea plants (Camellia sinensis L.) and confer high levels of L-theanine accumulation in a non-tea plant . Plant Biotechnol J. 2021; 19: 2395-7.

[16] Hazra A, Kumar R, Sengupta C et al. Genome-wide SNP discovery from Darjeeling tea cultivars-their functional impacts and application toward population structure and trait associations. Genomics. 2021; 113: 66-78.

[17] Zhang WY, Zhang Y, Qui H et al. Genome assembly of wild tea tree DASZ reveals pedigree and selection history of tea varieties. Nat Commun. 2020; 11: 3719.

[18] Liang XY, Liu S, Wang T et al. Metabolomics-driven gene mining and genetic improvement of tolerance to salt-induced osmotic stress in maize. New Phytol. 2021; 230: 2355-70.

[19] Chen W, Wang W, Peng M et al. Comparative and parallel genome-wide association studies for metabolic and agronomic traits in cereals. Nat Commun. 2016; 7: 12767.

[20] Chen W, Gao Y, Xie W et al. Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism. Nat Genet. 2014; 46: 714-21.

[21] Li K, Wang D, Gong L et al. Comparative analysis of metabolome of rice seeds at three developmental stages using a recombinant inbred line population. Plant J. 2019; 100: 908-22.

[22] Zhu GT, Wang S, Huang Z et al. Rewiring of the fruit metabolome in tomato breeding. Cell. 2018; 172: 249-61.

[23] Vernikos G, Medini D, Riley DR et al. Ten years of pan-genome analyses. Curr Opin Microbiol. 2015; 23: 148-54.

[24] Bayer PE, Golicz AA, Scheben A et al. Plant pan-genomes are the new reference. Nat Plants. 2020; 6: 914-20.

[25] LIU D, YANG L, ZHANG JZ et al. Domestication and breeding changed tomato fruit transcriptome. J Integr Agric. 2020; 19: 120-32.

[26] Hirsch CN, Foerster JM, Johnson JM et al. Insights into the maize pan-genome and pan-transcriptome. Plant Cell. 2014; 26: 121-35.

[27] Swanson-Wagner R, Briskine R, Schaefer R et al. Reshaping of the maize transcriptome by domestication. Proc Natl Acad Sci U S A. 2012; 109: 11878-83.

[28] Jin ML, Liu H, He C et al. Maize pan-transcriptome provides novel insights into genome complexity and quantitative trait variation. Sci Rep. 2016; 6: 18936.

[29] Won SY, Kwon SJ, Lee TH et al. Comparative transcriptome analysis reveals whole-genome duplications and gene selection patterns in cultivated and wild chrysanthemum species. Plant Mol Biol. 2017; 95: 451-61.

[30] Tettelin H, Masignani V, Cielewicz M. J. et al. genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome". Proc Natl Acad Sci U S A. 2005; 102: 13950-5.

[31] Lu F, Romay MC, Glaubitz JC et al. High-resolution genetic mapping of maize pan-genome sequence anchors. Nat Commun. 2015; 6: 6914.

[32] Qin P, Lu H, Du H et al. Pan-genome analysis of 33 genetically diverse rice accessions reveals hidden genomic variations. Cell. 2021; 184: 3542-58.

[33] Zhao Q, Feng Q, Lu H et al. Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice. Nat Genet. 2018; 50: 278-84.

[34] Liu YC, Du H, Li P et al. Pan-genome of wild and cultivated soybeans. Cell. 2020; 182: 162-76.

[35] Li YH, Zhou G, Ma J et al. De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits. Nat Biotechnol. 2014; 32: 1045-52.

[36] Lyu J . Crop genomics tomato pan-genome. Nat Plants. 2019; 5: 558-8.

[37] Tao YF, Luo H, Xu J et al. Extensive variation within the pan-genome of cultivated and wild sorghum. Nat Plants. 2021; 7: 766-73.

[38] Song JM, Guan Z, Hu J et al. Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus . Nat Plants. 2020; 6: 34-45.

[39] Song JM, Liu DX, Xie WZ et al. BnPIR: Brassica napus pan-genome information resource for 1689 accessions. Plant Biotechnol J. 2021; 19: 412-4.

[40] Petek M, Zagorscak M, Ramsak Z et al. Cultivar-specific transcriptome and pan-transcriptome reconstruction of tetraploid potato. Sci Data. 2020; 7: 249.

[41] Ma YL, Liu M, Stiller J et al. A pan-transcriptome analysis shows that disease resistance genes have undergone more selection pressure during barley domestication. BMC Genomics. 2019; 20: 12.

[42] Medina CA, Samac DA, Yu LX . Pan-transcriptome identifying master genes and regulation network in response to drought and salt stresses in alfalfa (Medicago sativa L.) . Sci Rep. 2021; 11: 17203.

[43] Zhang X, Chen S, Shi L et al. Haplotype-resolved genome assembly provides insights into evolutionary history of the tea plant Camellia sinensis . Nat Genet. 2021; 53: 1250-9.

[44] Grabherr MG, Haas BJ, Yassour M et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29: 644-52.

[45] Jiang XL, Huang K, Zheng G et al. CsMYB5a and CsMYB5e from Camellia sinensis differentially regulate anthocyanin and proanthocyanidin biosynthesis. Plant Sci. 2018; 270: 209-20.

[46] Wei K, Wang L, Zhang Y et al. A coupled role for CsMYB75 and CsGSTF1 in anthocyanin hyperaccumulation in purple tea . Plant J. 2019; 97: 825-40.

[47] Zheng XT, Chen YL, Zhang XH et al. ANS-deficient Arabidopsis is sensitive to high light due to impaired anthocyanin photoprotection. Funct Plant Biol. 2019; 46: 756-65.

[48] Qi Y, Zhao D, Lv L . Cloning of CsANS gene from tea plant (Camellia sinensis) and its functional analysis in transgenic tobacco (Nicotiana tabacum) . Journol of Agricultural Biotechnology. 2019; 27: 636-44.

[49] Mei Y, Xie H, Liu S et al. Metabolites and transcriptional profiling analysis reveal the molecular mechanisms of the anthocyanin metabolism in the "Zijuan" tea plant (Camellia sinensis var. assamica) . J Agric Food Chem. 2021; 69: 414-27.

[50] Zhu MZ, Zhou F, Ran LS et al. Metabolic profiling and gene expression analyses of purple-leaf formation in tea cultivars (Camellia sinensis var. sinensis and var. assamica) . Front Plant Sci. 2021; 12: 606962.

[51] Joshi R, Rana A, Kumar V et al. Anthocyanins enriched purple tea exhibits antioxidant, immunostimulatory and anticancer activities. J Food Sci Technol. 2017; 54: 1953-63.

[52] Mato M, Onozaki T, Ozeki Y et al. Flavonoid biosynthesis in pink-flowered cultivars derived from ’William Sim’ carnation (Dianthus caryophyllus) . Journal of the Japanese Society for Hortic Sci. 2001; 70: 315-9.

[53] Shimizu K, Ohnishi N, Morikawa N et al. A 94-bp deletion of anthocyanidin synthase gene in acyanic flower lines of lisianthus Eustoma grandiflorum (Raf.) Shinn. Journal of the Japanese society for Hortic Sci. 2011; 80: 434-42.

[54] Yang CL, Zhang J, Fang SZ et al. Cloning of ANS gene from oriental lily and phenotypic analysis of overexpression in Arabidopsis thaliana . Molecular Plant Breeding. 2021; 19: 6741-6.

[55] Wei W, Li S, Wang Y et al. Metabolome-based genome-wide association study provides genetic insights into the natural variation of foxtail millet . Front Plant Sci. 2021; 12: 665530.

[56] Chen J, Hu X, Shi T et al. Metabolite-based genome-wide association study enables dissection of the flavonoid decoration pathway of wheat kernels. Plant Biotechnol J. 2020; 18: 1722-35.

[57] Chen J, Xue MY, Liu HB et al. Exploring the genic resources underlying metabolites through mGWAS and mQTL in wheat: from large-scale gene identification and pathway elucidation to crop improvement. Plant Commun. 2021; 2: 100216.

[58] Fang CY, Luo J . Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism. Plant J. 2019; 97: 91-100.

[59] Ritchie MD, Holzinger ER, Li RW et al. Methods of integrating data to uncover genotype-phenotype interactions. Nat Rev Genet. 2015; 16: 85-97.

[60] Holzinger ER, Ritchie MD . Integrating heterogeneous high-throughput data for meta-dimensional pharmacogenomics and disease-related studies. Pharmacogenomics. 2012; 13: 213-22.

[61] Lu WJ, Du Q, Xiao L et al. Multi-omics analysis provides insights into genetic architecture of flavonoid metabolites in Populus . Ind Crop Prod. 2021; 168: 113612.

[62] Xia E-H, Tong W, Wu S et al. Tea plant genomics: achievements, challenges and perspectives. Hortic Res. 2020;7:7.

[63] Zhao J, Li PH, Xia T et al. 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. 2020; 40: 667-88.

[64] Zhu JY, Xu Q, Zhao S 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. 2020; 290: 110306.

[65] Li PH, Chen B, Zhang G et al. Regulation of anthocyanin and proanthocyanidin biosynthesis by Medicago truncatula bHLH transcription factor MtTT8. New Phytol. 2016; 210: 905-21.

[66] Xu HF, Wang N, Liu J et al. The molecular mechanism underlying anthocyanin metabolism in apple using the MdMYB16 and MdbHLH33 genes . Plant Mol Biol. 2017; 94: 149-65.

[67] Li MZ, Li Y, Guo L et al. Functional characterization of tea (Camellia sinensis) MYB4a transcription factor using an integrative approach . Front Plant Sci. 2017; 8: 943.

[68] Jin JQ, Liu YF, Ma CL et al. A novel F3’5’H allele with 14 bp deletion is associated with high catechin index trait of wild tea plants and has potential use in enhancing tea quality. J Agric Food Chem. 2018; 66: 10470-8.

[69] Kong W, Sun T, Zhang C et al. Comparative transcriptome analysis reveals the mechanisms underlying differences in salt tolerance between indica and japonica rice at seedling stage . Front Plant Sci. 2021; 12: 725436.

[70] Rymaszewski W, Vile D, Bediee A et al. Stress-related gene expression reflects morphophysiological responses to water deficit. Plant Physiol. 2017; 174: 1913-30.

[71] Pertea M, Kim D, Pertea GM et al. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nature Protoc. 2016; 11: 1650-67.

[72] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[73] Minh BQ, Schmidt HA, Chernomor O et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020; 37: 1530-4.

[74] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64.

[75] Pertea M, Petea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[76] Niknafs YS, Pandian B, Iyer HK et al. TACO produces robust multisample transcriptome assemblies from RNA-seq. Nat Methods. 2017; 14: 68-70.

[77] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and samtools. Bioinformatics. 2009; 25: 2078-9.

[78] Fu LM, Niu BF, Zhu ZW et al. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012; 28: 3150-2.

[79] Huerta-Cepas J, Szklarczyk D, Heller D et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019; 47: D309-14.

[80] Chen CJ, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[81] Browning BL, Tian X, Zhou Y et al. Fast two-stage phasing of large-scale sequence data. Am J Hum Genet. 2021; 108: 1880-90.

[82] Chi Z, Shan-Shan D, Jun-Yang X et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2018; 35: 1786-8.

[83] Yin L, Zhang H, Tang Z et al. rMVP: a memory-efficient, visualization-enhanced, and parallel-accelerated tool for genome-wide association study. Genomics Proteomics Bioinformatics. 2021; 19: 619-28.

[84] Pang JL, Fu J, Zong N et al. Kernel size-related genes revealed by an integrated eQTL analysis during early maize kernel development. Plant J. 2019; 98: 19-32.