不同氮形态处理条件下杨树根尖差异表达基因的特征分析

周 婧，李卓蓉，吴江婷

(林木遗传育种国家重点实验室，国家林业和草原局森林培育重点实验室，中国林业科学研究院林业研究所，北京 100091)

氮素是农林植物生长发育必不可少的大量元素之一，也是农林植物根系从土壤中吸收最多的矿质元素，对其生长发育具有重要意义[1-2]。但在土壤中无机氮含量往往较低[3]，从而导致农林植物生产力低下，限制了其经济效益的发挥。大量含硝态氮或铵态氮无机氮肥的施用，成为维持农林土壤肥力的主要手段之一。但盲目、过量施肥往往造成土壤养分失衡、氮素流失等诸多环境问题[4]。因此，深入挖掘农林植物响应氮素的关键基因，解析其分子调控机制，提高农林植物对氮素吸收同化能力，既可以实现农林生态系统氮素高效利用的经济效益，也能实现减少土壤氮肥施加，保护生态环境的环保效益。

杨树（PopulusL.）是我国主要的速生丰产用材林树种之一[5]。由于杨树生长迅速、对养分消耗大，对氮素有着很强的需求[3]。已有研究表明，不同氮形态处理，能够导致小黑杨（P.simonii×P.nigra）根系形态发生改变，从而影响其对氮素的吸收同化能力[6]。作者前期研究表明，硝态氮处理，能够改变灰杨（Populus × canescens）根尖形态，并影响根尖不同区段对硝态氮的吸收速率[7]。上述研究均表明，杨树根尖能够通过响应不同氮形态，改变其自身形态特征，从而影响根尖对不同氮形态的吸收同化能力[6-8]。

近年来，随着高通量测序技术的快速发展，转录组测序成为研究基因差异表达的重要手段。利用高通量测序技术对杨树根尖响应氮素的研究也取得了显著成果，但是研究大多集中在氮素亏缺方面[9-10]。例如，通过转录组重测序以及生物信息学分析表明，低氮处理导致灰杨根系中特异表达的PtaNAC1（NAC-domain protein）上调表达，且增加了杨树根系生物量，并显著改变了PtaNAC1下游基因的表达，从而影响杨树根系形态结构[9]。Dash 等的研究也表明，低氮处理条件下，在灰杨根系中，PtaHWS（Hawaiian Skirt），PtaNAC1和PtaRAP 2.11（subfamilies of ERF/AP2 TF family）显著差异表达，并影响杨树根系形态建成[10]。上述研究均通过转录组学和生物信息学方法，对杨树根系响应低氮胁迫的过程进行了研究。而不同氮形态处理条件下，杨树根尖存在怎样的差异表达模式尚不清楚，值得深入研究。

本研究以灰杨根尖为试材，利用Illumina 测序平台，筛选不同氮形态处理条件下，杨树根尖差异表达基因，并结合生物信息学分析，构建杨树根尖对不同氮形态响应过程的互作调控网络，挖掘关键基因。最终，阐明杨树根尖响应硝态氮或铵态氮，从而影响根尖生长发育过程的分子调控机制。研究成果为后续开发高氮素吸收利用效率的杨树新种质提供了科学依据。

1 材料和方法

1.1 试验材料

将生根的灰杨组培苗（4 周）转移至水培，隔天换水浇灌改良的LA（Long Ashton）营养液（0.5 mmol·L-1NH4NO3、0.5 mmol·L-1KCl、0.9 mmol·L-1CaCl2、0.3 mmol·L-1MgSO4、0.6 mmol·L-1KH2PO4、42 μmol·L-1K2HPO4、10 μmol·L-1Fe-EDTA、2 μmol·L-1MnSO4、10 μmol·L-1H3BO3、7 μmol·L-1Na2MoO4、0.05 μmol·L-1CoSO4、0.2 μmol·L-1ZnSO4和 0.2 μmol·L-1CuSO4,pH 5.5）。2 周后，选取具有相似高度的48 株植株， 成2 组。分别添加0.5 mmol·L-1硝态氮（NO3-）和0.5 mmol·L-1铵态氮（NH4+）代替LA 营养液中的NH4NO3，进行水培，培养时间为10 d[7]。对根系形态特征进行分析。同时，收获根尖0～4 cm样品，迅速置于液氮中，放于-80℃冰箱备用。将8 株植物的样本等量混合作为1 个重复，每种处理水平3 个生物学重复。

