Detection and Phylogenetic Analysis of Group 1 Coronaviruses in South American Bats (original) (raw)

• Journal List
• Emerg Infect Dis
• v.14(12); 2008 Dec
• PMC2634629

Emerg Infect Dis. 2008 Dec; 14(12): 1890–1893.

University of the West Indies, St. Augustine, Republic of Trinidad and Tobago (C.V.F. Carrington, J.E. Foster, N. Thompson, A.J. Auguste, V. Ramkissoon, A.A. Adesiyun)

University of Hong Kong, Hong Kong Special Administrative Region, People’s Republic of China (H.C. Zhu, J.X. Zhang, G.J.D. Smith, Y. Guan)

Corresponding author.

Address for correspondence: Christine V. F. Carrington, Biochemistry Unit, Bldg 36, Department of Preclinical Sciences, Faculty of Medical Sciences, University of the West Indies, Eric Williams Medical Sciences Complex, Mt. Hope, Republic of Trinidad and Tobago; email: ude.iwu.ats@notgnirrac.enitsirhc

This is a publication of the U.S. Government. This publication is in the public domain and is therefore without copyright. All text from this work may be reprinted freely. Use of these materials should be properly cited.

Abstract

Bat coronaviruses (Bt-CoVs) are thought to be the precursors of severe acute respiratory syndrome coronavirus. We detected Bt-CoVs in 2 bat species from Trinidad. Phylogenetic analysis of the RNA-dependent RNA polymerase gene and helicase confirmed them as group 1 coronaviruses.

Keywords: Bat, coronavirus, Bt-CoV, Trinidad, South America, dispatch

Bats are of particular interest as reservoirs for potentially emergent pathogens. Because of their abundance, wide distribution, and mobility, bats confer a greater risk for zoonotic transmission than other animals (1). Bats have long been known as the natural hosts for rabies virus and other lyssaviruses and were more recently identified as the reservoirs for emerging viruses such as Ebola, Hendra, and Nipah virus (reviewed in 1). The search for the animal reservoir of the severe acute respiratory syndrome coronavirus (SARS-CoV) led to extensive surveys of coronaviruses in wild and domestic animal populations in China, resulting in the detection of a wide variety of novel bat coronaviruses (Bt-CoVs) (2–5). The data suggest that the progenitor of the SARS-CoV, and all other coronaviruses in other animal hosts, originated in bats (4). Recent reports by Dominguez et al. (6) and Gloza-Rausch et al_._ (7) confirmed the existence of Bt-CoVs outside China, in the United States and Germany, respectively. Additionally antibodies reactive with SARS-CoV have been detected in African bat species (8). We report the detection and characterization of CoVs in bats from Trinidad, the southernmost island of the Caribbean archipelago, located 9 km (5.5 miles) off the northeastern coast of South America.

The Study

A total of 114 bats collected from their natural habitats from December 2006 through July 2007 (8 species from 10 locations; Table) were euthanized after being deeply anesthetized with 2% xylazine and 10% ketamine administered subcutaneously (in some cases after sedation with CO2). The bats were then taxonomically classified on the basis of morphology, and the carcasses were stored at –70°C until used. For sampling, the carcasses were thawed at 4°C for 3–4 hours, then oropharyngeal and anal samples were taken with Dacron-tipped swabs that were then placed in RNA_later_ (Ambion, Austin, TX, USA) and stored at –20°C until used.

Table

Number and location of bat species collected and tested from December2006 through July 2007, with bat coronavirus species status

*Indicates group from which coronavirus positive bats originated.

CoV detection and sequencing were conducted as previously described (4). Briefly, viral RNA was extracted from swabs by using the QIAamp viral RNA minikit (QIAGEN, Westburg, the Netherlands) and used as the template for reverse transcription–PCR (RT-PCR) detection of the CoV RNA–dependent RNA polymerase (RdRp) gene (9). Primers based on the RdRp gene, conserved for all known coronaviruses, were then used for RT-PCR detection. The RdRp PCR products were gel purified by using the QIAquick PCR purification kit (QIAGEN) and sequenced to confirm virus species. RNA from samples positive for coronavirus was then used for cDNA synthesis by using random hexamer, gene-specific, and oligo(dT) primers. The RdRp gene and 1b open reading frames, including the helicase (HEL) domain, were then sequenced, also as previously described (4). Sequences derived from this study were deposited in GenBank (accession nos. EU769557 and EU769558).

