key: cord-294698-mtfrbn87
authors: Kim, H. K.; Yoon, S.‐W.; Kim, D.‐J.; Koo, B.‐S.; Noh, J. Y.; Kim, J. H.; Choi, Y. G.; Na, W.; Chang, K.‐T.; Song, D.; Jeong, D. G.
title: Detection of Severe Acute Respiratory Syndrome‐Like, Middle East Respiratory Syndrome‐Like Bat Coronaviruses and Group H Rotavirus in Faeces of Korean Bats
date: 2016-05-23
journal: Transbound Emerg Dis
DOI: 10.1111/tbed.12515
sha: 
doc_id: 294698
cord_uid: mtfrbn87

Bat species around the world have recently been recognized as major reservoirs of several zoonotic viruses, such as severe acute respiratory syndrome coronavirus (SARS‐CoV), Middle East respiratory syndrome coronavirus (MERS‐CoV), Nipah virus and Hendra virus. In this study, consensus primer‐based reverse transcriptase polymerase chain reactions (RT‐PCRs) and high‐throughput sequencing were performed to investigate viruses in bat faecal samples collected at 11 natural bat habitat sites from July to December 2015 in Korea. Diverse coronaviruses were first detected in Korean bat faeces, including alphacoronaviruses, SARS‐CoV‐like and MERS‐CoV‐like betacoronaviruses. In addition, we identified a novel bat rotavirus belonging to group H rotavirus which has only been described in human and pigs until now. Therefore, our results suggest the need for continuing surveillance and additional virological studies in domestic bat.

Bats are considered reservoirs of several emerging infectious disease. Recently emerging human viruses, such as severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), Nipah virus, Hendra virus and Ebola virus, are thought be bat-borne viruses (Han et al., 2015) . In addition, new species of influenza A viruses, lyssaviruses, paramyxoviruses and coronaviruses continue to be discovered in bats around the world (Tong et al., 2013; Banyard et al., 2014; Yang et al., 2014; Mortlock et al., 2015) . Reassortant group A rotavirus was also detected in the faeces of a straw-coloured fruit bat (Esona et al., 2010) . Therefore, viruses that originate in bats may be a source of additional spill over from wildlife into domestic animals and humans (Plowright et al., 2015) .

Since the outbreak of MERS-CoV in Korea in 2015, the prediction and control of newly emerging viruses has gained urgency. However, in Korea, there have been a few studies on viruses in domestic bat species. Serological evidence of Japanese encephalitis virus (JEV) was reported in the early 1990s (Lee and Lee, 1992) . In 1995, there was a report on the isolation and genetic characterization on a hantavirus from bats (Jung and Kim, 1995) . More information of viruses in bats is needed to clarify the role of these animals in infectious diseases epidemiology.

Therefore, in this study, we investigated viruses in bat species in Korea, using 49 faecal samples collected from July to December 2015 in 11 sites in natural bat habitats.

From July to December in 2015, 49 faecal samples were collected at 11 sites in natural bat habitats, such as caves, an abandoned mine and under a bridge (Table 1) . Fresh bat faeces were placed into transport medium in a 10% suspension and were transported to the laboratory for further laboratory analysis.

The major bat species at the collection sites were determined based on morphology and on previous data from bats' roosting sites (Han et al., 2012) . In this study, the major bat species in the collecting sites were Rhinolophus ferrumequinum, Myotis macrodactylus, Myotis aurascens Kuzyakin, Myotis petax and Miniopterus schreibersii.

Consensus primer-based reverse transcriptase polymerase chain reactions (RT-PCRs) were performed to detect influenza A viruses, coronaviruses, lyssaviruses and paramyxoviruses (Fouchier et al., 2000; Poon et al., 2005; V azquez-Mor on et al., 2006; Tong et al., 2008) . The detailed information on the RT-PCRs is in Table S1 . A total of 49 bat faecal samples were tested. The amplicons from positive RT-PCRs were sequenced using target-specific forward and reverse primers synthesized by Cosmogenetech Co. Ltd (Daejeon, Korea). The nucleotide data obtained from Sanger sequencing of the PCR fragments were further analysed with related sequences in GenBank using BioEdit (Hall, 1999) and MEGA version 6 (Tamura et al., 2013) .

