Previous Article | Next Article ![]()
Journal of Virology, June 2008, p. 5408-5416, Vol. 82, No. 11
0022-538X/08/$08.00+0 doi:10.1128/JVI.00070-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Cincinnati Children's Hospital Medical Center,1 University of Cincinnati College of Medicine, Cincinnati, Ohio,2 Tulane National Primate Research Center, Covington, Louisiana,3 Tulane University School of Medicine, New Orleans, Louisiana4
Received 10 January 2008/ Accepted 21 March 2008
|
|
|---|
|
|
|---|
CVs are small, nonenveloped, icosahedral viruses with a positive-sense, single-stranded, polyadenylated RNA genome (of about 7.5 to 8.5 kb). The Lagovirus, SV, and Becovirus genomes are organized into two open reading frames (ORFs); ORF1 encodes the nonstructural (NS) proteins and the capsid protein, while the NV and Vesivirus genomes contain three ORFs, where the capsid protein is encoded by ORF2, a separate ORF. The translated proteins of all CVs contain several conserved amino acid (aa) motifs that are similar to those of picornaviruses (23): 2C helicase/NTPase (GXXGXGKS/T), 3C protease (GDCG), and RNA-dependent RNA polymerase (RdRp; GLPSG and YGDD). The ORF closest to the 3' end of the CV genome (ORF2 for Lagovirus, SV, and Becovirus and ORF3 for NV and Vesivirus) encodes a basic protein (pI,
10.0) that has been reported to be a minor structural protein for feline CV and the Norwalk virus (10, 34).
The only CV (Pan-1) that has been isolated from nonhuman primates belongs to the Vesivirus genus (31). Several studies have been conducted for the development of a primate animal model for studying HuCV disease. Chimpanzees inoculated with the Norwalk virus developed subclinical infections, with seroresponses and virus shedding but no vomiting or diarrhea (41). In a more recent study, Subekti et al. reported the development of clinical illness characterized by diarrhea, dehydration, vomiting, and viral shedding in newborn pigtail macaques inoculated with the Toronto virus (35). Recently, Rockx et al., evaluated the susceptibility of four different Old World monkey species to NVs (29) and found that rhesus macaques might be susceptible to NV infection, since one of the three animals inoculated with the Norwalk virus shed the virus for up to 19 days postinoculation and developed specific immunoglobulin M and immunoglobulin G responses. In the above-mentioned study, serological evidence of natural NV infection could not be demonstrated among different monkey species; however, Jiang et al. previously reported a high prevalence of anti-NV antibodies (both GI and GII) in serum samples collected from nonhuman primate species, including rhesus macaques (14).
In this study, we report the isolation and characterization of a novel, tissue culture propagable CV, the Tulane virus (TV), which was isolated from stool samples of rhesus macaques and represents a new genus within the Caliciviridae.
|
|
|---|
Detection of CV RNA in stool specimens by RT-PCR. Total RNA was extracted from 100 µl of 10% (wt/vol) stool specimens with the TRIzol reagent (Gibco BRL, Gaithersburg, MD) according to the manufacturer's instructions. Precipitated RNA was resolved in 20 µl of molecular-biology-grade water (Eppendorf AG, Hamburg, Germany), and 1 µl was used as the template in the RT-PCR. The remainder of the samples was stored at –70°C for repeated assays. RNA was amplified in 25-µl reaction mixtures by using the AccessQuick RT-PCR system (Promega, Madison, WI) according to the manufacturer's protocol. Primers P289H,I/P290H,I,J,K targeting conserved sequences in the RdRp region were used (9, 15). These primers are able to amplify both NVs and SVs. RT was performed at 48°C for 45 min, followed by a 3-min denaturation step at 94°C and 35 PCR cycles of denaturation for 30 s at 94°C, annealing for 1 min at 49°C, and elongation for 1 min at 72°C, followed by a final extension of 10 min at 72°C. Twelve microliters of each reaction mixture was analyzed on 1% agarose gels in the presence of ethidium bromide.
Amplification and sequence analysis of the TV genome. The primers used in this study are listed in Table 1. Based on the sequence data obtained from the short RdRp amplicon, a TV-specific primer (P776) was designed and used with an oligo(dT) primer to obtain an approximately 3-kb amplicon, including the 3' end of ORF1, the complete ORF2 and ORF3, and the 3' nontranslated region of the genome. Extended amplification of the genome toward the 5' end has been achieved by TV-specific reverse primers that were designed based on newly obtained sequences and degenerated forward primers targeting the GDCG (P774) and the GXXGXGKT amino acid motifs (P802) conserved in CV proteases and NTPases, respectively. The 5'-end sequence of the genome was determined by the 5' rapid amplification of cDNA ends system (Invitrogen, Carlsbad, CA) according to the manufacturer's recommendations by using TV-specific primers in combination with the abridged anchor primer and abridged universal primer and viral RNA extracted from tissue-culture-adapted, CsCl-purified TV as the template. Both homopolymeric deoxycytidine-tailed and deoxyadenine-tailed cDNAs were made to determine the exact 5' end. Additionally, the end ligation of proteinase-K-treated viral RNA with T4 RNA ligase I (New England Biolabs, Ipswich, MA) was also performed according to the method described previously by Mandl et al. (20), with slight modifications. Briefly, 1 µg extracted viral RNA was digested with proteinase K (100 µg/ml proteinase K, 30 mM Tris-Cl [pH 8.0], 30 mM EDTA, 0.5% sodium dodecyl sulfate [SDS]) at 55°C for 1 h, recovered by using an RNeasy mini column (Qiagen, Valencia, CA), and ligated with T4 RNA ligase I (New England Biolabs, Inc., Ipswich, MA) at 4°C overnight. The ligated RNA was recovered on an RNeasy mini column and used for cDNA synthesis with reverse primer P1080. First-round PCRs were performed with primers P1068 and P1080, and a second-round PCR was performed with primers P1035 and P1082, which flanked the ligated 5'-3' end junction of the genome. The PCR yielded a smeared product, with fragments of approximately 400 to 500 bp, which were gel purified and ligated into the pGEM-T vector (Promega, Madison, WI). Clones were PCR screened with M13 R/F primers, and clones with the longest inserts were selected for sequencing.
|
View this table: [in a new window] |
TABLE 1. List of primers used in the study
|
Phylogenetic analysis. Sequences (published in the GenBank) used in the phylogenetic analysis were of the following full-length CVs (accession nos. are shown in parentheses): for NV, Norwalk/68 (M87661), Southampton (L07418), MD145 (AY032605), and MNV-1 (DQ285629); for SV, Manchester/93 (X86560), and PEC-Cowden (AF182760); for Vesivirus, SMSV-1 (U15301), FCV-CFI68 (U13992), and FCV-Urbana (NC_001481); for Lagovirus, EBHSV-GD (NC_002615) and RHDV-FRG (M67473); and for Becovirus, BEC-NB (NC_004064) and Newbury-1 (DQ013304). Multiple alignments of nucleotides (nt) and amino acid sequences were created using the Omiga 2.0 software (Oxford Molecular Ltd, Oxford, United Kingdom). Alignments for individual NS proteins used in the analysis were created based on an initial alignment of the NS polyprotein and analysis for the presence of amino acid sequence motifs within the published proteolytic cleavage sites of other CVs (1, 19, 22, 24, 25, 32, 33). Individual NS proteins were realigned and edited to select partial alignments with relatively lower diversities. Bordering amino acid positions for the partial viral NTPase, protease, and polymerase sequences used in the final alignments are listed in Table 2. To retain the maximum number of conserved evolutionary sites, columns containing gaps were removed and the remaining consensus length alignments were used to calculate pairwise homology scores (percent) and construct phylogenetic trees. Dendrograms based on amino acid alignments were constructed by the unweighted pair group method with arithmetic mean and the neighbor-joining clustering methods of the Molecular Evolutionary Genetics Analysis (MEGA, version 3.1) software with Poisson correction distance calculations. The confidence values of the internal nodes were obtained by performing 125 and 1,025 bootstrap analyses, which yielded similar tree topologies and only slight differences (1 or 2%) in the branch bootstrap values. Pairwise distances and means within (intragenogroup) and between (intergenogroup) genogroups and clusters were also calculated.
|
View this table: [in a new window] |
TABLE 2. Partial NTPase, protease, and polymerase regions used in multiple sequence alignmentsa
|
Cells and medium were harvested 5 days postinoculation by scraping, transferred to sterile 15-ml tubes (Becton Dickinson and Company, Franklin Lakes, NJ), and freeze-thawed at –20°C to release cell-associated virus. Cell debris was removed by centrifugation at 2,000 rpm (600 x g) for 5 min, and 200 µl of the cell-free medium was transferred to fresh cell cultures, which were monitored daily for cytopathic effect by light microscopy. These steps were repeated for six passages. After the 4th and 6th passages, 100 µl medium from each cell line was tested for the presence of TV by RT-PCR using TV-specific primers. Later, additional human cell lines (293T, intestine [int] 407, FHs74 int.) were also tested for their permissibility.
Plaque purification. LLC-MK2 cell cultures grown to 80% confluence in six-well culture plates (Corning Incorporated Life Sciences, Lowell, MA) were incubated for 1 h with 1 ml of serially diluted chloroform-extracted and CsCl-purified TV obtained from the 4th passage. Following inoculation, cells were washed with medium and overlaid with 2 ml of agarose-medium containing 1 volume of 1.2% agarose in distilled water and 1 volume of 2x M199 medium supplemented with 10% fetal bovine serum. Plaque formation was evaluated daily by light microscopy. Single plaques were picked from wells with one to five isolated plaques, and agarose plugs were dissolved in 500 µl medium and used to inoculate new plates. The process was repeated six times. Viruses grown from nine final plaques were stocked at –70°C. In all subsequent experiments, plaque-purified TV M33 3/2 has been used.
Cesium chloride density gradient purification of TV. LLC-MK2 cultures grown to 90% confluence in roller bottles (BD Biosciences, San Jose, CA) and containing 100 ml culture medium were inoculated with 3 ml of TV inoculum (105.3 50% tissue culture infective doses per milliliter). At 72 h postinoculation, cultures were freeze-thawed once and attached cells were removed by moving the ice around the flask. Cell debris was separated from culture medium in a JA-10 rotor at 8,000 rpm (10,000 x g) for 20 min in a Beckman J2-21 centrifuge. The cell pellet was resuspended in culture medium, freeze-thawed three additional times to free entrapped virus, and centrifuged at 8,000 rpm (10,000 x g) for 20 min, and the supernatant was collected and added to the culture medium. Virus particles were precipitated by mixing 1 volume of a 30% polyethylene glycol 8000-0.5 M NaCl (pH 7.4) stock solution with 4 volumes of culture medium and stirring at room temperature for 2 h and at 4°C overnight. Polyethylene glycol precipitate was collected by centrifugation at 10,000 x g for 45 min and resuspended in 1/100 of the original culture volume in phosphate-buffered saline (PBS). The resulting material was mixed with CsCl to obtain a solution with a density of 1.34 g/ml by adding 5.11 g CsCl to 10 ml suspension. Gradients were formed by centrifugation at 35,000 rpm (150,000 x g) for 27 h in a SW 41Ti rotor using a Beckman Coulter Optima L-90 K ultracentrifuge. Gradients were fractionated by bottom puncture, and the density of each fraction was determined by a traditional Abbe refractometer. Viral RNA was extracted from 1 µl of each fraction and analyzed for the presence of TV RNA by RT-PCR. Peak fractions were diluted five times in PBS, and virus particles were collected by centrifugation at 35,000 rpm (150,000 x g) for 3 h in a SW 41Ti rotor, resuspended in PBS, and stored at –70°C.
Electron microscopy (EM). Drops of CsCl-gradient-purified TV were adsorbed onto Formvar-carbon 200-mesh copper grids (Electron Microscopy Sciences, Fort Washington, PA), stained with 1% ammonium-molybdate, and examined using a Hitachi-7600 transmission electron microscope equipped with an AMT digital camera. The mean diameter of TV particles was determined by measuring 25 virus particles at a direct magnification of x100,000 to x200,000.
|
|
|---|
Genomic organization of TV. The complete genome of the TV is 6,714 nt long without the poly(A) tail, making it the shortest known CV genome. It has a ribonucleoside composition of 27.46% A, 24.58% C, 22.07% G, and 25.87% U residues and starts with three G residues that are part of a 14-nt-long 5' untranslated region (UTR) (GGGTGACTAGAGCTATG). The genome contains three large ORFs. ORF1 spans from nt 15 to 4358, ORF2 from nt 4380 to 5984, and ORF3 from nt 5984 to 6640 (Table 3 and Fig. 1A). ORF1 and ORF2 are in frame +2 relative to the first nucleotide of the genome and separated by 21 nt, which is different from the case for NVs, where ORF1 and ORF2 are in different frames and overlap by 11 to 17 nt, but similar to the case for canine CV (GenBank accession no. NC_004542), a vesivirus, where ORF1 and ORF2 are in the same frame and are separated by 3 nt. The TV ORF3 is in frame +3 and overlaps with the last nucleotide of the ORF2 termination codon (TAATG; bold letters indicate the termination [stop] codons of ORF2; underlining the initiation [start] codons of ORF3), which is similar to NVs and different from vesiviruses and SVs (ATGA), lagovirus (ATGTCTGA), and becoviruses (TAATATG). For all three ORFs, the predicted initiation codons are in a strong consensus context for translation initiation of eukaryotic mRNA (17). The 5' UTR did not reveal significant homologies with the 5' region of the putative subgenomic RNA, which is different from the case for NVs; however, it aligned well to a region 21 nt downstream from the ORF2 initiation codon starting at nt 4404 (Fig. 1B). The TV genome ends with a 74-nt-long 3' UTR and a poly(A) tail of 26 to 62 adenine residues, as determined by the sequencing of 12 clones obtained from end-ligated genomic RNA.
|
View this table: [in a new window] |
TABLE 3. Comparison of the Tulane virus with representative caliciviruses of the other genera
|
![]() View larger version (15K): [in a new window] |
FIG. 1. (A) Genome organization of TV. Conserved amino acid motifs are indicated. (B) Alignment of the 5' ends of genomic and subgenomic TV RNAs. Initiation codons of the NS polyprotein and VP1 are boxed.
|
|
View this table: [in a new window] |
TABLE 4. Amino acid identity of the Tulane virus with representative caliciviruses of the other genera
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Predicted NS polyprotein cleavage sites of TV based on sequence alignments with published cleavage maps of other CVs. Since the alignments revealed the highest identity between TV and NVs, only published NV cleavage maps representing GI, GII, and GV are shown. Calculated protein molecular masses (in kilodaltons) and the positions of the amino acids at the cleavage sites are shown.
|
The TV VP2 is 218 aa long, with a calculated molecular mass of 22.8 kDa and an isoelectric point (pI) of 10.0, like the case for NVs. It also shows the highest amino acid identity (22 to 26%) with NVs (Table 4).
Phylogenetic analysis.
Pairwise homology scores from amino acid sequence alignments indicated that TV exhibits the highest amino acid identity with NVs in all regions analyzed (Table 4). This was confirmed by phylogenetic analyses conducted on the viral NS polyprotein, partial NTPase, protease, polymerase, and whole-capsid sequences (Fig. 3). VP2, Vpg, P16, and the N-terminal protein were not analyzed, because of the higher diversity of these proteins between the different CV genera. Phylogenetic trees obtained for each genomic region grouped viruses of the previously described genera in a manner similar to that in the trees reported by others, and TV was consistently placed on a branch rooting with NV. The mean distances between TV and the different genera also showed that TV is most closely related to NV (the mean distances between TV and NV were as follows: for ORF1 polyprotein,
was 1.459 and
M was 0.041; for NTPase,
was 0.791 and
M was 0.080; for protease,
was 1.317 and
M was 0.162; for polymerase,
was 1.124 and
M was 0.065; and for VP1,
was 1.442 and
M was 0.072.
indicates mean distance, and
M indicates standard error of the mean); however, these distances were longer than the distances between some of the established genera (the mean distances between SV and Vesivirus were as follows: for ORF1 polyprotein,
was 1.184 and
M was 0.030; for NTPase,
was 0.682 and
M was 0.068; for protease,
was 1.326 and
M was 0.153; for polymerase,
was 0.844 and
M was 0.047; and for VP1,
was 1.322 and
M was 0.067). Also the mean distances between TV and NV were 1.66 to 3.5 times longer than the mean distances within NV (intergenogroup distance) (the mean distances within NV were as follows: for ORF1 polyprotein,
was 0.550 and
M was 0.017; for NTPase,
was 0.444 and
M was 0.052; for protease,
was 0.367 and
M was 0.059; for polymerase,
was 0.392 and
M was 0.025; and for VP1,
was 0.864 and
M was 0.042). Taken together, these results indicate that TV represents a new genus within the Caliciviridae.
![]() View larger version (20K): [in a new window] |
FIG. 3. Unrooted phylogenetic trees based on amino acid sequence alignments of CV NTPase (A), polymerase (B), and VP1 (C). Trees were constructed by the neighbor-joining clustering method of MEGA 3.1, with Poisson distance calculations. The scale bars represent the phylogenetic distances expressed as units of amino acid substitutions per site. The confidence values of the internal nodes were obtained by performing 1,025 bootstrap analyses. Trees constructed for the NS polyprotein and protease exhibited similar topologies (data not shown).
|
![]() View larger version (39K): [in a new window] |
FIG. 4. Detection of TV-specific RNA at the 6th passage of a TV-positive stool sample in CaCo-2, MA-104, Vero, and LLC-MK2 cells. M, molecular size marker (1-kb DNA ladder; Invitrogen, Carlsbad, CA). Mock-infected cells were passed parallel with stool-inoculated cells for each cell line. RNA was extracted from 100 µl of cell-free tissue culture medium obtained from each mock or infected (stool) cell culture. RNA extracted from 100 µl of TV-positive stool that was used to inoculate the cells originally served as a positive control (+Contr). –Contr, negative control. Primers P885 and P886 (Table 1) were used to amplify a 523-bp TV-specific product.
|
Physicochemical characterization of TV. CsCl density gradient fractionation of plaque-purified TV revealed a clear peak, with a mean density (± standard deviation) of 1.3722 ± 0.014 g/ml containing small rounded virions with cup-like surface protrusions (Fig. 5). The mean diameter of TV particles is 35.8 ± 1.6 nm, which is similar to that of other CVs. SDS-polyacrylamide gel electrophoresis analysis of the purified TV revealed a single major protein of about 60 kDa, most likely the capsid protein (VP1) (Fig. 5).
![]() View larger version (79K): [in a new window] |
FIG. 5. Physicochemical properties of TV. (A) RT-PCR amplification of TV-specific RNA extracted from CsCl density gradient fractions revealed a peak with a mean density calculated as follows: = 1.37 g/ml. –Contr, negative control. (B) Virus particles with typical CV morphology and a mean diameter (calculated as = 35.8 nm) were visualized in the peak fractions. (C) A single protein band of about 60 kDa could be visualized by SDS-polyacrylamide gel electrophoresis of peak fractions.
|
|
|
|---|
The whole genome of the TV is 6,714 nt long, which is at least 600 nt shorter than the genomes of other CVs. To avoid the possibility of missing the exact 5' end, it was determined by several methods and all of them yielded the same result. Moreover, the transfection of LLC-MK2 cells with an in vitro-transcribed TV genomic RNA yielded infectious virions (data not shown), further confirming our data.
TV represents a new CV genus that is genetically more closely related to NVs than to other CVs. Phylogenetic analysis conducted on structural and NS proteins of TV and representatives of other CV genera consistently placed TV on a separate branch rooted together with NV (Fig. 3). This placement is in accordance with the pairwise homology and phylogenetic distance calculations that also indicated the closest relationship between TV and NVs. However, the low-homology scores observed in all regions and the equivalence of phylogenetic distances between TV and NV to distances between the established CV genera clearly place TV outside the NV genus. In accordance with the proposed nomenclature of Becovirus (26), we propose the name Recovirus (rhesus enteric Calicivirus) for this new genus represented by TV.
Additional features of TV with high similarity to those of NVs have been also noted. The lengths of the TV ORF2 and ORF3 and the Mrs of the encoded structural proteins VP1 and VP2 are comparable to the corresponding values for NVs (Table 3). Also, the junction between ORF2 and ORF3, which seems to be conserved for CVs of different genera, is similar to that of NVs (TAATG). The shorter genome length of TV is clearly due to a shorter ORF1 (Table 3). The TV ORF1 is 4,344 nt long, which is more than 1,000 nt shorter than the 5,370 nt of Norwalk virus ORF1. Based on alignments of deduced NS polyprotein sequences and the predicted protein cleavage map of TV (Fig. 2), the TV NS polyprotein is cleaved into at least six functional proteins, NH2-P26.4(N-terminal)-P38(NTPase)-P16-P10.8(Vpg)-P18(Prot)-P52.9(Pol)-COOH, all of which exhibit the highest amino acid identity with the corresponding proteins of NVs. The N-terminal protein, P16, and VPg are significantly (about 20 to 40%) smaller than those of NVs, while the remaining NS proteins have Mrs comparable to those of NVs.
One significant difference between TV and HuCVs is that TV could be readily grown in vitro. Whether this result is due to the difference in NS protein function or receptor usage or to other factors is unclear and needs to be studied. It also should be noted that although TV could be grown in LLC-MK2 cells (a rhesus macaque kidney cell line), growth attempts in other cell lines failed. Thus, investigating the nature of cellular and viral elements required for TV replication may provide important clues for identifying factors required for successful in vitro growth of HuCVs. The availability of a cell culture and a reverse genetics system makes TV a valuable model to address these issues.
Since TV was isolated from stool samples of animals without symptoms of diarrhea, the clinical disease linked to TV infection and its pathogenesis and whether TV can be developed to an animal model system for CV gastroenteritis remain to be established. A recent study reported the detection of NV-like particles in nonhuman primate stool samples, including rhesus macaques, by EM (38). Whether these particles represented TV or NVs is unknown. Taken together, these findings call for extended investigations of CVs in nonhuman primates. Previous attempts at the molecular detection of CVs in primate stool samples were unsuccessful, possibly due to primer selection. The detection of TV with the primer set P289H,I/P290H,I,J,K further emphasizes the usefulness of these primers in CV detection.
We thank Xiaoyun Deng, Weiming Zhong, and Mayra A. Cantu for technical support and Irene Hofmann for helping with the EM.
Published ahead of print on 2 April 2008. ![]()
|
|
|---|
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»