畜牧兽医学报  2020, Vol. 51 Issue (6): 1399-1407. DOI: 10.11843/j.issn.0366-6964.2020.023    PDF    
貉源阿留申病毒VP2和NS1蛋白的抗原性预测
吴艳虹, 赵永强, 韦韬, 章秀婷, 丛丽, 岳志刚, 肖家美, 邵西群     
中国农业科学院特产研究所, 长春 130112
摘要:旨在预测和分析貉源阿留申病毒(RFAV)VP2和NS1蛋白的抗原表位特征,筛选阿留申病毒属较保守的B、T细胞抗原表位。本研究对RFAV的近全长基因组进行克隆及测序,对其VP2和NS1基因编码蛋白的理化性质、二级结构、翻译后修饰位点和抗原表位进行预测,并将预测的修饰位点、抗原表位与其他阿留申病毒种的序列进行比较分析,筛选相对保守的修饰位点和抗原表位。结果显示,获得的RFAV基因组长4 327 bp,编码VP2蛋白的636个氨基酸,编码NS1蛋白的641个氨基酸。VP2和NS1蛋白均为亲水性蛋白,二级结构以无规则卷曲为主。在阿留申病毒属内,VP2蛋白有3个保守的B细胞抗原表位,3个保守的T细胞表位,8个保守的翻译后修饰位点;NS1蛋白有1个保守的B细胞抗原表位,2个保守的T细胞抗原表位和7个保守的翻译后修饰位点。本研究成功克隆了RFAV近全长基因组序列,全面对RFAV的VP2、NS1蛋白抗原表位、翻译后修饰位点进行预测,并分析其在阿留申病毒属内的保守和变异特征,为阿留申病毒免疫研究提供参考。
关键词阿留申病毒    VP2蛋白    NS1蛋白    抗原表位    翻译后修饰位点    
Cloning and Epitopes Analysis of VP2 and NS1 Proteins of Raccoon Dog and Arctic Fox Amdoparvovirus
WU Yanhong, ZHAO Yongqiang, WEI Tao, ZHANG Xiuting, CONG Li, YUE Zhigang, XIAO Jiamei, SHAO Xiqun     
Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun 130112, China
Abstract: This study aimed to analyze the epitope characteristics of VP2 and NS1 proteins of raccoon dog and arctic fox amdoparvovirus (RFAV), and further select the conserved B- and T-cell epitopes on the VP2 and NS1 proteins of RFAV in Amdoparvovirus genus. In the present study, the RFAV genome except the 5' and 3' end of the genome was cloned and sequenced, then the bioinformatics programs were used to predict the physical and chemical properties, Secondary structure, post-translational modifications and potential B- and T-cell epitopes of VP2 and NS1 proteins, and further compared the potential epitopes and post-translational modifications with the typical members within Amdoparvovirus genus. The 4 327 bp RFAV genomic sequence in length, encoding 636 amino acids of VP2 protein and 641 amino acids of NS1 protein, was obtained. Both VP2 and NS1 proteins are hydrophilic proteins and the secondary structure were mainly composed of random coils. We predicted there were 3 conserved B-cell epitopes, 3 conserved T-cell epitopes and 8 conserved post-translational modification sites on VP2 protein within the Amdoparvovirus genus; while NS1 protein had 1, 2 and 7 in the respective conserved items. In this study, we cloned the main genome sequence of RFAV, and predicted the potential epitopes and post-translational modifications on VP2 and NS1 proteins of RFAV. This study provides fundamental basis for the development of a candidate vaccine and disease diagnosis of amdoparvoviruses.
Key words: amdoparvovirus    VP2 protein    NS1 protein    epitope    post-translational modification sites    

