The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant

Abstract

Tea is one of the most popular nonalcoholic beverages due to its characteristic secondary metabolites with numerous health benefits. Although two draft genomes of tea plant (Camellia sinensis) have been published recently, the lack of chromosome-scale assembly hampers the understanding of the fundamental genomic architecture of tea plant and potential improvement. Here, we performed a genome-wide chromosome conformation capture technique (Hi-C) to obtain a chromosome-scale assembly based on the draft genome of C. sinensis var. sinensis and successfully ordered 2984.7 Mb (94.7%) scaffolds into 15 chromosomes. The scaffold N50 of the improved genome was 218.1 Mb, ~157-fold higher than that of the draft genome. Collinearity comparison of genome sequences and two genetic maps validated the high contiguity and accuracy of the chromosome-scale assembly. We clarified that only one Camellia recent tetraploidization event (CRT, 58.9–61.7 million years ago (Mya)) occurred after the core-eudicot common hexaploidization event (146.6–152.7 Mya). Meanwhile, 9243 genes (28.6%) occurred in tandem duplication, and most of these expanded after the CRT event. These gene duplicates increased functionally divergent genes that play important roles in tea-specific biosynthesis or stress response. Sixty-four catechin- and caffeine-related quantitative trait loci (QTLs) were anchored to chromosome assembly. Of these, two catechin-related QTL hotspots were derived from the CRT event, which illustrated that polyploidy has played a dramatic role in the diversification of tea germplasms. The availability of a chromosome-scale genome of tea plant holds great promise for the understanding of genome evolution and the discovery of novel genes contributing to agronomically beneficial traits in future breeding programs.

References

[1] Kanwar, J. et al. Recent advances on tea polyphenols. Front. Biosci. 4, 111-131 (2012).

[2] Pervin, M. et al. Beneficial effects of green tea catechins on neurodegenerative diseases. Molecules 23, 1297 (2018).

[3] Mancini, E. et al. Green tea effects on cognition, mood and human brain function: a systematic review. Phytomedicine 34, 26-37 (2017).

[4] Chen, L., Apostolides, Z. & Chen, Z. M. Global Tea Breeding: Achievements, Challenges and Perspectives (University Press-Springer, Hangzhou, Zhejiang; 2012).

[5] Chen, L., Yao, M. Z., Wang, X. C. & Yang, Y. J. Tea genetic resources in China. Int. J. Tea Sci. 8(2), 55-64 (2012).

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

[7] International Tea Committee . Annual Bulletin of Statistics (London, 2019).

[8] Rani, A., Singh, K., Sood, P., Kumar, S. & Ahuja, P. S. p-Coumarate:CoA ligase as a key gene in the yield of catechins in tea [Camellia sinensis (L.) O. Kuntze]. Funct. Integr. Genomics 9, 271-275 (2009).

[9] Singh, K. et al. An early gene of the flavonoid pathway, flavanone 3-hydroxylase, exhibits a positive relationship with the concentration of catechins in tea (Camellia sinensis). Tree Physiol. 28, 1349-1356 (2008).

[10] Rani, A., Singh, K., Ahuja, P. S. & Kumar, S. Molecular regulation of catechins biosynthesis in tea [Camellia sinensis (L.) O. Kuntze]. Gene 495, 205-210 (2012).

[11] Pang, Y. et al. Functional characterization of proanthocyanidin pathway enzymes from tea and their application for metabolic engineering. Plant Physiol. 161, 1103-1116 (2013).

[12] Deng, W. W., Ogita, S. & Ashihara, H. Ethylamine content and theanine biosynthesis in different organs of Camellia sinensis seedlings. Z. Naturforsch. C. J. Biosci. 64, 387-390 (2009).

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

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

[15] Xia, E. et al. The tea plant reference genome and improved gene annotation using long-read and paired-end sequencing data. Sci. Data 6, 122 (2019).

[16] Philippe, R. et al. A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat. Genome Biol. 14, R64 (2013).

[17] International Rice Genome Sequencing Project . The map-based sequence of the rice genome. Nature 436, 793-800 (2005).

[18] Schnable, P. S. et al. The B73 maize genome: complexity, diversity, and dynamics. Science 326, 1112-1115 (2009).

[19] Prochnik, S. et al. The cassava genome: current progress, future directions. Tropical Plant Biol. 5, 88-94 (2012).

[20] Jia, J. et al. Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation. Nature 496, 91-95 (2013).

[21] Wang, S. et al. Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes. Genome Biol. 16, 108 (2015).

[22] Ma, J. Q. et al. Large-scale SNP discovery and genotyping for constructing a high-density genetic map of tea plant using Specific-Locus Amplified Fragment Sequencing (SLAF-seq). PLoS ONE 10, e0128798 (2015).

[23] Kaplan, N. & Dekker, J. High-throughput genome scaffolding from in vivo DNA interaction frequency. Nat. Biotechnol. 31, 1143-1147 (2013).

[24] Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119-1125 (2013).

