畜牧兽医学报  2024, Vol. 55 Issue (2): 649-659. DOI: 10.11843/j.issn.0366-6964.2024.02.022    PDF    
河北省廊坊市牛主要病毒性腹泻病原感染状况检测及牛冠状病毒演化分析
李思远1,2, 付新成5, 袁雪松2,3, 毛立2,3, 蔡旭航1,2, 孙心如2, 黄金2,3, 谢玲玲6, 王府6, 周华7, 张琪1, 李基棕2,3,4, 李彬2,3,4     
1. 西北农林科技大学动物医学院, 杨凌 712100;
2. 江苏省农业科学院兽医研究所, 农业农村部兽用生物制品工程技术重点实验室, 南京 210014;
3. 南京农业大学动物医学院, 南京 210095;
4. 江苏大学生命科学学院 食品与生物工程学院, 镇江 212013;
5. 河北省廊坊市农业农村局, 廊坊 065000;
6. 贵州省种畜禽种质测定中心, 贵阳 550018;
7. 黔西市动物疫病预防控制中心, 黔西 551500
摘要:牛病毒性腹泻病原严重危害畜牧业生产, BVDV、BoAstV、BCoV、MRV、BKoV、BoNeV、BNoV、BRV、BToV被认为是引起牛群腹泻的病原, 不仅如此, BVDV、BCoV和BToV除了能引起消化道症状外还能引起呼吸道症状, BVDV不仅能引起消化道和呼吸道症状, 还会导致母牛生殖功能异常, 这些疾病对畜牧业造成巨大的经济损失。为了解牛主要病毒性腹泻病原流行情况, 本次试验通过PCR检测方法对我国河北地区323份腹泻牛粪样品进行常见9种病毒性牛腹泻病原进行检测。结果显示, BVDV阳性检出率1.85% (6/323)、BCoV阳性检出率14.24% (46/323)、BRV阳性检出率57.89% (187/323)、BoAstV阳性检出率0.31% (1/323)、BoNeV阳性检出率10.84% (35/323)、BNoV阳性检出率4.33%(14/323)、MRV阳性检出率2.48% (8/323)、BToV阳性检出率4.64% (15/323)、BKoV阳性检出率11.76% (38/323)。混合感染共有76份, 占比23.53%(76/323)。HBLF2302株的全基因序列、SHENME基因进行同源性和遗传进化分析发现, HBLF2302与BCoV/NMG1/2022基因同源性最高, 为GIIb亚型。基于氨基酸序列比对和HBLF2302的S蛋白三级结构预测发现, 在S蛋白中157号位突变为157T, 854号位点突变为854I。318V为独特突变, 改变了与侧链残基的作用力, 与631C形成氢键连接。并发现中国株区别于其它地区的GⅡb亚型, 位于ORF1a区域有三个独特位点。本研究丰富了牛腹泻病原的流行病学数据, 为牛腹泻病的防控奠定基础。
关键词病毒性腹泻    牛冠状病毒    流行病学调查    腹泻病原    
Detection of Bovine Viral Diarrhea Pathogens and Evolution Analysis of Bovine Coronavirus in Langfang, Hebei
LI Siyuan1,2, FU Xincheng5, YUAN Xuesong2,3, MAO Li2,3, CAI Xuhang1,2, SUN Xinru2, HUANG Jin2,3, XIE Lingling6, WANG Fu6, ZHOU Hua7, ZHANG Qi1, LI Jizong2,3,4, LI Bin2,3,4     
1. College of Veterinary Medicine, Northwest A&F University, Yangling 712100, China;
2. Institute of Veterinary Medicine, Jiangsu Academy of Agricultural Sciences, Key Laboratory of Veterinary Biological Engineering and Technology, Ministry of Agriculture and Rural Affairs, Nanjing 210014, China;
3. College of Veterinary Medicine, Nanjing Agricultural University, Nanjing 210095, China;
4. School of Food and Biological Engineering, School of Life Sciences, Jiangsu University, Zhenjiang 212013, China;
5. Langfang Municipal Bureau of Agriculture and Rural Affairs, Langfang 065000, China;
6. Guizhou Testing Center for Livestock and Poultry Germplasm, Guiyang 550018, China;
7. Qianxi Animal Disease Control Center, Qianxi 551500, China
Abstract: Cattle viral diarrhea pathogens have adverse effects on animal husbandry. BVDV, BoAstV, BCoV, MRV, BKoV, BoNeV, BNoV, BRV, BToV were well known as the cause of diarrhea. On the other hand, BVDV, BCoV and BToV also caused respiratory system symptoms. Among them, BVDV not only causes alimentary tract symptoms and respiratory symptoms, but also leads to abnormal reproductive function of cows causing great economic losses. In order to understand the prevalence of the main viral diarrhea pathogens in cattle, 323 fecal samples from diarrheal cow in Hebei province of China were detected by PCR to detect nine common viral diarrhea pathogens. Prevalence of 1.85% (6/323), 14.24% (46/323), 57.89% (187/323), 0.31 (1/323), 10.84% (35/323), 4.33% (14/323), 2.48% (8/323), 4.64% (15/323), 11.76% (38/323), was determined for BVDV, BCoV, BRV, BoAstV, BoNeV, BNoV, MRV, BToV, BKoV, respectively. There were 76 co-infection samples accounting for 23.53% (76/323). Homology and genetic evolution analysis of the complete genomic sequence, S, HE, N, M and E genes of HBLF2302 strain showed that HBLF2302 gene had high homology with BCoV/NMG1/2022 and was GIIb subtype. Based on the amino acid sequence alignment and the tertiary structure prediction of HBLF2302 S protein, it was found that the 157 site was mutated to 157T and the 854 site was mutated to 854I in the S protein. 318V is a unique mutation, which changes the force with the side chain residue and forms a hydrogen bond connection with 631C. It was found that the Chinese strain was different from the GⅡb subtype in other regions, and there were three unique sites in the ORF1a region. This study enriched the epidemiological data about bovine diarrhea pathogens for laying the foundation of bovine pathogens prevention and control.
Key words: viral diarrhea    bovine coronavirus    epidemiology investigation    diarrhea pathogens    

