Genetic and antigenic characterization of bovine viral diarrhea viruses isolated from cattle in Hokkaido,Japan

In our previous study, we genetically analyzed bovine viral diarrhea viruses (BVDVs) isolated from 2000 to 2006 in Japan and reported that subgenotype 1b viruses were predominant. In the present study, 766 BVDVs isolated from 2006 to 2014 in Hokkaido, Japan, were genetically analyzed to understand recent epidemics. Phylogenetic analysis based on nucleotide sequences of the 5′-untranslated region of viral genome revealed that 766 isolates were classified as genotype 1 (BVDV-1; 544 isolates) and genotype 2 (BVDV-2; 222). BVDV-1 isolates were further divided into BVDV-1a (93), 1b (371) and 1c (80) subgenotypes, and all BVDV-2 isolates were grouped into BVDV-2a subgenotype (222). Further comparative analysis was performed with BVDV-1a, 1b and 2a viruses isolated from 2001 to 2014. Phylogenetic analysis based on nucleotide sequences of the viral glycoprotein E2 gene, a major target of neutralizing antibodies, revealed that BVDV-1a, 1b and 2a isolates were further classified into several clusters. Cross-neutralization tests showed that BVDV-1b isolates were antigenically different from BVDV-1a isolates, and almost BVDV-1a, 1b and 2a isolates were antigenically similar among each subgenotype and each E2 cluster. Taken together, BVDV-1b viruses are still predominant, and BVDV-2a viruses have increased recently in Hokkaido, Japan. Field isolates of BVDV-1a, 1b and 2a show genetic diversity on the E2 gene with antigenic conservation among each subgenotype during the last 14 years.     

1株牛病毒性腹泻病毒分离毒株的基因组特征

旨在从宁夏某奶牛群持续感染牛分离牛源牛病毒性腹泻病毒（BVDV）,并解析其基因组特征,为研究我国不同地区BVDV分离株遗传演化规律提供理论依据。利用BVDV抗原检测试剂盒检测宁夏回族自治区银川市某示范区的240头高产奶牛间隔两周的双份抗凝血,筛选持续感染牛,分离血液淋巴细胞制备裂解液接种牛肾细胞（MDBK）,分离鉴定获得BVDV株,克隆测序获得全基因组序列,比较分析其遗传演化关系。从该示范区高产奶牛筛选获得2头持续感染牛,分离获得1株非致细胞病变型BVDV,命名为NX2019/01。测序获得基因组全序列（12 107 nt）,其中ORF长11 703 nt,编码3 898个氨基酸。在基因组水平,NX2019/01株与我国SD-15、ZM-95、XC、LN-1等1m亚型分离株相似性较高（92.17%~93.84%）,但E^rns、E1以及E2基因存在较大差异。示范区同群牛急性感染BVDV时,毒株E2蛋白N端编码区核苷酸突变可导致第9位或第67位氨基酸变异。重组分析表明,NX2019/01株E2基因179—288位核苷酸区段以及ZM-95株E1基因168位—E2基因332位核苷酸区段存在相似的重组信号,可能由主要亲本SD-15株与次要亲本LN-1株重组形成,表明NX2019/01株、ZM-95株在演化进程中与SD-15株以及LN-1株或早期流行的高度相似毒株存在密切关联。本研究从持续感染高产奶牛分离获得了牛源BVDV-1m亚型毒株,在基因组水平厘清了BVDV-1m亚型毒株的进化关系,并首次发现同亚型BVDV毒株基因同源重组,为进一步研究BVDV在我国的演化规律奠定了基础。     

Nucleotide sequence homology to bovine viral diarrhea virus 2 (BVDV 2) in the 5′ untranslated region of BVDVs from cattle with mucosal disease or persistent infection in Japan

Cytopathogenic and non-cytopathogenic bovine viral diarrhea viruses (BVDVs) were isolated from cattle with mucosal disease or persistent infection in Japan. These isolates were compared for antigenic properties by cross-neutralization tests with Japanese reference strains of BVDV belonging to classical type 1. Significantly low cross-reactivity to reference strains was noted, indicating the viruses to possibly represent a new serotype in Japan. Thus, to determine the genotype of the isolates, nucleotide sequences of the 5′ untranslated region were determined and compared with those of previously reported BVDV 1 and 2. The isolates were clearly shown to belong to BVDV 2, not to BVDV 1.     

Genetic and antigenic characterization of bovine viral diarrhoea virus type 2 isolated from cattle in India

