Screening and Expression Analysis of Different Genes for Oleoresin Production in the Specific Period of Pinus elliottii×P.caribaea by RNA-Seq Technology

doi:10.7525/j.issn.1673-5102.2021.03.012

Abstract

Abstract:

The experiment was conducted to explore the differential expression of genes which regulate oleoresin synthesis and transport in Pinus elliottii×P. caribaea. At the peak of oleoresin production， P. elliottii×P. caribaea with similar growth but significant differences in oleoresin production of plants were selected， and secondary xylem tissues were collected for high-throughput sequencing platform. The results showed that 594 108 382 clean reads were generated by sequencing. 190 514 unigenes were obtained by de novo splicing， and 304 unigenes were related to the synthesis of terpenoids. The 414 unigenes were obtained， differential expression analysis， including 183 unigenes up-regulated and 231 unigenes down regulated. The genes encoding CYPs， ABC transport proteins， PLR， Strictosidine synthase activity increased in high oleoresin-yielding， while the genes encoding prenyltransferase activity increased in low oleoresin-yielding. This study enriched the transcript information of P. elliottii× P. caribaea， and provided reference value for further study of the molecular mechanism of the difference in P. elliottii×P. caribaea and Pines plant resin production.

Key words: Pinus elliottii×P. caribaea, transcriptome sequencing, oleoresin production genes, qRT-PCR

CLC Number:

S791.244

Chen-Chen ZHOU, Yi-Liang LI, Zhe WANG, Fen-Cheng ZHAO, Sui-Ying ZHONG, Chang-Ming LIN, Zhi-Qiang TAN, Fu-Ming LI, Hui-Shan WU, Wen-Bing GUO, Fang-Yan LIAO, Jian-Ge WANG. Screening and Expression Analysis of Different Genes for Oleoresin Production in the Specific Period of Pinus elliottii×P.caribaea by RNA-Seq Technology[J]. Bulletin of Botanical Research, 2021, 41(3): 419-428.

Figures/Tables 12

Table 1

Table 2

Table 3

Table 4

Table 5

Table 6

Fig.1

Fig.2

Table 7

Table 8

Table 9

Fig.3

References 28

1	吴东山，杨章旗，黄永利.不同杂种松产脂力及松脂成分的变异研究［J］.生物质化学工程，2018，52（5）：6-12.
	Wu D S，Yang Z Q，Huang Y L.Comparative analysis of resin-producing capacity and resin composition of different hybrid pine［J］.Biomass Chemical Engineering，2018，52（5）：6-12.
2	沈熙环.澳大利亚的林木遗传改良工作［J］.世界林业研究，1996，9（6）：52-56.
	Shen X H.Forest tree improvement in Australia［J］.World Forestry Research，1996，9（6）：52-56.
3	Keeling C I，Bohlmann J.Genes enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens［J］.New Phytologist，2006，170（4）：657-675.
4	Boone C K，Aukema B H，Bohlmann J，et al.Efficacy of tree defense physiology varies with bark beetle population density：a basis for positive feedback in eruptive species［J］.Canadian Journal of Forest Research，2011，41（6）：1174-1188.
5	Coppen J J W，Hone G A.Gum naval stores，turpentine and rosin from pine resin［R］.Rome：Natural Resources Institute，1995.
6	Kelkar V M，Geils B W，Becker D R，et al.How to recover more value from small pine trees：essential oils and resins［J］.Biomass and Bioenergy，2006，30（4）：316-320.
7	刘彬，刘青华，周志春，等.基于高通量转录组测序筛选马尾松抗松材线虫病相关基因［J］.林业科学研究，2019，32（5）：1-10.
	Liu B，Liu Q H，Zhou Z C，et al.Identification of candidate constitutive expressed resistant genes of pine wilt disease in Pinus massoniana based on high-throughput transcriptome sequencing［J］.Forest Research，2019，32（5）：1-10.
8	Fox H，Doron-Faigenboim A，Kelly G，et al.Transcriptome analysis of Pinus halepensis under drought stress and during recovery［J］.Tree Physiology，2018，38（3）：423-441.
9	吴东山，杨章旗，黄永利.马尾松不同半同胞家系产脂力、松脂组分分析与评价［J］.北京林业大学学报，2019，41（2）：53-61.
	Wu D S，Yang Z Q，Huang Y L.Analysis and evaluation of resin productivity and resin component among different half sibling families of Pinus massoniana［J］.Journal of Beijing Forestry University，2019，41（2）：53-61.
10	李彦杰.湿地松家系产脂性状遗传分析与评价技术研究［D］.北京：中国林业科学研究院，2013.
	Li Y J.Genetic analysis and evaluation of slash pine families on resin production related traits［D］.Beijing：Chinese Academy of Forestry，2013.
