• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

2017年福建地区PRRSV分子检测及ORF5基因变异分析

韦应芳, 冀池海, 孙彦阔, 韩晓亮, 陈尧, 张桂红

韦应芳, 冀池海, 孙彦阔, 等. 2017年福建地区PRRSV分子检测及ORF5基因变异分析[J]. 华南农业大学学报, 2019, 40(3): 14-19. DOI: 10.7671/j.issn.1001-411X.201810027
引用本文: 韦应芳, 冀池海, 孙彦阔, 等. 2017年福建地区PRRSV分子检测及ORF5基因变异分析[J]. 华南农业大学学报, 2019, 40(3): 14-19. DOI: 10.7671/j.issn.1001-411X.201810027
WEI Yingfang, JI Chihai, SUN Yankuo, et al. Molecular detection of PRRSV and variation analysis of ORF5 gene in Fujian province in 2017[J]. Journal of South China Agricultural University, 2019, 40(3): 14-19. DOI: 10.7671/j.issn.1001-411X.201810027
Citation: WEI Yingfang, JI Chihai, SUN Yankuo, et al. Molecular detection of PRRSV and variation analysis of ORF5 gene in Fujian province in 2017[J]. Journal of South China Agricultural University, 2019, 40(3): 14-19. DOI: 10.7671/j.issn.1001-411X.201810027

2017年福建地区PRRSV分子检测及ORF5基因变异分析

基金项目: 国家重点研发计划项目(2018YFD0501200);国家生猪现代农业产业技术体系项目(CARS-35);国家自然科学基金(31272564)
详细信息
    作者简介:

    韦应芳(1994—),女,硕士研究生,E-mail:1764479177@qq.com

    通讯作者:

    张桂红(1968—),女,教授,博士,E-mail: guihongzh@scau.edu.cn

  • 中图分类号: S855.3