Previous studies have shown that bovine viral diarrhoea virus type 1 (BVDV-1) subtype b is predominantly circulating in Indian cattle. During testing for exotic pestiviruses between 2007 and 2010, BVDV-2 was identified by real time RT-PCR in two of 1446 cattle blood samples originating from thirteen states of India. The genetic analysis of the isolated virus in 5′ UTR, N^pro, entire structural genes (C, E^rns, E1 and E2), nonstructural genes NS2-3 besides 3′ UTR demonstrated that the nucleotide and amino acid sequences showed highest similarity with BVDV-2. The entire 5′ and 3′ UTR consisted of 387 and 204 nucleotides, respectively, and an eight nucleotide repeat motif was found twice within the variable part of 3′ UTR that may be considered as a characteristic of BVDV-2. The phylogenetic analysis revealed that the cattle isolate and earlier reported goat BVDV-2 isolate fall into separate clades within BVDV-2a subtype. Antigenic typing with monoclonal antibodies verified the cattle isolate also as BVDV-2. In addition, cross-neutralization tests using antisera raised against Indian BVDV strains circulating in ruminants (cattle, sheep, goat and yak) displayed significant antigenic differences only between BVDV-1 and BVDV-2 strains. This is the first identification of BVDV-2 in Indian cattle that may have important implications for immunization strategies and molecular epidemiology of BVD.     

甘肃、青海和宁夏地区牛病毒性腹泻病毒Erns基因的变异分析

分析2013—2019年中国西北部分省区不同基因亚型牛病毒性腹泻病毒（BVDV）抗原基因E^rns的分子特征,了解其遗传演化规律。从甘肃、青海、宁夏规模化牛场送检的疑似牛病毒性腹泻发病牛150份EDTA抗凝血提取总RNA,利用RT-PCR扩增病毒基因组E^rns-E1区,克隆测序后比对,构建系统进化树进行遗传演化关系分析。利用牛肾细胞MDBK对检出的不同基因亚型BVDV进行分离,并鉴定其生物型。RT-PCR扩增结果表明,BVDV总体阳性率为37.33%,其中甘肃省、青海省、宁夏回族自治区BVDV阳性率分别为37.68%、35.71%、40.00%。获得56份E^rns-E1 DNA,克隆测序获得33条不同的E^rns序列,长度均为681 bp,分析表明流行株分属10个BVDV基因亚型：BVDV-1a （2株）、BVDV-1b （5株）、BVDV-1c （1株）、BVDV-1d （3株）、BVDV-1m （11株）、BVDV-1o （1株）、BVDV-1p （4株）、BVDV-1q （4株）、BVDV-1v （1株）、BVDV-2a （1株）。分离获得BVDV-1a亚型、BVDV-1b亚型、BVDV-1v亚型、BVDV-2a亚型分离株各1株,BVDV-1 d亚型分离株2株,均为非致细胞病变型。各亚型株间E^rns基因核苷酸相似性以BVDV-1a～1d经典亚型株（79.8%～85.9%）或1m～1q及1v新亚型株（81.0%～87.3%）较高,以BVDV-1 m和BVDV-1p流行株亚型间相似性最高（87.3%）。各亚型株E^rns基因编码蛋白的RNA酶活性位点以及双链RNA作用基序（₁₃₉KKGK₁₄₂）保守,但E^rns第26位糖基化位点（26 NRSL）在1m～1q、1v亚型株移位（24 NVSR）。首次以E^rns核苷酸序列构建系统进化树,结果显示1m～1q及1v等亚型BVDV株在进化上关系较为密切。本研究首次选用E^rns靶标基因对甘肃、青海、宁夏部分省区牛源BVDV株进行同源性及系统进化分析,发现10个基因亚型流行株,以1m亚型株最为普遍,1m～1q及1v等亚型株亲缘关系密切。     

Genotypic characteristics of bovine viral diarrhea virus 2 strains isolated in northern Italy

Two strains of Bovine viral diarrhea virus 2 (BVDV-2) were isolated from calves in northern Italy. Variations in the 5'-untranslated region (UTR) of the genome were studied by primary structure alignment and neighbor-joining method based phylogenetic tree analyses and by palindromic nucleotide substitutions at the three variable loci in the 5'-UTR. Genetic analysis indicated their appurtenance to genovar BVDV-2a. Nucleotide sequence at the 5'-UTR of strain BS-95-II, one of the Italian isolates from healthy calves, showed 98% homology to that of the Japanese isolate OY89, a cytopathic strain derived from cattle with mucosal disease.     

Genoepidemiological evaluation of Bovine viral diarrhea virus 2 species based on secondary structures in the 5' untranslated region