11	刘青华，周志春，范辉华，等.马尾松产脂力与生长性状的家系变异及优良家系早期选择［J］.林业科学研究，2013，26（6）：686-691.
	Liu Q H，Zhou Z C，Fan H H，et al.Variation for oleoresin yielding capability and growth among families and early selection for superior families in Pinus massoniana［J］.Forest Research，2013，26（6）：686-691.
12	张长波，孙红霞，巩中军，等.植物萜类化合物的天然合成途径及其相关合酶［J］.植物生理学通讯，2007，43（4）：779-786.
	Zhang C B，Sun H X，Gong Z J，et al.Plant terpenoid natural metabolism pathways and their synthases［J］.Plant Physiology Communications，2007，43（4）：779-786.
13	王毅，原晓龙，陈伟，等.高产脂思茅松牻牛儿基牻牛儿基焦磷酸合成酶基因的克隆及表达［J］.福建农林大学学报：自然科学版，2018，47（6）：711-716.
	Wang Y，Yuan X L，Chen W，et al.Identification and gene expression analysis of geranylgeranyl diphosphate synthase gene from high resin yield Pinus kesiya var.langbianensis［J］.Journal of Fujian Agriculture and Forestry University：Natural Science Edition，2018，47（6）：711-716.
14	刘美佳，于怡琳，姜森，等.珠子参中法尼基焦磷酸合酶（FPS）对皂苷生物合成的影响研究［J］.植物研究，2018，38（4）：611-618.
	Liu M J，Yu Y L，Jiang S，et al.Effect of Farnesyl-Pyrophosphate Synthase（FPS） on the biosynthesis of saponins in Panax japonicus［J］.Bulletin of Botanical Research，2018，38（4）：611-618.
15	杨悦，杨长青，王勇，等.过表达丹参GGPPS研究丹参酮类的生物合成途径［J］.中国科学：生命科学，2017，47（5）：586-594.
	Yang Y，Yang C Q，Wang Y，et al.Investigation of the tanshinone biosynthesis pathway by over expressing GGPPS in Salvia miltiorrhiza［J］.Scientia Sinica Vitae，2017，47（5）：586-594.
16	Bai Q S，He B X，Cai Y L，et al.Transcriptomic and metabolomic analyses reveal several critical metabolic pathways and candidate genes involved in resin biosynthesis in Pinus massoniana［J］.Molecular Genetics and Genomics，2020，295（2）：327-341.
17	梁宗锁，方誉民，杨东风，等.植物萜类化合物生物合成与调控及其代谢工程研究进展［J］.浙江理工大学学报：自然科学版，2017，37（2）：255-264.
	Liang Z S，Fang Y M，Yang D F，et al.Biosynthesis，regulation and metabolic engineering of terpenoids in Plants［J］.Journal of Zhejiang Sci-Tech University：Natural Sciences，2017，37（2）：255-264.
18	Hamberger B，Ohnishi T，Hamberger B，et al.Evolution of diterpene metabolism：Sitka spruce CYP720B4 catalyzes multiple oxidations in resin acid biosynthesis of conifer defense against insects［J］.Plant Physiology，2011，157（4）：1677-1695.
19	Hamberger B，Bohlmann J.Cytochrome P450 mono-oxygenases in conifer genomes：discovery of members of the terpenoid oxygenase superfamily in spruce and pine［J］.Biochemical Society Transactions，2006，34（6）：1209-1214.
20	Geisler K，Jensen N B，Yuen M M S，et al.Modularity of conifer diterpene resin acid biosynthesis：P450 enzymes of different CYP720B clades use alternative substrates and converge on the same products［J］.Plant Physiology，2016，171（1）：152-164.
21	Westbrook J W，Resende Jr M F R，Munoz P，et al.Association genetics of oleoresin flow in loblolly pine：discovering genes and predicting phenotype for improved resistance to bark beetles and bioenergy potential［J］.New Phytologist，2013，199（1）：89-100.
22	Liu Q H，Zhou Z C，Wei Y C，et al.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［J］.PLoS One，2015，10（7）：e0132624.
23	陈晓明.马尾松产脂相关基因挖掘及表达规律研究［D］.北京：中国林业科学研究院，2018.
	Chen X M.Study on discovery and expression pattern of the genes related to rosin formation in Pinus Massoniana［D］.Beijing：Chinese Academy of Forestry，2018.
24	李海霞，张志毅，张谦，等.毛白杨TIR-NBS-LRR基因转化烟草的研究［J］.北京林业大学学报，2009，31（1）：73-78.
	Li H X，Zhang Z Y，Zhang Q，et al.Transformation of tobaccos with a TIR-NBS-LRR gene isolated from Populus tomentosa Carr.［J］.Journal of Beijing Forestry University，2009，31（1）：73-78.
25	张颖，李峰，刘崇怀，等.中国野生刺葡萄抗白腐病NBS-LRR类抗病基因同源序列的分离与鉴定［J］.中国农业科学，2013，46（4）：780-789.
	Zhang Y，Li F，Liu C H，et al.Isolation and identification of NBS-LRR resistance gene analogs sequences from Vitis davidi［J］.Scientia Agricultura Sinica，2013，46（4）：780-789.
