豬繁殖與呼吸系統綜合征病毒華南株GDgz的分離鑒定及遺傳進化分析

于林洋，董建國，2，張樂宜，劉燕玲，梁鵬帥，王 磊，宋長緒

（1.華南農業大學動物科學學院，廣東 廣州 510642；2.信陽農林學院牧醫工程學院，河南 信陽 464000）

豬繁殖與呼吸系統綜合征病毒華南株GDgz的分離鑒定及遺傳進化分析

于林洋1，董建國1，2，張樂宜1，劉燕玲1，梁鵬帥1，王 磊1，宋長緒1

（1.華南農業大學動物科學學院，廣東 廣州 510642；2.信陽農林學院牧醫工程學院，河南 信陽 464000）

為了解華南地區豬繁殖與呼吸綜合征病毒（PRRSV）流行株的變異特性，用Marc-145細胞從廣東省疑似豬繁殖與呼吸綜合征發病豬場采集的肺組織中分離得到1株PRRSV（命名為GDgz），該毒株能產生明顯的細胞病變。全基因組進化分析和同源性比對分析結果顯示，GDgz與中國高致病性毒株位于同一分支，與歐洲型毒株Lelystad-virus同源性最低、僅為60.3%，與美洲型經典毒株VR-2332的同源性為88.8%，與中國高致病性毒株JXA1和HuN4同源性為98.8% 和98.9%，與JXA1的疫苗株同源性為91.1%，與NADC30和NADC30-like毒株 CHsx1401和JL580同源性分別為84.5%、87.2%和83.6%。NSP2序列比對結果顯示，GDgz在高變區存在30個不連續氨基酸缺失。GP5序列比對結果顯示，GDgz在GP5抗原表位上有氨基酸突變，表明GDgz屬于美洲型高致病性毒株傳代致弱的疫苗株，且在抗原表位上有一定程度的突變。

豬；繁殖與呼吸綜合征病毒；分離鑒定；遺傳進化

豬繁殖與呼吸系統綜合征（P o r c i n e reproductive and respiratory syndrome，PRRS）是造成當今世界養豬業重大經濟損失的傳染性疾病之一，它主要造成母豬的繁殖系統障礙和各年齡豬群的呼吸系統疾病［1-2］。該病的病原是豬繁殖與呼吸系統綜合征病毒（PRRSV）。PRRSV為有囊膜的單股正鏈RNA病毒，屬于套式病毒目、動脈炎病毒科、動脈炎病毒屬［3］。

PRRSV基因組全長約15 kb，含有10個開放閱讀框（ORFs）、5′端UTR、3′端UTR和poly（A）尾巴。開放閱讀框分別為ORF1a、ORF1b、ORF2a、ORF2b、ORF3～7和ORF5a［4-5］。其中ORF1a和ORF1b進行12個非結構蛋白Nsp1～Nsp12的編碼，在病毒復制、轉錄和拮抗干擾素等方面起著重要作用［6-7］。ORF2～7進行結構蛋白GP2、GP3、GP4、GP5、M和N的編碼［8］。病毒中和抗體生成有關的抗原表位主要位于結構蛋白上［9］。5′端UTR用于調節病毒RNA功能及結構蛋白的表達，3′端UTR位于ORF7下游，主要進行病毒復制酶識別并結合以啟動負鏈RNA的合成［10-11］。根據PRRSV變異以及流行情況，PRRSV可分為以Lelystad-virus為代表的歐洲型（Ⅰ型）和以VR2332為代表的北美型（Ⅱ型）［12］。我國PRRSV主要以基因Ⅱ型為主。PRRS對我國養豬業造成巨大的經濟損失，尤其是1996年以低致病性毒株CH-1a為代表和2006年以高致病性毒株JXA1為代表的兩次PRRS暴發對我國養豬業造成了重大經濟損失［13］。2014年，有研究報道NADC30-like毒株在我國廣泛流行，進一步打擊了我國養豬業，也給PRRS的防控帶來了新挑戰［14-17］。

