Genetic diversity of epidemic enterovirus 71 strains recovered from clinical and environmental samples in Taiwan

Most enteroviruses excreted in human feces and urine are present in environmental water. In order to clarify the infection route of enterovirus, the detection of viruses in both clinical and environmental samples may contribute to understanding the mode of transmission of strains responsible for human infection. Thus, 21 epidemic enterovirus 71 strains from environmental water or stool samples were collected from HFMD children during 2005. Enterovirus genomic RNA was first amplified directly from clinical and environmental samples and then characterized by DNA sequencing and phylogenetic analysis. Results showed that these clinical strains share similar sequence identity (86.4-86.8%) to prototype BrCr based on the 5'-nontranslated region (NTR). However, environmental strains, except HME 77, share similar sequence identity (86.2-87.2%) to prototype BrCr. HME 77 showed higher sequence identity (90.1%). Results from phylogenetic analysis revealed that five environmental isolates were clustered as genogroup 3, which was closely related to a Taiwan outbreak in 1998 (AY055133). HME 77 was more closely related to a China epidemic isolate (AY895144), which belonged to genogroup 4. In contrast, all strains from clinical samples tested belonged to genogroup 3, which clustered with AY055133. In conclusion, there are two major epidemic clones (genogroups 3 and 4) prevalent in Taiwan since 2004 either in water or clinical patients.     

Genomic characterization of an enterovirus 97 strain isolated in Shandong,China

The genomic characterization of human enterovirus 97 (EV97) strain isolated from an acute flaccid paralysis case in Shandong province, China in 1999, is described. The strain, designated as 99188/SD/CHN/1999/EV97 (abbreviated as 99188), had a genome of 7394 nucleotides. Compared with other EV97 strains, it had 81.3–83.3% nucleotide similarity and 94.0–95.4% amino acid similarity in VP1 coding region, and it had 81.4% complete genomic similarity with prototype strain BAN99-10355. The most striking feature was the deletion of 18 nucleotides in the 3′ end of VP1 coding region, combined with two deletions and one insertion in 5′ and 3′ untranslated regions. All these findings demonstrated the strain 99188 had a distant genetic relationship with other EV97 strains. In the phylogenetic trees generated from VP1 and 3D sequences of human enterovirus species B (HEV-B), the lineages of strain 99188 were not congruent, suggesting the event of recombination. Similarity plot analysis further provided the evidence of recombination with other strains of HEV-B in P2 and P3 coding region. This is the first finding of EV97 in China and the third genomic sequence of EV97 reported.     

Complete genome analysis of human enterovirus B73 isolated from an acute flaccid paralysis patient in Shandong,China

Human enterovirus B73 (EV-B73) is a member of species Enterovirus B. To date, only one complete genome sequence of prototype strain CA55-1988 from California has been available. In this study, the complete genome analysis of an EV-B73 strain 088/SD/CHN/04 isolated from an acute flaccid paralysis case in Shandong Province, China in 2004 is conducted. It had 75.6 and 79.3 % nucleotide similarity with prototype strain CA55-1988 in the VP1 coding region and the complete genome, respectively. It had great VP1 nucleotide divergence (16.7–24.4 %) with EV-B73 strains from other parts of the world. Similarity plot and bootscanning analyses provided evidence of recombination with other EV-B viruses.     

Molecular analysis of an echovirus 3 strain isolated from an individual concurrently with appearance of islet cell and IA-2 autoantibodies

Growing evidence has implicated members of the genus Enterovirus of the family Picornaviridae in the etiology of some cases of type 1 diabetes (T1D). To contribute to an understanding of the molecular determinants underlying this association, we determined the complete nucleotide sequence of a strain of echovirus 3 (E3), Human enterovirus B (HEV-B) species, isolated from an individual who soon after virus isolation developed autoantibodies characteristic of T1D. The individual has remained positive for over 6 years for tyrosine phosphatase-related IA-2 protein autoantibodies and islet cell autoantibodies, indicating an ongoing autoimmune process, although he has not yet developed clinical T1D. The sequence obtained adds weight to the observation that recent enterovirus isolates differ significantly from prototype strains and provides further evidence of a role for recombination in enterovirus evolution. In common with most HEV-B species members, the isolate exhibits 2C and VP1 sequences suggested as triggers of autoimmunity through molecular mimicry. However, comparisons with the E3 prototype strain and previously reported diabetogenic and nondiabetogenic HEV-B strains do not reveal clear candidates for sequence features of PicoBank/DM1/E3 that could be associated with autoantibody appearance. This is the first time a virus strain isolated at the time of commencement of beta-cell damage has been analyzed and is an invaluable addition to enterovirus strains isolated previously at the onset of T1D in the search for specific molecular features which could be associated with diabetes induction.     

