Genetic basis of high aroma and stress tolerance in the oolong tea cultivar genome

Abstract

Tea plants (Camellia sinensis) are commercially cultivated in >60 countries, and their fresh leaves are processed into tea, which is the most widely consumed beverage in the world. Although several chromosome-level tea plant genomes have been published, they collapsed the two haplotypes and ignored a large number of allelic variations that may underlie important biological functions in this species. Here, we present a phased chromosome-scale assembly for an elite oolong tea cultivar, “Huangdan”, that is well known for its high levels of aroma. Based on the two sets of haplotype genome data, we identified numerous genetic variations and a substantial proportion of allelic imbalance related to important traits, including aroma- and stress-related alleles. Comparative genomics revealed extensive structural variations as well as expansion of some gene families, such as terpene synthases (TPSs), that likely contribute to the high-aroma characteristics of the backbone parent, underlying the molecular basis for the biosynthesis of aroma-related chemicals in oolong tea. Our results uncovered the genetic basis of special features of this oolong tea cultivar, providing fundamental genomic resources to study evolution and domestication for the economically important tea crop.

References

[1] Chen, F. et al. The sequenced angiosperm genomes and genome databases. Front. Plant Sci. 9, 418 (2018).

[2] Zhang, X. T., Wu, R. X., Wang, Y. B., Yu, J. X. & Tang, H. B. Unzipping haplotypes in diploid and polyploid genomes. Comput. Struct. Biotechnol. J. 18, 66-72 (2019).

[3] Zhang, X. T., Zhang, S. C., Zhao, Q., Ming, R. & Tang, H. B. Assembly of allele-aware, chromosomal-scale autopolyploid genomes based on Hi-C data. Nat. Plants 5, 833-845 (2019).

[4] Zhang, J. S. et al. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nat. Genet. 50, 1565-1573 (2018).

[5] Chen, H. T. et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat. Commun. 11, 2494 (2020).

[6] Chen, X. Q. et al. A high-quality Brassica napus genome reveals expansion of transposable elements, subgenome evolution and disease resistance. Plant Biotechnol. J. https://doi.org/10.1111/pbi.13493 (2020).

[7] Zhang, X. T. et al. Genomes of the banyan tree and pollinator wasp provide insights into fig-wasp coevolution. Cell 183, 875-889 (2020).

[8] Drew, L. The growth of tea. Nature 566, S2-S4 (2019).

[9] Chen, J. D. et al. The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant. Hortic. Res. 7, 63 (2020).

[10] Xia, E. H., Tong, W., Wu, Q., Wei, S. & Wan, X. C. Tea plant genomics: achievements, challenges and perspectives. Hortic. Res. 7, 7 (2020).

[11] Zhou, Y., Zeng, L. T., Hou, X. L., Liao, Y. Y. & Yang, Z. Y. Low temperature synergistically promotes wounding-induced indole accumulation by INDUCER OF CBF EXPRESSION-mediated alterations of jasmonic acid signaling in Camellia sinensis. J. Exp. Bot. 71, 2172-2185 (2020).

[12] Dodds, P. N. & Rathjen, J. P. Plant immunity: towards an integrated view of plant-pathogen interactions. Nat. Rev. Genet. 11, 539-548 (2010).

[13] Zeng, L. T., Watanabe, N. & Yang, Z. Y. Understanding the biosyntheses and stress response mechanisms of aroma compounds in tea (Camellia sinensis) to safely and effectively improve tea aroma. Crit. Rev. Food Sci. Nutr. 59, 2321-2334 (2018).

[14] Wei, C. L. 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).

[15] Xia, E. H. et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into genome evolution and adaptation of tea plants. Mol. Plant 13, 1013-1026 (2020).

[16] Zhang, Q. J. et al. The chromosome-level reference genome of tea tree unveils recent bursts of non-autonomous LTR retrotransposons to drive genome size evolution. Mol. Plant 13, 935-938 (2020).

