Human Betacoronavirus 2c EMC/2012–related Viruses in Bats, Ghana and Europe (original) (raw)

Abstract

We screened fecal specimens of 4,758 bats from Ghana and 272 bats from 4 European countries for betacoronaviruses. Viruses related to the novel human betacoronavirus EMC/2012 were detected in 46 (24.9%) of 185 Nycteris bats and 40 (14.7%) of 272 Pipistrellus bats. Their genetic relatedness indicated EMC/2012 originated from bats.

Keywords: Africa, Europe, Ghana, coronavirus, bats, human betacoronavirus, CoV-EMC, viruses

Coronaviruses (CoVs) are enveloped viruses with a positive-sense, single-stranded RNA genome (1). CoVs are classified into 4 genera: Alphacoronavirus, Betacoronavirus (grouped further into clades 2a–2d), Gammacoronavirus, and Deltacoronavirus. Two human coronaviruses (hCoVs), termed hCoV-OC43 and −229E, have been known since the 1960s and cause chiefly mild respiratory disease (2). In 2002–2003, an outbreak of severe acute respiratory syndrome (SARS) leading to ≈850 deaths was caused by a novel group 2b betacoronavirus, SARS-CoV (3). A likely animal reservoir for SARS-CoV was identified in rhinolophid bats (4,5). In the aftermath of the SARS pandemic, 2 hCoVs, termed hCoV-NL63 and -HKU1, and numerous novel bat CoVs were described.

In September 2012, health authorities worldwide were notified of 2 cases of severe respiratory disease caused by a novel hCoV (6,7). This virus, termed EMC/2012, was related to the 2c betacoronavirus clade, which had only been known to contain Tylonycteris bat coronavirus HKU4 and Pipistrellus bat coronavirus HKU5 (8).

We previously identified highly diversified alphacoronaviruses and betacoronaviruses, but not clade 2c betacoronaviruses, in bats from Ghana (9). We also identified sequence fragments from a 2c betacoronavirus from 1 Pipistrellus bat in Europe (10). In this study, we analyzed an extended sample of 4,758 bats from Ghana and 272 bats from 4 European countries.

The Study

Fecal specimens were collected from 10 bat species in Ghana and 4 Pipistrellus species in Europe (Table 1). Bats were caught during 2009–2011 with mist nets, as described (9), in 7 locations across Ghana and 5 areas in Germany, the Netherlands, Romania, and Ukraine (Figure 1). The species, age, sex, reproductive status, and morphologic measurements of the bats were recorded. Fecal pellets were collected and suspended in RNAlater Stabilization Reagent (QIAGEN, Hilden, Germany). RNA was purified as described (11). CoV was detected by using nested reverse transcription PCR (RT-PCR) targeting the RNA-dependent RNA polymerase (RdRp) gene (12) (see Table 1 for assay oligonucleotides).

Table 1. Overview of bats tested for 2c betacoronaviruses, Ghana and Europe.

Figure 1.

Location of bat sampling sites in Ghana and Europe. The 7 sites in Ghana (A) and the 5 areas in Europe (B) are marked with dots and numbered from west to east. a, Bouyem (N7°43′24.899′′ W1°59′16.501′′); b, Forikrom (N7°35′23.1′′ W1°52′30.299′′); c, Bobiri (N6°41′13.56′′ W1°20′38.94′′); d, Kwamang (N6°58′0.001′′ W1°16′0.001′′); e, Shai Hills (N5°55′44.4′′ E0°4′30′′); f, Akpafu Todzi (N7°15′43.099′′ E0°29′29.501′′); g, Likpe Todome (N7°9′50.198′′ E0°36′28.501′′); h, Province Gelderland, NED (N52°1′46.859′′ E6°13′4.908′′); i, Eifel area, federal state Rhineland-Palatinate, GER (N50°20′5.316′′ E7°14′30.912′′); j, Holstein area, federal state Schleswig-Holstein, GER (N54°14′51.271′′ E10°4′3.347′′); k, Tulcea county, ROU (N45°12′0.00′′ E29°0′0.00′′); l, Kiev region, UKR (N50°27′0.324′′ E30°31′24.24′′). NED, the Netherlands; GER, Germany; ROU, Romania; UKR, Ukraine.