Cloning and Nucleotide Analysis of Segment A Gene of Infectious Bursal Disease Virus Detected in Korea

A strain of infectious bursal disease virus (IBDV) was detected from bursal tissues of chicks which suffered from infectious bursal disease (IBD) in Chinju, Korea and provisionally named as Chinju strain. A full-length cDNA clone for segment A gene of the virus was constructed, and complete nucleotide sequence of the gene including noncoding region was determined and analyzed by comparison with that of other IBDV strains. The segment A gene of Chinju strain consisted of 3,269 nucleotides including 862 adenine (26.4%), 917 cytosine (28.0%), 854 guanine (26.1%) and 636 thymine (19.5%). There were regions for two open reading frames (ORFs), ORF1 encoding the VP5 with ATG codon at nucleotides 98–100 and ORF2 encoding the polyprotein of VP2, VP4 and VP3 in the nucleotides 132–3,170. In deduced translation the ORF2 encoded 1,012 amino acids. The full nucleotide sequence of segment A gene and amino acid sequence of ORF2 of the Chinju strain showed 98–99% homology with those of the very virulent IBDVs (vvIBDVs) such as HK46, OKYM, UK661, UPM97/61, D6948 and BD3/99. In phylogenetic analysis of nucleotide and amino acid sequences, the Chinju strain was also related closely to the vvIBDVs. Hence, it was suggested that the Chinju strain is a vvIBDV. The nucleotide and amino acid sequences of the Chinju strain with pertinent information can be useful for the development of genetically engineered vaccines and diagnostic reagents against vvIBDV.     

Molecular evidence of persistent echovirus 13 meningoencephalitis in a patient with relapsed lymphoma after an outbreak of meningitis in 2000

Enteroviral meningoencephalitis was diagnosed in a patient with an immunodeficiency syndrome acquired after treatment with rituximab for a relapsed primary B-cell lymphoma. A second meningoencephalitic episode was diagnosed 6 months later and was successfully treated with a combination of immunoglobulins and pleconaril. The infection was persistent since the enterovirus genome was detected in five sequential specimens of cerebrospinal fluid collected over 9 months. An echovirus 13 isolate was isolated in the first three samples. The viral sequence encoding the VP1 capsid protein of the three isolates was determined and was compared with that of four control viruses. The virus isolates recovered from the patient shared >99% nucleotide sequence similarity with one another. In a phylogenetic tree, they were directly related to a control virus obtained from a patient hospitalized in 2000 during an outbreak of enterovirus meningitis. The epidemiological origin of a chronic echovirus infection in a patient with immune deficiency suggests that the echovirus had been continuously circulating in the general population after the outbreak that had revealed its emergence.     

Complete nucleotide sequence of Alfalfa mosaic virus isolated from alfalfa (Medicago sativa L.) in Argentina

