Long-term surveillance of bat coronaviruses in Korea: Diversity and distribution pattern

Circulation characteristics of bat coronaviruses linked to bat ecological factors in Korea, 2021-2022

Considering that bat ecology alterations may be linked with pathogen spillover, research on bat coronaviruses, particularly on the infection and transmission pattern among bats in relation with their ecology, is essential. We captured bats distributed in Korea from 2021 to 2022, examined coronaviruses in oral swabs, feces, urine, and ectoparasites, and were able to detect alphacoronavirus. We investigated coronaviruses, but noted no substantial differences in the body condition index in the coronavirus-positive bats. Binary logistic regression analysis revealed that bat ecological factors that were significantly associated with coronavirus-positive were roost type, sample type, and bat species. Coronavirus-positive ectoparasite cases suggested additional study on the potential role of them as the viral transmission vectors or fomites. Reinfection of a different coronavirus in recaptured bats was evident, suggesting the possibility that coronavirus circulation can evade the potential protective immunity acquired from previous coronavirus infections. The present findings provide comprehensive information on the coronaviruses transmission dynamics within bat populations linked with bat ecology.

Differential prevalence and risk factors for infection with coronaviruses in bats collected from Yunnan Province, China

Coronaviruses (CoVs) pose a threat to human health globally, as highlighted by severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS) and the COVID-19 pandemic. Bats from the Greater Mekong Subregion (GMS) are an important natural reservoir for CoVs. Here we report the differential prevalence of CoVs in bats within Yunnan Province across biological and ecological variables. We also show the coexistence of CoVs in individual bats and identify an additional putative host for SARS-related CoV, with higher dispersal capacity than other known hosts. Notably, 11 SARS-related coronaviruses (SARSr-CoVs) were discovered in horseshoe bats (family Rhinolophidae) and a Chinese water myotis bat (Myotis laniger) by pan-CoV detection and Illumina sequencing. Our findings facilitate an understanding of the fundamental features of the distribution and circulation of CoVs in nature as well as zoonotic spillover risk in the One health framework.

Detection of coronaviruses in insectivorous bats of Fore-Caucasus, 2021

Coronaviruses (CoVs) pose a huge threat to public health as emerging viruses. Bat-borne CoVs are especially unpredictable in their evolution due to some unique features of bat physiology boosting the rate of mutations in CoVs, which is already high by itself compared to other viruses. Among bats, a meta-analysis of overall CoVs epizootiology identified a nucleic acid observed prevalence of 9.8% (95% CI 8.7-10.9%). The main objectives of our study were to conduct a qPCR screening of CoVs' prevalence in the insectivorous bat population of Fore-Caucasus and perform their characterization based on the metagenomic NGS of samples with detected CoV RNA. According to the qPCR screening, CoV RNA was detected in 5 samples, resulting in a 3.33% (95% CI 1.1-7.6%) prevalence of CoVs in bats from these studied locations. BetaCoVs reads were identified in raw metagenomic NGS data, however, detailed characterization was not possible due to relatively low RNA concentration in samples. Our results correspond to other studies, although a lower prevalence in qPCR studies was observed compared to other regions and countries. Further studies should require deeper metagenomic NGS investigation, as a supplementary method, which will allow detailed CoV characterization.

Opportunities and Limitations of Molecular Methods for Studying Bat-Associated Pathogens

Bats have been identified as reservoirs of zoonotic and potentially zoonotic pathogens. Significant progress was made in the field of molecular biology with regard to infectious diseases, especially those that infect more than one species. Molecular methods, sequencing and bioinformatics have recently become irreplaceable tools in emerging infectious diseases research and even outbreak prediction. Modern methods in the molecular biology field have shed more light on the unique relationship between bats and viruses. Here we provide readers with a concise summary of the potential and limitations of molecular methods for studying the ecology of bats and bat-related pathogens and microorganisms.

Characterization of replication and variations in genome segments of a bat reovirus, BatMRV/B19-02, by RNA-seq in infected Vero-E6 cells

