Distribution characteristics of soil organic carbon and its influencing factors in different vegetation types in loess region of northern Shaanxi Province, northwestern China
-
摘要:目的 探究陕北黄土区退耕还林（草）后形成的主要植被群落类型土壤有机碳空间分布特征及其影响因素，旨在为今后该地区人工林土壤碳汇管理以及生态效益评估提供参考依据，为我国北方森林土壤碳的相关研究积累基础数据。方法 以吴起县大吉沟森林公园内的油松林、沙棘林、草地、油松沙棘混交林为研究对象，选取典型样区，运用单因素方差分析与灰色关联法，探讨不同植被类型在0 ~ 100 cm土壤有机碳垂直变化规律及其主要影响因素。结果 （1）研究区土壤有机碳含量及储量具有明显表聚现象，且随土壤深度增加而降低。（2）不同植被类型下，土壤有机碳平均含量表现为沙棘林（7.03 g/kg） > 低坡度油松沙棘混交林（5.34 g/kg） > 草地（5.16 g/kg） > 高坡度油松沙棘混交林（3.87 g/kg） > 油松林（3 g/kg），沙棘林与油松林、高坡度油松沙棘混交林土壤有机碳平均含量呈显著性差异（P < 0.05）。（3）不同植被类型土壤有机碳储量介于41.11 ~ 74.76 t/hm2。（4）不同植被土壤剖面C/N在16.41 ~ 39.11之间，C/N均值由大到小表现为沙棘林（34.68） > 低坡度油松沙棘混交林（25.88） > 草地（25.82） > 油松林（23.08） > 高坡度油松沙棘混交林（22.71）。（5）不同植被类型土壤理化因子与有机碳含量关联度均在中等关联以上，与有机碳含量关系密切。结论 研究区在今后建设碳汇林时应充分考虑土壤有机碳影响因素，优先选择沙棘等优势树种。Abstract:Objective Our project objective is to explore the spatial distribution characteristics of soil organic carbon and its influencing factors of the main vegetation community types formed after returning farmland to forest land (grassland) in the loess region of northern Shaanxi Province, northwestern China and to evaluate soil carbon sink management and ecological benefits of planted forests in the future, in order to provide reference basis and accumulate basic data for related research on forest soil carbon in northern China.Method Taking the Pinus tabuliformis forest, Hippophae rhamnoides forest, grassland and mixed forest of Pinus tabuliformis and Hippophae rhamnoides in Dajigou Forest Park of Wuqi County as the research objects, the typical sample plots were selected, and the single factor analysis of variance and gray correlation method were used to explore the different types of vegetation in 0−100 cm soil organic carbon vertical change and its main influencing factors.Result (1) The soil organic carbon content and storage in the study area had obvious surface aggregation phenomenon, and decreased with increasing soil depth. (2) The average content of soil organic carbon under different vegetation types was shown as Hippophae rhamnoides forest (7.03 g/kg) > low-slope Pinus tabuliformis-Hippophae rhamnoides mixed forest (5.34 g/kg) > grassland (5.16 g/kg) > high-slope Pinus tabuliformis-Hippophae rhamnoides mixed forest (3.87 g/kg) > Pinus tabulacformis forest (3 g/kg). Significant differences in average soil organic carbon content were found between Hippophae rhamnoides forest and Pinus tabuliformis forest, high-slope Pinus tabuliformis-Hippophae rhamnoides mixed forest (P < 0.05). (3) Soil organic carbon reserve of different vegetation types ranged from 41.11 to 74.76 t/ha. (4) The C/N value of soil profiles of different vegetation ranged from 16.41 to 39.11, and the average C/N values varied from big to small as Hippophae rhamnoides forest (34.68) > low-slope Pinus tabuliformis-Hippophae rhamnoides mixed forest (25.88) > grassland (25.82) > Pinus tabuliformis forest (23.08) > High-slope Pinus tabuliformis-Hippophae rhamnoides mixed forest (22.71). (5) The correlations between soil physical and chemical factors and organic carbon content of different vegetation types are all above medium correlation, and closely related to organic carbon content.Conclusion In the future construction of carbon sink forests in the study area, the influencing factors of soil organic carbon should be fully considered, and preferentially select dominant tree species such as Hippophae rhamnoides forest.
-
类毒过敏原蛋白（venom-allergen proteins，VAPs）是一种富含半胱氨酸的分泌蛋白（cysteine-rich secretory protein，CRISP），属于精子包被蛋白（sperm coating protein）的Tpx-1 / Ag-5 / Pr-1 / Sc-7（SCP / TAPS）超家族[1]。动植物寄生线虫的栖息地不同，但它们的分泌物共享一组结构保守的蛋白质—VAPs。这些蛋白在线虫生命周期的寄生阶段上调表达，而且对线虫在植物和动物体内的持续侵染至关重要[2-3]。寄生线虫在侵染植物过程中使用口针刺穿植物细胞壁，将分泌的VAPs传递到细胞中，并从活植物细胞中获取低分子量化合物。目前，编码VAP蛋白的基因在许多植物寄生线虫中已被分离鉴定，如南方根结线虫（Meloidogyne incognita）[4-5]、大豆孢囊线虫（Heterodera glycines）[6]、马铃薯腐烂茎线虫（Ditylenchus destructor）[7]、松材线虫（Bursaphelenchus xylophilus）[8]、马铃薯金线虫（Globodera rostochiensi）[9]、禾谷胞囊线虫（Heterodera avenae）[10]和香蕉穿孔线虫（Radopholus similis）[11]。目前的研究主要集中在VAPs作为宿主免疫反应的激活剂和抑制剂的作用上。例如，马铃薯金线虫的Gr-VAP-1蛋白与番茄细胞外半胱氨酸蛋白酶Rcr3pim的结合会触发由表面定位受体Cf-2介导的防御相关超敏反应。另外，Gr-VAP-1基因在番茄（Solanum lycopersicum）中还可以抑制细胞外受体蛋白Cf-4及其同源诱导子Avr-4引起的细胞死亡。在转基因植物的细胞外基质中，异位表达的线虫基因VALs（Gr-VAP-1、Hs-VAP-1和Hs-VAP-2）的表达显著提高了其对植物寄生线虫的敏感性[9]。在植物寄生线虫中，VAPs基因的表达与入侵宿主及随后在宿主组织的迁移有关[2,6,12-13]。VAPs基因的RNA干扰（RNA interference，RNAi）会降低植物寄生线虫在宿主植物内迁移并成功建立永久性觅食位点的能力[9,13-14]。禾谷胞囊线虫的2个类毒素过敏原蛋白新基因HaVAP1和HaVAP2均参与了该线虫的寄生，但它们具有不同的功能。RNAi沉默HaVAP1基因增强了禾谷胞囊线虫的毒力，RNAi沉默HaVAP2基因阻碍了禾谷胞囊线虫的寄生[10]。这些研究表明，植物寄生线虫在宿主组织的迁移以及感染的早期阶段，口针分泌的VALs可能通过抑制表面定位受体介导的植物固有免疫力，从而增强宿主植物对线虫感染的总体敏感性，这对于调节宿主反应非常重要。根结线虫的Mi-VAP-1作为已知的效应子，对其沉默后的表达进行转录组分析，得出该基因的沉默对其他效应子表达产生干扰。这种对其他基因的级联效应也可有助于降低线虫的寄生能力[15]。香蕉穿孔线虫的类毒素过敏原蛋白RsVAP可能与LeRabA1d（Ras-related protein）相互作用以影响宿主防御反应，进而促进线虫感染[11]。因此，明确类毒素过敏原蛋白VAPs的结构与功能，有助于了解植物寄生线虫的致病机理，从而设计针对植物寄生线虫寄主的新型保护策略。
松材线虫病（pine wilt disease）即松树（Pinus spp.）萎蔫病，是一种以松材线虫为病原的国际检疫性林业病害[16]，已造成我国生态与林业经济的巨大损失。全国大面积的健康松林均面临着松材线虫入侵扩张的高风险和大流行的威胁[17-18]。2011年首次在松材线虫中发现一个VAP基因家族，包括3个以基因簇形式存在的功能性类毒过敏原蛋白基因（Bx-vap-1、Bx-vap-2、Bx-vap-3）和一个不表达的假基因（Bx-vap-P）[8]。Bx-vap-1与Bx-vap-2基因核苷酸序列不同，但蛋白序列相同。原位杂交显示这3个基因均在松材线虫的食道腺特异表达。免疫组化分析显示Bx-VAP1 蛋白特异存在于松材线虫的中食道球和食道腺中，该基因和蛋白的表达水平均于松材线虫在松树体内生长繁殖阶段达到最高，RNAi 沉默该基因会显著降低松材线虫的迁移率，但对繁殖率没有影响[14]。这些结果表明Bx-VAP1 蛋白可能经口针分泌至体外抑制松树的防卫反应，参与松材线虫的迁移过程。转录组分析发现松材线虫6个类毒过敏原蛋白[19]，为研究松材线虫 VAP 的多态性和功能提供了新的线索。后来，采用重叠延伸PCR 基因合成法人工合成松材线虫Bxvap2基因序列在大肠杆菌（Escherichia coli）内表达，获得原核表达的VAP2蛋白，纯化后注入新西兰白兔体内制备多克隆抗体[20]。利用杆状病毒表达系统在sf9昆虫细胞中表达Bx-vap-1基因[14,21]，获得VAP1蛋白接种的3年生马尾松（Pinus massoniana）后，寄主马尾松细胞发生质壁分离、细胞核降解等现象，表现出典型的细胞程序性死亡[22]。由此，本研究推测松材线虫VAPs作为关键的致病基因之一，在该线虫与寄主植物互作过程中发挥着至关重要的作用。
类毒素过敏原蛋白VAPs可能作为效应子或者病程相关蛋白参与松材线虫致病过程。虽然这类功能蛋白已被筛选出，但是其结构及具体功能尚不明确。因此，本文以松材线虫的4个VAPs 基因为例，检测不同发育阶段的基因表达量，构建原核表达载体，使其在大肠杆菌中得到快速高效表达，并以纯化的原核表达产物为抗原制备相应的多克隆抗体，分析此类蛋白的生物信息学，为进一步研究这类蛋白的组织定位及功能提供实验材料，为明确该基因在松材线虫与寄主植物互作中的作用机理奠定基础。
1. 材料与方法
1.1 供试虫源的准备
供试松材线虫虫株来自于本实验所保存的Nxy61 虫株，该虫株分离于浙江宁波2012 年 10 月当年发病的马尾松枝条，通过贝尔曼漏斗法分离后长期保存于长满灰葡萄孢（Botrytis cinerea）的玉米培养基中。试验接种线虫为混合虫态松材线虫（包括虫卵、2 ~ 4 龄幼虫以及成虫），在灰葡萄孢 PDA 培养基上25 ℃黑暗培养5 d后，贝尔曼漏斗法收集分离线虫，无菌水冲洗3 次后，再用0.5% 的硫酸链霉素溶液处理2 min 后，无菌水冲洗3 次，调整线虫悬浮液浓度为30 头/μL后，备用。
之后，获取不同发育阶段的松材线虫。收集松材线虫的卵，加入10 mL ddH2O于25 ℃处理24 h后，卵孵化为繁殖型2龄线虫（L2），收集L2接种在长满灰葡萄孢的PDA培养基上，24 h后获得繁殖型3龄线虫（L3）；L2培养52 h后获得繁殖型4龄线虫（L4），L2培养72 h后获得成虫。
将收集到的卵、L2、L3、L4、成虫经无菌水冲洗后13 000 rpm离心5 min，充分去除上清液，收集混合虫态松材线虫，液氮速冻后保存于−80 ℃冰箱，用于线虫RNA的提取。
1.2 松材线虫RNA的提取及cDNA的合成
利用RNA提取试剂盒RNeasy Plus Mini Kit（No. 74134；Qiagen，Hilden，Germany）提取松材线虫不同发育阶段卵、L2、L3、L4、成虫的RNA，通过NanoDrop 8000（Thermo，Waltham，MA，USA）检测其纯度和浓度。之后使用PrimeScript RT reagent Kit with gDNA Eraser（No. RR047Q，TaKaRa，Shiga，Japan）试剂盒进行松材线虫cDNA第一链合成，作为实时荧光定量PCR的模板。
1.3 不同发育阶段松材线虫Bx-VAPs基因的表达量
根据NCBI上松材线虫的VAPs基因的mRNA完整编码序列，利用Primer Premier 6.0软件设计实时荧光定量PCR引物（表1），以不同龄期的松材线虫cDNA 为模板，松材线虫β-actin 基因为内参（表1），根据KAPA Roche480 SYBRFast（北京，普凯瑞）操作步骤，使用Roche480 Real-time PCR 扩增仪进行实时荧光定量PCR分析，检测不同发育阶段的松材线虫Bx-VAP1、Bx-VAP2、Bx-VAP3和Bx-VAP4基因表达量。每个样品设3 个重复，以卵作为对照，采用2–ΔΔCt 法计算目的基因的相对表达量，采用SPSS软件进行单因素方差分析。
表 1 引物信息Table 1. Primer information引物名称
Primer name引物序列
Primer sequence引物用途
Primer usageBx-VAP1-F ATGTTCGGTCTACTGTTAGT 全长克隆
Full-length gene cloneBx-VAP1-R TCAGGCGCTACACAGACTGT Bx-VAP2-F ATGGTTCGAGTATTAGTTCT Bx-VAP2-R TCAGGCACTGCACAAACTGT Bx-VAP3-F ATGTTCCGAGTTCTCCTGAC Bx-VAP3-R CTATTCGGCACTGCAGAGTG Bx-VAP4-F ATGATAACTAAATTCCTGTG Bx-VAP4-R TTAAGCGACGCACAATCCCT EX-VAP1-F GTATTTTCAGGGATCCGAATTCGACAGGTTCAGCAACAGCCAA 添加双酶切位点
Addition of double cleavage sitesEX-VAP1-R GGTGGTGGTGGTGGTGCTCGAGTCAGGCGCTACACAGACTGT EX-VAP2-F GTATTTTCAGGGATCCGAATTCACCAAATTCAGCGAAAGCCAA EX-VAP2-R GGTGGTGGTGGTGGTGCTCGAGTCAGGCGCTACACAGACTGT EX-VAP3-F GTATTTTCAGGGATCCGAATTCGCCAGGTTAAGCGACGATGAG EX-VAP3-R GGTGGTGGTGGTGGTGCTCGAGCTATTCGGCACTGCAGAGTG EX-VAP4-F GTATTTTCAGGGATCCGAATTCGAAAAGTTGAATAGCGCTGAG EX-VAP4-R GGTGGTGGTGGTGGTGCTCGAGTTAAGCGACGCACAATCCCT RT-VAP1-F AATGTGATGGAGGCTGAC 实时荧光定量PCR
Real-time quantitative polymerase chain reactionRT-VAP1-R GCTACACAGACTGTCACT RT-VAP2-F GCTGTGCCATTCAGAACT RT-VAP2-R GCTCCCATCATGTTTCCA RT-VAP3-F TTGGACCTTCACCGTCTG RT-VAP3-R CTTCTTCTTCCTCGCACTG RT-VAP4-F CTCCGACCTGGAAGACAA RT-VAP4-R CGCTCTCATGGCAATGTT β-actin -F TCCGTACCCTGAAGTTGGCTAACC β-actin -R AAGTGGAGACGAGGGAATGGAACC 1.4 松材线虫Bx-VAPs基因的克隆
根据NCBI上松材线虫的VAPs基因的mRNA完整编码序列，利用Primer Premier 6.0软件设计克隆引物，由华大蛋白（Beijing，China）合成（表1），使用PrimeSTAR HS（TaKaRa，日本）目的基因全长PCR扩增体系：cDNA模板1 μL，2 × PrimeSTAR Max Premix12.5 μL，正反向引物各1 μL，补ddH2O至25 μL。PCR反应程序：94 ℃预变性5 min；35个循环的94 ℃变性30 s，55 ℃退火10 s，72 ℃延伸1 min；72 ℃再延伸10 min，保存温度为4 ℃。使用北京全式金生物技术有限公司的pEASY-Blunt Simple Cloning Kit （No. CB111-01；TransGen Biotech，Beijing，China）试剂盒进行连接并转化，将得到的阳性克隆菌液，送至北京六合华大基因科技有限公司进行质粒提取和双向测序，获得测序正确的Bx-VAPs基因全长质粒。
1.5 重组质粒pET-32b-VAPs的构建
将pET-32b空载体分别经EcoRⅠ/XhoⅠ双酶切，37 ℃处理3 h，回收相应目的片段。利用Primer Premier 5.0软件设计扩增引物（表1），添加双酶切位点EcoRⅠ（GAATTC）和XhoI（CTCGAG），引物由北京六合华大基因科技有限公司合成。以克隆的VAPs全长质粒为模板、EX-VAPn-F/ EX-VAPn-R为引物分别进行VAPs的扩增。25 µL PCR扩增体系和反应条件同1.4。产物经1%的琼脂糖凝胶电泳检测，在−20 ℃条件下保存。使用PCR 产物纯化试剂盒（EZNA cycle pure，kit vacuum protocal，OMEGA公司）纯化并回收目的片段。使用一步法克隆试剂盒（ClonExpress® II One Step Cloning Kit，Vazyme）将回收目的片段产物与线性化pET-32b载体进行同源重组，并将其转化至DH5α感受态细胞中，次日挑取单克隆菌落进行PCR鉴定，测序正确的菌液扩大培养并提取重组质粒，命名为pET-32b-VAPn。
1.6 重组Bx-VAPs蛋白的诱导表达与纯化
将测序正确的重组质粒pET-32b-VAPn分别转化大肠杆菌 BL21DE3感受态细胞，并涂布于相应LB平板，37 ℃倒置培养16 h。次日挑单菌落进行PCR鉴定正确后，送至六合华大基因科技有限公司测序，并将测序正确的目的菌株Bx-VAPn用15%甘油保存，用于目的蛋白表达。将验证正确的目的菌株接种至Amp+（100 μg/mL）抗性LB液体培养基中37 ℃ 220 rpm振荡培养至对数生长期（OD600介于0.4 ~ 0.6）。加入终浓度为1 mmol/L的IPTG，37 ℃ 180 rpm振荡培养5 h后，于 50 ml离心管中4000 rpm离心15 min，收集诱导前和诱导后菌体，与1 × SDS混合煮沸，进行SDS-聚丙烯酰胺凝胶电泳（SDS-PAGE）检测。诱导后的菌液以0.1 mmol/L PBS洗涤沉淀后再次重悬，冰浴条件下超声波破碎后，4 ℃ 12 000 rpm 离心 15 min，使上清和菌体碎片完全分离，至溶液澄清。收集上清液，经Ni亲和层析预装柱纯化后，分别用含10 mmol/L和250 mmol/L咪唑的缓冲液分别洗脱杂质蛋白和目的蛋白。使用 Ni2+-NTA 琼脂糖凝胶预装柱（全式金）对 VAPs蛋白进行纯化。纯化后的蛋白进行SDS-PAGE分析以判断其纯化效果。
1.7 重组Bx-VAPs蛋白多克隆抗体的制备与鉴定
1.7.1 动物免疫
将纯化后的Bx-VAPn蛋白分别作为免疫原与弗氏完全佐剂，取200 µL 1∶1混匀液（含100 µg的Bx-VAPn重组蛋白）在每只Balb/c小鼠进行背部皮下脊柱两侧均匀免疫（8 ~ 10个点）。初免14 d或21 d后，分别将Bx-VAPn蛋白与等体积的弗氏不完全佐剂混匀后进行二次免疫。二免14 d后，进行三免，方法同二免。三免7 d后，尾静脉采血，3000 ~ 10 000 rpm离心5 min，分离血清（10 ~ 20 µL）。
1.7.2 抗体效价
按《SOP-抗体鉴定》方法进行效价检测：
包被：用 pH9.6 CB的1 × 包被液，按1 ug/mL的抗原质量浓度，4 ℃包被过夜（100 µL/孔），洗板1次，用封闭液200 µL /孔37 ℃封闭2 h，拍干备用；
加样：将多抗血清（1 mg/mL）用样稀液依次稀释100、200、400、800、1 600、3 200、6 400、12 800倍，以阴性血清1 000倍稀释液做阴性对照，100 µL/孔37 ℃孵育30 min，洗板3次拍干；
二抗：取适当稀释的酶标二抗（羊抗鼠IgG-HRP），100 µL /孔37 ℃孵育30 min，洗板5次拍干；
显色：取显色液按100 µL /孔 加入酶标板，37 ℃避光显色15 min；
终止：各反应孔加入终止液 50 µL终止反应；
读值：酶标仪设定450 nm和630 nm双波长读值，采用OD450 ~ OD630 的值为1为抗体；记录并分析数据，1/2最大吸光度对应的稀释倍数为多抗血清效价。
1.8 Bx-VAPs蛋白的蛋白免疫印迹
纯化血清抗体anti-VAPs分别稀释至1 mg/mL后，留存备用。取50 µg 经过蛋白定量处理的松材线虫蛋白液（Bx）和原核表达的蛋白（Bx-VAPn），检测抗体的特异性，进行15% SDS-PAGE电泳，转膜、封闭后以1∶10 000的稀释比孵育对应的多克隆抗体，二抗采用羊抗鼠IgG-HRP（1∶4 000稀释），进行DAB显色压片曝光并拍照保存。
1.9 Bx-VAPs蛋白的生物信息学在线分析
利用ExPASyProtparam在线程序（https://web.expasy.org/protparam/）分析Bx-VAPs蛋白的理化性质，利用SOPMA在线程序预测Bx-VAPs蛋白的二级结构，利用SWISS-MODEL在线程序（https://swissmodel.expasy.org/）预测Bx-VAPs蛋白的三级结构，利用SingalP5.0在线程序（http://www.cbs.dtu.dk/services/SignalP/）预测Bx-VAPs蛋白是否存在信号肽，利用DNAStar软件中Protean软件对Bx-VAPs进行B细胞抗原表位预测与分析。
2. 结果与分析
2.1 松材线虫Bx-VAPs的基因全长克隆
扩增松材线虫的Bx-VAPs基因全长并进行琼脂糖凝胶电泳检测，条带单一。测序分析结果显示克隆所得的基因序列与NCBI和Wormbase上所提交的基因序列（表2）大小一致。Bx-VAP1基因编码区长度为615 bp，Bx-VAP2基因编码区长度为621 bp，Bx-VAP3基因编码区长度为612 bp，Bx-VAP4基因编码区长度为606 bp（图1）。
表 2 松材线虫的Bx-VAPs蛋白序列及理化性质Table 2. Protein sequence and physicochemical properties of Bx-VAPs in B. xylophilus蛋白名称
Protein name基因库登录号
GenBank ID氨基酸数
Number of
amino acid分子量
Molecular
mass/Da理论等电点
Theoretical pI信号肽的位置
Signal peptide
location/aa跨膜结构域
Transmembrane
domain磷酸化位点数量
Number of phosphorylation site丝氨酸
Serine苏氨酸
Threonine络氨酸
ComplexinBx-VAP1 ADG86237.1 204 22 434.86 4.93 17 ~ 18 无 No 12 17 4 Bx-VAP2 ADG86238.1 206 22 418.83 4.68 17 ~ 18 无 No 15 15 4 Bx-VAP3 ADG86239.1 202 22 446.27 4.31 18 ~ 19 无 No 16 11 8 Bx-VAP4 CAD5147609.1 201 22 038.58 5.19 17 ~ 18 无 No 14 17 7 ![]()
图 1 Bx-VAPs基因PCR扩增的电泳图Figure 1. Electrophoresis map of PCR amplification of Bx-VAPs gene2.2 不同发育阶段松材线虫Bx-VAPs基因的表达量
收集各个龄期松材线虫的cDNA，分析Bx-VAPs基因的表达量（图2），结果表明在成虫时期，松材线虫的Bx-VAP1、Bx-VAP2基因相对表达量达到150和120倍。Bx-VAP3的基因表达量在L2和L3时期略有升高，Bx-VAP4基因表达量在不同龄期变化不超过3倍。
![]()
图 2 不同发育阶段松材线虫Bx-VAPs的相对表达量E代表卵，L2 ~ L4 代表繁殖型2 ~ 4龄幼虫，A代表成虫。不同小写字母表示不同发育阶段基因表达量差异显著性(P < 0.05)。下同。E indicates egg stage, L2 − L4 indicate propagative 2nd to 4th instar larvae, A indicates adult. Different lowercase letters indicate significant differences in gene expression at different developmental stages (P < 0.05). The same below.Figure 2. Relative expression level of Bx-VAPs in B. ylophilus at various developmental stages2.3 松材线虫Bx-VAPs蛋白的原核表达
构建pET-32b-VAPn重组表达载体，并转化至DH5α宿主菌中，获得重组蛋白pET-32b-VAPs菌液，添加IPTG诱导后，分别在37 ℃条件下诱导表达重组蛋白5 h。经SDS-PAGE分析显示，pET-32b-VAP1和pET-32b-VAP2在22.0 kDa左右显示条带，pET-32b-VAP3和pET-32b-VAP4在22 ~ 31 kDa间显示条带，且重组蛋白大多以包涵体形式存在于沉淀中；约1 L菌液经诱导后，菌体经亲和层析纯化后蛋白质量浓度可达到 0.6 mg/mL，获得8 mg蛋白，经考马斯亮蓝染色，TBST脱色后，利用软件Image J对其进行灰度分析，纯度达90 %以上，可作为免疫抗原（图3）。
![]()
图 3 重组蛋白在大肠杆菌BL21DE3中表达的SDS-PAGE分析1. 无IPTG诱导的BL21DE3-BxVAPn全菌；2. IPTG诱导的BL21DE3-BxVAPn 全菌；3. 上清液；4. 包涵体；5. 纯化的BxVAPn蛋白；6. 低分子量蛋白标记。1, BL21DE3–BxVAP1 without IPTG; 2, BL21DE3- BxVAP1 induced with IPTG; 3, supernatant; 4, inclusion body; 5, purified BXVAP1 protein; 6, low molecular mass protein marker.Figure 3. SDS-PAGE analysis of the expression of recombinant protein in E. coli BL21DE32.4 anti-VAPs抗体的效价及蛋白免疫印记结果
间接ELISA法测定结果发现：anti-VAP1多抗血清效价达到1∶12 800，anti-VAP2多抗血清效价达到1∶51 200，anti-VAP3多抗血清效价达到1∶25 600，anti-VAP4多抗血清效价达到1∶51 200（图4）。
免疫印迹法Western blot验证抗体特异性结果如图5所示，Bx-VAP1、Bx-VAP2、Bx-VAP3的抗体分别与对应的原核表达的VAP1、VAP2、VAP3反应，同时均与Bx蛋白液产生反应；Bx-VAP4的抗体与原核表达的VAP4产生反应，但并未与Bx蛋白液发生反应。结果表明所制备的Bx-VAP1、Bx-VAP2、Bx-VAP3多克隆抗体特异性好，灵敏度高，具有与松材线虫的蛋白提取液结合的能力。
![]()
图 5 Western blot验证抗体的特异性M. 低分子量蛋白质标记物；Bx. 松材线虫蛋白液。M, low molecular mass protein marker; Bx, protein solution of B. xylophilus.Figure 5. Western blot verifying antibody specificity2.5 松材线虫Bx-VAPs蛋白的功能和结构预测
利用ExPASy Protparam在线程序分析Bx-VAPs蛋白的理化性质，结果显示：Bx-VAP1编码204个氨基酸，Bx-VAP2编码206个氨基酸，Bx-VAP3编码202个氨基酸；Bx-VAP4编码201个氨基酸（图6A），均为疏水性蛋白（图6B）。为进一步研究Bx-VAPs蛋白的进化关系，根据Bioedit多序列比对结果，应用MEGA 6.0软件和N–J邻接法构建Bx-VAPs蛋白的系统进化树。系统进化树（图6C）显示：松材线虫Bx-VAPs蛋白与拟松材线虫（Bursaphelenchus mucronatus）、豆伞滑刃线虫（Bursaphelenchus doui）、马铃薯腐烂茎线虫、南方根结线虫聚为同一小分支，遗传距离较小，亲缘关系较近。
![]()
图 6 Bx-VAPs蛋白生物信息学分析A. Bx-VAPs的氨基酸序列；B. Bx-VAPs蛋白的亲水/疏水性；C. Bx-VAPs蛋白系统发育树。A, amino acid sequence of Bx-VAPs; B, predicted hydrophilicity/ hydrophobicity of Bx-VAPs; C, phylogenetic tree analysis of Bx-VAPs protein.Figure 6. Bioinformatics information analysis of the Bx-VAPs protein利用SingalP5.0在线程序预测Bx-VAPs蛋白均存在信号肽（表2），其中Bx-VAP1的信号肽位置为1 ~ 18号的氨基酸，Bx-VAP2的位置为1 ~ 18号的氨基酸；Bx-VAP3的位置为1 ~ 19号的氨基酸；Bx-VAP4的位置为1 ~ 18号的氨基酸。通过TMHMM2.0在线程序预测跨膜结构域，结果显示（表2），Bx-VAP1、Bx-VAP2、Bx-VAP3、Bx-VAP4蛋白均无跨膜结构域。利用NetPhos3.1Server在线程序预测Bx-VAPs蛋白的磷酸化位点，结果显示（表2），Bx-VAP1蛋白共有33个磷酸化位点，其中丝氨酸位点12个，苏氨酸位点17个，酪氨酸位点4个；Bx-VAP2蛋白共有34个磷酸化位点，其中丝氨酸位点15个，苏氨酸位点15个，酪氨酸位点4个；Bx-VAP3蛋白共有35个磷酸化位点，其中丝氨酸位点16个，苏氨酸位点11个，酪氨酸位点8个；Bx-VAP4蛋白共有38个磷酸化位点，其中丝氨酸位点14个，苏氨酸位点17个，酪氨酸位点7个。
采用SMART软件进行松材线虫Bx-VAPs蛋白的功能结构域分析可知（图7），Bx-VAP1、Bx-VAP2、Bx-VAP3、Bx-VAP4均具有保守的SCP结构域，分别在各自的第23 ~ 180位氨基酸、第23 ~ 180位氨基酸、第23 ~ 176位氨基酸、第23 ~ 175位氨基酸。此外，Bx-VAP2具有一个低复杂性区域，在第190 ~ 205位氨基酸。
利用SWISS-MODEL在线程序对Bx-VAP1蛋白的三级结构进行预测（图8A），结果显示，所得结构有较多的α螺旋和无规则卷曲，与利用SOPMA在线程序预测的Bx-VAPs蛋白的二级结构（图8B）结果一致。Bx-VAP1蛋白二级结构包括α螺旋（h，45.59%）、无规则卷曲（c，39.71%）、延伸链（e，13.24%）、β转角（t，1.47%）4种类型。Bx-VAP2蛋白二级结构包括α螺旋（h，42.72%）、无规则卷曲（c，42.23%）、延伸链（e，12.14%）、β转角（t，2.91%）4种类型。Bx-VAP3蛋白二级结构包括α螺旋（h，40.39%）、无规则卷曲（c，47.29%）、延伸链（e，10.84%）、β转角（t，1.48%）4种类型。Bx-VAP4蛋白二级结构包括α螺旋（h，37.31%）、无规则卷曲（c，43.28%）、延伸链（e，16.92%）、β转角（t，2.49%）4种类型。
![]()
图 8 Bx-VAPs的三级结构（A）和二级结构（B）预测图h代表α-螺旋，c代表无规则卷曲，e代表延伸链，t代表β转角。h represents α-helix, c represents random coil, e represents extended strand, and t represents β turn.Figure 8. Secondary structure (B) and tertiary structure (A) prediction diagram of Bx-VAPs2.6 松材线虫Bx-VAPs蛋白的B细胞表位预测
利用Protean软件，对Bx-VAPs的二级结构、表面可及性、可塑性、B细胞抗原表位等进行了预测与分析（图9）。综合Gamier-Robson和Chou-Fasman 2种方法预测Bx-VAP1蛋白的二级结构（图9A）发现其主要以α螺旋为主，主要分布在24 ~ 31、36 ~ 45、50 ~ 53、58 ~ 80、92 ~ 96、122 ~ 127、135 ~ 139、178 ~ 185氨基酸残基；β转角主要形成于20 ~ 22、83 ~ 84、116 ~ 118、129 ~ 130、142 ~ 143、151 ~ 154、170 ~ 172、189 ~ 191、193 ~ 200氨基酸残基；无规则卷曲主要分布在80 ~ 83、88 ~ 90、103 ~ 111、118 ~ 199、131 ~ 135、155 ~ 158、173 ~ 177氨基酸残基，β折叠较少。利用Kyte-Doolittle方法对Bx-VAP1蛋白的亲水性进行分析（图9A），结果显示，Bx-VAP1亲水性区域主要分布在：17 ~ 31、32 ~ 53、55 ~ 88、104 ~ 145、151 ~ 166、172 ~ 179、196 ~ 197氨基酸残基。可及性较高区域位于Bx-VAP1蛋白分子表面，而可及性较低区域埋藏于分子内部的区域。Bx-VAP1蛋白的表面可及性较高区域是：17 ~ 26、33 ~ 36、42 ~ 47、65 ~ 73、106 ~ 112、114 ~ 131、161 ~ 164氨基酸残基。利用Karplus-Schulz方法对Bx-VAP1蛋白的可塑性进行分析（图9A），结果显示，柔韧性较高区域有14 ~ 27、34 ~ 50、54 ~ 58、66 ~ 72、100 ~ 112、114 ~ 124、127 ~ 131、140 ~ 143、151 ~ 165、186 ~ 200氨基酸残基。这些区域可能具有一定的柔韧性，形成表位的可能性较大，容易与抗体进行嵌合。利用Jameson-Wolf方法对Bx-VAP1蛋白的B细胞抗原表位进行预测（图9A），可看出Bx-VAP1存在10个潜在的B细胞抗原表位位点，预测的B细胞抗原表位区域为16 ~ 29、31 ~ 59、65 ~ 75、78 ~ 91、99 ~ 112、115 ~ 125、126 ~ 133、140 ~ 147、151 ~ 166、186 ~ 204氨基酸残基。结合该蛋白的二级结构、亲水性、可塑性、抗原性和表面可及性进行综合分析，Bx-VAP1蛋白的优势抗原表位可能存在于20 ~ 23、32 ~ 35、46 ~ 49、81 ~ 84、104 ~ 111、116 ~ 119、129 ~ 132、140 ~ 143氨基酸残基及其附近。
![]()
图 9 Bx-VAPs蛋白的B细胞抗原表位预测A. Bx-VAP1蛋白的B细胞表位预测；B. Bx-VAP2蛋白的B细胞表位预测；C. Bx-VAP3蛋白的B细胞表位预测；D. Bx-VAP4蛋白的细胞表位预测。A, B cell epitope prediction of Bx-VAP1 protein; B, B cell epitope prediction of Bx-VAP2 protein; C, B cell epitope prediction of Bx-VAP3 protein; D, B cell epitope prediction of Bx-VAP4 protein.Figure 9. B cell antigenic epitope prediction of Bx-VAPs protein预测Bx-VAP2蛋白的二级结构（图9B），发现其主要以α螺旋为主，主要分布在25 ~ 31、50 ~ 53、58 ~ 77、92 ~ 94、108 ~ 115、122 ~ 126、135 ~ 142氨基酸残基；β转角主要形成于35 ~ 36、84 ~ 86、117 ~ 118、151 ~ 155、164 ~ 165、170 ~ 172、189 ~ 191、193 ~ 196、199 ~ 202氨基酸残基；无规则卷曲主要分布在32 ~ 34、38 ~ 41、80 ~ 82、87 ~ 90、102 ~ 108、119 ~ 121、129 ~ 134、156 ~ 160、173 ~ 175氨基酸残基；β折叠较少，主要形成于182 ~ 185氨基酸残基。利用Kyte-Doolittle方法对Bx-VAP2蛋白的亲水性进行分析（图9B），结果显示，Bx-VAP2亲水性区域主要分布在：20 ~ 54、64 ~ 94、109 ~ 133、151 ~ 166、173 ~ 179、188 ~ 201氨基酸残基。Bx-VAP2蛋白的表面可及性较高区域是：19 ~ 25、32 ~ 38、44 ~ 48、69 ~ 73、80 ~ 84、114 ~ 129、140 ~ 143、154 ~ 155、163 ~ 165氨基酸残基。利用Karplus-Schulz方法对Bx-VAP2蛋白的可塑性进行分析（图9B），结果显示，柔韧性较高区域有19 ~ 26、33 ~ 51、53 ~ 58、68 ~ 72、81 ~ 90、101 ~ 107、14 ~ 123、141 ~ 143、151 ~ 164、175 ~ 176、186 ~ 201氨基酸残基。利用Jameson-Wolf方法对Bx-VAP2蛋白的B细胞抗原表位进行预测（图9B），可看出Bx-VAP2存在10个潜在的B细胞抗原表位位点，预测的B细胞抗原表位区域：19 ~ 26、32 ~ 74、81 ~ 89、100 ~ 107、113 ~ 125、140 ~ 146、150 ~ 167、170 ~ 179、186 ~ 205氨基酸残基。综合分析Bx-VAP2蛋白的B细胞优势抗原表位可能存在于35 ~ 37、83 ~ 86、162 ~ 165氨基酸残基及其附近。
预测Bx-VAP3蛋白的二级结构（图9C），发现其主要以α螺旋为主，主要分布在35 ~ 37、151 ~ 166、189 ~ 193氨基酸残基；β转角主要形成于46 ~ 47、82 ~ 84、86 ~ 88、99 ~ 101、111 ~ 113、123 ~ 126、137 ~ 140、151 ~ 166、158 ~ 161、170 ~ 171、180 ~ 183氨基酸残基；无规则卷曲主要分布在31 ~ 33、39 ~ 44、48 ~ 58、80 ~ 81、90 ~ 91、95 ~ 96、102 ~ 103、117 ~ 121、128 ~ 129、152 ~ 153、156 ~ 157、172 ~ 176、187 ~ 188氨基酸残基；β折叠较少，主要形成于59 ~ 60、130 ~ 133氨基酸残基。利用Kyte-Doolittle方法对Bx-VAP3蛋白的亲水性进行分析（图9C），结果显示，Bx-VAP3亲水性区域主要分布在：16 ~ 69、76 ~ 127、136 ~ 142、150 ~ 166、168 ~ 199氨基酸残基。Bx-VAP3蛋白的表面可及性较高区域是：20 ~ 26、31 ~ 39、42 ~ 51、61 ~ 67、80 ~ 90、95 ~ 102、105 ~ 122、138 ~ 140、159 ~ 162氨基酸残基。利用Karplus-Schulz方法对Bx-VAP3蛋白的可塑性进行分析（图9C），结果显示，柔韧性较高区域有7 ~ 9、21 ~ 26、31 ~ 56、63 ~ 66、82 ~ 90、98 ~ 101、103 ~ 107、109 ~ 120、123 ~ 127、137 ~ 140、147 ~ 159、169 ~ 197氨基酸残基。利用Jameson-Wolf方法对Bx-VAP3蛋白的B细胞抗原表位进行预测（图9C），可看出Bx-VAP3存在10个潜在的B细胞抗原表位位点，预测的B细胞抗原表位区域：18 ~ 29、32 ~ 59、62 ~ 71、80 ~ 91、98 ~ 128、138 ~ 163、170 ~ 202氨基酸残基。综合分析Bx-VAP3蛋白的优势抗原表位可能存在于31 ~ 33、48 ~ 58、80 ~ 81、117 ~ 121、184氨基酸残基及其附近。
预测Bx-VAP4蛋白的二级结构（图9D），发现其主要以α螺旋为主，主要分布在16 ~ 26、67 ~ 78、120 ~ 122、135 ~ 141、189 ~ 193氨基酸残基；β转角主要形成于33 ~ 34、47 ~ 49、82 ~ 84、86 ~ 88、159 ~ 161、169 ~ 170、184 ~ 186、195 ~ 197氨基酸残基；无规则卷曲主要分布在41 ~ 46、54 ~ 56、80 ~ 81、99 ~ 101、108 ~ 110、129 ~ 131氨基酸残基；β折叠，主要形成于5 ~ 12、57 ~ 61、142 ~ 150、155 ~ 157、162 ~ 166、179 ~ 182、198 ~ 201氨基酸残基。利用Kyte-Doolittle方法对Bx-VAP4蛋白的亲水性进行分析（图9D），结果显示，Bx-VAP4亲水性区域主要分布在：17 ~ 53、58 ~ 88、90 ~ 93、102 ~ 143、183 ~ 195氨基酸残基。Bx-VAP4蛋白的表面可及性较高区域是：18 ~ 24、32 ~ 39、61 ~ 73、81 ~ 84、116 ~ 121、140 ~ 141、188 ~ 190氨基酸残基。利用Karplus-Schulz方法对Bx-VAP4蛋白的可塑性进行分析（图9D），结果显示，柔韧性较高区域有15 ~ 26、32 ~ 49、53 ~ 56、62 ~ 72、81 ~ 85、87 ~ 89、99 ~ 109、114 ~ 123、128 ~ 132、136 ~ 143、152 ~ 154、175 ~ 178、182 ~ 192氨基酸残基。利用Jameson-Wolf方法对Bx-VAP4蛋白的B细胞抗原表位进行预测（图9D），可看出Bx-VAP4存在10个潜在的B细胞抗原表位位点，预测的B细胞抗原表位区域：15 ~ 24、32 ~ 49、53 ~ 58、63 ~ 75、77 ~ 89、100 ~ 110、113 ~ 123、137 ~ 144、169 ~ 171、175 ~ 178、184 ~ 198氨基酸残基。结合该蛋白的二级结构、亲水性、可塑性、抗原性和表面可及性进行综合分析，Bx-VAP4蛋白的优势抗原表位可能存在于33 ~ 36、82 ~ 84氨基酸残基及其附近。
3. 结论与讨论
本研究成功构建4个不同的VAPs原核表达载体，并在大肠杆菌中获得分子量为22 ~ 31 kDa的VAPs蛋白，其中原核表达的Bx-VAP1蛋白与真核表达的VAP1蛋白分子量23.3 kDa[21]和22 kDa[14]相近，充分证明了VAP1蛋白原核表达成功。Bx-VAP2原核表达蛋白与人工合成的VAP2基因在大肠杆菌中原核表达获得的VAP2蛋白分子量（21.8 KDa）相似[20]。Bx-VAP3和Bx-VAP4蛋白首次在大肠杆菌中表达成功。同时，Bx-VAP1、Bx-VAP2、Bx-VAP3、Bx-VAP4基因在卵、L2、L3、L4、成虫中均有表达，其中Bx-VAP1、Bx-VAP2和Bx-VAP3基因的表达与前人研究类似[8]，未发表的转录组数据中显示VAP1基因和VAP2基因在成虫表达量最高，分别为对照组的419.18倍和18.03倍，VAP3和VAP4基因在L3中表达量较高，表达量分别为8.5和6.4倍。推测VAP1和VAP2可能与成虫的发育有关，VAP3和VAP4可能与L3的发育有关，这些基因的持续性表达模式与松材线虫这类迁移性内寄生线虫整个生活史过程中的侵染、迁徙和取食行为之间是否存在必然联系需要深入研究。考虑到mRNA与蛋白质表达水平存在不一致性，我们有必要针对已获得的松材线虫VAPs蛋白进行更加详尽的表达和定位分析，从而为进一步探讨该类基因的功能提供一定的依据。
大肠杆菌表达系统具有遗传背景清晰、操作过程简单、生长繁殖快、产物稳定性好、成本低、产量高、不易污染等优点，是基因表达技术中最早发展且应用最广的经典表达系统[21]。目前，已得到了广泛应用并取得了巨大的科研价值和经济效益。原核表达中信号肽序列的高疏水性会抑制大肠杆菌细胞的生长和目的蛋白的表达，但是不会对表达蛋白的生物学活性造成影响[23-24]，因此本研究设计了去信号肽的酶切引物，以满足松材线虫VAPs的免疫学检测的实验需求。同时，本研究利用纯化后的Bx-VAPs分别进行了Balb/c小鼠免疫，获取了高效价的抗体anti-VAP1、anti-VAP2和anti-VAP3，并通过蛋白免疫印迹成功检测松材线虫的蛋白液，但anti-VAP4抗体未能与松材线虫蛋白液结合反应，这可能是由于Bx-VAP4基因在扩散型4龄LⅣ中大量表达（未发表转录组数据），而本研究中的松材线虫蛋白液是室内培养的繁殖型不同虫态的松材线虫混合液。VAPs是松材线虫的分泌蛋白，该类蛋白在寄主体内的分布和流动是值得关注的问题。得到的重组VAPs分别作为抗原，获得了相应的抗体，使用这些相应的抗体结合ELISA 或免疫印迹等方法对相关病害木进行检测，结合免疫荧光及激光共聚焦对该蛋白在寄主中位置进行鉴定，为后续体外结合实验验证松材线虫的分泌蛋白-寄主植物蛋白相互作用、质谱分析深入研究 VAPs的作用机制及蛋白间的相互作用提供了有力的实验材料基础。
生物信息学分析结果显示，这4个Bx-VAPs蛋白均无跨膜结构域，蛋白质的二级结构有α螺旋、β折叠、β转角、无规则卷曲4种，其中α螺旋和β折叠常位于蛋白质的内部，因此不易与抗体结合，β转角和无规则卷曲更容易突出在蛋白质的表面，与抗体嵌合可能性很大，成为候选抗原表位的可能性更大。运用DNAStar软件对Bx-VAPs蛋白的二级结构、亲水性、表面可及性、柔韧性、B细胞抗原表位等进行分析与预测，综合分析得出Bx-VAP1蛋白的优势B细胞抗原表位较多。该预测结果在后续研究中可通过人工合成表位多肽，或者将抗原表位多肽与其他抗原性较大的蛋白相偶联来增加其免疫原性，以进行免疫学验证。4个Bx-VAPs蛋白均具有SCP结构域，属于前人报道的SCP/TAPS家族蛋白，这类蛋白的功能各不相同，如免疫应答、生殖、病原体入侵等。松材线虫Bx-VAPs蛋白的功能还有待进一步研究。
综上所述，本研究首次通过原核表达体系，诱导表达和纯化多个Bx-VAPs蛋白，并对其基本理化性质、结构与功能特性、潜在B细胞抗原表位进行生物信息学分析，制备了多克隆抗体，对松材线虫蛋白也有一定的免疫反应，为进一步制备Bx-VAPs的单抗、开发松材线虫诊断试剂盒及深入认识 VAPs蛋白在松材线虫致病过程中的作用奠定了基础。
-
![]()
图 2 不同植被类型土壤有机碳储量
不同小写字母表示同种植被类型土壤有机碳储量在不同土层间差异显著(P < 0.05)。Different lowercase letters indicate that soil organic carbon storage of the same plantation type has significant differences among varied soil layers (P < 0.05).
Figure 2. Soil organic carbon storage under different vegetation types
表 1 样地基本情况
Table 1 Basic situation of the sample plots
样地编号
Sample plot No.植被类型
Vegetation type密度/(株·hm−2) Density/(tree·ha−1) 平均树高
Mean tree
height/m坡度
Slope/(°)海拔
Altitude/m油松
Pinus tabuliformis沙棘
Hippophae rhamnoidesPT 油松 Pinus tabuliformis 1 200 3.33 17 1 386 HR 沙棘 Hippophae rhamnoides 2 300 2.62 17 1 406 HPⅠ 油松沙棘(低坡度)
Pinus tabuliformis-Hippophae rhamnoides (low slope)800 1 500 1.88 12 1 396 HPⅡ 油松沙棘(高坡度)
Pinus tabuliformis-Hippophae rhamnoides (high slope)700 1 400 1.98 29 1 380 LL 达乌里胡枝子、赖草
Lespedeza davurica, Leymus secalinus28 1 398 
下载: 导出CSV
表 2 不同植被类型土壤有机碳含量统计特征
Table 2 Statistical characteristics of soil organic carbon content in different vegetation types
植被类型
Vegetation type平均值
Mean value/(g·kg−1)最大值
Max. value/(g·kg−1)最小值
Min. value/(g·kg−1)标准差
Standard deviation标准误
Standard error变异系数
Variation coefficient/%HR 7.03a 13.96 3.19 0.2 0.08 53 PT 3b 5.53 2.31 0.18 0.07 38 LL 5.16ab 9.44 3.19 0.17 0.06 40 HPⅠ 5.34ab 8.92 3.24 0.2 0.06 44 HPⅡ 3.87b 5.41 3.14 0.23 0.07 20 注：表中不同字母表示不同植被类型的土壤有机碳含量差异显著(P < 0.05)。Note: different letters indicate significant differences in soil organic carbon content of varied vegetation types (P < 0.05).
下载: 导出CSV
表 3 不同植被类型土壤理化因子与土壤有机碳的关联度
Table 3 Correlations between soil physical and chemical factors and soil organic carbon under different vegetation types
项目 Item HPⅠ HPⅡ PT LL HR 关联度
Correlation
degree排序
Rank关联度
Correlation
degree排序
Rank关联度
Correlation
degree排序
Rank关联度
Correlation
degree排序
Rank关联度
Correlation
degree排序
Rank土壤密度
Soil bulk density0.78 3 0.77 2 0.76 2 0.59 9 0.94 1 毛管孔隙度
Capillary porosity0.65 7 0.64 7 0.57 9 0.89 2 0.63 9 非毛管孔隙度
Non-capillary porosity0.68 5 0.76 3 0.63 8 0.78 3 0.78 5 饱和含水量
Saturated moisture0.58 12 0.63 8 0.66 6 0.90 1 0.62 10 毛管持水量
Capillary water capacity0.61 10 0.62 9 0.56 10 0.89 2 0.64 8 电导率
Conductivity0.78 3 0.77 2 0.77 1 0.61 8 0.73 6 pH 0.66 6 0.83 1 0.69 4 0.70 5 0.63 9 速效氮
Available nitrogen0.6 11 0.65 6 0.63 8 0.72 4 0.68 7 速效磷
Available phosphorus0.64 8 0.68 5 0.67 5 0.62 7 0.68 7 速效钾
Available potassium0.72 4 0.67 6 0.71 3 0.78 3 0.80 4 全氮
Total nitrogen0.62 9 0.7 4 0.65 7 0.70 5 0.83 2 全磷
Total phosphorus0.86 1 0.56 10 0.65 7 0.54 10 0.59 11 全钾
Total potassium0.84 2 0.76 3 0.51 11 0.63 6 0.81 3
下载: 导出CSV
-
[1] Batjes N H. Total carbon and nitrogen in the soils of the world[J]. European Journal of Soil Science, 2014, 65(1): 10−21. doi: 10.1111/ejss.12114_2.
[2] Li Q, Yu P, Li G, et al. Overlooking soil erosion induces underestimation of the soil C loss in degraded land[J]. Quaternary International, 2014, 349(28): 287−290.
[3] Chen D M, Lan Z C, Hu S J, et al. Effects of nitrogen enrichment on belowground communities in grassland: relative role of soil nitrogen availability vs. soil acidification[J]. Soil Biology and Biochemistry, 2015, 89(10): 99−108.
[4] Post W M, Kwon K C. Soil carbon sequestration and land-use change: processes and potential[J]. Global Change Biology, 2000, 6(3): 317−327. doi: 10.1046/j.1365-2486.2000.00308.x.
[5] 陈心桐, 徐天乐, 李雪静, 等. 中国北方自然生态系统土壤有机碳含量及其影响因素[J]. 生态学杂志, 2009, 38(4):1133−1140. Chen X T, Xu T L, Li X J, et al. Soil organic carbon concentrations and the influencing factors in natural ecosystems of northern China[J]. Chinese Journal of Ecology, 2009, 38(4): 1133−1140.
[6] 郭璐璐, 李安迪, 商宏莉, 等. 川西贡嘎山不同森林生态系统土壤有机碳垂直分布与组成特征与组成特征[J]. 中国农业气象, 2018, 39(10):636−643. doi: 10.3969/j.issn.1000-6362.2018.10.002. Guo L L, Li A D, Shang H L, et al. Total and labile organic carbon in soils of three subalpine forest types in Gongga Mountain, western Sichuan[J]. Chinese Journal of Agrometeorology, 2018, 39(10): 636−643. doi: 10.3969/j.issn.1000-6362.2018.10.002.
[7] 丁越岿, 杨劼, 宋炳煜, 等. 不同植被类型对毛乌素沙地土壤有机碳的影响[J]. 草业学报, 2012, 21(2):18−25. doi: 10.11686/cyxb20120203. Ding Y K, Yang J, Song B Y, et al. Effect of different vegetation types on soil organic carbon in Mu Us Desert[J]. Acta Prataculturae Sinica, 2012, 21(2): 18−25. doi: 10.11686/cyxb20120203.
[8] 李裕元, 邵明安, 郑纪勇, 等. 黄土高原北部草地的恢复与重建对土壤有机碳的影响[J]. 生态学报, 2007, 27(6):2279−2287. doi: 10.3321/j.issn:1000-0933.2007.06.017. Li Y Y, Shao M A, Zheng J Y, et al. Impact of grassland recovery and reconstruction on soil organic carbon in the northern Loess Plateau[J]. Acta Ecologica Sinica, 2007, 27(6): 2279−2287. doi: 10.3321/j.issn:1000-0933.2007.06.017.
[9] 张宏, 黄懿梅, 安韶山, 等. 黄土高原森林带植被群落下土壤活性有机碳研究[J]. 水土保持研究, 2013, 20(3):65−70,77. Zhang H, Huang Y M, An S S, et al. Soil active organic carbon with different plant communities on the Loess Plateau[J]. Research of Soil and Water Conservation, 2013, 20(3): 65−70,77.
[10] 孟国欣, 查同刚, 张晓霞, 等. 植被类型和地形对黄土区退耕地土壤有机碳垂直分布的影响[J]. 生态学杂志, 2017, 36(9):2447−2454. Meng G X, Zha T G, Zhang X X, et al. Effects of vegetation type and terrain on vertical distribution of soil organic carbon on abandoned farmlands in the Loess Plateau[J]. Chinese Journal of Ecology, 2017, 36(9): 2447−2454.
[11] 王鑫, 杨德刚, 熊黑钢, 等. 新疆不同植被类型土壤有机碳特征[J]. 干旱区研究, 2017, 34(4):782−788. Wang X, Yang D G, Xiong H G, et al. Characteristics of soil organic carbon under different vegetation types in Xinjiang[J]. Arid Zone Research, 2017, 34(4): 782−788.
[12] 武小钢, 郭晋平, 杨秀云, 等. 芦芽山典型植被土壤有机碳剖面分布特征及碳储量[J]. 生态学报, 2011, 31(11):3009−3019. Wu X G, Guo J P, Yang X Y, et al. Soil organic carbon storage and profile inventory in the different vegetation types of Luya Mountain[J]. Acta Ecologica Sinica, 2011, 31(11): 3009−3019.
[13] 弓文艳, 陈丽华, 郑学良. 基于不同林分类型下土壤碳氮储量垂直分布[J]. 水土保持学报, 2019, 33(1):152−157,164. Gong W Y, Chen L H, Zheng X L. Vertical distributions of soil carbon and nitrogen reserves in different forests[J]. Journal of Soil and Water Conservation, 2019, 33(1): 152−157,164.
[14] 刘小璐. 生态建设条件下坡面土壤水分消耗与补偿特征模拟研究[D]. 西安: 西安理工大学, 2019. Liu X L. Simulation study on soil water consumption and compensation characteristics of slope unoer ecological construction conditions[D]. Xi’an: Xi’an University of Technology, 2019.
[15] 王秋兵, 段迎秋, 魏忠义, 等. 沈阳市城市土壤有机碳空间变异特征研究[J]. 土壤通报, 2009, 40(2):252−257. Wang Q B, Duan Y Q, Wei Z Y, et al. Spatial varibility of urban soil organic carbon in Shenyang City[J]. Chinese Journal of Soil Science, 2009, 40(2): 252−257.
[16] Deng J L. Introduction to grey system theory[J]. The Journal of Grey System, 1989, 1(1): 1−24.
[17] 罗上华, 马蔚纯, 王祥荣, 等. 城市环境保护规划与生态建设指标体系实证[J]. 生态学报, 2003, 23(1):45−55. doi: 10.3321/j.issn:1000-0933.2003.01.007. Luo S H, Ma W C, Wang X R, et al. A case study on indicator system of urban environmental protection and ecological construction[J]. Acta Ecologica Sinica, 2003, 23(1): 45−55. doi: 10.3321/j.issn:1000-0933.2003.01.007.
[18] 王文静, 王百田, 吕钊, 等. 山西太岳山不同林分土壤有机碳储量研究[J]. 干旱区资源与环境, 2013, 27(1):81−85. Wang W J, Wang B T, Lü z, et al. Soil organic carbon reserve of different forests in Taiyue Mountain[J]. Journal of Arid Land Resources and Environment, 2013, 27(1): 81−85.
[19] 周莉, 李保国, 周广胜. 土壤有机碳主导影响因子及其研究进展[J]. 地球科学进展, 2005, 20(1):99−105. doi: 10.3321/j.issn:1001-8166.2005.01.016. Zhou L, Li B G, Zhou G S. Advances in controlling factors of soil organic carbon[J]. Advances in Earth Science, 2005, 20(1): 99−105. doi: 10.3321/j.issn:1001-8166.2005.01.016.
[20] 赵俊勇, 孙向阳, 李素艳, 等. 辽宁省老秃顶子不同林分类型土壤有机碳储量和碳氮垂直分布特征[J]. 东北林业大学学报, 2016, 44(10):65−68,78. doi: 10.3969/j.issn.1000-5382.2016.10.013 Zhao J Y, Sun X Y, Li S Y, et al. Soil organic carbon storage and vertical distribution of carbon and nitrogen under different forest types in the Laotudingzi Mountain of Liaoning[J]. Journal of Northeast Forestry University, 2016, 44(10): 65−68,78. doi: 10.3969/j.issn.1000-5382.2016.10.013
[21] 李龙, 姚云峰, 秦富仓, 等. 黄花甸子流域人工林土壤有机碳密度分布特征[J]. 西北农林科技大学(自然科学版), 2016, 44(2):77−82. Li L, Yao Y F, Qin F C, et al. Distribution characteristics of soil organic carbon density of different forests in Huanghuadianzi Watershed[J]. Journal of Northwest A&F University (Natural Science Edition), 2016, 44(2): 77−82.
[22] 权伟, 郑方东, 戎建涛. 浙江乌岩岭7种林分土壤碳密度及碳氮比分布特征[J]. 南京林业大学学报(自然科学版), 2019, 43(1):175−180. Quan W, Zheng F D, Rong J T. Soil carbon density and C/N distribution of seven forest types in Wuyanling Nature Reserve, Zhejiang Province[J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2019, 43(1): 175−180.
[23] 黄从德, 张健, 杨万勤, 等. 四川森林土壤有机碳储量的空间分布特征[J]. 生态学报, 2009, 29(3):1217−1225. doi: 10.3321/j.issn:1000-0933.2009.03.017. Huang C D, Zhang J, Yang W Q, et al. Spatial distribution characteristics of forest soil organic carbon stock in Sichuan Province[J]. Acta Ecologica Sinica, 2009, 29(3): 1217−1225. doi: 10.3321/j.issn:1000-0933.2009.03.017.
[24] 文伟, 彭友贵, 谭一凡, 等. 深圳市森林土壤主要类型有机碳分布特征[J]. 西南林业大学学报(自然科学), 2018, 38(6):106−113. Wen W, Peng Y G, Tan Y F, et al. Distribution characteristics of organic carbon in main forest soil types in Shenzhen City[J]. Journal of Southwest Forestry University (Natural Sciences), 2018, 38(6): 106−113.
[25] 朱丽琴, 黄荣珍, 段洪浪, 等. 红壤侵蚀地不同人工恢复林对土壤总有机碳和活性有机碳的影响[J]. 生态学报, 2017, 37(1):249−257. Zhu L Q, Huang R Z, Duan H L, et al. Effects of artificially restored forests on soil organic carbon and active organic carbon in eroded red soil[J]. Acta Ecologica Sinica, 2017, 37(1): 249−257.
[26] 周志文. 我国中东部不同气候带森林土壤有机碳氮分布特征研究[D]. 南京: 南京农业大学, 2015. Zhou Z W. A syudy on the distribution of forest soil organic carbon and total nitrogen in different climatic zones of the middle east China[D]. Nanjing: Nanjing Agricultural University, 2015.
[27] 艾宁, 魏天兴, 朱清科. 陕北黄土高原不同植被类型下降雨对坡面径流侵蚀产沙的影响[J]. 水土保持学报, 2013, 27(2):26−30. Ai N, Wei T X, Zhu Q K. The effect of rainfall for runoff-erosion-sediment yield under the different vegetation types in Loess Plateau of northern Shaanxi Province[J]. Journal of Soil and Water Conservation, 2013, 27(2): 26−30.
[28] 李龙, 周飞, 田杰, 等. 地形因子对半干旱地区土壤有机碳含量的影响[J]. 北方园艺, 2019, 43(16):104−109. Li L, Zhou F, Tian J, et al. Effects of topographic tactors on soil organic carbon content in semi-arid regions[J]. Northern Horticulture, 2019, 43(16): 104−109.
[29] 卓志清, 李勇, 兴安, 等. 东北旱作区土壤碳氮磷生态化学计量特征及其影响因素[J]. 农业机械学报, 2019, 50(10):259−268, 336. doi: 10.6041/j.issn.1000-1298.2019.10.030. Zhuo Z Q, Li Y, Xing A, et al. Characteristic of ecological stoichiometry of soil C, N and P and its influencing factors in dry farming region of northeast China[J]. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(10): 259−268, 336. doi: 10.6041/j.issn.1000-1298.2019.10.030.
[30] 吴晓玲, 张世熔, 蒲玉琳, 等. 川西平原土壤微生物生物量碳氮磷含量特征及其影响因素分析[J]. 中国生态农业学报(中英文), 2019, 27(10):1607−1616. Wu X L, Zhang S R, Pu Y L, et al. Distribution characteristics and impact factors of soil microbial biomass carbon, nitrogen and phosphorus in western Sichuan plain[J]. Chinese Journal of Eco-Agriculture, 2019, 27(10): 1607−1616.
[31] 韩磊, 庄涛, 周慧华, 等. 小清河滨岸带土壤碳氮变化及影响因素研究[J]. 环境科学与技术, 2019, 42(6):28−34. Han L, Zhuang T, Zhou H H, et al. Study on the change of soil carbon and nitrogen and their influencing factors in the riparian zone of Xiaoqing River[J]. Environmental Science & Technology, 2019, 42(6): 28−34.
[32] 王云琦, 王玉杰. 三峡库区林地土壤有机碳含量特征及效应[J]. 长江流域资源与环境, 2010, 19(12):1448−1455. Wang Y Q, Wang Y J. Content of soil organic carbon in forest soil and its effects in the Three Gorges Reservoir Area[J]. Resources and Environment in the Yangtze Basin, 2010, 19(12): 1448−1455.
[33] 赵维俊, 孟好军, 马剑, 等. 祁连山哈溪林区森林土壤电导率剖面变化特征[J]. 林业科技通讯, 2018(11):7−10. Zhao W J, Meng H J, Ma J, et al. Change features of forest soil conductivity profile in Haxi District, Qilian Mountain[J]. Forest Science and Technology, 2018(11): 7−10.
[34] 高艺宁, 许丽, 林凤友, 等. 矿区复垦地土壤有机碳分布及与土壤化学特性的关系[J]. 内蒙古农业大学学报(自然科学版), 2016, 37(1):54−60. Gao Y N, Xu L, Lin F Y, et al. Relationship between soil organic carbon distribution and soil chemical properties in mining area[J]. Journal of Inner Mongolia Agricultural University (Natural Science Edition), 2016, 37(1): 54−60.
[35] 张鹏, 张涛, 陈年来. 祁连山北麓山体垂直带土壤碳氮分布特征及影响因素[J]. 应用生态学报, 2009, 20(3):518−524. Zhang P, Zhang T, Chen N L. Vertical distribution patterns of soil organic carbon and total nitrogen and related affecting factors along northern slope of Qilian Mountains[J]. Chinese Journal of Applied Ecology, 2009, 20(3): 518−524.
-
期刊类型引用(6)
1. 莫崇杏，董明亮，李荣生，余纽，郑显澄，杨锦昌. 米老排杂交子代苗期生长性状遗传变异及选择. 森林与环境学报. 2023(05): 555-560 . 
百度学术
2. Shuchun Li，Jiaqi Li，Yanyan Pan，Xiange Hu，Xuesong Nan，Dan Liu，Yue Li. Variation analyses of controlled pollinated families and parental combining ability of Pinus koraiensis. Journal of Forestry Research. 2021(03): 1005-1011 . 必应学术
3. 潘艳艳，许贵友，董利虎，王成录，梁德洋，赵曦阳. 日本落叶松全同胞家系苗期生长性状遗传变异. 南京林业大学学报(自然科学版). 2019(02): 14-22 . 百度学术
4. 秦光华，宋玉民，乔玉玲，于振旭，彭琳. 旱柳苗高年生长与气象因子的灰色关联度. 东北林业大学学报. 2019(05): 42-45+51 . 百度学术
5. 李峰卿，陈焕伟，周志春，楚秀丽，徐肇友，肖纪军. 红豆树优树种子和幼苗性状的变异分析及优良家系的初选. 植物资源与环境学报. 2018(02): 57-65 . 百度学术
6. 张素芳，张磊，赵佳丽，张莉，张含国. 长白落叶松小RNA测序和其靶基因预测. 北京林业大学学报. 2016(12): 64-72 . 本站查看
其他类型引用(6)
计量
- 文章访问数: 1631
- HTML全文浏览量: 737
- PDF下载量: 112
- 被引次数: 12
