The reference genome of Camellia chekiangoleosa provides insights into Camellia evolution and tea oil biosynthesis

Abstract

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.

References

[1] Zhu M, Shi T, Chen Y et al. Prediction of fatty acid composition in camellia oil by 1H NMR combined with PLS regression. Food Chem. 2019; 279: 339-46.

[2] Zhou J, Ai Z, Wang H et al. Phosphorus alleviates aluminum toxicity in Camellia oleifera seedlings. Int J Agric Biol. 2019; 21: 237-43.

[3] Zhang D, Yu J, Zhang R et al. Teaoil Camellia-Eastern “olive” for the world. Acta Hortic. 2006; 769: 43-8.

[4] Sokoła-Wysoczańska E, Wysoczanski T, Wagner J et al. Polyunsaturated fatty acids and their potential therapeutic role in cardiovascular system disorders-a review. Nutrients. 2018; 10: 1561.

[5] He M, Qin CX, Wang X et al. Plant unsaturated fatty acids: biosynthesis and regulation. Front Plant Sci. 2020; 11: 390.

[6] Wang Y, Sun D, Chen H et al. Fatty acid composition and antioxidant activity of tea (Camellia sinensis L.) seed oil extracted by optimized supercritical carbon dioxide. Int J Mol Sci. 2011; 12: 7708-19.

[7] Qian J, Liu Y, Ma C et al. Positive selection of squalene synthase in Cucurbitaceae plants. Int J Genomics. 2019; 2019: 1.

[8] Shang Y, Huang S . Multi-omics data-driven investigations of metabolic diversity of plant triterpenoids. Plant J. 2019; 97: 101-11.

[9] Xie Y, Wang X . Comparative transcriptomic analysis identifies genes responsible for fruit count and oil yield in the oil tea plant Camellia chekiangoleosa. Sci Rep. 2018; 8: 6637.

[10] Wang X, Zeng Q, del Mar Contreras M et al. Profiling and quantification of phenolic compounds in Camellia seed oils: natural tea polyphenols in vegetable oil. Food Res Int. 2017; 102: 184-94.

[11] Guo H, Tan H, Zhou J . Proximate composition of Camellia chekiangoleosa Hu fruit and fatty acid constituents of its seed oil. Journal of Zhejiang University (Agriculture & Life Sciences). 2010; 36: 662-9.

[12] Liu ZW, Li H, Liu JX et al. Integrative transcriptome, proteome, and microRNA analysis reveals the effects of nitrogen sufficiency and deficiency conditions on theanine metabolism in the tea plant (Camellia sinensis). Hortic Res. 2020; 7: 65.

[13] Lu L, Chen H, Wang X et al. Genome-level diversification of eight ancient tea populations in the Guizhou and Yunnan regions identifies candidate genes for core agronomic traits. Hortic Res. 2021; 8: 190.

[14] Xia EH, Tong W, Wu Q et al. Tea plant genomics: achievements, challenges and perspectives. Hortic Res. 2020; 7: 7.

[15] Yin X, Li T, Huang B et al. Complete chloroplast genome of Camellia chekiangoleosa (Theaceae), a shrub with gorgeous flowers and rich seed oil. Mitochondrial DNA Part B. 2021; 6: 840-1.

[16] Wang Y, Chen F, Ma Y et al. An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis). Hortic Res. 2021; 8: 176.

[17] Chen J, Hao Z, Guang X et al. Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation. Nature Plants. 2019; 5: 18-25.

[18] Alix K, Gérard PR, Schwarzacher T et al. Polyploidy and interspecific hybridization: partners for adaptation, speciation and evolution in plants. Ann Bot. 2017; 120: 183-94.

[19] 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.

[20] Chen JD, Chao Z, Ma JQ et al. The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant. Hortic Res. 2020; 7: 63.

[21] Yukawa Y, Takaiwa F, Shoji K et al. Structure and expression of two seed-specific cDNA clones encoding stearoyl-acyl carrier protein desaturase from sesame, Sesamum indicum L. Plant Cell Physiol. 1996; 37: 201-5.

[22] Combs R, Bilyeu K . Novel alleles of FAD2-1A induce high levels of oleic acid in soybean oil. Mol Breed. 2019; 39: 79.

[23] Zheng Y, Chen C, Liang Y et al. Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing. Tree Physiol. 2019; 39: 356-71.

[24] Shimura K, Okada A, Okada K et al. Identification of a biosynthetic gene cluster in rice for momilactones. J Biol Chem. 2007; 282: 34013-8.

[25] King AJ, Brown GD, Gilday AD et al. Production of bioactive diterpenoids in the Euphorbiaceae depends on evolutionarily conserved gene clusters. Plant Cell. 2014; 26: 3286-98.

