Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2007, p. 2144-2150, Vol. 45, No. 7
0095-1137/07/$08.00+0 doi:10.1128/JCM.00174-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Clinical Investigation Facility, David Grant USAF Medical Center, 101 Bodin Circle, Travis Air Force Base, California 94535,1 Viral and Rickettsial Disease Laboratory, California Department of Health Services, 850 Marina Bay Parkway, Richmond, California 94804,2 Laboratory of Experimental Virology, Department of Medical Microbiology, Center for Infection and Immunity Amsterdam, Academic Medical Center, University of Amsterdam, 1105AZ Amsterdam, The Netherlands,3 D. I. Ivanovsky Institute of Virology, Russian Academy of Medical Sciences, 123098 Moscow, Russia4
Received 23 January 2007/ Returned for modification 20 February 2007/ Accepted 17 April 2007
|
|
|---|
|
|
|---|
Two species of viruses are distinguished within the Cardiovirus genus in the Picornaviridae family of viruses. One is the Theiler viruses, which includes Theiler's murine encephalomyelitis virus (TMEV), Theiler-like virus, and Vilyuisk virus (18-20). TMEV primarily infects mice, and Theiler-like virus was isolated from a rat (18). Vilyuisk virus, believed to be responsible for a neurodegenerative disease in the Yakuts people in the Vilyuy River area in Siberia, primarily infects humans (15). Interestingly, Vilyuisk virus, the only cardiovirus suspected to cause disease in humans, is not fully sequenced. The second species is the encephalomyocarditis viruses (EMCVs), which includes mengovirus, EMCV, Maus Elberfeld virus, and Columbia SK virus (4, 16). In addition to mice, these cardioviruses also infect humans, pigs, horses, elephants, squirrels, and monkeys (1, 6, 7, 9, 12).
In this study we sought to characterize a cardiovirus that was isolated from the stool of a patient who presented with FUO. Here, we used DNase sequence-independent single-primer amplification (DNase-SISPA) to obtain genetic information about the virus. Nearly full-genome sequencing and phylogenetic analysis was also performed.
|
|
|---|
Cell culture. The virus was first cultured in WI-38 cells and passed to HFDL-645 cells obtained from the State of California Viral and Rickettsial Disease Laboratory (VRDL) at the San Diego County Public Health Department Laboratory. L-645 is a line of fetal human diploid fibroblasts started and maintained by the VRDL. The sample was subsequently sent to the VRDL for further characterization. There, it was grown in in-house human fetal diploid kidney (HFDK) cells but failed to grow in primary monkey kidney, A-549, BSC, and RD cells. The virus grew in suckling mice. We did not harvest sera from any mouse for sequencing or any other analysis of Saffold virus (SAF-V).
Further characterization was carried out at the VRDL. By electron microscopy the agent had a diameter of 28 to 30 nm and appeared as a typical picornavirus. It was acid stable, was not neutralized by pooled human gamma globulin, and was temperature sensitive (no growth at 33°C).
DNase-SISPA. We screened cell-free medium that was cultured from HFDK cells. One hundred microliters of cell-free medium was filtered through a 0.22-µm-pore-size filter (Ultrafree MC; Millipore, Bedford, MA) and treated with 250 units of DNase I (Roche Diagnostics, Mannheim, Germany) to remove contaminating double-stranded DNA. DNase I-resistant nucleic acids were purified using a MagNA Pure LC DNA Isolation Kit I (Roche Diagnostics, Mannheim, Germany). Single-stranded DNA was created as follows: 50 µl of extracted nucleic acid was incubated with 240 pmol of primer A (5'-GTTTCCCAGTCACGATCNNNNNNNNN-3') and 40 units of RNase I for 5 min at 65°C. The mixture was cooled to 25°C, and the following components were added: 7.5 µl of Stratascript buffer, 7.5 µl of 0.1 M dithiothreitol, 1 µl of 10 mM deoxynucleoside triphosphate (dNTP) mix, and 800 units of Stratascript reverse transcriptase (Stratagene, La Jolla, CA). The reaction mixture was incubated at 42°C for 30 min. Double-stranded DNA was generated by incubating 44 µl of single-stranded DNA with 5 units of Klenow fragment (exo) (New England Biolabs, Beverly, MA) and 5 µl of 10x EcoPol buffer (New England Biolabs, Beverly, MA) at 10°C for 5 min, followed by 8 min at 37°C. Twenty-five microliters of double-stranded DNA was subsequently amplified by 5 units of Takara Ex Taq (Takara Bio, Shiga, Japan) with 100 pmol of primer B (5'-GTTTCCCAGTCACGATC-3') in 1x Takara buffer, 2 mM MgCl2, and 0.8 mM dNTP with the following cycling conditions: 40 cycles of 94°C for 21 s, 40°C for 30 s, 50°C for 30 s, and 72°C for 1 min. Next, 44 µl of PCR product was digested with 5 units of the restriction enzyme Csp6.I (Fermentas, Hanover, MD) in 1x B+ buffer (Fermentas, Hanover, MD). Twenty picomoles of adaptors composed of the hybridized oligonucleotides NBam24 5'-AGGCAACTGTGCTATCCGAGGGAG-3' and NCsp11 5'-TACTCCCTCGG-3' was then ligated to 15.5 µl of the restriction enzyme-digested DNA using 200 units of T4 DNA ligase (Invitrogen, Carlsbad, CA) in 1x T4 DNA ligase buffer for 10 min at room temperature. Two microliters of the ligation reaction mixture was used as a template in a 50-µl PCR mixture containing 25 pmol of NBam24, 1x Takara buffer, 0.8 mM dNTPs, and 2.5 units of Takara Ex Taq. The PCR mixture was heated to 72°C for 3 min before the addition of Takara Ex Taq. The mixture was then heated for an additional 5 min at 72°C to generate the primer-binding sites. Cycling conditions were as follows: 94°C for 1 min and 72°C for 3 min for 40 cycles.
Analysis of amplified DNA. Amplified PCR products were analyzed by polyacrylamide gel electrophoresis (PAGE). The distinct DNA bands seen by PAGE were excised and pooled. DNA from the pooled PAGE bands from each sample were purified and subcloned into pCR-II-TOPO (Invitrogen, Carlsbad, CA). Because excised DNA bands were pooled prior to subcloning, we were unable to relate each sequence to a particular PAGE band. Plasmids containing inserts were purified and sequenced using the T7 sequencing primer. We sequenced our cloned and uncloned PCR fragments by using a Big Dye terminator, version 3.1, sequencing kit (Amersham Pharmacia Biotechnology) and an ABI Prism 3130x genetic analyzer (Applied BioSystems) according to the manufacturers instructions.
Amplification of the SAF-V genome. To amplify regions of SAF-V flanking the sequences that we found, we designed primers based on an alignment of TMEV (X56019), Theiler-like virus (AB090161), mengovirus (L22089), and EMCV (M81861), using Primer Premier, version 5.0, software (Palo Alto, CA). All amplicons were then sequenced using primer walking. The 5' end of SAF-V was found using a kit for the 5' rapid amplification of cDNA ends (Invitrogen) according to the manufacturer's instructions.
Phylogenetic analysis of the SAF-V genome. Phylogenetic analysis was performed using sequences representing full-length genomes of virus species representing all genera of the Picornaviridae. Sequences of all recognized individual virus species were included in the analysis except for the members of the Enterovirus genus, for which participation in the analysis was limited to the two most divergent members of the genus for presentation purposes. Picornavirus sequences for echovirus 5 (accession number AF083069), human rhinovirus 1B (D00239), porcine enterovirus 8 (AF406813), human hepatitis virus type A (M20273), simian hepatitis A virus (D00924), avian encephalomyelitis virus (AJ225173), Ljungan virus (AF327921), human parechovirus 2 (AJ005695), echovirus 22 (L02971), echovirus 23 (AF055846), mengovirus (L22089), EMCV (X87335), TMEV (M20562), Theiler-like virus NGS910 (AB090161), foot-and-mouth disease virus (AF308157), equine rhinovirus type 1 (X96870), Aichi virus (AB040749), porcine teschovirus (AF231769), porcine enterovirus 1 (AJ011380), equine rhinovirus 3 (AF361253), and equine rhinovirus type 2 (X96871) were retrieved from the GenBank and aligned using the CLUSTALX program (24). A neighbor-joining tree based on nucleotide distances with the Jukes-Cantor correction and pair-wise gap deletions with bootstrap resampling (100 replicates) was constructed using the MEGA3 software (13). More detailed phylogenetic analysis was performed for each genomic region (5' and 3' untranslated regions [UTRs], leader peptide, and the proteins VP1 to VP4, 2A to 2C, and 3A to 3D) of cardioviruses. For the Vilyuisk virus, no full-length genome sequence is available. Therefore, we used sequence M94868 in the analysis of homologous regions. Noncoding regions of cardioviruses were aligned at the nucleotide level, and coding regions were aligned at the amino acid level and then back-translated to the nucleotide sequences using the CLUSTALW program as implemented in the BioEdit program (8). For each region, phylogenetic trees were built with bootstrap resampling (100 replications) using the MEGA3 program and based on nucleotide p-distances for noncoding regions and on amino acid p-distances for coding regions. Nucleotide and amino acid similarities between the sequences of the SAF-V and other cardioviruses were calculated by the MEGA3 program. To check whether the SAF-V is a recombinant of known viruses, recombination analysis was performed using the bootscan method of the SimPlot program (http://sray.med.som.jhmi.edu/SCRoftware/).
Nucleotide sequence accession number. The sequence of the SAF-V has been deposited in the GenBank database under accession number EF165067.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Growth of SAF-V in tissue culture and suckling mice
|
![]() View larger version (13K): [in a new window] |
FIG. 1. DNase-SISPA amplification of tissue culture sample containing SAF-V. PCR products were analyzed on a 6.5% acrylamide gel. Lane L contains molecular size markers in base pairs.
|
Molecular characteristics of SAF-V. The genetic organization of SAF-V is similar to that of other cardioviruses since it has one open reading frame (ORF) that has 13 cleavage sites in its genome, which should separate the viral polypeptide into the following proteins: L, VP4, VP2, VP3, VP1, 2A, 2B, 2C, 3A, 3B, 3C, and 3D (Fig. 2).
![]() View larger version (12K): [in a new window] |
FIG. 2. Schematic representation of SAF-V (A), Theiler virus (B), TMEV (C), and EMCV (D) genomes. The boxes represent the genes that encode proteins. Solid black horizontal lines represent the UTRs. Yellow boxes represent the genes that encode the L* proteins.
|
|
View this table: [in a new window] |
TABLE 2. Nucleotide identities between the sequences of the SAF-V and other cardioviruses
|
|
View this table: [in a new window] |
TABLE 3. Amino acid identities between the sequences of SAF-V and other cardioviruses
|
![]() View larger version (42K): [in a new window] |
FIG. 3. Alignment of the cardiovirus leader peptides. The zinc finger domain is highlighted in gray, the acidic domain is highlighted in aqua, and the Ser/Thr domain is highlighted in yellow.
|
The EF loop structure of cardioviruses, associated with virulence, is located in the VP2 protein. Specifically, there are two loops within the larger loop structure. The first part of the SAF-V EF loop structure is unique since it has only 2 amino acid identities to other cardioviruses (Fig. 4). In contrast, the second loop of SAF-V has 7 of 18 (39%) amino acid identities to the second loops of Theiler-like virus and Vilyuisk virus (Fig. 4). Interestingly, there is an insertion in the middle of the second loop of the SAF-V EF loop, which is not found in any other cardiovirus VP2 protein (Fig. 4). The CD loop structure of cardioviruses, also associated with virulence, is located in the VP1 protein. The first of the two SAF-V CD loops has no amino acid identity to other cardioviruses, and the second half of the loop is deleted (Fig. 5). The second CD loop of Theiler-like virus and Vilyuisk virus has 47% (8/17) and 39% (7/17) amino acid identity to the second CD loop of SAF-V, respectively. Taken together this suggests that the structure of SAF-V is unlike all known cardioviruses.
|
View larger version (17K): [in a new window] |
FIG. 4. Alignment of the EF loop structure amino acids in VP2 of cardioviruses SAF-V, TMEV, Theiler-like virus, Vilyuisk virus, EMCV, and mengovirus. Areas highlighted in gray are part of the EF loop structure. The area highlighted in yellow is an addition to the CD loop structure of SAF-V. Amino acids shown in blue are identical to the amino acids in SAF-V.
|
|
View larger version (20K): [in a new window] |
FIG. 5. Alignment of the CD loop amino acids in VP1 of cardioviruses SAF-V, TMEV, Theiler-like virus, Vilyuisk virus, EMCV, and mengovirus. Areas highlighted in yellow are part of the CD loop structure. Amino acids shown in red are identical to the amino acids in SAF-V.
|
![]() View larger version (14K): [in a new window] |
FIG. 6. Evolutionary relationships among Picornaviridae based on full-length genomes of picornaviruses (A) and based on sequences of each genomic region of cardioviruses (B). GenBank accession numbers are given in Materials and Methods. Abbreviations used for virus names: HEchoV 5, 22, and 23, human echoviruses 5, 22, and 23; HRV 1B, human rhinovirus 1B; PEV 1 and 8, porcine enterovirus 1 and 8; HAV, human hepatitis A virus; SHAV, simian hepatitis A virus; AEMV, avian encephalomyelitis virus; LjunganV, Ljungan virus; HPEchoV 2, human parechovirus 2; MengoV, mengovirus; T lV, Theiler-like virus; FMDV, foot-and-mouth disease virus; ERV 1, 2, and 3, equine rhinovirus 1, 2, and 3; AichiV, Aichi virus; PTV, porcine teschovirus.
|
![]() View larger version (36K): [in a new window] |
FIG. 7. Bootscan analysis of SAF-V based on full-length genomes. For virus names, see the legend of Fig. 6.
|
|
|
|---|
Molecular evolutionary analysis demonstrated that SAF-V is a new member of the Cardiovirus genus in the Picornaviridae family (Fig. 6 and 7). In the phylogenetic analysis of the full-length genome sequences, SAF-V appeared to be most closely related to Theiler-like virus, followed by TMEV of mice and rats (Fig. 6A). The only human virus related to the Theiler viruses identified so far was the Vilyuisk virus (15). This virus, for which only short genetic regions have been sequenced, was known to be the most divergent virus within the Theiler group of viruses (Fig. 6B). Our analysis of the genetic region for which amino acid sequence of the Vilyuisk virus is available demonstrated that SAF-V is closely related to the Theiler's group (Fig. 6B). The bootstrap analysis of SAF-V demonstrated that this virus is not a recombinant of other known picornaviruses (Fig. 7).
Currently, two virus species are recognized among cardioviruses: EMCVs (mengovirus and EMCV, for which full-genome sequences are available, as well as Columbia SK and Maus Elberfeld viruses, which are not sequenced) and the Theiler virus species (TMEV, Theiler-like virus, and Vilyuisk virus) (23). The demarcation criteria in the genus are that the members of a species share greater than 70% amino acid identity in the VP1, 2C, 3C, and 3D coding regions and that they share a natural host range. While our results demonstrated that the SAF-V is related to the Theiler virus species, its amino acid similarity in the VP1 region to the recognized members of the Theiler virus species is below 70% (Table 2).
The genetic organization of SAF-V is very similar to other cardioviruses. The SAF-V structural proteins VP1, VP2, VP3, and VP4 are in the same locations as those of TMEV (Fig. 2). In addition, they are approximately the same size. Specifically, the loci of genes in SAF-V are more similar to Theiler virus species than the EMCV species (Fig. 2), which is not surprising in light of the phylogenetic data. Since SAF-V does not have a poly(C) sequence in its 5' UTR, which is associated with virulence in EMCV (13), we would not expect SAF-V to be virulent. Conversely, the absence of a poly(C) sequence in the 5' UTR does not preclude SAF-V virulence (14, 26).
The four loop structures that form the exterior of the cardiovirus surface are two small loops within the EF loop structure of VP2 and the two loops of the CD loop structure of VP1 (10). The EF loop structure of SAF-V is unique since the second EF loop contains a 4-amino-acid insertion, and both loops combined contain no more than 26% amino acid identity to the EF loops of other cardioviruses (Fig. 4). Similarly, the CD loop structure of SAF-V also has low homology to all other cardioviruses (0% in the first loop and 47% in the second to Theiler-like virus) as well as a 4-amino-acid deletion in the first loop. In contrast, TMEV has 63 and 44% amino acid identity to the first and second CD loops of Theiler-like virus, respectively. This suggests that the viral surface of SAF-V has a novel structure and possibly a different mechanism for infection. Furthermore, the primary sequences of the EF and CD loop structures in SAF-V do not offer any clues as to whether SAF-V is virulent like the GDVII strain of TMEV or persistent like TMEV strain DA.
Transmission of TMEV and EMCV occurs by the fecal-oral route (3, 22). We do not know whether SAF-V caused the patient's fever; therefore, it is also difficult to hypothesize how or when the patient was infected with SAF-V. Since TMEV is pathogenic in mice, further investigation is under way to determine the prevalence of SAF-V in humans.
Because of the length of time between the illness and our molecular characterization, we do not know much about the patient's health around the time of initial infection (e.g., characteristics of the illness and/or length of symptoms beyond the 24th day) other than the fact that the stool sample was collected 24 days after the onset of FUO, which fits the clinical definition of FUO (21). In addition, we do not know what other pathogens were tested for. However, SAF-V was the only agent isolated from the stool of the study subject. Nevertheless, we cannot claim that SAF-V is the etiological agent because we have not analyzed fecal samples from patients with FUO of known and unknown etiology.
We identified and characterized a new human virus, SAF-V, which was cultured from the stool of a pediatric patient presenting with FUO. SAF-V represents a new member of the genus Cardiovirus in the Picornaviridae family. The prevalence and pathogenicity of this virus remain to be determined.
The views expressed in this material are those of the authors and do not reflect the official policy or position of the U.S. Government, the Department of Defense, or the Department of the Air Force.
Published ahead of print on 25 April 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»