key: cord-0761486-tp6mlqtw
authors: Li, Yingli; Chang, Jitao; Wang, Qian; Yu, Li
title: Isolation of two Chinese bovine enteroviruses and sequence analysis of their complete genomes
date: 2012-08-01
journal: Arch Virol
DOI: 10.1007/s00705-012-1424-6
sha: 9a09287031c2bfe6b8683f61181e7bf94fe3b2a2
doc_id: 761486
cord_uid: tp6mlqtw

In this study, RNA corresponding to bovine enterovirus (BEV) was detected in 24.6 % of faecal samples (17/69) from diarrheic and healthy cattle in six different areas in China by an RT-PCR screening method. Furthermore, two cytopathic agents, designated as BHM26 and BJ50, were isolated from the bovine diarrheic fecal samples. During passage in MA104 cells, ultrathin sections of virus-infected monolayers were examined using a transmission electron microscope, and a large number of symmetrical virus crystals were seen in the cytoplasm, with monomorphic small viral particles of 27-30 nm in diameter. The full-length RNA genomes were 7433 and 7416 nucleotides long, respectively, with a genome organization analogous to that of picornaviruses. Phylogenetic analysis of the VP1 and VP3 capsid protein coding sequences suggested that the viruses BHM26 and BJ50 belong to genotype 2 of the BEV cluster B (BEV-B). In addition, sequence comparisons of the 5′ and 3′ UTRs and P1, P2 and P3 subgenomic regions of the two isolates suggested that there were intergenotypic recombination events occurring during evolution of the BHM26 and BJ50 isolates.

Bovine enteroviruses (BEVs), which are members of the family Picornaviridae, genus Enterovirus, form two phylogenetic clusters (designated BEV-A and -B), with two and three genotypes, respectively [15] . BEV virions consist of a naked capsid that surrounds a core of positive-sense, single-stranded RNA. Hydrated native particles are 27-30 nm in diameter in electron micrographs, which reveal the virion to be a smooth sphere. The BEV genome is approximately 7.4 kb long and has a typical picornavirus genome organization, including VPg following the 5 0 UTR, structural (VP4, VP2, VP3 and VP1) regions, nonstructural (2A, 2B, 2C, 3A, 3B, 3C and 3D) regions, 3 0 UTR and a poly (A) tail [2] .

BEVs were first isolated in the late 1950s [3, 6] . They were isolated from (i) faeces of cattle with symptoms of pneumonia, respiratory disease, enteritis, dysentery and fertility disorders, (ii) fetal fluids from an aborted calf, (iii) the faeces of apparently healthy animals, and (iv) sewage, treated effluents, and materials from a farm environment. Difficulties in reproducing clinical symptoms following experimental infection of animals led to the conclusion that BEVs were of only minor veterinary medical importance. However, they have been shown to be widely detectable in a farm environment and have thus been proposed as markers for fecal contamination from cattle [7] . It is generally assumed that BEVs are avirulent, and if so, it is possible that they would make useful vectors for delivery of small molecules to the gut. Additionally, BEVs have been shown to have oncolytic properties and may be useful as therapeutic agents [12, 13] . For these various reasons, we were interested in learning about the epidemiology and molecular evolution of BEVs from China.

In this study, a preliminary epidemiological study was performed, the results of which suggested that BEVs are present in certain cattle in China. Also, the present report described the first isolation, by means of MA104 cell cultures, of two isolates of BEV (designated BHM26 and BJ50) obtained from rectal swabs from diarrheic cattle of dairy herds in China. These viruses were characterized molecularly through sequence analysis of complete genomes, and the BHM26 and BJ50 isolates were therefore classified into genotype 2 of the BEV-B cluster. The complete genome sequences of two Chinese BEV isolates not only enhance the database of sequences of members of the family Picornaviridae, providing data of use to picornavirus virologists, but they also can be used further for construction of a BEV infectious clone and development of BEV-vectored vaccines based on epitopes.

A total of 69 fecal samples were collected from cattle at six farms in China from October 2007 to July 2010. Out of a total of 69 fecal samples, 59 were from animals with diarrhea and 11 were from non-diarrheic animals. The samples were collected from Inner Mongolia (n = 11), Heilongjiang (n = 12), Anhui (non-diarrheic samples, n = 11), Liaoning (n = 12), Xinjiang (n = 12) and Shanghai (n = 11), which are the main production areas of cattle in China. Stool specimens were stored at -70°C until the testing was performed.