1.2 试验方法

1.2.1 根尖形态特征分析 试验苗进行不同氮形态处理10 d 后，对杨树根系构型进行观察测量，包括测量主根长度，以及测量可见侧根距根尖的距离。每个处理水平测量8 株植株，且每种处理水平3 个生物学重复。进一步，对数据统计进行分析，利用Statgraphics 软件，先对数据进行正态性检验，采用单因素方差分析（ANOVA），以P＜0.05 作为统计意义上的显著水平。

1.2.2 不同氮形态处理条件下转录组测序mRNAs 测序文库构建和测序按照Illumina 公司提供的标准步骤执行。将上述收获的2 种不同氮形态处理条件下的根尖（0～4 cm）样品，每个处理3个重复，分别利用RNA 提取试剂盒提取总RNAs（TRK1001,联川生物技术公司，杭州，中国）。将提取的总RNAs 经过DNase I 消化后，富集含有polyA 尾的mRNA，反转录成cDNA。合成第二链cDNA 后，用QIA quick PCR extraction 试剂盒进行纯化。经过回收纯化、PCR 富集后，获得cDNA文库。将构建好的文库，利用杭州联川生物技术公司IlluminaHiSeqTM4 000 高通量测序技术，进行转录组测序分析。

1.2.3 转录组测序结果分析 将测得的RNAs 序列，使用联川生物技术公司开发的软件进行分析。首先，将得到的原始序列进行去杂质处理，包括接头序列的读取片段、质量值Q ≤ 5 的低质量序列读取片段和未知碱基数比例＞10%的序列读取片段。经过筛选得到的序列读取片段使用SOAPaligner/soap2 将测序结果比对到灰杨基因组数据库（http://aspendb.uga.edu/index.php/databases/spta-717-genome）[11]。采用FPKM（Fragments Per Kilobase of exon per Million mapped）对基因的表达水平进行定量[12]。基于FPKM 值，使用Ballgown package计算mRNAs 的差异表达水平。利用mRNAs 在硝态氮处理条件下的FPKM 除以铵态氮处理条件下FPKM 来计算基因的差异倍数（Fold change，FC）。差异表达mRNAs 筛选阈值为log2（FC）≥1 或≤ -1，且P＜ 0.05。

1.2.4 荧光定量PCR 引物设计及合成 以上述提取的不同氮形态处理杨树根尖样品为模板，使用Takara 公司PrimeScriptTMRT reagent Kit with gDNA Eraser 试剂盒反转录合成cDNA 第一链，作为荧光定量PCR 反应的模板。从转录组数据中选取10 个差异表达的基因，根据灰杨数据库，利用Primer Premier 3.0 软件设计引物，并由上海生工生物工程有限公司合成相关引物。Actin为RTqPCR 的内参基因。序列见表1。使用SYBR Premix Ex TaqTM荧光定量试剂盒进行荧光定量检测。每个样品进行3 次重复。将qPCR 得到的Ct 值进行归一化，计算差异表达基因的相对表达量[7]。

表1 RT-qPCR 引物序列Table 1 Primers used for RT-qPCR

1.2.5 差异基因的功能注释与富集分析 通过Blast2GO 软件，对所有差异基因进行GO 功能注释[13]。随后进行GO 显著性富集分析，以上述得到目的GO 条目（term）为单位，通过P值Bonferroni校正（P＜ 0.05）来定义差异表达基因中显著富集的GO term。同时，结合生物信息学KEGG 软件，分析差异表达基因的生物学功能[14]。同样以correctedP＜ 0.05 为标准，筛选在差异表达基因中显著富集的代谢途径。利用MapMan 软件、Popgenie 数据库（https://popgenie.org/exnet）以及Cytoscape程序，构建杨树根尖响应不同氮形态相关差异表达基因的互作调控网络[15]。

2 结果与分析

