-
合理的林分密度能够显著促进林木的径向生长[1-2],从而提高木材产量,解析种植密度促进林木径向生长的内在机制意义重大。林木的径向生长(即次生木质部生成)来源于维管形成层的活动。维管形成层细胞向外分化产生次生韧皮部,向内分化产生次生木质部,其中,次生木质部细胞不断积累,使植株径向增粗[3]。由维管形成层分化产生成熟的次生木质部细胞需要经历木质部母细胞分化和分裂、木质部细胞增大、次生壁加厚、细胞程序化凋亡和心材形成多个阶段[4],各个阶段在分子水平上均受到一大批基因的调控。因此,研究木质部在不同种植密度下的基因表达模式,挖掘和鉴定响应种植密度的木质部生成关键基因并解析其作用机制是研究种植密度影响林木径向生长内在机制的重要内容。
以往以模式植物拟南芥(Arabidopsis thaliana (Linn.) Heynh.)和杨属(Populus)为研究对象,分离鉴定了一些参与木质部生成各个阶段的关键调控基因,在木质部生成方面取得了重要进展。诸多研究表明,调控因子如受体激酶和转录因子等在木质部生成过程中发挥重要的调控作用,如PTL(PETAL LOSS,花瓣缺失)转录因子作为抑制子参与调控拟南芥根系发育过程中的形成层细胞分裂和木质部细胞分化活动[5];类受体激酶PXY和其家族成员PXL1、PXL2参与调控拟南芥维管组织的细胞分裂活动[6-7];类受体激酶FEI和SHOU参与次生壁物质合成,fei突变体纤维素合成减少,而SHOU能够恢复fei突变体表型[8-9]。次生壁合成也需要一些转录因子参与,Chen等[10]在毛果杨(Populus trichocarpa Torr. & Gray)中绘制了一个4层级的细胞壁合成转录调控网络,其中包含17个转录因子和27个细胞壁成分合成基因,它们协同调控木质部生成;最近在杨属、松属(Pinus)等树种中又相继证实了多个NAC和MYB家族成员参与木材细胞壁物质合成过程[11-13]。由此可见,次生木质部生成的分子机制十分复杂,仍然有许多未知亟待解析。
桉属(Eucalyptus)与杨属为世界公认的速生人工林树种,但二者木材的生物量和质量存在明显差异。这2个树种在木质部生成调控机制上可能存在差异,但目前关于桉属植物木质部生成相关重要基因的解析及其调控机制的研究相对较少[14-15],同时,关于种植密度影响径向生长的分子基础研究尚未见报道。因此,本研究以不同种植密度下尾巨桉(Eucalyptus urophylla × Eucalyptus grandis W.Hill ex Maiden)的主茎为试验材料,联合PacBio Iso-Seq和RNA-Seq转录组技术,通过对尾巨桉木质部差异表达转录本(DETs)的表达模式展开分析,挖掘和鉴定桉树中响应种植密度的木质部生成关键基因,以期为科学阐明种植密度影响林木径向生长的分子机制奠定理论基础。
-
试验地位于广东省湛江市南方国家级林木种苗示范基地内(111°38′ E,21°30′ N)。试验材料为长势均一、主茎通直、无病虫害的尾巨桉DH32-29无性系组织培育苗。2019年7月14日,将幼苗移植至高0.5 m直径0.3 m的种植袋中,并按照10株·m−2高种植密度(株行距0.3 m × 0.3 m)和5株·m−2低种植密度(株行距0.6 m × 0.6 m)移放至长 × 宽 × 高为5.0 m × 4.5 m × 0.5 m的样地内,每个处理设置3个重复,共计6个样地。在幼苗生长过程中,水分和营养充分供应且保持一致。苗木种植后定期测定生长性状,于直径及其增长率表现出显著差异的时期进行取样(2019年10月15日)(表1)。每个样地选择1株平均木,取距离地面10 cm处主茎的木质部、韧皮部,共计12个样本。液氮冷冻后,带回试验室置于−80℃的超低温冰箱保存备用。样品命名为H-X1~3、L-X1~3、H-P1~3、L-P1~3,其中,H为高种植密度,L为低种植密度,X表示木质部,P表示韧皮部,数字1~3代表生物学重复。
表 1 种植密度对尾巨桉幼苗生长性状的影响
Table 1. Effect of planting density on growth traits for seedlings of E. urophylla × E. grandis
处理
Treatment平均直径/mm
Average diameter直径生长量/mm
Diameter increment平均株高/cm
Average height株高生长量/cm
Height increment低种植密度 Low planting density 6.84 ± 0.19 a 2.36 ± 0.28 a 65.72 ± 2.64 a 18.80 ± 2.51 a 高种植密度 High planting density 5.41 ± 0.20 b 1.41 ± 0.15 b 58.53 ± 2.21 a 13.97 ± 2.92 a 注:不同小写字母表示不同种植密度间的生长性状具有显著性差异(P=0.01)。
Note: The different letters indicate growth traits with significant differences in different planting densities(P=0.01). -
使用植物RNA提取试剂盒DP130508(天根生物公司)分别提取12个样本的RNA,等量混合后使用SMARTer™ PCR cDNA Synthesis Kit试剂盒合成mRNA的全长cDNA,通过PacBio平台进行全长转录组测序(百迈客公司);使用CD-HIT软件[16]去除转录本中的冗余序列;使用BLAST软件(version 2.2.26)[17]将得到非冗余转录本序列和NR、eggNOG、Pfam、GO、Swissprot、KOG、COG、KEGG数据库比对进行转录本功能注释;使用iTAK软件[18]进行调控因子预测。
对12个样本的RNA分别进行cDNA合成和文库构建,委托百迈客公司进行Illumina二代测序。使用RSEM[19]软件,通过比对到三代转录本上的位置信息,对转录本的表达水平进行定量,采用FPKM(Fragments Per Kb of transcript per Million fragments mapped)作为衡量指标。
-
使用DESeq2[20]对样本H-X1、H-X2、H-X3和L-X1、L-X2、L-X3进行分析,获得2个种植密度间DETs。筛选标准为差异倍数(Fold Change) ≥ 2且错误发现率(FDR)< 0.01;使用百迈客云平台热图绘制工具绘制DETs 热图,每个表达量的数值均取以10为底数的对数(Log10)[21];使用STRING数据库[22]对DETs 构建蛋白互作网络。
-
以转录组试验中2个种植密度下尾巨桉幼苗的木质部RNA为模板,使用TSK302S RT6 cDNA Synthesis Kit Ver 2(擎科生物公司)合成cDNA模板,从443个DETs中随机选取9个DETs,进行qRT-PCR验证转录组测序结果的准确性;选择9个DETs进行组织表达分析,材料为尾巨桉幼苗的叶片、韧皮部、第2节间木质部(木质部2)和第7节间木质部(木质部7)4个部位,RNA的提取采用植物RNA提取试剂盒V1.5(百菲特生物公司)。利用Primer 5.0设计qRT-PCR引物,以arp4(Actin-related protein 4)编码基因作为内参,引物见表2。使用2 × T5 Fast qPCR Mix(SYBR Green I)荧光定量试剂盒(擎科生物公司)进行荧光定量检测。每个样品设置4个技术重复,采用2−ΔΔCt计算相对表达量。使用SPSS 25.0处理数据,使用Excel 2016软件绘图,以P < 0.01表示差异显著。
表 2 qRT-PCR引物序列
Table 2. Primers for qRT-PCR
序号
Number转录本编号
Transcripts ID上游引物
Forward primer下游引物
Reverse primer1 MIX-PB_transcript_46917 GGAGATGAAGAATGGCGTGC CAGCGGAAGAACTCGCATAC 2 MIX-PB_transcript_39150 GGGACTAAAATTATCGGATTACGGT GCATAGTCTGTTGATGTTGGTAGCA 3 MIX-PB_transcript_72449 GTTCCTGGTCTCCGACAGTT AGGCAAGATACAGCACCCTG 4 MIX-PB_transcript_15209 CTCTGAAATCGGTGACTTCCTC AAGCATCTTTGCCACCCAAC 5 MIX-PB_transcript_17152 TCGTGGATGCGTTGCTTACT GGGCTCTCGGGTTCCTAATG 6 MIX-PB_transcript_61822 CAAACGAGTGGGAGCAAGGA CGGATCTCAATGCCATGCTG 7 MIX-PB_transcript_39624 TCGAGCACGCATTTTCTGTTT GGGTCGCAGGTGGTTCTATT 8 MIX-PB_transcript_6405 GTTGGATTTGGGTCGTGTCG GAGACTAGACACACCGGACAG 9 MIX-PB_transcript_6089 CCTTCGCTAAGAACAATCCACG AGCTTGCGGGTACCATGTCATA 10 MIX-PB_transcript_67597 CGACCAGATTCCAGGAGAGC GCCAGATCGAGACACCTTAAACT 11 MIX-PB_transcript_18228 GCTTTCTCCCCACAAGAGGA TCGTGTTTGGCGAGTGTCTA 12 MIX-PB_transcript_51497 GCCACCAAAGGACCGGAATA ACTTCTCTTGGCACTCACCG 13 MIX-PB_transcript_15184 GGCGGGAAAAGAAAAGCACG GCGAAATCAGAGACGCTCAC 14 MIX-PB_transcript_70069 CCTGATGATACGGGCTTCGT ATGATCCGTCACTCTTGCTGT 15 MIX-PB_transcript_73425 AGCTGGAATGTTGTAGTCGCT CTCAAGGTCCACATCTGCCAA 16 MIX-PB_transcript_12050 GTTTCCTGGCCGATAACCCA AGGAGGAGGGGTCAACACTT 17 MIX-PB_transcript_52911 GGTCATGGAGTGTGTCAAGC TTCCCACGCCATTGTAAGAGA 18 MIX-PB_transcript_30100 TGCTGATGCTTGACAGACAG AAGAAGCACCCACAATCCAC 19 arp4 GGAATTGTCGCGGATTGGGA GCCTTTTCTCGTTGCTGCTG -
对高、低2个种植密度下尾巨桉主茎木质部和韧皮部进行Pacbio全长转录组测序,共获得247 536条环形一致序列(CCS),其中,209 654条为全长非嵌合序列;对全长非嵌合序列进行聚类分析,获得81 317条一致性序列;进一步对一致性序列进行校正,最终获得高质量转录本(准确度>99%)79 965条,低质量转录本1 281条,高质量转录本达到98%以上;经CD-HIT去冗余后,获得非冗余转录本序列45 490条(表3)。
表 3 PacBio Iso-Seq测序数据评估统计表
Table 3. Statistical table of PacBio Iso-Seq data evaluation
分类
Classification统计
CountCCS数量
Number of CCS247 536 全长非嵌合序列数量
Number of full-length non-chimeric reads209 654 一致性序列数量
Number of consensus isoforms81 317 高质量转录本数量
Number of polished high-quality isoforms79 965 低质量转录本数量
Number of polished low-quality isoforms1 281 高质量转录本所占比例
Percent of polished high-quality isoforms/%98.34 非冗余转录本数
Number of non-redundant transcripts45 490 对高、低两个种植密度下尾巨桉主茎木质部进行RNA-Seq测序,共获得130 292 914条150 bp双端Reads。6个样本中,测序获得最多的Clean Reads为23 322 440条,最少为19 998 315条;获得Clean Data碱基数最多为6 988 132 460 bp,最少为5 990 895 036 bp;GC含量为50.30%~51.05%,%≥Q30(Clean Data质量值≥30的碱基所占的百分比)为93.88%~94.41%(表4),GC含量和%≥Q30均值分别为50.60%、94.19%。
表 4 RNA-Seq测序数据评估统计表
Table 4. Statistical table of RNA-Seq data evaluation
样品
SamplesClean Reads数量
Number of Clean ReadsClean Data碱基数/bp
Number of Clean DataGC含量/%
GC Content大于或等于Q30的百分比/%
≥ Q30H-X1 21 882 942 6 557 229 728 50.52 94.09 H-X2 22 886 461 6 843 134 076 50.81 94.41 H-X3 21 093 958 6 317 393 546 51.05 94.40 L-X1 19 998 315 5 990 895 036 50.42 94.21 L-X2 21 108 798 6 323 257 814 50.52 94.16 L-X3 23 322 440 6 988 132 460 50.30 93.88 共计/均值
Total/Mean130 292 914 39 020 042 660 50.60 94.19 -
将三代转录组测序获得的45 490个非冗余转录本与NR、eggNOG、Pfam、GO、Swissprot、KOG、COG、KEGG数据库比对并进行功能注释,有44 888条转录本至少比对到一个数据库中,注释率达98.68%(表5),其中,NR数据库注释率达98.47%,eggNOG、Pfam、Swissprot、GO、KOG数据库注释到的转录本数量均超过50%,KEGG数据库注释到的转录本数量最低。
表 5 转录本在各数据注释情况
Table 5. Statistics of BLAST annotation of different databases
数据库
Annotated databases注释到的转录本数目
Number of transcripts注释率/%
Annotation rateNR 44854 98.47 eggNOG 43433 95.47 Pfam 37496 82.43 GO 36631 80.52 Swissprot 33276 73.15 KOG 27767 61.04 COG 19986 44.93 KEGG 17093 37.58 All 44888 98.68 -
为了充分了解不同种植密度下尾巨桉木质部转录本的表达模式,以H-X为对照组,L-X为试验组,筛选获得差异表达转录本443个,其中,215个DETs在低种植密度条件下上调表达,228个DETs下调表达。注释结果显示,443个DETs中包含了60个假定的调控因子,归为37个家族(图1A、图2),其中,28个差异表达调控因子在低种植密度下尾巨桉的木质部中上调表达,编码21类调控因子。这些调控因子与多种生物学过程有关,GARP-ARR-B、Trihelix和RLK-Pelle_LRR-XI-1参与调控维管形成层活动;C3H、MYB、NAC和RLK-Pelle_LRR-XIIIa参与调控细胞壁成分合成;RLK-Pelle_LRR-III参与次生生长过程;RLK-Pelle_RLCK-V参与细胞周期调控;GARP-G2-like和WNK_NRBP参与节律调控。32个DETs下调表达,所对应调控因子的功能多与外界胁迫有关,如热激转录因子HSF、参与盐胁迫的bHLH和STE_STE7、参与干旱胁迫的RLK-Pelle_WAK_LRK10L-1、参与病虫害应答的RLK-Pelle_LRR-VI-2和RLK-Pelle_LRR-Xa、对脱落酸(ABA)敏感的bHLH和RLK-Pelle_CrRLK1L-1等。除编码调控因子的转录本外,一些转录本的表达模式也发生了显著的变化,在低种植密度植株的木质部中细胞壁物质合成相关蛋白(蔗糖合成酶,纤维素合成酶,漆酶,糖基转移酶,糖基水解酶)和木质素物质合成相关编码基因多为上调表达(图1B)。
-
通过STRING数据库对443个DETs进行蛋白互作网络分析,检测到一簇相互作用的基因簇,共包含30个转录本(图3,附表1)。这些基因主要与参与翻译后修饰、蛋白质转换和伴侣蛋白,细胞周期控制、细胞分裂和染色体分割,转录等功能。MIX-PB_transcript_67597作为网络的核心基因,与MIX-PB_transcript_30100和MIX-PB_transcript_52911占据了网络的中心位置(图3)。数据库注释表明:MIX-PB_transcript_67597编码类受体激酶家族RLK-Pelle_LRR-XI-1的PXL2,注释功能为转录。转录本MIX-PB_transcript_52911编码蛋白CUL1参与目标蛋白的泛素化;MIX-PB_transcript_30100编码蛋白T15D22.7参与硫醇氧化酶和硫氧还蛋白家族蛋白的合成,二者在数据库的注释功能均为细胞周期控制、细胞分裂和染色体分割。这3个基因在低种植密度植株的木质部中均上调表达,可能在尾巨桉主茎径向生长中发挥关键作用(附表1)。
图 3 不同种植密度下尾巨桉木质部中差异表达转录本蛋白互作网络
Figure 3. The interaction network of proteins of differentially expressed transcripts in xylem for E. urophylla × E. grandis under different planting densities
表 1 蛋白互作网络中30个参与种植密度响应的差异表达转录本注释
Table 1. The annotation of 30 differentially expressed transcripts responding to planting density in protein interaction network
转录本 Transcripts ID 转录本的功能 Transcript function 注释Annotation MIX-PB_transcript_67597 转录
Transcription富含亮氨酸的重复类受体蛋白激酶PXL2(前体)
Leucine-rich repeat receptor-like protein kinase PXL2 (Precursor)MIX-PB_transcript_52911 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning类枯灵素蛋白1
cullin-1-like ,CUL1MIX-PB_transcript_30100 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning巯基氧化酶1(前体)
Sulfhydryl oxidase 1 (Precursor), T15D22.7MIX-PB_transcript_31553 -- 推定的抗病蛋白RGA4 Putative disease resistance protein RGA4 MIX-PB_transcript_4630 -- Toll/白细胞介素-1受体-类蛋白
Toll/interleukin-1 receptor-like protein, TIRMIX-PB_transcript_41125 -- 可能的无活性富亮氨酸重复受体类蛋白激酶
Probably inactive leucine-rich repeat receptor-like protein kinase At3g28040MIX-PB_transcript_67116 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_66737 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_28715 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_32439 信号转导机制
Signal transduction mechanisms可能的蛋白磷酸酶2C60
Probable protein phosphatase 2CMIX-PB_transcript_41236 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones类硫氧还蛋白1-1,叶绿体(前体)
Thioredoxin-like 1-1, chloroplastic (Precursor)MIX-PB_transcript_53209 -- 金属内蛋白酶1(前体)
Metalloendoproteinase 1 (Precursor)MIX-PB_transcript_5308 -- 枯草杆菌蛋白酶样蛋白酶(前体)
Subtilisin-like protease (Precursor)MIX-PB_transcript_29586 脂质运输和代谢Lipid transport and metabolism C型1类尼曼-匹克蛋白
Niemann-Pick C1 protein isoform X1MIX-PB_transcript_48714 一般功能预测General function prediction only 锚定重复序列蛋白
Ankyrin repeat-containing protein, At5g02620MIX-PB_transcript_59534 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_8322 转录Transcription 热激蛋白83
Heat shock protein 83MIX-PB_transcript_76091 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones热激同源蛋白70
Heat shock cognate 70 kDa protein HSP70MIX-PB_transcript_63092 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones热激同源蛋白70
Heat shock cognate 70 kDa protein HSP70MIX-PB_transcript_12904 一般功能预测General function prediction only 多元醇转运蛋白5
Polyol transporter 5, PLT5MIX-PB_transcript_32159 转录Transcription 热应激转录因子A-2 T19L18.4
Heat stress transcription factor A-2 T19L18.4MIX-PB_transcript_26893 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_53117 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_44896 氨基酸转运和代谢
Amino acid transport and metabolism类酪氨酸脱羧酶1
tyrosine decarboxylase 1-likeMIX-PB_transcript_17956 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesDnaJ同源蛋白(前体)
DnaJ protein homolog (Precursor), DNAJ1MIX-PB_transcript_7819 一般功能预测
General function prediction onlyF框/亮氨酸精氨酸精氨酸重复蛋白3
F-box/LRR-repeat protein 3, FBL3MIX-PB_transcript_44207 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesE3泛素蛋白连接酶环1
E3 ubiquitin-protein ligase RING1MIX-PB_transcript_13679 一般功能预测General function prediction only 顶端向地性蛋白5 Protein SHOOT GRAVITROPISM 5 MIX-PB_transcript_10365 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning可能的丝氨酸/苏氨酸蛋白激酶
Probable serine/threonine-protein kinase, At1g54610MIX-PB_transcript_71371 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesE3泛素蛋白连接酶环1
E3 ubiquitin-protein ligase RING1 -
以转录组测序中2个种植密度下的木质部RNA为模板,随机选择9个DETs进行qRT-PCR验证,qRT-PCR结果与转录组测序的表达结果基本一致(图4),证明转录组测序数据准确可靠。
图 4 9个差异表达转录本的qRT-PCR验证
Figure 4. Confirmation of the expression profiles of 9 differentially expressed transcripts by qRT-PCR
为了验证木质部生成候选基因的功能,通过qRT-PCR对9个木质部生成候选基因进行分析,检测了它们在不同组织的基因表达情况。这9个DETs分别是转录因子基因PTL、MYB46、C3H14、NAC86,细胞壁合成酶基因CesA、LAC17以及PXL2、T15D22.7、CUL1。qRT-PCR结果(图5):这些编码基因在叶片中几乎不表达,在韧皮部、木质部均有不同程度的表达,其中,PXL2在木质部2中的表达量最高,在韧皮部和木质部7中表达水平也较高;CUL1和T15D22.7在韧皮部和木质部中高水平表达;PTL在韧皮部和木质部7中表达量最高,显著高于叶片和木质部2;NAC86在韧皮部中表达水平最高,但在木质部2和木质部7中也有较高表达水平;MYB46、C3H14、CesA和LAC17的表达模式相对一致,在木质部2的表达量最高,木质部7中略低,而韧皮部中表达量极低。这些结果进一步证明了测序结果的可靠性,同时也表明这些基因可能在木质部生成中发挥重要作用。
基于转录水平解析尾巨桉径向生长对种植密度的响应
Response of Stem Radial Growth of Eucalyptus urophylla × E. grandis to Planting Density Based on Transcriptome Analysis
-
摘要:
目的 挖掘和鉴定桉树中响应种植密度的木质部生成关键基因,阐释种植密度影响林木径向生长的分子机制。 方法 通过PacBio Iso-Seq和RNA-Seq联合分析,鉴定了在高、低2个种植密度水平下尾巨桉DH32-29主茎的差异表达转录本并利用qRT-PCR进行组织表达分析。 结果 分别获得45 490条在主茎中表达的非冗余全长转录本和443条在木质部中差异表达的转录本,并鉴定出60条差异表达的调控因子编码基因。在低种植密度条件下,植株直径生长量显著增加;参与细胞分裂调控的PXL2及其互作基因CUL1和T15D22.7、参与次生壁调控的MYB46、C3H14同源基因和其下游基因CesAs、LACs均上调表达,这些基因很可能是种植密度影响主茎径向生长过程中负责调控细胞分裂和次生壁合成的重要成员;参与筛管发育的NAC86同源基因及参与形成层发育和木质部分化的抑制子PTL同源基因也上调表达,它们可能促进木质部生成,这与草本植物中的功能不同。组织表达分析结果显示:PXL2、CUL1、T15D22.7、NAC86和PTL在韧皮部和木质部中优势表达;MYB46、C3H14、CesA和LAC17在木质部中优势表达。 结论 本研究鉴定了一批参与种植密度响应的木质部生成候选基因,构建了种植密度影响尾巨桉径向生长的分子调控网络模型,研究结果为深入解析种植密度影响林木径向生长的分子机制奠定基础。 Abstract:Objective To identify the key genes of secondary xylem development response to planting density for a well understanding of molecular mechanism of planting density affecting radial growth of eucalypts. Method By a combination of PacBio Iso-Seq and RNA-Seq analysis, the differentially expressed transcriptomes of xylem cells in Eucalyptus urophylla × E. grandis were identified under high and low planting densities. The tissues expression profiles of these key genes were analyzed via qRT-PCR. Result A total of 45 490 non-redundant full-length transcripts and 443 transcripts differentially expressed in xylem cells were obtained under high and low planting densities, and 60 transcripts encoding regulatory factors were obtained. Under low planting density, the diameters of trees increased significantly. The PXL2 and its interactional genes CUL1, T15D22.7 related to cell division, the MYB46, C3H14 with their downstream genes CesAs and LACs related to secondary wall regulation were preferentially expressed in the xylem cells. These genes might play key roles in the regulation of diameter growth under different densities. In addition, the NAC86 homologous genes involved in sieve element development and the inhibitor PTL homologous genes with dual functions in cambial cell proliferation and xylem differentiation were also up-regulated. They could promote the xylem development, which were different from the functions in herbaceous plants. The results of tissue expression analysis showed that PXL2, CUL1, T15D22.7, NAC86 and PTL were predominantly expressed in phloem and xylem, whereas MYB46, C3H14, CesA and LAC17 were predominantly expressed in xylem. Conclusion In this study, the candidate genes of xylem development related to planting density are identified and a model of molecular regulatory network that how the planting density affects radial growth of E. urophylla × E. grandis is proposed, which will benefit the intensive study of the molecular mechanism under different planting densities affecting radial growth for trees. -
表 1 种植密度对尾巨桉幼苗生长性状的影响
Table 1. Effect of planting density on growth traits for seedlings of E. urophylla × E. grandis
处理
Treatment平均直径/mm
Average diameter直径生长量/mm
Diameter increment平均株高/cm
Average height株高生长量/cm
Height increment低种植密度 Low planting density 6.84 ± 0.19 a 2.36 ± 0.28 a 65.72 ± 2.64 a 18.80 ± 2.51 a 高种植密度 High planting density 5.41 ± 0.20 b 1.41 ± 0.15 b 58.53 ± 2.21 a 13.97 ± 2.92 a 注:不同小写字母表示不同种植密度间的生长性状具有显著性差异(P=0.01)。
Note: The different letters indicate growth traits with significant differences in different planting densities(P=0.01).表 2 qRT-PCR引物序列
Table 2. Primers for qRT-PCR
序号
Number转录本编号
Transcripts ID上游引物
Forward primer下游引物
Reverse primer1 MIX-PB_transcript_46917 GGAGATGAAGAATGGCGTGC CAGCGGAAGAACTCGCATAC 2 MIX-PB_transcript_39150 GGGACTAAAATTATCGGATTACGGT GCATAGTCTGTTGATGTTGGTAGCA 3 MIX-PB_transcript_72449 GTTCCTGGTCTCCGACAGTT AGGCAAGATACAGCACCCTG 4 MIX-PB_transcript_15209 CTCTGAAATCGGTGACTTCCTC AAGCATCTTTGCCACCCAAC 5 MIX-PB_transcript_17152 TCGTGGATGCGTTGCTTACT GGGCTCTCGGGTTCCTAATG 6 MIX-PB_transcript_61822 CAAACGAGTGGGAGCAAGGA CGGATCTCAATGCCATGCTG 7 MIX-PB_transcript_39624 TCGAGCACGCATTTTCTGTTT GGGTCGCAGGTGGTTCTATT 8 MIX-PB_transcript_6405 GTTGGATTTGGGTCGTGTCG GAGACTAGACACACCGGACAG 9 MIX-PB_transcript_6089 CCTTCGCTAAGAACAATCCACG AGCTTGCGGGTACCATGTCATA 10 MIX-PB_transcript_67597 CGACCAGATTCCAGGAGAGC GCCAGATCGAGACACCTTAAACT 11 MIX-PB_transcript_18228 GCTTTCTCCCCACAAGAGGA TCGTGTTTGGCGAGTGTCTA 12 MIX-PB_transcript_51497 GCCACCAAAGGACCGGAATA ACTTCTCTTGGCACTCACCG 13 MIX-PB_transcript_15184 GGCGGGAAAAGAAAAGCACG GCGAAATCAGAGACGCTCAC 14 MIX-PB_transcript_70069 CCTGATGATACGGGCTTCGT ATGATCCGTCACTCTTGCTGT 15 MIX-PB_transcript_73425 AGCTGGAATGTTGTAGTCGCT CTCAAGGTCCACATCTGCCAA 16 MIX-PB_transcript_12050 GTTTCCTGGCCGATAACCCA AGGAGGAGGGGTCAACACTT 17 MIX-PB_transcript_52911 GGTCATGGAGTGTGTCAAGC TTCCCACGCCATTGTAAGAGA 18 MIX-PB_transcript_30100 TGCTGATGCTTGACAGACAG AAGAAGCACCCACAATCCAC 19 arp4 GGAATTGTCGCGGATTGGGA GCCTTTTCTCGTTGCTGCTG 表 3 PacBio Iso-Seq测序数据评估统计表
Table 3. Statistical table of PacBio Iso-Seq data evaluation
分类
Classification统计
CountCCS数量
Number of CCS247 536 全长非嵌合序列数量
Number of full-length non-chimeric reads209 654 一致性序列数量
Number of consensus isoforms81 317 高质量转录本数量
Number of polished high-quality isoforms79 965 低质量转录本数量
Number of polished low-quality isoforms1 281 高质量转录本所占比例
Percent of polished high-quality isoforms/%98.34 非冗余转录本数
Number of non-redundant transcripts45 490 表 4 RNA-Seq测序数据评估统计表
Table 4. Statistical table of RNA-Seq data evaluation
样品
SamplesClean Reads数量
Number of Clean ReadsClean Data碱基数/bp
Number of Clean DataGC含量/%
GC Content大于或等于Q30的百分比/%
≥ Q30H-X1 21 882 942 6 557 229 728 50.52 94.09 H-X2 22 886 461 6 843 134 076 50.81 94.41 H-X3 21 093 958 6 317 393 546 51.05 94.40 L-X1 19 998 315 5 990 895 036 50.42 94.21 L-X2 21 108 798 6 323 257 814 50.52 94.16 L-X3 23 322 440 6 988 132 460 50.30 93.88 共计/均值
Total/Mean130 292 914 39 020 042 660 50.60 94.19 表 5 转录本在各数据注释情况
Table 5. Statistics of BLAST annotation of different databases
数据库
Annotated databases注释到的转录本数目
Number of transcripts注释率/%
Annotation rateNR 44854 98.47 eggNOG 43433 95.47 Pfam 37496 82.43 GO 36631 80.52 Swissprot 33276 73.15 KOG 27767 61.04 COG 19986 44.93 KEGG 17093 37.58 All 44888 98.68 表 1 蛋白互作网络中30个参与种植密度响应的差异表达转录本注释
Table 1. The annotation of 30 differentially expressed transcripts responding to planting density in protein interaction network
转录本 Transcripts ID 转录本的功能 Transcript function 注释Annotation MIX-PB_transcript_67597 转录
Transcription富含亮氨酸的重复类受体蛋白激酶PXL2(前体)
Leucine-rich repeat receptor-like protein kinase PXL2 (Precursor)MIX-PB_transcript_52911 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning类枯灵素蛋白1
cullin-1-like ,CUL1MIX-PB_transcript_30100 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning巯基氧化酶1(前体)
Sulfhydryl oxidase 1 (Precursor), T15D22.7MIX-PB_transcript_31553 -- 推定的抗病蛋白RGA4 Putative disease resistance protein RGA4 MIX-PB_transcript_4630 -- Toll/白细胞介素-1受体-类蛋白
Toll/interleukin-1 receptor-like protein, TIRMIX-PB_transcript_41125 -- 可能的无活性富亮氨酸重复受体类蛋白激酶
Probably inactive leucine-rich repeat receptor-like protein kinase At3g28040MIX-PB_transcript_67116 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_66737 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_28715 -- 烟草花叶病毒抗性蛋白N
TMV resistance protein NMIX-PB_transcript_32439 信号转导机制
Signal transduction mechanisms可能的蛋白磷酸酶2C60
Probable protein phosphatase 2CMIX-PB_transcript_41236 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones类硫氧还蛋白1-1,叶绿体(前体)
Thioredoxin-like 1-1, chloroplastic (Precursor)MIX-PB_transcript_53209 -- 金属内蛋白酶1(前体)
Metalloendoproteinase 1 (Precursor)MIX-PB_transcript_5308 -- 枯草杆菌蛋白酶样蛋白酶(前体)
Subtilisin-like protease (Precursor)MIX-PB_transcript_29586 脂质运输和代谢Lipid transport and metabolism C型1类尼曼-匹克蛋白
Niemann-Pick C1 protein isoform X1MIX-PB_transcript_48714 一般功能预测General function prediction only 锚定重复序列蛋白
Ankyrin repeat-containing protein, At5g02620MIX-PB_transcript_59534 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_8322 转录Transcription 热激蛋白83
Heat shock protein 83MIX-PB_transcript_76091 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones热激同源蛋白70
Heat shock cognate 70 kDa protein HSP70MIX-PB_transcript_63092 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperones热激同源蛋白70
Heat shock cognate 70 kDa protein HSP70MIX-PB_transcript_12904 一般功能预测General function prediction only 多元醇转运蛋白5
Polyol transporter 5, PLT5MIX-PB_transcript_32159 转录Transcription 热应激转录因子A-2 T19L18.4
Heat stress transcription factor A-2 T19L18.4MIX-PB_transcript_26893 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_53117 转录Transcription 热激因子蛋白30
Heat shock factor protein HSF30MIX-PB_transcript_44896 氨基酸转运和代谢
Amino acid transport and metabolism类酪氨酸脱羧酶1
tyrosine decarboxylase 1-likeMIX-PB_transcript_17956 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesDnaJ同源蛋白(前体)
DnaJ protein homolog (Precursor), DNAJ1MIX-PB_transcript_7819 一般功能预测
General function prediction onlyF框/亮氨酸精氨酸精氨酸重复蛋白3
F-box/LRR-repeat protein 3, FBL3MIX-PB_transcript_44207 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesE3泛素蛋白连接酶环1
E3 ubiquitin-protein ligase RING1MIX-PB_transcript_13679 一般功能预测General function prediction only 顶端向地性蛋白5 Protein SHOOT GRAVITROPISM 5 MIX-PB_transcript_10365 细胞周期控制,细胞分裂,染色体分割
Cell cycle control, cell division, chromosome partitioning可能的丝氨酸/苏氨酸蛋白激酶
Probable serine/threonine-protein kinase, At1g54610MIX-PB_transcript_71371 翻译后修饰,蛋白质转换,伴侣蛋白
Posttranslational modification, protein turnover, chaperonesE3泛素蛋白连接酶环1
E3 ubiquitin-protein ligase RING1 -
[1] 李 洁, 列志旸, 许松葵, 等. 不同密度的银合欢林生长分析[J]. 中南林业科技大学学报, 2016, 36(6):70-74. [2] 楚秀丽, 王 艺, 金国庆, 等. 不同生境、初植密度及林龄木荷人工林生长、材性变异及林分分化[J]. 林业科学, 2014, 50(6):152-159. [3] Larson P R. The Vascular Cambium: Development and Structure[M]. Berlin: Springer-Verlag Berlin Heidelberg, 1994: 9-30. [4] Plomion C, LeprovostG, Stokes A. Wood formation in trees[J]. Plant Physiology, 2001, 127(4): 1513-1523. doi: 10.1104/pp.010816 [5] Zhang J, Eswaran G, Serra A J, et al. Transcriptional regulatory framework for vascular cambium development in Arabidopsis roots[J]. Nature Plants, 2019, 5(10): 1033-1042. doi: 10.1038/s41477-019-0522-9 [6] Fisher K, Turner S. PXY, a receptor-like kinase essential for maintaining polarity during plant vascular-tissue development[J]. Current Biology, 2007, 17(12): 1061-1066. doi: 10.1016/j.cub.2007.05.049 [7] Smit M E, Mcgregor S R, Sun H, et al. A PXY-mediated transcriptional network integrates signaling mechanisms to control vascular development in Arabidopsis[J]. The Plant Cell, 2019, 32(2): 319-335. [8] Xu S, Rahman A, Baskin I T. Two leucine-rich repeat receptor kinases mediate signaling, linking cell wall biosynthesis and ACC synthase in Arabidopsis[J]. The Plant Cell, 2008, 20(11): 3065-3079. doi: 10.1105/tpc.108.063354 [9] Polko K J, Barnes J W, Voiniciuc C, et al. SHOU4 proteins regulate traffificking of cellulose synthase complexes to the plasma membrane[J]. Current Biology, 2018, 28(19): 1-9. [10] Chen H, Wang J P, Liu H, et al. Hierarchical transcription factor and chromatin binding network for wood formation in Populus trichocarpa[J]. The Plant Cell, 2019, 31(3): 602-626. [11] Akiyoshi N, Nakano Y, Sano R, et al. Involvement of VNS NAC-domain transcription factors in tracheid formation in Pinus taeda[J]. Tree Physiology, 2019, 40(6): 1-13. [12] Sun Y, Ren S, Ye S, et al. Identification and functional characterization of PtoMYB055 involved in the regulation of the lignin biosynthesis pathway in Populus tomentosa[J]. International Journal of Molecular Sciences, 2020, 21(14): 4857. [13] Balmant K, Noble J, Alves F C, et al. Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populusdeltoides[J]. Genome Resesrch, 2020, 30(8): 1131-1143. doi: 10.1101/gr.261438.120 [14] Soler M, Plasencia A, Larbat R, et al. The Eucalyptus linker histone variant EgH1.3 cooperates with the transcription factor EgMYB1 to control lignin biosynthesis during wood formation[J]. New Phytologist, 2017, 213(1): 287-299. doi: 10.1111/nph.14129 [15] Ployet R, Soler M, Carocha V, et al. Long cold exposure induces transcriptional and biochemical remodelling of xylem secondary cell wall in Eucalyptus[J]. Tree Physiology, 2017, 38(3): 1-14. [16] Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences[J]. Bioinformatics, 2006, 22(13): 1658-1659. doi: 10.1093/bioinformatics/btl158 [17] Altschul S F, Madden T L, Schäffer A A, et al. Gapped BLAST and PSI BLAST: a new generation of protein database search programs[J]. Nucleic Acids Research, 1997, 25(17): 3389-3402. doi: 10.1093/nar/25.17.3389 [18] Zheng Y, Jiao C, Sun H, et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases[J]. Molecular Plant, 2016, 9(12): 1667-1670. doi: 10.1016/j.molp.2016.09.014 [19] Li B, Dewey C N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome[J]. BMC Bioinformatics, 2011, 12(1): 323. doi: 10.1186/1471-2105-12-323 [20] Robinson M D, McCarthy D J, Smyth G K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data[J]. Bioinformatics, 2010, 26(1): 139-140. doi: 10.1093/bioinformatics/btp616 [21] Becker R A, Chambers J M, Wilks A R. The New S Language[M].Monterey: Wadsworth and Brooks/Cole, 1988. [22] Franceschini A, Szklarczyk D, Frankild S, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration[J]. Nucleic Acids Research, 2013, 41: D808. [23] Zhong R, Mccarthy R L, Lee C, et al. Dissection of the transcriptional program regulating secondary wall biosynthesis during wood formation in poplar[J]. Plant Physiology, 2011, 157(3): 1452-1468. doi: 10.1104/pp.111.181354 [24] Ning K, Ding C, Huang Q, et al. Transcriptome profiling revealed diverse gene expression patterns in poplar (Populus×euramericana) under different planting densities[J]. PLoS One, 2019, 14(5): e0217066. doi: 10.1371/journal.pone.0217066 [25] Etchells J P, Mishra L S, Kumar M, et al. Wood formation in trees is increased by manipulating PXY-regulated cell division[J]. Currenrt Biology, 2015, 25(8): 1050-1055. doi: 10.1016/j.cub.2015.02.023 [26] Zhong R, Richardson E A, Ye Z. The MYB46 transcription factor is a direct target of SND1 and regulates secondary wall biosynthesis in Arabidopsis[J]. The Plant Cell, 2007, 19(9): 2776-2792. doi: 10.1105/tpc.107.053678 [27] Ko J H, Kim W C, Han K H, et al. Ectopic expression of MYB46 identifies transcriptional regulatory genes involved in secondary wall biosynthesis in Arabidopsis[J]. Plant Journal, 2009, 60(4): 649-665. doi: 10.1111/j.1365-313X.2009.03989.x [28] Kim W C, Ko J H, Han K H. Identification of a cis-acting regulatory motif recognized by MYB46, a master transcriptional regulator of secondary wall biosynthesis[J]. Plant Molecular Biology, 2012, 78(4-5): 489-501. doi: 10.1007/s11103-012-9880-7 [29] Kim W C, Ko J H, Kim J Y, et al. MYB46 directly regulates the gene expression of secondary wall-associated cellulose synthases in Arabidopsis[J]. Plant Journal, 2013, 73(10): 26-36. [30] Zhong R, Mccarthy R L, Haghighat M, et al. The poplar MYB master switches bind to the SMRE site and activate the secondary wall biosynthetic program during wood formation[J]. PloS One, 2013, 8(7): e69219. doi: 10.1371/journal.pone.0069219 [31] Chai G, Qi G, Cao Y, et al. Poplar PdC3H17 and PdC3H18 are direct targets of PdMYB3 and PdMYB21, and positively regulate secondary wall formation in Arabidopsis and poplar[J]. New Phytologist, 2014, 203(2): 520-534. doi: 10.1111/nph.12825 [32] Furuta K M, Yadav S R, Lehesranta S, et al. Arabidopsis NAC45/86 direct sieve element morphogenesis culminating in enucleation[J]. Science, 2014, 345(6199): 933-937. doi: 10.1126/science.1253736