為了進一步了解我國華南地區PRRSV的流行及變異情況，我們于2016年10月對廣東省廣州市某種豬場暴發疾病死亡的30日齡保育小豬進行檢測，結果呈PRRSV陽性，最后通過在Marc-145細胞上將病毒成功分離并對其全基因組進行測定和分析，將測序結果與其他代表PRRSV毒株進行序列比對和進化樹分析，進而分析所分離的毒株與目前國內PRRSV流行毒株的相似性，從而了解廣東省目前PRRSV毒株的流行情況。

1 材料與方法

1.1 試驗材料

1.1.1 病料 廣東省廣州市疑似暴發PRRS的自然發病豬，采集發病豬肺臟組織于-80℃凍存。

1.1.2 細胞、菌體及載體 本試驗所使用的Marc-145細胞、克隆宿主菌E. coliDH5α感受態細胞以及陽性對照病毒株GDzj株均由華南農業大學動物科學學院豬病防控研究室保存。pMD-19T克隆載體購自寶生物工程（大連）有限公司。

1.1.3 主要試劑 病毒DNA/RNA提取試劑盒購自Magen公司；反轉錄試劑盒、膠回收試劑盒購于Promega公司；RT-PCR相關試劑盒購自TOYOBO生物科技有限公司和Vazyme公司；DNA Marker DL2000等購自TaKaRa公司； EB替代染料購自廣州華奇盛生物科技有限公司；pEASY-Blunt Simple Cloning Vector和Trans-5α感受態細胞均購自北京全式金生物技術有限公司。

1.2 試驗方法

1.2.1 引物設計與合成 根據GenBank中發布的VR-2332等PRRSV毒株全基因組序列進行生物信息學分析，選取保守區域設計13對特異性引物（表1）用于擴增全基因組序列。引物由蘇州金維智生物公司合成。

表1 PRRSV全基因組擴增引物

1.2.2 病料處理及PCR檢測 將采集的發病豬肺臟組織剪碎，加入1 mL PBS研磨后收集研磨液于1.5 mL EP管中，12 000 r/min離心5 min，取上清用0.22 μm微孔濾器過濾除菌，收集濾液并取200 μL進行病毒DNA/RNA核酸提取及PRRSV NSP2、CSFV、PPV的 RT-PCR檢測與PRV、PCV2的常規PCR檢測，剩余濾液于-80℃保存。

1.2.3 病毒分離 使用12孔細胞培養板培養Marc-145細胞，待Marc-145細胞長滿單層后棄去上清，用PBS洗3遍，接種50 μL過濾過的含有病毒的病料上清，加入100 μL不含血清、含有2%雙抗的1 640培養液。37℃下孵育1 h后，棄去上清，加入含2%血清的1 640維持液1 mL，置于細胞培養箱中培養3～5 d，約80%細胞出現細胞病變效應CPE后即可收取樣品。將毒液置于-80℃冰箱保存。

1.2.4 病毒接毒 使用25 mL培養瓶培養Marc-145細胞，待Marc-145細胞長滿單層后棄去上清，用PBS洗3遍，接種200 μL過濾的含有病毒的病料上清，加入1 mL不含血清的含有2%雙抗的1 640培養液。37℃下孵育1 h，棄去上清，加入含2%血清的1 640維持液5 mL，置于細胞培養箱中培養3～5 d，約80%細胞出現CPE后即可收取樣品。將毒液置于-80℃冰箱保存。

1.2.5 核酸提取及反轉錄 取在Marc-145細胞上分離的第3代GDgz毒株，根據Magen公司產品試劑盒進行病毒核酸提取，然后根據Promega反轉錄試劑盒對所提取的病毒核酸進行反轉錄。