Bovine viral diarrhea virus 2 (BVDV-2) strains demonstrated in cattle, sheep, and adventitious contaminants of biological products have been evaluated by the palindromic nucleotide substitutions (PNS) method at the three variable loci (V1, V2 and V3) in the 5' untranslated region (UTR), to determine their taxonomical status. Variation in conserved genomic sequences was used as parameter for epidemiological evaluation of the species in relation with geographical distribution, animal host and virulence. Four genotypes, BVDV-2a, BVDV-2b, BVDV-2c, and BVDV-2d have been identified within the species. Taxonomical segregation corresponded to geographical distribution of genotype variants. Genotype 2a was present worldwide, and was the only circulating also in sheep, in addition to cattle. Genotypes 2b, 2c and 2d were restricted to South America. Contamination of biological products was related to genotypes 2a and 2d. Genetic variation could be related with chronological diffusion of the BVDV-2 species variants in different geographic areas. Chronologically, the species emerged in North America in 1978, spreading in UK and Japan, continental Europe, South America and New Zealand. Correlation between clinical features related with isolation of BVDV-2 strains and genetic variation indicated that subgenotype 1, variant 4 of genotype 2a was related with hemorrhagic syndrome. These observations suggest that evaluation of genomic secondary structure, by identifying markers for expression of virus biological activities and species evolutionary history, may be applied as useful tool for epidemiological evaluation of the BVDV-2 species, and possibly for other species of the genus Pestivirus.     

Genetic characterization of bovine viral diarrhea virus strains in Beijing,China and innate immune responses of peripheral blood mononuclear cells in persistently infected dairy cattle

To acquire epidemiological data on the bovine viral diarrhea virus (BVDV) and identify cattle persistently infected (PI) with this virus, 4,327 samples from Holstein dairy cows were screened over a four-year period in Beijing, China. Eighteen BVD viruses were isolated, 12 from PI cattle. Based on genetic analysis of their 5''-untranslated region (5''-UTR), the 18 isolates were assigned to subgenotype BVDV-1m, 1a, 1d, 1q, and 1b. To investigate the innate immune responses in the peripheral-blood mononuclear cells of PI cattle, the expression of Toll-like receptors (TLRs), RIG-I-like receptors, interferon-α (IFN-α), IFN-β, myxovirus (influenza virus) resistance 1 (MX1), and interferon stimulatory gene 15 (ISG15) was assessed by qPCR. When compared with healthy cattle, the expression of TLR-7, IFN-α, and IFN-β mRNA was downregulated, but the expression of MX1 and ISG-15 mRNA was upregulated in PI cattle. Immunoblotting analysis revealed that the expression of interferon regulatory factor 3 (IRF-3) and IRF-7 was lower in PI cattle than in healthy cattle. Thus, BVDV-1m and 1a are the predominant subgenotypes in the Beijing region, and the strains are highly divergent. Our findings also suggest that the TLR-7/IRF-7 signaling pathway plays a role in evasion of host restriction by BVDV.     

Genetic variety of bovine viral diarrhea virus 2 strains isolated from sheep

Bovine viral diarrhea virus 2 (BVDV-2) strains, isolated from sheep showing clinical symptoms of border disease, have been evaluated by the palindromic nucleotide substitution (PNS) method at the three variable loci (V1, V2 and V3) in the 5'-untranslated region (UTR) of genomic RNA. The characteristic two base-pairings common to the BVDV-2 species, a C-G pairing which was common to the V1 locus, and a G*U pairing common to the V2 locus, were observed in all tested strains. Strains BD-78 and C413 were identified by a unique C-G pairing at position 4 from the bottom of the V2 stem region, which is characteristic to BVDV-2b. BVDV-2d characteristic U-A pairing at position 18 of the V1 stem region was observed in five strains, Lees, 167 237, 168 149, 173 157 and 175 375. No strains have been assigned to the genotypes BVDV-2a or BVDV-2c. Furthermore, the investigation at the level of the 5'-UTR excluded the application in sheep of the proposed BVDV-2 genetic virulence markers described in cattle. The two specific positions of uracil and cytosine nucleotides related to low or high virulence where indifferently present in the ovine BVDV-2 strains responsible of border disease.     

Multiple diagnostic tests to identify cattle with Bovine viral diarrhea virus and duration of positive test results in persistently infected cattle

