Abstract
Masson pine (Pinus massoniana L.) is one of the most important species for oleoresin production worldwide. However, the mechanism of oleoresin synthesis in masson pine remains poorly studied. We revealed the key differentially expressed genes associated with terpenoid biosynthesis in trees with different yields of oleoresin in different tissues (trunk xylem and needle) during different seasons (spring and autumn) based on the combined application of the SMRT full-length transcriptome and Illumina RNA sequencing under field operations. A total of 76,754 full-length non-redundant transcripts were obtained in the study. Transcriptome data were annotated with KEGG and GO, and they have been shown to be of high quality by qRT-PCR. Chloroplastic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway and cytosolic mevalonate (MVA) pathway genes have different preferences under different conditions. Moreover, the phylogenetic analysis results revealed three classes of DXS genes with different functions. Among the tissue-specific genes, AACT1 was expressed in the trunk xylem at the highest levels, while MCS (1, 4) and HDR2 were significantly upregulated in the needle. Among stage-specific genes, DXS6 and HDS6 were expressed in the trunk xylem in August (autumn) at the highest levels. Furthermore, all of the terpene synthase (TPS) differentially expressed genes (DEGs) were clustered into the gymnosperm-specific TPS-d subfamily, including sesquiterpene synthase and diterpene synthase genes. According to the protein sequence alignment, diTPS1 and diTPS5 with three conserved active sites catalysed the synthesis of diterpene resin acids. In the low-yield oleoresin clone, diTPS (1, 5) and GGPPS (4, 6) were upregulated in the trunk xylem. The integration of the transcriptomic data indicated that the AACT, MCS, HDR, DXS, HDS, and TPS (with two conserved active sites) may play critical roles in the terpenoids biosynthesis of masson pine. The results of our study can be useful for developing an understanding of the coniferous genome and to significantly accelerate confer breeding.
Similar content being viewed by others
Abbreviations
- MEP pathway:
-
2-C-methyl-D-erythritol 4-phosphate pathway
- DXS:
-
1-deoxy-d-xylulose-5-phosphate synthase
- DXR:
-
1-deoxy-d-xylulose-5-phosphate reductoisomerase
- MCT:
-
2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase
- CMK:
-
4-diphosphocytidyl-2-C-methyl-D-erythritol kinase
- MCS:
-
2-C-methyl-D-erythritol 2, 4-cyclodiphosphate synthase
- HDS:
-
2-C-methyl-D-erythritol 2, 4-cyclodiphosphate synthase
- HDR:
-
2-C-methyl-D-erythritol 2, 4-cyclodiphosphate reductase
- MVA pathway:
-
cytosolic mevalonate pathway
- AACT:
-
acetyl-CoA acetyltransferase
- HMGS:
-
hydroxy methylglutaryl-CoA synthase
- HMGR:
-
hydroxy methylglutaryl-CoA reductase
- MK:
-
mevalonate kinase
- PMK:
-
phosphomevalonate kinase
- MDD:
-
diphosphomevalonate decarboxylase
- IDI:
-
isopentenyl-diphosphate delta-isomerase
- IPP:
-
isopentenyl diphosphate
- DMAPP:
-
dimethylallyl diphosphate
- GPPS:
-
geranyl diphosphate synthase
- FPPS:
-
farnesyl pyrophosphate synthase
- GGPPS:
-
geranylgeranyl diphosphate synthase
- P450:
-
cytochrome P450
- TPS:
-
terpenoid synthase
- diTPS:
-
diterpene synthase
- Nr:
-
NCBI non-redundant protein sequences
- COG:
-
Clusters of Orthologous Groups
- KOG:
-
EuKaryotic Ortholog Groups
- eggNOG:
-
evolutionary genealogy of genes: non-supervised orthologous groups of proteins
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- Swiss-Prot:
-
a manually annotated and reviewed protein sequence database
- Pfam:
-
protein family
- GO:
-
gene ontology
- SSR:
-
simple sequence repeats
- CDS:
-
coding region sequence
- qRT-PCR:
-
quantitative real-time PCR
- ROI:
-
reads of insert
- CCS:
-
circular consensus sequences
- FLNC:
-
full-length of non-chimeric
- DEGs:
-
differentially expressed genes
References
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25(17):3389–3402. https://doi.org/10.1093/nar/25.17.3389
Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106
Au KF, Underwood JG, Lee L, Wong WH (2012) Improving PacBio long read accuracy by short read alignment. PLoS One 7(10):e46679. https://doi.org/10.1371/journal.pone.0046679
Bohlmann J, Phillips M, Ramachandiran V, Katoh S, Croteau R (1999) cDNA cloning, characterization, and functional expression of four new monoterpene synthase members of the Tpsd gene family from grand fir (Abies grandis). Arch Biochem Biophys 368:232–243. https://doi.org/10.1006/abbi.1999.1332
Burke C, Croteau R (2002) Interaction with the small subunit of geranyl diphosphate synthase modifies the chain length specificity of geranylgeranyl diphosphate synthase to produce geranyl diphosphate. J Biol Chem 277:3141–3149. https://doi.org/10.1074/jbc.M105900200
Byun-McKay A, Godard KA, Toudefallah M, Martin DM, Alfaro R, King J, Bohlmann J, Plant AL (2006) Wound-induced terpene synthase gene expression in Sitka spruce that exhibit resistance or susceptibility to attack by the white pine weevil. Plant Physiol 140:1009–1021. https://doi.org/10.1104/pp.105.071803
Celedon JM, Yuen MMS, Chiang A, Henderson H, Reid KE, Bohlmann J (2017) Cell-type- and tissue-specific transcriptomes of the white spruce (Picea glauca) bark unmask fine-scale spatial patterns of constitutive and induced conifer defense. Plant J 92(4):710–726. https://doi.org/10.1111/tpj.13673
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, Turner SW, Korlach J (2013) Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 10(6):563–569. https://doi.org/10.1038/nmeth.2474
Rodrigues-Corrêa KC d S, de Lima JC, Fett-Neto AG (2012) Pine oleoresin: tapping green chemicals, biofuels, food protection, and carbon sequestration from multipurpose trees. Food and Energy Security 1:81–93. https://doi.org/10.1002/fes3.13
da Silva Rodrigues KC, Apel MA, Henriques AT, Fett-Neto AG (2011) Efficient oleoresin biomass production in pines using low cost metal containing stimulant paste. Biomass Bioenergy 35(10):4442–4448. https://doi.org/10.1016/j.biombioe.2011.08.021
de Oliveira Junkes CF, de Araújo Júnior AT, de Lima JC, et al (2019) Resin tapping transcriptome in adult slash pine (Pinus elliottii var. elliottii). Industrial Crops and Products. https://doi.org/10.1016/j.indcrop.2019.111545
Eid J, Fehr A, Gray J, Luong K, Lyle J, Otto G, Peluso P, Rank D, Baybayan P, Bettman B, Bibillo A, Bjornson K, Chaudhuri B, Christians F, Cicero R, Clark S, Dalal R, deWinter A, Dixon J, Foquet M, Gaertner A, Hardenbol P, Heiner C, Hester K, Holden D, Kearns G, Kong X, Kuse R, Lacroix Y, Lin S, Lundquist P, Ma C, Marks P, Maxham M, Murphy D, Park I, Pham T, Phillips M, Roy J, Sebra R, Shen G, Sorenson J, Tomaney A, Travers K, Trulson M, Vieceli J, Wegener J, Wu D, Yang A, Zaccarin D, Zhao P, Zhong F, Korlach J, Turner S (2009) Real-time DNA sequencing from single polymerase molecules. Science (New York, NY) 323(5910):133–138. https://doi.org/10.1126/science.1162986
Faldt J, Martin D, Miller B et al (2003) Traumatic resin defense in Norway spruce (Picea abies): methyl jasmonate-induced terpene synthase gene expression, and cDNA cloning and functional characterization of (+)-3-carene synthase. Plant Mol Biol 51(1):119–133. https://doi.org/10.1023/a:1020714403780
Franceschi VR, Krokene P, Christiansen E, Krekling T (2005) Anatomical and chemical defenses of conifer bark against bark beetles and other pests. New Phytol 167(2):353–375. https://doi.org/10.1111/j.1469-8137.2005.01436.x
Hackl T, Hedrich R, Schultz J, Förster F (2014) Proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 21:21–3011. https://doi.org/10.1093/bioinformatics/btu392
Hall DE, Zerbe P, Jancsik S, Quesada AL, Dullat H, Madilao LL, Yuen M, Bohlmann J (2013) Evolution of conifer diterpene synthases: diterpene resin acid biosynthesis in lodgepole pine and jack pine involves monofunctional and bifunctional diterpene synthases. Plant Physiol 161(2):600–616. https://doi.org/10.1104/pp.112.208546
He Y, Yan Z, Du Y et al (2017) Molecular cloning and expression analysis of two key genes, HDS and HDR, in the MEP pathway in Pyropia haitanensis. Sci Rep 7:17499. https://doi.org/10.1038/s41598-017-17521-9
Henry LK, Thomas ST, Widhalm JR, Lynch JH, Davis TC, Kessler SA, Bohlmann J, Noel JP, Dudareva N (2018) Contribution of isopentenyl phosphate to plant terpenoid metabolism. Nature Plants 4:721–729. https://doi.org/10.1038/s41477-018-0220-z
Hsieh MH, Goodman HM (2005) The Arabidopsis IspH homolog is involved in the plastid nonmevalonate pathway of isoprenoid biosynthesis. Plant Physiol 138(2):641–653. https://doi.org/10.1104/pp.104.058735
Huber DPW, Philippe RN, Godard KA, Sturrock RN, Bohlmann J (2005) Characterization of four terpene synthase cDNAs from methyl jasmonate-induced Douglas-fir, Pseudotsuga menziesii. Phytochemistry 66(12):1427–1439. https://doi.org/10.1016/j.phytochem.2005.04.030
Kawoosa T, Singh H, Kumar A, Sharma SK, Devi K, Dutt S, Vats SK, Sharma M, Ahuja PS, Kumar S (2010) Light and temperature regulated terpene biosynthesis: hepatoprotective monoterpene picroside accumulation in Picrorhiza kurrooa. Functional and Integrative Genomics 10(3):393–404. https://doi.org/10.1007/s10142-009-0152-9
Keeling CI, Bohlmann J (2006) Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens. New Phytol 170:657–675. https://doi.org/10.1111/j.1469-8137.2006.01716.x
KIM SM, Kuzuyama T, Chang Y-J, Kim S-U (2005) Functional identification of Ginkgo biloba 1-deoxy-D-xylulose 5-phosphate synthase (DXS) gene by using Escherichia coli disruptants defective in DXS gene. Agric Chem 48(2):101–104
Kim SM, Kuzuyama T, Chang YJ, Kim SU (2006) Cloning and characterization of 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase (MECS) gene from Ginkgo biloba. Plant Cell Rep 25(8):829–835. https://doi.org/10.1007/s00299-006-0136-3
Kim YB, Kim SM, Kang MK, Kuzuyama T, Lee JK, Park SC, Shin SC, Kim SU (2009) Regulation of resin acid synthesis in Pinus densiflora by differential transcription of genes encoding multiple 1-deoxy-d-xylulose 5-phosphate synthase and 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase genes. Tree Physiol 29(5):737–749. https://doi.org/10.1093/treephys/tpp002
Li B, Dewey CN (2014) RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. In: Bioinformatics: The Impact of Accurate Quantification on Proteomic and Genetic Analysis and Research pp 12(1):323–323. https://doi.org/10.1186/1471-2105-12-323
Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22(13):1658–1659. https://doi.org/10.1093/bioinformatics/btl158
Liu Q, Zhou Z, Wei Y, Shen D, Feng Z, Hong S (2015) Genome-wide identification of differentially expressed genes associated with the high yielding of oleoresin in secondary xylem of masson pine (Pinus massoniana Lamb) by transcriptomic analysis. PLoS One 10(7):e0132624. https://doi.org/10.1371/journal.pone.0132624
Mahmoud SS, Croteau RB (2002) Strategies for transgenic manipulation of monoterpene biosynthesis in plants. Trends Plant Sci 7(8):366–373. https://doi.org/10.1016/S1360-1385(02)02303-8
Mao J, He Z, Hao J et al (2019) Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses. PeerJ 7:e6124. https://doi.org/10.7717/peerj.6124
McKay SAB, Hunter WL, Godard KA et al (2003) Insect attack and wounding induce traumatic resin duct development and gene expression of (-)-pinene synthase in Sitka spruce. Plant Physiol 133(1):368–378. https://doi.org/10.1104/pp.103.022723
Phillips MA, Bohlmann J, Gershenzon J (2006) Molecular regulation of induced terpenoid biosynthesis in conifers. Phytochem Rev 5(1):179–189. https://doi.org/10.1007/s11101-006-0001-6
Phillips MA, Walter MH, Ralph SG, Dabrowska P, Luck K, Urós EM, Boland W, Strack D, Rodríguez-Concepción M, Bohlmann J, Gershenzon J (2007) Functional identification and differential expression of 1-deoxy-D-xylulose 5-phosphate synthase in induced terpenoid resin formation of Norway spruce (Picea abies). Plant Mol Biol 65(3):243–257. https://doi.org/10.1007/s11103-007-9212-5
Rather GA, Sharma A, Pandith SA, Kaul V, Nandi U, Misra P, Lattoo SK (2018) De novo transcriptome analyses reveals putative pathway genes involved in biosynthesis and regulation of camptothecin in Nothapodytes nimmoniana (Graham) Mabb. Plant Mol Biol 96(1):197–215. https://doi.org/10.1007/s11103-017-0690-9
Ro D-K, Bohlmann J (2006) Diterpene resin acid biosynthesis in loblolly pine (Pinus taeda): functional characterization of abietadiene/levopimaradiene synthase (PtTPS-LAS) cDNA and subcellular targeting of PtTPS-LAS and abietadienol/abietadienal oxidase (PtAO, CYP720B1). Phytochemistry 67(15):1572–1578. https://doi.org/10.1016/j.phytochem.2006.01.011
Sallaud C, Giacalone C, Töpfer R, Goepfert S, Bakaher N, Rösti S, Tissier A (2012) Characterization of two genes for the biosynthesis of the labdane diterpene Z-abienol in tobacco (Nicotiana tabacum) glandular trichomes. Plant J 72(1):1–17. https://doi.org/10.1111/j.1365-313X.2012.05068.x
Schmidt A, Gershenzon J (2007) Cloning and characterization of isoprenyl diphosphate synthases with farnesyl diphosphate and geranylgeranyl diphosphate synthase activity from Norway spruce (Picea abies) and their relation to induced oleoresin formation. Phytochemistry 68(21):2649–2659. https://doi.org/10.1016/j.phytochem.2007.05.037
Schmidt A, Gershenzon J (2008) Cloning and characterization of two different types of geranyl diphosphate synthases from Norway spruce (Picea abies). Phytochemistry 69:49–57. https://doi.org/10.1016/j.phytochem.2007.06.022
Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM (2015) BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 19:19–3212. https://doi.org/10.1093/bioinformatics/btv351
Tholl D, Croteau R, Gershenzon J (2001) Partial purification and characterization of the short-chain prenyltransferases, geranyl diphosphate synthase and farnesyl diphosphate synthase, from Abies grandis (grand fir). Arch Biochem Biophys 386:233–242. https://doi.org/10.1006/abbi.2000.2212
Vallinayagam R, Vedharaj S, Yang WM, Lee PS, Chua KJE, Chou SK (2014) Pine oil-biodiesel blends: a double biofuel strategy to completely eliminate the use of diesel in a diesel engine. Appl Energy 130:466–473. https://doi.org/10.1016/j.apenergy.2013.11.025
Verne S, Jaquish B, White R, Ritland C, Ritland K (2011) Global transcriptome analysis of constitutive resistance to the white pine weevil in spruce. Genome Biology and Evolution 5:851–867. https://doi.org/10.1093/gbe/evr069
Walter MH, Hans J, Strack D (2002) Two distantly related genes encoding 1-deoxy-D-xylulose 5-phosphate synthases: differential regulation in shoots and apocarotenoid-accumulating mycorrhizal roots. Plant J 31(3):243–254. https://doi.org/10.1046/j.1365-313X.2002.01352.x
Westbrook JW, Resende MFR, Munoz P et al (2013) Association genetics of oleoresin flow in loblolly pine: discovering genes and predicting phenotype for improved resistance to bark beetles and bioenergy potential. New Phytol 199:89–100. https://doi.org/10.1111/nph.12240
Wiberley AE, Donohue AR, Westphal MM, Sharkey TD (2009) Regulation of isoprene emission from poplar leaves throughout a day. Plant, Cell and Environment 32(7):939–947. https://doi.org/10.1111/j.1365-3040.2009.01980.x
Yang Y, Hou S, Fan W, Lilan L, Hui N, Xia W, Wei J (2019) Expression patterns of some genes involved in tanshinone biosynthesis in Salvia miltiorrhiza roots. Ind Crop Prod 130:606–614. https://doi.org/10.1016/j.indcrop.2019.01.001
Yu G, Wang LG, Han Y, He QY (2012) ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS A Journal of Integrative Biology 16:284–287. https://doi.org/10.1089/omi.2011.0118
Zheng Y, Jiao C, Sun H, Rosli HG, Pombo MA, Zhang P, Banf M, Dai X, Martin GB, Giovannoni JJ, Zhao PX, Rhee SY, Fei Z (2016) iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol Plant 9:1667–1670. https://doi.org/10.1016/j.molp.2016.09.014
Zhou K, Gao Y, Hoy JA, Mann FM, Honzatko RB, Peters RJ (2012) Insights into diterpene cyclization from structure of bifunctional abietadiene synthase from Abies grandis. J Biol Chem 287(9):6840–6850. https://doi.org/10.1074/jbc.M111.337592
Zulak KG, Bohlmann J (2010) Terpenoid biosynthesis and specialized vascular cells of conifer defense. J Integr Plant Biol 52(1):86–97. https://doi.org/10.1111/j.1744-7909.2010.00910.x
Funding
This research was funded by the National Key Research and Development Plan (grant number 2017YFD060030304); the Innovation Talent Program of Guizhou Province, People’s Republic of China (grant number 2016-4010); and the Postgraduate Research Project of Guizhou Province (grant number Qian JiaoHe Postgraduate Research Project Zi YJSCXJH [2018]048]).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Ethical approval
The research conducted was in compliance with the laws of the country in which the research was conducted.
Additional information
Communicated by F.P. Guerra
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Mei, L., Li, Z., Yan, Y. et al. Identification and functional study of oleoresin terpenoid biosynthesis-related genes in masson pine (Pinus massoniana L.) based on transcriptome analysis. Tree Genetics & Genomes 16, 53 (2020). https://doi.org/10.1007/s11295-020-01448-w
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11295-020-01448-w