牛冠状病毒(bovine coronavirus, BCoV)属于β冠状病毒,在全球范围广泛流行。主要编码5个结构蛋白,血凝素酯酶(hemagglutinin-esterase,HE)、刺突蛋白(spike protein,S)、包膜蛋白(envelope protein,E)、膜蛋白(membrane,M)和核衣壳蛋白(nucleocapsid,N)。BCoV最初被认为是肠道病毒, 可以造成新生犊牛腹泻[1-2]、牛冬季痢疾,但迄今为止,越来越多的临床案例表明该病毒可以从呼吸道检测出并导致各年龄段牛呼吸道症状[3-4],因此BCoV也是呼吸道病原的说法越来越频繁地被提及。BCoV具有很大的跨物种传播潜力,CoV-OC43被认为与BCoV相近,这两种病毒均能以9-O乙酰唾液酸作为糖结合底物[5-6],受人白细胞抗原Ⅰ亚群(human leukocyte antigen,HLA-Ⅰ)影响[5, 7-8],所以BCoV常被作于人呼吸道冠状病毒(人冠状病毒OC-43、SARS-CoV2、HKU1)的研究模型。除了感染牛以外,在其他动物(如羊驼、羊、长颈鹿等)中也检出类BCoV,其中更包括了德国的人呼吸道株4403[9-12], 具有溢出跨宿主传播风险。

