key: cord-0989681-fqxjzavm
authors: Anindita, Paulina Duhita; Sasaki, Michihito; Setiyono, Agus; Handharyani, Ekowati; Orba, Yasuko; Kobayashi, Shintaro; Rahmadani, Ibnu; Taha, Siswatiana; Adiani, Sri; Subangkit, Mawar; Nakamura, Ichiro; Sawa, Hirofumi; Kimura, Takashi
title: Detection of coronavirus genomes in Moluccan naked-backed fruit bats in Indonesia
date: 2015-02-04
journal: Arch Virol
DOI: 10.1007/s00705-015-2342-1
sha: f8f4ccd8f853366f0593b634dc6f4a49079a4a4b
doc_id: 989681
cord_uid: fqxjzavm

Bats have been shown to serve as natural reservoirs for numerous emerging viruses including severe acute respiratory syndrome coronavirus (SARS-CoV). In the present study, we report the discovery of bat CoV genes in Indonesian Moluccan naked-backed fruit bats (Dobsonia moluccensis). A partial RNA-dependent RNA polymerase gene sequence was detected in feces and tissues samples from the fruit bats, and the region between the RdRp and helicase genes could also be amplified from fecal samples. Phylogenetic analysis suggested that these bat CoVs are related to members of the genus Betacoronavirus.

Deltacoronavirus [10] . Members of the genera Alphacoronavirus and Betacoronavirus are known to cause human disease, whereas those of the genera Gammacoronavirus and Deltacoronavirus are causative agents of animal disease [11] . Bats are known to be reservoirs for many emerging infectious diseases that have zoonotic potential [12] and are believed to be the original source of SARS-CoV [13] and MERS-CoV [14] . In Southeast Asia, bats harboring coronaviruses have been discovered in a number of countries, including the Philippines [15, 16] and Thailand [17] . However, there have been no studies focused on CoVs in Indonesia. Here, we report the discovery of bat CoVs in Moluccan naked-backed fruit bats, Dobsonia moluccensis, in Indonesia.

We University. Fecal, tracheal swabs, and tissues samples including brain, lung, liver, spleen, and kidney from each bat were collected and placed into RNAlater (Life Technologies, Carlsbad, CA). All samples were exported to Japan with the permission of the Directorate General of Livestock and Animal Health Services, the Ministry of Agriculture, Republic of Indonesia.

Morphological features and nucleotide sequence analysis of both mitochondrial 16S rRNA and cytochrome b were conducted as described previously [18] to identify the species of each fruit bat. Feces and tissue samples were subjected to RNA extraction using TRIzol Reagent (Life Technologies) according to the manufacturer's instructions. The RNAs were then screened for genomes of CoVs using a nested RT-PCR with primer sets that specifically amplify the conserved region of RdRp, as described previously [19] (region A in Fig. 1 ). The helicase region was amplified using degenerate primers specific for members of the genus Betacoronavirus [5] (region B in Fig. 1 ). PCR products of both the RdRp and helicase regions were gel-purified, and nucleotide sequencing was performed using a BigDye Terminator v.1.1/3.1 Cycle Sequencing Kit (Life Technologies) and analyzed using an ABI PRISM 3130 Genetic Analyzer (Life Technologies). One-step RT-PCR was performed using a new primer set and a PrimeScript II High Fidelity OneStep RT-PCR Kit (TaKaRa, Tokyo, Japan) to amplify the region not covered in the first and second sets of PCR experiments (region C in Fig. 1 ). Based on alignment of the sequences obtained in this study, the forward primer 5 0 -CACGCAACTTGTTGTAATGCGTCAGAGA-3 0 and the reverse primer 5 0 -CACGTGCTTTTGCAGGCACTAT ACGAC-3 0 , corresponding to the RdRp gene at nucleotide position 14850-14877 and to the helicase gene at nucleotide position 16714-16740, respectively, of bat coronavirus (BatCoV) HKU9 (NC009021) were used. The fragments obtained, covering part of the RdRp and helicase genes (region D in Fig. 1 ), were sequenced.

PCR-amplified nucleotide sequences were aligned and compared to the sequences of known CoVs using the MEGA 6 program. Nucleotide sequence identity values were calculated using GENETYX software ver. 10 (GENETYX, Tokyo, Japan). Phylogenetic trees were constructed based on nucleotide sequences, using the maximum-likelihood method with the Tamura-Nei model and 1000 bootstrap replicates [20, 21] .

