Chromosome-level genome assembly and resequencing of camphor tree (Cinnamomum camphora) provides insight into phylogeny and diversification of terpenoid and triglyceride biosynthesis of Cinnamomum

Abstract

Cinnamomum species attract attentions owing to their scents, medicinal properties, and ambiguous relationship in the phylogenetic tree. Here, we report a high-quality genome assembly of Cinnamomum camphora, based on which two whole-genome duplication (WGD) events were detected in the C. camphora genome: one was shared with Magnoliales, and the other was unique to Lauraceae. Phylogenetic analyses illustrated that Lauraceae species formed a compact sister clade to the eudicots. We then performed whole-genome resequencing on 24 Cinnamomum species native to China, and the results showed that the topology of Cinnamomum species was not entirely consistent with morphological classification. The rise and molecular basis of chemodiversity in Cinnamomum were also fascinating issues. In this study, six chemotypes were classified and six main terpenoids were identified as major contributors of chemodiversity in C. camphora by the principal component analysis. Through in vitro assays and subcellular localization analyses, we identified two key terpene synthase (TPS) genes (CcTPS16 and CcTPS54), the products of which were characterized to catalyze the biosynthesis of two uppermost volatiles (i.e. 1,8-cineole and (iso)nerolidol), respectively, and meditate the generation of two chemotypes by transcriptional regulation and compartmentalization. Additionally, the pathway of medium-chain triglyceride (MCT) biosynthesis in Lauraceae was investigated for the first time. Synteny analysis suggested that the divergent synthesis of MCT and long-chain triglyceride (LCT) in Lauraceae kernels was probably controlled by specific medium-chain fatty acyl-ACP thioesterase (FatB), type-B lysophosphatidic acid acyltransferase (type-B LPAAT), and diacylglycerol acyltransferase 2b (DGAT 2b) isoforms during co-evolution with retentions or deletions in the genome.

References

[1] Cronquist A, ed. An Integrated System of Classification of Flowering Plants . Cambridge, England: Columbia University Press, 1981.

[2] Chaw SM, Liu YC, Wu YW et al. Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution. Nat Plants. 2019; 5: 63-73.

[3] Ma J, Sun P, Wang D et al. The Chloranthus sessilifolius genome provides insight into early diversification of angiosperms. Nat Commun. 2021; 12: 6929.

[4] Qin LY, Hu Y, Wang J et al. Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome. Nat Plants. 2021; 7: 1239-53.

[5] Rendón-Anaya M, Ibarra-Laclette E, Méndez-Bravo A et al. The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation. Proc Natl Acad Sci U S A. 2019; 116: 17081-9.

[6] Chen SP, Sun WH, Xiong YF et al. The Phoebe genome sheds light on the evolution of magnoliids. Hortic Res. 2020; 7: 146-59.

[7] Hu L, Xu Z, Wang M et al. The chromosome-scale reference genome of black pepper provides insight into piperine biosynthesis. Nat Commun. 2019; 10: 1-11.

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

[9] Chen YC, Li Z, Zhao YX et al. The Litsea genome and the evolution of the laurel family. Nat Commun. 2020; 11: 1675-89.

[10] Lv Q, Qiu J, Liu J et al. The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis. Plant J. 2020; 103: 1910-23.

[11] Chen F, Dong W, Zhang J et al. The sequenced angiosperm genomes and genome databases. Front Plant Sci. 2018; 9: 418.

[12] Rohwer JG . Toward a phylogenetic classification of the Lauraceae: evidence from matK sequences. Syst Bot. 2000; 25: 60-71.

[13] Li HW et al. Flora of China . Beijing, China and St. Louis, Missouri, USA: Science Press and Missouri Botanical Garden Press; 2008.

[14] Sun BX, Zhao HL . A new species of Cinnamomum from Yunnan. J Yunnan Univ. 1991; 13: 93-4.

[15] Huang JF, Li L, van der Werff H et al. Origins and evolution of cinnamon and camphor: a phylogenetic and historical biogeographical analysis of the Cinnamomum group (Lauraceae) . Mol Phylogenet Evol. 2016; 96: 33-44.

