Bats are natural reservoirs for potential zoonotic viruses. In this study, next-generation sequencing was performed to obtain entire genome sequences of picornavirus from a picornavirus-positive bat feces sample (16BF77) and to explore novel viruses in a pooled bat sample (16BP) from samples collected in South Korea, 2016. Fourteen mammalian viral sequences were identified from 16BF77 and 29 from 16BP, and verified by RT-PCR. The most abundant virus in 16BF77 was picornavirus. Highly variable picornavirus sequences encoding 3Dpol were classified into genera Kobuvirus, Shanbavirus, and an unassigned group within the family Picornaviridae. Amino acid differences between these partial 3Dpol sequences were ≥ 65.7%. Results showed that one bat was co-infected by picornaviruses of more than two genera. Retrovirus, coronavirus, and rotavirus A sequences also were found in the BP sample. The retrovirus and coronavirus genomes were identified in nine and eight bats, respectively. Korean bat retroviruses and coronavirus demonstrated strong genetic relationships with a Chinese bat retrovirus (RfRV) and coronavirus (HKU5-1), respectively. A co-infection was identified in one bat with a retrovirus and a coronavirus. Our results indicate that Korean bats were multiply infected by several mammal viruses.

Bats have been identified as a natural reservoir of several potentially zoonotic viruses, including Lyssavirus, Ebola virus, Marburg virus, Hendra virus, Nipah virus, as well as severe acute respiratory syndrome and Middle East respiratory syndrome coronavirus (CoV). Here, we performed a molecular epidemiological investigation of South Korean bat viruses. Genetic comparative analysis was performed on the spike glycoprotein gene of the detected MERS-related CoVs. Among 1640 samples (348 oral swabs, 1199 faecal samples, 83 urine samples and 10 bat carcass) collected across 24 South Korean provinces during 2017-2019, CoV was detected in 82 samples (75 faeces and seven oral swab samples) from 11 provinces. Surveillance over the 3 years during which samples were collected revealed significantly higher CoV detection rates between spring and autumn, and a high detection rate in Vespertillionidae and Rhinolophidae bats. Our phylogenetic analysis shows that Korean bat CoVs are genetically diverse regardless of their spatiotemporal distribution and their host species, and that the discovered bat CoVs belong to various subgenera within the Alpha-and Betacoronavirus genera. Twenty detected MERS-related CoVs belonging to the genus Betacoronavirus were similar to the Ia io bat CoV NL140422 and NL13845 strains. A comprehensive genetic analysis of two Korean bat MERS-related CoV spike receptor binding domain (RBDs)(176 and 267 strains) showed that the 18 critical residues that are involved in interactions with the human DPP4 receptor are most similar to the NL13845 strain, which is known to not bind with hDPP4. A deeper analysis of the interfacing residues in the Korean bat MERS-related CoVs RBD-hDPP4 complexes showed that the Korean bat CoVs has fewer polar contacts than the NL13845 strain. Although further study will be needed, these results suggest that Korean bat MERS-related CoVs are unlikely to bind with hDPP4. Nevertheless, these findings highlight the need for continuous monitoring to identifying the origin of new infectious diseases, specifically mutant CoV.