2.1 根系形态特征分析

在不同氮形态处理条件下，灰杨根系形态特征具有显著差别（图1）。硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍，这说明不同氮形态处理影响杨树根系生长发育过程。同时，在两种处理水平下，能够发挥吸收营养和水分作用的可见侧根发生的位置相对一致。因此选取侧根形成前的0～4 cm 根尖部分进行后续研究。

图1 不同氮形态处理后，灰杨根系表型Fig.1 Morphological parameters of P. × canescens roots with different nitrogen forms

2.2 转录组测序数据分析

基于上述形态特征分析的结果，以不同氮形态处理杨树根尖为材料，进行高通量转录组测序分析。测序结果显示，两个处理水平的测序分别获得57 001 409 和57 517 359 条有效序列（表2），占各自文库的比例均超过了99.20%（表2）。两个处理水平的测序数据，能够比对到灰杨数据库中的碱基数比例分别为81.59%和82.60%，说明测序所得序列与灰杨数据库比对性较好，可以用于后续分析。

表2 转录组文库测序数据统计Table 2 Distribution of RNA-seq in different categories

2.3 不同氮形态处理条件下，杨树根尖差异表达基因分析

在不同氮形态处理条件下，以铵态氮处理水平为对照，在硝态氮处理条件下，共筛选出2 207 个差异表达基因（p＜ 0.05）。相比较于铵态氮处理，硝态氮处理有1 414 个基因上调表达，793 个基因下调表达。有32 个基因显著差异表达（p＜ 0.001），见图2a。其中，Potri.013G102100基因差异变化最显著，log2（FC）为-6.56，该基因注释为II 型过氧化物酶（peroxiredoxin type 2）；Potri.004G138100基因上调倍数最显著，log2（FC）为6.18，然而该基因尚未被注释，功能未知。同时，筛选出71 个与植物生长发育相关的差异表达基因。其中，转录因子作为一大类调控基因，在杨树根尖通过响应不同氮形态而参与其生长发育过程中发挥着重要作用。例如，NFYA（nuclear transcription factor Y subunit alpha）、GRF（growthregulating factor）和ARF（auxin response factor）类转录因子，它们均能够通过调控硝态氮转运体NRT（nitrate transporters）基因的表达响应不同氮形态，从而调节植物根系生长发育过程[16-17]。进一步，挑选了10 个显著差异表达基因，利用RTqPCR 进一步证实了测序结果的可靠性（图2b）。

图2 不同氮形态处理条件下，杨树根尖显著差异表达的基因Fig.2 Significantly differentially expressed mRNAs in root tips of P. × canescens under different nitrogen form treatments

2.4 不同氮形态处理条件下，灰杨根尖差异表达基因的GO 功能划分

为进一步了解差异表达基因的功能，将差异表达基因映射GO 数据库，根据序列同源性，将不同氮形态处理条件下，2 207 个差异表达基因分为3 大类总共50 个功能组，包括生物学过程（biological process，BP）、细胞组分（cellular component，CC）和分子功能（molecular function，MF）（图3）。结果显示，生物学过程主要集中在转录调控、刺激反应和代谢过程。在细胞组分类别中，差异基因被注释到最多的亚类依次为细胞核、质膜和细胞质部分。在分子功能中，注释为ATP 结合和蛋白结合的差异基因最多，其次为转录因子活性和DNA 结合。以上注释较多的过程所涉及到的差异基因，可能参与杨树根尖对不同氮形态的响应，从而影响根尖生长发育过程。进一步挖掘这些基因的功能，研究其响应不同氮形态的作用机制，将有助于为培育高氮素吸收利用效率的杨树新种质奠定基础。

图3 不同氮形态处理条件下，灰杨根尖差异表达基因的GO 功能聚类注释Fig.3 Gene Ontology classification annotation of differentially expressed genes in root tips of P. × canescens under different nitrogen form treatments

2.5 不同氮形态处理条件下，灰杨根尖差异表达基因的KEGG 分析