[16] Song K, He M, Yu J et al. Characterization of the chloroplast genome of the family Lauraceae plant species, Cinnamomum cassia . Mitochondrial DNA B Resour. 2019; 4: 3906-7.

[17] Damasceno CSB, Fabri Higaki NT, Dias JFG et al. Chemical composition and biological activities of essential oils in the family Lauraceae: a systematic review of the literature. Planta Med. 2019; 85: 1054-72.

[18] Tholl D . Biosynthesis and biological functions of terpenoids in plants. Adv Biochem Eng Biotechnol. 2015; 148: 63-106.

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

[20] Liu G, Liu JJ, He ZR et al. Implementation of CsLIS/NES in linalool biosynthesis involves transcript splicing regulation in Camellia sinensis . Plant Cell Environ. 2017; 41: 176-86.

[21] Nagegowda DA . Plant volatile terpenoid metabolism: biosynthetic genes, transcriptional regulation and subcellular compartmentation. FEBS Lett. 2010; 584: 2965-73.

[22] Fan R, Li Y, Li C et al. Differential microRNA analysis of glandular trichomes and young leaves in Xanthium strumarium L. reveals their putative roles in regulating terpenoid biosynthesis. PLoS One. 2015; 10: e0139002.

[23] Xu X, Jiang Q, Ma X et al. Deep sequencing identifies tissue-specific microRNAs and their target genes involving in the biosynthesis of tanshinones in Salvia miltiorrhiza . PLoS One. 2014; 9: e111679.

[24] Li X, Xu Y, Shen S et al. Transcription factor CitERF71 activates the terpene synthase gene CitTPS16 involved in the synthesis of E-geraniol in sweet orange fruit . J Exp Bot. 2017; 68: 4929-38.

[25] Nieuwenhuizen NJ, Chen X, Wang MY et al. Natural variation in monoterpene synthesis in kiwifruit: transcriptional regulation of terpene synthases by NAC and ethylene-insensitive 3-like transcription factors. Plant Physiol. 2015; 167: 1243-58.

[26] Yang T, Li J, Wang H-X et al. A geraniol-synthase gene from. Phytochemistry. 2005; 66: 285-93.

[27] Yahyaa M, Matsuba Y, Brandt W et al. Identification, functional characterization, and evolution of terpene synthases from a basal dicot. Plant Physiol. 2015; 169: 1683-97.

[28] Lin YL, Lee YR, Huang WK et al. Characterization of S-(+)-linalool synthase from several provenances of Cinnamomum osmophloeum . Tree Genet Genomes. 2014; 10: 75-86.

[29] Gan Y, Song Y, Chen Y et al. Transcriptome analysis reveals a composite molecular map linked to unique seed oil profile of Neocinnamomum caudatum (Nees) Merr. BMC Plant Biol. 2018; 18: 303-18.

[30] Ibarra-Laclette E, Méndez-Bravo A, Pérez-Torres CA et al. Deep sequencing of the Mexican avocado transcriptome, an ancient angiosperm with a high content of fatty acids. BMC Genomics. 2015; 16: 599-617.

[31] Kilaru A, Cao X, Dabbs PB et al. Oil biosynthesis in a basal angiosperm: transcriptome analysis of Persea Americana mesocarp. BMC Plant Biol. 2015; 15: 203-22.

[32] Pollard MR, Anderson L, Fan C et al. A specific acyl-ACP thioesterase implicated in medium-chain fatty acid production in immature cotyledons of Umbellularia californica . Arch Biochem Biophys. 1991; 284: 306-12.

[33] Davies HM . Medium chain acyl-acp hydrolysis activities of developing oilseeds. Phytochemistry. 1993; 33: 1353-6.

[34] Nam JW, Yeon J, Jeong J et al. Overexpression of acyl-acp thioesterases, CpFatB4 and CpFatB5, induce distinct gene expression reprogramming in developing seeds of Brassica napus . Int J Mol Sci. 2019; 20: 3334.

[35] Yuan YJ, Gao L, Sun R et al. Seed-specific expression of an acyl-acyl carrier protein thioesterase CnFatB3 from coconut (Cocos nucifera L.) increases the accumulation of medium-chain fatty acids in transgenic arabidopsis seeds . Sci Hortic. 2017; 223: 5-9.