Several tests for Bovine viral diarrhea virus (BVDV) were applied to samples collected monthly from December 20, 2005, through November 27, 2006 (day 0 to day 342) from 12 persistently infected (PI) cattle with BVDV subtypes found in US cattle: BVDV-1a, BVDV-1b, and BVDV-2a. The samples included clotted blood for serum, nasal swabs, and fresh and formalin-fixed ear notches. The tests were as follows: titration of infectious virus in serum and nasal swabs; antigen-capture (AC) enzyme-linked immunosorbent assay (ELISA), or ACE, on serum, nasal swabs, and fresh ear notches; gel-based polymerase chain reaction (PCR) testing of serum, nasal swabs, and fresh ear notches; immunohistochemical (IHC) testing of formalin-fixed ear notches; and serologic testing for BVDV antibodies in serum. Of the 12 animals starting the study, 3 died with mucosal disease. The ACE and IHC tests on ear notches had positive results throughout the study, as did the ACE and PCR tests on serum. There was detectable virus in nasal swabs from all the cattle throughout the study except for a few samples that were toxic to cell cultures. The serum had a virus titer ≥ log₁₀ 1.60 in all samples from all the cattle except for 3 collections from 1 animal. Although there were several equivocal results, the PCR test most often had positive results. The BVDV antibodies were due to vaccination or exposure to heterologous strains and did not appear to interfere with any BVDV test. These findings illustrate that PI cattle may be identified by several tests, but differentiation of PI cattle from cattle with acute BVDV infection requires additional testing, especially of blood samples and nasal swabs positive on initial testing. Also, calves PI with BVDV are continual shedders of infectious virus, as shown by the infectivity of nasal swabs over the 11-mo study.     

Reactivity and prevalence of neutralizing antibodies against Japanese strains of bovine viral diarrhea virus subgenotypes

Bovine viral diarrhea virus (BVDV) field isolates show genetic and antigenic diversity. At least 14 subgenotypes of BVDV-1 and 4 of BVDV-2 have been identified in Artiodactyla worldwide. Of these, 6 subgenotypes of BVDV-1 and 1 of BVDV-2 have been isolated in Japan. Previously, we reported that each subgenotype virus expresses different antigenic characteristics. Here we investigated the reactivity of neutralizing antibodies against representative strains of Japanese BVDV subgenotypes using sera from 266 beef cattle to estimate the prevalence of this epidemic virus among cattle in Japan. Antibody titers at concentrations at least 4-fold higher than antibodies against other subgenotype viruses were considered subgenotype specific. Subgenotype-specific antibodies were detected from 117 (80.7%) of 145 sera samples (69.7% against BVDV-1a, 1.4% against BVDV-1b, 8.3% against BVDV-1c, and 1.4% against BVDV-2a). The results suggest that neutralization tests are useful in estimating currently epidemic subgenotypes of BVDV in the field.     

Genetic heterogeneity of bovine viral diarrhoea virus (BVDV) isolates from Turkey: identification of a new subgroup in BVDV-1

Genetic heterogeneity of Turkish ruminant pestiviruses was investigated by phylogenetic analysis of complete N(pro) encoding nucleotide sequences. A total of 30 virus isolates obtained from 15 provinces around the country between 1997 and 2005 were included in the phylogenetic analysis. Virus isolates mostly originated from cattle with one isolate from sheep. The bovine isolates all belonged to BVDV-1, the sheep isolate to BVDV-2. Fifteen isolates formed a new subgroup within BVDV-1, tentatively named BVDV-1l. The remaining bovine isolates were typed as BVDV-1a (n=4), BVDV-1b (n=4), BVDV-1d (n=3), BVDV-1f (n=2) and BVDV-1h (n=1). The isolates allocated to BVDV-1l originated from various geographical regions in different years. There was no correlation between genetic grouping and locations where isolates were obtained. Viruses originating from one farm in most cases belonged to the same subgroup (n=5). This study indicates that the newly detected subgroup BVDV-1l is predominant and widespread in Turkey. Moreover, an ovine virus isolate was identified as the first member of BVDV-2 reported in Turkey. A serological survey using samples from western Turkey indicated that BVDV-2 is also present in cattle.     

Variation in pestivirus growth in testicle primary cell culture is more dependent on the individual cell donor than cattle breed

The causes of bovine respiratory disease complex (BRDC) are multifactorial and include infection with both viral and bacterial pathogens. Host factors are also involved as different breeds of cattle appear to have different susceptibilities to BRDC. Infection with bovine pestiviruses, including bovine viral diarrhea virus 1 (BVDV1), BVDV2 and ‘HoBi’-like viruses, is linked to the development of BRDC. The aim of the present study was to compare the growth of different bovine pestiviruses in primary testicle cell cultures obtained from taurine, indicine and mixed taurine and indicine cattle breeds. Primary cells strains, derived from testicular tissue, were generated from three animals from each breed. Bovine pestivirus strains used were from BVDV-1a, BVDV-1b, BVDV-2a and ‘HoBi’-like virus. Growth was compared by determining virus titers after one passage in primary cells. All tests were run in triplicate. Virus titers were determined by endpoint dilution and RT-qPCR. Statistical analysis was performed using one way analysis of variance (ANOVA) followed by the Tukey’s Multiple Comparison Test (P?0.05). Significant differences in virus growth did not correlate with cattle breed. However, significant differences were observed between cells derived from different individuals regardless of breed. Variation in the replication of virus in primary cell strains may reflect a genetic predisposition that favors virus replication.     