1.2.6 PRRSV GDgz毒株目的片段分段擴增 使用TOYOBO生物科技有限公司高保真酶對目的片段進行擴增。擴增體系為50 μL，KOD-Plus-Neo 1 μL，10×PCR Buffer for KODPlus-Neo 5 μL，2 mmol/L dNTP 5 μL，25 mmol/L MgSO43 μL，cDNA模板4 μL，上游引物2 μL，下游引物2 μL，RNase Free H2O 28 μL。PCR擴增程序：94℃預變性2 min；98℃變性10 s，62℃退火30 s，68℃延伸30 s/kb，35個循環；72℃延伸10 min。取5 μL的PCR產物進行1%瓊脂糖凝膠電泳。

1.2.7 序列測定 PCR擴增產物通過1%瓊脂糖凝膠電泳鑒定后，使用膠回收試劑盒對目的片段進行膠回收。然后連接到北京全式金生物有限公司的pEASY-Blunt Simple Cloning Vector上。經過轉化挑菌鑒定后送往測序公司進行測序。

1.2.8 序列比對 利用Lasergene及Mega 5.0軟件將分離的毒株全基因組與國內外分離毒株進行序列比對并構建進化樹，并分別將NSP2基因序列和GP5氨基酸序列與其他國內外代表毒株的NSP2基因序列和GP5氨基酸序列進行比對分析。PRRSV代表毒株的基因序列均來源于GenBank。

2 結果與分析

2.1 病料的RT-PCR檢測

用NSP2基因檢測引物對處理好的病料樣品進行檢測，擴增出NSP2基因片段，瓊脂糖凝膠電泳分析結果表明，與陽性對照片段大小一致（圖1），而樣品中RNA病毒CSFV、PPV與DNA病毒PRV、PCV2等檢測均呈陰性（圖略）。

圖1 PRRSV NSP2基因片段擴增電泳結果

2.2 病毒分離

將經處理的病料接種到Marc-145細胞，培養48～96 h，每天觀察細胞情況，待細胞聚集、皺縮并脫落等有明顯CPE時，收取病毒液（圖2，封三）。

2.3 PRRSV全基因組測定

運用13對相互重疊的特異性引物，通過RT-PCR法擴增病毒的13個cDNA片段。瓊脂糖凝膠電泳分析顯示，擴增的13 個RT-PCR 產物大小與預期長度一致（圖3）。將擴增的片段膠回收后連接pEASY-Blunt Simple Cloning Vector上送公司測序，最后對13個片段測序結果進行分析，SeqMan軟件拼接序列，獲得全長的PRRSV基因組序列，全長為15 353 nt，命名為GDgz。

2.4 PRRSV全基因進化樹的構建

圖3 PRRSV基因片段的RT-PCR擴增結果

運用分子生物學軟件對GDgz毒株和其他國內外具有代表性的毒株進行進化樹構建，結果（圖4）顯示，整個遺傳進化樹可分為以Lelystad virus為代表的歐洲型和以VR2332為代表的美洲型兩大分支，以VR-2332為代表的美洲型分支又可進一步分為以美國NADC30株為代表的分支、以VR2332和我國低致病性毒株CH-1a為代表的分支以及以高致病性毒株JXA1、HUN4為代表的分支。結果表明，GDgz與國內的高致病性毒株JXA1、HUN4、JXwn06屬于同一分支，且與高致病性毒株的疫苗株JXA1-P80同源性最高。

圖4 基于PRRSV全基因組的進化樹分析

2.5 PRRSV全基因同源性分析

圖5 PRRSV全基因序列比對結果

用DNAStar軟件進行PRRSV全基因的進化分析，結果（圖5）顯示，GDgz與JXA1-P80同源性最高、為99.1%，與HP-PRRSV毒株JXA1、HUN4和JXwn06同源性分別為98.8%、 98.9%和98.9%，與NADC30和NADC30-like毒株 JL580、CHsx1401同源性分別為84.5%、87.2%和83.6%，與美洲型代表株VR-2332同源性為88.8%，與歐洲型代表株Lelystad virus同源性僅為60.2%。

2.6 PRRSV NSP2氨基酸序列比對