[36] Eccleston VS, Cranmer AM, Voelker TA et al. Medium-chain fatty acid biosynthesis and utilization in Brassica napus plants expressing lauroyl-acyl carrier protein thioesterase. Planta. 1996; 198: 46-53.

[37] Voelker T, Hayes TR, Cranmer AM et al. Genetic engineering of a quantitative trait: metabolic and genetic parameters influencing the accumulation of laurate in rapeseed. Plant J. 1996; 9: 229-41.

[38] Slabaugh MB, Leonard JM, Knapp SJ . Condensing enzymes from Cuphea wrightii associated with medium chain fatty acid biosynthesis. Plant J. 1998; 13: 611-20.

[39] Schutt BS, Brummel M, Schuch R et al. The role of acyl carrier protein isoforms from Cuphea lanceolata seeds in the de-novo biosynthesis of medium-chain fatty acids . Planta. 1998; 205: 263-8.

[40] Knutzon DS, Lardizabal KD, Nelsen JS et al. Cloning of a coconut endosperm cDNA encoding a 1-acyl- sn-glycerol-3-phosphate acyltransferase that accepts medium chain-length substrates . Plant Physiol. 1995; 109: 999-1006.

[41] Kim HJ, Silva JE, Iskandarov U et al. Structurally divergent lysophosphatidic acid acyltransferases with high selectivity for saturated medium chain fatty acids from Cuphea seeds. Plant J. 2016; 84: 1021-33.

[42] Iskandarow U, Silva JE, Kim HJ et al. A specialized diacylglycerol acyltransferase contributes to the extreme medium-chain fatty acid content of Cuphea seed oil. Plant Physiol. 2017; 174: 97-109.

[43] Ayme L, Jolivet P, Nicaud J-M et al. Molecular characterization of the Elaeis guineensis medium-chain fatty acid diacylglycerol acyltransferase DGAT1-1 by heterologous expression in Yarrowia lipolytica . PLoS One. 2015; 10: e0143113.

[44] Wan-yang S, Wei H, Guang-yu W et al. Study on chemical constituents of the essential oil and classification of types from Cinnamomum camphora . J Integr Plant Biol. 1989; 31: 209-14.

[45] An Q, Ren JN, Li X et al. Recent updates on bioactive properties of linalool. Food Funct. 2021; 12: 10370-89.

[46] Kulkarni M, Sawant N, Kolapkar A et al. Borneol: a promising monoterpenoid in enhancing drug delivery across various physiological barriers. AAPS PharmSciTech. 2021; 22: 145.

[47] Ma R, Su P, Jin B et al. Molecular cloning and functional identification of a high-efficiency (+)-borneol dehydrogenase from Cinnamomum camphora (L.) Presl. Plant Physiol Bioch. 2021; 158: 363-71.

[48] Hu JN, Zhang B, Zhu XM et al. Characterization of medium-chain triacylglycerol (MCT)-enriched seed oil from Cinnamomum camphora (Lauraceae) and its oxidative stability. J Agric Food Chem. 2011; 59: 4771-8.

[49] Vranová E, Coman D, Gruissem W . Network analysis of the MVA and MEP pathways for isoprenoid synthesis. Annu Rev Plant Bio. 2013; 64: 665-700.

[50] Yuan L, Voelker TA, Hawkins DJ . Modification of the substrate specificity of an acyl-acyl carrier protein thioesterase by protein engineering. Proc Natl Acad Sci U S A. 1995; 92: 10639-43.

[51] Fan LP, Liu J, Nie K et al. Synthesis of medium chain length fatty acid ethyl esters in engineered Escherichia coli using endogenously produced medium chain fatty acids. Enzyme Microb Tech. 2013; 53: 128-33.

[52] Wu ZY, ed. Flora of China . Beijing, China: Beijing Science Press; 2001.

[53] Green SA, Chen X, Nieuwenhuizen NJ et al. Identification, functional characterization, and regulation of the enzyme responsible for floral (E)-nerolidol biosynthesis in kiwifruit (Actinidia chinensis) . J Exp Bot. 2012; 63: 1951-67.