A novel CoV was detected in insectivorous Nycteris cf. gambiensis specimens (online Technical Appendix Figure, wwwnc.cdc.gov/EID/pdfs/12-1503-Techapp.pdf; GenBank accession nos. JX899382–JX899384). A real-time RT-PCR was designed to permit sensitive and quantitative detection of this CoV (Table 1). Only Nycteris bats were positive for CoV (46 [24.9%] of 185 specimens) (Table 1). Demographic factors predictive of CoV in captured Nycteris bats were assessed. Juvenile bats and lactating females were significantly more likely to be CoV-infected than were adult and nonlactating female bats, respectively (Table 2). Virus concentrations in feces from Nycteris bats were high (median 412,951 RNA copies/g range 323–150,000,000 copies/g).

Table 2. Possible factors predictive of 2c betacoronavirus detection in Nycteris cf. gambiensis bats, Ghana and Europe*.

The 398-bp CoV RdRp screening fragment was extended to 816 bp, as described (5), to enable more reliable taxonomic classification. We previously established _RdRp_-grouping units (RGU) as a taxonomic surrogate to enable prediction of CoV species on the basis of this 816-bp fragment when no full genome sequences could be obtained. According to our classification, the amino acid sequences in the translated 816-bp fragment of the tentative betacoronavirus species (RGU) differed from each other by at least 6.3% (5). The new Nycteris bat CoV differed from the 2c-prototype viruses HKU4 and HKU5 by 8.8%–9.6% and from EMC/2012 by 7.5% and thus constituted a novel RGU. A partial RdRp sequence fragment of a P. pipistrellus bat CoV from the Netherlands, termed VM314 (described by us in [10]), was completed toward the 816-bp fragment to refine the RGU classification of EMC/2012. EMC/2012 differed from VM314 by only 1.8%.

Because of the genetic similarity between EMC/2012 and VM314, we specifically investigated Pipistrellus bats from 4 European countries for 2c betacoronaviruses. We detected betacoronaviruses in 40 (14.7%) of 272 P. pipistrellus, P. nathusii, and P. pygmaeus bats from the Netherlands, Romania, and Ukraine (Table 1; GenBank accession nos. KC243390-KC243392) that were closely related to VM314. The VM314-associated Pipistrellus bat betacoronaviruses differed from EMC/2012 by 1.8%. The difference between EMC/2012 and HKU5 was 5.5%–5.9%. In summary, HKU5, EMC/2012, and the VM314-associated clade form 1 RGU according to our classification system, and the VM314-Pipistrellus bat clade contains the closest relatives of EMC/2012. HKU4 and the Nycteris CoV define 2 separate tentative species in close equidistant relationship.

We conducted a Bayesian phylogenetic analysis. In this analysis, the Nycteris bat CoV clustered as a phylogenetically basal sister clade with HKU4, HKU5, and EMC/2012 and the associated European Pipistrellus viruses (Figure 2).

Figure 2.