牛病毒性腹泻/黏膜病病毒(bovine viral diarrhea virus, BVDV)是全世界牛常见的疾病之一。BVDV可以导致牛生育能力下降、产奶量下降、呼吸道症状、腹泻、流产和产死胎等生殖功能异常[13-14]。在亚洲不同国家报道中,牛星状病毒(bovine astroviruses, BoAstV)被证明与牛腹泻有关,但报道的发病率不一致[15-16]。牛轮状病毒(bovine rotavirus virus, BRV)是牛腹泻最常见的原因,能导致显著的并发症和死亡率,在新生犊牛中尤为明显[17-18]。牛环曲病毒(bovine torovirus, BToV)主要引起小牛和成年牛腹泻,还具有双重组织嗜性,可以在呼吸道中检出,并导致牛腹泻,严重患病牛可至死亡[19-21]。牛诺瓦病毒(bovine norovirus, BoNoV)是引起动物胃肠炎的重要病原[22]。嵌杯病毒科的牛纽布病毒(bovine nebovirus, BoNeV)可导致患病牛小肠绒毛萎缩,引起腹泻[22-23]。哺乳动物正呼肠孤病毒(mammalian orthoreovirus, MRV)和牛嵴病毒(bovine kobuvirus, BKoV)也被认为是引起牛腹泻的病原之一[24-25]

本研究对我国河北省廊坊市以BCoV为主的腹泻病原进行流行病学调查,对BCoV遗传进化分析,为该地这9种病的发病、流行情况及疫病防控提供科学数据。

1 材料与方法 1.1 样品

从我国河北廊坊地区14个养殖场中采集了323份腹泻奶牛粪样。

1.2 主要试剂

RNA提取试剂盒、逆转录酶、Green Taq Mix、2×Phanta Max Master Mix (Dye Plus) 均购自南京诺唯赞生物科技有限公司,琼脂糖购自北京擎科生物科技有限公司,胶回收试剂盒购自Omega Bio-Tek公司, 磷酸盐缓冲液(PBS)购自塞维尔生物科技有限公司。

1.3 主要仪器

PCR扩增仪和离心机:德国Eppendorf公司;凝胶成像系统:上海天能公司。

1.4 样品处理

将采集的粪样每份转移部分至2 mL离心管中,其余冻存至-80 ℃冰箱。每个EP管中加入1 mL经高压灭菌的PBS溶液进行震荡充分混匀,离心后弃去沉淀取上清用0.22 mm滤器过滤后流液备用。

1.5 引物合成

参考文献[26-32]合成BCoV、BVDV、BoAstV、BRV、BToV、BNoV、BoNeV、MRV、BKoV相关片段的引物(表 12)。

表 1 病毒性腹泻病原检测引物信息 Table 1 The sequence of primers used in this study
表 2 BCoV SHE引物信息 Table 2 The sequence of primers to detect S and HE gene
1.6 病毒RNA提取及转录

参照说明进行提取总RNA后进行反转录,反转录产物cDNA冻存于-80 ℃冰箱备用。

1.7 RT-PCR扩增

根据表 1检测引物,采用20 μL体系:10 μL Mix,0.4 μL上游引物F,0.4 μL下游引物R,1 μL cDNA,8.2 μL ddH2O。使用降落PCR方法扩增片段,反应程序:预变性95 ℃ 5 min;95 ℃ 15 s,48~55 ℃ 15 s,72 ℃ 30 s,共35个循环;彻底延伸72 ℃ 10 min。PCR产物经1%琼脂糖凝胶电泳,观察结果。

根据SHE特异性引物扩增BCoV阳性样品,采用25 μL体系:12.5 μL Mix,1 μL上游引物F,1 μL下游引物R,2 μL cDNA,8.5 μL ddH2O。使用降落PCR方法扩增目的片段,程序为:预变性95 ℃ 5 min;95 ℃ 15 s,48~55 ℃ 15 s,72 ℃ 1 min共35个循环;彻底延伸72 ℃ 10 min。得到产物使用1%琼脂糖凝胶电泳,观察结果。

1.8 胶回收验证及进化树构建

按胶回收试剂盒说明书回收目的条带胶块,由上海生工生物工程股份有限公司测序,通过测序结果进行序列比对。同时将提取的RNA送上海探普生物公司进行病毒基因组下一代测序。

使用MEGA 11软件进行BCoV的HESNME、ORF1a和全序列基因核酸序列比对,并使用最大似然法构建进化树。

2 结果 2.1 RT-PCR检测

