Isolation and molecular characterization of bovine herpesvirus 4 from cattle in mainland China.

Bovine herpesvirus 4 (BoHV-4) is one of the most important of the known viral respiratory and reproductive pathogens of both young and adult cattle. However, BoHV-4 has not been isolated or detected in mainland China prior to this study. In 2019, BoHV-4 strain 512 was isolated from cattle in Heilongjiang Province, China, using MDBK cells, and characterized by PCR, nucleotide sequence analysis, and transmission electron microscopy. Two other unknown herpesvirus strains, BL6010 and J4034, which were isolated from cattle in 2009 in China and stored at -70℃, were also propagated in MDBK cells and identified as BoHV-4 by PCR. Phylogenetic analysis based on partial nucleotide sequences of the thymidine kinase (TK) gene and glycoprotein B (gB) gene for the three isolates indicated that these three Chinese strains belong to BoHV-4 genotype 1. A preliminary virus neutralization test revealed that 64% of the 70 bovine sera (45/70) collected from Inner Mongolia Autonomous Region, China, had anti-BoHV-4 antibodies and that natural BoHV-4 infection occurred in cattle in China. Here, we report for the first time the isolation and molecular characterization of BoHV-4 from cattle in mainland China.

Introduction

Bovine herpesvirus 4 (BoHV-4) is a member of the family Herpesviridae, subfamily Gammaherpesvirinae, and genus Rhadinovirus. BoHV-4 has been isolated from cattle with a variety of clinical signs, including respiratory disease [1, 2], postpartum metritis [3-6], abortion [7-12], mastitis [13-15], and vulvovaginitis [16]. It has also been isolated from apparently healthy animals [17]. Recently, BoHV-4 has also been detected in dairy cows with subclinical endoenteritis [6, 18]. BoHV-4 has no close biological or virological relationship to other known herpesviruses infecting members of the family Bovidae [19]. Cattle are the natural hosts of this virus, but African buffaloes are also susceptible to BoHV-4 infection [20]. It has also been isolated from calf kidney cell cultures [21], bovine arterial endothelial cell cultures [22], and from various hosts, including American bison [23], sheep [24], and cats [25, 26]. Goats, cats, guinea pigs, and rabbits can be infected experimentally [27-29].

The role of BoHV-4 in infections of the respiratory and genital tract has been studied by several research groups. The virus has been reported to be responsible for post-partum and chronic metritis, alone or in combination with other pathogens [3]. Viral genomic DNA and antigens of BoHV-4 have been detected in the lactiferous duct and sinus epithelium of mammary tissues collected from a cow with clinical mastitis [30]. BoHV-4 has also been isolated from three vaginal discharge samples from cows with post-partum metritis, some of which had BoHV-4-specific neutralizing antibodies [4]. Ex vivo explant models have shown that BoHV-4 can replicate and invade bovine genital mucosae and respiratory epithelial cells [31, 32]. Recently, a survey showed that cows with post-partum metritis shed BoHV-4 through several bodily secretions, including vaginal secretions, nasal exudates and milk, indicating that horizontal transmission may occur in BoHV-4-positive herds [6]. Like other herpesviruses, BoHV-4 has been isolated from persistently infected cattle [33]. Latency in lymphoid tissues and prolonged viremia associated with mononuclear cells has been demonstrated [34]. The latent virus can be reactivated by stress factors or glucocorticoid treatment [33].

The BoHV-4 genome consists of a linear dsDNA of approximately 145 kbp that contains a unique central sequence of 110 kbp with 41.4% G-C content, flanked at both ends by tandem repeats with higher G-C content, a noncoding region called polyrepetitive DNA (prDNA) [35, 36]. The central part of the genome is well conserved among BoHV-4 strains, but the prDNA varies in size depending on the number of repetitions of a 200-bp sequence [35]. All BoHV-4 strains analyzed so far exhibit, similar restriction profiles and a well-conserved central part of the genome [37]. BoHV-4 has been divided into two groups based on the restriction pattern. One is the Movar 33/63-like (Movar-like) group and the other is the DN-599-like group. Most of the Movar 33/63-like viruses were isolated in Europe, whereas most of the DN-599-like viruses were isolated in North America [38]. In 2012, analysis of BoHV-4 strains isolated from cervico-vaginal mucus of aborted cows in Argentina revealed a high genetic divergence among strains, which could be classified in three different groups: genotype 1 comprises Movar-like strains (European prototype), genotype 2 includes DN-599-like strains (American prototype), and the novel genotype 3 includes MGA1075 and 07 435 strains [37].

