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

State Key Laboratory of Emerging Infectious Diseases, Department of Microbiology, The University of Hong Kong, Hong Kong SAR, China,1 HKU-Pasteur Research Centre, Hong Kong SAR, China2
Received 23 April 2008/ Accepted 3 June 2008
|
|
|---|
|
|
|---|
![]() View larger version (7K): [in a new window] |
FIG. 1. Schematic diagrams of the mink astrovirus (MAstV) genome and bat astrovirus (BatAstV) AFCD337 partial genome. ORFs, protease motif (Pro), and RdRp motif are shown in the diagram. The unsequenced putative 5' end of the genome region of BatAstV AFCD337 is represented by a dotted line. The 2.5-kb ORF1a (partial) region and the ORF1b and ORF2 regions used for phylogenetic analysis in Fig. 2 are indicated by arrows. An, poly(A) tail.
|
Most of the surveillance studies of astroviruses focused on humans and domesticated animals, and relatively little is known about the prevalence of astroviruses in wildlife. The role of bats as the reservoirs for zoonotic diseases including rabies, Hendra, Nipah, and Ebola viruses has been highlighted in recent years (reviewed in reference2). Insectivorous bats have also been shown to harbor a range of novel coronaviruses including the precursor of severe acute respiratory syndrome coronavirus (15, 25). Some species of bats live in close proximity to human habitation, and thus it is important to have a better understanding of the virus ecology found in bats. The range and diversity of coronaviruses found in bats have led to the hypothesis that bat coronaviruses may be the precursors of most other mammalian group 1 and 2 coronaviruses (30). These findings highlight the importance of identifying novel viruses in wildlife in general and bats in particular. We used random primers to detect novel viruses in bat fecal specimens. Here, we report the discovery of novel astroviruses in bats in Hong Kong. The remarkably high prevalence and genetic diversity of astroviruses in various bat species found within a relatively small geographic area highlight the need for study in other species of bats and in other geographic locations.
|
|
|---|
Viral nucleic acid extraction and reverse transcription-PCR (RT-PCR). RNA from 140 µl of sample in transport medium was extracted by a QIAamp virus RNA mini kit (Qiagen) following the protocol provided by the manufacturer. Purified RNA was eluted in 60 µl of elution buffer provided in the extraction kit. cDNA was generated from RNA using Superscript III reverse transcriptase (Invitrogen) in a 20-µl reaction mixture containing 150 ng of random hexamers or a 0.5 µM concentration of a gene-specific reverse primer (5'-TTTGGTCCNCCNCTCCAAA-3') targeting the 3' end of ORF1b, 10 mM dithiothreitol, 0.5 mM deoxynucleoside triphosphate mix, 1x First-Strand buffer (Invitrogen), and 200 U of reverse transcriptase. Reaction mixtures were incubated at 25°C for 5 min, followed by 50°C for 60 min, and then the enzyme was inactivated by heating at 70°C for 15 min.
Random hexamer-generated cDNA was screened for the presence of astrovirus using heminested PCR targeting the RdRp gene. A 50-µl PCR mixture was set up containing 1 U of Accuprime Taq DNA polymerase in 1x reaction buffer (Invitrogen), a 2 µM concentration (each) of forward and reverse primers, and 2 µl of cDNA or 1 µl of the first PCR product as a template. First-round PCR was carried out with a mixture of two forward primers, 5'-GARTTYGATTGGRCKCGKTAYGA-3' and 5'-GARTTYGATTGGRCKAGGTAYGA-3', and reverse primer 5'-GGYTTKACCCACATNCCRAA-3'. After an initial incubation at 94°C for 1 min, 30 cycles of amplification were carried out consisting of denaturation at 94°C for 30 s, annealing at 50°C for 30 s, and extension at 68°C for 30 s. Heminested PCR was carried out with a mixture of two forward primers, 5'-CGKTAYGATGGKACKATHCC-3' and 5'-AGGTAYGATGGKACKATHCC-3', and the same reverse primers used in the first-round PCR; the thermocycling conditions were the same as those used for the first-round PCR, except that 40 cycles of amplification were performed. Water controls were included in each run of the RT-PCR assay. PCR products were analyzed by standard agarose gel electrophoresis. The expected product size of the second PCR was 422 bp. All positive results were verified by direct DNA sequencing of the PCR amplicons.
Cloning and sequencing of PCR products. PCR products were purified by a QIAquick PCR purification kit (Qiagen) following the manufacturer's instructions. Long PCR products (product sizes of >1,000 bp) were gel purified with a QIAquick gel extraction kit (Qiagen) and then cloned into pCR2.1-TOPO plasmids (Invitrogen) for DNA sequencing. Multiple clones of a PCR product were picked and sequenced by using a BigDye Terminator, version 3.1, Cycle Sequencing kit (Applied Biosystems). Sequencing products were analyzed by a PRISM 3700 DNA analyzer (Applied Biosystems).
One astrovirus-positive bat specimen, bat astrovirus AFCD337, was chosen for sequencing of the viral genome from PCR products derived from random or sequence-specific primed cDNA. The ORF2 region and the 3' end of the virus sample were amplified using a 3' rapid amplification of cDNA ends system (Invitrogen) following the protocol provided by the manufacturer (primers and PCR conditions available on request). ORF1a (partial) and ORF1b sequences were assembled from multiple overlapping sequences derived from PCR amplicons. Nine additional ORF1a (partial)/ORF1b astrovirus sequences from other samples were obtained using a similar method. The deduced sequences of these samples have at least threefold sequence coverage.
Phylogenetic analysis. Sequence editing and sequence identity calculations were done using BioEdit, version 7.0.4 (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Alignments of nucleotide sequences and amino acid sequences were done using Clustal W (29) with default parameters. Phylogenetic trees were constructed using Clustal X, version 2.0 (14), and Mega 4 (28) by the neighbor-joining method with the nucleotide substitution model of maximum composite likelihood and default parameters. Bootstrap values of the phylogenetic tree constructed were generated by doing 1,000 replicates.
Nucleotide sequence accession numbers. The sequence of bat astrovirus AFCD337 reported in this paper was deposited in the GenBank database under accession number EU847155 and the RdRp and the partial ORF1b sequences of other strains of bat astroviruses were deposited under accession numbers EU847144 through EU847154 and EU847156 through EU847220. For the genetic analysis in this paper, other astrovirus genomes were retrieved from the GenBank, including human astrovirus type 1 strain Oxford (L23513), human astrovirus type 1 strain Dresden (AY720892), human astrovirus type 4 strain Dresden (AY720891), human astrovirus type 4 isolate Goiania/GO/12/95/Brazil (DQ070852), human astrovirus type 4 isolate Guangzhou (DQ344027), human astrovirus type 5 isolate Goiania/GO/12/94/Brazil (DQ028633), human astrovirus type 8 (AF260508), mink astrovirus (NC_004579), ovine astrovirus (NC_002469) and turkey astrovirus type 1 (NC_002470).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Detection of astrovirus in bats by RT-PCR
|
These same specimens had previously been tested for bat coronaviruses (3, 4, 25). While 6% of bats were coinfected with both a bat astrovirus and a bat coronavirus, such coinfection appears to be randomly distributed, and there was no positive or negative statistical association between the presence of these two viruses (chi-square test with Yates correction, P = 0.82).
Longitudinal study of astrovirus in bats at a single habitat. The bats listed in Table 1 include 157 M. magnater and M. pusillus bats that were captured at four separate visits over a 2-year period at one habitat, an abandoned mine cave in Hong Kong. M. magnater bats were found throughout the period while M. pusillus bats were mainly found in two visits carried out in December 2005 and March 2006; only one was found in a visit in August 2005. Sixty-two (54%) out of 115 rectal swabs and 2 (2%) out of 116 throat swabs collected from M. magnater bats and 18 (55%) out of 33 rectal swabs and 5 (15%) out of 33 throat swabs collected from M. pusillus bats were positive for astrovirus. The overall positive rate in individual bats (either rectal or throat swab or both positive) for M. magnater ranged from 36 to 100% at each of the four visits which spanned the winter, spring, and summer seasons; for M. pusillus the positive rates were 50% and 70% in the two instances when adequate bats were sampled.
Genetic and phylogenetic analysis of a novel astrovirus from M. pusillus. The bat astrovirus AFCD337 detected in a rectal specimen from a M. pusillus bat collected in March 2006 was chosen as a representative virus for more extensive genome sequencing. Approximately 74% of the genome of this novel astrovirus was obtained by direct RT-PCR amplification from a rectal swab sample. The partial [excluding the poly(A) tail] 5,067-nucleotide (nt) genome was constructed by aligning sequences from multiple overlapping regions. Amino acid sequences deduced from the viral genome include part of ORF1a and the complete ORF1b, ORF2, and 3' untranslated region (UTR) followed by the poly(A) tail at the 3' end (Fig. 1).
Amino acid sequences encoded by ORFs of the novel bat astrovirus were compared with sequences of other astrovirus genomes including human astrovirus types 1, 4, 5, and 8 (17, 26) as well as mink astrovirus (22), ovine astrovirus, and turkey astrovirus type 1 (11) (Table 2). The findings show that the identified bat astrovirus AFCD337 is a novel mamastrovirus clearly distinct from other known astroviruses. It has <53% and <27% genetic similarity to other known astroviruses in the ORF1b and ORF2 regions, respectively.
|
View this table: [in a new window] |
TABLE 2. Amino acid sequence similarities between prototype bat astrovirus AFCD337 and subgroup A (AFCD68) and B (WCF140) bat astroviruses and other astroviruses
|
The putative ORF2 of the virus has a size of 2,553 nt, which is the largest astrovirus capsid gene known. The N-terminal half of the ORF2 protein, which was previously shown to be more conserved among astroviruses and proposed to be the core assembly domain of the viral capsid (12), was also found to be relatively conserved in this bat astrovirus. The amino acid sequence similarities of this N-terminal half of the bat astrovirus capsid protein to human astrovirus type 1 Oxford, ovine astrovirus, and mink astrovirus are 36.3%, 45.0%, and 39.5%, respectively, compared with <27% similarity for the ORF2 region overall (Table 2). Thus, the C-terminal half of this bat astrovirus protein was highly divergent compared with other astroviruses. This observation supports the speculation that the C-terminal half of the protein is located on the surface of the viral particle and constitutes a region of the capsid that contributes to the species-specific tropism of the virus (12).
Phylogenetic analysis of astrovirus in bats. Phylogenetic analyses of the ORF1a (partial) and ORF1b regions and of ORF2 confirm that bat astrovirus AFCD337 is a novel, distinct astrovirus (Fig. 2). To better define the genetic diversity within the bat astroviruses by species, time, and geographic location, 77 PCR RdRp amplicons (422 nt) obtained from the rectal swab samples in the screening PCR assay were selected for genetic sequence analysis. Bat astrovirus gene sequences were aligned with sequences of other astroviruses including mink astrovirus, ovine astrovirus, human astrovirus type 1 Oxford, human astrovirus type 1 Dresden, human astrovirus type 2, human astrovirus type 4 Dresden, human astrovirus type 4 Goiania, human astrovirus type 4 Guangzhou, human astrovirus type 5 Goiania, and human astrovirus type 8. An avian astrovirus, turkey astrovirus, was included as an outgroup. A phylogenetic tree was constructed from the sequence alignment (Fig. 3). The 72 astroviruses detected in bats cluster together to form a novel group of viruses within the cluster of mamastroviruses (Fig. 3). Within this group are found two subgroups of viruses, one that includes the majority of astroviruses detected from M. magnater, M. pusillus, and M. schreibersii (including bat astrovirus AFCD337) bats and another subgroup that includes most of the viruses detected in M. chinensis and M. ricketti bats. Other than these two major groups of bat astroviruses, a few astroviruses detected in M. magnater and M. pusillus bats and a virus detected in a P. abramus bat appear to have an outgroup relationship to the others, albeit with weak levels of statistical confidence in this phylogenetic topology (Fig. 3). To further investigate this, a 750-nt region of the RdRp and ORF1b (3' end) of these and other representative bat astroviruses were sequenced. A phylogenetic tree based on the aligned protein-encoding sequences was produced (Fig. 4). The phylogeny of the major group of viruses related to bat astrovirus AFCD337 is confirmed. Interestingly, however, AFCD11 from P. abramus and AFCD57 from M. magnater appear phylogenetically related to the human astroviruses although with modest bootstrap support. Their close relationship has also been confirmed by phylogenetic analysis using an alternative method, MrBayes (data not shown).
![]() View larger version (10K): [in a new window] |
FIG. 2. Phylogenetic analysis of the partial ORF1a ( 800 nt), ORF1b, and ORF2 nucleotide sequences comparing bat astrovirus AFCD337 (underlined) with astroviruses of other species. Nine other bat astrovirus sequences are included in the ORF1a (partial) and ORF1b phylogenetic trees. Alignment was based on the encoded amino acid sequences. Astroviruses are indicated as follows: BatAstV, bat; MAstV, mink; OAstV, ovine; HAstV-1, human type 1 (also types 4, 5, and 8); and TAstV-1, turkey type 1.
|
![]() View larger version (12K): [in a new window] |
FIG. 3. Phylogenetic tree constructed with RdRp gene sequences (422 nt) amplified by an RT-PCR screening assay. Sequences of 77 bat astroviruses (BatAstVs) and other respective sequences of different astroviruses isolated from human (HAstV types 1, 2, 4, 5, and 8), mink (MAstV), ovine (OAstV), and turkey type 1 (TAstV-1) viruses were included and aligned based on the nucleotide sequences. For the bat specimens collected during the phase 2 longitudinal study, the sampling dates are indicated in parentheses on the right according to the following code: 1, June 2005; 2, August 2005; 3, December 2005; 4, March 2006.
|
![]() View larger version (18K): [in a new window] |
FIG. 4. Phylogenetic tree constructed with 750-nt sequences of the RdRp gene and ORF1b (3' end) of representative astroviruses isolated from bat (BatAstV), human (HAstV types 1, 4, 5, and 8), mink (MAstV), ovine (OAstV), and turkey type 1 (TAstV-1) viruses. These sequences were aligned based on the encoded amino acid sequences and were reverse-translated back to nucleotides for the phylogenetic analysis.
|
|
|
|---|
The diversity of astroviruses in bats is remarkable. There is no significant phylogenetic clustering of viruses found within a single sampling occasion. The values of astrovirus RdRp amino acid pairwise similarity found within a single bat species, i.e., M. magnater, captured in a single habitat (mine cave) ranged between 51.1% and 100%, with 97.9% of these values lower than 90%. In contrast, the pairwise amino acid sequence similarities between the same gene region of human astroviruses from geographically diverse regions were estimated to range between 92.6% and 99.1%. It has been previously reported that the amino acid sequence identities of the RdRp gene (covering 80% of the RdRp gene regions analyzed in this report) between four groups of avastroviruses, i.e., turkey astrovirus type 1-like viruses, turkey astroviruses type 2-like viruses, avian nephritis virus-like viruses, and chicken-origin astroviruses, detected in different regions were also highly diverse, ranging from 50.1% to 73.8% identity (23). The high virus detection rates of bat astrovirus in our surveillance, taken together with the marked genetic diversity of viruses from bats within the same habitat, are reminiscent of our observations with bat coronaviruses (4). These findings may suggest that bats are persistently infected with astroviruses although mark-recapture studies are needed to confirm this contention. Other mammalian astrovirus infections tend to be short-lived in immunocompetent humans or other animals. However, type 3 human astrovirus has been associated with persistent gastroenteritis in immunocompetent children although the same virus serotype was not repeatedly demonstrated over the full period of clinical diarrhea (1).
Five bat astrovirus sequences from M. magnater, M. pusillus, and P. abramus failed to cluster with the subgroups A and B referred to above, and some of these (AFCD11 from P. abramus and AFCD57 from M. magnater) cluster, rather, with the human astroviruses although with weak statistical support (Fig. 4). Whether these bat viruses are related to the precursor of human astrovirus is yet to be further investigated. Further sequence data from these strains may help elucidate this phylogenetic association.
Evidence of recombination between astroviruses and also between coronaviruses is well documented (13, 24, 27). It has been reported that a stem-loop motif in the 3' UTR was found conserved in mamastroviruses, turkey astrovirus type 1, and avian infectious bronchitis virus, which is a group 3 coronavirus (10). Interestingly, this stem-loop motif was not recognized in bat astrovirus AFCD337. This 3' UTR stem-loop motif is also absent in turkey astrovirus type 2. However, a phylogenetic analysis of 3' UTRs did not indicate a close phylogenetic relationship between the two sequences of bat astrovirus AFCD337 (81 nt) and turkey astrovirus type 2 (accession no. NC_005790) (196 nt), which lack the 3' UTR stem-loop structure (data not shown). Recently, a report on the identification of a novel coronavirus from liver tissue of a whale with pulmonary disease and terminal acute liver failure showed that ORF6 of the novel coronavirus possessed significant amino acid similarity to human astrovirus capsid proteins (21). The high rates of infection of bats in the same mine cave habitat with coronaviruses and astroviruses imply frequent coinfection with both viruses. Therefore, we searched for sequences of bat astrovirus AFCD337 similar to bat coronaviruses using the BLAST program with algorithms allowing a word size down to seven bases. However, no sequences with significant similarity were detected between bat astroviruses and coronaviruses cocirculating within the same species within the same habitat.
The discovery of novel diverse astroviruses in bats and the genetic analysis of such viruses are likely to provide new insights into the ecology and evolution of astroviruses and reinforce the role of bats as a reservoir of viruses that sometimes pose a zoonotic threat to human health. More extensive surveillance for astroviruses in bats of different species and in different geographic areas is needed to further address these questions.
We thank K. S. Cheung, C. T. Shek, and C. S. M. Chan of the Department of Agriculture, Fisheries and Conservation, Hong Kong, for facilitating this study.
Published ahead of print on 11 June 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»