Temporal and spatial variations of soil bulk density by experimental forest fire in Daxing'an Mountains of northeastern China
-
摘要:目的研究不同强度火干扰下寒温带针叶林土壤密度的时空变化格局及其形成机制, 旨在更好地认识火干扰后森林土壤密度的时空动态规律, 为火烧迹地改造和生态恢复提供参考。方法在大兴安岭北部寒温带杜香-杜鹃-落叶松林地内, 进行严密监控下的秋季林火点烧试验。采用网格法布设固定样地, 在火烧前后连续定点观测土壤密度的时空动态。结果林火干扰后:(1)各时段重度火烧区土壤密度较轻、中度火烧区高5%~10%, 且差异达极显著水平。(2)时间序列上土壤密度的变化模式因火烧强度而异, 其中轻度火烧区土壤密度因火烧立即下降而后逐渐升高, 中度火烧区土壤密度未因火烧立即产生实质性变化而后持续升高, 重度火烧区(及3个火烧强度区综合分析)表现为火烧后土壤密度立即升高并持续增强。(3)土壤密度(或土壤密度变化)空间格局与火烧强度空间格局极显著正相关, 火烧强度的作用效果在火烧后立即显现, 其相关程度在融雪季、雨季/生长季进一步加深。结论兴安落叶松林火烧迹地恢复初期, 因火烧强度差异导致的土壤密度变化方向的初始差异在雨季/生长季消失(均表现为升高), 中、重度火烧区土壤密度显著升高。Abstract:ObjectiveThe temporal and spatial variations and its formation mechanism on soil bulk density were analyzed in a cold temperate coniferous forest under different intensity fire to understand the dynamic temporal and spatial law of soil bulk density after fire interference and to provide reference for post-fire forest improvement and ecological restoration.MethodAn autumn experimental fire was conducted in the boreal Ledum palustre-Rhododendron dauricum-Larix gmelinii forests in Daxing'an Mountains of northeastern China. Before and after the fire, the temporal and spatial variations of soil bulk density were examined continuously by latticed co-coupled sampling method.ResultThe results showed that: (1) At each time node, soil bulk density of severe burning area was 5%-10% higher than that of mild or moderate burning area after the fire. (2) In time series, soil bulk density variation pattern was depended on fire intensity, i.e., soil bulk density of mild burning area decreased immediately and then increased, that of moderate burning area was not changed immediately and then increased continually, that of severe burning area or of all the monitoring site was risen continually. (3) Significant positive correlation was between spatial pattern of soil bulk density (or soil bulk density change) and fire intensity, correlation appeared immediately after the fire, and correlation coefficient was ascended after snowmelt season and rainy season/growing season.ConclusionAt the initial stages of post-fire recovery, the original change direction of soil bulk density in Larix gmelini forest was different because of fire intensity differences and then all increased after rainy season/growing season, which increased significantly in moderate burning area and severe burning area.
-
蔗糖代谢在植物生长发育过程中发挥着重要的调控作用。作为植物中重要的碳源存储形式之一，蔗糖的合成和代谢过程倍受研究者们的关注。其中蔗糖合成过程主要由蔗糖磷酸合酶（sucrose phosphate synthase，SPS）和蔗糖磷酸化酶（sucrose phosphorylase，SPP）进行调控，而蔗糖合酶（sucrose synthase，SUS）调控蔗糖与UDP-葡萄糖和果糖间的可逆合成和分解过程，蔗糖转化酶（invertase，INV）则不可逆地调控蔗糖分解形成葡萄糖和果糖[1-2]。蔗糖转化酶作为蔗糖分解过程中的关键酶，按照其在植物细胞内分布部位的不同，可分为液泡转化酶（vacuolar invertase，VIN）、细胞壁转化酶（cell wall invertase，CWIN）和胞质转化酶（cytoplasmic invertase，CINV） 3类；按照其pH值的不同，细胞壁转化酶和液泡转化酶又被称为酸性转化酶（acid invertase，AIN），而胞质转化酶则被称为碱性/中性转化酶（alkaline/neutral invertase，NIN）[3]。因二者进化来源不一，碱性/中性转化酶与酸性转化酶的同源性较低[4]，其中碱性/中性转化酶类属于糖基水解酶家族100，主要定位于叶绿体、线粒体和胞质溶胶内，与植物的生长发育和胁迫响应相关[5-7]。
前人对于拟南芥（Arabidopsis thaliana）的研究表明：碱性/中性转化酶在植物的生长发育过程中发挥着重要的调控作用[1]。在拟南芥中，碱性/中性转化酶基因家族共包含9个成员，其中AtCINV1（cytoplasmic invertase 1）调控拟南芥根形态的正常发育，抑制其表达会抑制主根生长和促进侧根的发育[8-9]；AtINVE（alkaline/neutral invertase E）则参与拟南芥叶绿体的发育，抑制其表达不仅造成叶绿体发育畸形且提高了植株的氮同化速率[10]；AtINVH（alkaline/neutral invertase H）则在拟南芥再生组织部位中具有特异地高表达，其突变体表现出茎缩短、开花时间延迟等特性，且其根部中的活性氧（ROS）水平急剧下降[11]。此外，在百脉根（Lotus japonicus）中，LjINV1（alkaline/neutral invertase 1）基因的突变会导致植株根和茎长度的缩短，且无花粉产生，影响植株的生殖发育[12]。在巴西橡胶树（Hevea brasiliensis）中，转化酶活性则与乳胶的产量紧密关联[13]。水稻（Oryza sativa）中OsCYT-INV1（cytoplasmic invertase 1）基因的突变则会导致根长的缩短和开花时间的延长以及不育[14]。在杂交山杨（Populus tremula × tremuloides）中，转化酶的活性则与纤维素合成相关，其活性的下降会导致纤维素含量的下降，影响植物的次生生长[15]。此外，碱性/中性转化酶的活性还与植物的胁迫防御相关。在拟南芥中AtINVA（alkaline/neutral invertase A）和AtINVG（alkaline/oneutral invertase G）与氧化防御相关基因的表达紧密相关，过表达AtINVA和AtINVG会抑制抗坏血酸过氧化物酶基因AtAPX2（oxidative stress-responsive ascorbate peroxidase 2）的表达[16]。在小麦（Triticum aestivum）中，渗透胁迫和低温处理则会引起叶片中转化酶活性的增加[17]。在枸橘（Poncirus trifoliata）中，过表达PtrA/NINV（alkaline/neutral invertase）可降低植株体内的活性氧水平并增强其光合作用效率，以应对冷、高盐和干旱等不利环境条件的影响[18]。在小麦中，Ta-A/N-Inv1（Alkaline/neutral invertase 1）的表达则与条锈病相关，其可能通过提高细胞内己糖水平从而降低光合作用，减少活性氧的产生，从而阻止叶片坏死[19]。
杨树（Populus spp.）因其木材性质优良、速生以及基因组结构简单，在木本植物中作为模式植物而具有重要的研究价值。随着毛果杨（Populus trichocarpa）基因组数据的公布，以杨树为代表植物开展对木本植物新陈代谢的研究显得尤为方便[20]。基于基因组数据，前人针对杨树的转化酶基因家族也已开展了相关分析，并对其成员在不同组织部位的表达模式进行了研究[5, 21]，然而对于杨树中转化酶成员基因功能的研究却相对较少。本研究中，我们以毛白杨（Populus tomentosa）碱性/中性转化酶基因PtoNIN1为例，对其基因结构进行相关分析，并对其在不同组织部位的表达模式进行了探究；同时以模式植物拟南芥为例，对其进行过表达验证，初步探究PtoNIN1在植物生长发育过程中的作用，为研究杨树其他转化酶成员功能提供了借鉴，也为揭示杨树蔗糖代谢过程奠定了基础。
1. 材料和方法
1.1 植物材料
本研究中毛白杨根、茎、叶材料来源于实验室内部保存生长6个月的毛白杨优良无性系组培苗TC1521，同时参照Chen等[22-23]方法分别采集成年毛白杨植株8个不同发育时期下雌花芽和雄花芽材料，以及成年树的成熟叶。采用改良的CTAB法分别提取各样品的总RNA，通过Reverse Transcription System Ⅰ（Promega）系统反转录合成cDNA，用于PtoNIN1基因的克隆和组织特异性表达研究。
转基因受体材料为实验室保存的哥伦比亚生态型拟南芥。拟南芥种子经由5%次氯酸钠溶液（AR，Macklin）消毒15 min后移入1/2 MS培养基（添加20 g/L蔗糖和5.5 g/L琼脂），4 ℃冰箱低温储藏3 d后移至气候箱内培养。培养箱温度为（22 ± 1） ℃，采用16 h光照/8 h黑暗的长光照周期进行培养。培养7 d后移至含有蛭石的营养土内进行培养。
1.2 毛白杨同源PtoNIN1基因的克隆与测序
以毛果杨转化酶基因家族的同源基因PtrNIN1（Potri.008G101500）为索引，基于已公布的毛白杨基因组序列[24]设计并合成了克隆引物（表1）。以毛白杨组培苗叶片材料的cDNA为模版，进行PCR扩增。PCR反应体系为：2 × Phanta® Max Master Mix （Vazyme）10 μL，前后引物（1 μmol/L）各0.5 μL，模板cDNA（25 ~ 50 ng/μL）2 μL，ddH2O 0.8 μL，总体系为20 μL。PCR的反应条件为：预变性98 ℃，3 min；变性98 ℃，30 s；退火58 ℃，60 s；延伸72 ℃，1 min，共进行35个循环；再延伸72 ℃，5 min，最后12 ℃保存。PCR克隆后通过FastPure® Gel DNA Extraction Mini Kit （Vazyme）进行切胶回收。连接转化后送至公司进行测序，测序数据由北京睿博兴科生物科技有限公司提供。
表 1 本实验所用引物序列Table 1. Primers used in the study基因
Gene正向引物
Forward primer反向引物
Reverse primerC-PtoNIN1 ATGAACAGTAGTAGTTGTA TCAGACATGAATTTGAGATC O-PtoNIN1 AGAACACGGGGGACTCTTGACATGAACAGT
AGTAGTTGTATTGGTAGGGGAAATTCGAGCTGGTCACTCAGACATG
AATTTGAGATCTTGCPtoNIN1 AACAATGAAACCCTGTTGCAG CCGTATCTTGTGACCAACTACGG PtoActin GGATGGCTATATAGGTTGCAGGAA GAGATGCTGGGAAGGAATTGAAAG AtSUS2 AGGGTGTACCAAATCTCAT CATAGTGAAAGCTGTGTGG AtSUS3 GAGATACCGCAGGGAGAGTT CAGCATTTCAGTCTCAAGGG AtSUS4 CAAGTGTAAGCATGACCC AGGAACAGCTTGAGCCAG AtSPP2 AGAAGCTAGCAACTTCCCTGAG GCTAACCTTGTGTGGCCTCT AtSPS2 TGTTTCAACCTTGGCACAAA GTGGTTTCAACCCTCCTGAA AtEF1A4α TTGACCAGATCAACGAGCCCAAGA ACTCGTGGTGCATCTCAACAGACT 注：C-PtoNIN1为PtoNIN1基因CDS区克隆所用引物，O-PtoNIN1为PtoNIN1过表达载体构建所用引物，PtoNIN1为PtoNIN1基因表达量检测所用引物。Notes: C-PtoNIN1 is used for the gene colning of PtoNIN1, O-PtoNIN1 is used for the construction of over-expression vector, and PtoNIN1 is used for qRT-PCR. 1.3 毛白杨PtoNIN1基因的生物信息学分析
基于Phytozome数据库（https://phytozome-next.jgi.doe.gov/）的毛果杨基因组数据绘制PtoNIN1基因的染色体定位和基因结构，蛋白的理化性质采用在线软件ProtParam tool（https://web.expasy.org/ protparam/）进行分析，保守结构域分析采用在线软件InterProScan（http://www.ebi.ac.uk/interpro/search）和Pfam（http://pfam.xfam.org/ search）进行分析，二级结构分析则通过在线软件SOPMA（https://npsa-prabi.ibcp.fr/ cgi-bin/npsa_automat.pl?page=npsa_sopma.html）进行统计分析，蛋白三级结构采用软件SWISS-MODEL（https://swissmodel.expasy.org/）进行分析，跨膜结构域预测采用在线软件TMHMM 2.0（https://services.healthtech.dtu.dk/service.php?TM HMM-2.0）进行分析，信号肽的预测采用在线软件SignalP-5.0（https://services. healthtech.dtu.dk/service.php?SignalP-5.0）进行分析，亚细胞定位预测采用在线软件PSORT Prediction（http://psort1.hgc.jp/form.html）进行，分别利用在线软件NetPhos 3.1（https://services.healthtech.dtu.dk/service.php?NetPhos-3.1）和NetNGlyc 1.0（https://services.healthte ch.dtu.dk/service.php?NetNGlyc-1.0）对PtoNIN1蛋白的磷酸化位点和糖基化位点进行预测分析，多序列同源比对则采用软件DNAMAN（V9.0.1.116）进行。启动子区序列分析则采用在线软件PlantCARE（http://bioinformatics.psb.ugent. be/-webtools/plantcare/html/）进行分析。
为进一步分析PtoNIN1蛋白的进化关系，基于基因的氨基酸序列通过MEGA11.0构建系统进化树，采用Neighbor-Joining（NJ）算法的Poisson模型，Complete deletion模式进行建树，bootstrap值取1 000。其中，毛果杨包含12个成员，PtrNIN1（Potri.008G101500）、PtrNIN2（Potri.013G006600）、PtrNIN3（Potri.008G024100）、PtrNIN4（Potri.010G236100）、PtrNIN5（Potri.005G010800）、PtrNIN6（Potri.004G186500）、PtrNIN7（Potri.005G239400）、PtrNIN8（Potri.019G082000）、PtrNIN9（Potri.004G167500）、PtrNIN10（Potri.002G173600）、PtrNIN11（Potri.009G129000）、PtrNIN12（Potri.013G110800）；拟南芥包含9个成员，AtCINV1（AT1G35580）、AtCINV2（AT4G09510）、AtINVA（AT1G56560）、AtINVB（AT4G34860）、AtINVC（AT3G06500）、AtINVD（AT1G22650）、AtINVE（AT5G22510）、AtINVF（AT1G72000）、AtINVH（AT3G05820）；水稻包含8个成员，OsNIN1（LOC_Os03g20020）、OsNIN2（LOC_Os01g22900）、OsNIN3（LOC_Os02g32730）、OsNIN4（LOC_Os04g33490）、OsNIN5（LOC_Os02g03320）、OsNIN6（LOC_Os11g07440）、OsNIN7（LOC_Os04g35280）、OsNIN8（LOC_Os02g34560）；其余蛋白成员依次为LbNI（KR026955）、DlNI（KP769773）、DoNI2（KY794404）、HbNIN1（GU573728）、HbNIN2（GU573727）、HbNIN3（KC577600）、HfCIN3（MT379656）、DcNI4（XM_020835477.2）、Ta-A/N-Inv1（AM295169）。其中毛果杨和水稻基因数据来自Phytozome数据库，拟南芥基因数据来自TAIR数据库（https://www.arabidopsis.org/），其他基因数据来自NCBI数据库（https://www.ncbi.nlm.nih.gov/）。同时利用在线软件MEME（https://meme-suite.org/meme/tools/meme）对不同物种间转化酶成员的保守序列元件进行预测分析。
1.4 毛白杨PtoNIN1基因的组织特异性表达分析
实时荧光定量PCR（qRT-PCR）引物（表1）使用Primer3在线软件（http://fokker.wi.mit.edu）进行设计，采用StepOne PlusTM Real-Time PCR System with Tower（ABI，USA）检测系统和ChamQ Universal SYBR qPCR Master Mix（Vazyme）进行qRT-PCR扩增检测，以PtoActin基因为内参对照[22]，采用2−∆∆Ct法计算基因的相对表达量。毛果杨碱性/中性转化酶基因家族成员表达量数据来源于Phytozome数据库中已公布的转录组数据。
1.5 农杆菌介导的拟南芥遗传转化
以实验室所有的pCAMBIA1301载体质粒为骨架，采用无缝克隆试剂盒ClonExpress®ⅡOne Step Cloning Kit（c112-02，Vazyme）构建由35S启动子驱动的PtoNIN1过表达载体，无缝克隆引物如表1所示。载体经由测序验证后转化农杆菌感受态GV3101，通过农杆菌介导的花序浸染法获得拟南芥转基因植株。拟南芥转基因植株采用潮霉素进行筛选，潮霉素使用质量浓度为20 mg/L。
1.6 拟南芥转基因植株的表型鉴定
经由PCR验证后，对萌发后培养30 d的T3代转基因拟南芥植株进行表型测定，同时参照王静澄[25]的方法测定转基因拟南芥根、茎及叶中的果糖和蔗糖含量。采用葡萄糖测定试剂盒（A154-1-1，南京建成生物工程研究所）测定转基因植株中葡萄糖的含量，每个转基因株系分别选取5个单株作为生物学重复，以拟南芥野生型植株为对照，计算转基因植株中蔗糖、葡萄糖和果糖的相对含量。采用qRT-PCR的方法测定拟南芥内源蔗糖代谢相关基因的表达量变化，引物序列如表1所示，以拟南芥AtEF1A4α基因为内参对照[26]，采用2−∆∆Ct法计算基因的相对表达量。
2. 结果与分析
2.1 毛白杨PtoNIN1基因的克隆与生物信息学分析
以毛果杨同源基因PtoNIN1和已公布的毛白杨基因组序列为参考，克隆并获得了毛白杨同源基因序列PtoNIN1。测序结果显示，毛白杨PtoNIN1基因的CDS区序列长度为2 073 bp，共编码690个氨基酸（图1a）；染色体定位结果显示，PtoNIN1位于毛白杨第8号染色体，共包含6个外显子区和5个内含子区（图1b，c）。氨基酸理化性质分析结果显示：PtoNIN1蛋白分子量77.80 kDa，理论等电点为6.22，分子式为C3 474H5 376N948O1 027S32，不稳定系数38.58，为稳定性蛋白。亲疏水性分析结果显示，PtoNIN1蛋白的亲水性总平均值为−0.277，为亲水性蛋白。
![]()
图 1 毛白杨中PtoNIN1的序列（a）、染色体定位（b）和基因结构（c）分析Figure 1. Sequence (a), chromosomal location (b) and gene structure (c) of PtoNIN1 in Populus tomentosa蛋白的二级结构分析结果显示：无规则卷曲在PtoNIN1蛋白中所占比例最高，为39.42%，其次为α螺旋（36.81%）、伸展链（18.26%）和β折叠（5.51%）（图2a）。三级结构分析结果显示：其可形成同型六聚体结构，与鱼腥藻（Anabaena spp.）中碱性转化酶InvA蛋白结构最为相近（图2b）。结构域分析结果显示：其含有糖苷酶超家族100（glycosyl-hydrolase-100 superfamily）的特征结构域，属于碱性和中性转化酶基因家族。亚细胞定位预测结果显示，其可能位于线粒体基质空间中。跨膜结构域和信号肽的预测结果显示，PtoNIN1蛋白不具备跨膜结构域和信号肽。磷酸化位点预测结果显示：PtoNIN1蛋白中存在着丝氨酸、苏氨酸和酪氨酸磷酸化位点，其中丝氨酸磷酸化位点占比最高（图2c）。糖基化位点预测结果显示，PtoNIN1蛋白的2号、37号、44号、63号、123号、147号、176号和189号位置上的天冬酰胺为潜在的糖基化氨基酸位点（图2d）。
![]()
图 2 毛白杨PtoNIN1蛋白的二级结构（a）、三级结构（b）、磷酸化（c）和糖基化（d）位点分析图a中紫色代表无规则卷曲，蓝色代表α螺旋，红色代表伸展链，绿色代表β折叠。In Fig. a, purple represents irregular curls, blue represents α-spiral, red represents extension chain, and green represents β-fold.Figure 2. Secondary (a) and tertiary (b) structure, phosphorylation (c) and glycosylation site (d) of PtoNIN1 in P. tomentosa2.2 毛白杨PtoNIN1基因的多序列比对和系统进化分析
以Phytozome v13.0数据库中毛果杨基因组数据为参照，对毛白杨碱性/中性转化酶PtoNIN1与毛果杨的碱性/中性转化酶基因家族成员进行多序列比对分析（图3）。其中PtoNIN1与毛果杨中性转化酶基因家族中12个成员的氨基酸相似性依次为95.38%、69.13%、53.98%、54.48%、66.05%、55.07%、44.06%、42.46%、43.19%、40.84%、44.06%和42.32%。
![]()
图 3 PtoNIN1蛋白与毛果杨碱性/中性转化酶基因家族成员的多序列比对Figure 3. Multiple sequence alignment of PtoNIN1 with the alkaline/neutral invertase gene family in P. trichocarpa此外，为了进一步探究PtoNIN1的蛋白功能，我们联合了与其同源性较高的其他物种中研究较为透彻的碱性/中性转化酶蛋白构建了系统进化树（图4）。聚类分析结果显示：植物中性转化酶基因家族可分为2大亚类，一类主要定位于胞质内，包含AtCINV1、AtCINV2、Ta-A/N-Inv1、LbNI、HbNIN1/2等成员；另一类主要定位于细胞器内，按照其定位的细胞器类型又可分为2种类型，一类位于线粒体型（Ⅰ-Ⅰ），一类位于叶绿体型（Ⅰ-Ⅱ）。其中，PtoNIN1与毛果杨中同源蛋白PtrNIN1亲缘关系最近，与拟南芥中AtINVA/H/C、水稻中OsNIN1/2和铁皮石斛（Dendrobium officinale）中DoNI2处于同一分支，表明其类属于线粒体型转化酶成员，这也与亚细胞定位预测的结果相一致。MEME分析结果显示：植物碱性/中性转化酶家族中共含有11个不同的保守基序，除毛果杨中PtoNIN7、PtoNIN9、PtoNIN10和PtoNIN11、水稻中OsNIN5和OsNIN7以及AtINVB和LbNI外，其余碱性/中性转化酶成员中均含有这11个保守的基序。
![]()
图 4 植物碱性/中性转化酶基因家族的系统进化和保守基序分析Figure 4. Phylogenetic and conserved motif analysis of alkaline/neutral invertase gene families in plant基于已公布的毛白杨基因组序列，我们以PtoNIN1基因上游3 kbp区域序列为例，对其启动子区顺式作用元件进行预测分析（图5）。其中，光响应元件（13个）是PtoNIN1中分布最为丰富的顺式作用元件，表明PtoNIN1是潜在的光调控基因。此外，干旱诱导（2个）、厌氧诱导（4个）、防御和胁迫响应（2个）以及赤霉素响应元件（2个）也分布于PtoNIN1的启动子区域。
![]()
图 5 PtoNIN1启动子区顺式作用元件分析Figure 5. Analysis of cis-acting elements in the promoter region of PtoNIN12.3 毛白杨PtoNIN1基因的时空特异性表达模式分析
为了进一步探究PtoNIN1在毛白杨发育过程中的作用，我们以PtoActin基因为内参基因，对不同组织部位（根、茎、叶及成年树的成熟叶）和不同发育时期毛白杨雄花芽和雌花芽的表达模式进行分析。结果显示：PtoNIN1在毛白杨根中表达量较高，其次是成年树的成熟叶，在茎和叶中表达量则相对较低（图6a）。对不同发育时期下雌花芽和雄花芽的研究结果表明：PtoNIN1在雌雄花芽成花诱导期（时期S1）具有高表达，随后表达量出现明显下降，并在整个花芽发育时期均保持较为稳定的表达趋势（图6b，c）。
![]()
图 6 PtoNIN1的基因表达模式分析S1. 成花诱导期；S2. 花原基形成期；S3. 器官发生期；S4. 伸长期；S5. 孢子形成期；S6 ~ S7. 休眠期；S8. 小孢子发生期。S1, floral induction stage; S2, primordia formation stage; S3, organogenesis stage; S4, enlargement stage; S5, archespore formation stage; S6−S7, dormancy stages; S8, microsporogenesis stage.Figure 6. Gene expression pattern analysis of PtoNIN12.4 拟南芥中PtoNIN1基因的过表达分析
为进一步探究PtoNIN1的基因功能，我们以pCAMBIA1301载体质粒为骨架，为PtoNIN1构建过表达载体（图7a），并对模式植物拟南芥进行过表达研究。T3代植物经由PCR扩增验证后共获得6个转基因过表达株系（图7b）。
![]()
图 7 拟南芥过表达载体的构建（a）和转基因株系验证（b）图a中基因黑色区域部分为用于荧光定量检测和转基因株系验证所用片段；图b中1 ~ 6分别代表了PtoNIN1的6个表达株系，col则为拟南芥野生型，M为5000 bp的DNA Marker。In Fig. a, the black region of the gene is the fragment used for qRT-PCR and transgenic line validation; in Fig. b, the number of 1−6 represent the six overexpression lines of PtoNIN1, respectively, col represents the wild-type of arabidopsis, and M is a DNA Marker for 5000 bp.Figure 7. PtoNIN1-overexpression vector (a) and validation of transgenic lines (b) in arabidopsis拟南芥中过表达PtoNIN1显著提高了植株中生物量的积累。与野生型相比，转基因株系中仅Pro35S:PtoNIN1/col-2的株高出现显著提升，其余株系与野生型相比并无明显差异（图8a，b）。然而与野生型相比，转基因株系中Pro35S:PtoNIN1/col-2和Pro35S:PtoNIN1/col-4的茎及角果鲜质量显著高于野生型，Pro35S:PtoNIN1/col-2、Pro35S:PtoNIN1/col-3和Pro35S:PtoNIN1/col-4的莲座叶鲜质量显著高于野生型，表明PtoNIN1在拟南芥中的过表达可能影响了植株生物量的积累（图8b，c）。对表型差异明显的转基因株系Pro35S:PtoNIN1/col-2/3/4中蔗糖、葡萄糖和果糖的测定分析结果表明：PtoNIN1的过表达并未影响莲座叶中蔗糖、葡萄糖和果糖的含量；然而与野生型相比，转基因株系茎及角果中蔗糖的含量显著提升，而Pro35S:PtoNIN1/col-4中果糖含量相比于野生型出现显著下降，其余株系中果糖和葡萄糖的含量相对于野生型并未产生显著差异（图8d ~ f）。
![]()
图 8 拟南芥及其转基因株系的表型分析图a中左侧为培养40 d的拟南芥野生型植株，右侧为培养40 d的Pro35S:PtoNIN1/col-2转基因过表达植株；图b ~ g中横坐标col表示野生型，1 ~ 6分别表示转基因株系Pro35S:PtoNIN1/col-1/2/3/4/5/6。其中图e ~ g中黑色柱子代表莲座叶部位的糖含量比较，灰色柱子代表茎及角果部位的糖含量比较；**代表在P < 0.01下差异显著，*代表在P < 0.05下差异显著。下同。In Fig. a, the left side shows the wild-type plants in arabidopsis cultured for 40 d, and the right side shows the transgenic plants of Pro35S: PtoNIN1/col-2 cultured for 40 d; in Fig. b−g, col represents the wild-type plants in arabidopsis, and the number of 1−6 represent the transgenic plants of Pro35S: PtoNIN1/col-1/2/3/4/5/6, respectively; thereinto, the black column in Fig. e−g represents the sugar content in the rosette leaves, and the gray column represents the sugar content in the stems and siliques. ** means significant difference at P < 0.01 level, * means significant difference at P < 0.05 level. The same below.Figure 8. Phenotypic analysis of arabidopsis and its transgenic lines此外，以3个表型最为突出的拟南芥转基因株系为例，我们也对外源基因PtoNIN1及拟南芥内源蔗糖代谢相关基因的表达量变化进行了分析。其中外源基因PtoNIN1的表达量与拟南芥表型具有显著相关关系，表型最为显著的Pro35S:PtoNIN1/col-2株系中外源PtoNIN1的表达量相对较高（图9a）。参照前人研究[26-28]，我们分别选取了野生型拟南芥中表达量相对较高的蔗糖合酶成员（AtSUS2-4）、蔗糖磷酸合酶成员（AtSPS2）和蔗糖磷酸化酶成员（AtSPP2）对拟南芥转基因株系中蔗糖代谢过程的变化进行分析。相比于野生型，拟南芥叶片、茎及角果中的AtSUS2-4、AtSPS2和AtSPP2的表达量均出现了显著上升，且相比于叶片，茎及角果中基因表达量的上升更为显著（图9b ~ f）。
![]()
图 9 拟南芥及其转基因株系的基因表达模式分析黑色柱子代表莲座叶部位中基因表达模式分析，灰色柱子代表茎及角果部位的基因表达模式分析，2 ~ 4分别表示转基因株系Pro35S:PtoNIN1/col-2/3/4。The black column represents the gene expression patterns in the rosette leaves, and the gray column represents the gene expression patterns in the stems and siliques, and the number of 2−4 represent the transgenic plants of Pro35S: PtoNIN1/col-2/3/4, respectively.Figure 9. Gene expression patterns of arabidopsis and its transgenic lines3. 讨 论
蔗糖代谢作为植物生长发育过程中的重要代谢过程，可通过分解蔗糖为植物提供生长发育所需的能源，另一方面蔗糖分解所形成的单糖分子又可作为信号分子调控植物的生长发育过程[1]。转化酶作为植物蔗糖分解过程中2类关键酶之一，在植物的生长发育、信号转导、碳源分配和外界环境刺激响应等方面发挥了重要作用[29]。然而，前人的研究多集中于酸性转化酶家族，碱性/中性转化酶因其稳定性差、活性低而在前期被研究者们认为其仅起到酸性转化酶功能的补充作用[30]。随着测序数据的公布和基因组学的不断发展，越来越多研究者意识到碱性/中性转化酶家族具有其独特的功能，且在水稻（Oryza sativa）[31]、玉米（Zea mays）[32]、棉花（Gossypium barbadense）[33]、柑橘（Citrus reticulata）[34]、葡萄（Vitis vinifera）[35]、梨（Pyrus bretschneideri）[36]和苹果（Malus pumila）[37]等多个植物物种中开展了基因家族的遗传进化分析。然而对于中性转化酶基因功能的研究却相对较少，且多集中于模式植物拟南芥[4, 21, 25]和一些草本植物，如水稻[14]、小麦[19]、枸杞[38]等，木本植物中的研究仍十分匮乏。
杨树作为木本植物中的模式植物，随着其基因组数据的公布，研究者们也针对其中性转化酶家族开展了相关研究[5]，然而对于杨树碱性/中性转化酶基因功能的研究却鲜有报道。本研究中我们以毛果杨中PtrNIN1序列为索引，对已公布的毛白杨基因组序列进行同源比对分析，从中获得了毛白杨同源基因序列PtoNIN1。对其进行生物信息学分析显示：其具有碱性/中性转化酶家族特征结构域，且包含12个保守的氨基酸基序，这与前人关于杨树碱性/中性转化酶基因家族的研究相一致[5]。对其进行三级结构预测分析显示：其与鱼腥藻中InvA结构极为类似，该基因是植物碱性/中性转化酶基因家族的潜在祖先之一，可特异识别蔗糖的α-1, 2-糖苷键[39]。同源序列比对和进化树分析均表明PtoNIN1与PtrNIN1具有高度相似性，表明碱性/中性转化酶在杨属不同物种间具有保守性。此外，PtoNIN1在进化上与AtINVA/H/C、OsNIN1和DoNI2处于同一分支，该类蛋白均定位于线粒体内，这也与亚细胞定位预测的结果相一致。在拟南芥中，AtINVA/C在植物的根、茎和花均具有明显的表达[40]，OsNIN1在水稻的各个组织部位均具有显著的表达[31]，DoNI2在铁皮石斛中根、茎、叶也被检测出具有较高的表达量[41]，表明线粒体定位的转化酶成员在植株各个组织部位的发育中都具有重要的功能。本研究中PtoNIN1在毛白杨的根、茎、叶及成熟叶同样被检测出具有较高的表达量，表明了其在毛白杨生长发育过程中的重要性，然而不同于草莓（Fragaria × ananassa）[42]和水稻[43]中的研究，PtoNIN1在根中的表达量最高，高于成熟叶。此外，AtINVA/C/H在拟南芥中被证实参与了植物的花发育过程，其突变体均延迟了植物的开花时间[11, 40]。本研究中我们对PtoNIN1在毛白杨雌雄花序发育过程中的作用进行了探究，PtoNIN1在雌雄花序的发育过程中具有相似的表达趋势，其在成花诱导期具有显著的高表达，随后保持了一个相对稳定的表达趋势，我们推测PtoNIN1可能对于雌雄花芽的诱导和发育过程的维持具有重要作用。启动子区分析结果表明PtoNIN1受光照、厌氧环境、干旱、防御和胁迫等环境信号调控，此外赤霉素响应元件也存在于PtoNIN1的启动子区，与木薯 （Manihot esculenta）中MeNINV1[44]相类似。
此外，为进一步探究PtoNIN1的生物学功能，我们以模式植物拟南芥为例对其进行了过表达研究。相比于野生型，转基因植株的莲座叶、茎及角果的鲜质量均得到了显著提升。前人对拟南芥中PtoNIN1同源基因AtINVA/C/H的研究已表明，该类型基因的突变会造成了植株茎的缩短[11, 40]，然而在拟南芥中过表达PtoNIN1却并未显著增强其植株高度（除Pro35S:PtoNIN1/col-2外），我们推测这可能是由于PtoNIN1和拟南芥中的同源基因发生了功能冗余。转化酶是蔗糖分解过程中的关键酶，我们同样对野生型和转基因植株间的可溶性糖含量进行了分析比较，其中PtoNIN1的过表达显著增强了茎及角果中蔗糖含量的积累，而葡萄糖和果糖的含量却并未发生明显变化，我们推测PtoNIN1可能增强了植物整体的蔗糖代谢过程，蔗糖代谢相关基因表达量验证的结果也证实了这一结论，蔗糖合酶、蔗糖磷酸化酶和蔗糖磷酸磷酸化酶相关基因的表达量在转基因株系中都得到了显著提升。
本研究通过同源克隆的方式从毛白杨基因组数据中获取碱性/中性转化酶成员PtoNIN1，剖析了PtoNIN1的基因和蛋白序列结构特征，揭示了PtoNIN1在不同器官组织和不同雌雄花芽发育时期的表达模式，同时利用模式植物拟南芥异源过表达初步探究了PtoNIN1的基因功能，为进一步深入研究杨树碱性/中性转化酶基因的遗传调控机制，揭示杨树蔗糖代谢过程奠定了前期基础。
-
![]()
图 1 样地监测点位布设及林火强度格局
数字2~9表示火烧强度等级，其中，2、3表示轻度火烧，4、5、6表示中度火烧，7、8、9表示重度火烧。
Figure 1. Monitoring sites (small black dots) and the pattern of FI of sample plot
Number from 2 to 9 means fire intensity (FI) subclass, mild subclass includes 2 and 3; moderate subclass includes 4, 5 and 6; severe subclass includes 7, 8 and 9.
![]()
图 2 火烧时间序列土壤密度
T0为火烧前1d，T1为火烧后1d，T2为融雪后，T3为雨季/生长季后。下同。总体即3个火烧强度区综合分析。
Figure 2. Soil bulk density (SBD) in time series
T0 means the day before fire, T1 means the first day after fire, T2 means after snowmelt season, and T3 means after rainy and growing season. The same below. The overall means all the fire intensity area.
表 1 火烧后土壤密度点位变化率分级统计
Table 1 Monitoring site classification statistics on post-fire soil bulk density change (SBDC)
火烧强度
Fire intensity时间节点
Time node火烧后土壤密度相对于火烧前变化率的分级统计
Monitoring site classification statistics on post-fire soil bulk density change/%平均变化率及成对数据t检验
Average change rate and the result of paired t test/%-5%~0 0~5% 5%~10% 10%~15% 15%~20% 20%~25% 25%~30% 轻 T1 75.00 25.00 — — — — — -1.04** Mild T2 35.00 65.00 — — — — — -0.28ns n=20 T3 25.00 75.00 — — — — — 0.30ns 中 T1 74.07 25.93 — — — — — 0.24ns Moderate T2 3.70 88.89 3.70 3.70 — — — 2.61** n=27 T3 — 55.56 37.04 3.70 3.70 — — 5.41** 重 T1 — — 42.11 47.37 10.53 — — 10.93** Severe T2 — — 42.11 47.37 10.53 — — 11.13** n=19 T3 — — 5.26 31.58 42.11 10.53 10.53 16.87** 注：每个点位分别计算土壤密度变化率，即SBDCTn(%)=(SBDTn-SBDT0)/SBDT0×100，数据为各变化区间的点位分布频率%，即点位分布频率=(某变化区间的点位数/该火烧区点位数n)×100；点位平均变化率为Tn时段该火烧区各点位变化率的平均值；经成对数据t检验分析，**表示相对于T0时段土壤密度显著变化(P < 0.01)，ns表示相对于T0时段无显著差异。Notes: SBDC of each monitoring site was calculated as follows: SBDCT n(%)=(SBDTn-SBDT0)/SBD T0×100, data was distribution frequency of monitoring site in each statistics classification, and was calculated as follows: distribution frequency (%)=(number of monitoring site in statistics classification/number of monitoring site in fire intensity area)×100; average change rate was average of SBDC at Tn time; ** means SBDC significant difference (P < 0.01) between Tn and T0, and ns means no significant difference by paired-samples t test. 
下载: 导出CSV
表 2 土壤密度及其变化率与火烧强度相关性分析
Table 2 Correlation analysis between SBD or SBDC and FI
时间序列
Time series土壤密度与火烧强度相关性
Correlation analysis between SBD and FI土壤密度变化率与火烧强度相关性
Correlation analysis between SBDC and FI相关方程
Correlation equation相关系数
Correlation coefficient(r)相关方程
Correlation equation相关系数
Correlation coefficient(r)T0 SBD=-0.007FI+0.690 0.468** T1 SBD=0.005FI+0.642 0.361** SBDC=2.006FI-7.437 0.794** T2 SBD=0.005FI+0.652 0.370** SBDC=1.988FI-6.086 0.805** T3 SBD=0.011FI+0.637 0.670** SBDC=3.007FI-8.378 0.877** 注：**表示土壤密度(或土壤密度变化率)与火烧强度之间相关性显著(P < 0.01)。Note: ** means significant correlation (P < 0.01) between SBD (or SBDC) and fire intensity.
下载: 导出CSV
-
[1] Johnstone J F, Chapin Ⅲ F S. Effects of soil burn severity on post-fire tree recruitment in boreal forest[J]. Ecosystems, 2006, 9:14-31. doi: 10.1007/s10021-004-0042-x
[2] 宋雨, 胡海清, 孙龙, 等.大兴安岭不同坡位地表可燃物含水率的动态变化与建模[J].森林工程, 2017, 33(5):1-7. http://d.old.wanfangdata.com.cn/Periodical/slgc201705002 Song Y, Hu H Q, Sun L, et al. Dynamic change and modeling of moisture content of surface fuel in different slope positions of Daxing'anling[J].Forest Engineering, 2017, 33(5):1-7. http://d.old.wanfangdata.com.cn/Periodical/slgc201705002
[3] Greene R S B, Chartres C J, Hodgkinson K C. The effects of fire on the soil in a degraded semi-arid woodland (Ⅰ): cryptogam cover and physical and micromorphological properties[J].Australian Journal of Soil Research, 1990, 28(5):755-777. doi: 10.1071/SR9900755
[4] 刘道锟, 孙海龙, 甘秋妹, 等.大兴安岭干旱阳坡不同植被退化阶段土壤理化性质与物种多样性研究[J].森林工程, 2016, 32(2): 1-6. doi: 10.3969/j.issn.1001-005X.2016.02.001 Liu D K, Sun H L, Gan Q M, et al. Soil physicochemical properties and plant species diversity of different vegetation degradation stages in arid sunny-slope of Great Xing'an Mountains[J].Forest Engineering, 2016, 32(2):1-6. doi: 10.3969/j.issn.1001-005X.2016.02.001
[5] Certini G. Effects of fire on properties of forest soils: a review[J].Oecologia, 2005, 143(1):1-10. doi: 10.1007/s00442-004-1788-8
[6] 孙向阳.土壤学[M].北京:中国林业出版社, 2005:132. Sun X Y. Soil science[M].Beijing: China Forestry Publishing House, 2005:132.
[7] 崔晓阳, 方怀龙.城市绿地土壤及其管理[M].北京:中国林业出版社, 2001:67. Cui X Y, Fang H L. Urban green land and its management[M].Beijing: China Forestry Publishing House, 2001:67.
[8] 张喜, 朱军, 崔迎春, 等.火烧对黔中喀斯特山地马尾松林土壤理化性质的影响[J].生态学报, 2011, 31(19): 5809-5817. http://d.old.wanfangdata.com.cn/Periodical/stxb201119045 Zhang X, Zhu J, Cui Y C, et al. Influence of fire on a Pinus massoniana soil in a karst mountain area at the center of Guizhou Province, China[J]. Acta Ecologica Sinica, 2011, 31(19):5809-5817. http://d.old.wanfangdata.com.cn/Periodical/stxb201119045
[9] Badìa D, Martì C. Plant ash and heat intensity effects on chemical and physical properties of two contrasting soils[J]. Arid Land Research Management, 2003, 17(1):23-41. doi: 10.1080/15324980301595
[10] 张玉红, 孙铭隆, 刘彤.林火对大兴安岭典型植被土壤理化性质的影响[J].东北林业大学学报, 2012, 40(6):41-43. doi: 10.3969/j.issn.1000-5382.2012.06.011 Zhang Y H, Sun M L, Liu T. Effect of forest fire on soil physical and chemical properties of typical forests in Daxing'an Mountains[J]. Journal of Northeast Forestry University, 2012, 40(6):41-43. doi: 10.3969/j.issn.1000-5382.2012.06.011
[11] 孙龙, 赵俊, 胡海清.中度火干扰对白桦落叶松混交林土壤理化性质的影响[J].林业科学, 2011, 47(2):103-110. doi: 10.3969/j.issn.1006-1126.2011.02.006 Sun L, Zhao J, Hu H Q. Effect of moderate fire disturbance on soil physical and chemical properties of Betula platyphylla-Larix gmelinii mixed forest[J].Scientia Silvae Sinicae, 2011, 47(2):103-110. doi: 10.3969/j.issn.1006-1126.2011.02.006
[12] 胡海清, 刘洋, 孙龙, 等.火烧对不同林型下森林土壤水分物理性质的影响[J].水土保持学报, 2008, 22(2):162-165. doi: 10.3321/j.issn:1009-2242.2008.02.037 Hu H Q, Liu Y, Sun L, et al. Effect of fire on soil hydro-physical properties under different types of forest land[J]. Journal of Soil and Water Conservation, 2008, 22(2):162-165. doi: 10.3321/j.issn:1009-2242.2008.02.037
[13] 耿玉清, 周荣伍, 李涛, 等.北京西山地区林火对土壤性质的影响[J].中国水土保持科学, 2007, 5(5):66-70. doi: 10.3969/j.issn.1672-3007.2007.05.012 Geng Y Q, Zhou R W, Li T, et al. Influences of forest fire on soil properties in Xishan Area of Beijing[J].Science of Soil and Water Conservation, 2007, 5(5):66-70. doi: 10.3969/j.issn.1672-3007.2007.05.012
[14] 许鹏波, 屈明, 薛立.火对森林土壤的影响[J].生态学杂志, 2013, 32(6):1596-1606. http://d.old.wanfangdata.com.cn/Periodical/stxb201707005 Xu P B, Qu M, Xue L. Effects of forest fire on forest soils[J].Chinese Journal of Ecology, 2013, 32(6):1596-1606. http://d.old.wanfangdata.com.cn/Periodical/stxb201707005
[15] 王艳霞, 杨桂英, 张政.火干扰对森林土壤理化性质的影响[J].防护林科技, 2012(4):1-3. doi: 10.3969/j.issn.1005-5215.2012.04.002 Wang Y X, Yang G Y, Zhang Z. Effect of fire disturbance on soil physical & chemical properties of forest[J].Protection Forest Science and Technology, 2012(4):1-3. doi: 10.3969/j.issn.1005-5215.2012.04.002
[16] 金慧, 赵莹, 赵伟, 等.林火对长白山自然保护区风灾区土壤性质的影响[J].北华大学学报(自然科学版), 2013, 14(6): 698-703. http://d.old.wanfangdata.com.cn/Periodical/bhdxxb201306020 Jin H, Zhao Y, Zhao W, et al. Influences of forest fire on soil properties in wind disaster area of Changbai Mountain Nature Reserve[J]. Journal of Beihua University(Natural Science), 2013, 14(6):698-703. http://d.old.wanfangdata.com.cn/Periodical/bhdxxb201306020
[17] Rab M A. Soil physical and hydrological properties following logging and slash burning in the Eucalyptus regnans forest of southeastern Australia[J].Forest Ecology and Management, 1996, 84:159-176. doi: 10.1016/0378-1127(96)03740-1
[18] 李海洋, 胡海清, 孙龙.我国森林死可燃物含水率与气象和土壤因子关系模型研究[J].森林工程, 2016, 32(3):1-6. doi: 10.3969/j.issn.1001-005X.2016.03.001 Li H Y, Hu H Q, Sun L. Research on relational models of moisture content of dead forest fuel with meteorological factors and soil factors in China[J].Forest Engineering, 2016, 32(3):1-6. doi: 10.3969/j.issn.1001-005X.2016.03.001
[19] 杨光, 舒立福, 邸雪颖.气候变化影响下大兴安岭地区21世纪森林火险等级变化预测[J].应用生态学报, 2012, 23(12): 3236-3242. http://d.old.wanfangdata.com.cn/Periodical/yystxb201212002 Yang G, Shu L F, Di X Y. Prediction on the changes of forest fire danger rating in Great Xing'an Mountain Region of Northeast China in the 21st century under effects of climate change[J]. Chinese Journal of Applied Ecology, 2012, 23(12):3236-3242. http://d.old.wanfangdata.com.cn/Periodical/yystxb201212002
[20] Cui X Y, Gao F, Song J F, et al. Changes in soil total organic carbon after an experimental fire in a cold temperate coniferous forest: a sequenced monitoring approach[J]. Geoderma, 2014, 226-227:260-269. doi: 10.1016/j.geoderma.2014.02.010
[21] 张韫.土壤·植物·水理化分析教程[M].北京:中国林业出版社, 2011:22-23. Zhang Y. Soil, plants, water physical and chemical analysis tutorial[M].Beijing: China Forestry Publishing House, 2011:22-23.
[22] 张敏.林火对土壤环境影响的研究[D].哈尔滨: 东北林业大学, 2015: 25-26. Zhang M. The effect of forest fire on soil environment[D].Harbin: Northeast Forestry University, 2015: 25-26.
[23] 温美丽, 刘宝元, 魏欣, 等.冻融作用对东北黑土容重的影响[J].土壤通报, 2009, 40(3):492-495. http://d.old.wanfangdata.com.cn/Periodical/trtb200903006 Wen M L, Liu B Y, Wei X, et al. Impact of freezing and thawing on bulk density of black soil[J]. Chinese Journal of Soil Science, 2009, 40(3):492-495. http://d.old.wanfangdata.com.cn/Periodical/trtb200903006
[24] 王恩姮, 赵雨森, 夏祥友, 等.冻融交替后不同尺度黑土结构变化特征[J].生态学报, 2014, 34(21):6287-6296. http://d.old.wanfangdata.com.cn/Periodical/stxb201421028 Wang E H, Zhao Y S, Xia X Y, et al. Effects of freeze-thaw cycles on black soil structure at different size scales[J]. Acta Ecologica Sinica, 2014, 34(21):6287-6296. http://d.old.wanfangdata.com.cn/Periodical/stxb201421028
-
期刊类型引用(4)
1. 杨桦，李祥乾，王帆，方睿，杨伟. 长足大竹象信息素结合蛋白CbuqPBP2互作蛋白的筛选与验证. 西北农林科技大学学报(自然科学版). 2024(01): 87-97 . 
百度学术
2. 万超，张月，胡莉，伍炳华，袁媛. 茉莉花JsMYB305转录因子的原核表达及蛋白纯化. 福建农业学报. 2022(02): 164-169 . 百度学术
3. 武建颖，张燕，孙贺贺，赵玉兰，董金皋，申珅，郝志敏. 玉米大斑病菌蛋白激酶A催化亚基StPKA-C1/C2的表达与互作蛋白筛选. 农业生物技术学报. 2022(10): 1976-1986 . 百度学术
4. 胡莉，万超，张蕖，陈清西，伍炳华，袁媛. 茉莉花JsMYB305转录因子互作蛋白的筛选及验证. 福建农业学报. 2022(09): 1135-1144 . 百度学术
其他类型引用(9)
计量
- 文章访问数: 1233
- HTML全文浏览量: 306
- PDF下载量: 32
- 被引次数: 13