阿留申病毒(amdoparvovirus)归属于细小病毒科(Parvoviridae)阿留申病毒属(Amdoparvovirus genus)[1]。目前,该属成员包括水貂阿留申病毒(Aleutian mink disease virus, AMDV)、小熊猫阿留申病毒(Red panda amdoparvovirus, RpAPV)、臭鼬阿留申病毒(Skunk amdoparvovirus, SKAV)、貉狐阿留申病毒(Raccoon dog and arctic fox amdoparvovirus, RFAV)、赤狐阿留申病毒(Red fox fecal amdoparvovirus, RFFAV)和灰狐阿留申病毒(Gray fox amdoparvovirus, GFAV)[2-8]。阿留申病毒为单链DNA病毒,基因组长约5 kb,其基因组编码结构蛋白VP1、VP2和非结构蛋白NS1、NS2、NS3[9]。VP2是主要的衣壳蛋白,与病毒感染宿主范围、组织趋向性和毒力相关,是开发疫苗的候选蛋白[10]。NS1蛋白是最大的非结构蛋白,与病毒DNA复制、转录调控和衣壳组装密切相关,在病毒的复制过程中起重要作用[11]。AMDV可引起水貂的慢性疾病和持续性感染,以高蛋白血症、浆细胞增多和免疫复合物沉积为主要特征[12-14]。目前,阿留申病毒尚无有效疫苗和特异性药物[15],预测潜在的表位、选择保护性表位以及构建表位疫苗是提高疫苗保护性免疫的关键。

本试验选取3个不同地区的6株RFAV进行克隆测序,使用不同软件对其VP2和NS1蛋白的理化性质、二级结构、翻译后修饰位点和抗原表位进行预测,并将预测的修饰位点、抗原表位与本实验室之前获得的RFAV及NCBI上公布的其他阿留申病毒种的序列进行比较分析,以期筛选出较保守的T细胞和B细胞抗原表位。此外,将已经验证的AMDV优势抗原表位与RFAV进行比较,总结其在RFAV中的特征。为进一步研究阿留申病毒VP2、NS1蛋白的抗原特征与免疫表位的发掘奠定基础。本研究使用不同的软件进行预测和综合分析,克服了单一参数模型的局限性,通过与阿留申病毒属的病毒进行比较分析,显著提高了抗原表位预测的准确性,为阿留申病毒疫苗的研制提供理论依据。

1 材料与方法 1.1 样品的采集和预处理

乌苏里貉血清样品采自3个不同地区的乌苏里貉养殖场。将新鲜的血液样本500ⅹg (Eppendorf 5804/5804R, 德国)离心10 min或放于4 ℃冰箱中过夜,收集血清,-20 ℃冰箱中保存。用于RFAV检测和克隆。

1.2 主要试剂与仪器

血液DNA直接PCR试剂盒购自辽宁生物医药有限公司;Trans 2K DNA Marker、胶回收试剂盒、质粒小提试剂盒、大肠杆菌Trans感受态细胞、pEASY-T1载体购自北京全式金生物技术有限公司;生物安全柜购自上海博讯实业有限公司;NanoDrop核酸浓度测定仪ND2000购自Thermo公司;离心机购自Sigma公司。AMDV阳性、阴性血清由本实验室保存。

1.3 RFAV的克隆与测序