RNA-dependent RNA polymerase (RdRp) gene and Spike genephylogenies including the novel betacoronaviruses from bats in Ghana and Europe. A) Bayesian phylogeny of an 816-nt RdRp gene sequence fragment corresponding to positions 14781–15596 in severe acute respiratory syndrome coronavirus (SARS-CoV) strain Frankfurt 1 (GenBank accession no. AY291315). Data were analyzed with MrBayes version 3.1 (http://mrbayes.sourceforge.net/) by using a WAG amino acid substitution model and 4 million generations sampled every 100 steps. Trees were annotated by using a burn-in of 10,000 and visualized with FigTree version 1.6.1 from the BEAST package (www.beast.bio.ed.ac.uk). A whale gammacoronavirus was used as an outgroup. The novel Nycteris bat viruses are shown in boldface and red, the novel Pipistrellus bat viruses and other bat CoVs in the 2c clade are shown in boldface and cyan, and the novel human betacoronavirus EMC/2012 is shown in boldface. Values at deep nodes represent statistical support of grouping by posterior probabilities. CoV clades are depicted to the right of taxa. B) Phylogeny of the complete Spike gene of clade 2c CoVs determined by using the neighbor-joining method with an amino acid percentage distance substitution model and the complete deletion option in MEGA5 (www.megasoftware.net). The Nycteris CoV Spike gene was equidistant from other 2c-CoV Spike genes with 45.6%–46.8% aa divergence. Human coronavirus (hCoV)–OC43 was used as an outgroup. No complete Spike gene sequence was available for VM314 or the novel Pipistrellus bat CoVs. Scale bar represents percentage amino acid distance. The analysis comprised 1,731 aa residues. C) Phylogeny of the partial Spike gene of clade 2c CoVs, including the novel CoVs of Pipistrellus bats from Europe, determined by using a nucleotide distance substitution model and the complete deletion option in MEGA5. Scale bar represents percentage nucleotide distance. The analysis comprised 131 nt corresponding to positions 25378–25517 in hCoV-EMC/2012. Oligonucleotide sequences of primers used to amplify full and partial Spike gene sequences are available on request from the authors. Values at deep nodes in B and C represent statistical support of grouping by percentage of 1,000 bootstrap replicates. GenBank accession numbers for the complete and partial Spike genes correspond to those given in panel A for the RdRp gene.

To confirm the _RdRp_-based classification, we amplified the complete glycoprotein-encoding Spike gene and sequenced it for the novel Nycteris bat virus. The phylogenetically basal position of the novel Nycteris bat virus within the 2c clade resembled that in the CoV RdRp gene (Figure 2). Partial sequences that could be obtained from the 3′-end of the Spike gene of three 2c Pipistrellus bat betacoronaviruses confirmed their relatedness to EMC/2012 (Figure 2).

Conclusions

We detected novel clade 2c betacoronaviruses in Nycteris bats in Ghana and Pipistrellus bats in Europe that are phylogenetically related to the novel hCoV EMC/2012. All previously known 2c bat CoVs originated from vespertilionid bats: VM314 originated from a P. pipistrellus bat from the Netherlands and HKU4 and HKU5 originated from Tylonycteris pachypus and P. abramus bats, respectively, from the People’s Republic of China. The Nycteris bat virus in Africa extends this bat CoV clade over 2 different host families, Nycteridae and Vespertilionidae (Technical Appendix Figure). Detection of genetically related betacoronaviruses in bats from Africa and Eurasia parallels detection of SARS-CoV in rhinolophid bats from Eurasia and related betacoronaviruses in hipposiderid bats from Africa (9).

The relatedness of EMC/2012 to CoVs hosted by Pipistrellus bats at high prevalence across different European countries and the occurrence of HKU5 in bats of this genus from China highlight the possibility that Pipistrellus bats might indeed host close relatives of EMC/2012. This suspicion is supported by observations that tentative bat CoV species (RGUs) are commonly detected within 1 host genus (5). Within the Arabian Peninsula, the International Union for Conservation of Nature (www.iucn.org) lists 50 bat species, including P. arabicus, P. ariel, P. kuhlii, P. pipistrellus, P. rueppellii, and P. savii bats. Because of the epidemiologic link of EMC/2012 with the Arabian Peninsula (6,7), bats from this area should be specifically screened.

The genomic data suggest that EMC/2012, like hCoV-229E and SARS-CoV, might be another human CoV for which an animal reservoir of closely related viruses could exist in Old World insectivorous bats (4,9). Whether cross-order (e.g., chiropteran, carnivore, primate) host switches, such as suspected for SARS-CoV, have occurred for 2c clade bat CoVs remains unknown. However, we showed previously that CoVs are massively amplified in bat maternity colonies in temperate climates (13). This amplification also might apply to the Nycteris bat CoV because, as shown previously for vespertilionid bats from temperate climates (14), detection rates of CoV are significantly higher among juvenile and lactating Nycteris bats. In light of the observed high virus concentrations, the use of water from bat caves and bat guano as fertilizer for farming and the hunting of bats as wild game throughout Africa (15) may facilitate host switching events. To our knowledge, no CoV has been isolated directly from bats. Further studies should still include isolation attempts to obtain full virus genomes and to identify virulence factors that may contribute to the high pathogenicity of EMC/2012 (7).

Technical Appendix Figure

Bat evolutionary lineages and species in which novel group 2c betacoronaviruses were detected, Ghana and Europe.

Acknowledgments

We thank Sebastian Brünink, Tobias Bleicker, and Monika Eschbach-Bludau for technical assistance. We are grateful to Ioan Coroiu, Carsten Dense, Regina Klüppel-Hellmann, Anda Culisier, Danny Culisier, Sabrina Stölting, the volunteers at the Bonn Consortium for Bat Conservation, Andreas Kiefer, Manfred Braun, Isaac Mawusi Adanyeguh, Lucinda Kirkpatrick, Mac Elikem Nutsuakor, David Ofori Agyei, Sarah Koschnicke, Julia Morrison, Emmanual Asare, and Thomas Kruppa for their help during the organization and conduct of field work. We thank Anna Marie Corman for assistance with geographic information processing.

For all capturing, sampling, and exportation of bat specimens, we obtained permission from the respective countries’ authorities.

This study was supported by the European Union FP7 projects EMPERIE (contract number 223498) and ANTIGONE (contract number 278976) and by the German Research Foundation (DFG grant DR 772/3-1, KA1241/18-1).

Biography

Dr Annan is a scientist affiliated with the Kumasi Centre for Collaborative Research in Tropical Medicine, Kumasi, Ghana. Her primary research interest is the characterization of human and novel zoonotic viruses.

Footnotes

Suggested citation for this article: Annan A, Baldwin HJ, Corman VM, Klose SM, Owusu M, Nkrumah EE, et al. Human betacoronavirus 2c EMC/2012–related viruses in bats, Ghana and Europe. Emerg Infect Dis [Internet]. 2013 Mar [_date cited_]. http://dx.doi.org/10.3201/eid1903.121503

1

These authors contributed equally to this article.

References

• 1.Woo PC, Lau SK, Huang Y, Yuen KY. Coronavirus diversity, phylogeny and interspecies jumping. Exp Biol Med (Maywood). 2009;234:1117–27. 10.3181/0903-MR-94 [DOI] [PubMed] [Google Scholar]
• 2.Saif LJ. Animal coronaviruses: what can they teach us about the severe acute respiratory syndrome? Rev Sci Tech. 2004;23:643–60 . [DOI] [PubMed] [Google Scholar]
• 3.Drosten C, Gunther S, Preiser W, van der Werf S, Brodt HR, Becker S, et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N Engl J Med. 2003;348:1967–76. 10.1056/NEJMoa030747 [DOI] [PubMed] [Google Scholar]
• 4.Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH, et al. Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005;310:676–9. 10.1126/science.1118391 [DOI] [PubMed] [Google Scholar]
• 5.Drexler JF, Gloza-Rausch F, Glende J, Corman VM, Muth D, Goettsche M, et al. Genomic characterization of severe acute respiratory syndrome–related coronavirus in European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences. J Virol. 2010;84:11336–49. 10.1128/JVI.00650-10 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 6.Corman VMEI, Bleicker T, Zaki A, Landt O, Eschbach-Bludau M, van Boheemen S, et al. Detection of a novel human coronavirus by real-time reverse-transcription polymerase chain reaction. Euro Surveill. 2012;17:20285 . [DOI] [PubMed] [Google Scholar]
• 7.Zaki AM, van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N Engl J Med. 2012;367:1814–20. 10.1056/NEJMoa1211721 [DOI] [PubMed] [Google Scholar]
• 8.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 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 9.Pfefferle S, Oppong S, Drexler JF, Gloza-Rausch F, Ipsen A, Seebens A, et al. Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana. Emerg Infect Dis. 2009;15:1377–84. 10.3201/eid1509.090224 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 10.Reusken CB, Lina PH, Pielaat A, de Vries A, Dam-Deisz C, Adema J, et al. Circulation of group 2 coronaviruses in a bat species common to urban areas in western Europe. Vector Borne Zoonotic Dis. 2010;10:785–91. 10.1089/vbz.2009.0173 [DOI] [PubMed] [Google Scholar]
• 11.Drexler JF, Corman VM, Muller MA, Maganga GD, Vallo P, Binger T, et al. Bats host major mammalian paramyxoviruses. Nat Commun. 2012;3:796. 10.1038/ncomms1796 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 12.de Souza Luna LK, Heiser V, Regamey N, Panning M, Drexler JF, Mulangu S, et al. Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription–PCR and nonfluorescent low-density microarray. J Clin Microbiol. 2007;45:1049–52. 10.1128/JCM.02426-06 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 13.Drexler JF, Corman VM, Wegner T, Tateno AF, Zerbinati RM, Gloza-Rausch F, et al. Amplification of emerging viruses in a bat colony. Emerg Infect Dis. 2011;17:449–56. 10.3201/eid1703.100526 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 14.Gloza-Rausch F, Ipsen A, Seebens A, Gottsche 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 [DOI] [PMC free article] [PubMed] [Google Scholar]
• 15.Mickleburgh S, Waylen K, Racey P. Bats as bushmeat: a global review. Oryx. 2009;43:217–34. 10.1017/S0030605308000938 [DOI] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Technical Appendix Figure

Bat evolutionary lineages and species in which novel group 2c betacoronaviruses were detected, Ghana and Europe.