Sequences were aligned with previously published CoV sequences from GenBank by using ClustalX (http://bips.u-strasbg.fr/fr/Documentation/ClustalX) then manually aligned by using the Se-Al program (http://tree.bio.ed.ac.uk/software/seal). The GenBank accession numbers of all sequences used are noted in the taxon names in Figures 1 and ​2. The RdRp sequences were trimmed to equal length, which created 2 datasets of 780 bp (n = 40) and 378 bp (n = 45). The latter included Bt-CoV sequences from North America and Germany that were too short to be included in the first dataset. A third dataset comprised an alignment of the HEL domain (n = 46) trimmed to 1,797 bp. Maximum likelihood (ML) phylogenies were inferred under a General Time Reversible (GTR + Γ4 + I) model, which was identified as the best-fit model of nucleotide substitution using MODELTEST version 3.7 (10). Bootstrapping was performed to assess the robustness of tree topologies by using 1,000 replicate neighbor-joining (NJ) trees under the ML substitution model. All analyses were performed with PAUP* version 4.0b (Sinauer Associates, Inc., Sunderland, MA, USA).

Maximum likelihood trees of coronaviruses based on A) 780-bp fragment of the RNA-dependent RNA polymerase gene and B) 1,797 bp of the helicase (HEL) domain of open reading frame 1b. Trees were inferred under the General Time Reversible (GTR + Γ4 + I) model by using PAUP* version 4.0b (Sinauer Associates, Inc., Sunderland, MA, USA). Bootstrap support values >90% are indicated. Previously defined phylogenetic groups and a putative novel group (10) are delineated by the bars on the right. The numbering of these groups is as described in the eighth report of the International Committee on Taxonomy of Viruses with the alternative grouping proposed by Tang et al_._ (4) in brackets. Trinidadian bat coronavirus sequences are highlighted in red. GenBank accession numbers are noted in parentheses. Scale bars indicate number of nucleotide substitutions per site.

Maximum likelihood tree of coronaviruses based on 378-bp fragment of the RNA-dependent RNA polymerase gene. The tree was inferred under the General Time Reversible (GTR + Γ4 + I) by using PAUP* version 4.0b (Sinauer Associates, Inc., Sunderland, MA, USA). Trinidadian bat coronavirus (Bt-CoV) sequences are highlighted in red and North American Bt-CoV in blue. Previously defined phylogenetic groups and a putative novel group (10) are delineated by the bars on the right. The numbering of these groups is as described in the eighth report of the International Committee on Taxonomy of Viruses with the alternative grouping proposed by Tang et al_._ (4) in brackets. Bootstrap support values for groups 1a, 1b, 2a–c, 3, and the lineage containing Trinidadian Bt-CoVs are shown. GenBank accession numbers are noted in parentheses. Scale bar indicates number of nucleotide substitutions per site.

CoV RNA was detected in 1 of 21 Glossophaga soricina and 1 of 5 Carollia perspicillata bats tested**.** The latter, designated Bt-CoV/Trinidad/1FY2B, was an adult male from the area of Fyzabad, and virus was detected in anal and oropharyngeal swabs. The infected G. soricina specimen (Bt-CoV/Trinidad/1CO7B) was an adult female from the area of Couva, and the virus was detected in the anal swab only. Subsequent sequencing was performed on virus from the anal swabs. A total of 3,905 bp of RdRp from Bt-CoV/Tri/1CO7B and 5,160 bp from Bt-CoV/Tri/1FY2B were sequenced, in addition to 1,782 bp from the HEL domain of both samples (full sequence data not shown). Across the shared 3,798-bp region of RdRp, divergence was 28.7% at the nucleotide level and 15.5% at the amino acid level. Similarly, the HEL domains of the 2 viruses showed 27.3% divergence at the nucleotide level and 13.3% at the amino acid level.

ML phylogenies inferred from the RdRp genes of 40 viruses (780 bp) and the HEL domains of 46 viruses (1,797 bp) are shown in Figure 1. In both cases, the 5 groups proposed by Tang et al. (4) on the basis of NJ trees, and the lineage containing the recently reported novel CoV sequences (11), were strongly supported with bootstrap values >95% in all cases. In each case the Trinidadian sequences clustered with group 1 CoVs within a clade containing all other group 1 bat and human CoVs as well as porcine CoV. To determine the phylogenetic relationship between the Trinidadian Bt-CoVs and North American Bt-CoVs, for which only relatively short sequences for RdRp were available, we inferred a second RdRp ML phylogeny based on a 378-bp fragment (Figure 2). When this shorter fragment was used, 5 groups were again defined, but bootstrap support for the putative group 5 (4) was lower. The level of divergence between Trinidadian sequences was notably higher than among the North American (6) and German sequences (7).

Conclusions

Our detection of RNA from group 1 CoVs in Trinidadian bats shows that Bt-CoVs have a wider distribution than previously suspected and is added support for bats as the original host species for these viruses. Group 1 CoVs form 2 well-supported clades designated 1a and 1b (12). The Trinidadian bt-CoV clustered within the latter clade, which contains all other group 1 Bt-CoVs, including those from Germany and North America, and the 3 known group 1 human CoVs associated with respiratory illness (13–15). Despite the geographic proximity of the bats from which the Trinidadian Bt-CoV sequences were derived—Couva and Fyzabad are 28 km (17 miles) apart, and Trinidad has an area of only 4,769 km2 (1,864 square miles)—they are relatively highly divergent. This divergence might reflect virus adaptations to different host species; however, more data would be needed to confirm this. Given the mobility of bats, the possibility of the viruses having different geographic origins (perhaps even from outside Trinidad) cannot be ruled out. Further work on CoV diversity in Trinidad and the rest of the Americas, as well as on the ecology and behavior of susceptible bat species, is needed to understand the origins, evolution, and dispersal of these viruses.

