Snow hydrological characteristics of Larix gmelinii forest in northern Daxing'an Mountains of northeastern China
-
摘要:目的研究大兴安岭北部地区兴安落叶松林雪水文特征，为今后积雪蒸发测定和区域水资源调控提供更为科学的数据支撑和理论基础。方法对观测期内16场降雪的大气降雪量以及对林内积雪深度、积雪密度以及雪水当量进行了周期性观测与统计分析。结果(1) 随着降雪级别的减少，兴安落叶松林的截留率呈现逐渐增加的趋势，分别为6.50%(暴雪)、9.04%(大雪)、9.8%(中雪)、15.7%(小雪)。可见，兴安落叶松林降雪截留最大截留率出现在小雪，暴雪时截留率则最小。(2)兴安落叶松林内积雪深度和林外相比差异不大，其中落叶松林林内积雪深度最深为68.6cm，林外空地积雪深度最深为74.8cm。(3)林内和林外的积雪密度在观测初期会随降雪的输入而降低，无降雪期有相应升高。融雪期随着气温的升高，雪密度减少速度会加快。在4月24—29日达到最大值，减少量分别为0.07和0.11g/cm3。(4)雪水当量速率减少量在4月24—29日达到最大值，分别为30.2和46.4mm。结论和林外空地相比，兴安落叶松林对积雪深度、积雪密度及雪水当量影响不大，说明兴安落叶松林在雪水文过程中对积雪特征影响尽管存在，但并不明显。该区雪蒸发日变化呈单峰曲线变化规律，积雪期的日蒸发量和蒸发速率均值分别为0.04mm和0.2×10-3mm/h，日蒸发量波动幅度在0.02~0.14mm之间，在融雪期间，日蒸发量和蒸发速率的均值分别为0.38mm和1.51×10-3mm/h，采用灰色关联度对各因子进行分析，得出净辐射是影响兴安落叶松林内积雪蒸发的主要因素。Abstract:ObjectiveTo provide more scientific data support and theoretical basis for the future snow cover evaporation measurement and regional water resurces requlation, the snow hydrological characteristics in the natural forest of Larix gmelinii was studied in northern Daxing'an Mountains of northeastern China. The snowfall interception and snow cover characteristics were systematically studied from the end of October 2015 to the beginning of May 2016.MethodDuring the observation period, the snowfall of 16 snowfalls, as well as the snow depth, snow density and snow water equivalent in the forest were periodically observed and statistically analyzed.Result(1) With the decreasing of snowfall level, the interception rate of Larix gmelinii forest showed a gradual increasing trend, which was 6.50% (blizzard), 9.04% (heavy snow), 9.8% (moderate snow), 15.7% (light snow). It can be seen that the maximum interception rate of snowfall in Larix gmelinii forest occurred in light snow, and the smallest occurred in blizzard. (2) The difference between the depth of snow inside Larix gmelinii forest and that outside the forest was not significant. Among them, the deepest snow depth in larch forest was 68.6cm, and that in outer space was 74.8cm. (3) The density of snow inside and outside the forest decreased with the input of snowfall in the early stage of observation, and did not increase in the absence of snowfall. Snowmelt as the temperature rises, the snow density will decrease faster. (4) Reaching the maximum between April 24 and April 29, with a decrease of 0.07 and 0.11g/cm3, respectively; the reduction of snow equivalent rate reached the maximum between April 24 and April 29, 30.2 and 46.4mm, respectively.ConclusionCompared with the open space outside the forest, the Larix gmelinii forest has little effect on the snow depth, snow cover density and snow water equivalent, indicating that the effect of Larix gmelinii forest on the hydrological characteristics of snow cover is not obvious. The daily change of snow evapotranspiration in this area showed a single-peak curve variation. The daily evaporation and evaporation rates in the snow-covered period were 0.04mm and 0.2×10-3mm/h, respectively. Daily evaporation fluctuated between 0.02-0.14mm. During the snowmelt period, the average daily evaporation and evaporation rates were 0.38mm and 1.51×10-3mm/h, respectively, and the factors were analyzed using the grey correlation degree. It is concluded that net radiation is the main factor influencing snowcover evaporation in Larix gmelinii forest.
-
Keywords:
- Larix gmelinii forest /
- snowcover evaporation /
- snowfall interception /
- snowmelting
-
蔗糖代谢在植物生长发育过程中发挥着重要的调控作用。作为植物中重要的碳源存储形式之一，蔗糖的合成和代谢过程倍受研究者们的关注。其中蔗糖合成过程主要由蔗糖磷酸合酶（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的基因功能，为进一步深入研究杨树碱性/中性转化酶基因的遗传调控机制，揭示杨树蔗糖代谢过程奠定了前期基础。
-
![]()
图 3 兴安落叶松林林内及林外空地积雪深度特征
Figure 3. Snow depth characteristics inside and outside Larix gemelinii forest
![]()
图 5 兴安落叶松林雪水当量特征
Figure 5. characteristics of snow water equivalence of larix gmelinii forest
![]()
图 6 积雪融雪期日蒸发量动态变化特征
Figure 6. Dynamic variation characteristics of daily evaporation during snowcover and snowmelting period
![]()
图 7 积雪融雪期气象因子的动态变化特征
Figure 7. Dynamic variation characteristics of meteorological factors during snowcover and snowmelting period
表 1 兴安落叶松天然林林冠截留特征
Table 1 Canopy interception characteristics of Larix gmelinii forest
降雪级别
Snowfall level林外降雪量
Snowfall outsideforest/mm林内降雪量
Snowfall insideforest/mm林冠截留量
Canopy interception/mm林冠截留率
Canopy interception rate/%小雪Light snow 1.94 1.63 0.31 15.7 中雪Moderate snow 2.73 2.46 0.27 9.8 大雪Heavy snow 5.88 5.35 0.53 9.04 暴雪Blizzard 20.63 19.28 1.35 6.5 
下载: 导出CSV
表 2 不同气象因子与积雪蒸发量的灰色关联度
Table 2 Grey correlation degree of different meteorological factors and snowcover evaporation
因子Factor 环境温度
Environment temperature/℃环境湿度
Environment humidity/%水汽压
Water vapor pressure/kPa风速
Wind speed/(m·s-1)净辐射
Net radiation/(MJ·m-2·d-1)积雪期蒸发量Snow period evaporation 0.831 1 0.741 2 0.541 1 0.544 6 0.896 5 融雪期蒸发量Snowmelting period evaporation 0.827 8 0.686 6 0.533 4 0.582 7 0.856 6
下载: 导出CSV
-
[1] 刘世荣, 常建国, 孙鹏森.森林水文学:全球变化背景下的森林与水的关系[J].植物生态学报, 2007, 31(5):753-756. http://d.old.wanfangdata.com.cn/Periodical/zwstxb200705001 Liu S R, Chang J G, Sun P S. Forest hydrology: forest and water in a context of global change[J]. Journal of Plant Ecology, 2007, 31 (5):753-756. http://d.old.wanfangdata.com.cn/Periodical/zwstxb200705001
[2] 盛后财, 蔡体久, 琚存勇.小兴安岭白桦林降水转化过程元素特征分析[J].北京林业大学学报, 2015, 37(2):59-66. doi: 10.13332/j.cnki.jbfu.2015.02.009 Sheng H C, Cai T J, Ju C Y. Element characteristics in the precipitation conversion process in Betula platyphlla forest of Xiaoxing'an Mountains, northeastern China[J]. Journal of Beijing Forestry University, 2015, 37 (2):59-66. doi: 10.13332/j.cnki.jbfu.2015.02.009
[3] 王增艳, 车涛. 2002—2009年中国干旱区积雪时空分布特征[J].干旱区研究, 2012, 29(3): 464-467. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ghqyj201203014 Wang Z Y, Che T. Spatiotemporal distribution of snow cover in arid regions in China[J]. Arid Zone Research, 2012, 29(3): 464-467. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ghqyj201203014
[4] Wei X, Liu S, Zhou G, et al. Hydrological processes in major types of Chinese forest[J]. Hydrological Processes, 2005, 19(1):63-75. http://cn.bing.com/academic/profile?id=ccec2017c9ca1ef970d1dedffceb46a4&encoded=0&v=paper_preview&mkt=zh-cn
[5] Dozlier J. Snow, snowmelt, rain, runoff, and chemistry in a Sierra Nevada Watershed[Z]. Final Report to California Air Resources Board, 1989: 16-26.
[6] Martinet J. Rango A. Paramenter values for snowmelt runoff modeling[J]. Journa1 of Hydrology, 1986, 84(3-4):197-219. http://cn.bing.com/academic/profile?id=6a99879c59f02f9b9a5b642dcab50102&encoded=0&v=paper_preview&mkt=zh-cn
[7] Sokratov S A. Parameters influencing the recrystallization rate of snow[J]. Cold Regions Science and Technology, 2001, 33(2-3): 263-274. doi: 10.1016/S0165-232X(01)00053-2
[8] Pomeroy J W, Brun E. Physical properties of snow[M]//Jones H G, Pomeroy J W, Walker D A, et al. Snow ecology: an interdisciplinary examination of snow-covered ecosystems. Cambridge: Cambridge University Press, 2001: 45-126.
[9] 白重媛, 大畑哲夫.天山乌鲁木齐河源1号冰川夏季消融期内反射率的变化[J].冰川冻土, 1989, 11(4): 311-324. http://www.cnki.com.cn/Article/CJFDTotal-BCDT198904002.htm Bai Z Y, Ohata T. Variations of albedo on the glacier No.1 at the head water of Urumqi River, Tianshan Moutains, during the summer ablation period[J]. Journal of Glaciology and Geocryology, 1989, 11(4): 311-324. http://www.cnki.com.cn/Article/CJFDTotal-BCDT198904002.htm
[10] 谢应钦, 张金生.雪层内太阳的穿透辐射[J].冰川冻土, 1988, 10(2): 135-142. http://www.cnki.com.cn/Article/CJFDTOTAL-BCDT198802004.htm Xie Y Q, Zhang J S. Solar penetration radiation in snow layers[J]. Journal of Glaciology and Geocryology, 1988, 10(2):135-142. http://www.cnki.com.cn/Article/CJFDTOTAL-BCDT198802004.htm
[11] Kormos P R, Marks D, McNamara J P, et al. Snow distribution, melt and surface water inputs to the soil in the mountain rain-snow transition zone[J]. Journal of Hydrology, 2014, 519:190-204. doi: 10.1016/j.jhydrol.2014.06.051
[12] Schelker J, Kuglerová L, Eklöf K, et al. Hydrological effects of clear-cutting in a boreal forest-snowpack dynamics, snowmelt and streamflow responses[J]. Journal of Hydrology, 2013, 484: 105-114. doi: 10.1016/j.jhydrol.2013.01.015
[13] Michael A, Rawlins, Kyle C, et al. Remote sensing of snow thaw at the pan-Arctic scale using the seawinds scatterometer[J]. Journal of Hydrology, 2005, 312(1-4):294-311. doi: 10.1016/j.jhydrol.2004.12.018
[14] Bernier P Y. Microwave remote sensing of snowpack properties: potential and limitations[J]. Nordic Hydrology, 1987, 18(1): 1-20. http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_bac5a3a177fc5f9ac0af5e29a5f820e4
[15] 王雪芹, 张元明, 蒋进, 等.古尔班通古特沙漠南部沙垄水分动态:兼论积雪融化和冻土变化对沙丘水分分异作用[J].冰川冻土, 2006, 28(2):262-268. doi: 10.3969/j.issn.1000-0240.2006.02.017 Wang X Q, Zhang Y M, Jiang J, et al. Variation pattern of soil water content in longitudinal dune in the southern part of Gurbantêggêt Desert: how snowmelt and frozen soil change affect the soil moisture[J]. Journal of Glaciology & Geocryology, 2006, 28(2):262-268. doi: 10.3969/j.issn.1000-0240.2006.02.017
[16] 王贺, 蔡体久, 满秀玲, 等.小兴安岭不同类型人工林林内积雪特征[J].水土保持学报, 2012, 26(6):263-267, 273. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201206052 Wang H, Cai T J, Man X L, et al. Characteristics of snowpack in different planted forests of Xiaoxing'anling[J]. Journal of Soil and Water Conservation, 2012, 26(6):263-267, 273. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201206052
[17] 李奕, 蔡体久, 盛后财, 等.大兴安岭地区天然樟子松林降雪截留及积雪特征[J].水土保持学报, 2014, 28(5):124-128. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201405022 Li Y, Cai T J, Sheng H C, et al. Characteristics of the snow interception and the snowpack in Scotch pine forest in Great Xing'an Mountains[J]. Journal of Soil and Water Conservation, 2014, 28(5):124-128. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201405022
[18] 俞正祥, 蔡体久, 朱宾宾.大兴安岭北部主要森林类型林内积雪特征[J].北京林业大学学报, 2015, 37(12):100-107. doi: 10.13332/j.1000-1522.20150175 Yu Z X, Cai T J, Zhu B B. Characteristics of snowpack in major forest types of northern Daxing' anling Mountains, northeastern China[J]. Journal of Beijing Forestry University, 2015, 37(12):100-107. doi: 10.13332/j.1000-1522.20150175
[19] Zhang J, Ashjian C, Campbell R, et al. The influence of sea ice and snow cover and nutrient availiability on the formation of massive under-ice phytoplankton blooms in the Chukchi Sea[J]. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 2015, 118:122-135. doi: 10.1016/j.dsr2.2015.02.008
[20] Golding D L, Swanson R H. Snow distribution patterns in clearings and adjacent forest[J]. Water Resources Research, 1986, 22(13): 1931-1940. doi: 10.1029/WR022i013p01931
[21] Martinee J. Expected snow loads on structure from incomplete hydrological data[J]. Journal of Glaciology, 1977, 19(81):185-195. doi: 10.1017/S0022143000029270
[22] Wever N, Schmid L, Heilig A, et al. Verification of the multi-layer SNOWPACK model with different water transport schemes[J]. The Cryosphere, 2015, 9(6): 2271-2293. doi: 10.5194/tc-9-2271-2015
[23] 杨俊华, 秦翔, 吴锦奎, 等.祁连山老虎沟流域春季积雪属性的分布及变化特征[J].冰川冻土, 2012, 34(5):1092-1093. http://d.old.wanfangdata.com.cn/Periodical/bcdt201205010 Yang J H, Qin X, Wu J K, et al. Distribution and variation of spring snow cover in Laohugou Watershed of the Qilian Mountains[J]. Journal of Glaciology and Geocryology, 2012, 34(5):1092-1093. http://d.old.wanfangdata.com.cn/Periodical/bcdt201205010
[24] 高培, 魏文寿, 刘明哲.中国西天山季节性积雪热力特征分析[J].高原气象, 2012, 31(4): 1075-1078. http://d.old.wanfangdata.com.cn/Periodical/gyqx201204022 Gao P, Wei W S, Liu M Z. Characteristic analysis on temperature change in snow layer in Western Tianshan Mountain in China[J]. Plateau Meteorology, 2012, 31(4): 1075-1078. http://d.old.wanfangdata.com.cn/Periodical/gyqx201204022
[25] 陆恒, 魏文寿, 刘明哲, 等.季节性森林积雪融雪期雪层含水率垂直廓线与时间变化特征[J].地理研究, 2011, 30(7):1245-1248. http://d.old.wanfangdata.com.cn/Periodical/dlyj201107010 Lu H, Wei W S, Liu M Z, et al. A study on the vertical profile of forest snow liquid water content and variation with time in the snowmelt period[J]. Geographical Research, 2011, 30(7):1245-1248. http://d.old.wanfangdata.com.cn/Periodical/dlyj201107010
[26] 刘家霖, 满秀玲.降雨和非降雨日兴安落叶松天然林蒸腾及蒸散发特征[J].生态学报, 2017, 37(15):5059-5069. http://d.old.wanfangdata.com.cn/Periodical/stxb201715014 Liu J L, Man X L.Characteristics of transpiration and evapo-transpiration from natural Larix gmelinii forests on rainy and non-rainy days[J]. Acta Ecologica Sinica, 2017, 37(15):5059-5069. http://d.old.wanfangdata.com.cn/Periodical/stxb201715014
[27] 刘海亮, 蔡体久, 满秀玲, 等.小兴安岭主要森林类型对降雪、积雪和融雪过程的影响[J].北京林业大学学报, 2012, 34(2): 20-25. http://j.bjfu.edu.cn/article/id/9722 Liu H L, Cai T J, Man X L, et al. Effects of major forest types of Xiaoxing'an Mountains on the process of snowfall, snow cover and snow melting[J]. Journal of Beijing Forestry University, 2012, 34(2):20-25. http://j.bjfu.edu.cn/article/id/9722
[28] 李辉东, 关德新, 吴家兵, 等.长白山阔叶红松林冬季雪面蒸发特征[J].应用生态学报, 2013, 24(4):1039-1046. http://d.old.wanfangdata.com.cn/Periodical/yystxb201304022 Li H D, Guan D X, Wu J B, et al. Characteristics of evaporation over broadleaved Korean pine forest in Changbai Mountains, Northeast China during snow cover period in winter[J]. Chinese Journal of Applied Ecology, 2013, 24(4):1039-1046. http://d.old.wanfangdata.com.cn/Periodical/yystxb201304022
[29] Mcjannet D, Vertessy R. Effects of thinning on wood production, leaf area index, transpiration and canopy Interception of a plantation subject to drought[J]. Tree Physiology, 2011, 21 (12/13): 1001-1008. http://cn.bing.com/academic/profile?id=de705ad9c3839e0d39309fb425d6fb00&encoded=0&v=paper_preview&mkt=zh-cn
[30] Bulcock H H, Jewitt G P W. Spatial mapping of leaf area index using hyperspectral remote sensing for hydrological application with a particular focus on canopy interception[J]. Hydrology and Earth System Science, 2010, 14(2): 383-392. doi: 10.5194/hess-14-383-2010
[31] Rutter N, Essery R, Pomeroy J, et al. Evaluation of forest snow processes models (SnowMIP2)[J/OL]. Journal of Geophysical Research: Atmospheres, 2009, 114[2017-08-11]. http://onlinelibrary.wiley.com/doi/10.1029/2008JD011063/full.
[32] Varhola A, Coops N C, Weiler M, et al. Forest canopy effects on snow accumulation and ablation: an integrative review of empirical results[J]. Journal of Hydrology, 2010, 392(3): 219-233. http://cn.bing.com/academic/profile?id=e815be270bf9fff351eb43f73e52a547&encoded=0&v=paper_preview&mkt=zh-cn
[33] Boon S. Snow ablation energy balance in a dead forest stand[J]. Hydrological Processes, 2009, 23(18): 2600-2610. doi: 10.1002/hyp.v23:18
[34] Coughlan J C, Running S W. Regional ecosystem simulation: a general model for simulating snow accumulation and melt in mountainous terrain[J]. Landscape Ecology, 1997, 12(3):119-136. doi: 10.1023/A:1007933813251
[35] Lee Y H, Mahrt L. An evaluation of snowmelt and sublimation over short vegetation in land surface modeling[J]. Hydrological Processes, 2004, 18: 3543-3557. doi: 10.1002/hyp.5799
[36] Suzuki K, Nakai Y. Canopy snow influence on water and energy balances in a coniferous forest plantation in northern Japan[J]. Journal of Hydrology, 2008, 352(1): 126-138. http://cn.bing.com/academic/profile?id=483a6d8db5f9163138b01bbf53530c1b&encoded=0&v=paper_preview&mkt=zh-cn
[37] Lundberg A, Calder I, Harding R. Evaporation of intercepted snow: measurement and modeling[J]. Journal of Hydrology, 1998, 206(3-4): 151-163. doi: 10.1016/S0022-1694(97)00016-4
-
期刊类型引用(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)
计量
- 文章访问数: 1405
- HTML全文浏览量: 329
- PDF下载量: 22
- 被引次数: 13
