Genotype and Phenotype Characterization of
Rhinolophus sp. Sarbecoviruses from Vietnam:
Implications for Coronavirus Emergence

Sarah Temmam, Tran Cong Tu, Béatrice Regnault, Massimiliano Bonomi,
Delphine Chrétien, Léa Vendramini, Tran Nhu Duong, Tran Vu Phong,
Nguyen Thi Yen, Hoang Ngoc Anh, et al.

► To cite this version:

Sarah Temmam, Tran Cong Tu, Béatrice Regnault, Massimiliano Bonomi, Delphine Chrétien, et al. Genotype and Phenotype Characterization of Rhinolophus sp. Sarbecoviruses from Vietnam: Implications for Coronavirus Emergence. Viruses, 2023, 15 (9), pp.1897. 10.3390/v15091897. pasteur-04228724

HAL Id: pasteur-04228724

https://pasteur.hal.science/pasteur-04228724

Submitted on 4 Oct 2023

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Genotype and Phenotype Characterization of *Rhinolophus* sp. Sarbecoviruses from Vietnam: Implications for Coronavirus Emergence

Sarah Temmam<sup>1,2,†</sup> ID, Tran Cong Tu<sup>3,†</sup>, Béatrice Regnault<sup>1,2</sup>, Massimiliano Bonomi<sup>4</sup>, Delphine Chrétien<sup>1,2</sup>, Léa Vendramini<sup>1,2</sup>, Tran Nhu Duong<sup>3</sup>, Tran Vu Phong<sup>3</sup>, Nguyen Thi Yen<sup>3</sup>, Hoang Ngoc Anh<sup>3</sup>, Tran Hai Son<sup>3 ID</sup>, Pham Tuan Anh<sup>3</sup>, Faustine Amara<sup>5</sup>, Thomas Bigot<sup>1,6 ID</sup>, Sandie Munier<sup>5</sup>, Vu Dinh Thong<sup>7</sup>, Sylvie van der Werf<sup>8,9</sup>, Vu Sinh Nam<sup>3</sup> and Marc Eloit<sup>1,2,10,* ID</sup>

1. Pathogen Discovery Laboratory, Institut Pasteur, Université Paris Cité, 75015 Paris, France

2. Institut Pasteur, The OIE Collaborating Center for the Detection and Identification in Humans of Emerging Animal Pathogens, Université Paris Cité, 75015 Paris, France

3. National Institute of Hygiene and Epidemiology, Hanoi 100000, Vietnam

4. Structural Bioinformatics Unit, Institut Pasteur, CNRS UMR3528, Université Paris Cité, 75015 Paris, France

5. Institut Pasteur, G5 Evolutionary Genomics of RNA Viruses, Université Paris Cité, 75015 Paris, France

6. Institut Pasteur, Bioinformatics and Biostatistics Hub, Université Paris Cité, 75015 Paris, France

7. Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology (VAST), Hanoi 70072, Vietnam

8. Molecular Genetics of RNA Viruses Unit, Institut Pasteur, CNRS UMR 3569, Université Paris Cité, 75015 Paris, France

9. Institut Pasteur, National Reference Center for Respiratory Viruses, Université Paris Cité, 75015 Paris, France

10. Ecole Nationale Vétérinaire d'Alfort, University of Paris-Est, 77420 Maisons-Alfort, France

* Correspondence: marc.eloit@pasteur.fr

† These authors contributed equally to this work.

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Temmam, S.; Tu, T.C.; Regnault, B.; Bonomi, M.; Chrétien, D.; Vendramini, L.; Duong, T.N.; Phong, T.V.; Yen, N.T.; Anh, H.N.; et al. Genotype and Phenotype Characterization of *Rhinolophus* sp. Sarbecoviruses from Vietnam: Implications for Coronavirus Emergence. *Viruses* **2023**, *15*, 1897. https://doi.org/10.3390/v15091897

Academic Editor: Norbert Nowotny

Received: 13 July 2023

Revised: 11 August 2023

Accepted: 6 September 2023

Published: 8 September 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Bats are a major reservoir of zoonotic viruses, including coronaviruses. Since the emergence of SARS-CoV in 2002/2003 in Asia, important efforts have been made to describe the diversity of *Coronaviridae* circulating in bats worldwide, leading to the discovery of the precursors of epidemic and pandemic sarbecoviruses in horseshoe bats. We investigated the viral communities infecting horseshoe bats living in Northern Vietnam, and report here the first identification of sarbecoviruses in *Rhinolophus thomasi* and *Rhinolophus siamensis* bats. Phylogenetic characterization of seven strains of Vietnamese sarbecoviruses identified at least three clusters of viruses. Recombination and cross-species transmission between bats seemed to constitute major drivers of virus evolution. Vietnamese sarbecoviruses were mainly enteric, therefore constituting a risk of spillover for guano collectors or people visiting caves. To evaluate the zoonotic potential of these viruses, we analyzed in silico and in vitro the ability of their RBDs to bind to mammalian ACE2s and concluded that these viruses are likely restricted to their bat hosts. The workflow applied here to characterize the spillover potential of novel sarbecoviruses is of major interest for each time a new virus is discovered, in order to concentrate surveillance efforts on high-risk interfaces.

**Keywords:** Sarbecovirus; horseshoe bat; virus evolution; spillover

## 1. Introduction

Bats constitute a reservoir of many viruses, among which several have the potential to cross the species barrier from their natural hosts and to infect humans and domestic animals, directly or via an intermediate host [1]. In particular, bats are a reservoir of three of the ten virus groups that present pandemic potential, as designated by the World Health Organization: henipaviruses, filoviruses and coronaviruses [2]. Indeed, during recent decades, several coronaviruses emerged from their bat hosts and were responsible for multiple and
severe human and animal epidemics or pandemics. These include severe acute respiratory syndrome (SARS) [3], Middle East respiratory syndrome (MERS) [4] and the coronavirus disease (COVID-19) in humans; and swine acute diarrhea syndrome (SADS) [5] and porcine epidemic diarrhea virus (PEDV) in pigs [6]. More than 4000 coronaviruses from 14 different bat families have been identified so far [7], a number that is probably underestimated given the richness of bat species distributed worldwide and the recurrent discovery of novel bat-associated coronaviruses [8–14]. Since their distribution and spillover potential are usually unknown, some of these viruses may present public and animal health concerns, and a global health approach is needed to better evaluate the risk of emergence of such coronaviruses [15].

In particular, the recent emergence of SARS-CoV-2, which led to the COVID-19 pandemic, highlighted the importance of monitoring the circulation of bat-associated sarbecoviruses worldwide. In fact, serological evidence of prior exposure of humans to bat sarbecoviruses in rural China suggested that wildlife-to-human spillovers from bats may occur frequently but silently [16–18]. This evidence relies mainly on antigen/antibody binding methods with questionable specificity, especially when the viral nucleoprotein is targeted. By contrast, very specific methods such as seroneutralization did not reveal significant cross-species transmissions of bat sarbecoviruses to humans [19]. However, we cannot exclude that such events sometimes occur and could represent opportunities for virus adaptation that may lead to better human-to-human transmission. The necessity of increasing surveillance of bat sarbecoviruses has therefore emerged. Most novel bat-related sarbecoviruses have been associated with horseshoe bats from Asia, and to a lesser extent from Europe [20–22] and Africa [23,24]. For example, surveys of bats in China have revealed high diversity and prevalence (5–10%) of sarbecoviruses in different rhinolophid bat species [25,26]. Bat precursors of SARS-CoV-2 have been identified in *Rhinolophus* sp. bats from Laos [27], Thailand [28], Cambodia [29], South Korea [30] and Japan [31]. For the majority of these novel viruses, it is still unknown if they are able to infect various hosts (including humans) apart from their bat reservoir.

Sarbecoviruses use the cellular angiotensin-converting enzyme 2 (ACE2) receptor to enter and infect mammalian cells [32,33]. The spike viral sequence, and especially the receptor-binding domain (RBD), are major determinants of host range and tropism [34]. ACE2 is a ubiquitous receptor present in all vertebrate animals, and several studies demonstrated that it is an efficient *Sarbecovirus* receptor in both bat reservoirs and susceptible mammalian hosts [35], such as palm civets [36], pangolins [37] and mustelids [38]. Several bat-associated sarbecoviruses were shown to use the ACE2 receptor to enter mammalian cells [39–41], while most of them presented deletions in their RBDs, precluding their use of ACE2 as a cellular receptor. However, some bat-associated sarbecoviruses were shown to infect human and bat cells through an ACE2-independent entry [42]. Recently, we demonstrated that the bat-associated *Rh. marshalli* BANAL-236 sarbecovirus, phylogenetically close to SARS-CoV-2, was able to use ACE2 to infect mammalian cells but presented an enteric tropism in primates contrary to its SARS-CoV-2 relative, possibly due to its lack of a furin cleavage site [19,27]. Therefore, several factors other than the ability to bind to ACE2 may influence host tropism and the restriction of certain sarbecoviruses to specific host species [43–45].

In Vietnam, the prevalence of coronaviruses circulating among bats is high and was estimated to be between 22% [46] and 75% [47] depending on the species and location of sampling. In a previous study conducted in Southern Vietnam, thirty-four *Pedacovirus* sequences close to Scotophilus bat CoV-512 were identified in *Scotophilus kuhlii* bats roosting in bat guano farms [48]. A more recent study, also implemented around bat guano collection sites and pig farming sites, identified progenitors of pig-infecting *Pedacovirus* in bats roosting at the human–wildlife–domestic animal (i.e., pig) interface [49]. Several *Merbecovirus*-, *Tegacovirus*- and *Nobecovirus*-related viruses were also identified, but no *Sarbecovirus*. Here, we report the first identification and characterization of *Rhinolophus* sp. sarbecoviruses from Northern Vietnam. We performed molecular and genetic character-
ization of these viruses through broad-range phylogenetic analyses, and evaluated their phenotype regarding their spillover potential to humans and wild and domestic animals by testing in silico and in vitro the ability of their RBDs to bind to mammalian ACE2s, a major determinant of host tropism. We conclude that these viruses are likely restricted to their bat hosts and are prone to recombinations and interspecies transmissions than can favor the emergence of viruses with an extended host spectrum or tropism.