Acknowledgments

We are grateful to Kirk Armour and other staff of the Anti-Rabies Unit in the Ministry of Agriculture, Land and Marine Resources, Trinidad and Tobago, for assistance with field collection of bats; to Vijaykrishna Dhanasekaran for technical assistance with sequencing; and to Rajendra Mahabir for assistance with taxonomic classification of bats.

Funding was provided by the University of the West Indies St. Augustine Campus Research and Publications Fund and the Li Ka Shing Foundation.

Biography

•

Dr Carrington is a senior lecturer in the Department of Preclinical Sciences of the Faculty of Medical Sciences, University of the West Indies, Trinidad and Tobago. Her main research interests are ecologic and evolutionary factors involved in emergence, dispersal, and maintenance of viral pathogens, and in particular zoonotic and vector-borne RNA viruses.

Footnotes

Suggested citation for this article: Carrington CVF, Foster JE, Zhu HC, Zhang JX, Smith GJD, Thompson N, et al. Detection and phylogenetic analysis of group 1 coronaviruses in South American bats. Emerg Infect Dis [serial on the Internet]. 2008 Dec [_date cited_]. Available from http://www.cdc.gov/EID/content/14/12/1890.htm

References

1. Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006;19:531–45. 10.1128/CMR.00017-06 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

2. Lau SK, Woo PC, Li KS, Huang Y, Tsoi HW, Wong BH, et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci U S A. 2005;102:14040–5. 10.1073/pnas.0506735102 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

3. Poon LL, Chu DK, Chan KH, Wong OK, Ellis TM, Leung YH, et al. Identification of a novel coronavirus in bats. J Virol. 2005;79:2001–9. 10.1128/JVI.79.4.2001-2009.2005 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

4. Tang XC, Zhang JX, Zhang SY, Wang P, Fan XH, Li LF, et al. Prevalence and genetic diversity of coronaviruses in bats from China. J Virol. 2006;80:7481–90. 10.1128/JVI.00697-06 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

5. Woo PC, Lau SK, Li KS, Poon RW, Wong BH, Tsoi HW, et al. Molecular diversity of coronaviruses in bats. Virology. 2006;351:180–7. 10.1016/j.virol.2006.02.041 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

6. Dominguez SR, O’Shea TJ, Oko LM, Holmes KV. Detection of group 1 coronaviruses in bats in North America. Emerg Infect Dis. 2007;13:1295–300. [PMC free article] [PubMed] [Google Scholar]

7. Gloza-Rausch F, Ipsen A, Seebens A, Göttsche M, Panning M, Felix Drexler J, et al. Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. Emerg Infect Dis. 2008;14:626–31. 10.3201/eid1404.071439 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

8. Muller MA, Paweska JT, Leman PA, Drosten C, Grywna K, Kemp A, et al. Coronavirus antibodies in African bat species. Emerg Infect Dis. 2007;13:1367–70. [PMC free article] [PubMed] [Google Scholar]

9. Guan Y, Zheng BJ, He YQ, Liu XL, Zhuang ZX, Chueng CL, et al. Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China. Science. 2003;302:276–8. 10.1126/science.1087139 [PubMed] [CrossRef] [Google Scholar]

10. Posada D, Crandall KA. Modeltest: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8. 10.1093/bioinformatics/14.9.817 [PubMed] [CrossRef] [Google Scholar]

11. Dong BQ, Liu W, Fan XH, Vijaykrishna D, Tang XC, Gao F, et al. Detection of a novel and highly divergent coronavirus from Asian leopard cats and Chinese ferret badgers in Southern China. J Virol. 2007;81:6920–6. 10.1128/JVI.00299-07 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

12. Gorbalenya AE, Snijder EJ, Spaan WJ. Severe acute respiratory syndrome coronavirus phylogeny: toward consensus. J Virol. 2004;78:7863–6. 10.1128/JVI.78.15.7863-7866.2004 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

13. Hamre D, Procknow JJ. A new virus isolated from the human respiratory tract. Proc Soc Exp Biol Med. 1966;121:190–3. [PubMed] [Google Scholar]

14. van der Hoek L, Pyrc K, Jebbink MF, Vermeulen-Oost W, Berkhout RJ, Wolthers KC, et al. Identification of a new human coronavirus. Nat Med. 2004;10:368–73. 10.1038/nm1024 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

15. Fouchier RA, Hartwig NG, Bestebroer TM, Niemeyer B, de Jong JC, Simon JH, et al. A previously undescribed coronavirus associated with respiratory disease in humans. Proc Natl Acad Sci U S A. 2004;101:6212–6. 10.1073/pnas.0400762101 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from Emerging Infectious Diseases are provided here courtesy of Centers for Disease Control and Prevention