The samples were diluted 1:10 in 10 mM phosphate-buffered saline (PBS) (pH 7.4), clarified by low-speed centrifugation (3000 9 g for 10 min), and kept at -70°C until they were used to isolate viruses. Monolayers of MA104 cells grown for 1 or 2 days were used for virus isolation. The growth medium was Dulbecco 0 s minimum essential medium (DMEM) supplemented with 8 % fetal bovine serum (FBS) and antibiotics (100 U of penicillin per ml and 100 lg of streptomycin per ml). The centrifuged fecal samples were filtered through 0.22-lm filters (Millex, Millipore). Subsequently, 1 ml of filtrate was inoculated onto a MA104 cell monolayer that had been grown in flasks and washed three times with PBS. After adsorption for 60 min at 37°C, 4 ml of fresh medium (without FBS) was added, the cultures were incubated at 37°C in a CO 2 incubator. When cytopathic effect (CPE) was not observed one week after inoculation, or the cells showed 50 to 80 % CPE, the cultures were frozen and thawed three times for harvest of cell lysates. Subsequent passages were carried out in the same manner by using 0.1 ml of cell suspension from the previous passage as the inoculum.

Total RNA from fecal samples and virus-infected MA104 cells was extracted by using TRIzol LS Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer 0 s instructions. Extracted viral RNA was used as a template for cDNA synthesis using PrimeScript Reverse Transcriptase (TaKaRa). Oligonucleotide primers (Table 1 ) used in this study were designed according to the consensus sequences of the 12 fully sequenced BEV RNA genomes available from the GenBank database. RT-PCR was carried out with an initial reverse transcription step of 60 min at 42°C, followed by PCR for activation at 94°C for 5 min, 30 cycles of denaturation at 94°C for 0.5 min, annealing for 0.5 min at 51°C to 55°C depending on the T m values of the different primers, and extension at 72°C for 1.5 min using a thermal cycler (Bio-Rad, USA). The final extension step was done at 72°C for 10 min. The 5 0 ends of the virus genomes were amplified using a 5 0 -Full RACE Kit (TaKaRa, Dalian, China) following the manufacturer 0 s instructions. Sequences at the 3 0 end of the virus genomes were determined using RT-PCR, with first-strand cDNA synthesis carried out using an oligo (dT) primer, thus taking advantage of the poly (A) tail in BEV RNA genomes as a priming site. All of the PCR products were detected by electrophoresis in 1 % agarose gels and were purified using a DNA Gel Extraction Kit (Watson Biotechnologies, Inc., Shanghai, China). They were then sequenced directly in both directions by a commercial service (Invitrogen, China).

The complete nucleotide sequences of the isolates were assembled from the sequence fragments using the SeqMan module of the DNASTAR software package (DNASTAR Inc., Madison, WI). Amino acid sequences were deduced from the nucleotide sequences using the EditSeq module. Multiple sequence alignments and phylogenetic trees were constructed by the neighbor-joining method using Molecular Evolutionary Genetics Analysis (MEGA, version 4.0) software, and bootstrap resampling analysis of 1000 replicates was performed.

A total of 69 bovine diarrheic and healthy faecal specimens from six different areas were examined by RT-PCR using the primer pair 6424U24 and 7450L24, which cover the highly conserved partial 3D gene and 3 0 UTR of BEVs. The results showed that virus-specific RNA could be detected directly in faecal material in 17 of 69 samples (24.6 %). The prevalence rate of BEV in the six cattle farms located in Inner Mongolia, Heilongjiang, Anhui, Liaoning, Xinjiang, and Shanghai was 45.4 % (5/11), 25 % (3/12), 27.2 % (3/11), 16.7 % (2/12), 25 % (3/12) and 9 % (1/11), respectively, indicating that BEVs are present in certain domestic cattle of China. All of the BEV-specific-RNA-positive faecal samples from Inner Mongolia and Xinjiang were inoculated onto MA104 cells for isolation of BEV. After the centrifuged fecal samples were inoculated onto cell monolayers, CPE was observed in the second-passage culture on the third day in culture plates inoculated with stool specimens 26 (from Inner Mongolia) and 99050 (from Xinjiang). The cultivated virus derived from sample 26 was designated BHM26, and that from the sample 99050 was designated BJ50. Distinctly recognizable CPE was observed until the third passage at 12 to 24 h after infection. The CPE produced by the two isolates consisted of obscure cell borders, cell rounding, piling up of cells, and detachment of cells from the surface of plates (Fig. 1) . Ultrathin sections of the infected monolayers were stained with uranyl acetate and lead citrate and examined using a transmission electron microscope. A large number of symmetrical virus crystals were seen in the cytoplasm, and the monomorphic small viral particles with a diameter of 27-30 nm were consistent with viruses of the family Picornaviridae (Fig. 2) . The presence of BEVs in cell cultures was further detected by RT-PCR with the primer set 6424U24 and 7450L24, indicating that the two isolates were BEV positive. PCR products were sequenced directly in both directions, and sequence alignments were conducted by using MEGA version 4.0. The results showed that the 3D/3 0 UTR (1026 nt) of BHM26 and BJ50 shared 86 % and 88 % nucleotide sequence identity, respectively, to that of BEV strain Jena38/02 (GenBank accession number: DQ092789), indicating that the two isolates are bovine enterovirus.