[17] Wang, X., Feng, H., Chang, Y., Ma, C. & Yang, Y. Population sequencing enhances understanding of tea plant evolution. Nat. Commun. 11, 4447 (2020).

[18] 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).

[19] Zhang, W. Y. et al. Genome assembly of wild tea tree DASZ reveals pedigree and selection history of tea varieties. Nat. Commun. 11, 3719 (2020).

[20] Fuchinoue, Y. Analysis of self-incompatibility alleles of major varieties of tea. Jpn. Agr. Res. Q. 13, 43-48 (1979).

[21] Huang, X. H. et al. Genomic analysis of hybrid rice varieties reveals numerous superior alleles that contribute to heterosis. Nat. Commun. 6, 6258 (2015).

[22] Shao, L. et al. Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis. Proc. Natl Acad. Sci. USA 116, 5653-5658 (2019).

[23] Zeng, L., Zhou, X., Su, X. & Yang, Z. Chinese oolong tea: An aromatic beverage produced under multiple stresses. Trends Food Sci. Technol. 106, 242-253 (2020).

[24] Zheng, Y. C. et al. Transcriptome and metabolite profiling reveal novel insights into volatile heterosis in the tea plant (Camellia Sinensis). Molecules 24, 3380 (2019).

[25] Cheng, H., Concepcion, G. T., Feng, X., Zhang, H. & Li, H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods 18, 170-175 (2021).

[26] Guan, D. et al. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics 36, 2896-2898 (2020).

[27] Simão, F. A., Waterhouse, R. M., Panagiotis, I., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[28] Ou, S., Chen, J. & Jiang, N. Assessing genome assembly quality using the LTR Assembly Index (LAI). Nucleic Acids Res. 46, e126 (2018).

[29] Cantarel, B. L. et al. MAKER: An easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 18, 188-196 (2008).

[30] Abrusán, G., Grundmann, N., DeMester, L. & Makalowski, W. TEclass-a tool for automated classification of unknown eukaryotic transposable elements. Bioinformatics 25, 1329-1330 (2009).

[31] Ou, S. J. & Jiang, N. LTR_retriever: a highly accurate and sensitive program for identification of long terminal-repeat retrotransposons. Plant Physiol. 176, 1410-1422 (2017).

[32] Alonge, M., Wang, X., Benoit, M., Soyk, S. & Pereira, L. Major impacts of widespread structural variation on gene expression and crop improvement in tomato. Cell 182, 145-161 (2020).

[33] Chen, F., Tholl, D., Bohlmann, J. R. & Pichersky, E. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 66, 212-229 (2011).

[34] Liu, G. et al. Implementation of CsLIS/NES in linalool biosynthesis involves transcript splicing regulation in Camellia sinensis. Plant Cell Environ. 41, 176-186 (2017).

[35] Zhou, Y. et al. Formation of (E)-nerolidol in tea (Camellia sinensis) leaves exposed to multiple stresses during tea manufacturing. Food Chem. 231, 78-86 (2017).

[36] Wang, X. W. et al. Formation of α-Farnesene in tea (Camellia sinensis) leaves induced by herbivore-derived wounding and its effect on neighboring tea plants. Int. J. Mol. Sci. 20, 4151 (2019).

[37] Xu, Q. S. et al. Unraveling a crosstalk regulatory network of temporal aroma accumulation in tea plant (Camellia sinensis) leaves by integration of metabolomics and transcriptomics. Environ. Exp. Bot. 149, 81-94 (2018).

[38] Stockinger, E. J., Gilmour, S. J. & Thomashow, M. F. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc. Natl Acad. Sci. USA 94, 1035-1040 (1997).

[39] Liu, Q. et al. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell 10, 1391-1406 (1998).

[40] Gilmour, S. J. et al. Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J. 16, 433-442 (1998).

[41] Wang, P. J. et al. Identification of CBF transcription factors in tea plants and a survey of potential CBF target genes under low temperature. Int. J. Mol. Sci. 20, 5137 (2019).