The complete nucleotide sequence of an Alfalfa mosaic virus (AMV) isolate infecting alfalfa (Medicago sativa L.) in Argentina, AMV-Arg, was determined. The virus genome has the typical organization described for AMV, and comprises 3,643, 2,593, and 2,038 nucleotides for RNA1, 2 and 3, respectively. The whole genome sequence and each encoding region were compared with those of other four isolates that have been completely sequenced from China, Italy, Spain and USA. The nucleotide identity percentages ranged from 95.9 to 99.1 % for the three RNAs and from 93.7 to 99 % for the protein 1 (P1), protein 2 (P2), movement protein and coat protein (CP) encoding regions, whereas the amino acid identity percentages of these proteins ranged from 93.4 to 99.5 %, the lowest value corresponding to P2. CP sequences of AMV-Arg were compared with those of other 25 available isolates, and the phylogenetic analysis based on the CP gene was carried out. The highest percentage of nucleotide sequence identity of the CP gene was 98.3 % with a Chinese isolate and 98.6 % at the amino acid level with four isolates, two from Italy, one from Brazil and the remaining one from China. The phylogenetic analysis showed that AMV-Arg is closely related to subgroup I of AMV isolates. To our knowledge, this is the first report of a complete nucleotide sequence of AMV from South America and the first worldwide report of complete nucleotide sequence of AMV isolated from alfalfa as natural host.     

Analysis of the full-length genome of a subgenotype IIIB hepatitis A virus isolate: primers for broadly reactive PCR and genotypic analysis

Among six known subgenotypes (IA, IB, IIA, IIB, IIIA, and IIIB) of human hepatitis A virus (HAV), the complete genomic sequence has not been determined for IIIB. In this study, the full-length genomic sequence of a IIIB HAV isolate (HA-JNG06-90F) recovered from a Japanese patient who contracted sporadic hepatitis A in 1990, was determined. The HA-JNG06-90F genome, which comprised 7462 nt excluding the poly(A) tail, was related most closely to NOR-21 of subgenotype IIIA with an identity of 89.1%, and was only 82.6-83.4% similar to human HAV isolates of genotypes I and II over the entire genome. Comparison of full-length genomic sequences of 20 reported isolates and HA-JNG06-90F generated optimal results for separation of different levels: the nucleotide identities were 80.7-86.6% at the genotype level, 89.1-91.9% at the subgenotype level, and 94.6-99.7% at the isolate level. Similar ranges of nucleotide identity were observed when comparing partial nucleotide sequences of the VP1-2B (481 nt; primer sequences at both ends excluded) and 3C/3D (590 nt) regions, which were amplifiable by PCR with primers designed from well-conserved areas of the HAV genome. All 66 samples with IgM-class HAV antibodies tested positive for HAV RNA by both VP1-2B (481 nt)-PCR and 3C/3D (590 nt)-PCR: subgenotype assignment was concordant in all samples tested (IA [n = 61], IB [n = 1], IIIA [n = 2] and IIIB [n = 2]). These results suggest that two broadly reactive PCRs using primers derived from the VP1-2B and 3C/3D regions, respectively, may be applicable to universal detection and phylogenetic analysis of various HAV strains.     

Influence of the genetic heterogeneity of the ISDR and PePHD regions of hepatitis C virus on the response to interferon therapy in chronic hepatitis C

Molecular typing of enteroviruses should ideally focus on regions encoding determinants for neutralization. Mapping of monoclonal neutralizing antibodies has shown the VPI protein, in particular its aminoterminal part, encompassing the B-C loop, to be one major antigenic region. We therefore sequenced 570 nucleotides from the 5'-end of the VP1 region of the genome for all 28 echovirus prototypes, and for 61 clinical isolates representing all different echovirus types. An analysis of 133 sequences, including 39 sequences retrieved from GenBank, classified all echoviruses in enterovirus group B confirming results from sequencing within the VP2 region. The nucleotide and amino acid divergence of VP1 sequences of homotypic strains varied from 7.5-23.0% and from 0.0-5.3%, respectively, when compared to their corresponding prototypes, whereas strains belonging to different serotypes these divergences were 22.1-38.9 % and 4.9-16.4 %, respectively. Despite these minimal overlaps, the VP1 sequence was always more similar to that of the homotypic prototype than to that of any heterotypic strain. For 13 out of 14 echovirus types, where multiple isolates were available, the corresponding VP1 sequences diverged more from those of the prototype than from the other homotypic sequences as a reflection of genetic drift. Because there was a complete concordance between the sequences of the region encoding the VP1 aminoterminus and the serotype (P< 0.00001) sequence analysis of this region might complement typing by neutralization, and classify correctly echovirus isolates that may not be typed conveniently by the antisera in hand.     

