Abstract
Bat coronavirus (BtCoV) is assumed to be a progenitor of severe acute respiratory syndrome (SARS)-related coronaviruses. To explore the distribution of BtCoVs in the Philippines, we collected 179 bats and detected viral RNA from intestinal or fecal samples by RT-PCR. The overall prevalence of BtCoVs among bats was 29.6 %. Phylogenetic analysis of the partial RNA-dependent RNA polymerase gene suggested that one of the detected BtCoVs was a novel alphacoronavirus, while the others belonged to the genus Betacoronavirus. Western blotting revealed that 66.5 % of bat sera had antibodies to BtCoV. These surveys suggested the endemic presence of BtCoVs in the Philippines.
Recently, bats have been recognized as a major reservoir of emerging viral infections that cause serious diseases in humans and other mammals [1]. Several viruses that cause severe diseases in humans have been detected in bats, such as lyssaviruses related to rabies virus, Hendra virus, Nipah virus, Ebola virus, Marburg virus and severe acute respiratory syndrome (SARS) coronavirus [1]. Hence, there is a growing focus on bats as zoonotic hosts.
SARS was caused by a previously unknown coronavirus [2]. From 2002 to 2003, the SARS epidemic affected mainly China, Hong Kong and Taiwan, afflicting about 8000 people and killing over 600 [3]. Investigations of bats revealed that horseshoe bats (Rhinolophus) in southern China were natural reservoirs of a coronavirus related to SARS-CoV [4, 5].
This prompted surveys of bat coronaviruses (BtCoVs) worldwide, and novel viruses were detected in North and South America, Africa, Europe and Asia [6–12], revealing a broad diversity of BtCoVs. Molecular phylogenetic analysis suggested that all mammalian coronaviruses including SARS-CoV were likely to have evolved from BtCoVs [13].
We previously reported that two groups of BtCoVs were detected in the Philippines: one belonging to the genus Alphacoronavirus and the other belonging to the genus Betacoronavirus [14]. To explore the diversity of BtCoVs in the Philippines, we carried out more surveys in this study. We investigated 179 intestinal or fecal samples to detect coronaviral RNA and performed phylogenetic analysis on isolated coronavirus sequences. We also detected virus-specific antibodies, indicating how widely distributed BtCoVs are in the Philippines.
A total of 179 bats (8 species) were captured at five different sites in the Philippines during the months of March 2009, August 2010 and May 2011. Each sampling was performed after receiving permission from the Government of the Philippines. For anesthesia, 15 mg/kg tiletamine and zolazepam (Virbac, Carros, France) were injected intraperitoneally into captured bats, and the anesthetized bats were sacrificed by cardiac exsanguination. Liver, lienal, renal, lung, fecal and intestinal specimens were collected from each captured bat and were shared with the coworkers. All of the experiments using animals were conducted under the supervision of, and with approval from the Animal Care and Use Committee of the Faculty of Agriculture, the University of Tokyo.
To determine the presence of viral RNA in bat samples, reverse transcriptase polymerase chain reaction (RT-PCR) was performed with a pair of consensus primers targeted to a highly conserved region of the RNA-dependent RNA polymerase (RdRp) gene [10]. Viral RNA was extracted from feces or intestinal tissue of the bats using a Total RNA Isolation Mini Kit (Agilent Technologies, Santa Clara, CA, USA). Complementary DNA was synthesized using a SuperScript III Kit (Invitrogen, Carlsbad, CA, USA). PCR was performed using GoTaq PCR Master Mix (Promega, Madison, WI, USA) and 0.2 μM primers (primer 1, 5′-GGTTGGGACTATCCTAAGTGTGA-3′; primer 2, 5′-CCATCATCAGATAGAATCATCATA-3′). PCR products were analyzed by agarose gel electrophoresis. Amplicons (440 nucleotides [nt] long) were gel–purified and directly sequenced with an ABI 3130 XL DNA Analyzer (Applied Biosystems, Foster City, CA, USA).
To detect anti-BtCoV antibodies in bat sera, western blot (WB) analysis was performed. Recombinant BtCoV nucleocapsid (N) protein produced by prokaryotic expression was used as an antigen. To amplify the full-length BtCoV N protein–coding region (1389 nt; AB543561), PCR was performed using KOD–Plus (Toyobo, Osaka, Japan) and 0.3 μM primers (forward primer, 5’-GAACAGATTGGAGGTATGTCTGGACGGAATAAGC-3’; reverse primer, 5’-TGCTCGAGTGCGGCCTTAGGATCTCTCATTAGCA-3’). The amplified fragment was cloned into the pE-SUMOstar vector (LifeSensors, Malvern, PA, USA). E. coli strain BL21 (DE3) was transformed with the recombinant plasmid. N protein expression was induced by the addition of 1 mM isopropyl–β–D–thiogalactopyranoside. The purity of proteins was checked by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), Coomassie brilliant blue (CBB) staining and WB analysis.
Recombinant BtCoV N protein was separated by 8 % SDS-PAGE. After blotting, nitrocellulose membranes were cut into strips and incubated for 1 h with bat serum samples (diluted 1:1000) or horseradish peroxidase (HRP)-conjugated anti-6x His tag antibody (Abcam, Cambridge, UK) at a dilution of 1:5000. The strips were rinsed and further incubated for 1 h with rabbit anti-bat IgG antibody (diluted 1:1000) [15] and subsequently incubated for 1 h with HRP-conjugated donkey anti-rabbit IgG (diluted 1:2000; GE Healthcare, Milwaukee, WI, USA). Signals were detected using an Enhanced Chemiluminescence (ECL) Detection Kit (GE Healthcare).
From 2009 to 2011, we sampled a total of 179 bats at five sites in the Philippines. To detect viral RNA, RT-PCR was performed using a pair of primers for a conserved 440-nt region within the coronavirus RdRp gene. Intestinal samples containing feces obtained in 2009 and feces collected in 2010 and 2011 were examined. As shown in Table 1, the overall prevalence of BtCoV was 53/179 (29.6 %). Amplified fragments were sequenced directly for characterization. A BLAST search (http://blast.ncbi.nlm.nih.gov/Blast) suggested that most of the sequences obtained in this study were similar to the sequence we reported previously (BtCoV/Philippines/Diliman1525G2/2008). These sequences were >98 % similar to one another and were obtained from lesser dog-faced fruit bats (Cynopterus brachyotis) and greater musky fruit bats (Ptenochirus jagori). Sequences from greater musky fruit bats (5/16 specimens positive by RT-PCR) that were similar to BtCoV/Philippines/Diliman1525G2/2008 were also found (87 % identical at the nucleotide level and 95 % identical in their translated amino acid sequences). Sequences that were >98 % identical were also amplified from Philippine sheath-tailed bats (Emballonura alecto), large rufous horseshoe bats (Rhinolophus rufus), Geoffroy’s rousette bats (Rousettus amplexicaudatus), and lesser long-tongued fruit bats (Macroglossus minimus). A novel sequence was detected in a diadem roundleaf bat (Hipposideros diadema), which had 82 % similarity to BtCoV/A977/2005, found in China [16], and 75 % to the BtCoV we reported previously (BtCoV/Philippines/Diliman1552G1/2008) [14]. Phylogenetic reconstruction from isolated sequences was performed using MEGA5 software [17] (Fig. 1). This suggested that one of the detected sequences (from the diadem roundleaf bat) was related to members of the genus Alphacoronavirus (BtCoV2231/Philippines/2010), while most of the BtCoVs belonged to the genus Betacoronavirus (BtCoV/Philippines/Diliman1525G2/2008 and BtCoV2265/Philippines/2010). The two sequences (BtCoV2265/Philippines/2010 and BtCoV2231/Philippines/2010) were deposited into the DNA Data Bank of Japan (DDBJ; http://www.ddbj.nig.ac.jp/) with the accession numbers AB683970 and AB683971, respectively.
WB analysis was performed to detect antibodies to the CoV N protein. The 6x His-SUMO–fused full–length N protein was expressed in BL21 (DE3) cells and used as antigen. The expression of the ~63-kDa recombinant protein was confirmed by CBB staining and WB using an anti-6x His antibody (Fig. 2). As the recombinant N protein was present in the pellet after centrifugation, this pellet was resuspended in phosphate-buffered saline (PBS) containing 0.5 % Triton X-100. It was confirmed that the recombinant protein was pure and showed a single band in WB using an anti-6x His antibody (Fig. 2b). The secondary supernatant (isolated after two rounds of centrifugation of the E. coli cells) was used as antigen to detect antibodies to BtCoV. Bat serum samples that detected a 63-kDa band at a dilution of 1:1000 were considered positive (Fig. 2c). Sera from 10 Leschnault rousette bats (Rousettus leschnaulti), which were obtained from a zoo in Japan, were used as negative controls; none showed bands at 63 kDa. As a positive control, serum from coronavirus infected Leschnault rousette bat was used [14]. Table 1 displays the prevalence of anti-BtCoV antibodies; 119/179 (66.5 %) serum specimens had antibody against N protein. Total WB positivity correlated with total RT-PCR positivity (p < 0.01, Fisher’s exact test).
Our surveys revealed that diverse populations of BtCoVs are continuously circulating in bats in the Philippines. Interestingly, CoVs in two different lineages were detected in greater musky fruit bats (BtCoV/Philippines/Diliman1525G2/2008 and BtCoV2265/Philippines/2010). Because this bat species is endemic to the Philippines [18], BtCoVs may have adapted to the local bat species and developed a unique ecology. While the bat-CoV/Philippines/Diliman1525G2/2008 lineage was only present in lesser dog-faced fruit bats and greater musky fruit bats, the BtCoV2265/Philippines/2010 lineage was detected in various bat species, including both insectivorous and frugivorous bats. These observations imply that interspecies jumping has occurred with BtCoVs in the Philippines. Given the wide host range of mammalian coronaviruses, it would not be surprising if BtCoVs crossed species barriers between frugivorous and insectivorous bats as well as barriers between different frugivorous bats. Besides SARS-CoV, cross-species transmission events have been implicated in the emergence of novel CoVs [19]. Recombination commonly occurs among CoVs [20] and also seems to happen in BtCoVs [21]. These two features could enable novel BtCoVs to emerge, leading to diversity and complexity among BtCoVs.
While numerous surveys to detect BtCoV RNA have been conducted, serological studies to find anti-CoV antibodies in bats are limited. It is possible that this lack of information has made our understanding of BtCoVs obscure. Only a few studies have investigated antibodies to SARS-related coronavirus in bats [4, 22, 23]. Our survey revealed that antibodies to BtCoV N protein were highly prevalent in bat species (66.5 %). This observation also supports the idea of continuous circulation of BtCoVs in the Philippines and implies that repetitive infection enables BtCoVs to circulate among bat colonies.
We previously reported two lineages of BtCoVs in the Philippines. One was an alphacoronavirus and the other was a betacoronavirus. In this study, we found two more distinct lineages: a betacoronavirus related to another previously described virus [14] and a novel alphacoronavirus. Our genetic and serological surveys suggest the widespread distribution and diversity of BtCoVs in the Philippines.
Abbreviations
- BtCoV:
-
Bat coronavirus
- CoV:
-
Coronavirus
- SARS:
-
Severe acute respiratory syndrome
- RT-PCR:
-
Reverse transcriptase polymerase chain reaction
- RdRp:
-
RNA–dependent RNA polymerase
- nt:
-
Nucleotide
- N:
-
Nucleocapsid
- PBS:
-
Phosphate-buffered saline
- WB:
-
Western blotting
- SDS-PAGE:
-
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
- CBB:
-
Coomassie brilliant blue
- HRP:
-
Horseradish peroxidase
References
Cc H, Hk V, Ds R, Tony S, Paul C (2008) Bats prove to be rich reservoirs for emerging viruses. Microbe 3:521–528
Drosten C, Günther S, Preiser W, van der Werf S, Brodt HR, Becker S, Rabenau H, Panning M, Kolesnikova L, Fouchier RA, Berger A, Burguière AM, Cinatl J, Eickmann M, Escriou N, Grywna K, Kramme S, Manuguerra JC, Müller S, Rickerts V, Stürmer M, Vieth S, Klenk HD, Osterhaus AD, Schmitz H, Doerr HW (2003) Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N Engl J Med 348:1967–1976
WHO (2003) Cumulative number of reported probable cases of Severe Acute Respiratory Syndrome (SARS). http://www.who.int/csr/sars/country/2003_05_20/en/. Accessed 27 November 2011
Lau SK, Woo PC, Li KS, Huang Y, Tsoi HW, Wong BH, Wong SS, Leung SY, Chan KH, Yuen KY (2005) Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci USA 102:14040–14045
Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH, Wang H, Crameri G, Hu Z, Zhang H, Zhang J, McEachern J, Field H, Daszak P, Eaton BT, Zhang S, Wang LF (2005) Bats are natural reservoirs of SARS-like coronaviruses. Science 310:676–679
Tong S, Conrardy C, Ruone S, Kuzmin IV, Guo X, Tao Y, Niezgoda M, Haynes L, Agwanda B, Breiman RF, Anderson LJ, Rupprecht CE (2009) Detection of novel SARS-like and other coronaviruses in bats from Kenya. Emerg Infect Dis 15:482–485
Carrington CV, Foster JE, Zhu HC, Zhang JX, Smith GJ, Thompson N, Auguste AJ, Ramkissoon V, Adesiyun AA, Guan Y (2008) Detection and phylogenetic analysis of group 1 coronaviruses in South American bats. Emerg Infect Dis 14:1890–1893
Gloza-Rausch F, Ipsen A, Seebens A, Göttsche M, Panning M, Felix Drexler J, Petersen N, Annan A, Grywna K, Müller M, Pfefferle S, Drosten C (2008) Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. Emerg Infect Dis 14:626–631
Dominguez SR, O’Shea TJ, Oko LM, Holmes KV (2007) Detection of group 1 coronaviruses in bats in North America. Emerg Infect Dis 13:1295–1300
Poon LL, Chu DK, Chan KH, Wong OK, Ellis TM, Leung YH, Lau SK, Woo PC, Suen KY, Yuen KY, Guan Y, Peiris JS (2005) Identification of a novel coronavirus in bats. J Virol 79:2001–2009
Donaldson EF, Haskew AN, Gates JE, Huynh J, Moore CJ, Frieman MB (2010) Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat. J Virol 84:13004–13018
Li L, Victoria JG, Wang C, Jones M, Fellers GM, Kunz TH, Delwart E (2010) Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses. J Virol 84:6955–6965
Vijaykrishna D, Smith GJ, Zhang JX, Peiris JS, Chen H, Guan Y (2007) Evolutionary insights into the ecology of coronaviruses. J Virol 81:4012–4020
Watanabe S, Masangkay JS, Nagata N, Morikawa S, Mizutani T, Fukushi S, Alviola P, Omatsu T, Ueda N, Iha K, Taniguchi S, Fujii H, Tsuda S, Endoh M, Kato K, Tohya Y, Kyuwa S, Yoshikawa Y, Akashi H (2010) Bat coronaviruses and experimental infection of bats, the Philippines. Emerg Infect Dis 16:1217–1223
Omatsu T, Ishii Y, Kyuwa S, Milanda EG, Terao K, Yoshikawa Y (2003) Molecular evolution inferred from immunological cross-reactivity of immunoglobulin G among Chiroptera and closely related species. Exp Anim 52:425–428
Tang XC, Zhang JX, Zhang SY, Wang P, Fan XH, Li LF, Li G, Dong BQ, Liu W, Cheung CL, Xu KM, Song WJ, Vijaykrishna D, Poon LL, Peiris JS, Smith GJ, Chen H, Guan Y (2006) Prevalence and genetic diversity of coronaviruses in bats from China. J Virol 80:7481–7490
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28:2731–2739
Ong P, Rosell-Ambal G, Tabaranza B, Heaney L, Pedregosa M, Paguntalan LM, Cariño AB, Ramayla S, Duya P, Warguez D, Alcala E, Garcia H, Pamaong R, Gonzalez JC, Lorica RP (2008) IUCN RedList of Thresatened Species. http://www.iucn.org. Accessed 1 February 2012
Vijgen L, Keyaerts E, Moës E, Thoelen I, Wollants E, Lemey P, Vandamme AM, Van Ranst M (2005) Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event. J Virol 79:1595–1604
Masters PS (2006) The molecular biology of coronaviruses. Adv Virus Res 66:193–292
Lau SK, Li KS, Huang Y, Shek CT, Tse H, Wang M, Choi GK, Xu H, Lam CS, Guo R, Chan KH, Zheng BJ, Woo PC, Yuen KY (2010) Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. J Virol 84:2808–2819
Müller MA, Paweska JT, Leman PA, Drosten C, Grywna K, Kemp A, Braack L, Sonnenberg K, Niedrig M, Swanepoel R (2007) Coronavirus antibodies in African bat species. Emerg Infect Dis 13:1367–1370
Lau SKP, Poon RWS, Wong BHL, Wang M, Huang Y, Xu H, Guo R, Li KSM, Gao K, Chan K-H, Zheng B-J, Woo PCY, Yuen K-Y (2010) Coexistence of different Genotypes in the same bat and sCharacterization of rousettus bat coronavirus HKU9 belonging to a novel betacoronavirus subgroup. J Virol 84:11385–11394
Acknowledgments
We greatly thank Mr. Edison Cosico and Mr. Eduardo Eres for the arrangements for bat collection. This study was supported by a grant from a Ministry of Education, Culture, Sports, Science and Technology, Japan.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Tsuda, S., Watanabe, S., Masangkay, J.S. et al. Genomic and serological detection of bat coronavirus from bats in the Philippines. Arch Virol 157, 2349–2355 (2012). https://doi.org/10.1007/s00705-012-1410-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00705-012-1410-z