[42] Hu, Z., Ban, Q., Hao, J., Zhu, X. & Li, Y. Genome-wide characterization of the c-repeat binding factor (CBF) gene family involved in the response to abiotic stresses in tea plant (Camellia sinensis). Front. Plant Sci. 11, 921 (2020).

[43] O’Malley, R. C. et al. Cistrome and epicistrome features shape the regulatory DNA landscape. Cell 166, 1598 (2016).

[44] Schuh, C. & Schieberle, P. Characterization of the key aroma compounds in the beverage prepared from Darjeeling black tea: quantitative differences between tea leaves and infusion. J. Agric. Food Chem. 54, 916-924 (2006).

[45] Chen, H. et al. Combinatorial evolution of a terpene synthase gene cluster explains terpene variations in Oryza. Plant Physiol. 182, 480-492 (2020).

[46] Zhou, H. C., Shamala, L. F., Yi, X. K., Yan, Z. & Wei, S. Analysis of terpene synthase family genes in Camellia sinensis with an emphasis on abiotic stress conditions. Sci. Rep. 10, 933 (2020).

[47] Jin, J. et al. Scenarios of genes-to-terpenoids network led to the identification of a novel α/β-Farnesene/β-ocimene synthase in Camellia sinensis. Int. J. Mol. Sci. 21, 655 (2020).

[48] Yang, Z. Y., Baldermann, S. & Watanabe, N. Recent studies of the volatile compounds in tea. Food Res. Int. 53, 585-599 (2013).

[49] Lin, J., Shi, Y. X., Dong, C. W. & Wang, X. C. Headspace volatiles influenced by infusion matrix and their release persistence: a case study of oolong tea. Food Sci. Biotechnol. 28, 1349-1358 (2019).

[50] Dolezel, J., Greilhuber, J. & Suda, J. Estimation of nuclear DNA content in plants using flow cytometry. Nat. Protoc. 2, 2233-2244 (2007).

[51] Wenger, A. M., Peluso, P., Rowell, W. J., Chang, P. C. & Hunkapiller, M. W. Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome. Nat. Biotechnol. 37, 1155-1162 (2019).

[52] Servant, N. et al. HiC-Pro: An optimized and flexible pipeline for Hi-C data processing. Genome Biol. 16, 259 (2015).

[53] Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 3, 95-98 (2016).

[54] Dudchenko, O. et al. De novo assembly of the aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92-95 (2017).

[55] Ranallo-Benavidez, T. R., Jaron, K. S. & Schatz, M. C. GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat. Commun. 11, 1432 (2020).

[56] Chen, C. J. et al. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194-1202 (2020).

[57] Bao, Z. R. & Eddy, S. R. Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 12, 1269-1276 (2002).

[58] Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, 351-358 (2005).

[59] Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573-580 (1999).

[60] Haas, B. J. et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494-1512 (2013).

[61] Haas, B. J. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 31, 5654-5666 (2003).

[62] Lomsadze, A., Ter-Hovhannisyan, V., Chernoff, Y. O. & Borodovsky, M. Gene identification in novel eukaryotic genomes by self-training algorithm. Nucleic Acids Res. 33, 6494-6506 (2005).

[63] Korf, I. Gene finding in novel genomes. BMC Bioinforma. 5, 59 (2004).

[64] Stanke, M., Schöffmann, O., Morgenstern, B. & Waack, S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinforma. 7, 62 (2006).

[65] Pertea, M., Kim, D., Pertea, G. M., Leek, J. T. & Salzberg, S. L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650-1667 (2016).

[66] Huerta Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309-D314 (2018).

[67] Wu, T. D. & Watanabe, C. K. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

[68] Tang, H. B. et al. Synteny and collinearity in plant genomes. Science 320, 486-488 (2008).

[69] Kurtz, S. et al. Versatile and open software for comparing large genomes. Genome Biol. 5, R12 (2004).

[70] Maria, N. & Schatz, M. C. Assemblytics: a web analytics tool for the detection of variants from an assembly. Bioinformatics 32, 3021-3023 (2016).