Complete nucleotide sequence analysis of the VP1 genomic region of Echoviruses 6 isolated from sewage in Greece revealed 98% similarity with Echoviruses 6 that were characterized from an aseptic meningitis outbreak 1 year later

The molecular characterization of two enterovirus strains (LR51A5 and LR61G3) isolated from the sewage treatment plant unit in Larissa, Greece, in May and June 2006 and the investigation of their relationship with enteroviruses of the same serotype isolated in Greece in 2001 and 2007 were performed by complete VP1 sequence analysis of the isolates. The close phylogenetic relationship and the high nucleotide similarity (98%) led to the conclusion that the virus isolated from sewage in 2006 was associated with that isolated from an aseptic meningitis outbreak 1 year later. Bootscan analysis of the VP1 genomic region revealed that intraserotypic multi-recombination events might have been involved in the evolutionary past history of the LR51A5 and LR61G3 isolates.     

Isolation and characterization of a Chinese strain of human enterovirus 74 from a healthy child in the Tibet Autonomous Region of China

Human enterovirus 74 (HEV74) is a recently described serotype within the species Human enterovirus B (HEV-B). Few nucleotide sequences of HEV74 are available, and only one complete genome sequence (the prototype strain) has been published. In this study, we report the complete genome sequence of an HEV74 strain isolated from a healthy child during a stool survey in the Tibet Autonomous Region of China. The results indicated that HEV74 may be a prevalent and common enterovirus type, and that HEV74 is globally distributed, especially in Asia. Sequence analysis revealed high variability among HEV74 strains and indicated frequent recombination within HEV-B.     

Nucleotide sequence of the gene encoding the RNA polymerase and the 3′ non-coding region of a bovine enterovirus Japanese isolate: rapid synonymous substitutions between European and Japanese strains

Summary.  The nucleotide sequences of the genome RNA encoding the RNA polymerase and the 3′ non-coding region (NCR) of bovine enterovirus (BEV) serotype I Japanese isolate, MZ468, were determined. The genetic distance between the two BEV serotype I strains, MZ468 and VG-5-27, was calculated by pairwise comparison of nucleotide sequences. The synonymous substitution rate was high (1.40 × 10⁻²/site/year), and of the same order as those of influenza virus HA, HIV-1 gag and env, and enterovirus 70 VP1 genes. Accepted December 1, 1997 Received June 2, 1997     

A new picornavirus isolated from bank voles (Clethrionomys glareolus)

A previously unknown picornavirus was isolated from bank voles (Clethrionomys glareolus). Electron microscopy images and sequence data of the prototype isolate, named Ljungan virus, showed that it is a picornavirus. The amino acid sequences of predicted Ljungan virus capsid proteins VP2 and VP3 were closely related to the human pathogen echovirus 22 (approximately 70% similarity). A partial 5' noncoding region sequence of Ljungan virus showed the highest degree of relatedness to cardioviruses. Two additional isolates were serologically and molecularly related to the prototype.     

Molecular characterization of human enteroviruses in clinical samples: comparison between VP2, VP1, and RNA polymerase regions using RT nested PCR assays and direct sequencing of products

Three nested RT-PCR assays were developed to permit sensitive typing of enteroviruses directly from clinical samples. These assays amplified short fragments from different genomic regions codifying for three proteins: VP2, VP1, and RNA polymerase. Given that enteroviruses have a high rate of degeneration within target codons among serotypes, the primers used consisted of mixed base and deoxyinosine residues. These techniques detected at 0.03-0.003 TCID50 of prototype Poliovirus 1 and Echovirus 30. They were used to characterize the enteroviral RNA detected in 18 CSF, stool, and throw swab samples and in 8 enterovirus isolates from patients with several syndromes. Phylogenetic analysis in each independent sequenced region grouped the enterovirus into four clusters, enabling genetic classification. A comparative study was performed among the 26 sequences obtained after direct sequencing of products with those available in the nucleotide databases. The efficiency of each assay for enterovirus identification was evaluated by both distance (Clustal) and similarity (M-NW) indices. Comparative results obtained independently in the three regions showed the highest yield of correlation between nucleotide sequences of all prototype serotypes and the analyzed genotypes in the VP1 region (26/26, 100% Clustal; 22/26, 85% M-NW). Conversely, the VP2 region failed to identify some of the circulating enteroviruses (17/26, 65% Clustal; 16/26, 62% M-NW). Using the RNA polymerase region, sequences from samples and isolates were associated with prototype strains whenever these were available (20/21, 95% Clustal; 12/21, 57% M-NW). These assays were useful for molecular identification of enterovirus directly from samples even when isolation was not possible.     

