作者简介:董波(1983—),男,天津人,博士,讲师,主要从事动物病毒学研究。E-mail: 381289930@qq.com
为了解2015—2017年闽西地区猪流行性腹泻病毒(PEDV)遗传进化规律,本研究收集了闽西地区12份PEDV阳性样品,对其M基因进行RT-PCR扩增,克隆至pMD18-T载体测序,并与国内外已知参考毒株序列进行比对及遗传进化分析。结果表明,12株闽西毒株M基因氨基酸序列同源性为100.0%;遗传进化分析表明,12株闽西毒株与经典毒株CV777、2010年前中国分离毒株LCZ、以及韩国毒株SM98、virulent DR13的亲缘性较远,而与2010年后国内分离毒株以及2013年美国分离毒株亲缘关系最近,说明闽西PEDV毒株属于目前国内外流行毒株。分子结构特征分析表明,PEDV M基因具备高度保守的特性,可以作为设计疫苗的候选基因。
To elucidate evolutionary characteristics of the memberane (M) genes of porcine epidemic diarrhea viruses in Fujian, China, 12 samples collected from Fujian province during 2015 to 2017 were screened by multiplex RT-PCR,M genes of PEDV viruses from positive samples were amplified by RT-PCR and cloned into pMD18-T vectors for sequencing respectively. The sequence data were done alignment and phylogenetic analysis with sequences from known reference strains. The M genes from 12 representative positive samples shared 100.0% identity at amino acid level. 12 representative positive samples are remote to classical strain CV777, SM98 and virulent DR13. 12 representative strains belong to the genotype variant strain. Molecular structure characteristics analysis showed that M gene is highly conservative and could be an optimal candidate protein for vaccine development.
猪流行性腹泻(porcine epidemic diarrhea, PED)是由猪流行性腹泻病毒(porcine epidemic diarrhea virus, PEDV)引起猪的一种急性、高度接触性肠道传染病, 以水样腹泻、呕吐和脱水为临床特征, 主要感染哺乳期仔猪, 感染率和病死率极高, 给养猪业造成极大损失[1, 2, 3]。该病第一次被报道于1971年的英国, 之后在欧洲、亚洲相继有该病的出现[4, 5]。我国第一次出现该病是在1976年, 继而呈现大面积流行趋势[6]。自2010年冬季以来, 我国10多个省市暴发PED, 对国内养猪业带来了严重影响, 甚至出现以前接种过PED疫苗的猪也未能幸免, 提示该流行毒株毒力较强, 而且与现有疫苗毒株存在较大的差异, 导致现有疫苗难以提供免疫保护[7, 8]。2011年, 张世忠等[9]对福建省各个地区的规模化养猪的顽固性腹泻病例进行检测, PEDV感染率高达52.2%。由此可见, 该病一旦暴发流行, 其危害十分严重, 对养猪业的潜在危害不容忽视。2013年, 美国多个洲先后有暴发PEDV, 病猪表现出水样腹泻、急性肠炎、呕吐等典型症状, 死亡率极高, 给美国养猪业造成重创, 引起了美国乃至全球养猪业的高度重视[10, 11]。
PEDV为单股正链RNA病毒, 全长大约28 kb, 由4种结构蛋白组成, 分别为纤突蛋白(spike, S)、膜蛋白(membrane, M)、核衣壳蛋白(nucleocapsid, N)和小膜蛋白(small membrane, E)[12]。其中, M基因全长681 bp, 编码226个氨基酸, 位于病毒表面, 是病毒囊膜中含量最多的基因[13]。M基因在PEDV的结构基因中, 相对较为保守, 然而, 近些年的相关文献显示, M基因也呈现了进化和变异的趋势。闽西地区地处闽粤赣三省交界, 一旦出现病毒变异并迅速蔓延, 对三省乃至全国的养猪业将造成不小的影响。因此, 本研究收集了2015年1月— 2017年6月闽西地区规模化猪场送检PEDV阳性病例, 采用RT-PCR技术, 扩增阳性样品全长M基因, 旨在从分子水平了解PEDV闽西地区流行毒株的特点, 进一步阐明闽西地区流行毒株的遗传变异规律, 为我国PEDV地方毒株的遗传演化分析提供基础资料, 并为下一步制备PEDV疫苗提供参考依据。
PEDV弱毒疫苗购自福州大北农公司; 野外阳性样品为2015年1月— 2017年6月年期间采自闽西地区规模化养猪场送检腹泻症状仔猪的粪便和小肠组织样品12份; 大肠埃希菌DH5α 感受态细胞, 由本实验室制备、保存; pMD18-T Vector购自宝生物工程(大连)有限公司。
DL2000 Marker、总RNA提取液(TaKaRa RNAiso Reagent)、病毒RNA/DNA提取试剂盒(TaKaRa MiniBEST Viral DNA/RNA Extraction Kit Ver.3.0)、AMV反转录酶、PrimeSTAR® HS DNA Polymerase 高保真酶、Oligo(dT)18、dNTP、DNA Marker DL 2000、RNA酶抑制剂购自宝生物工程(大连)有限公司; DNA凝胶回收试剂盒(DNA Gel Extraction Kit)及质粒提取试剂盒(Plasmid Miniprep Kit)购自AXYGEN公司产品; 感受态细胞制备试剂购自宝生物工程(大连)有限公司; 氯仿、异丙醇采用分析纯试剂。
按照试剂盒提供的操作手册, 自粪便和小肠组织以及灭活的PEDV样品(阳性对照)中提取总RNA。取洁净EP管, 加入RNA 11 μ L、Oligo(dT)18 1 μ L, 12 000 g离心40 s, 于70 ℃水浴10 min, 加入5× M-MLV Buffer 5 μ L、M-MLV 1 μ L、dNTP 1 μ L、DEPC水 6 μ L, 混匀后, 42 ℃水浴30 min, 95 ℃水浴10 min。cDNA于-80 ℃下保存。
参考GenBank上发表的PEDV CV777株M基因序列, 使用Premier Premier 5.0生物软件设计合成一对特异性引物, PEDV-M-F: TTCTATTCCCGTTGATGAGG; PEDV-M-R: CATGAAGCACTTTCTCACTATC。以cDNA为模板, 扩增病毒M基因。反应条件为:95 ℃预变性5 min; 94 ℃变性30 s, 56 ℃退火30 s, 72 ℃延伸45 s, 35个循环; 最后72 ℃延伸10 min。产物经电泳分离, 凝胶切割纯化, 凝胶成像系统下观察结果。
取5 μ L纯化产物与0.5 μ L pMD18-T载体(50 ng· μ L-1)进行连接反应, 16 ℃连接4 h后, 将连接产物转化入DH5α 感受态细胞, 涂板过夜, 挑取单个白色菌落, 每个样品挑选3个菌落, 以保证测序准确性。过夜培养, 提取质粒, 经鉴定正确后, 送上海生工生物工程技术服务有限公司测序。
使用DNAman、DNAstar、MEGA5.2分子生物学软件, 对M基因序列进行同源性及系统发育分析。
12份阳性样品编号为: FJLY1501、FJLY1502、FJLY1503、FJLY1504、FJLY1601、FJLY1602、FJLY16031、FJLY1604、FJLY1701、FJLY1702、FJLY1703、FJLY1704, 对M基因进行PCR扩增。M基因PCR产物经扩增后, 使用凝胶电泳进行鉴定, 结果显示获得大小约为681 bp基因片段(图1)。将纯化后的PCR产物连接到pMD18-T载体, 将阳性克隆送往上海生工生物工程技术服务有限公司测序, 序列结果正确。
将12份阳性样品的氨基酸序列与国内外已知代表毒株M基因序列进行比对及同源性分析。所引用的已知代表毒株信息见表1。12份阳性样品M基因与其他代表性毒株M基因进行氨基酸同源性比较后发现, 闽西毒株M基因之间的同源性为100.0%, 与经典毒株CV777和LCZ的同源性分别为98.7%和98.6%, 与2014年福建分离毒株FJ-ZP和FJ-FQ的同源性分别为99.6%和99.7%。与2010年以后中国分离毒株的同源性均在99%以上, 与2013年美国毒株的同源性也均在99%以上。
遗传进化分析表明, 42株PEDV可以分成两个大组。其中, Ⅰ 组包括经典毒株CV777, 3株韩国毒株, 6株中国毒株; Ⅱ 组则包括1株韩国毒株, 4株美国毒株, 22株2010年后中国毒株以及12株闽西分离毒株。系统发育分析显示, 闽西分离毒株与经典毒株CV777、2010年前中国分离毒株LCZ、SM98以及韩国毒株virulent DR13的亲缘性较远, 与2010年以后中国分离毒株的亲缘性较高, 与2013年美国毒株亲缘性也较高。见图2。
PEDV是PED的致病原, 主要感染对象为哺乳期仔猪, 致死率高, 对养猪业的危害非常严重。研究证实, PEDV只存在一个血清型, 然而, 分子流行病学调查显示, 近十年来我国流行毒株与欧洲流行毒株、原始毒株比较, 我国流行毒株的结构蛋白基因(M、N、S)均与欧洲流行毒株及原始毒株亲缘关系较远, 可能出现了新的基因型[14]。2011年, 刘孝珍等[15]将17株分离得到的PEDV S基因进行遗传进化分析研究, 证实我国出现了PEDV新的基因型。M基因能诱导机体产生针对病毒的中和抗体, 在病毒粒子的组装和出芽过程中具有重要作用[16]。M基因是高度保守的基因, 能介导机体产生干扰素, 常用作RT-PCR诊断靶基因和疫苗候选基因[17]。在遗传进化的分析中, 选择M基因, 有助于更多地了解病毒的流行趋势与进化规律, 对于病毒的防控和疫苗的研制具有重要意义。
2011年冬末, 广东、福建等地区持续暴发仔猪腹泻, 引起了学者的关注, 学者们普遍认为, 造成这次PEDV暴发的毒株为变异毒株。自此以后, 该变异毒株成为我国主要流行毒株。2013年, 张志等[18]将分离到的PEDV阳性毒株M基因进行遗传进化分析发现, 新分离毒株M基因与2010年后我国新分离毒株亲缘性较高, 为国内变异毒株。本研究中, 闽西毒株M基因氨基酸序列与2010年以后中国分离毒株以及2013年美国毒株的同源性均在99%以上, 证实了闽西流行毒株与国内其他地区流行毒株关系密切, 属于2010年后国内分离变异毒株。
2013年, 美国暴发了猪流行性腹泻, 给美国养猪业造成了严重影响, 经济损失巨大。这次暴发的PEDV, 与2010年我国暴发的PEDV的流行方式极为相似, 并且学者通过进化分析证实, 2013年在美国暴发的PEDV与2010年在中国流行的PEDV, 在亲缘性上很相近, 可以认为来自同一祖先。本研究中, 42株病毒可以分成2个大组, Ⅰ 组和Ⅱ 组。其中, Ⅰ 组主要包括CV777、LZC和SM98为代表的的经典株, 以及韩国DR13弱毒株和国内早期毒株; 而Ⅱ 组包含毒株均为2010年后中国分离毒株以及2013年后美国分离毒株, 这也证实了闽西毒株与2010年后国内分离毒株和2013年后美国分离毒株进化关系更为紧密, 为当前国内外流行毒株。同时, 亲缘比较结果表明, 闽西毒株与疫苗研发株CV777亲缘关系较远, 提示我们现在猪场所使用的疫苗, 对于防治PED的流行所起到的作用并不显著, 这也成为控制PEDV暴发的主要隐患。
M基因氨基酸结构特征比较显示, 与经典毒株CV777相比, 闽西毒株无特有的片段缺失和插入, 12份毒株存在3处相同突变, 分别为12E-Q、42A-V、213A-S, 这些突变位置与2010年后分离毒株一致, 证实了闽西毒株属于2010年后国内分离变异毒株。闽西毒株N-糖基化位点预测结果与经典毒株CV777一致, 且穿膜肽、亲水性、pI值等方面与经典毒株CV777也无差异, 这些结果进一步证明了PEDV M基因具备高度保守的特性, 可以作为设计疫苗的候选基因。
综上, 通过比较分析2015— 2017年闽西地区PEDV M基因氨基酸序列发现, 闽西毒株M基因与2010年国内分离毒株以及2013年美国分离毒株同源性较高, 属于国内外流行毒株。并且, 蛋白预测结果表明, 虽然闽西毒株M基因与经典毒株CV777相比存在3处基因改变, 但是大部分性质没有明显变化, 证明M基因是高度保守的。这些结果的获得, 有助于从分子水平了解闽西地区PEDV M基因分子结构及遗传变异特征, 为进一步研制有效的生物制剂和诊断方法提供资料支持, 同时为猪流行性腹泻的防控提供病原学和流行病学参考依据。
The authors have declared that no competing interests exist.