## 2. Materials and Methods

### 2.1. Ethical and Regulatory Statements

This study, including bat sampling and handling, was approved by the Ministry of Health (Vietnam) on the 21 September 2020 (reference 4026/QD-BYT). The detailed protocol was submitted to the Institutional Review Board (IRB) of the National Institute of Hygiene and Epidemiology (NIHE), which stated that the study was exempted from IRB request, on the 10 September 2020. All animals were captured, handled, sampled and released following previously published protocols and ASM guidelines [50]. All participants were protected by personal protective equipment that comprised protective clothing with long sleeves and covering the head, neck, arms, shoulders and upper body parts. Exportation from Vietnam and importation to France were conducted according to national regulations.

### 2.2. Bat Capture and Sample Collection

Two bat surveys were conducted from 22 February to 3 March and from 4 July to 13 July 2021 in 9 caves from two districts of Sơn La Province, Vietnam: Sốp Cộp from 22 to 26 February and Mộc Châu from 27 February to 3 March and from 9 to 13 July; and Sơn La City from 4 to 8 July (Figure 1, Table S1).

<figure><img src="image_4.png" /></figure>

**Figure 1.** Bat sampling in Northern Vietnam, 2021. Maps were constructed with Quantum GIS software (https://qgis.org/en/site/ accessed on 6 September 2023). Sampling performed during February–March 2021 is presented by a black dot (Mission 1), and sampling performed during July 2021 is presented by a red dot (Mission 2). The size of the pie charts is proportional to the number of bats sampled in each location. Bats were captured using mist nets and hand nets. Each captured bat was weighted using a Pesola spring scale, then the bats were individually kept in porous humidified cotton bags and placed in a cool dry place for sampling within six hours. The following basic set of samples was collected from each bat, and included saliva (oropharyngeal swab), feces (fresh fecal sample) or rectal swab, blood (serum; rbc/wbc pellet) and urine (free catch method or urogenital swab). After sampling, bats were released at the collection site.

Bats were morphologically identified as follows: photographs and selected external measurements of every captured bat were taken for identification based on morphological diagnoses in comparison with published descriptions and confirmed features of *Rhinolophus thomasi* from Cuc Phuong National Park, Northern Vietnam. Morphologically, *Rhinolophus*
*thomasi* is distinguished from other horseshoe bat species in Vietnam and surrounding countries by a complex of external characteristics including a short lancet with a rudimentary tip; a rounded connecting process; horseshoe does not cover the muzzle completely and has well-developed lateral leaflets; and a lower lip that has three mental grooves [51,52].

To confirm the morphological identifications, we took advantage of the concomitant sequencing of bat transcriptome and used the Barcode of Life Data Systems (BOLD) as previously described [53]. Briefly, all trimmed reads were mapped onto the Cytochrome Oxydase subunit I (COI) database, followed by de novo assembly of mapped reads. The resulting contigs were uploaded to the BOLD Identification System (https://www.boldsystems.org/index.php/IDS_OpenIdEngine) (accessed on 7 August 2023) for species identification. When a discrepancy between BOLD-based and morphology-based species identification was observed, we determined the most probable species by positioning the COI sequence of the bat individual into a phylogenetic tree based on 103 sequences of the N-terminal region of the Cytochrome Oxidase I gene (region of 657 bp), representative of 35 *Rhinolophus* species.

## 2.3. RNA Extraction and Pan-Coronavirus PCR Screening

Total RNA was extracted from 100 µL of the samples using the NucleoSpin® 8 Virus kit (Macherey-Nagel, Düren, North Rhine-Westphalia, Germany) according to the manufacturer's instructions. Extracted RNA was stored at -80 °C until further analysis. Reverse transcription of RNA extracts was performed with the SuperScript IV RT kit (Invitrogen, Carlsbad, CA, USA) using random hexamers and according to the manufacturer's recommendations. Samples were screened for the presence of coronaviruses with the pan-coronavirus nested PCR system targeting the viral polymerase, as previously described [54], and confirmed via Sanger sequencing.

## 2.4. Whole Genome Sequencing of Bat Sarbecoviruses

Complete *Sarbecovirus* genomes were obtained via *Betacoronavirus*-enriched next-generation sequencing, as previously described [27]. Briefly, specific *Betacoronavirus* primers were mixed with random hexamers and used in a reverse transcription reaction, followed by a MALBAC random amplification. NGS libraries were prepared for each individual positive bat sample using the DNA Prep (M) Tagmentation® kit (Illumina), and sequenced onto NextSeq500 or NextSeq2000 Illumina sequencers, at an average depth of 80 million 150 bp or 100 bp single reads per sample. When needed, conventional PCR and Sanger sequencing were used to fill the gaps in the genomes (primers are available on request).

## 2.5. Molecular and Phylogenetic Analyses

A comprehensive database of complete *Sarbecovirus* genomes was first built based on two recent studies [23,55] and completed with other sequences manually picked from GenBank and GISAID (Table S2) to cover non-human (bat, pangolin, palm civet and ferret badger) sarbecovirus diversity. A total of 221 complete genomes were selected. Several human prototype strains of SARS-CoV and SARS-CoV-2 were also included in the final dataset. For each complete genome, 5'-3' UTRs and intergenic regions were manually removed and genes were concatenated before being translated with C2A.A2C [56]. Protein sequences were then aligned with MAFFT (in auto mode, global alignment [57]) and backtranslated in nucleotides using the original sequences. Given the high variability in *Sarbecovirus* sequences, this method ensured the correct alignment of selected whole genomes. The same strategy was applied to align Vietnamese *Sarbecovirus* genomes.

Phylogenetic analyses were performed at the gene level, either for Vietnamese isolates or for the complete *Sarbecovirus* dataset, with W-IQ-TREE [58], using the maximum likelihood reconstruction method and the aBayes branch support parameter, after determination of the best substitution model fitting the data implemented in W-IQ-TREE. Trees were visualized with FigTree [59].
Identity matrices (expressed in %) were calculated at the whole genome and at each gene nucleotide and amino-acid level using CLC Main Workbench 21.0.4 (Qiagen, Hong Kong, China).

Recombination analyses were conducted on the Vietnamese or on the complete *Sarbecovirus* datasets using the "full exploratory recombination scan" implemented in the Recombination Detection Program (RDP v.5.36) [60], which includes RDP [61], GENECONV [62], MaxChi [63], Chimaera [64] and 3Seq [65] tools. Only recombination events detected with at least three of the above cited methods were considered true recombination events.

## 2.6. *In Silico Evaluation of RBD-ACE2 Binding Free Energy*

Bat ACE2 sequences missing from the databases at the time of analysis (i.e., *Rh. malayanus*, *Rh. thomasi* and *Chaerephon plicatus*) were reconstructed from our whole host genome sequences and completed with public data collected from NCBI/SRA. Briefly, a custom database comprising ACE2 sequences from Chiroptera or *Rhinolophus* sp. was created and queried using NCBI/SRA and our data using Diamond BlastX to identify ACE2-associated sequences. The assembly of identified sequences was then performed with Megahit.

The ACE2 genes from *Rhizomys pruinosus* (hoary bamboo rat) and *Hystrix brachyura* (Malayan porcupine) were similarly retrieved from genome assembly (respectively Gen-Bank GCA_009823505.1 and GCA_016801275.1) with BlastN against a custom database of Rodentia ACE2 genes.

To estimate the affinity between each of the *Sarbecovirus* RBD and ACE2 receptors studied here, we used the following in silico pipeline:

1. We used AlphaFold2 [66] to build structural models of ACE2 and RBDs. For each sequence, we generated 5 different models and selected the best model based on the AlphaFold2 pLDDT score.

2. For each potential RBD-ACE2 complex, we created 200 models of the complex using HADDOCK [67]. ACE2 and RBDs were first docked into a complex with the help of inter-subunit distance restraints. These restraints were defined between each pair of atoms at the interface of the SARS-CoV-2 RBD and human ACE2 based on the X-ray structure of this complex (PDB code 6M0J [68]). The generated models were then refined via molecular dynamics simulations in explicit solvent during the last step of the HADDOCK modeling pipeline.

3. The FoldX v. 5 scoring function [69] was used to estimate the RBD-ACE2 binding free energy from the ensemble of models generated with HADDOCK.

## 2.7. *Generation of Lentiviral Pseudoviruses and ACE2-Dependent Entry Assays*

Human codon-optimized sequences corresponding to the spike ectodomains of RtVN21-29, RtVN21-192, RtVN21-193 and RtVN21-201 were synthesized (GenScript) and cloned into a pVAX1 vector using a 2-fragment In-Fusion reaction with a sequence corresponding to the transmembrane and the C-terminal domains, with the 19 last amino acids deleted (positions 1222-1240 for RtVN21-29, positions 1224-1242 for RtVN21-192, positions 1216-1234 for RtVN21-193 and positions 1223-1241 for RtVN21-201).

Spike-pseudotyped lentiviral particles were produced via co-transfection in HEK-293T cells of the spike expression plasmid (5 µg of pVAX1-S-RtVN21-29, -192, -193, -201) together with a lentiviral backbone expressing the firefly luciferase (10 µg of pHAGE-CMV-Luc2-IRES-ZsGreen-W) and lentiviral helper plasmids encoding Gag-Pol, Tat and Rev (3.3 µg each of HDM-Hgpm2, HDM-tat1b and RC-CMV-Rev1b) using calcium phosphate precipitation [70]. Lentiviral particles were recovered 48 h post-transfection, clarified and tested for entry in HEK-293T expressing the human ACE2 receptor. Briefly, spike-pseudotyped lentiviral particles were serially diluted 3-fold and 50 µL was mixed with 50 µL of a HEK-293T-hACE2 cell suspension containing 4.10⁵ cells/mL and incubated in a 96-well white plate for 72 h. Firefly luciferase activity was then assessed by adding 100 µL
of BrightGlo substrate (Promega, Madison, WI, USA), and luminescence was measured using a Berthold Centro XS luminometer.

## 3. Results

As part of a global study assessing the viral diversity carried by horseshoe bats from Northern Vietnam, we report in this study the identification and the molecular characterization of seven novel sarbecoviruses.

### 3.1. *Sarbecovirus Detection Rate in Vietnamese Bats*

A total of 214 insectivorous bat rectal swabs were collected in Son La Province between February and July 2021 (22 February to 3 March and 4–13 July). They belonged to 19 bat species: *Aseliscus stoliczkanus* (N = 15), *Barbastrella* cf. *darjelingensis* (N = 1), *Hipposideros armiger* (N = 47), *Hipposideros larvatus* (N = 10), *Hipposideros gentilis* (N = 14), *Hypsugo* sp. (N = 8), *Megaderma spasma* (N = 3), *Miniopterus fuliginosus* (N = 10), *Myotis horsfieldii* (N = 1), *Myotis* sp. (N = 1), *Pipistrellus abramus* (N = 20), *Rhinolophus affinis* (N = 3), *Rhinolophus macrotis* (N = 1), *Rhinolophus marshalli* (N = 4), *Rhinolophus microglossus* (N = 9), *Rhinolophus pearsonii* (N = 3), *Rhinolophus pusillus* (N = 7), *Rhinolophus siamensis* (N = 1), *Rhinolophus thomasi* (N = 9) and *Taphozous* sp. (N = 48). Four *Rousettus* sp. frugivorous bat rectal swabs were also collected.

The pan-coronavirus PCR screening identified 21 coronaviruses, among which *Sarbecovirus*-related sequences were identified in seven of the eight *Rhinolophus thomasi* and in one of the *Rhinolophus siamensis* individuals. No sarbecovirus was identified in the other insectivorous or frugivorous bat species, suggesting a strong association of Vietnamese sarbecoviruses with their *Rhinolophus* sp. insectivorous hosts. One sarbecovirus originated from Khẩy lầu hamlet (Sốp cộp district, 220 km west from Hanoi), while the others were identified 80 km away, in two caves in Pha nhê hamlet (Mộc Châu district, 140 km west from Hanoï). Interestingly, two *Rhinolophus* species (i.e., *Rh. thomasi* and *Rh. siamensis*) were sampled in the same cave at the same time, highlighting the fact that different bat species live in sympatry in the same cave (Table 1). We tested oral swabs and urine samples collected from the same bat individuals for the presence of sarbecoviruses, and only two (RtVN21-193 and RtVN21-194) of the eight rectal swab-positive individuals also tested positive with oral swabs. No urine sample was positive, which suggests enteric shedding of these bat sarbecoviruses.

### 3.2. *Genomic Characterization of Rhinolophus* sp. *Sarbecoviruses*

The complete genomes of seven of the eight *Sarbecovirus* rectal swab strains were obtained after *Betacoronavirus*-enriched NGS. We were not able to obtain the complete genome of the oral swab-positive sarbecoviruses.

We first compared the genetic relationships at the gene level between these strains and identified different patterns of sequence clustering depending on the gene considered. For example, the clade that comprises RtVN21-197/-200 and -201 clustered with RtVN21-29 in ORF1a and with RtVN21-192 in ORF1b and in the spike, while they clustered with a subclade formed by RtVN21-29 and -192 in the nucleoprotein (Figure 2A). These changes in tree positioning suggest that recombinations occurred during RtVN21-29 and -192 evolution, while RtVN21-197/-200/-201 and RtVN21-193/RsVN21-195 always clustered together, suggesting higher stability of their genomes. To confirm that RtVN21-29 and RtVN21-192 were subject to multiple recombination events, we performed a full exploratory recombination analysis of the seven strains of Vietnamese sarbecoviruses, and observed that only RtVN21-29 and -192 presented recombinations (Figure 2B): RtVN21-29 presented three insertions for which the donor sequences were close to RtVN21-193 and RsVN21-195 (one in the ORF1ab and two in the S1 domain of the spike), and RtVN21-192 presented four insertions originating from strains close to RtVN21-197 and -201.
<table><caption>**Table 1.** Characteristics of *Sarbecovirus*-positive rectal swabs.</caption><thead><tr><th>Virus Strain</th><th>Province</th><th>District</th><th>Commune</th><th>Hamlet</th><th>Latitude (DDdddd')</th><th>Longitude (DDdddd')</th><th>Elevation (meters)</th><th>Bat Species</th><th>Bat COI Accession No.</th><th>Bat Weight (g)</th><th>Genome Coverage</th><th>CoV Accession No.</th></tr></thead><tbody><tr><td>RtVN21-29</td><td>Sơn La</td><td>Sốp cộp</td><td>Mường và</td><td>Khẩy lầu</td><td>20.882820</td><td>103.695280</td><td>936</td><td>Rhinolophus thomasi</td><td>OR427971</td><td>10.7</td><td>100%</td><td>OR261262</td></tr><tr><td>RtVN21-192</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 1</td><td>20.736283</td><td>104.504630</td><td>1059</td><td>Rhinolophus thomasi</td><td>OR435217</td><td>8.0</td><td>100%</td><td>OR261263</td></tr><tr><td>RtVN21-193</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus thomasi</td><td>OR427972</td><td>9.3</td><td>100%</td><td>OR261264</td></tr><tr><td>RtVN21-194</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus thomasi</td><td>OR427973</td><td>8.2</td><td>10.35%</td><td>nd</td></tr><tr><td>RsVN21-195</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus siamensis</td><td>OR427974</td><td>7.1</td><td>100%</td><td>OR261265</td></tr><tr><td>RtVN21-197</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus thomasi</td><td>OR427975</td><td>7.9</td><td>100%</td><td>OR261266</td></tr><tr><td>RtVN21-200</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus thomasi</td><td>OR427976</td><td>8.0</td><td>100%</td><td>OR261267</td></tr><tr><td>RtVN21-201</td><td>Sơn La</td><td>Mộc Châu</td><td>Lóng Sập</td><td>Pha nhê/cave 2</td><td>20.734502</td><td>104.505343</td><td>1073</td><td>Rhinolophus thomasi</td><td>OR427977</td><td>7.3</td><td>100%</td><td>OR261268</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Genomic characterization of Vietnamese sarbecoviruses. (A) Phylogenetic analysis of ORF1a, ORF1b, spike, ORF6, ORF7a and nucleoprotein genes of Vietnamese sarbecoviruses. Strains are represented by symbols and colored according to their position in the different clades. (B) Recombination analysis of Vietnamese sarbecoviruses. Recombinant fragments and corresponding breakpoints are indicated below each backbone and colored according to the donor genome sequence. Breakpoint coordinates correspond to the position in the alignment available at 10.6084/m9 figshare.23821428 (accessed on 6 September 2023). SARS-CoV-2 genome (NC_045512) served as reference for protein domains.</figcaption></figure>

Unsurprisingly, the distance matrices performed at both the gene and at the protein levels revealed that the less conserved proteins are the spike, the ORF8, the ORF3a and the ORF6 proteins, which exhibited 78%, 74%, 81% and 89% nucleotide identity between distant Vietnamese isolates, respectively (Figure S1).

Interestingly, different sarbecoviruses that belong to different clades were sampled in the same cave at the same time from two different *Rhinolophus* species (*thomasi* and *siamensis*), suggesting that at least two different lineages (one apparently restricted to
*Rh. thomasi* and one shared by *Rh. thomasi* and *Rh. siamensis*) of sarbecoviruses circulate concomitantly in the area. Conversely, viruses (RtVN21-29 and -192) that shared a common ancestor (at least for the nucleoprotein) are found in bats 80 km distant from each other despite the fact that insectivorous bats are not prone to long migrations, which suggests that bats living in caves not sampled in this study should harbor other close viral species (Figure 2).

At a more global scale, Vietnamese sarbecoviruses are clustered together but are quite distant from the SARS-CoV clade, while three distinct Vietnamese clades are clearly observed, both at the whole genome and at the spike gene level (Figure 3), which is consistent with previous observations (Figure 2). Similar conclusions could be drawn when looking at the RNA-dependent RNA polymerase (RdRP), at the nucleoprotein or at the N-terminal domain (NTD) and receptor-binding domains (RBDs) of the spike (Figure S2). Interestingly, African and European *Sarbecovirus* isolates clustered together and are distinct from Asian isolates, in addition to several bat species-specific clusters that are observed for *Rh. stheno*, *Rh. sinicus* and *Rh. hipposideros*, for example. This suggests that several sarbecoviruses may be specific to a given bat species, and that the geographical specificity observed between African/European and Asian *Sarbecovirus* isolates may only be explained by the presence of different bat populations in these continents.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Heatmap of nucleotide identities observed at the whole genome (A) and at the spike (B) levels between Vietnamese and 221 representative *Sarbecovirus* genomes. Aligned sequences were sorted by similarity. Low-to-high nucleotide identities are represented from blue to yellow. SARS-CoV- and SARS-CoV-2-related sequences are highlighted by white squares, and bat-specific clusters are highlighted by red squares. Identity matrices were calculated with CLC Main Workbench 21.0.4 (Qiagen) and the heatmap was created using GraphPad Prism 9.5.1 (GraphPad software). AS = Asian isolates; AF/EU = African/European isolates; Rh. = *Rhinolophus*. The complete identity matrices are presented in Table S3.</figcaption></figure>

Phylogenetic reconstructions of the evolution history of Vietnamese sarbecoviruses, along with representative SARS-CoV- and SARS-CoV-2-related animal viruses, placed Vietnamese strains differently when looking at the RdRP or at the spike genes, suggesting that recombinations may have occurred during Vietnamese *Sarbecovirus* evolution (Figure 4). For example, RtVN21-197/-200/-201 clustered together in the RdRP, with RtVN21-29 at a more basal position, but in the same clade. In addition, RtVN21-193/RsVN21-195 again clustered together but are placed in a sister clade of the one formed by RtVN21-197/-200/-
201/-192/-29. Conversely, these clusters are clearly distinct and distant from each other in the spike gene.

To evidence that recombinations have occurred during Vietnamese *Sarbecovirus* evolution, which could explain these different tree positionings, we performed a full exploratory recombination analysis, and identified 24 recombination events (Table S4). Some events were shared by two or more Vietnamese sarbecoviruses while other were unique. For example, the recombinant fragment #343 (starting at position 22,460, which corresponds to the S1 domain of the spike) was detected in four *Rh. thomasi* sequences (RtVN21-192, -197, -200 and -201), while fragment #758 (starting at position 3318, which corresponds to the ORF1a) was restricted to RtVN21-29. Donor sequences were sometimes unique (for example, RtVN21-192 for event #31 impacting RtVN21-197, -200 and -201) or multiple, meaning that the donor sequence was not clearly identified, possibly due to the lack of an exhaustive inventory of sarbecoviruses circulating in the area.

Interestingly, two Chinese *Sarbecovirus* strains, Rhinolophus sinicus Rs4231 and Rhinolophus affinis GD2017F, identified respectively in 2013 and 2017, clustered with Vietnamese strains at the RdRP level, suggesting that similar sarbecoviruses are shared by different rhinolophid bat species living in remote areas (Figure S3). While identified in the same *Rh. thomasi* bat colony, when looking at the spike gene, the RtVN21-197/-200/-201/-192 clade was closely related to several *Rh. sinicus*, *Rh. affinis* and *Aselliscus stoliczkanus* Chinese strains while the RtVN21-193/RsVN21-195 clade was more closely related to *Rh. pusillus* and *Rh. ferrumequinum* strains from China (Figure S4). More globally, even if local clusters of bat species-specific sarbecoviruses are observed, it seems that no strict host specificity was observed during *Sarbecovirus* evolution, as illustrated by many *Rh. sinicus*-associated clusters that frequently include sarbecoviruses from other bat species. This suggests that exchanges of related sarbecoviruses are frequent among different bat species that live in sympatry in the same colonies.

### 3.3. Prediction of Zoonotic Potential of Vietnamese Sarbecoviruses

#### 3.3.1. In Silico Evaluation of RBD-ACE2 Binding Free Energy

To study the potential risk of infection with Vietnamese sarbecoviruses in different vertebrate hosts, we used an in silico pipeline to estimate the binding free energy of the RBDs of representative sarbecoviruses (including SARS-CoV and SARS-CoV-2) and the ACE2 receptors of various animal species. Representative sarbecoviruses were selected from each of the three major clades in the RBD phylogenetic tree (Figure S5A). Animal species were selected to cover the diversity of rhinolophid bats that carry sarbecoviruses, and various domestic (pigs, cattle, dogs, cats, horses, chickens) and wild (rodents, pangolins, civets, porcupines, mustelids) animals that have been demonstrated or suspected to act as intermediate hosts for certain sarbecoviruses were selected.

Although SARS-CoV and SARS-CoV-related bat and civet sarbecoviruses showed predicted good affinities for human, pig, cat, raccoon dog, horse, porcupine, mustelid, rabbit and alpaca ACE2s, more distant SARS-CoV-related bat sarbecoviruses (such as *Rh. malayanus* BANAL-20-116 or Vietnamese *Rh. thomasi* and *Rh. siamensis* viruses) showed weak to very weak affinity for any human or animal ACE2 tested (Figure 5). Interestingly, we identified two large deletions in the receptor-binding motif of the RBDs of Vietnamese sarbecoviruses involving known amino-acid contact residues between virus RBDs and ACE2 receptors (Figure S5B), which are likely responsible for the apparent inability of our viruses to bind efficiently to any animal ACE2.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Phylogenetic analysis of Vietnamese and 221 representative sarbecoviruses. Trees were inferred for the RdRP and spike genes, and midpoint-rooted for clarity. Bootstraps above 50 are represented by black dots. Sarbecoviruses from *Rhinolophus* sp. bat species are represented by colored squares while *Chaerephon plicatus* and *Aselliscus stoliczkanus* sarbecoviruses are represented by dashed lines. Complete phylogenetic reconstructions are available in Figures S3 (RdRP) and S4 (spike).</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Estimation of the binding free energy of representative sarbecovirus RBDs for different animal ACE2 receptors. Scale was set up based on the binding free energy of human SARS-CoV-2 with human ACE2 (−10.0; white). Strong predicted binding free energy (below −10.0) is represented by shades of green while weak binding free energy (above −10.0) is represented by shades of red. “R” represents ACE2 sequence that was manually reconstructed from public and private datasets (see Section 2).</figcaption></figure>

### 3.3.2. ACE2-Dependent Entry Assay of Vietnamese Sarbecoviruses

To confirm these in silico predictions, we generated lentiviral particles pseudotyped either with the spike protein of one representative strain of Vietnamese sarbecoviruses (i.e., RtVN21-192, 193, -201 or -29) or with SARS-CoV-2 Wuhan (strain BetaCoV/France/IDF0372/2020, GISAID accession number EPI_ISL_406596). We tested the ability of these pseudotypes to enter into 293T-hACE2 cells (Figure 6) and observed that none of the Vietnamese strains were able to enter cells, unlike SARS-CoV-2. This result confirmed the in silico predictions and suggests that the deletion observed in the receptor-binding motif (Figure S5B) influences the ability of Vietnamese sarbecoviruses to use the human ACE2 receptor to infect humans.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** ACE2-dependent entry assay of *Rhinolophus* sp. sarbecovirus pseudotypes. Results of a single experiment performed in triplicate are expressed in relative luminescence units (RLUs) produced by the firefly luciferase expressed from the lentiviral vector. Error bars indicate SD. Statistical significance was assessed using one-way ANOVA and was compared with the RLU signal of an empty vector as a reference: ns: non-significant; ****: *p* < 0.0001.</figcaption></figure>

**4. Discussion**

Southeast Asia is a major hotspot for coronavirus emergence due to the huge richness in bat species that live and interact with high densities of human and domestic animal populations. This proximity is a key component of the processes of emergence of bat-borne zoonotic viruses. Direct (bat-to-human) or indirect (via an intermediate host, also known as a bridge host) spillovers could occur, and may lead to epidemics or pandemics, as observed in recent decades. Monitoring the circulation of coronaviruses in their bat reservoirs before their emergence in human or in domestic animal populations and understanding the molecular and epidemiological factors that shape spillover events remain challenges for public health.

In a work aiming at modeling the risk of *Sarbecovirus* spillover to humans, Sánchez et al. identified Northern Vietnam as one of the areas of Southeast Asia that harbors the highest diversity of bat-associated SARS-related coronaviruses [71]. Similarly, Hassanin et al. predicted that bats from Northern Vietnam could harbor SARS-CoV-like viruses [72]. Apart from these predictions, sarbecoviruses were never identified from Vietnam before this study, although extensive sampling efforts have been made. Here, we report the first identification and characterization of *Rhinolophus* sp. sarbecoviruses in Northern
Vietnam. Numerous insectivorous bat species were tested for the presence of coronaviruses, including horseshoe bats, which are known to be a major reservoir of sarbecoviruses [7], and only *Rhinolophus* bats (i.e., *Rh. thomasi* and *Rh. siamensis*) were found to be positive for sarbecoviruses. This finding revealed a strong association of Vietnamese sarbecoviruses with their *Rhinolophus* hosts, and indeed it has been proposed that such strong association reflects the long coevolution of bat-associated sarbecoviruses with their bat reservoirs [7].

Phylogenetic characterization of Vietnamese sarbecoviruses revealed three major find-
ings. First, the changes in tree positioning and the recombination events observed between
the different strains suggest that recombination was a major driver of *Sarbecovirus* evolution
and a possible source of novel viruses, as supposed for SARS-CoV-2 [73]. Second, the ob-
servation that different *Rh. thomasi* sarbecoviruses (RtVN21-193 vs. RtVN21-197, -200, -201)
belonging to different clades were sampled in the same cave at the same time shows that
different lineages of viruses are circulating concomitantly in the same area, which favors
recombination events. Finally, the fact that the same virus (presenting among others the
same receptor-binding domain) is found in different *Rhinolophus* species (RtVN21-193 from
*Rh. thomasi* and RsVN21-195 from *Rh. siamensis*) suggests that cross-species transmission
among bats is a major trait of sarbecoviruses (Figure 4, [25,72]). This observation suggests
that, even if sarbecoviruses coevolved with their bat hosts, transmission of viruses across
different bat species is frequent, which favors the emergence of novel viruses through
recombinations or other evolutionary processes. This in part likely reflects bat roosting
behavior and propensity to share the same or close habitats, which again favor the ex-
change of viruses and recombination. Understanding the drivers of the emergence of new
pathogens is one of the major challenges in the ecology of infectious diseases. It has been
proposed that Southeast Asia (comprising, among others, Southern China, Northern Laos
and Northern Vietnam) constitutes a hotspot for coronavirus emergence because of a favor-
able bioclimatic environment and a high diversity of bat species [61,72]. The presence of
this hotspot is due in part to the climate, the landscape and the human population density
of the region [74,75]. In a recent study, Beyer et al. reported, for example, the impact of
climate change on the geographical distribution of bats, which has resulted in an increase
in bat species richness in specific locations that could have been the origin of *Sarbecovirus*
emergence [76].

In line with previous observations that showed that coronavirus carriage in bats was
mainly enteric [27,46], we observed that few (25%) positive individuals presented a positive
oral swab, and none presented a positive urine sample. Indeed, a recent study showed that
the ACE2 mRNA expression level within bats greatly varies between organs and is higher
in the intestine than in the kidney, for example [77]. Enteric shedding of bat sarbecoviruses
therefore constitutes a risk of spillover, especially when considering people frequently
exposed to bat guano such as guano collectors, people working in guano farms or people
visiting caves who are exposed to guano dust [47,49]. Knowing the ability of Vietnamese
*Rhinolophus* sp. sarbecoviruses to infect humans is therefore crucial. There is currently no
evidence to suggest that these viruses pose a significant threat to public health. Therefore,
to decipher the zoonotic potential of these sarbecoviruses, we evaluated their ability to
enter mammalian cells using ACE2 receptors.

First, we evaluated in silico the interaction strength between the RBDs of Vietnamese
sarbecoviruses and different mammalian ACE2 receptors. We did not restrict the analysis to
the bat/human interface because, even if the intermediate host for SARS-CoV-2 emergence
is yet unknown, past emergence of coronaviruses frequently implies the presence of a
bridge host [6], such as domestic (dogs, pigs, dromedary camels) or wild animals (palm
civets, raccoon dogs) that are bred as food resources for humans in several Southeast
Asian countries, including Vietnam. The analysis indicated no robust interaction between
Vietnamese *Sarbecovirus* RBDs and human or wild or domestic animal ACE2. The predictive
value of these results is supported by the accurate prediction of the positive interaction
of human ACE2 with SARS-CoV-2 RBDs, along with the strong affinity of SARS-CoV
for swine ACE2, as previously reported [78,79]. When looking at the receptor-binding
motif (RBM) of Vietnamese sarbecoviruses, we observed two important deletions that may directly impact the affinity of their RBDs to mammalian ACE2 receptors by changing the structure of the spike protein. In nature, hundreds of sarbecoviruses like those described here have been found in bats, but only a small fraction of them have the ability to infect cells using ACE2. It has been proposed that bat sarbecoviruses be classified according to their RBD sequence and its subsequent ACE2 utilization [42]. Guo et al. reported that among the four clades of sarbecoviruses that can be identified, only Clade I viruses (harboring no deletion in the RBM) can utilize endogenous ACE2 for cell entry, and only these viruses are able to replicate in vitro or in vivo. Our viruses fall within Clade II of this classification, and initial attempts to isolate such viruses from field samples or propagate in vitro from reverse genetics were unsuccessful. Authors have also showed that in the presence of high levels of trypsin, some of the spike proteins of Clade II viruses were capable of mediating infection and replication in human (including Caco-2 intestine cell line) and bat intestine primary cells while TMPRSS2 and furin endogenous proteases had no effect on Clade II virus entry [42].

Apart from the rich trypsin environment, other factors such as the temperature may impact the infection and replication abilities of sarbecoviruses. Partridge et al. observed that the whole spike (and not the S1 or RBD domains) of SARS-CoV-2 exhibits a temperature-dependent cellular interaction; however, the underlying mechanisms are not fully understood [80]. Bats are the only mammalian species that is able to fly, which required evolutionary adaptations for them to survive the elevated body temperatures occurring during flight [81,82]. Therefore, it can be proposed that the unique temperature regulation of bats may impact the ability of bat sarbecoviruses to bind efficiently to bat cells, thus preventing other homeothermic species from being infected by bat sarbecoviruses that coevolved over long periods with their hosts. In line with the presence of the two deletions in the receptor-binding motif, such findings suggest that these viruses are likely unable to cross the species barrier and to pose a risk to human or animal health.

In addition to the presence of two deletions in the RBM, several amino-acid substitutions were observed at the contact residues. The plasticity of the spike protein of sarbecoviruses, which allows them to adapt to a broad range of ACE2 receptors, is now well known, and such substitutions at contact residues allow several viruses to increase their affinity for ACE2 [27]. Molecular evolution analyses indicated that these key residues are under positive selection, suggesting that the spike protein of bat sarbecoviruses and mammalian ACE2 may have coevolved over a long time and adapted to selection pressure from each other [83,84]. The substitutions observed in Vietnamese contact residues may therefore be a signature of the inability of a given sarbecovirus to interact efficiently with ACE2, as observed, for example, for Rhinolophus stheno RsYN09 [85], Rhinolophus affinis YN2020A [55], Rhinolophus sinicus Rs4247, Rhinolophus ferrumequinum Rf4092 [26] and Rhinolophus pusillus F46 [86] sarbecoviruses. These simple features may serve as preliminary markers of the ability or inability of a new sarbecovirus to bind to mammalian host cells via ACE2, and possibly infect them, as studied in depth when SARS-CoV-2 emerged [87].

Since the RBD alone cannot predict the behavior of the whole spike in the context of a mature virion, especially given the presence of the two deletions in the receptor-binding motif, we constructed pseudotyped lentiviruses to measure the interaction of the spike protein of Vietnamese coronaviruses with human ACE2. The results of entry assays confirmed that Vietnamese sarbecoviruses were not able to efficiently interact with hACE2, preventing their entry into 293T human cells. Indeed, hundreds of sarbecoviruses have been found in bats, but only a fraction of them have the ability to use ACE2 [42]. It can be hypothesized that other proteins act as alternative receptors, in complement to or in place of ACE2 [88]. For example, an ACE2-independent entry of SARS-CoV-2 into different human cells via the neuropilin-1 [89] or glycosaminoglycans [90], which are much more expressed at the surface of certain cells (including intestine and kidney cells) than hACE2, has been reported [91].
While apparently lacking the ability to interact with human (and animal) ACE2, and hence unlikely to be zoonotic without the acquisition of several sequence changes, Vietnamese sarbecoviruses enrich the pool of viruses circulating in the area and present opportunities for homologous recombination that could lead, for example, to a new virus able to cross the species barrier. Also, recombination could lead to increased virulence compared with the parental strains [92], to the capability to escape neutralization by antibodies able to neutralize the parental strain [93] or to changes in host or tissue tropism [94]. Therefore, identifying novel coronaviruses should not be restricted to the (bat) reservoir but also should include any potential bridge hosts that could be infected by emerging coronaviruses. The identification of such hosts could easily be implemented by applying a straightforward strategy that includes in silico binding analyses (conducted over large RBD/ACE2 datasets) in order to capture a global picture of the host range and risk of cross-species transmission of emerging sarbecoviruses, and to improve the surveillance of specific interfaces at higher risk of emergence.

**5. Conclusions**

Southeast Asia hosts a great diversity of wildlife and is undergoing important climate and land-use changes that can increase contact between wildlife and humans. Continued and expanded surveillance of bat populations and other key wild and domestic animals in Southeast Asia, and especially Vietnam, is a crucial component of future pandemic preparedness and prevention. Such programs should include the surveillance of sentinel interfaces and the use of novel and easy-to-implement molecular and bioinformatics tools to identify high-risk wildlife-to-human interfaces.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v15091897/s1: Figure S1: Amino-acid (lower) and nucleotide (upper) identity matrices of Vietnamese bat sarbecoviruses performed at the gene level. Figure S2: Heatmap of nucleotide identities observed at the RNA-dependent RNA polymerase (A), nucleoprotein (B), N-terminal domain (C) and receptor-binding domain (D) levels between Vietnamese and 221 representative Sarbecovirus genomes. Figure S3: Extended phylogenetic tree of Vietnamese and 221 representative sarbecoviruses performed in the RNA-dependent RNA polymerase. Figure S4: Extended phylogenetic tree of Vietnamese and 221 representative sarbecoviruses performed in the spike protein. Figure S5: (A) Phylogenetic reconstruction of the evolutionary history of the RBDs of Vietnamese and 221 representative sarbecoviruses in light with ACE2 usage (yellow). Sequences selected for in silico RBD/ACE2 interactions are noted with a red asterisk; (B) Alignments of the receptor-binding motif (RBM) of Vietnamese and representative sarbecoviruses. SARS-CoV (top) and SARS-CoV-2 (bottom) contact residues are highlighted. Deletions in the RBM are presented in gray. Table S1: Detailed sampling of bats. Table S2: Acknowledgment of sharing of SARS-CoV-2-related genome sequences from the GISAID database. Table S3: Complete identity matrices used to generate the heatmaps presented in Figure 3 (whole genome and spike) and in Figure S2 (RdRP, N, S-NTD and S-RBD). Table S4: Detailed results of the recombination analysis conducted between Vietnamese and 221 representative Sarbecovirus genomes.

**Author Contributions:** Conceptualization: S.T., T.C.T., B.R., M.B., S.M., V.D.T., S.v.d.W., V.S.N. and M.E.; methodology: all authors; software: S.T., M.B. and T.B.; validation: S.T., T.C.T., B.R., M.B., S.M., V.D.T., S.v.d.W., V.S.N. and M.E.; field investigation: T.C.T., T.N.D., T.V.P., N.T.Y., H.N.A., T.H.S., P.T.A. and V.S.N.; writing—original draft preparation: S.T.; writing—review and editing: all authors; supervision: M.B., S.M., S.v.d.W., V.S.N. and M.E.; project administration: M.E.; funding acquisition: S.v.d.W. and M.E. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Institut Pasteur "Covid Taskforce", by the French Agence Nationale de la Recherche ANR (BAT-COR-ASIA project, grant no. ANR-22-CE35-0004) and in part by the H2020 project 101003589 (RECOVER) and the Labex IBEID (ANR-10-LABX62-IBEID). Vu Dinh Thong was supported by the Vietnam Academy of Science and Technology through the project ID VAST04.03/21-22.

**Institutional Review Board Statement:** The animal study was approved by the Ministry of Health (Vietnam) (protocol code 4026/QD-BYT on the 21 September 2020).
**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The complete genome sequences of Vietnamese sarbecoviruses and bat-associated COIs were deposited in GenBank (accession numbers OR261262-68 and OR427971-77, respectively). GISAID identifier: EPI_SET_230623sh (DOI: 10.55876/gis8.230623sh). All genome sequences and associated metadata in this dataset are published in GISAID's EpiCoV database. The alignments used to reconstruct phylogenies presented in this study and the complete RDP5 recombination analysis are openly available in FigShare (DOI: 10.6084/m9 figshare.23821428 and 10.6084/m9 figshare.23821719, respectively).

**Acknowledgments:** The authors want to acknowledge Amal Bennouna for her help with Q-GIS, and all Vietnamese colleagues for their valuable help with field work. Vu Dinh Thong was supported by the Vietnam Academy of Science and Technology through the project ID VAST04.03/21-22.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Han, H.J.; Wen, H.L.; Zhou, C.M.; Chen, F.F.; Luo, L.M.; Liu, J.W.; Yu, X.J. Bats as reservoirs of severe emerging infectious diseases. *Virus Res.* **2015**, *205*, 1–6. [CrossRef] [PubMed]

2. Letko, M.; Seifert, S.N.; Olival, K.J.; Plowright, R.K.; Munster, V.J. Bat-borne virus diversity, spillover and emergence. *Nat. Rev. Microbiol.* **2020**, *18*, 461–471. [CrossRef] [PubMed]

3. Li, W.; Shi, Z.; Yu, M.; Ren, W.; Smith, C.; Epstein, J.H.; Wang, H.; Crameri, G.; Hu, Z.; Zhang, H.; et al. Bats are natural reservoirs of SARS-like coronaviruses. *Science* **2005**, *310*, 676–679. [CrossRef]

4. Anthony, S.J.; Gilardi, K.; Menachery, V.D.; Goldstein, T.; Ssebide, B.; Mbabazi, R.; Navarrete-Macias, I.; Liang, E.; Wells, H.; Hicks, A.; et al. Further Evidence for Bats as the Evolutionary Source of Middle East Respiratory Syndrome Coronavirus. *mBio* **2017**, *8*, e00373-17. [CrossRef] [PubMed]

5. Zhou, P.; Fan, H.; Lan, T.; Yang, X.L.; Shi, W.F.; Zhang, W.; Zhu, Y.; Zhang, Y.W.; Xie, Q.M.; Mani, S.; et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. *Nature* **2018**, *556*, 255–258. [CrossRef] [PubMed]

6. Keusch, G.T.; Amuasi, J.H.; Anderson, D.E.; Daszak, P.; Eckerle, I.; Field, H.; Koopmans, M.; Lam, S.K.; Das Neves, C.G.; Peiris, M.; et al. Pandemic origins and a One Health approach to preparedness and prevention: Solutions based on SARS-CoV-2 and other RNA viruses. *Proc. Natl. Acad. Sci. USA* **2022**, *119*, e2202871119. [CrossRef]

7. Ruiz-Aravena, M.; McKee, C.; Gamble, A.; Lunn, T.; Morris, A.; Snedden, C.E.; Yinda, C.K.; Port, J.R.; Buchholz, D.W.; Yeo, Y.Y.; et al. Ecology, evolution and spillover of coronaviruses from bats. *Nat. Rev. Microbiol.* **2022**, *20*, 299–314. [CrossRef]

8. Rahman, S.; Ullah, S.; Shinwari, Z.K.; Ali, M. Bats-associated beta-coronavirus detection and characterization: First report from Pakistan. *Infect. Genet Evol.* **2023**, *108*, 105399. [CrossRef]

9. Speranskaya, A.S.; Artiushin, I.V.; Samoilov, A.E.; Korneenko, E.V.; Khabudaev, K.V.; Ilina, E.N.; Yusefovich, A.P.; Safonova, M.V.; Dolgova, A.S.; Gladkikh, A.S.; et al. Identification and Genetic Characterization of MERS-Related Coronavirus Isolated from Nathusius' Pipistrelle (Pipistrellus nathusii) near Zvenigorod (Moscow Region, Russia). *Int. J. Environ. Res. Public Health* **2023**, *20*, 3702. [CrossRef]

10. Wang, J.; Pan, Y.F.; Yang, L.F.; Yang, W.H.; Lv, K.; Luo, C.M.; Wang, J.; Kuang, G.P.; Wu, W.C.; Gou, Q.Y.; et al. Individual bat virome analysis reveals co-infection and spillover among bats and virus zoonotic potential. *Nat. Commun.* **2023**, *14*, 4079. [CrossRef]

11. Van Brussel, K.; Mahar, J.E.; Ortiz-Baez, A.S.; Carrai, M.; Spielman, D.; Boardman, W.S.J.; Baker, M.L.; Beatty, J.A.; Geoghegan, J.L.; Barrs, V.R.; et al. Faecal virome of the Australian grey-headed flying fox from urban/suburban environments contains novel coronaviruses, retroviruses and sapoviruses. *Virology* **2022**, *576*, 42–51. [CrossRef] [PubMed]

12. Muzeniek, T.; Perera, T.; Siriwardana, S.; Bas, D.; Kaplan, F.; Öruc, M.; Becker-Ziaja, B.; Perera, I.; Weerasena, J.; Handunnetti, S.; et al. Full Genome of batCoV/MinFul/2018/SriLanka, a Novel Alpha-Coronavirus Detected in *Miniopterus fuliginosus*, Sri Lanka. *Viruses* **2022**, *14*, 337. [PubMed]

13. Wang, N.; Luo, C.M.; Yang, X.L.; Liu, H.Z.; Zhang, L.B.; Zhang, W.; Li, B.; Zhu, Y.; Peng, C.; Shi, Z.L.; et al. Genomic Characterization of Diverse Bat Coronavirus HKU10 in *Hipposideros Bats. Viruses* **2021**, *13*, 1962. [CrossRef]

14. Hardmeier, I.; Aeberhard, N.; Qi, W.; Schoenbaechler, K.; Kraettli, H.; Hatt, J.M.; Fraefel, C.; Kubacki, J. Metagenomic analysis of fecal and tissue samples from 18 endemic bat species in Switzerland revealed a diverse virus composition including potentially zoonotic viruses. *PLoS ONE* **2021**, *16*, e0252534. [CrossRef] [PubMed]

15. Daszak, P.; Keusch, G.T.; Phelan, A.L.; Johnson, C.K.; Osterholm, M.T. Infectious Disease Threats: A Rebound to Resilience. *Health Aff.* **2021**, *40*, 204–211. [CrossRef]

16. Wang, N.; Li, S.Y.; Yang, X.L.; Huang, H.M.; Zhang, Y.J.; Guo, H.; Luo, C.M.; Miller, M.; Zhu, G.; Chmura, A.A.; et al. Serological Evidence of Bat SARS-Related Coronavirus Infection in Humans, China. *Virol. Sin.* **2018**, *33*, 104–107. [CrossRef]
17. Li, H.; Mendelsohn, E.; Zong, C.; Zhang, W.; Hagan, E.; Wang, N.; Li, S.; Yan, H.; Huang, H.; Zhu, G.; et al. Human-animal interactions and bat coronavirus spillover potential among rural residents in Southern China. *Biosaf. Health* **2019**, *1*, 84–90.
[CrossRef]

18. Evans, T.S.; Tan, C.W.; Aung, O.; Phyu, S.; Lin, H.; Coffey, L.L.; Toe, A.T.; Aung, P.; Aung, T.H.; Aung, N.T.; et al. Exposure to diverse sarbecoviruses indicates frequent zoonotic spillover in human communities interacting with wildlife. *Int. J. Infect. Dis.* **2023**, *131*, 57–64.
[CrossRef]

19. Temmam, S.; Montagutelli, X.; Herate, C.; Donati, F.; Regnault, B.; Attia, M.; Baquero Salazar, E.; Chretien, D.; Conquet, L.; Jouvion, G.; et al. SARS-CoV-2-related bat virus behavior in human-relevant models sheds light on the origin of COVID-19. *EMBO Rep.* **2023**, *24*, e56055.
[CrossRef]

20. Tan, C.C.S.; Trew, J.; Peacock, T.P.; Mok, K.Y.; Hart, C.; Lau, K.; Ni, D.; Orme, C.D.L.; Ransome, E.; Pearse, W.D.; et al. Genomic screening of 16 UK native bat species through conservationist networks uncovers coronaviruses with zoonotic potential. *Nat. Commun.* **2023**, *14*, 3322.
[CrossRef]

21. Apaa, T.; Withers, A.J.; Staley, C.; Blanchard, A.; Bennett, M.; Bremner-Harrison, S.; Chadwick, E.A.; Hailer, F.; Harrison, S.W.R.; Loose, M.; et al. Sarbecoviruses of British horseshoe bats; sequence variation and epidemiology. *J. Gen. Virol.* **2023**, *104*, 001859.
[CrossRef]

22. Alkhovsky, S.; Lenshin, S.; Romashin, A.; Vishnevskaya, T.; Vyshemirsky, O.; Bulycheva, Y.; Lvov, D.; Gitelman, A. SARS-like Coronaviruses in Horseshoe Bats (*Rhinolophus* spp.) in Russia, 2020. *Viruses* **2022**, *14*, 113.
[CrossRef]

23. Wells, H.L.; Letko, M.; Lasso, G.; Ssebide, B.; Nziza, J.; Byarugaba, D.K.; Navarrete-Macias, I.; Liang, E.; Cranfield, M.; Han, B.A.; et al. The evolutionary history of ACE2 usage within the coronavirus subgenus *Sarbecovirus*. *Virus Evol.* **2021**, *7*, veab007.
[CrossRef] [PubMed]

24. Tao, Y.; Tong, S. Complete Genome Sequence of a Severe Acute Respiratory Syndrome-Related Coronavirus from Kenyan Bats. *Microbiol. Resour. Announc.* **2019**, *8*, e00548-19.
[CrossRef]

25. Latinne, A.; Hu, B.; Olival, K.J.; Zhu, G.; Zhang, L.; Li, H.; Chmura, A.A.; Field, H.E.; Zambrana-Torrelio, C.; Epstein, J.H.; et al. Origin and cross-species transmission of bat coronaviruses in China. *Nat. Commun.* **2020**, *11*, 4235.
[CrossRef]

26. Hu, B.; Zeng, L.P.; Yang, X.L.; Ge, X.Y.; Zhang, W.; Li, B.; Xie, J.Z.; Shen, X.R.; Zhang, Y.Z.; Wang, N.; et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. *PLoS Pathog.* **2017**, *13*, e1006698.
[CrossRef] [PubMed]

27. Temmam, S.; Vongphayloth, K.; Baquero, E.; Munier, S.; Bonomi, M.; Regnault, B.; Douangboubpha, B.; Karami, Y.; Chrétien, D.; Sanamxay, D.; et al. Bat coronaviruses related to SARS-CoV-2 and infectious for human cells. *Nature* **2022**, *604*, 330–336.
[CrossRef]

28. Wacharapluesadee, S.; Tan, C.W.; Maneeorn, P.; Duengkae, P.; Zhu, F.; Joyjinda, Y.; Kaewpom, T.; Chia, W.N.; Ampoot, W.; Lim, B.L.; et al. Evidence for SARS-CoV-2 related coronaviruses circulating in bats and pangolins in Southeast Asia. *Nat. Commun.* **2021**, *12*, 972.
[CrossRef]

29. Delaune, D.; Hul, V.; Karlsson, E.A.; Hassanin, A.; Ou, T.P.; Baidaliuk, A.; Gámbaro, F.; Prot, M.; Tu, V.T.; Chea, S.; et al. A novel SARS-CoV-2 related coronavirus in bats from Cambodia. *Nat. Commun.* **2021**, *12*, 6563.
[CrossRef] [PubMed]

30. Kim, Y.; Son, K.; Kim, Y.S.; Lee, S.Y.; Jheong, W.; Oem, J.K. Complete genome analysis of a SARS-like bat coronavirus identified in the Republic of Korea. *Virus Genes* **2019**, *55*, 545–549.
[CrossRef]

31. Murakami, S.; Kitamura, T.; Suzuki, J.; Sato, R.; Aoi, T.; Fujii, M.; Matsugo, H.; Kamiki, H.; Ishida, H.; Takenaka-Uema, A.; et al. Detection and Characterization of Bat Sarbecovirus Phylogenetically Related to SARS-CoV-2, Japan. *Emerg. Infect. Dis.* **2020**, *26*, 3025–3029.
[CrossRef]

32. Starr, T.N.; Zepeda, S.K.; Walls, A.C.; Greaney, A.J.; Alkhovsky, S.; Veesler, D.; Bloom, J.D. ACE2 binding is an ancestral and evolvable trait of sarbecoviruses. *Nature* **2022**, *603*, 913–918.
[CrossRef] [PubMed]

33. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C.; et al. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. *Nature* **2003**, *426*, 450–454.
[CrossRef] [PubMed]

34. Letko, M.; Marzi, A.; Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. *Nat. Microbiol.* **2020**, *5*, 562–569.
[CrossRef]

35. Liu, Y.; Hu, G.; Wang, Y.; Ren, W.; Zhao, X.; Ji, F.; Zhu, Y.; Feng, F.; Gong, M.; Ju, X.; et al. Functional and genetic analysis of viral receptor ACE2 orthologs reveals a broad potential host range of SARS-CoV-2. *Proc. Natl. Acad. Sci. USA* **2021**, *118*, e2025373118.
[CrossRef] [PubMed]

36. Li, W.; Zhang, C.; Sui, J.; Kuhn, J.H.; Moore, M.J.; Luo, S.; Wong, S.K.; Huang, I.C.; Xu, K.; Vasilieva, N.; et al. Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2. *EMBO J.* **2005**, *24*, 1634–1643.
[CrossRef]

37. Xiao, K.; Zhai, J.; Feng, Y.; Zhou, N.; Zhang, X.; Zou, J.J.; Li, N.; Guo, Y.; Li, X.; Shen, X.; et al. Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins. *Nature* **2020**, *583*, 286–289.
[CrossRef]

38. Ren, W.; Lan, J.; Ju, X.; Gong, M.; Long, Q.; Zhu, Z.; Yu, Y.; Wu, J.; Zhong, J.; Zhang, R.; et al. Mutation Y453F in the spike protein of SARS-CoV-2 enhances interaction with the mink ACE2 receptor for host adaption. *PLoS Pathog.* **2021**, *17*, e1010053.
[CrossRef]
39. Guo, H.; Hu, B.; Si, H.R.; Zhu, Y.; Zhang, W.; Li, B.; Li, A.; Geng, R.; Lin, H.F.; Yang, X.L.; et al. Identification of a novel lineage bat SARS-related coronaviruses that use bat ACE2 receptor. *Emerg. Microbes Infect.* **2021**, *10*, 1507–1514. [CrossRef]

40. Mou, H.; Quinlan, B.D.; Peng, H.; Liu, G.; Guo, Y.; Peng, S.; Zhang, L.; Davis-Gardner, M.E.; Gardner, M.R.; Crynen, G.; et al. Mutations derived from horseshoe bat ACE2 orthologs enhance ACE2-Fc neutralization of SARS-CoV-2. *PLoS Pathog.* **2021**, *17*, e1009501. [CrossRef]

41. Hu, Y.; Liu, K.; Han, P.; Xu, Z.; Zheng, A.; Pan, X.; Jia, Y.; Su, C.; Tang, L.; Wu, L.; et al. Host range and structural analysis of bat-origin RshSTT182/200 coronavirus binding to human ACE2 and its animal orthologs. *EMBO J.* **2023**, *42*, e111737. [CrossRef]

42. Guo, H.; Li, A.; Dong, T.Y.; Su, J.; Yao, Y.L.; Zhu, Y.; Shi, Z.L.; Letko, M. ACE2-Independent Bat Sarbecovirus Entry and Replication in Human and Bat Cells. *mBio* **2022**, *13*, e0256622. [CrossRef] [PubMed]

43. Aicher, S.M.; Streicher, F.; Chazal, M.; Planas, D.; Luo, D.; Buchrieser, J.; Nemcova, M.; Seidlova, V.; Zukal, J.; Serra-Cobo, J.; et al. Species-Specific Molecular Barriers to SARS-CoV-2 Replication in Bat Cells. *J. Virol.* **2022**, *96*, e0060822. [CrossRef]

44. Fujima, S.; Kosuke, Y.; Kimura, I.; Tukunaga, K.; Ito, J.; Sato, K. Determination of the factors responsible for host tropism of SARS-CoV-2-related bat coronaviruses. *bioRxiv* **2023**, preprint.

45. Yan, H.; Jiao, H.; Liu, Q.; Zhang, Z.; Xiong, Q.; Wang, B.J.; Wang, X.; Guo, M.; Wang, L.F.; Lan, K.; et al. ACE2 receptor usage reveals variation in susceptibility to SARS-CoV and SARS-CoV-2 infection among bat species. *Nat. Ecol. Evol.* **2021**, *5*, 600–608. [CrossRef] [PubMed]

46. Berto, A.; Anh, P.H.; Carrique-Mas, J.J.; Simmonds, P.; Van Cuong, N.; Tue, N.T.; Van Dung, N.; Woolhouse, M.E.; Smith, I.; Marsh, G.A.; et al. Detection of potentially novel paramyxovirus and coronavirus viral RNA in bats and rats in the Mekong Delta region of southern Viet Nam. *Zoonoses Public Health* **2018**, *65*, 30–42. [CrossRef]

47. Huong, N.Q.; Nga, N.T.T.; Long, N.V.; Luu, B.D.; Latinne, A.; Pruvot, M.; Phuong, N.T.; Quang, L.T.V.; Hung, V.V.; Lan, N.T.; et al. Coronavirus testing indicates transmission risk increases along wildlife supply chains for human consumption in Viet Nam, 2013–2014. *PLoS ONE* **2020**, *15*, e0237129. [CrossRef] [PubMed]

48. Phan, M.V.T.; Ngo Tri, T.; Hong Anh, P.; Baker, S.; Kellam, P.; Cotten, M. Identification and characterization of *Coronaviridae* genomes from Vietnamese bats and rats based on conserved protein domains. *Virus Evol.* **2018**, *4*, vey035. [CrossRef]

49. Latinne, A.; Nga, N.T.T.; Long, N.V.; Ngoc, P.T.B.; Thuy, H.B.; Long, P.T.; Phuong, N.T.; Quang, L.T.V.; Tung, N.; Nam, V.S.; et al. One Health Surveillance Highlights Circulation of Viruses with Zoonotic Potential in Bats, Pigs, and Humans in Viet Nam. *Viruses* **2023**, *15*, 790. [CrossRef]

50. Sikes, R.S.; Animal Care and Use Committee of the American Society of Mammalogists. Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. *J. Mammal.* **2016**, *97*, 663–688. [CrossRef]

51. Wilson, D.E.; Mittermeier, R.A. (Eds.) Handbook of the Mammals of the World. In *Bats*; Lynx Edicions: Barcelona, Spain, 2019; Volume 9.

52. Thong, V.D. *Taxonomy and Echolocation of Vietnamese Bats*; Publishing House for Science and Technology: Hanoi, Vietnam, 2021; 258p. (In Vietnamese)

53. Ratnasingham, S.; Hebert, P.D.N. BOLD: The Barcode of Life Data System (www.barcodinglife.org). *Mol. Ecol.* **2007**, *7*, 355–364.

54. Chu, D.K.; Leung, C.Y.; Gilbert, M.; Joyner, P.H.; Ng, E.M.; Tse, T.M.; Guan, Y.; Peiris, J.S.; Poon, L.L. Avian coronavirus in wild aquatic birds. *J. Virol.* **2011**, *85*, 12815–12820.

55. Wu, Z.; Han, Y.; Wang, Y.; Liu, B.; Zhao, L.; Zhang, J.; Su, H.; Zhao, W.; Liu, L.; Bai, S.; et al. A comprehensive survey of bat sarbecoviruses across China in relation to the origins of SARS-CoV and SARS-CoV-2. *Natl. Sci. Rev.* **2022**, *10*, nwac213.

56. C2A.A2C. Available online: https://gitlab.pasteur.fr/GIPhy/C2A.A2C (accessed on 27 June 2023).

57. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. *Brief Bioinform.* **2019**, *20*, 1160–1166. [PubMed]

58. Trifinopoulos, J.; Nguyen, L.T.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. *Nucleic Acids Res.* **2016**, *44*, W232–W235.

59. FigTree. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 27 June 2023).

60. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. *Virus Evol.* **2020**, *7*, veaa087. [PubMed]

61. Martin, D.; Rybicki, E. RDP: Detection of recombination amongst aligned sequences. *Bioinformatics* **2000**, *16*, 562–563. [CrossRef] [PubMed]

62. Padidam, M.; Sawyer, S.; Fauquet, C.M. Possible emergence of new geminiviruses by frequent recombination. *Virology* **1999**, *265*, 218–225.

63. Smith, J.M. Analyzing the mosaic structure of genes. *J. Mol. Evol.* **1992**, *34*, 126–129. [CrossRef]

64. Posada, D.; Crandall, K.A. Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 13757–13762. [PubMed]

65. Boni, M.F.; Posada, D.; Feldman, M.W. An exact nonparametric method for inferring mosaic structure in sequence triplets. *Genetics* **2007**, *176*, 1035–1047. [PubMed]
66. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. *Nature* **2021**, *596*, 583–589. [PubMed]

67. Van Zundert, G.C.P.; Rodrigues, J.P.G.L.M.; Trellet, M.; Schmitz, C.; Kastritis, P.L.; Karaca, E.; Melquiond, A.S.J.; van Dijk, M.; de Vries, S.J.; Bonvin, A.M.J.J. The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. *J. Mol. Biol.* **2016**, *428*, 720–725. [PubMed]

68. Lan, J.; Ge, J.; Yu, J.; Shan, S.; Zhou, H.; Fan, S.; Zhang, Q.; Shi, X.; Wang, Q.; Zhang, L.; et al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. *Nature* **2020**, *581*, 215–220.

69. Delgado, J.; Radusky, L.G.; Cianferoni, D.; Serrano, L. FoldX 5.0: Working with RNA, small molecules and a new graphical interface. *Bioinformatics* **2019**, *35*, 4168–4169. [CrossRef]

70. Crawford, K.H.D.; Eguia, R.; Dingens, A.S.; Loes, A.N.; Malone, K.D.; Wolf, C.R.; Chu, H.Y.; Tortorici, M.A.; Veesler, D.; Murphy, M.; et al. Protocol and Reagents for Pseudotyping Lentiviral Particles with SARS-CoV-2 Spike Protein for Neutralization Assays. *Viruses* **2020**, *12*, 513.

71. Sánchez, C.A.; Li, H.; Phelps, K.L.; Zambrana-Torrelio, C.; Wang, L.F.; Zhou, P.; Shi, Z.L.; Olival, K.J.; Daszak, P. A strategy to assess spillover risk of bat SARS-related coronaviruses in Southeast Asia. *Nat. Commun.* **2022**, *13*, 4380.

72. Hassanin, A.; Tu, V.T.; Curaveau, M.; Csorba, G. Inferring the ecological niche of bat viruses closely related to SARS-CoV-2 using phylogeographic analyses of Rhinolophus species. *Sci. Rep.* **2021**, *11*, 14276.

73. Lytras, S.; Hughes, J.; Martin, D.; Swanepoel, P.; de Klerk, A.; Lourens, R.; Kosakovsky Pond, S.L.; Xia, W.; Jiang, X.; Robertson, D.L. Exploring the Natural Origins of SARS-CoV-2 in the Light of Recombination. *Genome Biol. Evol.* **2022**, *14*, evac018.

74. Muylaert, R.L.; Kingston, T.; Luo, J.; Vancine, M.H.; Galli, N.; Carlson, C.J.; John, R.S.; Rulli, M.C.; Hayman, D.T.S. Present and future distribution of bat hosts of sarbecoviruses: Implications for conservation and public health. *Proc. Biol. Sci.* **2022**, *289*, 20220397.

75. Wang, Z.; Huang, G.; Huang, M.; Dai, Q.; Hu, Y.; Zhou, J.; Wei, F. Global patterns of phylogenetic diversity and transmission of bat coronavirus. *Sci. China Life Sci.* **2023**, *66*, 861–874.

76. Beyer, R.M.; Manica, A.; Mora, C. Shifts in global bat diversity suggest a possible role of climate change in the emergence of SARS-CoV-1 and SARS-CoV-2. *Sci. Total Environ.* **2021**, *767*, 145413.

77. Li, L.; Zhang, L.; Zhou, J.; He, X.; Yu, Y.; Liu, P.; Huang, W.; Xiang, Z.; Chen, J. Epidemiology and Genomic Characterization of Two Novel SARS-Related Coronaviruses in Horseshoe Bats from Guangdong, China. *mBio* **2022**, *13*, e0046322. [PubMed]

78. Chen, W.; Yan, M.; Yang, L.; Ding, B.; He, B.; Wang, Y.; Liu, X.; Liu, C.; Zhu, H.; You, B.; et al. SARS-associated coronavirus transmitted from human to pig. *Emerg. Infect. Dis.* **2005**, *11*, 446–448. [CrossRef]

79. Weingartl, H.M.; Copps, J.; Drebot, M.A.; Marszal, P.; Smith, G.; Gren, J.; Andova, M.; Pasick, J.; Kitching, P.; Czub, M. Susceptibility of pigs and chickens to SARS coronavirus. *Emerg. Infect. Dis.* **2004**, *10*, 179–184. [PubMed]

80. Partridge, L.J.; Urwin, L.; Nicklin, M.J.H.; James, D.C.; Green, L.R.; Monk, P.N. ACE2-Independent Interaction of SARS-CoV-2 Spike Protein with Human Epithelial Cells Is Inhibited by Unfractionated Heparin. *Cells* **2021**, *10*, 1419. [CrossRef] [PubMed]

81. Chionh, Y.T.; Cui, J.; Koh, J.; Mendenhall, I.H.; Ng, J.H.J.; Low, D.; Itahana, K.; Irving, A.T.; Wang, L.F. High basal heat-shock protein expression in bats confers resistance to cellular heat/oxidative stress. *Cell Stress Chaperones* **2019**, *24*, 835–849.

82. O'Shea, T.J.; Cryan, P.M.; Cunningham, A.A.; Fooks, A.R.; Hayman, D.T.; Luis, A.D.; Peel, A.J.; Plowright, R.K.; Wood, J.L. Bat flight and zoonotic viruses. *Emerg. Infect. Dis.* **2014**, *20*, 741–745.

83. Tian, J.; Sun, J.; Li, D.; Wang, N.; Wang, L.; Zhang, C.; Meng, X.; Ji, X.; Suchard, M.A.; Zhang, X.; et al. Emerging viruses: Cross-species transmission of coronaviruses, filoviruses, henipaviruses, and rotaviruses from bats. *Cell Rep.* **2022**, *39*, 110969.

84. Guo, H.; Hu, B.J.; Yang, X.L.; Zeng, L.P.; Li, B.; Ouyang, S.; Shi, Z.L. Evolutionary Arms Race between Virus and Host Drives Genetic Diversity in Bat Severe Acute Respiratory Syndrome-Related Coronavirus Spike Genes. *J. Virol.* **2020**, *94*, e00902-20.

85. Zhou, H.; Ji, J.; Chen, X.; Bi, Y.; Li, J.; Wang, Q.; Hu, T.; Song, H.; Zhao, R.; Chen, Y.; et al. Identification of novel bat coronaviruses sheds light on the evolutionary origins of SARS-CoV-2 and related viruses. *Cell* **2021**, *184*, 4380–4391.e14. [CrossRef]

86. Wang, L.; Fu, S.; Cao, Y.; Zhang, H.; Feng, Y.; Yang, W.; Nie, K.; Ma, X.; Liang, G. Discovery and genetic analysis of novel coronaviruses in least horseshoe bats in southwestern China. *Emerg. Microbes Infect.* **2017**, *6*, e14. [CrossRef]

87. Damas, J.; Hughes, G.M.; Keough, K.C.; Painter, C.A.; Persky, N.S.; Corbo, M.; Hiller, M.; Koepfli, K.P.; Pfenning, A.R.; Zhao, H.; et al. Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates. *Proc. Natl. Acad. Sci. USA* **2020**, *117*, 22311–22322. [CrossRef] [PubMed]

88. Lim, S.; Zhang, M.; Chang, T.L. ACE2-Independent Alternative Receptors for SARS-CoV-2. *Viruses* **2022**, *14*, 2535. [CrossRef]

89. Cantuti-Castelvetri, L.; Ojha, R.; Pedro, L.D.; Djannatian, M.; Franz, J.; Kuivanen, S.; van der Meer, F.; Kallio, K.; Kaya, T.; Anastasina, M.; et al. Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity. *Science* **2020**, *370*, 856–860. [CrossRef]

90. Hao, W.; Ma, B.; Li, Z.; Wang, X.; Gao, X.; Li, Y.; Qin, B.; Shang, S.; Cui, S.; Tan, Z. Binding of the SARS-CoV-2 spike protein to glycans. *Sci. Bull.* **2021**, *66*, 1205–1214. [CrossRef] [PubMed]

91. Li, Y.; Xu, Q.; Ma, L.; Wu, D.; Gao, J.; Chen, G.; Li, H. Systematic profiling of ACE2 expression in diverse physiological and pathological conditions for COVID-19/SARS-CoV-2. *J. Cell Mol. Med.* **2020**, *24*, 9478–9482. [CrossRef]

92. Schroeder, S.; Mache, C.; Kleine-Weber, H.; Corman, V.M.; Muth, D.; Richter, A.; Fatykhova, D.; Memish, Z.A.; Stanifer, M.L.; Boulant, S.; et al. Functional comparison of MERS-coronavirus lineages reveals increased replicative fitness of the recombinant lineage 5. *Nat. Commun.* **2021**, *12*, 5324. [CrossRef]
93. Terada, Y.; Matsui, N.; Noguchi, K.; Kuwata, R.; Shimoda, H.; Soma, T.; Mochizuki, M.; Maeda, K. Emergence of pathogenic coronaviruses in cats by homologous recombination between feline and canine coronaviruses. *PLoS ONE* **2014**, *9*, e106534.
[CrossRef]

94. Laude, H.; Godet, M.; Bernard, S.; Gelfi, J.; Duarte, M.; Delmas, B. Functional domains in the spike protein of transmissible gastroenteritis virus. *Adv. Exp. Med. Biol.* **1995**, *380*, 299–304.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.