Molecular characterization of <Emphasis Type="Italic">Bovine virus diarrhea viruses</Emphasis> species 2 (BVDV-2) from cattle in Turkey

Five BVDV species 2 (BVDV-2) isolates were detected from cattle in Turkey. Phylogenetic analysis was performed based on the 5′ untranslated regions (UTRs) and E2 coding gene regions, respectively. The isolates were closely related to BVDV-2a strains from North America and Canada used as references. This is the first report of the detection of BVDV-2 in naturally infected Turkish cattle. It is important to consider BVDV-2 for planning future BVDV control and vaccination programs in Turkey.     

The prevalent genotypes of bovine viral diarrhea virus in Japan, Germany and the United States of America

Genotypes and subgenotypes of bovine viral diarrhea virus (BVDV) field isolates from Japan, Germany and the United States of America (USA) were identified, and the prevalent pattern of BVDV in individual countries was estimated genetically. Subgenotypes were determined based on phylogenetic analyses of nucleotide sequences of a part of the E2-coding gene of BVDV. Forty-five, 61 and 56 BVDV strains were isolated from naturally infected cattle in Japan, Germany and USA, respectively, between 1980 and 2003. The most prevalent BVDV in these three countries was BVDV-1b. The second most prevalent BVDV strains were 1a, 1d and BVDV-2 in Japan, Germany and USA, respectively. The most prevalent subgenotype 1b in each country constructed individual small clusters in the subgenotype 1b branch in the phylogenetic tree. Although cattle and/or cattle products were moving among the three countries as part of international trade, the distribution of BVDV in the field in each country showed long-standing individual patterns.     

Genomic analyses of bovine viral diarrhea viruses isolated from cattle imported into Japan between 1991 and 2005

Thirty-one isolates of bovine viral diarrhea virus (BVDV) isolated within the past 15 years from imported cattle by the Japanese Animal Quarantine Service (AQS) were used in this study in which a 5'-untranslated region of each isolate was genetically analyzed. Twenty-six of the 31 isolates were classified as BVDV1 and the remainder as BVDV2. Phylogenetic analysis of the RT-PCR fragments amplified from the isolates showed the presence of viruses belonging to the BVDV1a, BVDV1b, BVDV1c, unclassified BVDV1 genotypes, and BVDV2. From the cattle of Australian origin, 16 of 17 isolates were classified as BVDV1c. This result was in agreement with a report showing that BVDV1c was a predominant subgenotype in Australia. From the cattle of North American origin, BVDV1 and BVDV2 species were both found. BVDV2 from the North American cattle was identified as the same cluster as the BVDV 890 strain, which is the prototype of BVDV2. These results suggest that the BVDVs isolated from exported cattle at the AQS reflect the predominant genotypes of BVDVs found in the exporting countries. The unclassified BVDV1 genotype of Chinese origin was in the same cluster as the ZM-95 strain, which was isolated from pigs in China. In this study, the genomic properties of 31 isolates of BVDV collected in the AQS were investigated. We concluded that isolates are genetically heterogeneous but geographically restricted. The information obtained from this report will be useful when carrying out epidemiological surveys of BVDV isolated in Japan.     

牛病毒性腹泻病毒BVDV-JL分离株全基因序列测定与分析

为了解中国农业科学院特产研究所分子生物学重点实验室前期分离的牛病毒性腹泻病毒BVDV-JL毒株完整的基因序列信息,本试验对BVDV-JL F6代毒株进行了完整基因序列测定。利用反转录—聚合酶链式反应(RT-PCR)方法分段扩增了BVDV-JL分离毒株的7段cDNA片段,分别克隆于pMD18-T载体并进行测序。BVDV-JL株基因组序列全长为12276 bp,编码3901个氨基酸。基因组两端为非编码区5'UTR和3'UTR。基因比对及进化树分析结果表明BVDV-JL与BVDV CP7同源性最高,核苷酸同源性为93.2%,归类于BVDV-1b2基因亚型。BVDV-JL是第1个在中国报道的BVDV-1b2亚型毒株。了解BVDV-JL完整基因序列有利于中国BVDV流行病学调查。