26	王猛.马尾松抑制消减文库的构建及抗病性相关基因的克隆［D］.长沙：中南林业科技大学，2010.
	Wang M.Construction of suppression subtracted cDNA Library of Pinus massoniana and cloning of disease resistance genes related to pine wilt diseas［D］.Changsha：Central South University of Forestry and Technology，2010.
27	祝钦泷，郭铁英，眭顺照，等.野生金荞麦高黄酮含量愈伤组织中一个类异黄酮还原酶（IRL）基因的克隆与分析［J］.药学学报，2009，44（7）：809-819.
	Zhu Q L，Guo T Y，Sui S Z，et al.Molecular cloning and characterization of a novel isoflavone reductase-like gene（FcIRL）from high flavonoids-producing callus of Fagopyrum cymosum［J］.Acta Pharmaceutica Sinica，2009，44（7）：809-819.
28	张美萍，王义，孙春玉，等.辐照西洋参培养物皂苷次生代谢调控的研究［J］.核农学报，2003，17（3）：207-211.
	Zhang M P，Wang Y，Sun C Y，et al.Study on the secondary metabolic regulation of callus of Panax quinquefolium L.［J］.Acta Agriculturae Nucleatae Sinica，2003，17（3）：207-211.

	产脂量 Yield of oleoresin(g)	树高 Tree height (m)	胸径 DBH (cm)	材积 Volume (m³)
产脂量 Yield of oleoresin(g)	1
树高 Tree height(m)	0.204^*	1
胸径 DBH(cm)	0.389^**	0.791^**	1
材积 Volume(m³)	0.425^**	0.835^**	0.980^**	1

		平方和SS	自由度df	均方MS	F	显著性p
材积 Volume(m³)	组间Between groups	0.021	1	0.021	3.909	0.076
	组内Within groups	0.054	10	0.005
	总计Total	0.075	11
产脂量 Oleoresin(g)	组间Between groups	3 290 302.413	1	3 290 302.413	176.704	0.000
	组内Within groups	186 203.927	10	186 203.927
	总计Total	3 476 506.340	11

样品 Sample	Raw Reads	Clean Reads	Clean Bases	错误率 Error (%)	Q20含量 Q20(%)	Q30含量 Q30(%)	GC含量 GC content(%)
H1	53680890	51782426	7.77G	0.03	97.65	93.34	44.79
H2	53249060	51725948	7.76G	0.03	97.75	93.56	44.71
H3	60645330	58618334	8.79G	0.03	97.69	93.43	44.86
H4	47981110	46300100	6.95G	0.03	97.66	93.37	44.56
H5	48513474	46527430	6.98G	0.03	97.60	93.20	44.57
H6	46606048	44244760	6.64G	0.03	97.60	93.24	44.22
L1	51182036	49598284	7.44G	0.03	97.75	93.56	45.05
L2	48064268	46743400	7.01G	0.03	97.50	93.04	44.85
L3	44098940	42197964	6.33G	0.03	97.46	92.96	45.09
L4	58865662	54092934	8.11G	0.03	97.55	93.16	44.54
L5	52851222	51255740	7.69G	0.03	97.55	93.14	44.95
L6	52687308	51021062	7.65G	0.03	97.48	92.97	44.74

注释类别 Annotation category	基因数 Number of genes	百分比 Percentage(%)
注释到NR数据库 Annotated in NR	87 608	45.98
注释到NT数据库 Annotated in NT	110 155	57.81
注释到KO数据库 Annotated in KO	38 361	20.13
注释到Swiss-Prot数据库 Annotated in Swiss-Prot	86 442	45.37
注释到PFAM数据库 Annotated in PFAM	91 373	47.96
注释到GO数据库 Annotated in GO	91 373	47.96
注释到KOG数据库 Annotated in KOG	23 382	12.27
注释到所有数据库 Annotated in all Databases	10 265	5.38
注释到至少一个数据库 Annotated in at least one Database	144 023	75.59
总基因数 Total unigenes	190 514	100.00

基因编号 Gene ID	FC对数值 log FC	q值 q value	描述 Description
Cluster-37595.96613	5.153 8	0.0091454	ABC transporter C family member
Cluster-37595.128007	-9.396 4	0.0099982	ABC transporter
Cluster-37595.77181	22.432	0.000000000092327	ATP binding
Cluster-37595.89400	24.203	0.0000000000010724	ABC transporter
Cluster-37595.20113	22.805	0.000000000048995	ABC transporter B family member
Cluster-37595.127923	-22.481	0.000000000084192	ADP binding
Cluster-36047.0	25.034	0.00000000000064934	Cytochrome P450
Cluster-31181.0	22.569	0.000000000074589	Cytochrome P450
Cluster-14555.0	4.897	0.0018662	Cytochrome P450