萬(wàn)春和 黃 瑜

（1.福建省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所/福建省禽病防治重點(diǎn)實(shí)驗(yàn)室/福建省畜禽疫病防治工程技術(shù)研究中心 福州 350013；2.龍巖學(xué)院/預(yù)防獸醫(yī)學(xué)與生物技術(shù)福建省高校重點(diǎn)實(shí)驗(yàn)室 福建龍巖 364012）

根據(jù)國(guó)際病毒分類委員會(huì) （International Committee on Taxonomy of Viruses，ICTV） 最 新 分 類（https://talk.ictvonline.org/taxonomy/）， 鵝 圓 環(huán) 病 毒（Goose circovirus，GoCV）屬于圓環(huán)病毒科（Circoviridae）圓環(huán)病毒屬（Circovirus）成員。

GoCV 病毒顆粒為無(wú)囊膜、二十面體對(duì)稱、單股負(fù)鏈環(huán)狀閉合DNA 病毒， 直徑約15 nm。 GoCV 于1999 年由Soike 等從生長(zhǎng)發(fā)育不良、 體重下降的患鵝中發(fā)現(xiàn)[1]。 隨后 ，GoCV 在我國(guó)臺(tái)灣[2]、浙江[3]、江西[4]、廣東[5]、山東[6]、福建[7]、重慶[8]、廣西[9]等多地鵝群中發(fā)現(xiàn)。 GoCV 基因組全長(zhǎng)為1 821 nt 左右，其5'-端具有典型的頸環(huán)結(jié)構(gòu)， 其頂部有1 個(gè)保守的9 堿基序列“TATTATT①AC”（①position 1）；其 3'-端基因間區(qū)存在正向或反向的重復(fù)序列。 與其他圓環(huán)病毒屬成員基因組結(jié)構(gòu)類似，GoCV 還編碼有兩個(gè)大的開(kāi)放閱讀框（ORF），位于病毒左側(cè)與病毒復(fù)制相關(guān)的非結(jié)構(gòu)蛋白（Rep）和位于右側(cè)抗原相關(guān)的結(jié)構(gòu)蛋白（Cap）[1-4]。

近年來(lái)，在鴨圓環(huán)病毒（Duck circovirus，DuCV）基因組中除ORF-V1 和ORF-C1 兩個(gè)主要ORF 之外，新鑒定的位于ORF-V1 互補(bǔ)鏈上的第三個(gè)功能性O(shè)RF，命名為 ORF3[10]。 對(duì) DuCV 的 ORF3 基因分析發(fā)現(xiàn)，ORF3 基因在鴨圓環(huán)病毒兩大基因型（DuCV-G1 和 DuCV-G2） 存 在 特 征 性 差 異 ，即DuCV-G1 的 ORF3 基因大小均為 237 bp， 編碼78 個(gè)氨基酸； 而DuCV-G2 的 ORF3 基因大小均為297 bp，編碼98 個(gè)氨基酸；同型之間氨基酸的同源性為 89.9%～100%， 而兩型之間氨基酸的同源性只有72.2%～81.0%。進(jìn)一步對(duì)DuCV-G1 和DuCV-G2 ORF3 基因分析比對(duì)發(fā)現(xiàn)， 由于DuCV-G1 的ORF3基因發(fā)生T236A 突變（由TTG 突變?yōu)門AG，產(chǎn)生新的終止子（Premature stop condon）導(dǎo)致翻譯提前終止， 使 DuCV-G1 比 DuCV-G2 的 ORF3 基因 3'-末端縮短 60 個(gè)核苷酸（即 20 個(gè)氨基酸）[11]。

目前，尚未見(jiàn)GoCV ORF3 基因的研究報(bào)告，本研究通過(guò)分析本團(tuán)隊(duì)前期鑒定的GoCV JX1 株（GenBank 登錄號(hào)：GU320569）[4]基因組 ORF3 基因特征， 為后續(xù)開(kāi)展GoCV ORF3 基因功能研究提供參考和借鑒。

1 材料與方法

1.1 毒株 朗德鵝源 GoCV（JX1 株，GenBank 登錄號(hào)：GU320569）[4]由福建省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所鑒定并保存。

1.2 JX1 株ORF3 基因基本理化特性分析 經(jīng)序列分析，JX1 株ORF3 基因位于ORF-V1 互補(bǔ)鏈上127 至 426 位，其基因全長(zhǎng)為 300 bp，編碼 99 個(gè)氨基酸，分子量為11 583.22 Da，理論等電點(diǎn)（Theoretical pI） 為 8.84； 總原子數(shù)為 1 619， 其組成為C515H802N154O144S4，帶負(fù)電荷氨基酸（Asp+Glu）為 8 個(gè)，帶正電荷氨基酸（Arg+Lys）為 10 個(gè)；不穩(wěn)定指數(shù)為80.95，是一個(gè)不穩(wěn)定蛋白；脂肪系數(shù)89.70，總平均疏水性指數(shù)為-0.472。