为进一步了解上述差异表达基因的功能，对识别出的差异表达基因进行KEGG 通路分析。共获得1 192 个差异表达基因，涉及到20 条KEGG 通路（图4）。其中，有9 条与代谢相关，并且大多与氮代谢和氨基酸代谢相关，有4 条与生物合成相关。KEGG 分析从功能的角度锁定基因，直观的显示了杨树根尖响应不同氮形态差异基因的代谢过程，有助于更好地研究杨树根尖响应不同氮形态的分子机制。

图4 KEGG 通路分析鉴定差异表达基因Fig.4 KEGG pathway analysis of significantly differentially expressed target genes

2.6 与氮代谢相关差异表达基因的MapMan分析

利用MapMan 软件，筛选出36 个与氮代谢通路过程、各类氨基酸的生物合成以及代谢过程相关的差异表达基因（表3）。其中，有8 个基因参与到氮代谢过程，26 个基因参与到氨基酸代谢与合成过程，还有2 个基因参与到硝态氮转运过程。已有研究表明，硝态氮处理能够改变NRTs 的表达水平，进一步影响植物根尖生长发育过程[18-19]。

表3 氮代谢相关差异表达基因Table 3 Differentially expressed genes related to nitrogen metabolism

2.7 不同氮形态处理条件下，灰杨根尖差异表达基因的互作网络分析

为了进一步研究这些基因的功能，利用Popgenie数据库以及Cytoscape 程序，构建了灰杨根尖对不同氮形态响应过程的互作调控网络。在基因作用网络中，检测到一组相互作用的与氮代谢过程相关的基因簇，共包含13 个基因（图5）。其中，硝酸还原酶（Potri.005G172400）基因作为硝态氮同化过程中的关键酶，占据了网络的核心位置。已有研究表明，硝酸还原酶作为植物氮素同化的关键酶之一，能够调控根系生长素水平，从而影响根尖生长发育过程[20-21]。进一步，挑选了调控网络中7 个显著差异表达基因，选取不同氮形态处理10 天的灰杨根尖材料进行RT-qPCR 实验。结果显示：这7 个基因在RT-qPCR 的结果与转录组分析的趋势一致（图6）。

图5 不同氮形态处理条件下，灰杨根尖差异表达基因的互作网络Fig.5 The interaction network of differentially expressed genes in root tips of P. × canescens under different nitrogen forms treatments

图6 不同氮形态处理条件下，灰杨根尖差异基因表达分析Fig.6 Validation of significantly differentially expressed genes under different nitrogen forms treatments by sRNA-seq and RT-qPCR

3 讨论

氮素是植物生长发育所必需的大量矿质营养元素之一，是蛋白质、核酸、叶绿素等重要物质的组成成分[22]。土壤中的氮素都需要在适宜的温度、水分和通气条件下，在土壤微生物和酶的作用下，将其水解或氧化为硝态氮或铵态氮，才能够被植物根系直接吸收利用[23-24]。而土壤中的不同氮形态可能导致杨树根尖氮代谢相关基因表达模式发生改变，从而影响植物根尖生长和发育过程[25]。本研究在不同氮形态处理10 d 后，对杨树根尖进行观察发现，硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍。类似的现象在小黑杨根系中（P.simonii×P.nigra）也有报道[6]。为了进一步了解木本植物根尖对不同氮形态响应的过程，本研究利用转录组测序技术，对不同氮形态处理条件下，杨树根尖差异表达基因进行了分析。