[54] Qiu FY, Wang X, Zheng Y et al. Full-length transcriptome sequencing and different Chemotype expression profile analysis of genes related to monoterpenoid biosynthesis in Cinnamomum porrectum . Int J Mol Med. 2019; 20: 6230.

[55] Marcais G, Kingsford C . A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011; 27: 764-70.

[56] Vurture GW, Sedlazeck FJ, Nattestad M et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics. 2017; 33: 2202-4.

[57] Koren S, Walenz BP, Berlin K et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017; 27: 722-36.

[58] Chin CS, Peluso P, Sedlazeck FJ et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016; 13: 1050-4.

[59] Boetzer M, Henkel CV, Jansen HJ et al. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011; 27: 578-9.

[60] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

[61] Mount DW. Using the basic local alignment search tool (BLAST). CSH Protoc. 2007; 2007: pdb.top17.

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

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

[64] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[65] Servant N, Varoquaux N, Lajoie BR et al. HiC-pro: an optimized and flexible pipeline for hi-C data processing. Genome Biol. 2015; 16: 259-70.

[66] Burton JN, Adey A, Patwardhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[67] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; 4.10.

[68] Jurka J, Kapitonov VV, Pavlicek A et al. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005; 110: 462-7.

[69] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[70] Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999; 27: 573-80.

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

[72] Slater GS, Birney E . Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005; 6: 31.

[73] Wu TD, Watanabe CK . GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005; 21: 1859-75.

[74] Stanke M, Keller O, Gunduz I et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006; 34: W435-9.

[75] Korf I. Gene finding in novel genomes. Bmc Bioinformatics. 2004; 5: 59.

[76] Ashburner M, Ball CA, Blake JA et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000; 25: 25-9.

[77] Li L Jr, Stoeckert CJ, Roos DS . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003; 13: 2178-89.

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

[79] Suyama M, Torrents D, Bork P . PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006; 34: W609-12.

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

[81] Nguyen LT, Schmidt HA, von Haeseler A et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015; 32: 268-74.

[82] Kalyaanamoorthy S, Minh BQ, Wong TKF et al. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017; 14: 587-9.

[83] Mirarab S, Reaz R, Bayzid MS et al. ASTRAL: genome-scale coalescent-based species tree estimation. Bioinformatics. 2014; 30: i541-8.

[84] Emms DM, Kelly S . STAG: species tree inference from all genes. BioRxiv. 2018. (preprint: not peer reviewed).

[85] Yang Z . PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24: 1586-91.

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

[87] Wang Y, Tang H, DeBarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49.

[88] Huang S, Li R, Zhang Z et al. The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009; 41: 1275-81.

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

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

[91] Cingolani P, Platts A, Wang LL et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff . Fly. 2012; 6: 80-92.

[92] Ronquist F, Teslenko M, van der Mark P et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012; 61: 539-42.

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

[94] Paradis E, Claude J, Strimmer K . APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004; 20: 289-90.

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

[96] Francis RM . Pophelper: an R package and web app to analyse and visualize population structure. Mol Ecol Resour. 2017; 17: 27-32.

[97] Price AL, Patterson NJ, Plenge RM et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38: 904-9.

[98] Langmead B, Salzberg SL . Fast gapped-read alignment with bowtie 2. Nat Methods. 2012; 9: 357-9.

[99] Dai X, Zhuang Z, Zhao PX . psRNATarget: a plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018; 46: W49-54.

[100] Eddy SR . Accelerated profile HMM searches. PLoS Comput Biol. 2011; 7: e1002195.

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

[102] Livak KJ, Schmittgen TD . Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001; 25: 402-8.

[103] Li-Beisson Y, Shorrosh B, Beisson F et al. Acyl-lipid metabolism. Arabidopsis Book. 2010; 8: e0133.

[104] Daniel V, Malapi WM, Anne CJ . SimpleSynteny: a web-based tool for visualization of microsynteny across multiple species. Nucleic Acids Res. 2016; 44: W41-5.

[105] Lai JS . Canoco 5: a new version of an ecological multivariate data ordination program. Biodivers Sci. 2013; 21: 765-8.