Virus isolation was attempted by inoculation of African green monkey kidney (Vero E6), Yaeyama flying fox kidney (FBKT), and Leschenault's rousette kidney (Dem-KT1) cells. Vero E6 cells were obtained from Prof. Ayato Takada (Research Center for Zoonosis Control, Hokkaido University, Hokkaido, Japan) [22] , and FBKT and Dem-KT1 cells were obtained from Prof. Ken Maeda (Yamaguchi University, Yamaguchi, Japan) [23] . Fecal and tissue homogenates that were positive for CoV RNA were used for inoculation of cells after clearance by centrifugation at 10009g for 5 min. Cells were incubated at 37°C for 30 min with fecal samples and for 1 h with tissue homogenates. Inocula were replaced by fresh culture medium supplemented with 2 % FBS, 2 % antibioticantimycotic solution containing penicillin (10,000 U/ml), streptomycin (10 lg/ml) and amphotericin B (25 lg/ml) (Sigma, St. Louis, MO), and gentamycin (25 lg/ml). The supernatants of cultured cells were transferred to fresh cell monolayers each week. RNAs were extracted from culture supernatants and subjected to nested RT-PCR targeting the RdRp gene up until the fifth passage. The presence of cytopathic effect (CPE) was also checked daily until the end of the experiment. The species of the captured bats were shown to be Dobsonia moluccensis, Acerodon celebensis, Pteropus sp.

(genetically close to Pteropus hypomelanus), and Pteropus vampyrus (Table 1) . Nested RT-PCR with three out of a total of 74 fecal samples from Dobsonia species produced a fragment of 455 bp in length after removal of the primer region, which was expected size of the RdRp fragment (region A in Fig. 1 and Table 1 ). The samples were IFB2012-8F, IFB2012-13F, and IFB2012-17F. The RdRp gene is the one most commonly used for PCR amplification in CoV surveillance, since it contains the most regions that are conserved in all CoVs [24] . Sequence analysis demonstrated that all samples contained the sequence of the RdRpgene-conserved motifs A (DYPKCD) and C (XSDD), which form the polymerase catalytic active site [24] . The partial RdRp sequences of IFB2012-13F and IFB2012-17F were identical. There were two nucleotide substitutions when the RdRp sequences of IFB2012-13F and IFB2012-17F were compared with the RdRp sequence of IFB2012-8F.

In addition to the partial RdRp gene, a portion of the helicase gene was also amplified, producing a fragment of 355 bp (region B in Fig. 1 ), from the three samples that were positive for PCR targeting the RdRp region. In the helicase gene, no nucleotide differences were found between IFB2012-13F and IFB2012-17F. However, two nucleotide substitutions were observed between IFB2012-8F and IFB2012-13F or IFB2012-17F.

After sequencing the partial RdRp and partial helicase regions, we designed new primers to amplify the sequence between these regions (region C in Fig. 1) . A 1335-bp fragment was detected in two samples (IFB2012-8F and IFB2012-17F). The sequence of a fragment (region D in Fig. 1 ), referred to as ''partial RdRp-helicase'', was then constructed by incorporating overlapping sequences of partial RdRp, partial helicase, and the region in between. The nucleotide sequence comparison of partial RdRp-helicase from IFB2012-8F and IFB2012-17F showed 13 nucleotide substitutions. BLAST searches showed that IFB2012-8F and IFB2012-17F had the highest nucleotide sequence identity to BatCoV HKU9-10-2 and BatCoV KY06, at 86 % and 85 %, respectively. Following removal of primer sequences from the partial sequence of the RdRp-helicase gene at the 5 0 and 3 0 end, the length of the sequence was 2,116 bp, and a phylogenetic tree comparing the sequence of the samples with known CoVs, including SARS-CoV and MERS-CoV, was constructed (Fig. 2) . The phylogenetic tree showed that IFB2012-8F and IFB2012-17F formed a distinct branch that was closely related to BatCoV HKU9 (NC009021), HKU9-2 (EF605514), HKU9-5-2 (HM211099), and HKU9-10-2 (HM211101) from China [25] as well as BatCoV KY06 from Kenya (HQ728483) [26] . All of these CoVs belong to the genus Betacoronavirus.

Nucleotide sequences of partial RdRp-helicase genes derived from fecal samples IFB2012-17F and IFB2012-8F were deposited under accession numbers AB918718 and AB918719, respectively, in the DNA Data Bank of Japan (DDBJ, Tokyo) nucleotide database.

Various tissues and tracheal swab samples collected from the bats where the feces yielded a positive result were also subjected to nested RT-PCR for the partial RdRp region. Spleen tissue of IFB2012-13F, brain tissue of IFB2012-17F, and lung tissue of IFB2012-17F showed positive signals for RdRp gene amplification. Sequence analysis of tissue samples showed sequences identical to the partial RdRp sequence obtained from fecal samples (IFB2012-13F and IFB2012-17F). The respiratory epithelium is known to be a target for CoVs [27] , and therefore, detection of bat CoV genome fragments in lung samples might be expected. However, the detection of CoV genome fragments in spleen and brain samples suggested that bat CoV infection in these species may not be restricted to the respiratory or alimentary tract.