在不同氮形态处理10 d 后，经过高通量测序分析，获得了2 207 个差异表达基因。其中，Potri.013G102100基因编码一个II 型过氧化物酶（peroxiredoxin type 2），其下调倍数最为显著，前人研究表明，II 型过氧化物酶是植物中分布最广的一类非典型2-半胱氨酸过氧化物酶（2-Cysteine Peroxiredoxin），其在对抗氧化能力和共生固氮中发挥着重要作用[26]。另有研究表明，2-半胱氨酸过氧化物酶可以与氮同化过程中亚硝酸还原酶（nitrite reductase）和谷氨酸合成酶（glutamate synthase）相互作用，从而影响氮代谢过程[27]。Potri.004G138100基因上调倍数最显著，但该基因功能未知，推测其可能参与杨树根尖对不同氮形态的响应，从而影响杨树根尖生长发育过程。同时，在本实验中，部分差异表达转录因子通过响应不同氮形态，影响杨树根尖生长发育过程，这与在草本植物中的研究结果相似[6,8,28]。例如，在水稻（Oryza sativaL.）和小麦（Triticum aestivum）中，过表达NFYAs 能够调控NRTs的表达[18,29]。同时，增加了转基因小麦根尖对硝态氮的吸收，促进了侧根生长[18]。在本研究中，相比较于铵态氮处理，硝态氮处理导致灰杨根尖PcNFYA7（Potri.011G 101000）下调表达，同时促进杨树主根伸长。在拟南芥（Arabidopsis thaliana）根中柱鞘细胞中，5mM 硝态氮可以促进ARF6/8的表达水平，促进了拟南芥侧根起始和随后萌发的过程[17,30]。在本研究中，相比较于铵态氮处理，硝态氮处理导致灰杨根尖中PcARF8/16（Potri.004G078200和Potri.010G223200）上调表达，同时促进了主根伸长。上述研究表明，灰杨根尖部分差异表达的基因，能够通过响应不同氮形态而参与调控根尖生长发育过程。

对差异基因进行GO 和KEGG 通路分析。结果表明，一些差异表达基因参与氮代谢和氨基酸代谢途径。例如，氮代谢、色氨酸代谢、丙酮酸代谢和氨基甲酸代谢途径等。同时，利用MapMan筛选出36 个氮代谢相关差异表达基因。其中，值得关注的是两个硝态氮转运体基因PcNRT1.5（Potri.014G179400和Potri.003G088800）。相比较于铵态氮处理，硝态氮处理导致灰杨根尖中PcNRT1.5显著下调表达。在拟南芥中，AtNRT1.1作为第一个被鉴定出的NRT 在根中有转运和感应硝态氮两个作用[31]。同时，在小麦和水稻中，硝态氮处理均能改变NRT 基因的表达水平，从而影响植物根尖生长发育过程[18-19]。该结果说明PcNRT1.5在杨树根尖响应不同氮形态过程中发挥着重要的调控作用。进一步，通过Popgenie 数据库以及Cytoscape程序，构建了灰杨根尖响应不同氮形态的互作调控网络。在基因作用网络中，硝酸还原酶（Potri.005G 172400）作为硝态氮同化过程中第一步的关键酶，在杨树根尖对不同氮形态响应过程中发挥着关键作用。这与在草本植物中的研究结果类似[20]。在拟南芥根中，硝酸还原酶作为氮素同化和根系构型重塑的关键酶之一，能够调控根系生长素水平，从而影响根尖生长发育过程[20]。同时，在水稻中，硝酸还原酶生成的一氧化氮，通过诱导水稻侧根形成和无机氮的吸收，提高水稻对氮的吸收能力[21]。而杨树作为多年生木本植物，其根系比浅根系草本植物发达得多，且构型更为复杂[9]。同时，木本植物往往遭受季节性、物候性土壤氮素反复变化，面临的土壤氮素环境复杂程度远远大于草本植物[4,32-33]。因此，杨树根系中硝酸还原酶可能存在着更为复杂的氮适应性变化的作用机制，值得深入研究。同时，通过基因互作网络的分析，可以逐步了解参与互作网络的基因在响应不同氮形态过程中的作用，其相互作用关系也有待进一步确认。

4 结论

本研究在不同氮形态处理10 d 后，对杨树根尖进行观察发现，硝态氮处理下的主根长度几乎是铵态氮处理条件下的一倍。同时，通过转录组测序分析，筛选出2 207 个差异表达基因。通过GO 功能和KEGG 通路分析，理清了差异表达基因富集的分子功能与代谢通路。利用MapMan筛选出36 个氮代谢相关差异表达基因，并构建了这些基因的互作网络，得到了关键基因。综上所述，本研究为挖掘和寻找杨树根尖响应不同氮形态关键调控基因提供了参考，并阐述了这些关键基因可能通过响应不同氮形态，影响杨树根尖生长发育过程，这为进一步培育出高效吸收利用氮素的杨树新种质奠定了基础。