Bats (order Chiroptera) comprise more than 20% of all living mammalian species; their high diversity and broad geographic distribution distinguishes them from other mammals, after rodents. Bats have been identified as a natural reservoir of several potentially zoonotic viruses, including DNA viruses (such as adenovirus, herpesvirus, and parvovirus) and RNA viruses (such as coronavirus, picornavirus, astrovirus, retrovirus, and reovirus) [1] .
Astroviruses (AstVs) are non-enveloped, positive-sense, single-stranded RNA viruses with a genome ranging from approximately 6.2 to 7.7 kb in size [2] . The Astroviridae family is made up of two genera: Mamastrovirus and Avastrovirus [3] . Mamastrovirus viruses have been identified in numerous mammalian species such as human, pig, mink, sheep, and bat [4, 5] . Avastrovirus viruses have been isolated from avian species including turkey, duck, chicken, fowl, goose, and many others [3, 6] . Astrovirus infections cause gastroenteritis in humans and other mammals [7] . Recently, astrovirus strains associated with encephalitis and ganglionitis were reported in humans, pigs, cattle, and sheep [8] [9] [10] [11] [12] [13] .
Bat astroviruses (BtAstVs) have been discovered in several locations in many countries across Europe, China, South-East Asia, and Africa [14, 15] . In recent studies, diverse coronaviruses and novel rotaviruses were discovered in the feces of Korean bats [16, 17] ; bat astrovirus, however, has not been identified in Korea. Current surveillance studies on bats are insufficient in Korea. Therefore, in this study, we performed a molecular and epidemiological investigation 
A total of 363 bat samples (sixty-one oral swabs, 244 fecal samples, and fifty-eight bat carcasses) were collected at twenty-nine sites across sixteen provinces in Korea from January to September 2016. Bats were captured by net for sample collection, and released immediately after sampling. Oral swabs and feces were transported in viral transport medium (Copan Diagnostics, Murrieta, CA, USA) after being collected from individual bats, and were stored at -70 °C until use. Bat carcasses were autopsied for the collection of internal organs including the brain, lung, intestine, heart, and liver; all collected tissues were immediately ground in phosphate-buffered saline (PBS) containing 1% antibioticantimycotic solution (Corning Inc., Corning, NY, USA). The supernatant of the ground tissues was also stored at -70 °C until use. RNA was extracted within 1 week of the date the samples were collected. 
Phylogenetic analysis was carried out based on the partial 422-bp RdRp sequence of bat astrovirus. Available sequences were retrieved from GenBank. Multiple alignments of astroviral sequences were generated by the Clustal W algorithm in BioEdit v. 7.0.9.0 [19] . Phylogenetic trees were constructed using the neighbor-joining (NJ) method in MEGA 7 [20] . The topologies of the NJ trees were evaluated by the maximum likelihood method and bootstrap values were determined by 1,000 iterations using MEGA 7.
Bat astroviruses were identified from thirty-four samples collected between January and September 2016 ( Fig. 1 ). They were collected from the Sunchang, Yeongcheon, and Danyang regions of Korea. The highest detection rate for astroviral genes was found in Sunchang at 61.5% (8/13 bats), followed by 53.8% in Yeongcheon (7/13 bats). In two areas of Danyang (DY1 and DY4), the detection rate of astrovirus genes was found to be 7.4% (7/94 bats) and 48% (12/25 bats), respectively (Table 1 ). In addition, the detection of astroviral genes by PCR was confirmed in samples collected in January, March, April, and July (Fig. 1 ). In particular, the highest detection rate (50%, 20/40 bats) was found in the samples collected during April, followed by January (25%, 2/8 bats). The detection rates in March and July were 10.5%
(2/19 bats) and 7.2% (10/138 bats), respectively (Fig. 1b) .
Astrovirus was detected consistently in Danyang during January (2/2 bats, 100%), March (2/14 bats, 14.3%), April (12/19 bats, 63.2%), and July (3/13 bats, 23.1%). Additionally, the eight astroviruses detected in Sunchang were only found during April, and the seven astroviruses detected in Yeongcheon were only found during July (Fig. 1a) . Bat astroviruses were discovered in the following four bat species of the thirteen analyzed: Miniopterus fuliginosus, Myotis macrodactylus, M. petax, and Rhinolophus ferrumequinum. The highest positive rate was identified in M. fuliginosus at 25.3% (25/99 bats), followed by M. petax at 25% (4/16 bats). The presence of astroviral RNA sequences in M. macrodactylus and R. ferrumequinum were confirmed at 10% (1/10 bats) and 9.3% (4/43 bats), respectively. The detection rate of bat astrovirus was the greatest in tissue samples (21/57, 36.8%), followed by feces (13/244, 5.3%). No bat astroviruses were detected in oral swabs.
A total of ninety-three astroviral sequences were fully identified from multiple clones isolated from thirty-four bats (accession No. MG840952-MG841043 and MH254918), and newly found astroviruses from Korean bats were tentatively denominated in the order of BtAstV-sample type and number-clone number. The amino acid identity within each clone from the same bat ranged from 50.8-100%. The lowest amino acid identity within each clone from individual bats was in BT34. The amino acid identity between BtAstV-BT34-1 and BtAstV-BT34-2 was 50.8%, and the amino acid identity between BtAstV-BT34-1 and BtAstV-BT34-3 was 55.8%. However, the amino acid identity between BtAstV-BT34-2 and BtAstV-BT34-3 was 75.8% (violet squares in Fig. 2 ). In addition, the amino acid identity of BtAstV-BT42-1 was 57.5% when compared to both BtAstV-BT42-2 and BtAstV-BT42-3 (orange triangles in Fig. 2) .
The amino acid identity among the Korean bat astrovirus sequences ranged from 46.6-100%, and the amino acid identity within single bat species captured from the same habitat also had a wide range of 53.4-100%. The amino acid identity among bat astrovirus sequences from group B in Fig. 2 was 100%, and the amino acid identity of bat astrovirus sequences within group E in Fig. 2 was ≥ 99.1% (Fig. 2) . However, the pairwise amino acid similarity among Korean bat astroviruses from group F (BtAstV-BF73-1, BtAstV-BF73-3, and BtAstV-BT34-1) and D (BtAstV-BT54-1, BtAstV-BF67-2, and BtAstV-BF67-3) in Fig. 2 were the lowest, at 46.6-48.3%. Additionally, the pairwise amino acid similarity among the astroviruses from group F and B in Fig. 2 were confirmed to be 47.5-51.6%. In group C in Fig. 2 , the bat astrovirus detected from R. ferrumequinum (BtAstV-BF67-1) and the bat astrovirus detected from M. fuliginosus (BtAstV-BT53-1) were closely clustered, with an amino acid identity of ≥ 99.1%. Furthermore, in group A in Fig. 2 , bat astroviruses identified from M. macrodactylus (BtAstV-BT40-1, BtAstV-BT40-2, and BtAstV-BT40-3) were clustered with the astroviruses confirmed from M. fuliginosus (BtAtV-BT53-2 and BtAstV-BF78-2); the average amino acid identity among the bat astroviruses in group A in Fig. 2 was ≥ 94.1%.
The detection of partial ORF2 (capsid protein) sequences from several astroviral RdRp positive samples-BT6 and BT27 of R. ferrumquinum, BT51 of M. petax, and BF69 and BF80 of M. fuliginosusi-was performed. Two partial ORF2 (632 bp) sequences (BT6 and BT27) were obtained and the phylogenetic distance was calculated. The ORF2 (Fig. 3 ).
In the present study, we identified the infection of four Korean bat species-M. fuliginosus, M. macrodactylus, M. petax, and R. ferrumequinum-with various astroviruses. These bats are distributed across a wide geographical area that includes Korea, Japan, and Russia [21] [22] [23] [24] . The discovery of astroviruses from these bat species has been reported in China and Hungary, but not Korea [25] [26] [27] . Bat astroviruses are known to have an exceedingly high genetic diversity. In a previous study, the amino acid similarity of the RdRp gene within bat astroviruses in Asia and Europe was verified to be ~49.2% [26, 28] . In this study, the amino acid identity of the partial RdRp gene in Korean bat astroviruses ranged from 46.6-100%, and the amino acid identity of astroviruses within a single bat species captured from the same habitat also ranged from 53.4-100%. Moreover, the amino acid identity of multiple astroviral strains from within an individual bat was 50.8-100%. This amino acid identity range within multiple clones of bat astroviruses indicated the existence of astrovirus haplotypes in Korean bats; the very low amino acid identity may indicate the possibility of infection by two or more astrovirus species in individual bats, or hyperpolymorphisms within identical astrovirus strains; however, more systematic studies and more genomic sequences are needed to verify this. Additionally, astroviruses are known to widely spread among Chiroptera and show no host restriction [15, 18, 25, 26] . Our result showed that astrovirus sequences in multiple clones from one M. fuliginosus bat were closely clustered with bat astroviruses from different host species: M. macrodactylus and R. ferrumequinum. In group D, the clustering between M. fuliginosus (BtAstV-BT54-1) and R. ferrumequinum (BtAstV-BF67-2) showed an amino acid similarity of ≥ 99%; this was supported by high posterior probabilities (pp = 0.99) from ▸ phylogenetic topology. These phylogenetic results suggest that genetic variants were generated at high rates within individual bats, and that active transmission events between bat species or families occurred between the Vespertilionidae and Rhinolophidae families. Furthermore, of the available astrovirus sequences, the partial amino acid sequence of ORF2 from two bat astroviruses was confirmed to have the highest identity (73.3-74.8%) with the bat astrovirus LC03 (FJ571074) and LS11 (FJ571067) sequences from China. It is unclear whether these represent novel astrovirus species according to ICTV guidelines [29] , as we could not characterize the entire ORF2 sequence. Therefore, further studies on the characterization of the whole genomes of bat astroviruses is necessary, as Korean bat astroviruses have not been comprehensively analyzed until now.
Various emerging viruses, such as nipah virus, filovirus, coronavirus, and astrovirus in bats, have been known to display seasonal variations in concentration or prevalence [28, 30, 31] . For astrovirus persistence, recent studies have shown that dry roosting at the beginning of the rainy season may increase the rate of transmission [32] . Accordingly, a high concentration and high prevalence of astrovirus was found in the samples collected during spring and autumn, indicating that they were most susceptible in the periods before and after parturition [28] . Our results showed that the prevalence of astrovirus in Korean bats was at its maximum during spring, the period before parturition, and before the rainy season. Knowing the seasonal fluctuations of bat astrovirus infections is required to understand the transmission and dispersal of astroviruses within bat populations.
Korean bat astroviruses were distributed across three regions: Sunchang, Danyang, and Yeongcheon. A recent study showed that no direct transmission occurred across different geographical locations through the analysis of astrovirus haplotypes in individual bats [33] ; however, our study identified different phylogenetic topologies. The indiscriminate phylogenetic clustering of astroviral sequences obtained from different locations implies the direct (or indirect) transmission of bat astroviruses with no geographical restrictions in Korea. Additionally, astrovirus was highly prevalent in feces because astroviruses replicate in the gastrointestinal tract, and transmission has been found to occur through the oral-fecal route [34] . Interestingly, Korean bat astroviruses only showed high detection rates in tissue (including colon) and fecal samples, but not in oral swabs.
In bats, the co-infection and recombination of viruses has been occasionally reported [18, 35] . This suggests that the ecological characteristics of bats may increase the opportunities for recombination in various viruses. The movement of bats due to a decrease in the size of their habitats causes an increase in contact opportunities with domestic animals as well as humans; zoonotic pathogens in bats may thus pose a potential risk of transmission to humans and domestic animals. Therefore, further surveillance studies of bat viruses are required due to the lack of clear knowledge on bat pathogens in South Korea.