在检测323份样品中BCoV总体阳性检出率为14.24% (46/323),BVDV为1.85% (6/323),BRV为57.89% (187/323),BoAstV为0.31% (1/323),BoNeV为10.84% (35/323),BoNoV为4.33% (14/323),MRV为2.48% (8/323),BToV为4.64% (15/323),BKoV为11.76% (38/323)。该地区BCoV检出率并不高,而BRV高达57.89%,检出率最低的是BoAstV,仅有0.31%。各个牛场不同病原阳性检出率差异较大,BVDV检出率为0.00%~17.65%,BRV检出率为0.00%~100%,BCoV检出率为0.00%~45.45%,BKoV检出率为0.00%~42.47%,BNeV检出率为0.00%~62.50%,BNoV检出率为0.00%~30.00%,BToV检出率为0.00%~13.33%,MRV检出率为0.00%~11.76%,BoAstV检出率为0.00%~3.33%。其中BRV在每个牛场中均能检出,在7个奶牛场中检出率为100.00%。而BCoV也在12个奶牛场中检出,仅次于BRV。

2.2 混合感染检测

总共有22种混合感染类型,最多同时感染5种病毒。以BRV和BoNeV混合感染检出率最高,高达28.95% (22/76),BCoV混合感染比率为36.84%(28/76) (表 3)。

表 3 混合感染状况分析结果 Table 3 Analysis result of co-infection status
2.3 BCoV的HES基因序列扩增

阳性样本采用特异性引物扩增HES片段,可分段扩增出1株HES全长序列,命名为HBLF2302株(图 1)。

M.DL2000 DNA相对分子质量标准; 1. HBLF2302检测片段;2. HA;3. HB;4. HE;5. S1;6. 35;7. 3536;8. 36;9. 3637;10. S2;11. S3;12. DL5000 DNA相对分子质量标准 M.DL2000 DNA maker; 1. HBLF2302 detected segment; 2. HA; 3. HB; 4. HE; 5. S1;6. 35;7. 3536;8. 36;9. 3637;10. S2;11. S3;12. DL5000 DNA maker 图 1 BCoV HES基因扩增结果 Fig. 1 Amplification consequence of HE and S genes of BCoV
2.4 BCoV结构蛋白遗传进化分析

根据NCBI数据库及DNAstar 7.0中的Megalign比对,HBLF2302 S基因与我国华东地区MZ711 357.1株相似性最高,达到99.5%。HBLF2302 HE基因,与我国西南地区获得的MN982 191.1序列相似性最高,可达99.53%,序列中无12个碱基插入。HBLF2302 EMN基因与我国其它地区测得的序列相似性均在99.7%以上,差异较小。

从NCBI数据库获得的部分序列S序列进行比对,S进化树结果表明HBLF2302属于GⅡb基因型,为国内BCoV流行的基因型,值得注意的是,HY24株为我国牦牛源毒株,p95为我国台湾省于HRT-18细胞上测得,属于GⅠa基因型。不仅如此BCoV-GX-NN190313和BJ232犬冠状病毒处于独立分支,提示国内BCoV可能存在新基因型。韩国流行株主要以GⅡa基因型为主,但也出现有GⅡb基因型。表明不同基因型跨地区传播的风险在增加。可能存在新基因型的出现,有待进一步测序验证。HEEMN基因在不同基因型中未见有明显分化,表明在不同基因型中除了S差异较大具有分型价值外,其余四个结构蛋白可能不可作为分型依据。

2.5 BCoV S蛋白三级结构分析

通过SWISS-MODEL同源建模及PyMDL分析,发现氨基酸157号区别于其他参考序列157V和157I,854I较为少见,318号位为独特突变位点,在HBLF2302株中处于选择。其中157 T和318V位于S1,854I位于S2。318号位点在HBLF2302上由318I突变为318V时,与631C形成氢键(图 2)。

A. 318I局部结构;B. 318V局部结构 A. Local structure of 318I of S; B. Local stucture of 318V of S 图 2 BCoV S蛋白三级结构图 Fig. 2 Protein tertiary structure diagram of S protein from BCoV
2.6 BCoV ORF1a序列及分析