Molecular detection of PRRSV and variation analysis of ORF5 gene in Fujian province in 2017

  • 摘要:
    目的 

    研究2017年福建地区中小规模猪场猪繁殖与呼吸综合征病毒(PRRSV)毒株的分子流行情况。

    方法 

    从福建地区中小规模猪场收集83份疑似蓝耳阳性病料,对编码GP5蛋白的ORF5基因进行了RT-PCR检测和遗传进化分析。

    结果 

    采集的83份疑似蓝耳病料中阳性率为73%,测得的序列之间ORF5的核苷酸相似性达80.3%~100.0%,其中所测的各个分支毒株与其所属分支的代表毒株NADC30、GM2和JXA1的核苷酸相似性分别为92.6%~93.9%、89.8%~93.0%和93.4%~99.3%,检测到新分支中新出现的毒株与JXA1的核苷酸相似性为91.0~93.0%。GP5氨基酸比对结果显示不同亚群GP5信号肽区、诱导表位、主要中和表位及高变区等功能区发生了显著变异。

    结论 

    福建地区出现了PRRSV新分支,使福建PRRSV变异种类更加多样化,建议加强PRRSV流行及变异的监测。

    Abstract:
    Objective 

    To study the molecular prevalence of the porcine peproductive and respiratory syndrome virus (PRRSV) strain in small and medium-sized pig farms of Fujian province in 2017.

    Method 

    We collected 83 suspected blue ear positive disease materials from small and medium-sized pig farms in Fujian province. RT-PCR detection and genetic evolution analysis of ORF5 gene which encoded GP5 protein was performed.

    Result 

    The positive rate of 83 suspected blue ear disease samples was 73%, and the nucleotide similarity among the measured ORF5sequences was 80.3%−100.0%. The nucleotide similarities between tested strains of different branches and the representative strains of those branches, NADC30, GM2 and JXA1, were 92.6%−93.9%, 89.8%−93.0% and 93.4%−99.3%, respectively. The new identified strain in the new branch had nucleotide similarity of 91.0%−93.0% to JXA1. The GP5 amino acid alignment showed significant variation in functional regions such as GP5 signal peptide region, induction epitope, principal neutralizing epitope and hypervariable region in different subgroups.

    Conclusion 

    A new branch of PRRSV has emerged in Fujian area, making the PRRSV variants in Fujian more diverse. It is necessary to strengthen the monitoring of PRRSV epidemics and mutations.

  • 丛枝菌根真菌(Arbuscular mycorrhizal fungi,AMF)是属于球囊霉门Glomeromycota的专性共生微生物,可以与陆地80%以上的植物根系共生。AMF自身生长所需的碳源完全依赖宿主植物,作为回馈,它们可以帮助宿主吸收矿质营养,提高宿主植物对生物及非生物胁迫的抗性[1]。研究发现AMF对植物种群、群落甚至是草地生态系统都有重要的调控作用[2-4]

    高寒草甸是青藏高原东南部主要的草地类型[5],其主要建群植物为莎草科嵩草属的小嵩草Kobresia pygmaea、矮嵩草K. humilis、线叶嵩草K. capillifolia和藏嵩草K. tibetica[6-7],因此也称作高寒嵩草草甸。莎草科植物一直被认为是非丛枝菌根植物,但有研究发现在青藏高原AMF可以与嵩草属植物共生[8-9],并对土壤团聚体的形成起重要作用[10]。基于AMF在草地生态系统中的重要性以及高寒嵩草草甸在青藏高原草地中所占的比例,研究高寒嵩草草甸生态系统中AMF的群落构建及动态对深入理解青藏高原草地土壤−植被相互作用具有重要意义。

    近几十年来,受全球气候变暖和人类活动加剧的影响,青藏高原约1/3的草地发生了不同程度的退化[11]。施肥是退化草地恢复的有效方式之一,目前关于养分添加(主要是氮和磷)能否及如何影响高寒草甸AMF群落的研究相对稀少,并且仅限于海拔不超过3 500 m的区域[12-14],不能代表海拔为3 200~5 200 m的青藏高原高寒嵩草草甸[7]。目前基于传统形态鉴别[15]和分子测序[16]的研究表明,随海拔提升,AMF侵染率、孢子数量和物种丰度下降,AMF群落组成也随之变化。不同海拔高度AMF群落对施肥的响应可能不同,为了更全面地认识高寒草甸,有必要选择海拔更高的区域进行研究。

    本试验选取海拔约4 500 m的高寒嵩草草甸施肥样地,研究不同氮、磷添加对青藏高原高寒嵩草草甸AMF群落的影响及潜在驱动因子。本研究的核心科学问题为探究高寒嵩草草甸根系AMF群落对氮、磷添加的响应机制,及在添加氮、磷条件下主导AMF群落变化的因素。

    研究样地位于中国科学院青藏高原研究所那曲生态环境观测站(E91°12ʹ~93°02ʹ,N30°31ʹ~31°55ʹ),属于西藏那曲东南部,海拔4 480 m,年平均气温−2.1 ℃,年平均降水406.2 mm,高原山地气候,土壤类型为始成土,主要植物种类有小嵩草、矮嵩草和藏薹草Carex thibetica等。

    2013年设置了48个5 m × 5 m的小区,完全随机设计,3个施氮梯度,4个施磷梯度,4次重复。施用的氮肥为尿素,各处理每年施用量分别为0、7.5、15.0 g·m−2,分别记为N0、N1和N2。施用的磷肥为过磷酸钙,有效成分为P2O5,各处理每年P2O5施用量分别为0、7.5、15.0、30.0 g·m−2,分别记为P0、P1、P2和P3。每年7月末施肥1次。

    2016年8月18日,在每个小区随机挑选1个植被覆盖度相似的点作为样点。用土钻从每个样点均取深20 cm、直径8 cm的土块。将土块过2 mm筛,去除土和石块,将保留的根系装入自封袋作为1份样品,−20 ℃保存。另取部分过筛的土装入自封袋,4 ℃保存。共计取得48份根系样品和48份土壤样品,冷藏运送至北京林业大学草地资源与生态实验室。在进行后续试验分析前,将所有根系样品室温解冻,然后用清水洗净。

    取适量土壤样品进行土壤化学成分分析。土壤有机碳用重铬酸钾滴定法[17]测定;土壤全氮用凯氏定氮法[18]测定;土壤铵态氮和硝态氮均用1 mol·L−1的氯化钾溶液浸提,然后用紫外分光光度法[19]测定;土壤全磷和有效磷分别用高氯酸−硫酸和碳酸氢钠浸提,并用钼−锑比色法[20-21]测定;pH在水土质量比为2.5︰1的条件下测量。

    根据Brundrett等[22-23]的方法,先从每份根系样品中剪取若干长1.5 cm的幼嫩根段,用台盼蓝染色。染色后在每份样品中挑出30个根段,切成1 cm根段,然后将每10个根段平行地放置在载玻片上,滴加乳酸甘油制片,每份样品制作3个玻片。在200×显微镜下根据McGonigle等[24]的方法测定AMF侵染率。

    在每份根系样品中,剪取40根长1 cm的幼嫩根段,用CTAB法[25]提取其中的DNA。提取出的DNA用无菌去离子水稀释20倍作为DNA原液,然后进行巢式PCR扩增。2轮PCR的反应体系均为25 μL:2 × Pfu PCR MasterMix (KP201) 12.5 μL,无菌去离子水8.5 μL,DNA模板2 μL以及2个引物各1 μL。第1轮PCR使用的引物为NS31(5ʹ-TTGGAGGGCAAGTCTGGTGCC-3ʹ)和AML2(5ʹ- GAACCCAAACACTTTGGTTTCC-3ʹ),反应条件:94 ℃预热3 min;94 ℃变性45 s,53 ℃退火45 s,72 ℃延伸45 s,35个循环;72 ℃延伸10 min。第1轮PCR的产物稀释至1/10作为第2轮PCR的模板。第2轮PCR使用的引物为AMV4.5NF(5ʹ-AAGCTCGTAGTTGAATTTCG-3ʹ)和AMDGR(5ʹ-CCCAACTATCCCTATTAATCAT-3ʹ)[26],其中AMDGR的5ʹ端带有由12个碱基组成的barcode。第2轮PCR反应条件与第1轮PCR相同,所得的产物用10 g·L−1 2 × TAE的琼脂糖凝胶进行电泳(电压120 V,时间30 min)检测。若样品有目标DNA条带(约280 bp)检出,取4 μL第2轮PCR模板用2倍的反应体系进行第2轮PCR。若样品无目标条带检出,用该样品的DNA原液重复巢式PCR,若2次重复均无目标条带,则去除该样品。根据此方法N1P1处理中的1个样品无条带,去除后剩余47个样品。切取含目标条带的凝胶,使用DNA胶纯化试剂盒(Axygen,美国)进行DNA纯化。纯化后的DNA使用Nanodrop 8000超微量分光光度计测定浓度,根据测定的浓度用移液枪从每个样品中吸取含100 ng DNA的溶液混合至1个离心管,送至中国科学院成都生物研究所在Illumina Miseq平台进行高通量测序。

    测序得到的原始数据使用QIIME[27]进行分析。首先对数据进行质量控制,最大错误期望=0.5,最短序列长度=200,然后依次去除重复、chimeras和singletons。使用UPARSE-OTU算法和Silva数据库,依据97%的相似度阈值划分代表OTU并制作OTU表。在Excel中除去不属于球囊霉门、出现的样本数少于3以及序列数小于总序列量0.01%的OTU,并删除序列数小于总序列1%的样品以降低误差。对所有样品根据最小的样本序列量进行重新抽样,以减小样本量不同造成的差异。此外,将代表OTU序列上传至NCBI,与GenBank的已有序列进行比对,获得相近的参考序列。使用代表序列和参考序列在MEGA 7中使用p-distance模型,在bootstrap值为1 000下构建邻接树,直观查看代表OTU的系统发生情况。本研究获得的代表OTU序列均上传至欧洲核苷酸档案库,序列号为LR535994~LR536029。

    统计OTU数量作为OTU丰度,并计算AMF各科在各样品和各处理中的相对丰度。在R(3.2.2)中使用Vegan包中diversity功能计算每个样品的Shannon多样性指数。

    进行方差分析前,将数据进行转换以符合数据分布的正态性和方差同质性。其中,对土壤化学成分、OTU丰度和Shannon多样性指数数据进行对数转换,侵染率和AMF科相对丰度数据作反正弦平方根转换。转换后的数据在JMP 11中进行双因素方差分析,分析氮、磷添加是否显著影响以上数据,具有显著差异(P<0.05)的结果使用Tukey's HSD检验进行多重比较。

    在群落分析前,先应用R对OTU−样品矩阵进行Hellinger转换以减少稀有OTU的影响,然后用vegdist功能将矩阵转化为Bray-Curtis矩阵。在R中使用Vegan包中的adonis功能进行PerMANOVA分析,研究施氮、磷和其交互作用对AMF群落组成的影响。将土壤化学成分数据作为环境变量导入R并转化为Euclidean矩阵,与之前的群落矩阵一起进行Mantel分析,以研究土壤化学成分与AMF群落组成的关系。

    最后,进行排序分析以研究不同处理AMF的群落分布情况。首先用原始的OTU−样品矩阵进行去趋势对应分析,根据排序轴长结果进一步使用典范对应分析,并用envfit功能进行Monte Carlo检验找出与AMF群落分布有显著相关的环境因子。最终使用ggplot包进行画图,直观显示各处理下AMF群落分布及其和环境因子的关系。

    方差分析得出,不同施氮处理铵态氮和硝态氮含量差异均极显著(P<0.001);不同施磷处理全磷和有效磷含量也有极显著(P<0.001)差异;氮、磷添加的交互作用对硝态氮和有效磷含量有显著(P=0.021,P=0.040)影响(表1)。其中,N2处理铵态氮和硝态氮含量高于N0和N1处理;P3处理全磷和有效磷量高于其余施磷处理,P2处理有效磷含量高于P0处理(表2)。氮、磷添加及其交互作用未对土壤有机碳、全氮和pH产生显著影响(表1)。

    表  1  不同施肥处理土壤化学成分含量差异显著性分析
    Table  1.  Significance analyzes of soil chemical component content differences in different fertilization treatments
    指标 Index 氮添加
    Nitrogen addition
    磷添加
    Phosphorus addition
    氮、磷添加交互作用
    Interaction between nitrogen and phosphorus additions
    F P F P F P
    有机碳 Organic carbon 0.610 0.549 0.940 0.431 1.849 0.117
    全氮 Total nitrogen 0.737 0.486 0.343 0.794 0.674 0.672
    铵态氮 Ammonium nitrogen 15.662 < 0.001 1.094 0.367 1.939 0.105
    硝态氮 Nitrate nitrogen 21.780 < 0.001 2.711 0.061 2.964 0.021
    全磷全磷 Total phosphorus 0.304 0.740 17.135 < 0.001 0.389 0.881
    有效磷 Available phosphorus 0.502 0.610 33.642 < 0.001 2.509 0.040
    pH 1.336 0.276 0.329 0.804 0.948 0.473
    下载: 导出CSV 
    | 显示表格
    表  2  不同施肥处理试验样地的土壤化学成分含量1)
    Table  2.  Soil chemical component contents of the experiment field in different fertilization treatments
    处理 Treatment w/% w/(mg·kg−1) pH
    N P 有机碳
    Organic carbon
    全氮
    Total nitrogen
    全磷
    Total phosphorus
    铵态氮
    Ammonium nitrogen
    硝态氮
    Nitrate nitrogen
    有效磷
    Available phosphorus
    N0 P0 4.79±0.45a 0.26±0.01a 0.040±0.001cd 2.53±0.77b 7.49±1.43bc 2.72±0.49d 6.74±0.08a
    P1 3.61±0.79a 0.21±0.06a 0.040±0.006cd 3.91±2.16b 5.23±1.25bc 6.39±2.34cd 6.77±0.11a
    P2 3.25±0.58a 0.23±0.04a 0.044±0.007abcd 2.06±0.27b 6.26±1.90bc 13.64±7.19bcd 6.66±0.06a
    P3 3.69±0.65a 0.27±0.09a 0.057±0.012abc 2.15±1.37b 4.67±1.14c 21.78±4.00abc 6.65±0.12a
    N1 P0 3.59±0.42a 0.22±0.01a 0.038±0.001d 1.83±0.73b 15.32±5.27bc 2.17±1.35d 6.69±0.08a
    P1 3.35±0.40a 0.24±0.01a 0.046±0.003abcd 2.25±1.07b 6.96±0.95bc 9.73±0.83bcd 6.47±0.15a
    P2 3.77±0.47a 0.22±0.04a 0.042±0.004bcd 6.07±4.38b 9.00±0.49bc 5.49±2.61d 6.73±0.19a
    P3 3.83±0.73a 0.24±0.04a 0.060±0.010ab 6.28±5.37b 5.79±0.91bc 26.02±13.21ab 6.67±0.23a
    N2 P0 2.98±0.32a 0.19±0.01a 0.038±0.003d 27.07±11.00a 35.66±15.70ab 2.28±0.12d 6.64±0.11a
    P1 3.57±0.94a 0.23±0.08a 0.044±0.004abcd 13.72±9.81ab 24.76±22.05bc 4.54±1.96d 6.61±0.15a
    P2 3.29±0.39a 0.24±0.05a 0.047±0.003abcd 7.13±0.85b 14.94±3.58bc 11.30±5.63bcd 6.61±0.05a
    P3 4.24±1.50a 0.23±0.03a 0.060±0.010a 15.41±11.90ab 59.55±23.70a 35.21±5.05a 6.60±0.07a
     1)表中数据为平均值±标准差;同列数据后不同小写字母表示差异显著(P < 0.05,Tukey′s HSD法);N0、N1、N2指氮施用量为0、7.5、15.0 g·m −2的处理;P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理
     1) Data in the table were mean value ± standard deviation; Different lowercase letters in the same column indicated significant differences (P<0.05, Tukey′s HSD test); N0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2
    下载: 导出CSV 
    | 显示表格

    图1所示,根系样品AMF侵染率平均值为73.1%。方差分析结果表明,氮、磷的施加及其交互作用对AMF侵染率无显著(P= 0.350,P = 0.119,P = 0.562)影响。

    图  1  不同施肥处理下丛枝菌根真菌侵染率
    N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理;P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理
    Figure  1.  Arbuscular mycorrhizal fungal colonization rates in different fertilization treatments
    N0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2

    测序数据除杂后,剩余43个样品,共获得234 041条AMF序列。每个样品序列数从297至15 397不等。将所有样品重新抽样至样本量为297后,获得36个OTU,每个样品的OTU数量从1至12个不等,平均值为5.6个。根据稀释曲线(图2)可以看出,重新抽样后N0P2、N1P0、N1P3、N2P2和N2P1对应的曲线末端趋于平滑,说明测序深度较为饱和。经比对,所有AMF的OTU分属7个科,其中球囊霉科有27个OTU,10 285条序列,占总序列的80.5%;多孢囊霉科Diversisporaceae有3个OTU,666条序列,占总序列的5.2%;碎球囊霉科Claroideoglomeraceae有2个OTU,146条序列,占总序列的1.1%;无梗囊霉科Acaulosporaceae、原囊霉科Archaeosporaceae、巨孢囊霉科Gigasporaceae和和平囊霉科Pacisporaceae各有1个OTU,序列数分别为441、532、351和350条,共占总序列的13.1%(图3)。根据OTU代表序列和参考序列构建的邻接树如图4所示。

    图  2  不同施肥处理下分子测序的稀释曲线
    N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理;P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理
    Figure  2.  Rarefaction curve of sequencing samples in different fertilization treatments
    N0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2
    图  3  不同施肥处理下丛枝菌根真菌各科的相对丰度
    N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理;P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理
    Figure  3.  Relative abundance of different arbuscular mycorrhizal fungal families in different fertilization treatments
    N0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2
    图  4  36个丛枝菌根真菌OTU代表序列及其参考序列构建的邻接树
    模型:p-distance;Boostrap值:1 000;DQ846895作为outgroup
    Figure  4.  Neighbor-joining tree constructed based on representative sequences of 36 arbuscular mycorrhizal fungal OTUs and their reference sequences
    Model: p-distance;Boostrap value: 1 000; DQ846895 was used as an outgroup

    方差分析结果显示,氮、磷添加及其交互作用对OTU丰度和Shannon多样性指数无显著影响(图5A5B)。但施氮对球囊霉科的相对丰度有显著影响(P<0.001),N2处理球囊霉科的相对丰度显著低于N1(图5C)。

    图  5  不同施肥处理丛枝菌根真菌OTU丰度、Shannon多样性指数及球囊霉科相对丰度
    不同柱子上不同大写字母表示在P<0.001水平差异显著(Tukey’ s HSD检验)
    Figure  5.  OTU richness, Shannon diversity index of arbuscular mycorrhizal fungi and relative abundance of Glomeraceae in different fertilization treatments
    Different capital letters on different columns indicated significant differences at P<0.001 level (Tukey’ s HSD test)

    PerMANOVA结果表明氮、磷添加处理及其交互作用对AMF群落组成无显著影响(P = 0.680,P= 0.473,P= 0.589)。Mantel分析结果显示AMF群落组成与有机碳含量和硝态氮含量有显著正相关(r=0.176, P=0.110;r=0.142, P=0.041)关系。此外,Monte Carlo结果表明有机碳、硝态氮、全磷和有效磷含量均对AMF群落结构有显著影响(r=0.04, P=0.001;r=0.327, P=0.013;r=0.185, P=0.025;r=0.188, P=0.020)。这4个环境因子在典范对应分析中对排序结果的解释量为11.70%,第1轴和第2轴的解释量分别为6.13%和2.61%。(图6)。

    图  6  不同施肥处理丛枝菌根真菌群落以及显著环境变量的典范对应分析图
    Figure  6.  Canonical correspondence analysis plot of arbuscular mycorrhizal fungal community distribution and significant environmental variables among different fertilization treatments

    在农田和天然草地生态系统中,施加氮和磷会降低AMF侵染率、OTU丰度和多样性,改变根系和根际土壤中AMF的群落结构[28-29]。研究认为出现这种现象的原因是土壤养分的增加使植物更多依靠自身根系去吸收营养,降低对帮助其吸收养分的微生物的依赖[30]。本研究发现,施氮显著提高了土壤中铵态氮和硝态氮的含量,施磷显著增加了土壤中全磷和有效磷的含量,但氮、磷的施加对AMF侵染率、OTU丰度、多样性以及群落组成均无显著影响。因此,本研究的结果与此前在其他生态系统下的研究结果存在差异[28-29],该差异可能是由以下原因造成的。

    首先,差异可能是海拔不同造成的。目前相关研究[12-14]的海拔不超过3 500 m,低于本试验样地的海拔4 480 m。海拔会显著影响AMF的群落和功能,Gai等[15]和Shi等[31]通过形态学方法发现不同海拔条件下AMF侵染率、孢子密度以及菌丝密度都会有显著变化;Liu等[16]和Li等[32]通过分子测序手段发现不同海拔条件下AMF群落组成有显著差异。本研究与此前研究的海拔不同,气候、植被及土壤条件都会有不同,进而导致AMF对氮、磷添加的差异性响应。

    第二,在高寒地区AMF不仅能帮助植物进行养分吸收,同时也可帮助植物应对环境胁迫。因此即使在施肥条件下植物不需要AMF来帮助其获取养分,还是要与AMF保持共生关系应对胁迫。Xiang等[13]在海拔3 220 m的青藏高原高寒草地施肥3年,发现同时添加氮和磷时AMF的OTU丰度和多样性均显著高于对照,也高于氮、磷单独添加时的水平。Chen等[33]研究认为,虽然在养分充足的条件下,植物在营养吸收上降低了对AMF的依赖,然而,植物依然需要AMF来帮助其对抗环境胁迫,尤其是球囊霉科可以提高宿主抗寒性。在本研究的高寒草甸生态系统中,一方面氮、磷的添加使得植物群落降低了对AMF吸收养分的依赖,另一方面植物群落却增加了依靠AMF来对抗胁迫的需要。因此,2个效应综合使得在氮、磷添加的情况下AMF的OTU丰度、多样性和群落组成未受显著影响。

    另一个可能引起差异的原因是样品内及样品间的差异较大,影响了氮、磷添加的效应。从稀释曲线可以看出,个别样品的曲线在终点处仍在上升,说明其测序深度未达到饱和,此时增加取样量可能会有更多的OTU出现。同时,通过典范对应分析图可以直观地发现,N2P3处理3次重复的AMF群落分布距离较大,该处理3次重复的AMF群落的差异是所有处理的重复间差异最大的。基于样本内的误差,本试验结果存在一定限制性。

    本研究中,PerMANOVA分析发现,施加氮、磷不影响AMF群落组成,但通过Mantel分析发现AMF群落的组成与有机碳和硝态氮含量有显著的正相关关系。典范对应分析的结果表明有机碳、硝态氮、全磷和有效磷含量与AMF群落的分布有显著的正相关关系,其中有机碳与AMF群落分布的相关性最大。有机碳主要对低磷条件(P0和P1处理)AMF群落有影响,硝态氮、全磷和有机碳主要对高磷条件(P3处理)AMF群落有影响。相关研究也表明,土壤有机质、硝态氮和有效磷均是对高寒草甸AMF群落有显著影响的土壤成分[34-35]

    Zheng等[34]发现施氮能改变AMF群落,但不是通过直接作用,而是间接地通过改变土壤其他成分(如有机碳)和植物群落来影响AMF群落。土壤成分也可以直接或间接地影响微生物群落[36-37]。因此可以解释本研究中施肥未对AMF群落产生影响,但特定的土壤因子与AMF群落分布有显著相关性。

    综上所述,基于青藏高原4 500 m海拔高寒草甸的研究发现,氮、磷添加对AMF的侵染率、OTU丰度和多样性无显著影响。该发现回答了本研究要明确的第1个核心科学问题,即青藏高原高寒嵩草草甸根系中的AMF群落不受氮、磷添加的影响,与海拔较低的AMF群落对氮、磷的响应不同;对第2个核心科学问题,影响青藏高原高寒嵩草草甸根系AMF群落的因子主要是有机碳和硝态氮含量,全磷和有效磷含量是次要因子。未来研究还需要系统选取不同海拔的样地,以获取青藏高原不同海拔高寒草甸AMF群落对施肥的响应。另外,在全面剖析AMF群落变化规律及驱动机制的基础上,亟需研究AMF群落变化如何反馈影响地上植被的个体生长与群落动态,探索高寒草地退化及恢复演替中的地下生态过程及机理。

  • 图  1   PRRSV GP5氨基酸序列比对分析结果

    PS:信号肽;D:诱导表位;PNE:主要中和表位;HVR1、HVR2:高变区;T1、T2:T细胞表位;B:B细胞表位;TM1、TM2、TM3:跨膜区域;黄色高亮部分表示糖基化位点

    Figure  1.   Comparison analysis of amino acid locus variation in PRRSV GP5

    PS: Peptide signal; D: Decoy epitope; PNE: Principal neutralizing epitope; HVR1, HVR2: Hypervariable region; T1, T2: T cell epitope; B: B cell epitope; TM1, TM2, TM3: Transmembrane regions; Yellow highlighted region refers to glycosylation site

    图  2   ORF5基因遗传进化分析

    采用NJ法建树,自举值检验为1 000

    Figure  2.   Genetic evolution analysis of the ORF5 gene

    The tree is built using the NJ method, and the bootstrap value test is 1 000

    表  1   引物和扩增片段长度

    Table  1   Primer and amplified fragment length

    引物名称
    Primer name
    引物序列(5'→3')
    Primer sequence
    长度/bp
    Length
    ORF5-1-F TGAGACCATGAGGTGGGC 726
    ORF5-1-R GAAAACGCCAAAAGCACC
    ORF5-2-F TGCTCCATTTCATGACAC 974
    ORF5-2-R GCATCTGGAGGTGATGAAT
    下载: 导出CSV
  • [1]

    SNIJDER E J, KIKKERT M, FANG Y. Arterivirus molecular biology and pathogenesis[J]. J Gen Virol, 2013, 94(10): 2141-2163.

    [2]

    FENG Y, ZHAO T, NGUYEN T, et al. Porcine respiratory and reproductive syndrome virus variants, Vietnam and China, 2007[J]. Emerg Infect Dis, 2008, 14(11): 1774-1776. doi: 10.3201/eid1411.071676

    [3]

    TIAN K, YU X, ZHAO T, et al. Emergence of fatal PRRSV variants: Unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark[J]. PLoS One, 2007, 2(6): e526. doi: 10.1371/journal.pone.0000526

    [4]

    TONG G Z, ZHOU Y J, HAO X F, et al. Highly pathogenic porcine reproductive and respiratory syndrome, China[J]. Emerg Infect Dis, 2007, 13(9): 1434-1436. doi: 10.3201/eid1309.070399

    [5]

    ZHOU Y J, HAO X F, TIAN Z J, et al. Highly virulent porcine reproductive and respiratory syndrome virus emerged in China[J]. Transbound Emerg Dis, 2008, 55(3/4): 152-164.

    [6] 宁宜宝, 郑杰, 张纯萍, 等. 我国南方猪高热病的研究: Ⅱ: 猪繁殖与呼吸障碍综合征病毒的分离、鉴定和致病性测定[J]. 中国兽药杂志, 2007(1): 14-18. doi: 10.3969/j.issn.1002-1280.2007.01.005
    [7]

    LUNNEY J K, FANG Y, LADINIG A, et al. Porcine reproductive and respiratory syndrome virus (PRRSV): Pathogenesis and interaction with the immune system[J]. Annu Rev Anim Biosci, 2016, 4: 129-154. doi: 10.1146/annurev-animal-022114-111025

    [8] 邓蓉, 郑文堂, 肖红波. 中国生猪生产区域布局与经营规模分析[J]. 现代化农业, 2016(10): 44-46. doi: 10.3969/j.issn.1001-0254.2016.10.024
    [9]

    ZHAO K, YE C, CHANG X, et al. Importation and recombination are responsible for the latest emergence of highly pathogenic porcine reproductive and respiratory syndrome virus in China[J]. J Virol, 2015, 89(20): 10712-10716. doi: 10.1128/JVI.01446-15

    [10]

    ZHANG Q, JIANG P, SONG Z, et al. Pathogenicity and antigenicity of a novel NADC30-like strain of porcine reproductive and respiratory syndrome virus emerged in China[J]. Vet Microbiol, 2016, 197: 93-101. doi: 10.1016/j.vetmic.2016.11.010

    [11]

    LIU J K, WEI C H, DAI A L, et al. Complete genomic characterization of two European-genotype porcine reproductive and respiratory syndrome virus isolates in Fujian province of China[J]. Arch Virol, 2017, 162(3): 823-833. doi: 10.1007/s00705-016-3136-9

    [12]

    BAI X, WANG Y, XU X, et al. Commercial vaccines provide limited protection to NADC30-like PRRSV infection[J]. Vaccine, 2016, 34(46): 5540-5545. doi: 10.1016/j.vaccine.2016.09.048

    [13] 张金辉. 生猪禁养政策持续加码, 养猪企业何去何从?[J]. 猪业科学, 2016(11): 142-143. doi: 10.3969/j.issn.1673-5358.2016.11.051
    [14] 王亚辉, 彭华, 李娟. 我国21个地区公布畜禽禁养时间表2016年环保政策将更严[J]. 中国猪业, 2016(5): 13-18. doi: 10.3969/j.issn.1673-4645.2016.05.002
    [15] 刘岳. 猪繁殖与呼吸综合征病毒GP5蛋白分段表达及鉴定[D]. 泰安: 山东农业大学, 2010.
    [16]

    XIE J, CUI T, CUI J, et al. Epidemiological and evolutionary characteristics of the PRRSV in Southern China from 2010 to 2013[J]. Microb Pathog, 2014, 75: 7-15. doi: 10.1016/j.micpath.2014.08.001

    [17]

    SHI M, LAM T T, HON C C, et al. Molecular epidemiology of PRRSV: A phylogenetic perspective[J]. Virus Res, 2010, 154(1/2): 7-17. doi: 10.1016/j.virusres.2010.08.014

    [18]

    LU W H, TUN H M, SUN B L, et al. Re-emerging of porcine respiratory and reproductive syndrome virus (lineage 3) and increased pathogenicity after genomic recombination with vaccine variant[J]. Vet Microbiol, 2015, 175(2/3/4): 332-340. doi: 10.1016/j.vetmic.2014.11.016

    [19] 申思. HP-PRRSV囊膜糖蛋白5糖基化位点突变及功能初探[D]. 泰安: 山东农业大学, 2016.
  • 期刊类型引用(3)

    1. 杨磊,袁斌,郑锷,叶晨朔,王思行,何贤俊,张潇潇,黄山,胡伟伟,邵敏. 珠江三角洲秋季生物质燃烧对有机气溶胶的贡献. 中国环境科学. 2023(01): 20-28 . 百度学术
    2. 黄巧义,于俊红,黄建凤,黄旭,李苹,付弘婷,唐拴虎,刘一锋,徐培智. 广东省主要农作物秸秆养分资源量及替代化肥潜力. 生态环境学报. 2022(02): 297-306 . 百度学术
    3. 郭畅,刘剋. 基于MODIS数据京津冀的秸秆焚烧火点监测与分析. 科技风. 2019(31): 125-126 . 百度学术

    其他类型引用(4)

图(2)  /  表(1)
计量
  • 文章访问数:  1733
  • HTML全文浏览量:  0
  • PDF下载量:  2002
  • 被引次数: 7
出版历程
  • 收稿日期:  2018-10-17
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2019-05-09

目录

/

返回文章
返回