使用DANStar對GDgz和國內外代表毒株進行NSP2氨基酸的序列比對，結果（圖6）顯示，GDgz與高致病性毒株HUN4、JXA1、JXwn06和高致病疫苗株JXA1-P80一樣，NSP2在511、534～562位存在30個不連續氨基酸的缺失，在其他位點沒有氨基酸的缺失或插入。

圖6 PRRSV NSP2氨基酸序列比對結果

2.7 PRRSV GP5氨基酸序列比對

運用DNAStar軟件對GDgz和國內外參考毒株進行GP5蛋白的氨基酸序列比對，已有研究表明美洲型毒株GP5蛋白的表位有3個，2個（aa27～aa30 和 aa180～aa197）為非中和表位，1個（aa37～aa45）為中和表位（圖7），與美洲型代表毒株VR-2332相比，GDgz在非中和表位中有1個A29V的變異，與其他代表毒株相似；在185位有V185A突變，該突變與高致病性毒株及其疫苗株相似；在189位有I189L突變，該突變與國內經典毒株和高致病性毒株及其疫苗株相似。中和表位中，在39位有1個L39I突變，該突變與高致病性毒株及其疫苗株相似。

圖7 PRRSV GP5 氨基酸序列比對結果

3 討論

自從美國分離了第1株PRRSV毒株VR-2332以及1996年中國分離了第1株PRRSV CH-1a以來，PRRS已經成為嚴重影響養豬業發展的疾病之一，每年導致巨大的經濟損失。我國分離的毒株以北美型PRRSV毒株為主［18-19］，目前豬場主要通過疫苗免疫該病，國內用于PRRSV防控的疫苗主要有經典毒株和高致病性毒株減毒活疫苗以及滅活苗。由于我國PRRSV毒株種類眾多，研發生產疫苗的廠家眾多，致使疫苗使用混亂，一個豬場可能存在免疫多種疫苗的情況，導致不同毒株間發生遺傳重組的可能。加上RNA病毒復制酶低保真性和免疫系統壓力的影響，進而加速PRRSV病毒不斷發生變異［20-21］。2006年出現的高致病性豬繁殖與呼吸綜合征病毒和2014年出現的NADC30-like毒株說明我國PRRSV變異情況越來越嚴重，致使現有疫苗免疫效果差，不能很好地提供免疫保護，給疾病防控造成了巨大困難。本研究通過細胞培養分離病毒和RT-PCR方法對GDgz毒株全基因組進行分段擴增，并克隆、測序，通過相關軟件對序列進行拼接后得到GDgz毒株的全基因序列，并對其與其他PRRSV毒株同源性進行比較和構建進化樹。結果顯示GDgz毒株與中國高致病性毒株JXA1疫苗株JXA1-P80同源性最高，表明廣東豬場由于廣泛使用減毒活疫苗，導致疫苗株在豬場廣泛流行，給PRRSV的變異增加了風險。

NSP2基因是PRRSV基因組上變異最大的基因之一，遺傳進化分析表明，NSP2區域存在廣泛的基因缺失、插入和突變。有研究證實HPPRRSV毒株與經典毒株相比在NSP2區域存在30個不連續氨基酸的缺失［13］。對比國內外NSP2氨基酸序列可知，GDgz的NSP2存在30個不連續氨基酸的缺失，符合高致病性PRRSV毒株的特征。

抗原表位是決定病毒免疫特性的重要因素之一。目前，已有研究者證實美洲型 PRRSV GP5 蛋白抗原表位有1個中和表位和2個非中和表位［22］。GP5 蛋白由PRRSV ORF5 基因編碼，是病毒主要結構蛋白之一，為囊膜糖蛋白，具有較高的免疫原性和中和活性［23］。與美洲型代表毒株VR-2332相比，GDgz在非中和表位中有1個A29V的變異，在185位有V185A突變，在189位有I189L突變；中和表位中，在39位有1個L39I突變。這些突變很可能影響蛋白的免疫原性，進而影響疫苗的免疫效果，導致豬場免疫失敗，PRRS暴發和流行。