通过氨基酸比对结果表明,在目前上传至NCBI中的BCoV中国株的开放阅读框1a(Open reading frame ORF1a)基因中834、857和889号位为区别于其它地区的GⅡb株独特的位点。在其他序列中834W、857 L和889R号位均为一致,而在中国株中834W为834R,857 L为857 S,而在889R为终止密码子(图 3)。从遗传进化树上可以看出,除了牦牛源MH810 163.1 HY24仍与GⅠa型传统株处于同一分支外,中国株独立于同为GⅡb型的日本株和美洲株独立成簇(图 4)。由此可见,在ORF1a上,同一基因型的BCoV毒株存在一定差异性,并且中国株呈现出独特特征。

图 3 ORF1a 834、857和889号位氨基酸序列比对结果 Fig. 3 Amino acid sequence alignment of ORF1a 834, 857 and 889
图 4 基于BCoV ORF1a的进化树分析 Fig. 4 Phylogenetic tree analysis based on ORF1a of BCoV
2.7 BCoV全基因遗传进化分析

根据同源性分析结果,HBLF2302与从内蒙古地区获得的序列OP924 545.1相似性最高,在进化树上也与OP924 545.1进化距离接近。全基因组序列FJ938 067.1 4408株与基于S基因的进化树具有差异,在基于S基因进化树中与GⅠb基因型归为一簇,而在全基因组序列中与GⅠa基因型归为一簇(图 5)。

图 5 基于全基因组的BCoV遗传演化分析 Fig. 5 Genetic evolution analysis of BCoV based on whole genome
3 讨论

牛群腹泻的病原多样,呈现出多重病原混合感染的复杂情况,在临床中难以证实某种病原与腹泻的直接相关性。本研究表明,廊坊地区BCoV检出率为14.24%,其他9种腹泻病原的检出率为0.52%~57.89%。混合感染情况占35.37%。近年来,BCoV流行率有着上升趋势,不同地区流行率各有差异,可能存在地方性流行的特点[33-34]。在吉林部分地区报道牛冠状病毒阳性检出率为21.10%[31],云南部分地区为14.29%[32],在北疆地区检出率为0%~38.89%[33],在四川部分地区为16.80%[34],而在高原地区,腹泻牛样本检出率高达69.05%,患呼吸道疾病的牛群中检出率为72.45%[35]。由此可见,在不同地区不同样本牛冠状病毒发病率存在差异,且在腹泻牛群中检出率较高。不仅如此,在国外报道中,巴西为68.7%[36],日本为57.7%[37],韩国为39.8%[38],土耳其为21.9%[39],伊朗为7.2%[40],在不同国家腹泻牛群中检出率也有所差异。本文中均为腹泻样本,牛冠状病毒的检出率为14.24%,可能是由于存在多种病原微生物混合感染,非单一因素或病原引起的牛群腹泻或存在地区流行趋势,因此无法确切表明,牛冠状病毒与腹泻具有直接关联。

按照目前BCoV的基因型划分,不同地区存在着不同基因型的BCoV,美洲、东南亚和东亚流行株(除韩国外)主要为GⅡb,韩国主要为GⅡa,欧洲主要为GⅠb。此外还存少数以经典株为主的为GⅠa和GⅡc基因型。但由于对牛冠状病毒基因型划分依据仍不统一,按目前划分S基因与全基因序列存在差异。更值得注意的是BCoV全基因组序列差异和变异较小,不同地区流行的BCoV基因型主要可能是由贸易原因为主导,因此即使在地处亚洲的土耳其,在S基因上趋向于欧洲的GⅠb基因型。

不仅如此,关键的蛋白的突变和重组可能导致BCoV的易感性增加。S中的S1具有双受体结合基序[37],S1的N端结构域与9-O乙酰唾液酸结合[6],S1的受体结合结构域(receptor-binding domain,RBD)被认为与受体结合,并且大量选择位点位于RBD处,可能改变其宿主谱进而导致溢出。本研究获取S基因的157号位点发生突变,该突变区别于157V和157I,318V为独特突变位点, 位于S的NTD端,318号位点的突变改变了侧链残基间相互作用力,是否改变糖嗜性仍需进一步研究。除此之外,本文通过序列比对,发现在参考序列的中国株及HBLF2302株中的ORF1a上的834、857及889 aa存在独特突变位点,区别于其他地区的GⅡb基因型毒株并在889号位点出现终止密码子,是否在后续报道毒株中存在普遍性用于更加细化解析牛冠状病毒,仍需要大量测得序列来验证,并进一步分析。