In 2019, we used the non-cytopathic-effect (NCPE) strain 3877 of bovine viral diarrhea virus (BVDV) to inoculate two four-month old calves to establish an animal infection model of BVDV for challenge infection. Nasal swabs and buffy coat fractions were collected and used to inoculate cultured MDBK cells for virus isolation. A cytopathic effect (CPE) similar to that caused by bovine herpesvirus was observed in cultured MDBK cells inoculated with the buffy coat fractions from one calf. An attempt was made to identify the virus isolate using a pair of primers to detect glycoprotein B (gB) of BoHV-4 [39]. A 615-bp (base pair) fragment was obtained and sequenced. A BLAST search revealed that the sequence was more than 99% identical to those of the gB regions of most BoHV-4 strains. The isolate was designated as BoHV-4 strain 512. Two other unknown herpesvirus strains, B6010, isolated from nasal swabs of aborted cows from Heilongjiang Province, and J4034, from nasal swabs from six-month-old calves with respiratory diseases from Shandong Province, China, in 2009 that were stored at -70℃, were also propagated in MDBK cells and amplified by PCR and identified as BoHV-4. So far, the isolation and detection of BoHV-4 has not been reported in mainland China. The three Chinese BoHV-4 isolates were further characterized to examine their biological properties, their genomes were partially sequenced, and the resulting sequences were used for phylogenetic analysis. The intention of the present study was to genetically characterize the three wild-type BoHV-4 strains isolated from calves and a dairy cow showing different clinical signs and to provide a foundation for future epidemiological investigations and control of BoHV-4 infections in mainland China.

Materials and methods

Samples, cells and virus isolation

Anticoagulated blood samples were collected from two four-month-old calves experimentally infected with the NCPE strain 3877 of BVDV. The buffy coat fractions were separated from the blood samples using lymphocyte separation medium and resuspended in 1 ml of minimum essential medium (MEM, GIBCO) supplemented with 3% fetal bovine serum (Biochrom AG, Germany). MDBK cells were cultivated in MEM that was supplemented with 10% fetal bovine serum. The separated buffy coat fractions (200 µl) were inoculated onto MDBK cells cultured in 24-well cell culture plates and incubated for five days at 37°C. The cells were harvested by freezing and thawing three times, and three blind passages were carried out until a cytopathic effect (CPE) typical of BoHV-4 was observed. Another two unknown herpesvirus isolates, B6010 and J4034, were also propagated in a monolayer of MDBK cells.

PCR detection of the gB gene of BoHV-4 and nucleotide sequence analysis

A pair of primers, gB1 and gB2, designed based on gB region of BoHV-4 strain Movar 33/63 were used to amplify parts of the gB gene [39]. An approximately 615-bp product was amplified from the BoHV-4 genomic DNA using the primers gB1 (5'-CCC TTC TTT ACC ACC ACC TAC A-3') and gB2 (5'-TGC CAT AGC AGA GAA ACA ATG A-3').

Each virus isolate was inoculated onto a monolayer of MDBK cells, and the cell cultures were harvested by freezing and thawing when CPE was observed. Viral genomic DNA was extracted from 500 µl of the harvested cell culture and detected by PCR as described previously [40]. The conditions for PCR amplification were 45 s at 95 °C, 50 s at 58 °C, and 90 s at 72 °C, with a final extension step at 72 °C for 10 min.

A band of the expected size was obtained by electrophoresis in a 1.5% agarose gel in Tris-acetate EDTA buffer and recovered from the gel using a gel extraction kit (Thermo Fisher Scientific Inc., USA). Direct DNA sequencing was performed using an ABI A3730 automated sequencer and the primers gB1 and gB2. The BLAST search program (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi) was used to identify the sequences.

Electron microscopy

MDBK cells infected with the virus isolates were harvested and negative-stain preparations for transmission electron microscopy were made as described previously [41]. The virions were observed and photographed.

Virus neutralization test