最近國內新出現了NADC30-like毒株，有研究證實該類毒株具有高致病性，是由中國高致病性毒株和美國NADC30毒株發生重組所致，商品化疫苗對該類型毒株缺乏免疫保護作用。我們分離的毒株與商品化疫苗JXA1-P80高度同源，這些結果表明弱毒苗的廣泛使用導致中國 PRRSV 毒株變異情況越來越嚴重，有必要開發新型疫苗預防疾病發生，早日解決PRRS感染與流行的問題。

［1］ Dekkers J，Rowland R R，Lunney J K，et al. Host genetics of response to porcine reproductive and respiratory syndrome in nursery pigs［J/OL］. Vet Microbiol，2017，http：//www.sciencedirect. com/science/article/pii/S0378113516306782，2017-03-22.

［2］ Ke H，Yoo D. The viral innate immune antagonism and an alternative vaccine design for PRRS virus［J］. Vet Res，1997，28：305-352.

［3］ Snijder E J，Meulenberg J M. The molecular biology of arteriviruses［J］. Vet Microbiol，2017，81（Pt 4）：853–879.

［4］ Han J，Zhou L，Ge X，et al. Pathogenesis and control of the Chinese highly pathogenic porcine reproductive and respiratory syndrome virus［J/ OL］. Vet Microbiol，http：//www.sciencedirect. com/science/article/pii/S0378113516305120，2017-03-02.

［5］ Aeksiri N，Jantafong T. Structural insights into typeⅠ and typeⅡ of nsp4 porcine reproductive and respiratory syndrome virus （nsp4 PRRSV）by molecular dynamics simulations［J］. J Mol Graph Model，2017，74：125-134.

［6］ 劉準，王鳳雪，武華，等. 豬繁殖與呼吸綜合征病毒非結構蛋白研究進展［J］. 動物醫學進展，2010，31（12）：120-124.

［7］ Fang Y，Snijder E L. The PRRSV replicase：Exploring the multifunctionality of an intriguing set of nonstructural proteins［J］. Virus Res，2010，154（1-2）：61-76.

［8］ Wang Q，Li Y，Dong H，et al. Identification of host cellular proteins that interact with the M protein of a highly pathogenic porcine reproductive and respiratory syndrome virus vaccine strain ［J］. Virol J，2017，14（1）：39.

［9］ 李國娟，田永強，李建強，等. 豬繁殖與呼吸綜合征病毒分子生物學研究進展［J］. 生物技術通報，2009（12）：37-41.

［10］ 劉光清，薛 強，仇華吉，等. 豬繁殖與呼吸綜合征病毒CH-1a株非結構基因的分子克隆及其基因特征的研究［J］. 中國預防獸醫學報，2002，24（2）：81-87.

［11］ Leng C L，Tian Z J，Zhang W C，et al. Characterization of two newly emerged isolates of porcine reproductive and respiratory syndrome virus frome Northeast China in 2013［J］. Vet Microbiol，2014，171（1-2）：41-52.

［12］ Won I K，Jac J K，Sang H C，et al. Different biological characteristics of wild-type porcine respiratory syndrome viruses and vaccine viruses and identification of the corresponding gentic determinants［J］. J Clin Microbiol，2008，46（5）：1758-1768.

［13］ Bi C，Shao Z，Zhang Y，et al. Identification of a linear B-cell epitope on non-structural protein 12 of porcine reproductive and respiratory syndrome virus，using a monoclonal antibody［J/OL］. Arch Virol，http：//www.doc88.com/p-1611301780273. html，2017-04-01.

［14］ Li X，Bao H，Wang Y，Tian K Widespread of NADC30-like PRRSV in China：Another Pandora's box for Chinese pig industry as the outbreak of highly pathogenic PRRSV in 2006［J］? Infect Genet Evol. 2017，49：12-13.

［15］ 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.

［16］ Li X，Wu J，Tan F，et al. Genome characterization of two NADC30-like porcine reproductive and respiratory syndrome viruses in China［J］.Springerplus，2016，5（1）：1677.