[25] Dekker, J., Marti-Renom, M. A. & Mirny, L. A. Exploring the three-dimensional organization of genomes: interpreting chromatin interaction data. Nat. Rev. Genet. 14, 390-403 (2013).

[26] Xie, T. et al. De novo plant genome assembly based on chromatin interactions: a case study of Arabidopsis thaliana. Mol. Plant 8, 489-492 (2015).

[27] Jibran, R. et al. Chromosome-scale scaffolding of the black raspberry (Rubus occidentalis L.) genome based on chromatin interaction data. Hortic. Res. 5, 8 (2018).

[28] Qiao, X. et al. Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants. Genome Biol. 20, 38 (2019).

[29] Jiao, Y. et al. Ancestral polyploidy in seed plants and angiosperms. Nature 473, 97-100 (2011).

[30] Soltis, D. E. et al. Polyploidy and angiosperm diversification. Am. J. Bot. 96, 336-348 (2009).

[31] Wendel, J. F., Jackson, S. A., Meyers, B. C. & Wing, R. A. Evolution of plant genome architecture. Genome Biol. 17, 37 (2016).

[32] Moore, M. J., Bell, C. D., Soltis, P. S. & Soltis, D. E. Using plastid genome-scale data to resolve enigmatic relationships among basal angiosperms. Proc. Natl Acad. Sci. USA 104, 19363-19368 (2007).

[33] Paterson, A. H. et al. Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 492, 423-427 (2012).

[34] Wang, J. et al. An overlooked paleotetraploidization in Cucurbitaceae. Mol. Biol. Evol. 35, 16-26 (2018).

[35] Wang, J. P. et al. Two likely auto-tetraploidization events shaped kiwi fruit genome and contributed to establishment of the Actinidiaceae family. iScience 7, 230-240 (2018).

[36] Wu, H. et al. A high-quality Actinidia chinensis (kiwifruit) genome. Horticulture Res. 6, 117 (2019).

[37] Freeling, M. Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition. Annu. Rev. Plant Biol. 60, 433-453 (2009).

[38] Li, X. et al. Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 141, 1167-1184 (2006).

[39] Fan, C., Chen, Y. & Long, M. Recurrent tandem gene duplication gave rise to functionally divergent genes in Drosophila. Mol. Biol. Evol. 25, 1451-1458 (2008).

[40] Hanada, K., Zou, C., Lehti-Shiu, M. D., Shinozaki, K. & Shiu, S. H. Importance of lineage-specific expansion of plant tandem duplicates in the adaptive response to environmental stimuli. Plant Physiol. 148, 993-1003 (2008).

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

[42] Stein, J. C. et al. Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza. Nat. Genet. 50, 285-296 (2018).

[43] Fu, Y. et al. NBS-encoding genes in Brassica napus evolved rapidly after allopolyploidization and co-localize with known disease resistance loci. Front. Plant Sci. 10, 26 (2019).

[44] Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 126, 485-493 (2001).

[45] Guo, F., Guo, Y., Wang, P., Wang, Y. & Ni, D. Transcriptional profiling of catechins biosynthesis genes during tea plant leaf development. Planta 246, 1139-1152 (2017).

[46] Jin, J. Q., Ma, J. Q., Yao, M. Z., Ma, C. L. & Chen, L. Functional natural allelic variants of flavonoid 3’,5’-hydroxylase gene governing catechin traits in tea plant and its relatives. Planta 245, 523-538 (2017).

[47] Guo, L. et al. Functional analysis of flavonoid 3’-hydroxylase and flavonoid 3’,5’-hydroxylases from tea plant (Camellia sinensis), involved in the B-ring hydroxylation of flavonoids. Gene 717, 144046 (2019).

[48] Balentine, D. A., Wiseman, S. A. & Bouwens, L. C. The chemistry of tea flavonoids. Crit. Rev. Food Sci. Nutr. 37, 693-704 (1997).

[49] Shahidi, F. Antioxidants in food and food antioxidants. Nahrung 44, 158-163 (2000).

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

[51] Bickhart, D. M. et al. Single-molecule sequencing and chromatin conformation capture enable de novo reference assembly of the domestic goat genome. Nat. Genet. 49, 643-650 (2017).

[52] Wallace, J. G. et al. Association mapping across numerous traits reveals patterns of functional variation in maize. PLoS Genet. 10, e1004845 (2014).

[53] Rao, S. S. et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 159, 1665-1680 (2014).

[54] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

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

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

[57] Durand, N. C. et al. Juicebox provides a visualization system for Hi-C contact maps with unlimited zoom. Cell Syst. 3, 99-101 (2016).

[58] Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).

[59] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

[60] Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49 (2012).

[61] Myburg, A. A. et al. The genome of Eucalyptus grandis. Nature 510, 356-362 (2014).

[62] Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284-287 (2012).

[63] Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114-2120 (2014).

[64] Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[65] Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139-140 (2010).

[66] Marchler-Bauer, A. et al. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 39, D225-D229 (2011).

[67] McDonnell, A. V., Jiang, T., Keating, A. E. & Berger, B. Paircoil2: improved prediction of coiled coils from sequence. Bioinformatics 22, 356-358 (2006).