A novel sapelovirus-like virus isolation from wild boar

A novel sapelovirus-like virus was isolated from a wild boar (Sus scrofa). In this study, partial viral genomic nucleotide sequences were determined using the rapid determination system of viral nucleic acid sequences (RDV) ver. 3.1, which we recently developed for discovering novel viruses. Phylogenetic analysis of VP1 and 3A proteins and their encoding nucleotide sequences of enteroviruses and sapeloviruses indicated that the isolated virus was closely related to porcine sapelovirus. RT-PCR detected viral sequences in six of 48 wild boar fecal samples.     

Complete genomic sequence of a Tobacco rattle virus isolate from Michigan-grown potatoes

Tobacco rattle virus (TRV) causes stem mottle on potato leaves and necrotic arcs and rings in potato tubers, known as corky ringspot disease. Recently, TRV was reported in Michigan potato tubers cv. FL1879 exhibiting corky ringspot disease. Sequence analysis of the RNA-1-encoded 16-kDa gene of the Michigan isolate, designated MI-1, revealed homology to TRV isolates from Florida and Washington. Here, we report the complete genomic sequence of RNA-1 (6,791 nt) and RNA-2 (3,685 nt) of TRV MI-1. RNA-1 is predicted to contain four open reading frames, and the genome structure and phylogenetic analyses of the RNA-1 nucleotide sequence revealed significant homologies to the known sequences of other TRV-1 isolates. The relationships based on the full-length nucleotide sequence were different from than those based on the 16-kDa gene encoded on genomic RNA-1 and reflect sequence variation within a 20–25-aa residue region of the 16-kDa protein. MI-1 RNA-2 is predicted to contain three ORFs, encoding the coat protein (CP), a 37.6-kDa protein (ORF 2b), and a 33.6-kDa protein (ORF 2c). In addition, it contains a region of similarity to the 3′ terminus of RNA-1, including a truncated portion of the 16-kDa cistron. Phylogenetic analysis of RNA-2, based on a comparison of nucleotide sequences with other members of the genus Tobravirus, indicates that TRV MI-1 and other North American isolates cluster as a distinct group. TRV M1-1 is only the second North American isolate for which there is a complete sequence of the genome, and it is distinct from the North American isolate TRV ORY. The relationship of the TRV MI-1 isolate to other tobravirus isolates is discussed.     

Detection of enterovirus genome sequence from diarrheal feces of goat

Goat diarrheal feces were subjected to metagenome analysis by the next-generation sequencing. Nucleotide sequences with homology to enteroviruses were obtained. Primers for RT-PCR were designed based on the nucleotide sequence of these sequences at the 5′-untranslated region, and we determined 563 bp nucleotide sequences that showed homology to bovine-like and ovine enteroviruses (77–87 %). We named the virus detected in this study goat enterovirus G1 (GEV-G1). In the phylogenetic analysis, GEV-G1 belonged to a cluster containing ovine enteroviruses. To our knowledge, this is the first report on nucleotide sequences of an enterovirus infecting Japanese goats.     

Complete genome sequence of arracacha mottle virus

Arracacha mottle virus (AMoV) is the only potyvirus reported to infect arracacha (Arracacia xanthorrhiza) in Brazil. Here, the complete genome sequence of an isolate of AMoV was determined to be 9,630 nucleotides in length, excluding the 3′ poly-A tail, and encoding a polyprotein of 3,135 amino acids and a putative P3N-PIPO protein. Its genomic organization is typical of a member of the genus Potyvirus, containing all conserved motifs. Its full genome sequence shared 56.2 % nucleotide identity with sunflower chlorotic mottle virus and verbena virus Y, the most closely related viruses.