[26] Matsuba Y, Zi J, Jones AD et al. Biosynthesis of the diterpenoid lycosantalonol via nerylneryl diphosphate in Solanum lycopersicum. PLoS One. 2015; 10: e0119302.

[27] Naoumkina MA, Modolo LV, Huhman DV et al. Genomic and coexpression analyses predict multiple genes involved in triterpene saponin biosynthesis in Medicago truncatula. Plant Cell. 2010; 22: 850-66.

[28] Krokida A, Delis C, Geisler K et al. A metabolic gene cluster in Lotus japonicus discloses novel enzyme functions and products in triterpene biosynthesis. New Phytol. 2013; 200: 675-90.

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

[30] Wei C, 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.

[31] Pu X, Dong X, Li Q et al. An update on the function and regulation of methylerythritol phosphate and mevalonate pathways and their evolutionary dynamics. J Integr Plant Biol. 2021; 63: 1211-26.

[32] Alicandri E, Paolacci AR, Osadolor S et al. On the evolution and functional diversity of terpene synthases in the Pinus species: a review. J Mol Evol. 2020; 88: 253-83.

[33] Wang JR, Lin JF, Guo LQ et al. Cloning and characterization of squalene synthase gene from Poria cocos and its up-regulation by methyl jasmonate. World J Microbiol Biotechnol. 2014; 30: 613-20.

[34] Laranjeira S, Amorim-Silva V, Esteban A et al. Arabidopsis Squalene Epoxidase 3 (SQE3) complements SQE1 and is important for embryo development and bulk squalene epoxidase activity. Mol Plant. 2015; 8: 1090-102.

[35] Xue Z, Duan L, Liu D et al. Divergent evolution of oxidosqualene cyclases in plants. New Phytol. 2012; 193: 1022-38.

[36] Chen F, Tholl D, Bohlmann J et al. 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. 2011; 66: 212-29.

[37] Zhou HC, Shamala LF, Yi XK et al. Analysis of terpene synthase family genes in Camellia sinensis with an emphasis on abiotic stress conditions. Sci Rep. 2020; 10: 933.

[38] Aubourg S, Lecharny A, Bohlmann J . Genomic analysis of the terpenoid synthase (AtTPS) gene family of Arabidopsis thaliana. Mol Gen Genomics. 2002; 267: 730-45.

[39] Busquets A, Keim V, Closa M et al. Arabidopsis thaliana contains a single gene encoding squalene synthase. Plant Mol Biol. 2008; 67: 25-36.

[40] Liu Y, Zhou J, Hu T et al. Identification and functional characterization of squalene epoxidases and oxidosqualene cyclases from Tripterygium wilfordii. Plant Cell Rep. 2020; 39: 409-18.

[41] Rao SSP, Huntley MH, Durand NC et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014; 159: 1665-80.

[42] Chen S, Zhou Y, Chen Y et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics (Oxford, England). 2018; 34: i884-90.

[43] Cheng H, Concepcion GT, Feng X et al. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021; 18: 170-5.

[44] Dudchenko O, Batra SS, Omer AD et al. De novo assembly of the genome using Hi-C yields chromosome-length scaffolds. Science. 2017; 356: 92-5.

[45] Simão FA, Waterhouse RM, Ioannidis P et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics (Oxford, England). 2015; 31: 3210-2.

[46] Ou S, Su W, Liao Y et al. Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. Genome Biol. 2019; 20: 275.

[47] Kim D, Paggi JM, Park C et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15.

[48] Cantarel BL, Korf I, Robb SMC et al. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008; 18: 188-96.

[49] Lowe TM, Eddy SR . tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997; 25: 955-64.

[50] Kalvari I, Argasinkska J, Quinones-Olvera N et al. Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families. Nucleic Acids Res. 2018; 46: D335-42.

[51] Emms DM, Kelly S . OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015; 16: 157.

[52] Katoh K, Misawa K, Kuma KI et al. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002; 30: 3059-66.

[53] Stamatakis A . RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics (Oxford, England). 2014; 30: 1312-3.

[54] Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T . trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25: 1972-3.

[55] Hedges SB, Dudley J, Kumar S . TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics (Oxford, England). 2006; 22: 2971-2.

[56] Yang Z . PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997; 13: 555-6.

[57] De Bie T, Cristianini N, Demuth JP et al. CAFE: a computational tool for the study of gene family evolution. Bioinformatics (Oxford, England). 2006; 22: 1269-71.

[58] Yu G, Smith DK, Zhu H et al. Ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017; 8: 28-36.

[59] Zhang Z, Li J, Zhao XQ et al. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinformatics. 2006; 4: 259-63.

[60] Tang H, Bowers JE, Wang X et al. Synteny and collinearity in plant genomes. Science. 2008; 320: 486-8.