［17］ Sun Z，Wang J，Bai X，et al. Pathogenicity comparison between highly pathogenic and NADC30-like porcine reproductive and respiratory syndrome virus［J］. Arch Virol，2016，161（8）：2257-2261.

［18］ Zhou Y J，Yu H，Tian Z J，et al. Genetic diversity of the ORF5 gene of porcine reproductive and respiratory syndrome virus isolates in China from 2006 to 2008［J］. Virus Res，2009，144（1-2）：136-144.

［19］ Zhou L，Yang H. Porcine reproductive and respiratory syndrome in China［J］. Virus Res，2010，154（1-2）：31-37.

［20］ Yuan S，Nelsen C J，Murtaugh M P，et al. Recombination between North American strains of porcine reproductive and respiratory syndrome virus［J］. Virus Res，1999，61：87-98.

［21］ Leng X L，Li Z，Xia M，et al. Mutations in the genome of the highly pathogenic porcine reproductive and respiratory syndrome virus potentially related to attenuation［J］. Vet Microbiol，2012，157（1-2）：50-60.

［22］ Kijona F K，Gholamreza H，Denis K G，et al. Genetic variationand phylogentic analyses of the ORF5 gene of acute porcine respiratory and reproductive syndrome viruse isolate［J］. Vet Microbiol，2001，83（3）：249-263.

［23］ Zhang Q Y，Xu X J，You S M，et al. Emerging of two new subgenotypes of porcine reproductive and respiratory syndrome viruses in Southeast China［J］. Microb Pathog，2016，97：27-33.

［24］ 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.

［25］ Zhao K，Ye C，Chang XB，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.

（責任編輯 崔建勛）

Isolation，identification and genetic variation analysis of South China strain GDgz of porcine reproductive and respiratory syndrome virus

YU Lin-yang1，DONG Jian-guo1，2，ZHANG Le-yi1，LIU Yan-ling1，LIANG Peng-shuai1，WANG Lei1，SONG Chang-xu1
（1. College of Animal Science，South China Agricultural University，Guangzhou 510642，China；2. College of Animal Husbandry and Veterinary，Xinyang College of Agriculture and Forestry，Xinyang 464000，China）

To understand the evolution of porcine reproductive and respiratory syndrome virus（PRRSV）in South China，we isolated a PRRSV strain in Marc-145 cell from lung of suspected PRRSV-positive pig in South China and named GDgz. GDgz strain could cause typical cenotaphic effect. The results of phylogenetic tree and homology alignment analysis showed that GDgz belonged to the same cluster with high pathogenetic PRRSV（HPPRRSV） strains in China；GDgz had a far distance with European strain Lelystad virus，the nucleotide homology was only 60.3%，with American classic strain VR-2332 was 88.8%，with Chinese HP-PRRSV strains JXA1 and HuN4 were 98.8% and 98.9%，with the vaccine strain JXA1-P80 of JXA1 was 99.1%，with NADC30 and NADC30-like strains CHsx1401 and JL580 were 84.5%，87.2 and 83.6%，respectively. The result of sequence alignment showedthat there were 30 discontinuous amino acids deletion in the high variable region of NSP2. These results showed that GDgz belonged to vaccine strain of high pathogenetic North American type PRRSV and had some mutations in the epitope region.

porcine；reproductive and respiratory syndrome virus；isolation and identification；genetic variation

S855.3

A

1004-874X（2017）04-0138-08

于林洋，董建國，張樂宜，等. 豬繁殖與呼吸系統綜合征病毒華南株GDgz的分離鑒定及遺傳進化分析［J］.廣東農業科學，2017，44（4）：138-145.

2017-02-19

國家科技支撐計劃項目（2015BAD12B02-5）

于林洋（1992-），男，在讀碩士生，E-mail: 396329755@qq.com

宋長緒（1965-），男，博士，研究員，E-mail: cxsong2004@163.com