key: cord-334628-axon4jdc authors: Lee, Saemi; Jo, Seong-Deok; Son, Kidong; An, Injung; Jeong, Jipseol; Wang, Seung-Jun; Kim, Yongkwan; Jheong, Weonhwa; Oem, Jae-Ku title: Genetic Characteristics of Coronaviruses from Korean Bats in 2016 date: 2017-07-19 journal: Microb Ecol DOI: 10.1007/s00248-017-1033-8 sha: doc_id: 334628 cord_uid: axon4jdc Bats have increasingly been recognized as the natural reservoir of severe acute respiratory syndrome (SARS), coronavirus, and other coronaviruses found in mammals. However, little research has been conducted on bat coronaviruses in South Korea. In this study, bat samples (332 oral swabs, 245 fecal samples, 38 urine samples, and 57 bat carcasses) were collected at 33 natural bat habitat sites in South Korea. RT-PCR and sequencing were performed for specific coronavirus genes to identify the bat coronaviruses in different bat samples. Coronaviruses were detected in 2.7% (18/672) of the samples: 13 oral swabs from one species of the family Rhinolophidae, and four fecal samples and one carcass (intestine) from three species of the family Vespertiliodae. To determine the genetic relationships of the 18 sequences obtained in this study and previously known coronaviruses, the nucleotide sequences of a 392-nt region of the RNA-dependent RNA polymerase (RdRp) gene were analyzed phylogenetically. Thirteen sequences belonging to SARS-like betacoronaviruses showed the highest nucleotide identity (97.1–99.7%) with Bat-CoV-JTMC15 reported in China. The other five sequences were most similar to MERS-like betacoronaviruses. Four nucleotide sequences displayed the highest identity (94.1–95.1%) with Bat-CoV-HKU5 from Hong Kong. The one sequence from a carcass showed the highest nucleotide identity (99%) with Bat-CoV-SC2013 from China. These results suggest that careful surveillance of coronaviruses from bats should be continued, because animal and human infections may result from the genetic variants present in bat coronavirus reservoirs. Coronaviruses are large, pleomorphic, enveloped viruses containing a single linear, positive-sense single-stranded RNA molecule. Coronavirus genomes are approximately 27-32 kb in length, the largest continuous RNA genomes among mammalian viruses [1] . Coronaviruses are classified into four genera: alpha, beta, gamma, and delta [2] . Coronaviruses are the second most prevalent cause of the common cold in humans [3] . In livestock and poultry, coronaviruses are recognized causes of enteric and respiratory infections that are often fatal in young animals [4] . An outbreak of severe acute respiratory syndrome (SARS) in 2002 and 2003 resulted in infection of 8096 people worldwide, with 774 (9.5%) of them dying from this novel human coronavirus [5] . SARS coronavirus is believed to have been acquired from an animal species, and the Chinese horseshoe bat may have been the original source of SARS infection [6] [7] [8] [9] . The emergent Middle East respiratory syndrome coronavirus (MERS-CoV) may also have originated in bats, with fragments of MERS-CoV genes being identified in bats from both Saudi Arabia [10] and Africa [11] . More than 200 viruses have been isolated from and detected in bats [12] . The ability of bats to fly and migrate, as well as the large sizes of their social groups, predisposes bats to acquire and maintain viruses [13] . Given the import of MERS into South Korea [14] and the presence of SARS in the relatively close geographic location of China [9] (Fig. 3) , together with the fact that bats are a reservoir for coronaviruses, the prevalence of coronavirus infection in Korean bat species should provide valuable information. Although SARS-like and MERS-like bat coronaviruses have already been reported in the feces of Korean bats [14] , the prevalence and presence of coronaviruses in other bat samples have not been described. Therefore, we investigated the prevalence of coronaviruses in Korean bat species using various kinds of samples, including oral swabs, fecal and urine samples, and carcasses. The coronaviruses found were characterized for their relationships to each other and to coronaviruses isolated in South Korea and in other countries. A total of 672 samples (332 oral swabs, 245 fecal samples, 38 urine samples, and 57 bat carcasses) were collected at 33 sites of natural bat habitat from January 2016 to September 2016 (Table 1 ). Bats were captured by net for collection of oral swabs, fecal samples, and urine samples, and released immediately after sampling. Fecal samples were also collected from guano. The bat carcasses collected were naturally dead bats found at the site where the other samples were taken. Swab samples were kept in viral transport medium at 4°C, and organs collected from carcasses and other samples were stored at −80°C before processing. All fecal and tissue samples were resuspended in 1% antibiotic-antimycotic solution (Corning, USA) diluted in phosphate-buffered saline (PBS), and clarified by centrifugation at 3500×g for 10 min. RNA from 200 μL sample was extracted with the QIAamp® Viral RNA Mini Kit (Qiagen, Germany) and eluted in 60 μL RNase-free water. cDNA was synthesized by PrimeScript First Strand cDNA Synthesis Kit (TaKaRa, Japan) following the manufacturer's instructions. Bat-CoV screening was performed by a pancoronavirus PCR method based on primers used by Poon et al. (Corona 1 forward, 5′-GGTTGGGACTATCC TAAGTGTGA-3′, and Corona 2 reverse, 5′-CCAT CATCAGATAGAATCATCATA-3′) [15] , followed by sequencing of the amplified product to confirm Bat-CoV identification. The pan-coronavirus primers were used to amplify and sequence a 440-bp segment of the highly conserved RNA-dependent RNA polymerase (RdRp) gene. Cycle sequencing reactions were performed using the BigDye® Terminator Cycle Sequencing kit version 1.1 (Applied Biosystems, Foster City, CA, USA). Reactions were purified using SigmaSpin™ Post-reaction Clean-Up Columns (Sigma, St. Louis, MO, USA) and sequenced on an ABI PRISM 3130 automated capillary DNA sequencer (Applied Biosystems) according to the manufacturer's instructions. Primer sequences were excluded from the raw sequence data, and the sequences we analyzed were 392bp long. All sequences obtained in this study were submitted to GenBank (accession numbers KY432454-432471). The nucleotide sequences were aligned and compared to 52 selected human and animal CoV sequences available from the GenBank database using ClustalW software implemented in BioEdit version 7.0.9.0. The phylogenetic trees were drawn using the neighbor joining method using the maximum composite likelihood model with MEGA 7 software. Coronaviruses were detected in 13 oral swabs, one carcass (intestine), and four fecal samples from a total of 672 samples, giving an overall detection rate of 2.7% (18/672). Four of the 14 bat species tested were found to harbor coronavirus. Thirteen positive oral swabs were detected from one species of the family Rhinolophidae, which accounted for 24.3% (163/672) of the samples, and the other positive samples were from Vespertilionidae bats that comprised 55.5% (373/672) of the samples. The 13 coronavirus-positive oral swabs were collected from bats at an abandoned mine in Jeonbuk province ( Fig. 1 ). Oral swabs and other samples (n = 60) were obtained from three species of bats, Rhinolophus ferrumequinum, Miniopterus schreibersii, and Myotis macrodactylus, but coronaviruses were only detected in samples from R. ferrumequinum (Table 2 ). One carcass was collected from a cave in Chungbuk province (Fig. 1 ). In total, 94 samples were obtained from six bat species (M. schreibersii, Myotis petax, Vespertilio sinensis, Murina leucogaster, R. ferrumequinum, and M. macrodactylus), and coronavirus was detected from a V. sinensis sample. Four coronaviruspositive fecal samples were collected in four different regions, in Gyeongbuk province (Andong, Yeongdeok, and Gyeongju) and in the metropolitan city Gwangju (Fig. 1 ). In the forest of Andong, 43 samples in total were collected from three species of bats (Eptesicus serotinus, M. petax, and Pipistrelus abramus), and coronavirus was detected from one P. abramus bat. In Yeongdeok, a total of five samples was collected from P. abramus living in the forest, and one sample was positive for coronavirus. In the forest of Gyeongju, five samples in total were collected only from P. abramus bats, and one of the samples was positive for coronavirus. In the forest of Gwangju, 18 samples in total were collected from two species of bats (E. serotinus and Myotis aurascens) and coronavirus was detected from one E. serotinus bat ( Table 2) . To determine the genetic relationships between the 18 bat coronaviruses obtained in this study and previously known coronaviruses, the 392-nt RdRp sequences were analyzed phylogenetically. All the sequences obtained in this study belonged to the betacoronavirus genus (Fig. 2) . The sequences (Fig. 2) . The 13 sequences from oral swabs were 97.1-100.0% identical with each other, and the closest known CoV strain, with 97.1-99.7% nucleotide identity, was Bat-CoV-JTMC15 that was recently isolated from R. ferrumequinum in China. The nucleotide identity with SARS-CoV was only 86.7-89.0%. One sequence from a carcass showed the highest identity (99%) with Bat-CoV-SC2013 that was isolated from Vespertilio superans (synonym of V. sinensis), also from China [16] . The nucleotide identity with MERS-CoV was 83.9%. Four sequences from fecal samples shared 96.9-100.0% identity with each other and showed 94.1-95.1% identity with Bat-CoV HKU5 [17] . The nucleotide identity with MERS-CoV was 84.1-84.4%. Due to shorter reported sequences from Kim et al. [14] , another phylogenetic tree was constructed from the 355-bp RdRp sequences of Korean bat coronaviruses obtained by Kim et al. [14] and in this study to analyze the genetic relation between Bat-CoVs in South Korea (Fig. 3) . Thirteen sequences from oral swabs were clustered with Bat-CoV B15-21, which was detected in fecal bat samples collected from an abandoned mine in Gangwon province. The major bat species present there was R. ferrumequinum, but the coronaviruspositive bat species were not specified [14] . Nucleotide identity between Bat-CoV B15-21 and the 13 clustered sequences was 97.7-100.0%. The other three Korean bat coronaviruses (Bat-CoV B15-8, B15-40, and B15-41) belonged to the alphacoronavirus genus, but we did not detect coronaviruses belonging to this group in our study. There are approximately 24 species of bats in South Korea, and most are insectivorous and relatively small [18] . Bats in South Korea usually do not migrate over distances greater than 1000 km [19], but apart from two species found only in South Korea, they are generally distributed in neighboring countries such as China, East Siberia, Japan, and Taiwan. Although it is not clear if bats go between China and other countries, some species that migrate over 500 km can carry pathogens to South Korea. The fact that the 18 bat coronaviruses (Bat-CoVs) detected in this study were all In this study, the prevalence of coronavirus in Korean bat samples was 2.7% overall. The prevalence by sample type was 5.4% (13/332) for oral swabs, 1.8% (1/57) for carcasses, 1.6% (4/245) for fecal samples, and 0.0% (0/38) for urine samples. Reported rates from other countries for coronavirus prevalence in bats were 6.5% (64/985) and 5.3% (50/951) for China, 8.2% (16/195) for Italy, and 5.3% (32/606) for Mexico [20] [21] [22] [23] , not greatly different from our figure of 2.7% for prevalence in South Korea. Fourteen out of 24 Korean bat species were screened for CoVs in this study. The species positive for Bat-CoVs were R. ferrumequinum, V. sinensis, P. abramus, and E. serotinus. The Bat-CoVs that showed the highest nucleotide identity with our isolates were obtained from the same species that we identified as positive, except for E. serotinus. To our knowledge, this is the first time a Bat-CoV has been detected from E. serotinus in Asia. In Europe (Italy), a lineage C betacoronavirus was detected from E. serotinus in 2014 [24] . Bat-CoVs were usually detected at higher rates in alimentary samples (fecal samples, rectal swabs, and intestine) than in respiratory samples (oral and throat swabs) [7, 15, 17, 22, 25, 26] . We collected alimentary and respiratory samples and urine samples. However, Bat-CoVs were not detected in urine samples. Animal coronaviruses usually produce either enteric or respiratory infections. Animal models show similar clinical features to SARS and MERS patients; the majority of whom presented with respiratory symptoms but in some cases also suffered from enteric complications [27] . Unlike other human coronavirus infections, a number of MERS cases were associated with acute renal failure. MERS-CoV replication in the human kidney suggests the potential of shedding viruses in urine [28] . Bat-CoV isolate Bat SL-CoV-WIV1 (KF367457) showed infectivity to Rhinolophus sinicus kidney cell lines, which means it may be possible to detect Bat-CoV in urine [6] . The current study revealed that Korean bats have SARSand MERS-related Bat-CoVs, consistent with earlier work [ 1 4 ] . T h e p r e v i o u s s t u d y d e t e c t e d a l p h a -a n d betacoronaviruses only from fecal samples in four regions, two in Gangwon, one in Chungbuk, and one in Gyeongbuk [14] . The only betacoronavirus RdRp sequence Bat-CoV B15-21 (KU528590.1) was clustered with 13 Bat-CoVs β-CoV Lineage B β-CoV Lineage C α-CoV Fig. 3 Phylogenetic tree of Korean bat coronaviruses detected by Kim et al. [13] and in this study, based on 355-nt sequences. Phylogenetic trees were constructed by using the neighbor joining method and bootstrap values were determined by 3000 replicates. Scale bar the estimated genetic distance of these viruses. Black circle sequences from oral swabs. Black up-pointing triangle sequences from fecal samples. White uppointing triangle sequences from carcass (16BO122, 129, 132, 133, 134, 135, 139, 141, 143, 145, 150 , 154, and 155) detected from Danyang, Chungbuk province. The identity between the 14 sequences was very high at 97.7-100%, even though the sample types and collection sites were different. The species of bat infected with Bat-CoV B15-21 was not specified in the paper, but the major bat species of the sample collection site was R. ferrumequinum; the same species in which we detected 13 Bat-CoVs. The same Bat-CoVs were detected from fecal samples in the prior study [14] and from oral swabs in our study. Spike proteins, which are translated from spike genes, define viral tropism by receptor specificity and by membrane fusion activity during entry into cells [29] . Although the spike genes were not identified in this study, full spike gene analysis of Bat-CoV B15-21 suggested a low potential for human emergence. One MERS-like betacoronavirus, Bat-CoV B15-1-6 (KU528586.1, KU528585.1), was identified by Kim et al. [14] , but the similarity of this isolate to our isolates could not be analyzed because the RdRp sequence was not obtained. High nucleotide identity between Bat-CoV B15-1-6 and Bat-CoV-SC2013 (KJ473821) suggests that Bat-CoV B15-1-6 and 16BT3 might be very similar. Bat-CoVs have been identified in many countries including the Australasian region [10, 11, 14-16, 20-26, 30, 31] . But Bat-CoVs identified in Korea were genetically close to Bat-CoVs from China, a close neighbor of Korea. According to the results of Bat-CoV research in China, the degree of host restriction for coronaviruses in bat populations was high. For example, similar viruses were detected in the same species of bat in different regions, approximately 1600 km apart, while two different bat species in the same habitat had different coronaviruses [22] . In our results, three Bat-CoVs (16BF109, 211, and 244) from P. abramus were clustered with each other, even though samples were collected from different habitats in different regions (Andong, Yeongdeok, and Gyeongju). The Bat-CoV detected from E. serotinus bat in the metropolitan city Gwangju (16BF104) was also clustered with these three Bat-CoVs mentioned above. However, there was no matching virus for 16BF104 found in the same species of bat in other locations, nor in other species in the region where samples were collected. Our data could support the possibility of relationships between genetic variation, geographic locations where Bat-CoVs were isolated, and bat host species, but more data on Bat-CoVs in Korea would be needed to address this. In summary, the prevalence of coronavirus in Korean bats was determined and found to be comparable to or lower than the rates in other countries. We detected SARS-related and MERS-like Bat-CoVs that were genetically very similar to Bat-CoVs identified in China. For the study of Bat-CoV seasonal occurrence and transmission, continuous monitoring during all seasons is required, and ideally alimentary, respiratory, and other samples need to be collected from each bat. Future work to determine the complete genome sequences of Bat-CoVs from South Korea would give a more complete picture of their pathogenicity and the possibility of infectious spread to other animals and humans. Bat CoV-JTMC15(KU182964) Bat CoV-273(DQ648856) SARS Bat CoV-Lingbao(KF294456) SARS Bat CoV-Rf1(DQ412042) Bat CoV-Rp(JX993987) Bat CoV HKU3-1(DQ022035) Bat CoV-Rp3(DQ071615) Bat CoV RsSHC014(KC881005) SARS CoV-MA15(JF292915) SARS CoV-Tor2(JX163928) SARS CoV-Urbani(AY278741) SARS-CoV(NC004718) Bat CoV BM48-31/BGR/2008(GU190251) Bat CoV SLO1A0082(GQ404797) PHEV(NC007732) MHV(NC001846) HCoV-HKU-1(NC006577) Bat CoV-HKU9(NC009021) TCV(NC010800) IBV(NC001451) Bat CoV-HKU4(NC009019) Kbat CoV-16BF104 Kbat CoV-16BF109 Kbat CoV-16BF221 Kbat CoV-16BF244 Bat CoV HKU5-1(DQ249217) Bat CoV-SC2013(KJ473821) Kbat Cov-16BT3 Bat CoV(KC243390) Bat CoV BtCoV/UKR-G17(KC243392) MERS-CoV(NC019843) Human beta-CoV(JX869059) Human beta-CoV(KC667074) Human beta-CoV(KC776174) Coronavirus genome structure and replication A comparative sequence analysis to revise the current taxonomy of the family Medical microbiology Animal coronaviruses: what can they teach us about the severe acute respiratory syndrome? Ecology, evolution and classification of bat coronaviruses in the aftermath of Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor Bats are natural reservoirs of SARSlike coronaviruses Science Veterinary microbiology and microbial disease A review of studies on animal reservoirs of the SARS coronavirus Middle East respiratory syndrome coronavirus in bats Close relative of human Middle East respiratory syndrome coronavirus in bat Bats and zoonotic viruses: can we confidently link bats with emerging deadly viruses? Bats: important reservoir hosts of emerging viruses Detection of severe acute respiratory syndrome-like, Middle East respiratory syndrome-like bat coronaviruses and group H rotavirus in Faeces of Korean bats Identification of a novel coronavirus in bats MERS-related betacoronavirus in Vespertilio superans bats Genetic characterization of betacoronavirus lineage C viruses in bats reveals marked sequence divergence in the spike protein of pipistrellus bat coronavirus HKU5 in Japanese pipistrelle: implications for the origin of the novel Middle East respiratory syndrome coronavirus Bats. Moonji publishing, Seoul 19 Detection of coronaviruses in bats of various species in Italy Viruses Prevalence and genetic diversity of coronaviruses in bats from Detection and characterization of diverse alpha-and betacoronaviruses from bats in Alpha and lineage C betaCoV infections in Italian bats Coronaviruses in bent-winged bats (Miniopterus spp Detection of alpha and betacoronaviruses in multiple Iberian bat species Microbiology, ecology, and natural history of coronaviruses In-vitro renal epithelial cell infection reveals a viral kidney tropism as a potential mechanism for acute renal failure during Middle East respiratory syndrome (MERS) coronavirus infection Coronavirus spike proteins in viral entry and pathogenesis Coronavirus infection and diversity in bats in the Australasian region EcoHealth Detection of group 1 coronaviruses in bats in North America Emerg Acknowledgements We thank Dr. C.W. Jeong and his colleagues for their efforts in the collection of wild bat carcasses and samples. This research was supported by grant no. 2016-01-01-033 from the NIER of the Republic of Korea. The funders had no roles in the study design, data collection and analysis, decision to publish, or the preparation of the manuscript.