1.3 GoCV 序列分析

1.3.1 全基因組遺傳進(jìn)化分析 由于DuCV 不同基因型ORF3 基因存在特征性差異， 為探討 GoCV ORF3 基因是否存在類型現(xiàn)象， 本研究先對(duì)Gen-Bank 中登錄的 GoCV 全基因組序列 （共 67 株GoCV）進(jìn)行遺傳分析。 利用 MEGA6.0 軟件[12]繪制GoCV 全基因組遺傳進(jìn)化樹(shù)， 繪制遺傳進(jìn)化樹(shù)采用NJ 法 （Neighbor-Joining Methods）， 重 復(fù) 計(jì) 算1 000 次（Bootstrap=1 000）。

1.3.2 ORF3 基因遺傳進(jìn)化分析 將上述 67 株GoCV 的 ORF3 基因利用 MEGA6.0 軟件繪制 ORF3基因遺傳進(jìn)化樹(shù)， 繪制遺傳進(jìn)化樹(shù)采用 NJ 法（Neighbor-Joining Methods）， 重 復(fù) 計(jì) 算 1 000 次（Bootstrap=1 000）。

1.3.3 ORF3 基因核苷酸同源性比較 利用Lasergene MegAlign 軟件（By Clustal W Method），將 JX1株與其他66 株GoCV 的ORF3 基因進(jìn)行核苷酸同源性分析， 比較JX1 株與其他ORF3 基因不同遺傳進(jìn)化分支毒株的核苷酸同源性。

2 結(jié) 果

2.1 全基因組遺傳進(jìn)化分析 從GoCV 全基因組遺傳進(jìn)化分析可見(jiàn) （見(jiàn)圖1）， 可分為5 個(gè)基因型（Genotype 1 至 Genotype 5）。 其 中 Genotype 3 型（7 株）、Genotype 4 型（1 株）和 Genotype 5 型（2 株）均為歐洲波蘭分離株， 其余57 株均處于經(jīng)典的Genotype 1 型和Genotype 2 型。需要指出的是，波蘭地區(qū)的GoCV 全基因組遺傳進(jìn)化均為復(fù)雜[13]，Genotype 2a1 亞型和Genotype 2a2 亞型也有分布， 且歐洲其他地區(qū)流行的GoCV 也均處于Genotype 2 型內(nèi)的Genotype 2a1 亞型和Genotype 2a2 亞型。

2.2 ORF3 基因組遺傳進(jìn)化分析 從ORF3 基因組遺傳進(jìn)化分析可見(jiàn)（見(jiàn)圖2），同樣可分為5 個(gè)遺傳進(jìn)化分支（Clade 1 至Clade 5），但其與全基因組基因分型存在很大不同，Clade 1 分支GoCV 均為中國(guó)株，其中中國(guó)臺(tái)灣省GoCV 均處于Clade 1a 亞分支，中國(guó)其他地區(qū)GoCV 則分別處于Clade 1b 亞分支、Clade 1c 亞分支和 Clade 1d 亞分支； 歐洲地區(qū)GoCV 則分處于 Clade 2 分支、Clade 3 分支、Clade 4分支和Clade 5 分支。

2.3 ORF3 基因分析 對(duì)所有 67 株 GoCV 的ORF3 基因進(jìn)行分析發(fā)現(xiàn)， 其中 60 株 GoCV 的ORF3 基因長(zhǎng)度與 JX1 株一致， 均為 300 bp。 但是6 株 GoCV （G1 株、G4 株、G9 株、G15 株、G16 株、GoCV-Hun2 株）的 ORF3 基因長(zhǎng)度為 309 bp；需要指出的是， 該6 株GoCV 在ORF3 基因進(jìn)化上均處于Clade 2 分支， 但處于Clade 2 分支的其他GoCV的ORF3 基因長(zhǎng)度均為300 bp， 關(guān)于該差異的生物學(xué)意義還需要作更進(jìn)一步深入研究。

2.4 ORF3 基因核苷酸同源性比較 從核苷酸同源性分析比較結(jié)果可見(jiàn)，JX1 株與 yk1 株（AY633653）、yk3 株（DQ192280）的核苷酸同源性最高，均為100%；與DG1 株的核苷酸同源性最低，僅為89.0%。

3 討 論

圖1 GoCV 全基因組遺傳進(jìn)化分析

圖2 ORF3 基因組遺傳進(jìn)化分析

