Volume 19, Number 8—August 2013
Group C Betacoronavirus in Bat Guano Fertilizer, Thailand
To the Editor: Bats play a critical role in the transmission and origin of zoonotic diseases, primarily viral zoonoses associated with high case-fatality rates, including those caused by Nipah virus (NiV) and severe acute respiratory syndrome (SARS)–like coronavirus (CoV) infections (1). Recently, the World Health Organization (WHO) reported 44 confirmed cases of human infection with Middle East respiratory syndrome CoV, resulting in 22 deaths. Full-genome and phylogenetic analyses of these Middle East respiratory syndrome CoVs have been published elsewhere (2). The identified viruses from 2 patients (previously referred to as England/Qatar/2012 and EMC/2012) are genetically related and belong to group C betacoronavirus, which is most related to CoVs from Nycteris bats in Ghana and Pipistrellus bats in Europe (2,3). In addition, bat CoVs HKU4 and HKU5 originated from Tylonycteris pachypus and Pipistrellus abramus bats, respectively, in the People’s Republic of China (4). Bats are also known to harbor and transmit nonviral zoonotic pathogens, including the fungal pathogen Histoplasma capsulatum, which causes histoplasmosis in humans (5).
Bat guano is sold for use as a fertilizer in several countries, including Thailand, Indonesia, Mexico, Cuba, and Jamaica. The practice of collecting and harvesting bat guano may pose a considerable health risk because guano miners have a high level of contact and potential exposure to bat-borne pathogens. To assess pathogens in bat guano, we examined bat guano from a cave in the Khao Chong Phran Non-hunting Area (KCP-NHA) in Ratchaburi Province, Thailand, where bat guano was sold as agricultural fertilizer, for the presence of NiV, CoV, and H. capsulatum fungi. Bats from 14 species in 7 families have been found roosting within this area. Tadarida plicata bats are the most abundant species (2,500,000 bats), and 3 other species of bats found at the site each had thousands of members: Taphozous melanopogon, Taphozous theobaldi, and Hipposideros larvatus.
A random sample of dry bat guano, ≈100 g, was collected in a sterile plastic bag weekly from the main cave at KCP-NHA from September 2006 through August 2007. The specimens were sent for analysis by express mail (at room temperature within 2–3 days) to the WHO Collaborating Centre for Research and Training in the Viral Zoonoses Laboratory at Chulalongkorn University. Samples were frozen immediately at –80°C until nucleic acids were extracted and PCR assays were run. A total of 52 collected bat guano specimens were examined in this study.
Two aliquots of feces from each weekly specimen (104 samples total) were screened for CoV, NiV, and H. capsulatum by PCR. RNA was extracted from 10 mg of fecal pellet by using the QIAamp Viral RNA Mini Kit (QIAGEN, Hilden, Germany). CoV RNA was detected by using nested reverse transcription PCR with the degenerated primers to amplify the RNA-dependent RNA polymerase (RdRp) gene (6). NiV RNA was detected by duplex nested reverse transcription PCR (7). To detect H. capsulatum and other fungi, we extracted genomic DNA directly from bat guano by using the silica-guanidine thiocyanate protocol, NucliSense Isolation Reagent (bioMérieux, Boxtel, the Netherlands), according to the manufacturer’s protocol. We tested for fungal ribosomal DNA (rDNA) in extracted total nucleic acid specimens by using the PCR protocol designed to amplify all rDNA from 4 major fungus phyla at the internal transcribed spacer 1 and 2 regions (8).
Four (3.8%) of 104 samples were positive for CoV. They were collected on September 2, 2006 (KCP9), October 26, 2006 (KCP12), November 14, 2006 (KCP15), and March 4, 2007 (KCP31). Three of the 4 positive CoV sequences (KCP9, KCP12, and KCP15) were identical at 152 nt of the RdRp region (ATCGTGCTATGCCTAATATGTGTAGGATTTTTGCATCTCTCATATTAGCTCGTAAACACAATACTTGTTGTAGTGTTTCAGACCGCTtTtATAGACTTGCaAACGAGTGTGCGCAAGTCTTGAGTGAGTATGTGCTATGTGGTGGTGGCTAT) and phylogenetically clustered with the group C betacoronavirus (Figure), with 76%, 80%, and 77% nt identity to bat CoV HKU4, bat CoV HKU5, and human CoV EMC and England1_CoV, respectively. The other CoV sequence (KCP31: ATCGTGCACTTCCCAATATGATACGCATGATTTCCGCCATGATTTTGGGATCAAAGCATGTTACTTGCTGTGACACATCTGATAAGTATTACCGTCTTTGTAATGAGCTtGCACAAGTTTTGACAGAGGTTGTTTATTCTAATGGTGGTTTC) showed 82% nt identity with bat CoV HKU8, an alphacoronavirus. Although we recognize that longer sequences or full genomes may alter the topology of the phylogeny slightly and give stronger branch support, we expect that the overall topology and placement of these CoVs would remain consistent. Samples from particular bat species could not be identified because bats of different species roost in this cave, and samples were pooled during collection for bat guano fertilizer. The detection of CoVs in bat guano from the KCP-NHA cave in Ratchaburi was consistent with the previous finding of alphacoronavirus from Hipposideros armiger bats from the same province in 2007, but those researchers tested fresh bat feces (9).
All bat guano samples screened by PCR were negative for NiV and Histoplasma spp. but were positive for group C betacoronavirus. The natural reservoir and complete geographic distribution of this CoV are currently unknown. Although we did not isolate live virus from these samples, the detection of nucleic acid and previous isolation of viruses from bat feces and urine (10) warrants some concern that guano miners might be exposed to bat pathogens in fresh excreta as well as in soil substances. We suggest that guano miners use preventive measures of personal hygiene and improved barrier protection to reduce the possibility of exposure to zoonotic pathogens.
This research was supported by the Thailand Research Fund (DBG5180026, RDG5420089); the Higher Education Research Promotion and National Research University Project of Thailand, Office of the Higher Education Commission (HR1160A-55); the Thai Red Cross Society, Thailand; the Naval Health Research Center BAA-10-93 under Cooperative Agreement no. W911NF-11-2-004; the Centre National de la Recherche Scientifique; a nonbiodefense Emerging Infectious Diseases Research Opportunities Award from the National Institute of Allergy and Infectious Diseases of the National Institutes of Health (1 R01 AI079231-01); and the United States Agency for International Development Emerging Pandemic Threats Program, PREDICT project, under the terms of Cooperative Agreement no. GHN-A-OO-09-00010-00.
- Olival KJ, Epstein JH, Wang LF, Field HE, Daszak P. Are bats unique viral reservoirs? In: Aguirre AA, Ostfeld RS, Daszak P, editors. New directions in conservation medicine: applied cases of ecological health. 2nd ed. Oxford: Oxford University Press; 2012. p. 195–212.
- Cotten M, Lam TT, Watson SJ, Palser AL, Petrova V, Grant P, Full-genome deep sequencing and phylogenetic analysis of novel human betacoronavirus. Emerg Infect Dis. 2013;19:736–42.
- Annan A, Baldwin HJ, Corman VM, Klose SM, Owusu M, Nkrumah EE, Human betacoronavirus 2c EMC/2012–related viruses in bats, Ghana and Europe. Emerg Infect Dis. 2013;19:456–9 .
- Woo PC, Lau SK, Li KS, Poon RW, Wong BH, Tsoi HW, Molecular diversity of coronaviruses in bats. Virology. 2006;351:180–7.
- Lyon GM, Bravo AV, Espino A, Lindsley MD, Gutierrez RE, Rodriguez I, Histoplasmosis associated with exploring a bat-inhabited cave in Costa Rica, 1998–1999. Am J Trop Med Hyg. 2004;70:438–42 .
- Tong S, Conrardy C, Ruone S, Kuzmin IV, Guo X, Tao Y, Detection of novel SARS-like and other coronaviruses in bats from Kenya. Emerg Infect Dis. 2009;15:482–5.
- Wacharapluesadee S, Hemachudha T. Duplex nested RT-PCR for detection of Nipah virus RNA from urine specimens of bats. J Virol Methods. 2007;141:97–101.
- White TJ, Bruns T, Lee S, Taylor JW. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods and applications. New York: Academic Press, Inc.; 1990. p. 315–22.
- Gouilh MA, Puechmaille SJ, Gonzalez JP, Teeling E, Kittayapong P, Manuguerra JC. SARS-coronavirus ancestor’s foot-prints in South-East Asian bat colonies and the refuge theory. Infect Genet Evol. 2011;11:1690–702.
- Chua KB, Koh CL, Hooi PS, Wee KF, Khong JH, Chua BH, Isolation of Nipah virus from Malaysian Island flying-foxes. Microbes Infect. 2002;4:145–51.