A hyperimmune serum pool against bovine herpesvirus 1 (BoHV-1) was prepared by immunizing rabbits with purified and concentrated BoHV-1. A total of 70 bovine sera were collected from cattle in the Inner Mongolia Autonomous Region, China. A volume of 100 μl of serially diluted serum samples in MEM was mixed with an equal volume of a suspension containing 100 TCID50 of BoHV-4 strain 512. After incubation at 37℃ for 1 hr, the mixture was added to an MDBK cell monolayer, and the residual virus was allowed to adsorb for 1 h. The cells were then kept in a 5% CO2 incubator at 37℃ and examined daily under an optical microscope for one week to observe the appearance of CPE. The neutralizing antibody titer was defined as the highest serum dilution that resulted in reduction of CPE by at least 50%. Negative control serum samples were collected from mock-inoculated rabbits.

Partial sequence analysis of the TK gene of BoHV-4 isolates and phylogenetic analysis

In addition to the gB sequences, partial thymidine kinase (TK) gene sequences of the three Chinese isolates were also determined. A 576-bp product encompassing the TK region was amplified from the three isolates using the primers TK1 (5'-GTT GGG CGT CCT GTA TGG TAG C-3') and TK2 (5'-TGT ATG CCC AAA ACT TAT AAT ATG ACC AG-3') as described previously [9].