本研究对中国河北地区牛冠状病毒及其他病毒性腹泻病原进行了流行病学调查,丰富了牛冠状病毒及其他常见病毒性腹泻病原的流调数据,为进一步研究牛冠状病毒及牛腹泻性病原提供了数据。

4 结论

廊坊地区牛冠状病毒检出率为14.24%与其他不同地区存在差异,检出率较低。测得的HBLF2302属于GⅡb亚型流行株。在ORF1a发现区别于中国以外地区GⅡb型独特的位点,并在S蛋白NTD端上有独特突变位点。

参考文献
[1]
CLARK M A. Bovine coronavirus[J]. Br Vet J, 1993, 149(1): 51-70. DOI:10.1016/S0007-1935(05)80210-6
[2]
TSUNEMITSU H, SAIF L J. Antigenic and biological comparisons of bovine coronaviruses derived from neonatal calf diarrhea and winter dysentery of adult cattle[J]. Arch Virol, 1995, 140(7): 1303-1311. DOI:10.1007/BF01322757
[3]
RAHE M C, MAGSTADT D R, GROELTZ-THRUSH J, et al. Bovine coronavirus in the lower respiratory tract of cattle with respiratory disease[J]. J Vet Diagn Invest, 2022, 34(3): 482-488. DOI:10.1177/10406387221078583
[4]
SOULES K R, RAHE M C, PURTLE L, et al. Bovine coronavirus infects the respiratory tract of cattle challenged intranasally[J]. Front Vet Sci, 2022, 9: 878240. DOI:10.3389/fvets.2022.878240
[5]
COLLINS A R. Human coronavirus OC43 interacts with major histocompatibility complex class I molecules at the cell surface to establish infection[J]. Immunol Invest, 1994, 23(4-5): 313-321. DOI:10.3109/08820139409066826
[6]
PENG G Q, XU L Q, LIN Y L, et al. Crystal structure of bovine coronavirus spike protein lectin domain[J]. J Biol Chem, 2012, 287(50): 41931-41938. DOI:10.1074/jbc.M112.418210
[7]
COLLINS A R. HLA class I antigen serves as a receptor for human coronavirus OC43[J]. Immunol Invest, 1993, 22(2): 95-103. DOI:10.3109/08820139309063393
[8]
SZCZEPANSKI A, OWCZAREK K, BZOWSKA M, et al. Canine respiratory coronavirus, bovine coronavirus, and human coronavirus OC43:receptors and attachment factors[J]. Viruses, 2019, 11(4): 328. DOI:10.3390/v11040328
[9]
ERLES K, TOOMEY C, BROOKS H W, et al. Detection of a group 2 coronavirus in dogs with canine infectious respiratory disease[J]. Virology, 2003, 310(2): 216-223. DOI:10.1016/S0042-6822(03)00160-0
[10]
VLASOVA A N, SAIF L J. Bovine coronavirus and the associated diseases[J]. Front Vet Sci, 2021, 8: 643220. DOI:10.3389/fvets.2021.643220
[11]
ZHANG X M, HERBST W, KOUSOULAS K G, et al. Biological and genetic characterization of a hemagglutinating coronavirus isolated from a diarrhoeic child[J]. J Med Virol, 1994, 44(2): 152-161. DOI:10.1002/jmv.1890440207
[12]
AMER H M. Bovine-like coronaviruses in domestic and wild ruminants[J]. Anim Health Res Rev, 2018, 19(2): 113-124. DOI:10.1017/S1466252318000117
[13]
DEREGT D, TESSARO S V, BAXI M K, et al. Isolation of bovine viral diarrhoea viruses from bison[J]. Vet Record, 2005, 157(15): 448-450. DOI:10.1136/vr.157.15.448
[14]
BROCK K V. The many faces of bovine viral diarrhea virus[J]. Vet Clin North Am: Food Anim Pract, 2004, 20(1): 1-3. DOI:10.1016/j.cvfa.2003.12.002
[15]
CAFFARENA R D, CASAUX M L, SCHILD C O, et al. Causes of neonatal calf diarrhea and mortality in pasture-based dairy herds in Uruguay: a farm-matched case-control study[J]. Braz J Microbiol, 2021, 52(2): 977-988. DOI:10.1007/s42770-021-00440-3
[16]
ZHU J, QI M P, JIANG C W, et al. Prevalence of bovine astroviruses and their genotypes in sampled Chinese calves with and without diarrhoea[J]. J Gen Virol, 2021, 102(8): 001640.
[17]
PAPP H, LÁSZLÓ B, JAKAB F, et al. Review of group A rotavirus strains reported in swine and cattle[J]. Vet Microbiol, 2013, 165(3-4): 190-199. DOI:10.1016/j.vetmic.2013.03.020
[18]
BÁNYAI K, ESONA M D, MIJATOVIC S, et al. Zoonotic bovine rotavirus strain in a diarrheic child, Nicaragua[J]. J Clin Virol, 2009, 46(4): 391-393. DOI:10.1016/j.jcv.2009.08.005
[19]
DRAKER R, ROPER R L, PETRIC M, et al. The complete sequence of the bovine torovirus genome[J]. Virus Res, 2006, 115(1): 56-68. DOI:10.1016/j.virusres.2005.07.005
[20]
HOET A E, CHO K O, CHANG K O, et al. Enteric and nasal shedding of bovine torovirus (Breda virus) in feedlot cattle[J]. Am J Vet Res, 2002, 63(3): 342-348. DOI:10.2460/ajvr.2002.63.342
[21]
ITO T, OKADA N, OKAWA M, et al. Detection and characterization of bovine torovirus from the respiratory tract in Japanese cattle[J]. Vet Microbiol, 2009, 136(3-4): 366-371. DOI:10.1016/j.vetmic.2008.11.014
[22]
SMILEY J R, CHANG K O, HAYES J, et al. Characterization of an enteropathogenic bovine calicivirus representing a potentially new calicivirus genus[J]. J Virol, 2002, 76(20): 10089-10098. DOI:10.1128/JVI.76.20.10089-10098.2002
[23]
CHO E H, SOLIMAN M, ALFAJARO M M, et al. Bovine nebovirus interacts with a wide spectrum of histo-blood group antigens[J]. J Virol, 2018, 92(9): e02160-17.
[24]
LI H P, TANG C, YUE H. Molecular detection and genomic characteristics of bovine kobuvirus from dairy calves in China[J]. Infect Genet Evol, 2019, 74: 103939. DOI:10.1016/j.meegid.2019.103939
[25]
WANG L Y, FREDRICKSON R, DUNCAN M, et al. Bovine kobuvirus in calves with diarrhea, United States[J]. Emerg Infect Dis, 2020, 26(1): 176-178. DOI:10.3201/eid2601.191227
[26]
WU Q, LI J Z, WANG W, et al. Next-generation sequencing reveals four novel viruses associated with calf diarrhea[J]. Viruses, 2021, 13(10): 1907. DOI:10.3390/v13101907
[27]
LI H, ZHANG B, YUE H, et al. First detection and genomic characteristics of bovine torovirus in dairy calves in China[J]. Arch Virol, 2020, 165(7): 1577-1583. DOI:10.1007/s00705-020-04657-9
[28]
GUO Z J, HE Q F, YUE H, et al. First detection of Nebovirus and Norovirus from cattle in China[J]. Arch Virol, 2018, 163(2): 475-478. DOI:10.1007/s00705-017-3616-6
[29]
吕珽, 陈虹吟, 汤承, 等. 牦牛源正呼肠孤病毒2型的检测和分离鉴定[J]. 畜牧兽医学报, 2021, 52(8): 2361-2368.
LV T, CHEN H Y, TANG C, et al. Detection, isolation and identification of mammalian orthoreovirus serotype 2 in Yaks[J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(8): 2361-2368. (in Chinese)
[30]
REUTER G, BOLDIZSÁR Á, PANKOVICS P. Complete nucleotide and amino acid sequences and genetic organization of porcine kobuvirus, a member of a new species in the genus Kobuvirus, family Picornaviridae[J]. Arch Virol, 2009, 154(1): 101-108. DOI:10.1007/s00705-008-0288-2
[31]
孙飞雁, 叶京飞, 魏宇, 等. 吉林省肉牛冠状病毒感染状况调查及分析[J]. 畜牧兽医学报, 2023, 54(2): 673-82.
SUN F Y, YE J F, WEI Y, et al. Epidemiological investigation and analysis of bovine coronavirus in beef cattle in Jilin province[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(2): 673-682. (in Chinese)
[32]
张松梅, 王鑫文, 代琦, 等. 云南省部分地区牛冠状病毒检测及其N基因序列分析[J]. 中国动物检疫, 2022, 39(11): 35-40.
ZHANG S M, WANG X W, DAI Q, et al. Detection of bovine coronavirus and analysis on its N gene sequence in some areas of Yunnan province[J]. China Animal Health Inspection, 2022, 39(11): 35-40. (in Chinese)
[33]
张坤, 剡根强, 王静梅, 等. 新疆北疆部分地区致犊牛腹泻冠状病毒的分子流行病学调查[J]. 中国兽医科学, 2015, 45(12): 1270-1276.
ZHANG K, YAN G Q, WANG J M, et al. Molecular epidemiological investigation of calf diarrhea coronavirus in parts of northern Xinjiang[J]. Chinese Veterinary Science, 2015, 45(12): 1270-1276. (in Chinese)
[34]
邓飞, 陈弟诗, 陈斌, 等. 2020年四川省部分地区腹泻牛群牛冠状病毒感染情况调查[J]. 中国动物检疫, 2021, 38(9): 13-16.
DENG F, CHEN D S, CHEN B, et al. Investigation on infection status of bovine coronavirus in diarrheal cattle in some regions of Sichuan province in 2020[J]. China Animal Health Inspection, 2021, 38(9): 13-16. (in Chinese)
[35]
何琪富. 牦牛冠状病毒的分子检测及基因组研究[D]. 成都: 西南民族大学, 2019.
HE Q F. Molecular detection and genomic characterization of bovine coronavirus in yak[D]. Chengdu: Southwest Minzu University, 2019. (in Chinese)
[36]
COURA F M, FREITAS M D, RIBEIRO J, et al. Longitudinal study of Salmonella spp., diarrheagenic Escherichia coli, Rotavirus, and Coronavirus isolated from healthy and diarrheic calves in a Brazilian dairy herd[J]. Trop Anim Health Prod, 2015, 47(1): 3-11. DOI:10.1007/s11250-014-0675-5
[37]
MAWATARI T, HIRANO K, IKEDA H, et al. Surveillance of diarrhea-causing pathogens in dairy and beef cows in Yamagata Prefecture, Japan from 2002 to 2011[J]. Microbiol Immunol, 2014, 58(9): 530-535. DOI:10.1111/1348-0421.12174
[38]
PARK J, HAN D G, KIM S, et al. Prevalence of coronavirus from diarrheic calves in the Republic of Korea[J]. Asian Pac J Trop Biomed, 2018, 8(1): 1-6. DOI:10.4103/2221-1691.221037
[39]
SEVINC TEMIZKAN S, ALKAN F. Bovine coronavirus infections in Turkey: molecular analysis of the full-length spike gene sequences of viruses from digestive and respiratory infections[J]. Arch Virol, 2021, 166(9): 2461-2468. DOI:10.1007/s00705-021-05147-2
[40]
LOTFOLLAHZADEH S, MADADGAR O, MOHEBBI M R, et al. Bovine coronavirus in neonatal calf diarrhoea in Iran[J]. Vet Med Sci, 2020, 6(4): 686-694. DOI:10.1002/vms3.277

(编辑   白永平)