key: cord-0929754-thraojqm authors: Dharmayanti, Ni Luh Putu Indi; Nurjanah, Diana; Nuradji, Harimurti; Maryanto, Ibnu; Exploitasia, Indra; Indriani, Risa title: Molecular detection of bat coronaviruses in three bat species in Indonesia date: 2021-08-09 journal: J Vet Sci DOI: 10.4142/jvs.2021.22.e70 sha: c4c72bdb7f1ff324480b2fb524dcdc54f0321a71 doc_id: 929754 cord_uid: thraojqm Bats are an important reservoir of several zoonotic diseases. However, the circulation of bat coronaviruses (BatCoV) in live animal markets in Indonesia has not been reported. Genetic characterization of BatCoV was performed by sequencing partial RdRp genes. Real-time polymerase chain reaction based on nucleocapsid protein (N) gene and Enzyme-linked immunosorbent assay against the N protein were conducted to detect the presence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral RNA and antibody, respectively. We identified the presence of BatCoV on Cynopterus brachyotis, Macroglossus minimus, and Rousettus amplexicaudatus. The results showed that the BatCoV included in this study are from an unclassified coronavirus group. Notably, SARS-CoV-2 viral RNA and antibodies were not detected in the sampled bats. Bats have important roles as natural reservoir hosts of zoonotic diseases and act as agents of zoonosis transmission, including the transmission of coronaviruses that led to severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) [1] . Bats have good adaptability skills, which support their resistance against several viral infections [2] . Studies in Indonesia have revealed that bats from several regions harbor Hendra and Nipah viruses [3, 4] . There are more than 1,200 species of bats, and approximately 20% (239 species) have been reported in Indonesia [5] [6] [7] . This high diversity of bat species provides a wide variety of cell types and receptors that may facilitate the bat species becoming a potential source for coronavirus transmission [8] . of viruses into the environment [11] . Wild animal meat ("bushmeat") comes from various sources, including bats. Transmission of bushmeat-associated pathogen occurs through direct contact with the body fluids and feces of the wild animals. Ebola virus, human immunodeficiency virus-1, monkeypox virus, and SARS-CoV are infectious pathogens associated with bushmeat [12] [13] [14] . SARS-CoV and MERS-CoV are zoonotic coronaviruses that believed originating from bats and have been transmitted to other species (civets and camels, respectively), which act as intermediate hosts before being transmitted to humans [15] [16] [17] . The SARS-CoV-2, which recently caused a worldwide pandemic, was shown to have a nucleotide similarity of up to 96.1% with BatCoV RaTG13 originating from Rhinolophus affinis and first isolated in 2013 [18, 19] . Various types of interactions between bats, other animals, and humans may contribute to the interspecies transmission of coronaviruses, including interactions that may occur in live animal markets [15, 20, 21] . Live animal markets are located in most regions of Indonesia, 2 of which are the Tomohon market in Manado, North-Sulawesi and the Depok market in Surakarta, Central Java, where people sell wild animals for consumption and for other purposes [22] . Live animal markets are potential meeting points for the interspecies transmission of several diseases from animals to human or vice versa (i.e., zoonotic diseases) [12] [13] [14] 23] . In 2015, Anindita et al. [24] reported a bat betacoronavirus (Bat betaCoV) in Dobsonia moluccensis acquired in Paguyaman, Gorontalo Province, Indonesia. Furthermore, Febriani et al. [25] reported that Pteropus alecto in Gorontalo Province carried Bat betaCoV. This study reports for the first time the presence of bat coronaviruses (BatCoV) in 3 species of bats sold by traders at live animal markets in Central Java Province, and bat collectors in West Java Province and Yogyakarta Province, Indonesia. The results of this study are expected to contribute to elucidating BatCoV ecology in Indonesia and identify the species of bats that are potential hosts of BatCoV. Identification of bat species was performed using photo-documentation and external morphological (morphometric) measurements based on various indicators, such as forearm length, tibia length, hind leg length, ear length and shape, body length, body color, presence of claw on the second finger of the wing, tail length, body weight, shape of the muzzle and tongue, and color on the edge of the ear [7] . Enzyme-linked immunosorbent assay (ELISA) for the detection of antibodies directed against the nucleocapsid protein of SARS-CoV-2 ELISA was carried out using ID Screen SARS-CoV-2 Double Antigen Multi-species kits. Briefly, 25 μL of dilution buffer 13 was added to each well, 25 μL of the negative control solution was added to wells A1 and B1, 25 μL of the positive control solution was added to wells C1 and D1, and 25 μL of each sample to be tested were added to the remaining wells. The plates were covered and incubated for 45 min ± 5 min at 37°C (± 2°C). The wells were then emptied and washed 3 times with at least 300 μL of the kit's wash solution without drying the wells between washes. Conjugate 1X was prepared by diluting the concentrated conjugate 10-fold to 1:10 in dilution buffer 13, and 100 μL of conjugate 1X was added to each well. The plate was again covered and incubated for 30 min ± 3 min at 21°C (± 5°C). The wells were then emptied and washed 3 times with at least 300 μL of wash solution, avoiding drying of the wells between washes. Substrate solution (100 µL) was added to each well, and the plate was then covered and incubated for 20 min ± 2 min at 21°C (± 5°C) in the dark. Stop solution (100 µL) was added to each well, in the same order as described above to stop the reaction. Optical density was then read and recorded at 450 nm and the S/P percentage (S/P%) was calculated as: Where OD PC is optical density of positive control, OD NC is optical density of negative control. Samples presenting a S/P% less than or equal to 50% were considered negative. Those between 50% and 60% considered as doubtful, and those greater than or equal to 60% considered positive (interpretation: ≤ 50%, negative; 50%-60%, doubtful; ≥ 60%, positive). A 25-μL reaction mixture was established containing 5 μL of RNA and 12.5 μL of 2 × reaction buffer as provided with the Superscript III one-step RT-PCR system with Platinum Taq Polymerase (Invitrogen, USA) containing 0.4 mM of each deoxyribonucleotide triphosphates and 3.2 mM magnesium sulfate, 1 μL of reverse transcriptase/Taq mixture from the kit, and 1 μL for each primer probe (N1, N2, and RP). Thermal cycling was performed at 50°C for 15 min for reverse transcription, followed by 95°C for 2 min and 45 cycles of 95°C for 15 sec, 55°C for 30 sec using an Applied Biosystems 7300 Real-Time PCR System (Applied Biosystems, USA). The isolation of BatCoV genetic material obtained from individual field samples was carried out using The QIAamp Viral RNA Mini Kit (Qiagen, Germany) following the kit instructions. The extracted RNA was stored at −20 o C for RT-PCR. Genetic material identification was performed using primers designated for partial RdRp genes as described by Woo et al. [13] . The partial RdRp gene was amplified by performing RT-PCR using the AB9700 thermal cycler. The PCR mixture (25 µL) contained 12.5 µL of 2X reaction mix, 1 µL of each forward and reverse primer (20 pmol/µL), 0.5 µL SuperScript III RT/Platinum Taq Mix, and 10 µL of RNA template. The initial stage of RT-PCR was reverse transcription at 42°C for 45 min, followed by enzyme inactivation at 95°C for 3 min. Amplification included 35 cycles of denaturation for 30 sec at 95 o C, annealing for 40 sec at 48 o C, and elongation for 40 sec at 72 o C. Final elongation was carried out at 72 o C for 10 min. Amplicon products were visualized using 2% gel agarose, with gel electrophoresis run for 30 min at 100 V. Partial DNA sequencing PCR products were separated in 1% agarose by electrophoresis and the amplicon was excised and purified using a QIAquick gel purification kit (QIAGEN, Hilden, Germany). The sequencing method used was direct sequencing using a cycle sequencing kit (BigDye Terminator version 3.1; Applied Biosystems) on a Genetic Analyzer 3130 (Applied Biosystems). The nucleotide sequencing data obtained in this study were analyzed together with the genetic data available in the National Center for Biotechnology Information for the coronavirus RdRp gene. The production of multiple alignments of the gene and the analysis of the amino acids were carried out using BioEdit version 7 (http://www.mbio.ncsu.edu/ BioEdit). Phylogenetic trees were generated by maximum likelihood (1,000 replicates) using the Tamura-Nei algorithm in MEGA version 5.2 (http://www.megasoftware.net). All isolates in this study were submitted to GenBank (www.ncbi.nlm.nih.gov) and given accession numbers MW652309-MW652323 and MZ451148-MZ451156. A total of 182 rectal swab samples from 4 species of bats were acquired from bat collectors who were going to sell the bats to restaurants for consumption and from animal markets in several regencies/cities in Central Java Province (Surakarta City, Magelang Regency), Yogyakarta Province, and West Java Province (Bogor City, Cianjur Regency) in 2020 and 2021. Seventy-two BatCoV-positive bats (39.56%), based on PCR targeting the RdRp region, were detected among the 182 samples from 3 bat species: C. brachyotis, M. minimus, and R. amplexicaudatus ( Table 1) Related to close contact with humans, serological and rRT-PCR investigations were also conducted to investigate the possible circulation of SARS-CoV-2-related coronaviruses in bats in Indonesia. The samples tested produced negative results in both the ELISA tests to detect antibodies against the N protein of the SARS-CoV-2 and the rRT-PCR assays targeting the nucleocapsid protein (N) gene of SARS-CoV-2 (Tables 3 and 4) . Bats are widely distributed globally, although each species has a specific geographic area [26] . Woo et al. [27] reported that only alpha-and betaCoV have been identified in bats. Bat betaCoV was previously reported on D. moluccensis collected from the Paguyaman Regency, Gorontalo Province, Indonesia. However, D. moluccensis is a native bat species from Moluccas, Indonesia [6] . The phylogenetic tree analysis indicated that the BatCoVs in the previous study formed distinct branches but were closely related to BatCoV HKU9, HKU9-2, HKU9-5-2, and HKU9-10-2 originating from China and BatCoV KY06 from Kenya [24] . In 2018, BatCoVs were detected in P. alecto from the Gorontalo Province, Indonesia. Those isolates (INDSWBT101, INDSWBT102, INDSWBT103) have 98% nucleotide similarity with Indonesian isolates, namely Bat Coronavirus Indonesia IFB 2012-8F D. moluccensis [24] , meanwhile the 5 other isolates (INDSWBT110, INDSWBT192, INDSWBT198, INDSWBT180, INDSWBT195) have 85%-88% nucleotide similarity with Bat Coronavirus BtCoV/B55762/S. hea/CB/Tha/6/2012 isolated from Scotophilus heathii in Thailand [25] . Indonesia is reported to have 81 species of bats that belong to suborder Megachiroptera (fruit nectar-eating bats) and 158 species in the suborder Microchiroptera (non-fruit-nectar-eating bats). As many as 10 families of bats have been reported in Indonesia, including Pteropodidae, which belongs to the Megachiroptera suborder, and Rhinopomatidae, Emballonuridae, Nycteridae, Megadermatidae, Rhinolophidae, Hiposideridae, Vespertilionidae, Minioptereridae, and Molossidae families that belong to the Microchiroptera suborder. As an archipelago country, Indonesia has the highest number of fruit-and nectar-eating bat species in the world. Fruit-nectar bats are very important for the process of pollination, fertilization, and seed dispersal. There are at least 21 genera and 81 species that belong to Pteropodidae family [6, 7, 28] , and, in this study, 3 species were observed to carry BatCoV (C. brachyotis, R. amplexicaudatus, and M. minimus). The discovery of BatCoV in several species of bats in Indonesia suggests that bats are potential natural reservoir hosts for coronavirus and may transmit that virus to other species, such as occurred in previous infectious disease outbreaks such as SARS-CoV and MERS-CoV [23, 29] . Bats are reported to possess efficient and varied antiviral responses associated with adaptations in their immune system and their ability to evolve. The adaptive immune mechanism in bats can suppress the pathological effects of the inflammation caused by viral infection. However, various factors, such as stress, may contribute to unbalancing the mechanism, resulting in increased viral replication and shedding, and potentially becoming a source of cross-species virus transmission, including human transmission [11] . Live animal markets are considered a major area for the spread of zoonotic diseases. Several zoonotic diseases have originated from live animal markets, including SARS-CoV and avian influenza virus [30] . Transportation, cages, and environmental conditions at live animal markets may trigger stress responses in various animals, including bats. Under stress conditions, animals tend to be more vulnerable to viruses and infections. In addition, cage systems, which are often stacked, may facilitate efficient transmission of the virus among animals [31] . In Indonesia, bats and other animals are widely traded in several live animal markets. During studies that we conducted in several animal markets in Indonesia, we observed that bats are kept close to other animals, including various species of poultry, other birds, as well as various mammals and reptiles. Collectors capture bats in forests or caves and sell them to traders in the animal markets and to restaurant owners. Bats may be held together with other animals for several days to months prior to being sold (Fig. 3) . Bat meat is commonly consumed among local communities due to traditional beliefs that such meat can increase stamina and cure particular diseases, such as asthma [22] . Bushmeat trade-related activities are considered to be connected to the transmission of zoonotic diseases. In addition to the people who consume bushmeat, people who slaughter, cut the body or organs, or have direct contact with the blood, urine, and feces, and those who may be scratched or bitten by wild animals have a high potential of being infected by zoonotic disease agents. Of all emerging diseases, zoonoses have been identified as the most significant increasing threat to global health [14] . 3 . Description of how bats in this study were captured, transported, and sold by bat collectors. Collectors captured various species of bats from forests and caves, transported the bats to a markets or their home by motorcycle. People/buyers purchased bats from traders at live animal markets or directly from the collector's house. Bats in the anthropocene Bat flight and zoonotic viruses Seroepidemiologi Nipah virus pada kalong dan ternak babi di beberapa wilayah di Indonesia Status infeksi virus hendra pada kalong (Pteropus Spp.) Di Pontianak Global patterns in coronavirus diversity Checklist of The Mammals of Indonesia: Scientific, English, Indonesia Name and Distribution Area Table in Indonesia Including CITES, IUCN and Indonesian Category for Conservation Global epidemiology of bat coronaviruses Response to the Ebola-related bushmeat consumption ban in rural Côte d'Ivoire A study on use of animals as traditional medicine by Sukuma Tribe of Busega District in North-western Tanzania Immune system modulation and viral persistence in bats: understanding viral spillover Serologic and behavioral risk survey of workers with wildlife contact in China Molecular diversity of coronaviruses in bats Bushmeat and emerging infectious diseases: lessons from Africa Human infection with MERS coronavirus after exposure to infected camels, Saudi Arabia Evidence for zoonotic origins of Middle East respiratory syndrome coronavirus Bats, civets and the emergence of SARS Genomic characterization of the 2019 novel humanpathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan A pneumonia outbreak associated with a new coronavirus of probable bat origin Organ distribution of severe acute respiratory syndrome (SARS) associated coronavirus (SARS-CoV) in SARS patients: implications for pathogenesis and virus transmission pathways A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster Quantifying the bat bushmeat trade in North Sulawesi, Indonesia, with suggestions for conservation action Review of bats and SARS Detection of coronavirus genomes in Moluccan naked-backed fruit bats in Indonesia Bat coronavirus of Pteropus alecto from Gorontalo Province, Indonesia Field identification key and guide for bats of the United States of America Discovery of seven novel mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus Kelelawar Sulawesi Jenis dan Peranannya. Jakarta: Penerbit Badan Penelitian dan Pengembangan Kesehatan Genomic characterization of a newly discovered coronavirus associated with acute respiratory distress syndrome in humans Wildlife, exotic pets, and emerging zoonoses Illicit wildlife trade, wet markets, and COVID-19: preventing future pandemics The authors express their appreciation and thanks to Teguh S, Ace ES, Suraida Meisari DVM, Zaki Aminullah DVM for technical assistance and the agricultural services of the respective regencies/cities staff for field activity assistance. Ni Luh Putu Indi Dharmayanti and Diana Nurjanah are the main authors and contributed equally to this article. NegativeNegative This study shows that 3 species of bats, C. brachyotis, M. minimus, and R. amplexicaudatus, that were collected by bat traders/bat collectors and examined in this study were found to carry viruses of an unclassified Bat coronavirus phylogenetic group. The ELISA results did not reveal antibodies against the nucleocapsid protein of the SARS-CoV-2 and the rRT-PCR assay results did not show the presence of SARS-CoV-2 viral RNA. Close contacts between humans, bats, and other animals have a high potential for transmitting zoonotic diseases. Further studies regarding coronaviruses carried by bats or other animals and the possible effects of environmental conditions are needed to identify possible novel virus transmission routes, particularly in live animal markets. Early detection of pathogen transmission and the application of appropriate control measures may minimize the destructive impacts on global health. Thus, more surveillance studies are needed to investigate the potentially important role of bats as natural reservoir hosts in the interspecies transmission of coronaviruses.