A 212-nucleotide region of TK and a 435-nucleotide region of gB for BoHV-4 were used to construct phylogenetic trees for genetic typing of the three Chinese isolates [9]. Reference sequences from previously identified BoHV-4 isolates were also included. The phylogenetic analysis was conducted using Clustal W program and Molecular Evolutionary Genetics Analysis (MEGA) version X (https://megasoftware.net/) with bootstrap values calculated from 1000 replicates. The neighbor-joining algorithm was used for construction of the trees. To confirm the genotyping results, the deduced amino acid sequences of the corresponding proteins of genotypes 1 and 2 were also aligned with the MEGA version X.

Results

The three isolates were cultured and passaged in MDBK cells. All three isolates uniformly produced CPE characteristic of many herpesvirus isolates, with many scattered, rounded, refractory cells and hole formation in the monolayer of MDBK cells. The gB gene fragments of BoHV-4 with the expected size of 615 bp were detected in the three Chinese isolates by PCR (data not shown). The recovered fragments were sequenced and found in a BLAST search to be more than 99% identical to the corresponding sequences of other BoHV-4 strains. Therefore, the three Chinese isolates were identified as BoHV-4. Typical virions of approximately 100-150 nm in diameter were observed in negative-stain preparations of MDBK cells inoculated with the three isolates (Fig. 1). None of the three isolates could be neutralized by a pool of rabbit hyperimmune serum against BoHV-1, showing that there was no cross-neutralizing between BoHV-4 and BoHV-1. Forty-five of the 70 bovine sera (64%) had neutralization antibodies against BoHV-4 in a virus neutralization test performed using strain 512, and indicating that natural BoHV-4 infections had occurred in tested cattle from Inner Mongolia Autonomous Region, China. Two BoHV-4 strains, 512 and B6010, were isolated from cattle in Heilongjiang Province, and one strain, J4034, was isolated from cattle in Shandong Province, China. However, the prevalence of this virus is still unknown, and further research is needed.

Fig. 1

Electron micrographs of the three Chinese BoHV-4 strains propagated in MDBK cell cultures exhibiting typical herpesvirus morphology. A. strain 512. B. strain J4034. C. strain B6010

The partial TK and gB sequences of the three Chinese BoHV-4 isolates were deposited in the GenBank database with the accession numbers MN735172-MN735177. Using the computer program DNAStar (DNAStar Inc., Madison, WI), it was found that the level of conservation among the three Chinese isolates was very high, with 99.0 to 100% nucleotide sequence identity in the gB gene, and 99.8 to 100.0% identity in the TK gene. The partial TK genes amplified using the primers TK1 and TK2 were found by BLAST search to be 95.77 to 100% identical to those of other BoHV-4 strains. Alignment of the amino acid sequences encoded within the amplified gB fragments showed a higher level of conservation than in the nucleotide sequences between the three Chinese isolates. The predicted amino acid sequences of the partial gB genes of the three Chinese isolates were 100.0% identical to each other, and 84.7 to 100.0% identical to the corresponding sequences of other BoHV-4 isolates.

Phylogenetic analysis based on a 212-nucleotide region of the TK gene (Fig. 2A) showed that the three Chinese isolates clustered with BoHV-4 genotype 1 isolates. The phylogenetic analysis based on partial sequences of the gB region yielded similar results. The three Chinese isolates clustered together within the BoHV-4 genotype 1 lineage and formed a separate branch with a bootstrap value of 95% in the gB phylogenetic tree (Fig. 2B). The phylogenetic analysis based on the gB region was also in concordance with that based on the TK region for the BoHV-4 genotype 3 isolates. However, three Argentinean BoHV-4 isolates of genotype 2 (07 568, 08 476, 09 759) did not form a separate branch, showing that the phylogenetic tree constructed using gB sequences was different from that constructed using TK sequences (Fig. 2A and B). Therefore, the TK gene appears to be more suitable than the gB gene for genotyping BoHV-4 isolates. To further evaluate the use of the gB gene for genetic typing of BoHV-4 genotypes 1 and 2, the deduced amino acid sequences of the genotypes 1 and 2 were aligned, and a 9-amino-acid insertion was found at positions 25 to 33 of genotype 1 using MEGA version X (Fig. 3). This was observed previously with Argentinean BoHV-4 isolates [42], revealing significant differences between the gB gene sequences of genotypes 1 and 2. A BLAST search showed that genotypes 1 and 3 also had an 8-amino-acid insertion at positions 44 to 51 (data not shown). Therefore, the differences in gB amino acid sequences support the use of the gB gene for genetic typing of BoHV-4 genotypes 1, 2 and 3. The data obtained in this study also support use of the 212-bp fragments of the TK gene for genetic typing of BoHV-4 genotypes 1, 2 and 3.

Fig. 2

Phylogeny of the three Chinese BoHV-4 isolates. The TK phylogenetic tree was constructed using TK gene nucleotide sequences of the three Chinese isolates obtained in this work and sixteen BoHV-4 reference isolates retrieved from GenBank (A TK). The gB phylogenetic tree was constructed using the gB gene nucleotide sequences of the three Chinese isolates detected in this work and eight BoHV-4 reference isolates retrieved from GenBank (B gB). The trees were generated using the Neighbor-Joining method and bootstrap testing. Numbers at branches indicate the percentage of 1000 bootstrap replicates that support each phylogenetic branch. The BoHV-4 reference isolates and their GenBank accession numbers are as follows: B11-41 (AB035515), Movar 33/63 (AB035516), 86 068 (AB035517), V.test (JN133502), 66-p-37 (AF318573, Z15044), 07 759 (JQ838048), 08 209 (JQ838049), 08 263 (JQ838050), 8 362 (JQ838051), 08 433 (JQ838054), granulose (KP209014), DN-599 (JQ838062), 07 435 (JQ838046, KP209016), 09 227 (JQ838056, KP209024), MGArom (EU244697), MGA1075 (EU244700). The isolates 09 508 (KU180394), 07 568 (KP209017), 08 476 (KP209023), and 09 759 (KP209027) were only used in the gB tree

Fig. 3

Alignment of partial gB amino acid sequences of genotypes 1 (66-p-347, 512, B6010, J4034, V.test and 09_508) and 2 (07_568, 08_476 and 09_759) of BoHV-4 (MEGA Version X). Missing bases in the sequence are indicated by dashes

Discussion

Recently, analysis of BoHV-4 strains isolated from aborted cows in Argentina revealed a high degree of genetic divergence among strains, which could be classified into three different groups: genotypes 1, 2, and 3 [34, 42, 43]. Our preliminary results revealed that the three Chinese isolates clustered within BoHV-4 genotype 1 in both the TK and gB phylogenetic trees. The phylogenetic analysis using a portion of the gB region was in concordance with that using TK gene sequences for genotypes 1 and 3. This indicated that the gB region might be a candidate target for typing BoHV-4 genotypes 1 and 3, but three Argentinean BoHV-4 isolates of genotype 2 did not form a separate branch (Fig. 2B). Therefore, the phylogenetic analysis results presented here support the use of the TK gene as a suitable target for genotyping BoHV-4 but do not fully support the use of the gB gene as a candidate target for genotype 2. A 9-amino-acid insertion was found at positions 25 to 33 of gB of BoHV-4 genotypes 1 and 2 (Fig. 3), and an 8-amino-acid insertion was found at positions 44 to 51 of gB of BoHV-4 genotypes 1 and 3. In general, the data obtained from the gB phylogenetic tree also supported the use of a 212-bp fragment of the TK gene as a valid target for genetic typing of BoHV-4.

In order to evaluate whether the level of genetic variation presented here is typical of BoHV-4 genotypes 1, 2 and 3, more BoHV-4 genome sequences might be needed. It is notable that BoHV-4 genotypes 2 and 3 have only been detected in America, and the novel genotype 3 has only been detected in South America. Prior to this study, one BoHV-4 strain was isolated from bovine arterial endothelial cells in Taiwan, China [22]. Another strain of BoHV-4, B11-40 was isolated from the spinal cord of a cow with astasia in Japan [38]. Three BoHV-4 strains, 512, B6010, and J4034 were isolated from cattle in mainland China in this study. To date, five Asian strains have been isolated, and four of them clustered in BoHV-4 genotype 1 with the Movar-like strains (European prototype). Therefore, the results presented here suggest that the differences between genotypes 2 and 3 might be a result of geographic isolation. Sequence analysis of more BoHV-4 field isolates from different geographic regions is needed, particularly for genotypes 2 and 3, to investigate differences between genotypes. Further immunological investigations might be required to determine the biological relevance of the three proposed BoHV-4 genotypes and their importance for future vaccine development. The identification of sequence variability is also valuable for selecting appropriate primers and excluding false-negative PCR results for BoHV-4.

BoHV-4 has been identified as one of the most important of the known viral pathogens of both young and adult cattle with respiratory or reproductive problems [2, 6, 12, 44–47]. Recently, a strong association was found between BoHV-4 infection and lactation number, lactation stage, and postpartum metritis [6, 48]. Respiratory and reproductive problems continue to cause substantial economic losses for the global cattle industry. BoHV-1, BVDV, bovine respiratory syncytial virus (BRSV), bovine parainfluenza virus type 3 (BPIV3), and bovine adenovirus type 3 (BAV-3) are the main recognized viral pathogens involved in the bovine respiratory disease complex (BRDC) [49, 50], but the only recognized virus involved in BRDC in China is BoHV-1, which was first isolated from cattle imported from New Zealand in 1980 [51]. BVDV infection is very common in China; however, its role in BRDC is not clear [52]. Other important viruses, such as BRSV, BPIV3 and BVA-3, have been isolated or detected in China [40, 41, 53]. So far, the only viruses known to be involved in bovine reproductive problems in China are BoHV-1 and BVDV. The isolation and detection of BoHV-4 was not reported prior to this study. This is the first report about the isolation of BoHV-4 in mainland China. Two Chinese BoHV-4 strains, B6010 and J4034, were isolated from cattle showing different clinical signs. The B6010 strain was the only virus isolated from an aborted cow. This is notable because BoHV-1 and BVDV are the main known viral pathogens involved in abortion in pregnant cows in China. This indicated that the BoHV-4 might be involved in abortion in pregnant cows. Another strain, J4034, was isolated from calves showing respiratory diseases, and it might be involved in calf respiratory disease. The BoHV-4 strain 512 was isolated from a four-month-old calf experimentally infected with an NCPE BVDV strain, but no BVDV was isolated from the calf after inoculation with BVDV. A previous infection with BoHV-4 might be a good explanation for the failure of BVDV isolation from this calf and suggests that BoHV-4 might interfere with experimental infection with BVDV. However, the role of BoHV-4 infection in respiratory or reproductive problems in cattle in China is still not known and needs further investigation.

In conclusion, partial gB gene sequences of the three Chinese BoHV-4 isolates from this study were 99.0% to 100.0% identical to each other, and all three appeared to belong to genotype 1. A virus neutralization test using hyperimmune sera against BoHV-1 showed no cross-neutralization between BoHV-4 and BoHV-1. About 64% of 70 bovine sera (45/70) collected from Inner Mongolia Autonomous Region had neutralizing antibodies against BoHV-4 strain 512. This suggests that natural BoHV-4 infection in cattle might be more widespread than expected in China and that an extensive monitoring program is needed. A serological investigation will be useful to assess the host immune response to the virus. On the basis of this research, detection methods could be established for measuring antibodies against BoHV-4 and implementing epidemiological investigations of BoHV-4 in China. The new isolates could be evaluated for future BoHV-4 vaccine development. Extended research on the impact of BoHV-4 on BRDC and reproductive problems in China will be beneficial to the cattle industry.