1.3.1 引物设计   根据NCBI上公布的阿留申病毒的序列信息,在保守区域选择引物,并通过Oligo 7.0 (Molecular Biology Insights, Inc.,Cascade, CO, USA, http://www.oligo.net)进行验证。引物由大连TaKaRa生物技术公司合成。本试验引物序列见表 1

表 1 引物信息 Table 1 Primer information

1.3.2 克隆及测序   用本实验室设计的阿留申病毒保守引物和血液DNA直接PCR试剂盒对乌苏里貉血清样品进行检测,AMDV阳性、阴性血清样本分别作为阳性和阴性对照。3个地区各筛选2份阳性样品进行克隆,分别命名为DA-RD17、DA-RD15、ZJ-L3、ZJ-L5、HB-Y4和HB-Y5。PCR总反应体系为40 μL,3.6 μL核酸释放剂加入到2.0 μL血清样本中裂解10 min。PCR反应液:32 μL PCR reaction buffer,上、下游引物(表 1)各1.5 μL,1 μL Taq酶,加RNase-free H2O至总反应体系为40 μL。PCR反应条件:95 ℃预变性5 min;94 ℃变性40 s,52~53 ℃ (表 1)退火35 s,72 ℃延伸1~2 min (表 1),35个循环;72 ℃ 10 min。PCR产物用1.5%的琼脂糖电泳验证。阳性产物用胶回收试剂盒纯化,连接pEASY-T1载体,转化Trans-T1感受态细胞,在LB固体培养基中过夜培养,挑取单菌落于LB液体培养基中扩大培养,菌液PCR鉴定正确后送公司测序(库美生物,长春)。用Lasergene软件(DNAStar)对序列进行组装分析,用MEGA 7.0将VP2和NS1基因序列翻译成蛋白序列[16]

1.4 RFAV VP2和NS1蛋白序列预测分析

利用ExPASy’s预测服务器ProtParam (https://web.expasy.org/protparam/)预测VP2、NS1蛋白基本理化性质[17];GOR4(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_gor4.html)预测蛋白二级结构[18];Motif Scan(https://myhits.isb-sib.ch/cgi-bin/motif_scan)预测翻译后修饰位点[19];Protean(DNAStar)分析其亲水性指数、柔韧性指数、可及性参数以及抗原指数,推测B细胞抗原表位的可能位置,并结合ABCpred(https://webs.iiitd.edu.in/raghava/abcpred/ABC_submission.html)和IEDB(http://www.iedb.org/)方案预测结果筛选一致性B细胞线性表位[20-22];联合Protean和IEDB在线软件预测T细胞抗原表位,筛选一致性T细胞表位。

以AMDV-G株(JN040434.1)的VP2、NS1蛋白序列作为残基定位参考。对预测的RFAV的B、T细胞抗原表位与NCBI上公布的RFAV、RpAPV、SKAV、RFFAV、AMDV的VP2、NS1蛋白序列进行比对,分析各表位氨基酸差异,排除突变集中区域,综合评价VP2和NS1蛋白的抗原表位区域,最终筛选相对保守的抗原表位。并对翻译后修饰位点的保守性进行分析,筛选完全保守的翻译后修饰位点,可进一步研究翻译后修饰位点对蛋白抗原性的影响。

2 结果 2.1 RFAV近全长序列的扩增

RFAV血清样品PCR分段扩增后,产物用1.5%琼脂糖凝胶电泳检测,结果如图 1,片段大小与目的条带基本符合,切胶回收并加载到pEASY-T1载体测序。

A、B.NS-1、NS-2片段;M. Trans 2K DNA分子质量标准;1~6.RFAV样品DA-RD17、DA-RD15、ZJ-L3、ZJ-L5、HB-Y4和HB-Y5 A, B. PCR product of NS-1, NS-2; M. Trans 2K DNA marker; 1-6. RFAV sample of DA-RD17, DA-RD15, ZJ-L3, ZJ-L5, HB-Y4 and HB-Y5 图 1 RFAV PCR产物 Fig. 1 PCR amplification of RFAV
2.2 RFAV蛋白基本理化性质分析

本研究获得的RFAV基因组长4 327 bp,编码VP2蛋白的636个氨基酸,编码NS1蛋白的641个氨基酸。RFAV样品DA-RD17、DA-RD15、ZJ-L3、ZJ-L5、HB-Y4和HB-Y5的VP2蛋白相似性为96.5%~99.4%,NS1蛋白相似性为92.5%~97.7%。以DA-RD7为代表进行RFAV的抗原性分析。PortParam分析得到RFAV的VP2蛋白理论等电点为5.49,带正、负电荷数氨基酸分别为58、77个,总平均亲水性为-0.730,为亲水性蛋白质;NS1蛋白理论等电点为8.64,带正、负电荷数氨基酸分别为84、75个,总平均亲水性为-0.635,为亲水性蛋白质。B细胞抗原表位应具有亲水性高的特征,有利于与抗体的结合,因此,亲水性与蛋白的B细胞抗原表位密切相关。

2.3 RFAV蛋白二级结构分析

GOR4预测结果显示,VP2蛋白二级结构中,α螺旋占17.77%,β片层占25.47%,无规则卷曲占56.76%。NS1蛋白α螺旋占28.97%,β片层占18.54%,无规则卷曲形成占52.49%。说明VP2和NS1蛋白二级结构以无规则卷曲为主。α螺旋、β片层结构规则不易形变,与抗体嵌合能力较差,且经常处于蛋白质内部,因而很少成为抗原表位。无规则卷曲结构易于扭曲、盘旋并位于蛋白质分子表面,有利于抗体嵌合,为B细胞优势抗原表位的可能性较大[23]

2.4 RFAV翻译后修饰位点分析

利用在线软件MotifScan(https://myhits.isb-sib.ch/cgi-bin/motif_scan)预测RFAV蛋白翻译后修饰位点,结果如表 2

表 2 RFAV蛋白翻译后修饰位点预测 Table 2 Post-translational modification sites of deduced RFAV proteins
2.5 RFAV蛋白B、T细胞抗原表位分析

B细胞抗原表位根据Protean、IEDB、ABCpred软件综合分析,VP2蛋白筛选出6个B细胞抗原表位,分别位于296—303、330—338、415—424、519—530、560—568和610—619位氨基酸残基附近。利用同样的方法对NS1蛋白进行综合分析,筛选出4个B细胞抗原表位,分别位186—211、333—338、528—541和581—586位氨基酸残基附近(表 3)。以上B细胞抗原表位均位于蛋白的亲水区和无规则卷曲区。

表 3 RFAV主要蛋白抗原表位 Table 3 Antigenic epitopes of RFAV proteins

T细胞抗原表位分别用AMPHI、Rothbard-Taylor和IEDB软件预测,预测结果有差异,也有部分重合,选择重合的部分作为优势T细胞抗原表位。VP2蛋白优势T细胞抗原表位有4个:53—59、197—211、275—310和396—412,而NS1蛋白有2个:225—239和588—602(表 3)。

VP2蛋白部分翻译后修饰位点在预测的RFAV抗原表位内,但多为磷酸化位点,如酪蛋白激酶Ⅱ磷酸化位点304—307、421—424,蛋白激酶C磷酸化位点283—285、297—302,酪氨酸激酶磷酸化位点85—93;RFAV的NS1蛋白的翻译后修饰位点有部分位点在预测的RFAV的抗原表位内,如N-糖基化位点205—206、591—594、599—602,酪蛋白激酶Ⅱ磷酸化位点200—203、338,N-酰基化位点527—532,蛋白激酶C磷酸化位点596—598(表 2)。翻译后修饰对于蛋白的功能发挥具有重要意义,尤其是糖基化修饰,经糖基化修饰的蛋白可以引起特定的免疫反应、降低蛋白抗原的免疫原性,也有研究表明蛋白的糖基化修饰可能在病毒的免疫逃逸过程起着重要的作用[24]

2.6 阿留申病毒蛋白保守抗原表位、翻译后修饰位点分析

以AMDV-G(JN040434.1)株的VP2、NS1蛋白序列为残基位点定位序列,确定RFAV的VP2、NS1蛋白抗原表位、翻译后修饰位点所在氨基酸残基位置。将筛选的VP2、NS1蛋白B/T细胞抗原表位和翻译后修饰位点与本研究克隆的其他5份RFAV样品及阿留申病毒属内病毒相应蛋白进行比较分析。结果显示VP2蛋白B细胞抗原表位296—303、330—338、415—424在阿留申病毒属内变异位点较多,519—530、560—568、610—619相对保守(图 2A);T细胞抗原表位275—310、396—412在阿留申病毒属内变异位点较多,其中275—298保守;而53—59、197—211相对保守(图 2B)。NS1蛋白B细胞抗原表位186—211、333—338、581—586在阿留申病毒属内变异位点较多,528—541相对保守;T细胞抗原表位225—239、588—602在阿留申病毒属内变异位点较多,但区段231—239、590—600相对保守(图 2C)。

A. VP2蛋白B细胞抗原表位;B.VP2蛋白T细胞抗原表位;C. NS1蛋白B、T细胞抗原表位 A. Potential B cell epitopes of VP2 protein; B. Potential T cell epitopes of VP2 protein; C. Potential B and T cell epitopes of NS1 protein 图 2 VP2、NS1蛋白的抗原表位分析 Fig. 2 The comparison of the antigenic epitopes on VP2 and NS1 proteins of the typical members within Amdoparvovirus

VP2蛋白翻译后修饰位点117SPK119147TVTE150283STR285358GLQGSY363390STR392447SWEE450459THFE462553NLTD556在阿留申病毒属中完全保守,NS1蛋白翻译后修饰位点38YT39318SKK320340ED341399TL400423GK424435GTGKTL440456TTS458在阿留申病毒属内完全保守(图略)。

3 讨论

本研究克隆并测序了来自不同地区的6株RFAV,为进一步开展RFAV分子研究提供材料。本实验室已收集RFAV丰富序列信息,为RFAV的主要蛋白VP2和NS1的免疫相关抗原性系统比较分析奠定基础。目前阿留申病毒疫苗的研制主要集中在VP2和NS1蛋白,但均为B细胞抗原表位的预测和验证,T细胞抗原表位尚未解析。已有研究表明AMDV-G的VP2蛋白386QQNFSTRYIYDRN398428HTAKQPKLRTPPIHHSKID446459THFEDEV-IYLDYFNF473为B细胞优势抗原表位[25-28],NS1蛋白302TKVTSKKEVANPVQQPSKKL321为B细胞优势抗原表位[29]。以AMDV非结构蛋白NS1接种水貂,只对接种水貂产生部分保护[30]。将上述优势抗原表位与RFAV进行比较分析,其中,459THFEDE-VIYLDYFNF473在RFAV的VP2蛋白中完全保守。其他优势表位在RFAV中有部分位点发生突变,如VP2蛋白Y395→395F、R436→436S的突变,NS1蛋白K308→308R/G、N312→312E、V314→A314、Q315→315K、L321→321M的突变,但通常是性质相同、相近或种类相同的氨基酸之间发生相互替换。这些位置的氨基酸可能在抗原抗体的识别、结合和相互作用中起重要作用,这些表位也可能是RFAV的优势抗原表位。

本研究结合RFAV的VP2、NS1蛋白的二级结构、亲疏水性特征对其B、T细胞抗原表位进行筛选及分析。结果表明,预测的B细胞抗原表位位于蛋白的无规则卷曲区域,具有亲水性高、柔韧性好、表面可及性大、抗原性指数高的特征。预测的T细胞抗原表位是基于MHC-Ⅱ类分子与抗原肽的结合的分值高低进行筛选,主要是线性短肽。部分B、T细胞抗原表位覆盖翻译后修饰位点。研究表明参与免疫反应的大多数分子是经糖基化修饰后的糖蛋白,例如免疫球蛋白、细胞因子、补体系统等[31-32]。人免疫缺陷病毒(HIV)的gp120分子高度糖基化使机体免疫系统误将gp120分子识别为自身分子或通过空间位阻作用使病毒更易于逃避中和抗体的中和作用[31]。丙型肝炎病毒(HCV)的包膜糖蛋白上的N-糖苷可影响某些抗原表位而抑制适应性免疫应答,使病毒逃逸机体的免疫攻击[32]。说明蛋白的糖基化修饰可能在病毒的免疫逃逸过程起着重要的作用。而磷酸化修饰的病毒蛋白参与调控病毒复制、病毒增殖和病毒粒子装配等一系列病毒的代谢活动,同时,也影响宿主细胞内的信号转导,抑制宿主基因组复制和表达[33-36]。位于RFAV的VP2、NS1蛋白B、T细胞抗原表位区内的翻译后修饰位点是否与蛋白抗原性相关,尚需进一步验证。目前,研究抗原糖基化与免疫之间的关系, 通过糖基化的改变使病毒抗原成为有效的抗病毒疫苗是一种新的研究思路。

在阿留申病毒属内,RFAV的VP2蛋白B细胞抗原表位519—530、560—568、610—619相对保守,这与之前试验确定的AMDV的VP2蛋白B细胞抗原表位386—398、428—446和459—473不同[25-28],但通过ABCpred软件预测的RFAV的VP2蛋白B细胞抗原表位422—438、442—458可覆盖AMDV的428—448表位,IEDB软件预测的RFAV的B细胞抗原表位382—388、435—444、462—471与AMDV的386—398、428—446、459—473部分覆盖。NS1蛋白B细胞抗原表位528—541相对保守。这与之前试验确定的AMDV的NS1蛋白B细胞抗原表位302—321不同,但通过ABCpred、IEDB、Protean软件预测的RFAV的NS1蛋白B细胞抗原表位298—304、319—324与AMDV的302—321部分覆盖。说明软件预测有一定的可靠性,推测阿留申病毒尚有未发现或验证的抗原表位。VP2蛋白T细胞抗原表位53—59、197—211、275—298相对保守,NS1蛋白T细胞抗原表位231—239、590—600相对保守。但哪些表位能诱导机体发生特异性细胞免疫和体液免疫还需进一步验证。

4 结论

成功克隆了貉源阿留申病毒近全长基因组序列(4 327 bp),全面对RFAV的VP2、NS1蛋白抗原表位、翻译后修饰位点进行预测,并分析其在阿留申病毒属内的保守和变异特征,为阿留申病毒免疫研究提供参考。

参考文献
[1] KING A M Q, LEFKOWITZ E J, ADAMS M J, et al. Virus taxonomy:ninth report of the international committee on taxonomy of viruses[M]. Amsterdam: Elsevier Academic Press, 2011: 405-425.
[2] PORTER D D, LARSEN A E, PORTER H G. The pathogenesis of Aleutian disease of mink:I. In vivo viral replication and the host antibody response to viral antigen[J]. J Exp Med, 1969, 130(3): 575–593. DOI: 10.1084/jem.130.3.575
[3] KENYON A J, GANDER J E, LOPEZ C, et al. Isolation of Aleutian mink disease virus by affinity chromatography[J]. Science, 1973, 179(4069): 187–189. DOI: 10.1126/science.179.4069.187
[4] ALEX C E, KUBISKI S V, LI L L, et al. Amdoparvovirus infection in red pandas (Ailurus fulgens)[J]. Vet Pathol, 2018, 55(4): 552–561. DOI: 10.1177/0300985818758470
[5] CANUTI M, DOYLE H E, BRITTON A P, et al. Full genetic characterization and epidemiology of a novel amdoparvovirus in striped skunk (Mephitis mephitis)[J]. Emerg Microbes Infect, 2017, 6(5): e30.
[6] SHAO X Q, WEN Y J, BA H X, et al. Novel amdoparvovirus infecting farmed raccoon dogs and arctic foxes[J]. Emerg Infect Dis, 2014, 20(12): 2085–2088. DOI: 10.3201/eid2012.140289
[7] BODEWES R, RUIZ-GONZALEZ A, SCHAPENDONK C M, et al. Viral metagenomic analysis of feces of wild small carnivores[J]. Virol J, 2014, 11: 89. DOI: 10.1186/1743-422X-11-89
[8] LI L L, PESAVENTO P A, WOODS L, et al. Novel amdovirus in gray foxes[J]. Emerg Infect Dis, 2011, 17(10): 1876–1878. DOI: 10.3201/eid1710.110233
[9] HUANG Q F, LUO Y, CHENG F, et al. Molecular characterization of the small nonstructural proteins of parvovirus Aleutian mink disease virus (AMDV) during infection[J]. Virology, 2014, 452-453: 23–31. DOI: 10.1016/j.virol.2014.01.005
[10] COSTELLO, STEENFOS, JENSEN, et al. Epitope mapping of Aleutian mink disease parvovirus virion protein VP1 and 2[J]. Scand J Immunol, 1999, 49(4): 347–354. DOI: 10.1046/j.1365-3083.1999.00499.x
[11] CANUTI M, WHITNEY H G, LANG A S. Amdoparvoviruses in small mammals:expanding our understanding of parvovirus diversity, distribution, and pathology[J]. Front Microbiol, 2015, 6: 1119.
[12] MCKENNA R, OLSON N H, CHIPMAN P R, et al. Three-dimensional structure of Aleutian mink disease parvovirus:implications for disease pathogenicity[J]. J Virol, 1999, 73(8): 6882–6891. DOI: 10.1128/JVI.73.8.6882-6891.1999
[13] GONG Q L, LI D, DIAO N C, et al. Mink Aleutian disease seroprevalence in China during 1981-2017:A systematic review and meta-analysis[J]. Microb Pathog, 2020, 139: 103908. DOI: 10.1016/j.micpath.2019.103908
[14] FARID A H, HUSSAIN I. A comparison between intraperitoneal injection and intranasal and oral inoculation of mink with Aleutian mink disease virus[J]. Res Vet Sci, 2019, 124: 85–92.
[15] KOWALCZYK M, JAKUBCZAK A, HORECKA B, et al. A comparative molecular characterization of AMDV strains isolated from cases of clinical and subclinical infection[J]. Virus Genes, 2018, 54(4): 561–569.
[16] KUMAR S, STECHER G, TAMURA K. MEGA7:molecular evolutionary genetics analysis version 7. 0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870–1874. DOI: 10.1093/molbev/msw054
[17] WILKINS M R, GASTEIGER E, BAIROCH A, et al. Protein identification and analysis tools in the ExPASy server[J]. Methods Mol Biol, 1999, 112: 531–552.
[18] XIA F, DOU Y, LEI G Q, et al. FPGA accelerator for protein secondary structure prediction based on the GOR algorithm[J]. BMC Bioinformatics, 2011, 12(1): S5.
[19] NIKITINA T V, TISHCHENKO L I. Computational search for potential posttranslational modification sites in human RNA polymerase Ⅲ subunits[J]. Mol Biol (Mosk), 2005, 39(3): 437–444.
[20] PLASTERER T N. Protean:protein sequence analysis and prediction[J]. Mol Biotechnol, 2000, 16(2): 117–125. DOI: 10.1385/MB:16:2:117
[21] COSTA J G, FACCENDINI P L, SFERCO S J, et al. Evaluation and comparison of the ability of online available prediction programs to predict true linear B-cell epitopes[J]. Protein Pept Let, 2013, 20(6): 724–730. DOI: 10.2174/0929866511320060011
[22] VITA R, MAHAJAN S, OVERTON J A, et al. The Immune Epitope Database (IEDB):2018 update[J]. Nucleic Acids Res, 2019, 47(D1): D339–D343. DOI: 10.1093/nar/gky1006
[23] 黄莉, 谢芝勋, 谢丽基, 等. 禽呼肠孤病毒σB和σC蛋白二级结构及其细胞抗原表位预测[J]. 中国畜牧兽医, 2016, 43(11): 2880–2885.
HUANG L, XIE Z X, XIE L J, et al. Prediction of secondary structure and B-and T-cell epitopes of σB and σC proteins of avian reovirus[J]. China Animal Husbandry & Veterinary Medicine, 2016, 43(11): 2880–2885. (in Chinese)
[24] 王晓龙, 王秀然, 卢天成. 蛋白质糖基化修饰的研究进展[J]. 基因组学与应用生物学, 2017, 36(10): 4380–4384.
WANG X L, WANG X R, LU T C. The research progress of protein glycosylation modification[J]. Genomics and Applied Biology, 2017, 36(10): 4380–4384. (in Chinese)
[25] BLOOM M E, BEST S M, HAYES S F, et al. Identification of Aleutian mink disease parvovirus capsid sequences mediating antibody-dependent enhancement of infection, virus neutralization, and immune complex formation[J]. J Virol, 2001, 75(22): 11116–11127. DOI: 10.1128/JVI.75.22.11116-11127.2001
[26] YI L, CHENG Y N, ZHANG M, et al. Identification of a novel Aleutian mink disease virus B-cell epitope using a monoclonal antibody against VP2 protein[J]. Virus Res, 2016, 223: 39–42. DOI: 10.1016/j.virusres.2016.06.014
[27] LU T F, WANG Y Z, GE J W, et al. Identification and characterization of a novel B-cell epitope on Aleutian Mink disease virus capsid protein VP2 using a monoclonal antibody[J]. Virus Res, 2018, 248: 74–79. DOI: 10.1016/j.virusres.2017.12.008
[28] CHEN X W, SONG C L, LIU Y, et al. Development of an enzyme-linked immunosorbent assay based on fusion VP2332-452 antigen for detecting antibodies against Aleutian mink disease virus[J]. J Clin Microbiol, 2016, 54(2): 439–442. DOI: 10.1128/JCM.02625-15
[29] 马凡舒, 张蕾, 王洋, 等. 水貂阿留申病病毒NS1蛋白B细胞抗原表位的分析与鉴定[J]. 中国兽医科学, 2016, 46(2): 210–216.
MA F S, ZHANG L, WANG Y, et al. Analysis and identification of B-cell epitopes on Aleutian mink disease virus NS1 protein[J]. Chinese Veterinary Science, 2016, 46(2): 210–216. (in Chinese)
[30] AASTED B, ALEXANDERSEN S, CHRISTENSEN J. Vaccination with Aleutian mink disease parvovirus (AMDV) capsid proteins enhances disease, while vaccination with the major non-structural AMDV protein causes partial protection from disease[J]. Vaccine, 1998, 16(11-12): 1158–1165.
[31] 向田, 章晓联. 病毒与宿主细胞的糖基化修饰及相关功能[J]. 生物化学与生物物理进展, 2017, 44(10): 898–907.
XIANG T, ZHANG X L. Glycosylation modification and related functions of virus and host cells[J]. Progress in Biochemistry and Biophysics, 2017, 44(10): 898–907. (in Chinese)
[32] REN Y S, MIN Y Q, LIU M, et al. N-glycosylation-mutated HCV envelope glycoprotein complex enhances antigen-presenting activity and cellular and neutralizing antibody responses[J]. Biochim Biophys Acta Gen Subj, 2016, 1860(8): 1764–1775. DOI: 10.1016/j.bbagen.2015.08.007
[33] 李芒芒, 周倩, 王鹏, 等. 磷酸化病毒蛋白的生物学功能及形成机制[J]. 微生物学通报, 2015, 42(6): 1133–1140.
LI M M, ZHOU Q, WANG P, et al. Biological function and formation mechanism of phosphorylated viral protein[J]. Microbiology China, 2015, 42(6): 1133–1140. (in Chinese)
[34] NUESCH J P F, CORBAU R, TATTERSALL P, et al. Biochemical activities of minute virus of mice nonstructural protein NS1 are modulated in vitro by the phosphorylation state of the polypeptide[J]. J Virol, 1998, 72(10): 8002–8012. DOI: 10.1128/JVI.72.10.8002-8012.1998
[35] GANAIE S S, ZOU W, XU P, et al. Phosphorylated STAT5 directly facilitates parvovirus B19 DNA replication in human erythroid progenitors through interaction with the MCM complex[J]. PLoS Pathog, 2017, 13(5): e1006370. DOI: 10.1371/journal.ppat.1006370
[36] NüESCH J P F, ROMMELAERE J. A viral adaptor protein modulating casein kinase II activity induces cytopathic effects in permissive cells[J]. Proc Natl Acad Sci U S A, 2007, 104(30): 12482–12487. DOI: 10.1073/pnas.0705533104