Mammalian orthoreoviruses (MEVs) that can cause enteric, respiratory, and encephalitic infections have been identified in a wide variety of mammalian species. Here, we report a novel MRV type 1 strain detected in Miniopterus schreibersii that may have resulted from reassortment events. Using next-generation RNA sequencing (RNA-seq), we found that the ratios of the RNA levels of the 10 reovirus segments in infected cells were constant during the late stages of infection. We also discovered that the relative abundance of each segment differed. Notably, the relative abundance of M2 (encoding the µ1 protein) and S4 (encoding the σ3 protein) RNAs was higher than that of the others throughout the infection. Additionally, massive junctions were identified. These results support the hypothesis that defective genome segments are generated and that cross-family recombination occurs. These data may further the study of gene function, viral replication, and virus evolution.

Testing bats in rehabilitation for SARS-CoV-2 before release into the wild

Several studies have suggested SARS-CoV-2 originated from a viral ancestor in bats, but whether transmission occurred directly or via an intermediary host to humans remains unknown. Concerns of spillover of SARS-CoV-2 into wild bat populations are hindering bat rehabilitation and conservation efforts in the United Kingdom and elsewhere. Current protocols state that animals cared for by individuals who have tested positive for SARS-CoV-2 cannot be released into the wild and must be isolated to reduce the risk of transmission to wild populations. Here, we propose a reverse transcription-quantitative polymerase chain reaction (RT-qPCR)-based protocol for detection of SARS-CoV-2 in bats, using fecal sampling. Bats from the United Kingdom were tested following suspected exposure to SARS-CoV-2 and tested negative for the virus. With current UK and international legislation, the identification of SARS-CoV-2 infection in wild animals is becoming increasingly important, and protocols such as the one developed here will help improve understanding and mitigation of SARS-CoV-2 in the future.

Genomic Comparisons of Alphacoronaviruses and Betacoronaviruses from Korean Bats

Coronaviruses are well known as a diverse family of viruses that affect a wide range of hosts. Since the outbreak of severe acute respiratory syndrome, a variety of bat-associated coronaviruses have been identified in many countries. However, they do not represent all the specific geographic locations of their hosts. In this study, full-length genomes representing newly identified bat coronaviruses in South Korea were obtained using an RNA sequencing approach. The analysis, based on genome structure, conserved replicase domains, spike gene, and nucleocapsid genes revealed that bat Alphacoronaviruses are from three different viral species. Among them, the newly identified B20-97 strain may represent a new putative species, closely related to PEDV. In addition, the newly-identified MERS-related coronavirus exhibited shared genomic nucleotide identities of less than 76.4% with other Merbecoviruses. Recombination analysis and multiple alignments of spike and RBD amino acid sequences suggested that this strain underwent recombination events and could possibly use hDPP4 molecules as its receptor. The bat SARS-related CoV B20-50 is unlikely to be able to use hACE2 as its receptor and lack of an open reading frame in ORF8 gene region. Our results illustrate the diversity of coronaviruses in Korean bats and their evolutionary relationships. The evolution of the bat coronaviruses related ORF8 accessory gene is also discussed.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses. Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eight following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin sarbecoviruses, (iv) three bat sarbecoviruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat sarbecoviruses from Yunnan, (vi) the bat sarbecovirus from Japan, (vii) the bat sarbecovirus from Bulgaria, and (viii) the bat sarbecovirus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes have less cytosines and more uracils at third codon-positions than other sarbecoviruses, whereas the genomes of pangolin sarbecoviruses show more adenines at third codon-positions. I suggest that taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the eight groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats and their latitudinal distribution is also discussed. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat sarbecovirus.

Using Environmental Sampling to Enable Zoonotic Pandemic Preparedness

The current pandemic caused by the SARS CoV-2, tracing back its origin possibly to a coronavirus associated with bats, has ignited renewed interest in understanding zoonotic spillovers across the globe. While research is more directed towards solving the problem at hand by finding therapeutic strategies and novel vaccine techniques, it is important to address the environmental drivers of pathogen spillover and the complex biotic and abiotic drivers of zoonoses. The availability of cutting-edge genomic technologies has contributed enormously to preempt viral emergence from wildlife. However, there is still a dearth of studies from species-rich South Asian countries, especially from India. In this review, we outline the importance of studying disease dynamics through environmental sampling from wildlife in India and how ecological parameters of both the virus and the host community may play a role in mediating cross-species spillovers. Non-invasive sampling using feces, urine, shed hair, saliva, shed skin, and feathers has been instrumental in providing genetic information for both the host and their associated pathogens. Here, we discuss the advances made in environmental sampling protocols and strategies to generate genetic data from such samples towards the surveillance and characterization of potentially zoonotic pathogens. We primarily focus on bat-borne or small mammal-borne zoonoses and propose a conceptual framework for non-invasive strategies to tackle the threat of emerging zoonotic infections.

