Estimation of PM2.5 emissions from forest fires in Fujian Sanming City of southern China during 2000-2011
-
摘要:目的福建三明地区属于森林火灾多发地区, 森林火险等级较高, 森林火灾在燃烧过程中释放大量可吸入颗粒物(PM2.5), 不仅对大气环境造成污染, 还会对人类的身体健康造成不良影响, 因此, 研究福建三明地区森林火灾PM2.5排放量具有重要意义。方法本文基于福建省三明市2000—2011年森林火灾数据, 采用蓄积量-生物量模型计算林火过程消耗生物量, 其中, 森林损失蓄积量和过火面积数据由当地防火部门提供。运用排放因子法计算得出PM2.5排放量, 各林型PM2.5排放因子由历年文献统计得出。针对福建三明地区4种主要林型:马尾松林、杉木林、阔叶混交林和针阔混交林对森林火灾中PM2.5排放量分别进行估算, 为该地区森林火灾污染物排放研究提供依据。结果结果表明, 福建三明地区2000—2011年间森林火灾PM2.5估算总排放量为458.789~2958.339t。其中杉木林燃烧PM2.5排放量估算为217.117~1072.489t, 阔叶混交林, 马尾松林, 针阔混交林PM2.5估算排放量分别为94.123~903.154t、93.379~478.634t、54.169~504.063t。结论综合分析, 按不同森林火灾等级分析, 由于重大火灾和较大火灾发生次数较多, PM2.5排放主要集中在重大火灾和较大火灾; 按不同森林类型分析, PM2.5排放主要集中在杉木林和阔叶混交林两种林型, 杉木林和阔叶混交林由于其地上生物量相对丰富, 且易发生森林火灾, 是今后重点抚育管理对象, 相关部门应在防火期到来之前做好修枝抚育间伐工作。Abstract:ObjectiveForest fire frequently occurs in Fujian Sanming area of southern China, and has high forest fire danger class. Forest fire releases large amount of particulate matter (PM2.5), which not only causes air pollution, but also has bad effect on human health. Therefore, the research on area forest fire PM2.5 emissions in Fujian Sanming City is of great significance.MethodBased on forest fire data in Fujian Sanming area during 2000-2011, we used volume-biomass model to calculate fire process consume biomass, volume loss and forest fire area data provided by the local forestry authorities. Using emission factor method can calculate PM2.5 emissions, and the forest PM2.5 emission factor was obtained from historical literature statistics. PM2.5 emissions of four main forest types:Pinus massoniana forests, Cunninghamia lanceolata forests, mixed broadleaved forests and coniferous broad leaved mixed forests were assessed individually in Fujian Sanming area. These results will provide theory basis for forest fire prevention work in this area.ResultThe results showed that forest fire PM2.5 emission estimates were 458.789-2958.339t in Fujian Sanming area during 2000-2011. PM2.5 emission of Cunninghamia lanceolate forests was 217.117-1072.489t, PM2.5 emission of temperate mixed broadleaved forests, Pinus massoniana forest, coniferous and broadleaved mixed forests was 94.123-903.154t, 93.379-478.634t, 54.169-504.063t, respectively.ConclusionFrom comprehensive analysis, as a result of more times of disastrous fire and large fire, PM2.5 emission mainly concentrated in the disastrous fire and large fire; about forest types, PM2.5 emissions are mainly concentrated in Cunninghamia lanceolata forests and mixed broadleaved forests. Because of its high biomass aboveground and prone to forest fire, these two forest types will be the key management object in future work, and the relevant departments should take good thinning of pruning work before the arrival of fireproofing period.
-
Keywords:
- forest fire /
- PM2.5 /
- emission factor method /
- fuel load /
- forest type
-
蔗糖代谢在植物生长发育过程中发挥着重要的调控作用。作为植物中重要的碳源存储形式之一，蔗糖的合成和代谢过程倍受研究者们的关注。其中蔗糖合成过程主要由蔗糖磷酸合酶（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 不同森林类型不同火强度下损失的生物量
Figure 1. Biomass loss of different forest types under varied fire intensities
![]()
图 2 不同森林类型不同火烧等级下PM2.5排放量
Figure 2. PM2.5 emission under different fire size classes for varied forest types
表 1 森林类型蓄积量-生物量模型参数
Table 1 Model parameters between volume and biomass of forest type
森林类型Forest type BEF (Biomass expansion factor) 杉木林Cunninghamia lanceolata forest BEF=2.4381 Age-0.3293(R2=0.746) 马尾松林Pinus massoniana forest BEF=0.9267 Age-0.1317(R2=0.447) 针阔混交林Coniferous and broadleaved mixed forest a=0.8019, b=12.2799 阔叶混交林Mixed broadleaved forest BEF=2.5887 Age-0.1821(R2=0.746) 
下载: 导出CSV
表 2 不同林型各器官生物量分配比
Table 2 Biomass allocation of various organs for different forest types
% 森林类型Forest type 生物量分配Biomass allocation 干Trunk 枝Branch 叶Leaf 马尾松林Pinus massoniana forest 73.2 16.3 4.5 杉木Cunninghamia lanceolata forest 67.0 11.0 9.2 针阔混交林Coniferous and broadleaved mixed forest 60.7~70.0 14.4~18.4 6.3~10.1 阔叶混交林Mixed broadleaved forest 65.4~72.9 12.3~19.8 7.5~11.4
下载: 导出CSV
表 3 乔木在不同等级森林火灾中各器官燃烧效率
Table 3 Combustion efficiency of each tree organs under different fire size classes
% 树木器官Tree organ 一般火灾General fire 较大火灾Larger fire 重大火灾Disastrous fire 特大火灾Serious fire 树枝Branch 14.1~25.0 29.8~65.0 37.3~85.0 53.7~100 树叶Leaf 22.5~32.5 64.0~89.5 89.5~95.5 95.5~100 树干Trunk 21.6~25.0 39.8~65.0 54.8~85.0 61.2~100
下载: 导出CSV
-
[1] 陈媛, 岑况, NORRA S, 等.北京市区大气气溶胶PM2.5污染特征及颗粒物溯源与追踪分析[J].现代地质, 2010, 24(2): 345 -354. doi: 10.3969/j.issn.1000-8527.2010.02.020 Chen Y, Cen K, NORRA S, et al. Study on pollution characteristics of PM2.5 in the aerosol, tracing and tracking atmospheric particulates in Beijing City[J]. Geoscience, 2010, 24(2):345-354. doi: 10.3969/j.issn.1000-8527.2010.02.020
[2] Fuzzi S, Baltensperger U, Carslaw K, et al. Particulate matter, air quality and climate: lessons learned and future needs[J]. Atmospheric Chemistry and Physics, 2015, 15(1):8217-8299. http://cn.bing.com/academic/profile?id=b0b75b06267a5b5af234541e6c4ab44d&encoded=0&v=paper_preview&mkt=zh-cn
[3] Don S, Ricardo C, Samuel T, et al. Using National Ambient Air Quality Standards for fine particulate matter to assess regional wildland fire smoke and air quality management[J]. Journal of Environmental Management, 2017, 201:345-356. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=2da9c8b8ed8d7f48b16aa5f7604f265c
[4] Dennekamp M, Straney L D, Erbas B, et al.Forest fire smoke exposures and out-of-hospital cardiac arrests in Melbourne, Australia:a case-crossover study[J].Environmental Health Perspectives, 2015, 123(10):959-964. doi: 10.1289/ehp.1408436
[5] Neal F, Breanna A, Richard A B, et al. The health impacts and economic value of wildland fire episodes in the US: 2008-2012[J]. Science of the Total Environment, 2018, 610-611:802-809. doi: 10.1016/j.scitotenv.2017.08.024
[6] Song C B, He J J, Wu L, et al. Health burden attributable to ambient PM2.5 in China[J]. Environmental Pollution, 2017, 223:575-586. doi: 10.1016/j.envpol.2017.01.060
[7] Tara S, Narasimhan L, Miriam R, et al.PM2.5 measurements in wildfire smoke plumes from fire seasons 2005-2008 in the Northwestern United States[J].Journal of Aerosol Science, 2011, 42(3):143-155. doi: 10.1016/j.jaerosci.2010.09.001
[8] Saunders R O, Waugh D W. Variability and potential sources of summer PM2.5 in the Northeastern United States[J]. Atmospheric Environment, 2015, 117:259-270. doi: 10.1016/j.atmosenv.2015.07.007
[9] Cesari D, De Benedetto G E, Dinoi A, et al. Seasonal variability of PM2.5 and PM10 composition and sources in an urban background site in Southern Italy[J]. Science of the Total Environment, 2018, 612:202-213. https://www.sciencedirect.com/science/article/pii/S0048969717322283
[10] Bondur V G. Satellite monitoring of trace gas and aerosol emissions during wildfires in Russia[J]. Lzvestiya Atmospheric & Oceanic Physics, 2016, 52(9):1078-1091. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c25a977f774c53b7d8f1a163f8318c7e
[11] Zhang L B, Liu Y Q, Hao L. Contributions of open crop straw burning emissions to PM2.5 concentrations in China[J].Environmental Research Letters, 2016, 11(1):1-9. http://cn.bing.com/academic/profile?id=386e629ccc026ddc45cfcb7e787d6fad&encoded=0&v=paper_preview&mkt=zh-cn
[12] 祝斌, 朱先磊, 张元勋, 等.农作物秸秆燃烧PM2.5排放因子的研究[J].环境科学研究, 2005, 18(2):29-34. doi: 10.3321/j.issn:1001-6929.2005.02.006 Zhu B, Zhu X L, Zhang Y X, et al. Emission factor of PM2.5 from crop straw burning[J]. Research of Environmental Sciences, 2005, 18(2):29-34. doi: 10.3321/j.issn:1001-6929.2005.02.006
[13] 唐喜斌, 黄成, 楼晟荣, 等.长三角地区秸秆燃烧排放因子与颗粒物成分谱研究[J].环境科学, 2014, 35(5):1623-1632. http://d.old.wanfangdata.com.cn/Periodical/hjkx201405001 Tang X B, Huang C, Lou S R, et al. Emission factors and PM chemical composition study of biomass in the Yangtze River Delta Region[J].Environmental Science, 2014, 35(5):1623-1632. http://d.old.wanfangdata.com.cn/Periodical/hjkx201405001
[14] 刘晓东, 王博.森林燃烧主要排放物研究进展[J].北京林业大学学报, 2017, 39(12):118-124. doi: 10.13332/j.1000--1522.20170307 Liu X D, Wang B. Review on the main emission products released by forest combustion[J]. Journal of Beijing Forestry University, 2017, 39(12):118-124. doi: 10.13332/j.1000--1522.20170307
[15] 靳全锋, 马祥庆, 王文辉, 等.中国亚热带地区2000—2014年林火排放颗粒物时空动态变化[J].环境科学学报, 2017, 37(6):2238-2247. http://www.cnki.com.cn/Article/CJFDTOTAL-HJXX201706028.htm Jin Q F, Ma X Q, Wang W H, et al. Temporal and spatial characteristics of particulate matter emission from forest fires in subtropical China during 2000-2014[J]. Acta Scientiae Circumstantiae, 2017, 37(6):2238-2247. http://www.cnki.com.cn/Article/CJFDTOTAL-HJXX201706028.htm
[16] 靳全锋, 王文辉, 马祥庆, 等.福建省2000—2010年林火排放污染物时空动态变化[J].中国环境科学, 2017, 37(2):476-485. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGHJ201702010.htm Jin Q F, Wang W H, Ma X Q, et al. Temporal and spatial dynamics of pollutants emission from forest fires in Fujian during 2000-2010[J]. China Environmental Science, 2017, 37(2):476-485. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGHJ201702010.htm
[17] 彭徐剑, 郭福涛, 胡海清, 等.大兴安岭森林火灾NO和SO2释放量的估算[J].林业科技开发, 2015, 29(5):134-138. http://d.old.wanfangdata.com.cn/Periodical/lykjkf201505033 Peng X J, Guo F T, Hu H Q, et al. Estimation of NO and SO2 emission by forest fires in the Da Hinggan Moutains[J]. China Forestry Science and Technology, 2015, 29(5):134-138. http://d.old.wanfangdata.com.cn/Periodical/lykjkf201505033
[18] 胡海清, 魏书精, 孙龙.大兴安岭2001—2010年森林火灾碳排放的计量估算[J].生态学报, 2012, 32(17):5373-5386. http://www.cnki.com.cn/Article/CJFDTOTAL-STXB201217011.htm Hu H Q, Wei S J, Sun L. Estimating carbon emissions from forest fires during 2001 to 2010 in Daxinganling Mountain[J]. Acta Ecologica Sinica, 2012, 32(17):5373-5386. http://www.cnki.com.cn/Article/CJFDTOTAL-STXB201217011.htm
[19] 鲁斯唯.海峡西岸经济区大气污染物排放清单研究[D].厦门: 厦门大学, 2014. Lu S W. Research on air pollutant emission inventory in the western Taiwan straits economic zone[D]. Xiamen: Xiamen University, 2014.
[20] 郭怀文.福建三明地区森林火灾碳排放研究[D].北京: 北京林业大学, 2013. Guo H W. Study on forest fires carbon emissions in Sanming area of Fujian Province[D]. Beijing: Beijing Forestry University, 2013.
[21] Seiler W, Crutzen P J. Estimates of gross and net fluxes of carbon between the biosphere and the atmosphere from biomass burning[J]. Climate Change, 1980, 2(3):207-247. doi: 10.1007/BF00137988
[22] 方精云, 陈安平, 赵淑清, 等.中国森林生物量的估算:对Fang等Science一文(Science, 2001, 291:2320-2322)的若干说明[J].植物生态学报, 2002, 26(2): 243-249. doi: 10.3321/j.issn:1005-264X.2002.02.018 Fang J Y, Chen A P, Zhao S Q, et al. Estimating biomass of China's forests: supplementary notes on report published in Science (2001, 291:2320-2322) by Fang et al[J]. Chinese Journal of Plant Ecology, 2002, 26(2): 243-249. doi: 10.3321/j.issn:1005-264X.2002.02.018
[23] 吕新双.黑龙江大兴安岭主要乔木树种火灾碳释放研究[D].哈尔滨: 东北林业大学, 2006. Lü X S.Carbon emissions from forest fire on main tree species in Daxing'an Mountain in Heilongjiang Province[D]. Harbin: Northeast Forestry University, 2006.
[24] Urbanski S P. Combustion fficiency and emission factors for wildfire-season fires in mixed conifer forests of the northern rocky mountains, US[J]. Infection and Immunity, 2008, 76(11): 5341-5349. doi: 10.1128/IAI.00786-08
[25] Alves C, Gon Alves C, Fernandes A P. Fireplace and woodstove fine particle emissions from combustion of western Mediterranean wood types[J]. Atmospheric Research, 2011, 101(3):692-700. doi: 10.1016/j.atmosres.2011.04.015
[26] Fernandes A P, Alves C A, Tarelho L, et al. Emission factors from residential combustion appliances burning Portuguese biomass fuels[J]. Journal of Environmental Monitoring, 2011, 13(11):3196-3206. doi: 10.1039/c1em10500k
[27] 张坤洪, 梁远楠, 杨俊伟, 等.西江中下游针阔混交林优势种群的结构及分布格局研究[J].广东林业科技, 2009, 25(4):1-6. doi: 10.3969/j.issn.1006-4427.2009.04.001 Zhang K H, Liang Y N, Yang J W, et al. Structure and spatial paterns of dominant populations in the conifer and broadleaf mixed forest in middle-lower reaches of Xijiang River[J]. Guangdong Forestry Science and Technology, 2009, 25(4):1-6. doi: 10.3969/j.issn.1006-4427.2009.04.001
[28] 赵凤君, 舒立福, 姚树人.森林火灾碳排放估算方法与研究进展[J].森林防火, 2012(1):25-29. doi: 10.3969/j.issn.1002-2511.2012.01.015 Zhao F J, Shu L F, Yao S R. Forest fire estimation method and research progress of carbon emissions[J]. Forest Fire Prevention, 2012(1):25-29. doi: 10.3969/j.issn.1002-2511.2012.01.015
[29] 陈戈萍, 张思玉, 王红干.马尾松凋落物燃烧烟气成分比较分析[J].消防科学与技术, 2014, 33(9):999-1002. doi: 10.3969/j.issn.1009-0029.2014.09.006 Chen G P, Zhang S Y, Wang H G. The compare of the components from combust emission of Pinus massoniana litterfall[J]. Fire Science and Technology, 2014, 33(9):999-1002. doi: 10.3969/j.issn.1009-0029.2014.09.006
-
期刊类型引用(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)
计量
- 文章访问数: 1428
- HTML全文浏览量: 352
- PDF下载量: 35
- 被引次数: 13