In total, 16 bat faecal samples that were collected in four sites were used for high-throughput sequencing ( Table 1) . The supernatants from suspensions of one, four, five and six faecal samples collected in Danyang, Wonju, Taebaek and Uljin, respectively, were pooled by sampling site. Each of the four pooled samples was filtered through a 0.2-lm filter and ultracentrifuged at 252 000 g for 1 h. Each pellet was resuspended in 500 ll of 19 digestion buffer (Turbo DNA Free Kit; Ambion, Darmstadt, Germany), and a total of 10 ll of Turbo DNase was added. After incubating for 30 min at 37°C, suspensions were transferred to a 1.5-ml reaction tube and centrifuged at 1500 g for 3 min. The supernatants were used for RNA extraction. RNA was extracted using Trizol LS (Invitrogen Corp., Carlsbad, CA, USA), following the manufacturer's manual. The extracted RNA was submitted to Macrogen (Seoul, Korea) for high-throughput sequencing in a HiSeq 2000 sequencing system based on the transcriptome de novo sequencing platform.

The obtained viral contigs were analysed and annotated by the MG-RAST server (Meyer et al., 2008) . The cut-off for the annotation was 10 À5 , 60% and 15 for maximum e-value, minimum percentage identity and minimum alignment length, respectively. The contigs annotated as mammalian viruses were validated through BLASTn (http:// www.ncbi.nlm.nih.gov). The GenBank accession numbers for the nucleotide sequences in this study are KU528584-KU528594.

There were no samples positive for influenza A viruses, lyssaviruses or paramyxoviruses.

Among 14 cases of the submitted samples, three cases were positive for coronavirus by RT-PCR, showing the expected 440 nt target. These included the following: (Case no. 1) one of four samples from HO cave in Wonju where Rhinolophus ferrumequinum and Myotis macrodactylus were the major species present, (Case no. 6) one of six samples BM abandoned mine in Inje where Rhinolophus ferrumequinum was the major species found and (Case no. 12) two of four samples BGS (Case no. 12) cave in Munkyung where Miniopterus schreibersii was found to be the major species (Table 1) .

In the HO cave, one of four faecal samples (B15-8) was positive by coronavirus-specific RT-PCR. When the PCRpositive fragment was sequenced, it was found to be a bat coronavirus (Bat CoV) High-throughput sequencing from the pooled bat faecal samples As shown in Table 1 , cases 1, 2, 3 and 4 were further analysed by high-throughput sequencing. By Hiseq 2000 sequencing, 4341, 162 559, 7690 and 16 657 contigs were obtained from cases 1, 2, 3 and 4, respectively. Based on the annotation by the MG-RAST system, several viral sequences that are related to mammalian viruses, insect viruses, plant viruses, fungal viruses and bacteriophages were identified.

The mammalian viral sequences annotated by the MG-RAST system were further validated by BLASTn (http:// www.ncbi.nlm.nih.gov). The BLASTn result is presented in Table 2 . Fifteen viral sequences from case no. 1 were similar to porcine rotavirus SKA-1 (69-73% nucleotide identity), Bat CoV SC2013 (70-91% nucleotide identity), Bat CoV HKU5-5 (82% nucleotide identity), Alphacoronavirus Eptesicus fuscus-related strains (95-96% nucleotide identity) and Rhinolophus pearsoni bunyavirus Shaanxi 2011 (72% nucleotide identity). Three viral sequences from case no. 2 were similar to Banna virus 02VN078b (90% nucleotide identity), Bat CoV CDPHE15 (80% nucleotide identity) and Bat CoV Neixiang-31 (82% nucleotide identity). There was no mammalian viral sequences in cases 3 and 4.

Partial RNA-dependent RNA polymerase (RDRP) sequences obtained from B15-8 (sample in case no. 2), B15-21 (sample in case no. 6), B15-40 and B15-41 (sample in case no. 12) were phylogenetically compared with other reference sequences (Fig. 1a) . Three of these sequences clustered with the alphacoronaviruses (B15-8, 40 and 41 Because an RDRP sequence from betacoronavirusrelated Bat CoV B15-1-6 (detected from the pooled sample in case no. 1) was not obtained, a partial spike gene sequence (c10151) was compared with the spike genes of reference betacoronaviruses (Fig. 1b) . In addition, the full spike gene sequence obtained from Bat CoV B15-21 was also analysed. Bat CoV B15-1-6 grouped with MERS-CoVlike Bat CoVs (Bat CoV strains HKU5, A434, SC2013 and HKU4), showing 80.9-82.2% nucleotide identities, and was found to be closely related to MERS-CoV EMC, showing 76.6% of nucleotide identity.

The amino acids sequence encoded by the full spike gene of B15-21 was aligned and compared with human SARS-CoV and SARS-CoV-like Bat CoV sequences ( Figure S1 ). The S2 domain is more conserved than the S1 domain in the spike protein. The receptor binding domain (RBD) of the SARS-CoV-like Bat CoVs, except for RaSHC014, had two major deletion sites, TGNYN (446-450) and PFSPDGKPCTPPA (472-484), compared to human SARS coronaviruses sequences. The receptor binding motif (RBM) and human ACE2 (hACE2)-interacting residue were highly variable between human SARS-CoVs and SARS-CoV-like Bat CoVs but RaSHC014 strain.

Porcine rotavirus-related sequences from pooled bat faeces in case no. 1 were analysed with reference sequences of rotaviruses in human and animals (Fig. 2) . Based on the partial sequences of the VP1, VP3 and VP4 regions, the rotavirus-related sequences were found to belong to rotavirus group H. The group H bat rotavirus in this study 

According to a recent report, there are 23 species of bats in Korea (Han et al., 2012) . Most Korean bats are thought to be insectivores, as no fruit bats have been found in Korea to date, except for some fruit bats as pets or indoor exhibitions (Han et al., 2010 (Chu et al., 2008) . As Bat CoVs B15-40 and 41 were also found in faeces in a cave where Miniopterus schreibersii was dominant, Bat CoV strains in the same lineage with Bat CoV 1A and 1B may be adapted to the Miniopterus spp. It was interesting that Bat CoV B15-21 found in this study was closely related to SARS-CoV-like Bat CoVs, such as Bat CoV Rf1, based on the partial RDRP and S genes, respectively. The newly identified Bat CoV B15-21 clustered with the SARS-CoV-like Bat CoV group in the phylogenetic analysis. The major bat species in the site where Bat CoV B15-21 was collected was Rhinolophus ferrumequinum, which is consistent with a previous report that SARS-CoVlike bat CoVs were mostly found in Rhinolophus spp. (Yuan et al., 2010) . The amino acids of its full spike gene showed that the B15-21 strain is closely related to SARS-CoV-like Bat CoVs Rp3, Rf1, HKU3-2 and 273 rather than human SARS-CoVs and Bat CoV RaSHC014, which was reported to have high potential for human emergence (Menachery et al., 2015) and showed two specific deletion sites in the RBD of the spike protein. One additional betacoronavirus (B15-1-6), which was not detected by the consensus primer-based RDRP-specific RT-PCR (Poon et al., 2005) , was detected by high-throughput sequencing as shown in Table 2 . Three contigs were similar to Bat CoV SC2013, which is a MERS-related betacoronavirus that was found in bats in China (Yang et al., 2014) . One contig was also similar to Bat CoV HKU5-5, which clustered with MERS-CoV (Woo et al., 2007; Annan et al., 2013) . Phylogenetic analysis based on partial spike gene sequences of Bat CoV B15-1-6 showed that it belonged to a MERS-CoV-like clade with previously known MERS-CoV and MERS-CoV-like Bat CoVs.

Although information of the associated bat species was not clearly available, this study first report that SARS-CoV and MERS-CoV-like bat CoVs exist in Korea. This may imply that other strains or types of those coronaviruses may circulate in bat species in Korea. Therefore, we cannot ignore appearance of novel or variant coronavirus through local emergence in Korea, like SARS and MERS in China and Middle East, respectively. In addition to coronaviruses, this study is the first to determine the existence of group H rotaviruses in bat faeces collected in Korea. Group H rotavirus is a recently proposed group of rotaviruses that include strains ADRV-N and B219, which infect human adults (Alam et al., 2007; Matthijnssens et al., 2012) . So far, group H rotaviruses have only been reported in human and pigs (Molinari et al., 2015) , but this study provides evidence that bat species may be a host of group H RVs. To confirm that, there should be follow-up studies including virus isolation and characterization, genomic analysis, continuous surveillance and VP6-based classification (Matthijnssens et al., 2012) to find its prevalence, epidemiology and zoonotic potential.

Localized emergences of zoonotic diseases are usually due to expansion of the wildlife-human interface (Morse et al., 2012) . The human-bat interface in Korea has not been fully described. People usually recognize that bats are flying and roosting in the deep forest or caves. However, the sites where the bat faecal samples were collected showed some evidence of human-bat contact ( Figure S2 ). Some caves or abandoned mine where bats are roosting are located near human habitats. People living there frequently visit into the caves or the abandoned mine to take a rest or to ferment foods. Some caves were also used for some religious ceremony by Korean traditional shamans.

In this study, SARS-CoV-like and MERS-CoV-like bat CoVs and group H rotavirus were detected for this first time in Korea, which may be of interest because of their zoonosis potential. Several kinds of evidence for humanbat contact were also found. Therefore, the potential for novel or variant human viruses from bats in Korea should not be ignored. Continuous surveillance and additional virological research should be a priority to prevent future emergences of zoonotic viral diseases. Intraspecies diversity of SARS-like coronaviruses in Rhinolophus sinicus and its implications for the origin of SARS coronaviruses in humans. J. Gen. Virol. 91, 1058-1062.

Additional Supporting Information may be found in the online version of this article: Table S1 . Information on targets and primers for viral detection. Figure S1 . Amino acid-based multiple alignment of spike proteins from B15-21 and other reference strains. Figure S2 . Possible evidence of human-bat interaction in the bat habitats in this study.

Genetic analysis of an ADRV-N-like novel rotavirus strain B219 detected in a sporadic case of adult diarrhea in Bangladesh

SKA-1(Porcine)-AB576629

J19(Human)-DQ113897

B219(Human)-EF453355

03V0567(Chicken)-NC_021590

CAL-1(Human)-EU490414

Bang373(Human)-NC_021541

03V0568(Chicken)-NC_021625

05V0049(Chicken)-NC_014511

OSU(Porcine)-GU199514

J19(Human)-DQ113900

B219(Human)-EF453357

SKA-1(Porcine)-AB576631

03V0567(Chicken)-NC_021589

Cal-1(Human)-EU490417

Bang-373(Human)-NC_021551

OSU(Porcine)-AY277921

Human betacoronavirus 2c EMC/2012-related viruses in bats

Lyssaviruses and bats: emergence and zoonotic threat

Genomic characterizations of bat coronaviruses (1A, 1B and HKU8) and evidence for co-infections in Miniopterus bats

Reassortant group A rotavirus from straw-colored fruit bat

Detection of influenza A viruses from different species by PCR amplification of conserved sequences in the matrix gene

BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT Nucl

Isolation of a zoonotic pathogen Kluyvera ascorbata from Egyptian fruit-bat Rousettus aegyptiacus

Sounds of the Bats in Korea. National Institute of Biological Resources Press

Bats as reservoirs of severe emerging infectious diseases

Genomic characterization of M and S RNA segments of hantaviruses isolated from bats

Detection of antibodies in Korean bats to Japanese encephalitis virus

VP6-sequence-based cutoff values as a criterion for rotavirus species demarcation

2015: A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence

The metagenomics RAST server -a public resource for the automatic phylogenetic and functional analysis of metagenomes

Genetic variability of VP6, VP7, VP4, and NSP4 genes of porcine rotavirus group H detected in Brazil

2012: Prediction and prevention of the next pandemic zoonosis

Novel paramyxoviruses in bats from Sub-Saharan Africa

Ecological dynamics of emerging bat virus spillover

Identification of a novel coronavirus in bats

MEGA6: molecular evolutionary genetics analysis version 6.0

Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses

RT-PCR for detection of all seven genotypes of Lyssavirus genus

Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features

MERS-related betacoronavirus in Vespertilio superans bats