Genomic Characterization of a Novel Alphacoronavirus Isolated from Bats, Korea, 2020

Coronavirus, an important zoonotic disease, raises concerns of future pandemics. The bat is considered a source of noticeable viruses resulting in human and livestock infections, especially the coronavirus. Therefore, surveillance and genetic analysis of coronaviruses in bats are essential in order to prevent the risk of future diseases. In this study, the genome of HCQD-2020, a novel alphacoronavirus detected in a bat (Eptesicus serotinus), was assembled and described using next-generation sequencing and bioinformatics analysis. The comparison of the whole-genome sequence and the conserved amino acid sequence of replicated proteins revealed that the new strain was distantly related with other known species in the Alphacoronavirus genus. Phylogenetic construction indicated that this strain formed a separated branch with other species, suggesting a new species of Alphacoronavirus. Additionally, in silico prediction also revealed the risk of cross-species infection of this strain, especially in the order Artiodactyla. In summary, this study provided the genetic characteristics of a possible new species belonging to Alphacoronavirus.

Genetic diversity of bat coronaviruses and comparative genetic analysis of MERS-related coronaviruses in South Korea

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.

Dispersal history of Miniopterus fuliginosus bats and their associated viruses in east Asia

In this study, we examined the role of the eastern bent-winged bat (Miniopterus fuliginosus) in the dispersion of bat adenovirus and bat alphacoronavirus in east Asia, considering their gene flows and divergence times (based on deep-sequencing data), using bat fecal guano samples. Bats in China moved to Jeju Island and/or Taiwan in the last 20,000 years via the Korean Peninsula and/or Japan. The phylogenies of host mitochondrial D-loop DNA was not significantly congruent with those of bat adenovirus (m2XY = 0.07, p = 0.08), and bat alphacoronavirus (m2XY = 0.48, p = 0.20). We estimate that the first divergence time of bats carrying bat adenovirus in five caves studied (designated as K1, K2, JJ, N2, and F3) occurred approximately 3.17 million years ago. In contrast, the first divergence time of bat adenovirus among bats in the 5 caves was estimated to be approximately 224.32 years ago. The first divergence time of bats in caves CH, JJ, WY, N2, F1, F2, and F3 harboring bat alphacoronavirus was estimated to be 1.59 million years ago. The first divergence time of bat alphacoronavirus among the 7 caves was estimated to be approximately 2,596.92 years ago. The origin of bat adenovirus remains unclear, whereas our findings suggest that bat alphacoronavirus originated in Japan. Surprisingly, bat adenovirus and bat alphacoronavirus appeared to diverge substantially over the last 100 years, even though our gene-flow data indicate that the eastern bent-winged bat serves as an important natural reservoir of both viruses.

Long-term surveillance of bat coronaviruses in Korea: Diversity and distribution pattern

Bats harbour diverse coronaviruses (CoVs), some of which are associated with zoonotic infections, as well as inter‐species transmission. In this study, a total of 512 bat faecal samples from the bat habitats at different geographical locations in South Korea were investigated between 2016 and 2019. Seventy‐eight samples were positive for coronaviruses (15.2%), comprising 68 alphacoronaviruses (13.3%) and 10 betacoronaviruses (2.0%). The positive rates tended to increase during the awakening (April) period. Notably, betacoronaviruses were only found in the site where Rhinolophus ferrumequinum was the major species of bats, and were related to SARS‐ and MERS‐related CoVs identified in China and South Korea, respectively. No betacoronaviruses were closely related to SARS‐CoV‐2 in this study. Alphacoronaviruses were detected in the sites where Hypsugo alaschanicus, Miniopterus fuliginosus, Miniopterus schreibersii, Rhinolophus ferrumequinum, Myotis bombinus, Myotis macrodactylus and Myotis petax were found to be the major bat species. Furthermore, alphacoronaviruses had higher genetic diversity than betacoronaviruses and had a wider distribution in Korea. Considering that different bat species are co‐roosting in crowded conditions in the same habitat, the diverse coronaviruses in Korean bats are likely to undergo cross‐species transmission events due to the richness in host species. Therefore, continuous monitoring should be performed, especially at the awakening time of the hibernating bats in the habitats where diverse bat species co‐roost, to better understand the evolution of coronaviruses in bats.