Nucleotide sequences of the complete BHM26 and BJ50 genomes Several overlapping PCR amplifications spanning the entire genomes of the BHM26 and BJ50 isolates were obtained, and their complete nucleotide sequences were determined. The BHM26 and BJ50 genomes consist of 7,433 bp and 7,416 bp, including 42 and 21 A bases, respectively, in the poly (A) tail. The complete genomes of the two isolates have a typical picornavirus genome organization, including a 5 0 UTR, a large single ORF, a 3 0 UTR 

Low nucleotide sequence identity of the genes encoding the capsid proteins 1C and 1D is the main criterion for the identification of genotypes of BEV [15] . Thus, phylogenetic analysis of the 1C and 1D encoding regions of the two isolates and other representative BEV strains was performed to determine the phylogenetic characteristics of the two BEV isolates. It was found that both BHM26 and BJ50 were classified into genotype 2 of the BEV-B cluster and were most closely related to BEV strain Jena38/02 (Fig. 3) , an isolate from Germany (GenBank accession number: DQ092788). Furthermore, similar results were obtained when the 1C and 1D encoding regions were used to analyze phylogenetic relatedness among the BEV strains, although the branch orders were somewhat different (Fig. 3) , showing that the 1C and 1D coding regions can be used simultaneously to determine the genotypes of the BEV strains.