Attempts at virus isolation using Vero E6, FBKT, and DemKT1 cells were unsuccessful. No CPE was observed in cell cultures, and viral RNA was not detected by RT-PCR from the culture supernatants after the fifth passage.

Nucleotide sequencing and phylogenetic analysis of partial RdRp-helicase regions of the samples in the present study and these of known CoVs suggest that the CoV genome fragments obtained in this study are related to [28] . Therefore, the CoVs detected in this study could not be defined at the species level, since only certain regions of conserved replicase domains were determined. However, the number of nucleotide substitutions over the partial RdRp-helicase sequence between IFB2012-8F and IFB2012-17F and phylogenetic analysis of partial RdRp-helicase suggested that the fruit bats harbor variants of a single coronavirus strain. In this study, CoV genomes containing regions resembling those of members of the genus Betacoronavirus were detected in Dobsonia moluccensis, known as Moluccan naked-backed fruit bat, in Indonesia. Dobsonia moluccensis belongs to family Pteropodidae, which is known to harbor a number of viruses, including CoVs [25, 29] , paramyxoviruses [18, 30] and rhabdoviruses [29, 31] . Nipah and Hendra viruses have also been detected in Dobsonia moluccensis [32, 33] . Fig. 2 Phylogenetic tree of partial RdRp-helicase using a 2,116-bp nucleotide sequence following removal of primer sequences at the 5 0 and 3 0 end of the partial RdRp-helicase gene obtained from fruit bats samples (IFB2012-8F and IFB2012-17F) and various coronaviruses derived from GenBank. The tree was constructed by the maximum-likelihood method using MEGA 6 software. Bootstrap values were calculated for 1000 replicates and are shown next to the branches. Bat coronaviruses detected in this study are indicated by black squares. SARS-CoV, severe acute respiratory syndrome coronavirus; avian IBV, avian infectious bronchitis virus; CoV, coronavirus This report highlights the fact that Indonesian fruit bats are a natural reservoir for various viruses, which have the potential to be zoonotically transmitted, since fruit bats are also consumed in certain parts of Indonesia.

What we are watching-five top global infectious disease threats, 2012: a perspective from CDC's Global Disease Detection Operations Center

Detection of coronaviruses in bats of various species in Italy

Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China

The Genome sequence of the SARS-associated coronavirus

Middle East respiratory syndrome coronavirus in bats, Saudi Arabia

Identification of a novel coronavirus in patients with severe acute respiratory syndrome

A novel coronavirus associated with severe acute respiratory syndrome

Characterization of a novel coronavirus associated with severe acute respiratory syndrome

Coronaviruses as DNA wannabes: a new model for the regulation of RNA virus replication fidelity

Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS

From SARS coronavirus to novel animal and human coronaviruses

Bats and their virome: an important source of emerging viruses capable of infecting humans

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor

Bat origins of MERS-CoV supported by bat coronavirus HKU4 usage of human receptor CD26

Genomic and serological detection of bat coronavirus from bats in the Philippines

SARS-Coronavirus ancestor's foot-prints in South-East Asian bat colonies and the refuge theory

Molecular detection of a novel paramyxovirus in fruit bats from Indonesia

Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription-PCR and nonfluorescent low-density microarray

Evolutionary trees from DNA sequences: a maximum likelihood approach

MEGA6: molecular evolutionary genetics analysis version 6.0

Characterization of the envelope glycoprotein of a novel filovirus, lloviu virus

Isolation of novel adenovirus from fruit bat

Molecular model of SARS coronavirus polymerase: implications for biochemical functions and drug design

Coexistence of different genotypes in the same bat and serological characterization of Rousettus bat coronavirus HKU9 belonging to a novel Betacoronavirus subgroup

Detection of novel SARS-like and other coronaviruses in bats from Kenya

Differential sensitivity of bat cells to infection by enveloped RNA viruses: coronaviruses, paramyxoviruses, filoviruses, and influenza viruses

Virus taxonomy: classification and nomenclature of viruses: Ninth Report of the International Committee on Taxonomy of Viruses

Coronavirus genomes in fruit bats in Indonesia 1117

Bats: important reservoir hosts of emerging viruses

Bats host major mammalian paramyxoviruses

Bat lyssaviruses, northern Vietnam

Nipah virus and bats

A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special?

Acknowledgements We thank Prof. Ayato Takada for providing the Vero E6 cell line, and Prof. Ken Maeda for providing the FBKT and DemKT1 cell lines for the current study. This study was partly supported by a grant from the Japanese Initiative of Global Research Network on Infectious Diseases (J-GRID) and a grant of the Ministry of Education, Culture, Sports, and Technology (MEXT), Japan.