目前， 對(duì)于圓環(huán)病毒科圓環(huán)病毒屬豬圓環(huán)病毒1 型和 2 型（PCV1 和 PCV2）的 ORF3 基因功能研究較多，且研究發(fā)現(xiàn)，PCV1 和 PCV2 與 DuCV 類似，其中PCV1 （對(duì)豬群不致?。?的ORF3 基因大小均為237 bp， 編碼 78 個(gè)氨基酸；PCV2 （對(duì)豬群致病）的ORF3 基因大小均為 315 bp， 編碼 104 個(gè)氨基酸[14]。Liu 等通過(guò)構(gòu)建ORF3 缺失突變體病毒與野生型病毒比較發(fā)現(xiàn)，ORF3 缺失突變株不影響PCV2 復(fù)制，但參與了外源性細(xì)胞凋亡。 通過(guò)測(cè)定caspase 3、caspase 8 和 caspase 9 基因表達(dá)發(fā)現(xiàn)，PCV2 誘導(dǎo)的細(xì)胞凋亡需要激活caspase 8 （而不是caspase 9），caspase 8 的激活導(dǎo)致caspase 3 的激活，使感染細(xì)胞中caspase 底物的裂解。 研究還發(fā)現(xiàn)，PCV2 的ORF3與 P53 競(jìng)爭(zhēng)結(jié)合 pPirh2，通過(guò)干擾 pPirh2 對(duì) P53 泛素化， 從而誘發(fā)了P53 介導(dǎo)的凋亡通路。 PCV2 的ORF3 蛋白通過(guò)與胞漿中G 蛋白信號(hào)調(diào)節(jié)因子（RGS16）相互作用，導(dǎo)致RGS16 發(fā)生泛素介導(dǎo)的蛋白酶降解，從而促進(jìn)NF-kB 通過(guò)ERK1/2 信號(hào)通路移位胞核內(nèi)，增強(qiáng)IL-6 和IL-8 的轉(zhuǎn)錄，因而在細(xì)胞周邊誘發(fā)嚴(yán)重的炎性反應(yīng)和白細(xì)胞浸潤(rùn)， 這可能也與PCV2 致病性相關(guān)[15-16]。研究還發(fā)現(xiàn)，過(guò)表達(dá)PCV2的ORF3 蛋白可誘導(dǎo)豬外周血淋巴細(xì)胞（PBMCs）凋亡[17]。 但有意思的是，對(duì) PCV1 的 ORF3 蛋白誘導(dǎo)的細(xì)胞凋亡研究結(jié)果表明，PCV1 的 ORF3 蛋白比PCV2 的ORF3 蛋白可誘導(dǎo)更多細(xì)胞發(fā)生凋亡[18]。

對(duì)于 DuCV 的 ORF3 基因功能研究發(fā)現(xiàn)，DuCV-G1 的 ORF3 比 DuCV-G2 的 ORF3 對(duì) 禽 源DF-1 細(xì)胞表現(xiàn)出更高的凋亡能力； 且缺失DuCVG2 的 ORF3 末端20 個(gè)氨基酸后的 ORF3 （長(zhǎng)度與DuCV-G1 的ORF3 一致）表現(xiàn)出更高的促宿主細(xì)胞凋亡能力[19-20]。 但是， 通過(guò)對(duì) 67 株 DuCV-G1 的ORF3 基因分析發(fā)現(xiàn)，不同遺傳進(jìn)化分支的ORF3 基因長(zhǎng)度卻與 DuCV 兩個(gè)基因型 （DuCV-G1 和DuCV-G2）、PCV 兩個(gè)基因型（PCV1 和 PCV2）存在特征性差異及核定位信號(hào)差異明顯不同， 雖然有6株GoCV 的 ORF3 基因?yàn)?309 bp，但是與其他 61 株GoCV 的ORF3 基因編碼區(qū)氨基酸序列一樣，均不存在典型的核定位信號(hào)[21]。

研究還發(fā)現(xiàn)，GoCV 全基因組與ORF3 基因的遺傳進(jìn)化分析結(jié)果存在明顯不同， 雖然絕大部分GoCV 從全基因組遺傳進(jìn)化分析也可以同樣分為兩個(gè)大的基因型， 但是卻與ORF3 基因的遺傳進(jìn)化分析結(jié)果存在明顯不同， 尤其是ORF3 基因的遺傳進(jìn)化表現(xiàn)為明顯的地域性，但并不表現(xiàn)出時(shí)間性；上述結(jié)果表明，ORF3 基因的進(jìn)化頻率并不完全與GoCV全基因組一致， 具體機(jī)制還需要作更進(jìn)一步的深入研究。