The complete genomes of BHM26 and BJ50 were divided into 5 0 UTR and 3 0 UTR and P1, P2 and P3 subgenomic regions for sequence comparisons with those of the referenced BEV strains. For these subgenomic regions, the two sequenced isolates shared 88.6 % (5 0 UTR), 94.1 % (3 0 UTR), 80.6 % (P1), 80.9 % (P2) and 82.2 % (P3) nucleotide sequence identity ( 

In the present study, we first conducted a preliminary epidemiological survey for detecting BEVs from diarrheic and healthy cattle in of China. A total of 69 bovine fecal samples collected from the six different areas of China were tested by RT-PCR, and 17 samples (24.6 %) in the six areas were positive for BEV, indicating that BEVs are widespread in dairy farms in China. Subsequently, two BEV isolates (designated BHM26 and BJ50) were isolated from bovine fecal specimens, and they were then characterized morphologically and genotypically. In addition, the complete sequences of the BHM26 and BJ50 genomes were determined and were submitted to the GenBank database under accession nos. HQ917060 and HQ917061, respectively. By phylogenetic analysis of the 1C and 1D capsid protein genes, both BHM26 and BJ50 were classified into genotype 2 of the BEV-B cluster. Sequence comparisons of the 5 0 UTR and 3 0 UTR and P1, P2 and P3 subgenomic regions suggested that BHM26 and BJ50 are intergenotypic recombinant BEV isolates. In general, incongruence between phylogenies of different genome regions is considered indicative of recombination events in enteroviruses [8, 11] . For example, phylogenetic analysis of German BEV isolate D 14/3/96 yielded incongruent results [15] in which comparisons of 1B and 1C indicated a similarity with genotype 2 of the BEV-A cluster; however, in the 1D comparison, the D 14/3/96 strain clustered with genotype 1. Another example of intragenotypic recombination is shown by the LC-R4 strain, which was grouped with the Vir 404/03 strain in the 1B phylogeny but with the VG(5)27 strain in the 1C and 1D phylogenies [15] . These findings are an indication of intergenotypic or intragenotypic recombination in the evolution of the BEVs. In this study, for the China BEV isolates BHM26 and BJ50, sequence analysis based on the 5 0 UTR and 3 0 UTR and P1, P2 and P3 subgenomic regions yielded some incongruent results, suggesting that there were intergenotypic recombination events occurring during evolution of the BHM26 and BJ50 isolates.

The earliest reports concerning the virulence of BEVs suggested that there might be an association between BEV infections and a range of diseases in cattle, as viruses were isolated from cattle with clinical signs that varied from respiratory to enteric to reproductive disease and infertility [1, 14] . However, bovine enterovirus has also been isolated from apparently healthy cattle and wildlife [5, 6, 9] . Interpretation and comparison of these reports are difficult because in many cases other infectious agents were identified or were not investigated. As reviewed by Jiménez-Clavero and co-workers [4] , enteroviruses are highly stable in the digestive tract, are shed in large amounts, and persist in the environment for prolonged times, and thus serve as indicators of fecal pollution. In susceptible individuals, enteroviruses cause subclinical infection or mild disease, although occasional infections may cause serious disease, such as the poliomyelitis caused by poliovirus in humans [10] . The BHM26 and BJ50 were isolated from bovine diarrheic fecal samples, but their pathogenicity in cattle experimentally infected with the BHM26 and BJ50 has not been conducted. It is well known that bovine rotavirus and coronavirus are the most important cause of diarrheal diseases in cattle, especially in neonatal calves. If one animal is simultaneously infected by rotavirus and coronavirus, it can cause significant morbidity and mortality. In this study, BEV-specific RNA not only was detected in diarrheic faecal samples (14 of 17 BEV-positive samples collected from diarrheic cattle) but also was detected in healthy faecal samples (3 of 17 BEV-positive samples collected from healthy cattle). At the same time, two other RNA viruses that are transmitted by the fecal-oral route (bovine coronavirus and rotavirus) were detected from diarrheic faecal samples. Although little is actually known regarding the pathogenesis of the Chinese BEV isolates within the enteric tract, according to the epidemiological survey in this study, the BEVs in the fecal specimens from China are only considered to be a minor cause of diarrhea in these dairy farms.

Several characteristics of the BEVs make them good candidates for development of a vaccine vector. First, they are physically very stable, e.g., BEV particles are resistant to lipid solvents (e.g., chloroform) and are stable at pH 2 to 10. Second, BEV is currently considered avirulent in cattle, making it safe to use as a vaccine vector. Third, since enteroviruses enter their host via the enteric tract, they could be administered orally as vaccine vectors to stimulate mucosal immunity. Thus, we are interested in their potential as vaccine vectors for expressing small peptides. We have recently constructed an infectious cDNA clone of the BEV BHM26 isolate on the basis of these sequencing data (Yingli Li et al., unpublished data). Our previous studies have finely mapped a conserved type O FMDV neutralizing epitope recognized by a specific monoclonal antibody [16] . Engineering of recombinant BEVs displaying the FMDV neutralizing epitope in our laboratory is underway for development of a BEV-vectored and epitopebased vaccine. 

Bovine viral diarrhea virus type 2-induced meningoencephalitis in a heifer

The Complete nucleotide sequence of a bovine enterovirus

Isolation and classification of viruses from cattle during outbreaks of mucosal or respiratory disease and from herds with reproductive disorders

Survey of bovine enterovirus in biological and environmental samples by a highly sensitive Real-Time Reverse Transcription-PCR

Fields virology, 4th edn

A serological classification of bovine enteroviruses

Bovine Enteroviruses as indicators of fecal contamination

Evolution of the genome of human enterovirus B: incongruence between phylogenies of the VP1 and 3CD regions indicates frequent recombination within the species

Molecular characterisation of Australian bovine enteroviruses

One hundred years of poliovirus pathogenesis

Evidence of recombination among enteroviruses

Therapeutic effects of bovine enterovirus infection on rabbits with experimentally induced adult T cell leukaemia

Bovine enterovirus as an oncolytic virus: Foetal calf serum facilitates its infection of human cells

Bovine enterovirus-1: characterization, replication and cytopathogenic effects

Molecular-based reclassification of the bovine enteroviruses

Identification of a conserved linear epitope on the VP1 protein of serotype O foot-and-mouth disease virus by neutralising monoclonal antibody 8E8