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

HAL
open science

# Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats

Feng Zhu, Veasna Duong, Xiao Fang Lim, Vibol Hul, Tanu Chawla, Lucy Keatts, Tracey Goldstein, Alexandre Hassanin, Vuong Tan Tu, Philippe Buchy, et al.

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

► To cite this version:

Feng Zhu, Veasna Duong, Xiao Fang Lim, Vibol Hul, Tanu Chawla, et al., Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats. Viruses, 2022, 14 (2), pp.176. 10.3390/v14020176. pasteur-04096662

HAL Id: pasteur-04096662

https://riip.hal.science/pasteur-04096662

Submitted on 13 May 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_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Brief Report

# Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats

Feng Zhu<sup>1,†</sup> <u>D</u>, Veasna Duong<sup>2,†</sup> <u>D</u>, Xiao Fang Lim<sup>1,‡</sup>, Vibol Hul<sup>2,3</sup> <u>D</u>, Tanu Chawla<sup>1</sup>, Lucy Keatts<sup>4,5</sup> <u>D</u>, Tracey Goldstein<sup>5,§</sup>, Alexandre Hassanin<sup>6</sup> <u>D</u>, Vuong Tan Tu<sup>7</sup> <u>D</u>, Philippe Buchy<sup>2,||</sup>, October M. Sessions<sup>1,¶,**</sup> <u>D</u>, Lin-Fa Wang<sup>1,*</sup> <u>D</u>, Philippe Dussart<sup>2,††,*</sup> <u>D</u> and Danielle E. Anderson<sup>1,8,*</sup> <u>D</u>

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

check for
updates

**Citation:** Zhu, F.; Duong, V.; Lim, X.F.; Hul, V.; Chawla, T.; Keatts, L.; Goldstein, T.; Hassanin, A.; Tu, V.T.; Buchy, P.; et al. Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats. *Viruses* **2022**, *14*, 176. https://doi.org/10.3390/v14020176

Academic Editor: Stefano Aquaro

Received: 31 December 2021

Accepted: 15 January 2022

Published: 18 January 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 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 have been recognized as an exceptional viral reservoir, especially for coronaviruses. At least three bat zoonotic coronaviruses (SARS-CoV, MERS-CoV and SARS-CoV-2) have been shown to cause severe diseases in humans and it is expected more will emerge. One of the major features of CoVs is that they are all highly prone to recombination. An extreme example is the insertion of the P10 gene from reoviruses in the bat CoV GCCDC1, first discovered in *Rousettus leschenaultii* bats in China. Here, we report the detection of GCCDC1 in four different bat species (*Eonycteris spelaea*, *Cynopterus sphinx*, *Rhinolophus shameli* and *Rousettus* sp.) in Cambodia. This finding demonstrates a much broader geographic and bat species range for this virus and indicates common cross-species transmission. Interestingly, one of the bat samples showed a co-infection with an Alpha CoV most closely related to RsYN14, a virus recently discovered in the same genus (*Rhinolophus*) of bat in Yunnan, China, 2020. Taken together, our latest findings highlight the need to conduct active surveillance in bats to assess the risk of emerging CoVs, especially in Southeast Asia.

**Keywords:** bats; coronavirus; GCCDC1; zoonosis; recombination; co-infection; cross-species transmission

<footnote n="1">1 Programme in Emerging Infectious Diseases, Duke-NUS Medical School, Singapore 169857, Singapore; feng.zhu@duke-nus.edu.sg (F.Z.); lxiaofan@dso.org.sg (X.F.L.); tanu.chawla@duke-nus.edu.sg (T.C.); october.sessions@nus.edu.sg (O.M.S.)</footnote>

<footnote n="2">2 Virology Unit, Institut Pasteur du Cambodge, Pasteur Network, Phnom Penh 120210, Cambodia; dveasna@pasteur-kh.org (V.D.); hvibol@pasteur-kh.org (V.H.); philippe.x.buchy@gsk.com (P.B.)</footnote>

<footnote n="3">3 Unité des Virus Émergents, (UVÉ: Aix-Marseille Univ-IRD 190-INSERM 1207), 13005 Marseille, France</footnote>

<footnote n="4">4 Wildlife Conservation Society, Health Program, Bronx, NY 10460, USA; lkeattsconmed@gmail.com</footnote>

<footnote n="5">5 One Health Institute, School of Veterinary Medicine, University of California, Davis, CA 95616, USA; tgoldstein@ucdavis.edu</footnote>

<footnote n="6">6 Institut de Systématique, Évolution, Biodiversité, Sorbonne Université, MNHN, CNRS, EPHE, UA, 75005 Paris, France; alexandre.hassanin@mnhn.fr</footnote>

<footnote n="7">7 Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology, No. 18, Hoang Quoc Viet Road, Cau Giay District, Hanoi 10072, Vietnam; vttu@iebr.ac.vn</footnote>

<footnote n="8">8 Peter Doherty Institute for Infection and Immunity, Department of Microbiology and Immunology, University of Melbourne, Melbourne 3000, Australia</footnote>

<footnote>* Correspondence: linfa.wang@duke-nus.edu.sg (L.-F.W.); pdussart@pasteur.mg (P.D.); danielle.anderson@unimelb.edu.au (D.E.A.)</footnote>

<footnote>† These authors contributed equally to this work.</footnote>

<footnote>‡ Present address: DSO National Laboratories, Singapore 118225, Singapore.</footnote>

<footnote>§ Present address: University of Illinois at Urbana-Champaign, Chicago, IL 61820, USA.</footnote>

<footnote>|| Present address: GlaxoSmithKline Vaccines R&D Greater China & Intercontinental, Singapore 139234, Singapore.</footnote>

<footnote>¶ Present address: Saw Swee Hock School of Public Health, National University of Singapore and National University Health System, Singapore 117549, Singapore.</footnote>

<footnote>** Present address: Department of Pharmacy, National University of Singapore, Singapore 117559, Singapore.</footnote>

<footnote>†† Present address: Virology Unit, Institut Pasteur de Madagascar, Pasteur Network, Antananarivo 101, Madagascar.</footnote>
# 1. Introduction

Coronaviruses (CoV) have been known to infect humans for decades, but the severity of disease attributed to this group of viruses was only first realized when severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2002 [1,2]. Bats were found to be the reservoir host of SARS-CoV, with zoonotic transmission to humans occurring through palm civets, an intermediate host [3]. In 2012, a decade after the SARS-CoV outbreak, Middle East respiratory syndrome coronavirus (MERS-CoV) emerged and this time the intermediate host was dromedary camels [4]. The third and undoubtedly most devastating CoV outbreak began in 2019, when SARS coronavirus 2 (SARS-CoV-2) emerged in Wuhan, China [5,6]. SARS-CoV-2 is the causative agent of the COVID-19 pandemic, and there is a global race to end the outbreak that, at the end of 2021, has led to over 286.8 million infections and 5.4 million deaths [7].

CoVs are single-stranded, positive-sense, enveloped RNA viruses, contain the largest ssRNA genome (>30 kb) of all known RNA viruses, and belong to the family Coronaviridae. CoVs are taxonomically grouped into four genera: *Alphacoronavirus*, *Betacoronavirus*, *Gammacoronavirus*, and *Deltacoronavirus*. Alpha and betacoronaviruses are known to infect humans and of the species currently classified in these genera, 54% are bat viruses [8].

It is well documented that replication of large DNA genomes is less error prone due to proofreading capacity of the DNA polymerase. In contrast to DNA polymerases, RNA dependent RNA polymerases (RdRp) lack replication fidelity [9]. CoVs are no exception and are not only prone to mutation during viral replication, but recombination. Recombination primarily occurs in the spike region [10]. Virus diversity and evolution increase the risk of zoonotic transmission as well as impacting both diagnostics and treatment.

One extreme example of a recombination event is the cross-family recombination of a Rousettus bat coronavirus (most related to HKU9) with a non-enveloped dsRNA Pteropine orthoreovirus creating a new CoV, named Rousettus bat coronavirus GCCDC1 (RoBat-CoV-GCCDC1), which incorporated a part of the reovirus p10 gene in between N and NS7a. GCCDC1 was first discovered in bats in China [11] in 2016 and then in Singapore bats in 2020 [12]. Genes p10 and NS7c are specifically present in GCCDC1 strains from China and Singapore but not in the HKU9 strain from Hong Kong.

In this study, we performed a metagenomic study to analyze bat samples from different times and geographic locations in Cambodia to investigate the presence of RoBat-CoV-GCCDC1.

# 2. Materials and Methods

## 2.1. *Sample Collection*

The bat sampling was performed in 2010 by the Muséum national d'Histoire naturelle (MNHN; Paris, France) and the Institut Pasteur du Cambodge (IPC) in cooperation with Cambodian Government partners. The MNHN was mandated by UNESCO and the National Authority of Preah Vihear to conduct a mammal survey in northern Cambodia. Bats were captured using mist nets and harp traps in three provinces, Preah Vihear, Steung Treng and Ratanakiri, to compare bat diversity on the two sides of the Mekong River. Selected bats captured were humanely euthanized in full compliance with local ethical and legal guidelines. Bat species were identified to genus or species level following morphological criteria and DNA barcoding.

Oral and rectal swabs were collected and stored in viral transport medium solution (VTM; containing tryptose phosphate Broth 2.95%, 145 mM NaCl, 5% gelatin, 54 mM Amphotericin B, 106 U/L penicillin-streptomycin, 80 mg/L gentamycin [Sigma-Aldrich]). All specimens were immediately transferred into liquid nitrogen containers before being transported to the Institut Pasteur du Cambodge laboratory where they were stored at -80 °C prior to testing.

Virus isolation was attempted by passaging samples on Vero-E6 (African green monkey kidney), BHK (baby hamster kidney) and LLC-MK2 (monkey kidney) and bat primary cell lines.
## 2.2. RNA Extraction

Following the manufacturer's instructions, 750 µL TRIzol reagent (Invitrogen, Grand Island, NY, USA) was added to 250 µL of sample, followed by an addition of 200 µL chloroform and shaking for 15 s. The sample was centrifuged for 15 min at 12,000 rpm at 4 °C. The clear upper aqueous layer containing RNA, was transferred to a new 1.5 mL tube, and 0.5 mL of isopropanol per mL of initial TRIzol was added. Gentle mixing by inverting 5 times was performed before incubating for 10 min at room temperature. The sample was centrifuged for 15 min at 12,000 rpm at 4 °C. The supernatant was discarded, and the remaining pellet was resuspended in 1000 µL of 80% ethanol. The sample was again centrifuged for 5 min at 10,000 rpm at 4 °C. The ethanol was removed, RNA pellet was dried and dissolved in 100% ethanol and stored at -80 °C until shipment on dry ice to Duke-NUS Medical School for sequencing.

## 2.3. Preparation of Illumina DNA Libraries from RNA

Illumina libraries were generated from total RNA with NEBNext Ultra Directional RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA) and NEB-Next Multiplex Oligos for Illumina (New England Biolabs) according to the manufacturer's instructions with minor modifications. Briefly, 5 µL of total RNA was added to first-strand synthesis buffer and random primers followed by incubation at 94 °C for 5 min to generate RNA fragments larger than 500 nucleotides (nt). The double stranded cDNA was purified using Mag-Bind RxnPure Plus beads (Omega Bio-Tek, Norcross, GA, USA) and eluted in 60 µL nuclease-free water after first-strand and second strand cDNA synthesis. A library size between 400-600 nt was size-selected with Mag-Bind RxnPure Plus beads (Omega Bio-Tek) in a two-step selection by adding 35 µL and subsequently 15 µL of beads to the reaction. The library was eluted in 20 µL of nuclease-free water and amplified by PCR (20 cycles). Libraries were purified using the MinElute PCR Purification kit (Qiagen, Germantown, MD, USA) and eluted in 25 µL nuclease-free H₂O. Libraries were visualized on a 1.5% agarose gel and quantified in a Bioanalyzer High Sensitivity DNA Assay (Agilent, Santa Clara, CA, USA).

## 2.4. Virus Enrichment

Targeted CoV genome enrichment was performed using a custom designed library [13] of biotinylated, 120-mer xGen Lockdown Probes (Integrated DNA Technologies, Singapore). Prior to capture of viral sequences, 1 µL of xGen Universal Blocker- TS Mix (Integrated DNA Technologies) matched according to the library index, were added to 20 µL of library DNA. A total of 0.5 µL of 5 µg Cot-1 DNA (Invitrogen) was added to block binding of probes to non-viral regions of library fragments. Blocked libraries were ethanol precipitated and resuspended in 2.5 µL of nuclease-free H₂O, 3 µL Nimblegen hybridization solution and 7.5 µL Nimblegen 2× hybridization buffer (Roche, Mannheim, Germany) following a 10 min incubation at room temperature. The resuspended libraries were denatured at 95 °C for 10 min and cooled on ice. A total of 3 pmol of probe mix was added and hybridized to the libraries for 4 h at 65 °C. To capture virus-specific library fragments, 100 µL of magnetic M-270 streptavidin Dynabeads (Life Technologies, Waltham, MA, US) were added to the hybridization reaction and the mix was incubated for a further 45 min at 65 °C, with shaking at 2000 rpm on a ThermoMixer C (Eppendorf, Hamburg, Germany) shaker. Streptavidin beads were washed to remove unbound DNA using the SeqCap EZ hybridization and wash kit (Roche) according to the manufacturer's instructions. A post-capture PCR amplification with P1 and P2 primers (Illumina, San Diego, CA, USA) was performed using the following conditions: initial denaturation at 95 °C for 2 min; 20 cycles of denaturation for 20 s at 95 °C, annealing for 20 s at 65 °C and extension for 15 s at 72 °C; and a final extension step for 3 min at 72 °C. The enriched library was purified using Agencourt AMPure XP beads (Beckman Coulter Genomics, Indianapolis, IN, USA) and eluted in 10 µL nuclease-free H₂O, visualized on a 1.5% agarose gel and quantified in a Bioanalyzer High Sensitivity DNA Assay (Agilent).
## 2.5. Data Analysis

Raw reads were first imported into Geneious Prime (version 2022.0.1) for downstream analysis and trimmed of adapters with BBDuk (version 38.84). Bat host species was validated by mapping the reads to whole mitochondrial genome or specific mitochondrial genes such as CYTB or COX (depending on sequence availability) by minimap2 and identifying the best hit with highest identity and coverage. De novo assembly was conducted with clean reads by SPAdes (version 3.13.0, http://cab.spbu.ru/software/spades/, accessed on 30 December 2021) in Metagenome mode. The longest contig for each sample was then blasted against GCCDC1 reference strain 346 (NC_030886) and other GCCDC1 strains in China and Singapore to evaluate the completeness of the genome. The name RK074 was assigned to the best contig (30,161 nt) for GCCDC1. Each sample was then individually mapped to the reference RK074 genome using Geneious assembler. Annotation of RK074 was performed by comparing and transferring the annotation of China and Singapore GCCDC1 after nucleotide sequence alignment conducted by MAFFT in Geneious Prime software. Individual gene alignment was generated by Geneious alignment and used to plot the phylogeny tree by the maximum-likelihood method with the general-time-reversible (GTR) model and 1000 bootstrap replicates in PHYML 3.0 software. Phylogeny trees of different genes were exported into R, integrated with sample information and plotted by the R package "ggtree" (version 2.2.4). SNP analysis was performed in Geneious Prime and fraction of specific SNPs in different countries or host bat genus was calculated and plotted by the R package "ggseqlogo" (version 0.1). Four different bat species distribution data across the world were derived from spatial data on terrestrial mammals from the International Union for the Conservation of Nature (IUCN 2021. The IUCN Red List of Threatened Species. January 2019 [version 6.2]. https://www.iucnredlist.org; Downloaded on 18 May 2021) and exported into R (version 4.0.2). For the 4 specific species, the data were extracted with R package "tidyverse" (version 1.3.0). World map data were retrieved and cropped using the R package "rnaturalearth" (version 0.1.0) and "sf" (version 0.9-6). Distribution data were plotted on the cropped world map using r package "ggplot2" (3.3.2). Written permission from IUCN Red List was obtained for the publication of spatial data used in this study.

# 3. Results

## 3.1. Sample Collection, Virus Isolation, RNA Extraction and NGS Analysis

A total of 10 healthy bats from four different genera were captured in the province of Steung Treng and Ratanakiri, North Cambodia during December 2010 (Table 1, Figure 1, Table S1). Oral and rectal swabs were collected and subjected to virus isolation and RNA extraction. Virus isolation was attempted on Vero-E6 (African green monkey kidney), BHK (baby hamster kidney) and LLC-MK2 (monkey kidney) and bat primary cell lines. However, the virus isolation was not successful. RNA collected from swabs was sequenced by NGS and targeted CoV genome enrichment as previously described [13]. After mapping the de novo assembly of each individual sample against known beta coronaviruses, 100% (10 of 10) returned GCCDC1 as the closest hits. Several rounds of mapping and refining of the contigs of GCCDC1 were carried out and eight of ten near-complete genomes (<10% for gaps) were assembled with two partial genomes (RK072 and RK022) covering 55% and 82.7%, respectively (Table 2). The best quality genome with 100% coverage was named RK074. Each genome was compared against the GCCDC1 strains found in China and Singapore. All genomes from Cambodian bats shared highest sequence identities with GCCDC1 346 among GCCDC1 356 from China and GCCDC1 from Singapore. It is interesting to note that GCCDC1 was found in four different bat species in Cambodia, covering the species diversity previously found in the Chinese and Singaporean studies [11,12], making Cambodia the overlapping "center" region for GCCDC1 bat hosts.

Depending on the genomic regions used for phylogenetic analysis, the trees displayed slightly different topologies (Figure 2). However, it is clear that the CoV strains detected in
Cambodian bats, albeit from different bat species, are clustered more closely than those detected in China or Singapore.

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

**Figure 1.** Sample location of GCCDC1 with host bat distribution map in East and Southeast Asia. Spatial distribution of four bat species were retrieved from IUCN Red List of Threatened Species version 1.18 https://www.iucnredlist.org. Downloaded on 18 May 2021. The map was further zoomed to show the two provinces as GCCDC1 sample locations in Cambodia.
<figure><img src="image_7.png" /></figure>

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

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

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

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

<figure><img src="image_12.png" /><figcaption>**Figure 2.** Phylogenetic analysis of GCCDC1. Phylogenetic trees based on the (A) whole genome, (B) RdRp, (C) spike, (D) NS3, (E) p10, and (F) N gene sequences. Beta coronavirus HKU9 was used as the outgroup. Trees were generated using PhyML with general-time-reversible (GTR) substitution model and 1000 bootstrap replicates. Numbers above or below the branches are percentage bootstrap values for the associated nodes. The x axis represents the number of substitutions per site.</figcaption></figure>
<table><caption>Table 1. Bat sampling.</caption><thead><tr><th>Animal ID</th><th>Country</th><th>Province</th><th>Village</th><th>Family</th><th>Genus</th><th>Species</th><th>Type of Sample</th><th>Date Sample Taken</th></tr></thead><tbody><tr><td>PH162</td><td>Cambodia</td><td>Steung Treng</td><td>Chhab Village</td><td>Pteropodidae</td><td>Cynopterus</td><td>sp.</td><td>Rectal swab</td><td>5–6-December-2010</td></tr><tr><td>PH178</td><td>Cambodia</td><td>Steung Treng</td><td>Chhab Village</td><td>Pteropodidae</td><td>Rousettus</td><td>sp.</td><td>Oral swab</td><td>6–7-December-2010</td></tr><tr><td>PH201</td><td>Cambodia</td><td>Steung Treng</td><td>Chhab Village</td><td>Rhinolophidae</td><td>Rhinolophus</td><td>shameli</td><td>Rectal swab</td><td>6–7-December-2010</td></tr><tr><td>RK022</td><td>Cambodia</td><td>Ratanakiri</td><td>Pouy Village</td><td>Pteropodidae</td><td>Cynopterus</td><td>sp.</td><td>Rectal swab</td><td>14–15-December-2010</td></tr><tr><td>RK058</td><td>Cambodia</td><td>Ratanakiri</td><td>Ta Kouy Village</td><td>Pteropodidae</td><td>Eonycteris</td><td>spelaea</td><td>Rectal swab</td><td>15–16-December-2010</td></tr><tr><td>RK059</td><td>Cambodia</td><td>Ratanakiri</td><td>Ta Kouy Village</td><td>Pteropodidae</td><td>Cynopterus</td><td>sphinx</td><td>Rectal swab</td><td>15–16-December-2010</td></tr><tr><td>RK072</td><td>Cambodia</td><td>Ratanakiri</td><td>Veung Seing Village</td><td>Pteropodidae</td><td>Cynopterus</td><td>sp.</td><td>Rectal swab</td><td>17–18-December-2010</td></tr><tr><td>RK074</td><td>Cambodia</td><td>Ratanakiri</td><td>Veung Seing Village</td><td>Pteropodidae</td><td>Eonycteris</td><td>spelaea</td><td>Rectal swab</td><td>17–18-December-2010</td></tr><tr><td>RK075</td><td>Cambodia</td><td>Ratanakiri</td><td>Pan Porng Village</td><td>Pteropodidae</td><td>Cynopterus</td><td>sp.</td><td>Rectal swab</td><td>18–19-December-2010</td></tr><tr><td>RK091</td><td>Cambodia</td><td>Ratanakiri</td><td>Pan Porng Village</td><td>Pteropodidae</td><td>Eonycteris</td><td>spelaea</td><td>Rectal swab</td><td>18–19-December-2010</td></tr></tbody></table>

<table><caption>Table 2. Genomic features for each of the ten GCCDC1 strains characterization by NGS.</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">Length</th><th rowspan="2">x<sup>#</sup></th><th rowspan="2">Genome</th><th rowspan="2">Genome</th><th rowspan="2">Genome</th><th colspan="2">ORF1ab</th><th colspan="2">S</th><th colspan="2">NS3</th><th colspan="2">E</th><th colspan="2">M</th><th colspan="2">N</th><th colspan="2">p10</th><th colspan="2">NS7a</th><th colspan="2">NS7b</th><th colspan="2">NS7c</th></tr><tr><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th><th>nt</th><th>aa</th><th>x<sup>#</sup></th></tr></thead><tbody><tr><td>PH162</td><td>30101</td><td>90.9</td><td>90.49</td><td>92</td><td>93.81</td><td>90.9</td><td>90.5</td><td>92</td><td>93.8</td><td>94.8</td><td>97</td><td>97.5</td><td>97</td><td>100</td><td>27.9</td><td>8.31</td><td>3.9</td><td>86.7</td><td>84.8</td><td>87</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>66.3</td><td>57.1</td><td>56</td></tr><tr><td>PH178</td><td>30101</td><td>93</td><td>92.84</td><td>94.4</td><td>92.07</td><td>93</td><td>92.8</td><td>94</td><td>92.1</td><td>93.1</td><td>94</td><td>97.5</td><td>97.4</td><td>100</td><td>28.9</td><td>8.31</td><td>5.2</td><td>97.4</td><td>97.7</td><td>100</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>73.6</td><td>60.1</td><td>69</td></tr><tr><td>PH201</td><td>30134</td><td>94.2</td><td>94.8</td><td>95.6</td><td>93.53</td><td>94.2</td><td>94.8</td><td>96</td><td>93.5</td><td>95.2</td><td>96</td><td>97.6</td><td>97.4</td><td>100</td><td>69.2</td><td>47.3</td><td>75</td><td>97.3</td><td>97.7</td><td>100</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>97.3</td><td>96.8</td><td>98</td></tr><tr><td>RK022</td><td>30113</td><td>82.8</td><td>77.79</td><td>81.8</td><td>83.93</td><td>82.8</td><td>77.8</td><td>82</td><td>83.9</td><td>81.1</td><td>83</td><td>89.4</td><td>85.7</td><td>90</td><td>27.9</td><td>8.31</td><td>3.9</td><td>82.2</td><td>78</td><td>78</td><td>96.4</td><td>97.5</td><td>100</td><td>97.5</td><td>96.8</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>78.1</td><td>68.9</td><td>75</td></tr><tr><td>RK058</td><td>30139</td><td>97.4</td><td>99.12</td><td>100</td><td>96.67</td><td>97.4</td><td>99.1</td><td>100</td><td>96.7</td><td>99.4</td><td>100</td><td>97.8</td><td>98.3</td><td>100</td><td>100</td><td>100</td><td>100</td><td>97.4</td><td>97.7</td><td>100</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>86.2</td><td>83.3</td><td>82</td></tr><tr><td>RK059</td><td>30101</td><td>93.8</td><td>94.31</td><td>95.9</td><td>91.01</td><td>93.8</td><td>94.3</td><td>96</td><td>91</td><td>91.7</td><td>93</td><td>92.8</td><td>91.9</td><td>94</td><td>25</td><td>4.55</td><td>0</td><td>92</td><td>90.5</td><td>91</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>97.5</td><td>96.1</td><td>98</td></tr><tr><td>RK072</td><td>29961</td><td>63.1</td><td>51.22</td><td>55.4</td><td>51.64</td><td>63.1</td><td>51.2</td><td>55</td><td>51.6</td><td>36.7</td><td>41</td><td>53.7</td><td>39.6</td><td>42</td><td>25</td><td>4.55</td><td>0</td><td>82.6</td><td>70.1</td><td>86</td><td>64.4</td><td>52.5</td><td>60</td><td>94.8</td><td>89.3</td><td>100</td><td>93.8</td><td>93.8</td><td>96</td><td>79.8</td><td>87.7</td><td>92</td><td>50.3</td><td>35.7</td><td>39</td></tr><tr><td>RK074</td><td>30145</td><td>97.6</td><td>99.43</td><td>100</td><td>96.45</td><td>97.6</td><td>99.4</td><td>100</td><td>96.5</td><td>99.2</td><td>100</td><td>97.2</td><td>97.4</td><td>100</td><td>100</td><td>100</td><td>100</td><td>97.3</td><td>97.7</td><td>100</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>99.6</td><td>100</td><td>100</td></tr><tr><td>RK075</td><td>30147</td><td>94.9</td><td>96.2</td><td>97.2</td><td>95.29</td><td>94.9</td><td>96.2</td><td>97</td><td>95.3</td><td>97.4</td><td>98</td><td>97.3</td><td>97</td><td>100</td><td>28.9</td><td>8.31</td><td>5.2</td><td>97.3</td><td>97.7</td><td>100</td><td>96.5</td><td>97.7</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>97.9</td><td>98.5</td><td>100</td><td>85.4</td><td>82.1</td><td>81</td></tr><tr><td>RK091</td><td>30147</td><td>94.3</td><td>94</td><td>96.1</td><td>94.04</td><td>94.3</td><td>94</td><td>96</td><td>94</td><td>95.1</td><td>98</td><td>96.1</td><td>92</td><td>100</td><td>98.3</td><td>100</td><td>100</td><td>96.4</td><td>96.8</td><td>100</td><td>98.1</td><td>99.6</td><td>100</td><td>98.2</td><td>98.9</td><td>100</td><td>97.3</td><td>100</td><td>100</td><td>99.5</td><td>98.8</td><td>100</td><td>52.3</td><td>39.8</td><td>36</td></tr></tbody></table>

<footnote>#: X represents coverage, shown in percentage.</footnote>
## 3.2. Bat Species-Specific Features

To investigate the bat species-specific features of GCCDC1, we performed whole-genome SNP analysis for all previous and current GCCDC1 strains in this study. A total of 33 sites in ORF1ab, Spike and NS3 genes were analyzed to compare the frequencies of non-synonymous mutations between different countries and host bat genus (Figure 3, Table S2). These 33 sites were chosen as they were non-synonymous substitutions occurring more than once in the examined sequences. There were no unique SNPs in GCCDC1 Chinese or Singaporean strains, while Cambodian strains had many. Both geographic region and host bat genus favoured some SNPs over others. For example, at Position 4, while the proportion of a T was quite similar in *Cynopterus*, *Eonycteris* and *Rousettus* host bats, there were none in China or Singapore strains. At Position 7, C was a hallmark for Cambodian strains, but not in the *Cynopterus* host. The *Cynopterus* and *Rhinolophus* host carried more unique SNPs than the *Eonycteris* and *Rousettus* strains.

<figure><img src="image_13.png" /><figcaption>**Figure 3.** Analysis on non-synonymous mutations in GCCDC1 in countries and host bat genus. A total of 33 sites were selected to calculate the frequencies of non-synonymous mutations in (A) different countries and (B) host bat genus. Locations of the sites and effect on protein sequence can be found on Sup Table 1. Due to mixture of reads in some samples, non-ATGC nucleotides were transferred into A or T or G or C in equal fractions.</figcaption></figure>

## 3.3. RsYN14-Like AlphaCoV in PH201

For bat sample PH201, an alphaCoV was found to be present in addition to GCCDC1, indicating a co-infection of two CoVs in two different genera. During sequence analysis and after removing the reads mapped to GCCDC1, the remaining reads could also be assembled into a near full genome (96.7% coverage) of a novel alphaCoV. The closest virus is RsYN14 which was discovered in the same genus of bat in Yunnan, China, 2020 [14]. The nucleotide sequence identity between the full genome PH201_AlphaCoV detected in Cambodia and RsYN14 is 83.2% (Table 3). Closer investigation at the gene level revealed
that PH201_AlphaCoV shares an average of 89.5% protein sequence identity with RsYN14 for all genes, ranging from ORF4a at 76.3% to M at 98.4%. Phylogenetic analysis based on the whole genome or individual genes (Figure 4) inferred the closeness of RsYN14 and PH201_AlphaCoV, with the second-closest virus being either *Rousettus* coronavirus HKU10 or *Rhinolophus* coronavirus Bt-CoV/Rh/YN2012, depending on the genes.

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

**Figure 4.** Phylogenetic analysis of PH201 AlphaCoV. Phylogenetic trees were based on the (A) whole genome or (B) RdRp, (C) spike and (D) N gene sequences. Beta coronavirus HKU9 was used as the outgroup. Trees were generated using PhyML with general-time-reversible (GTR) substitution model and 1000 bootstrap replicates. Numbers above or below the branches are percentage bootstrap values for the associated nodes. The x bar represents the number of substitutions per site. The PH201 AlphaCoV is marked with a star symbol.
<table><caption>**Table 3.** Genomic characteristics of AlphaCoV in PH201 bat.</caption><thead><tr><th>ID</th><th>Length</th><th>X#</th><th>MZ08395<br>(RsYN14)</th><th>nt</th><th>ORF1ab<br>aa</th><th>x#</th><th>nt</th><th>S<br>aa</th><th>x#</th><th>nt</th><th>ORF3a<br>aa</th><th>x#</th><th>nt</th><th>ORF4a<br>aa</th><th>x#</th><th>nt</th><th>E<br>aa</th><th>x#</th><th>nt</th><th>M<br>aa</th><th>x#</th><th>nt</th><th>N<br>aa</th><th>x#</th><th>nt</th><th>ORF8<br>aa</th><th>x#</th><th>nt</th><th>ORF9<br>aa</th><th>x#</th></tr></thead><tbody><tr><td>PH201</td><td>28807</td><td>96.7</td><td>83.23</td><td>82.67</td><td>90.37</td><td>96.9</td><td>80.36</td><td>91.09</td><td>96.9</td><td>85.67</td><td>90.61</td><td>100</td><td>77.94(early<br>termination)</td><td>76.29(early<br>termination)</td><td>100</td><td>90.89</td><td>93.55</td><td>93.8</td><td>95.22</td><td>98.35</td><td>100</td><td>86.75</td><td>91.48</td><td>100</td><td>89.89</td><td>90.91</td><td>100</td><td>83.77</td><td>80.79</td><td>100</td></tr></tbody></table>

<footnote>#: X represents coverage, shown in percentage.</footnote>
**4. Discussion**

Better understanding of potentially zoonotic viruses, such as bat CoVs, is highly important for risk assessment of potential spillover events which may lead to human-to-human transmission. Due to the high frequency of genome-level recombination, CoV evolution is dynamic and evolving virus variants poses a risk to public health.

Here we report the presence of multiple strains of the cross-family recombinant RoBat-CoV GCCDC1 in Cambodian bats. The authenticity of GCCDC1 found in Cambodia is confirmed by the close genetic relationship with those found in China and Singapore, as well as the presence and sequence conservation of the reovirus p10 gene located between N and NS7a.

This study highlights several significant observations. First, it indicates a wider geographic distribution of GCCDC1 in Asia than previously known. It will be interesting to investigate whether GCCDC1 is present in bats outside of Asia. Second, our data suggest cross-bat species infection by GCCDC1 can happen quite frequently as bats from all four genera sampled in this study were positive for GCCDC1. Third, SNP analysis indicated some bat species-specific associations could possibly be confounded by bat spatial distribution. Fourth, the detection of GCCDC1 in Cambodia in all bat species previously found in China or Singapore underscore not only the "central" location of Cambodia, but also the richness of both bat and virus diversity in Cambodia. Finally, the detection of a novel alphaCoV co-infected with GCCDC1 further highlights the propensity for co-evolution or spill-over of CoVs among different bat species. A recently discovered SARS-CoV-2 related bat coronavirus in Cambodia [15] was found in the same area, in the same host species (*Rhinolophus shameli*) and at the same time of sampling, suggesting the urgency of more detailed investigation and monitoring of co-evolution.

In conclusion, this study demonstrates the wide geographic location of an evolutionally dynamic recombinant bat CoV. Enhanced surveillance in a broad geographical area, in both bats and humans is required to monitor the evolution of this group of viruses and to assess potential zoonotic transmission into human populations.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14020176/s1, Table S1: bat species identification; Table S2: selected SNP sites for analysis.

**Author Contributions:** Collection of samples and lab testing in Cambodia V.H., V.D., A.H., V.T.T. Sample collection coordinator—L.K. Project coordination—P.B., L.K., T.G. Sampling and training on bat capture—A.H. and V.T.T. Project coordination in Singapore—L.-F.W., D.E.A., O.M.S. Lab testing and sequencing—D.E.A., X.F.L., T.C. Data analysis—F.Z. Preparation of manuscript—F.Z., V.D., D.E.A., P.D., with input from all authors. All authors have read and agreed to the published version of the manuscript.

**Funding:** DEA and L-FW were supported by Grants NRF2012NRFCRP001-056 and NRF2016NRF-NSFC002-013 from the Singapore National Research Foundation. This study was made possible by the generous support of the American people through the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT project (cooperative agreement number GHN-A-OO-09-00010-00 and AID-OAA-A-14-00102). The contents of this paper are the responsibility of the authors and do not necessarily reflect the views of the US Agency for International Development or the US Government. The fieldwork was supported in part by the National Authority for Preah Vihear, UNESCO, "Société des amis du Muséum et du Jardin des Plantes", and the Muséum national d'Histoire naturelle.

**Institutional Review Board Statement:** Field studies were carried out in compliance with the AR-RIVE guidelines.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Genomic sequences reported in this study have been deposited in GenBank under the accession number OM219639 to OM219649. Raw sequence reads reported in this study have been uploaded into SRA under BioProject ID PRJNA795576.
**Acknowledgments:** The investigators thank everyone involved in zoonotic virus surveillance and investigation in the Kingdom of Cambodia, especially the Ministry of Agriculture, Forestry and Fisheries and the Government of Cambodia for permission to conduct this work. We greatly thank the field team and diagnostic team in the Virology Unit, Institut Pasteur du Cambodge and laboratory members at Duke-NUS Medical School for all their hard work and dedication.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Ksiazek, T.G.; Erdman, D.; Goldsmith, C.S.; Zaki, S.R.; Peret, T.; Emery, S.; Tong, S.; Urbani, C.; Comer, J.A.; Lim, W.; et al. A novel coronavirus associated with severe acute respiratory syndrome. *N. Engl. J. Med.* **2003**, *348*, 1953–1966. [CrossRef] [PubMed]

2. Rota, P.A.; Oberste, M.S.; Monroe, S.S.; Nix, W.A.; Campagnoli, R.; Icenogle, J.P.; Penaranda, S.; Bankamp, B.; Maher, K.; Chen, M.H.; et al. Characterization of a novel coronavirus associated with severe acute respiratory syndrome. *Science* **2003**, *300*, 1394–1399. [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] [PubMed]

4. Haagmans, B.L.; Al Dhahiry, S.H.; Reusken, C.B.; Raj, V.S.; Galiano, M.; Myers, R.; Godeke, G.J.; Jonges, M.; Farag, E.; Diab, A.; et al. Middle East respiratory syndrome coronavirus in dromedary camels: An outbreak investigation. *Lancet Infect. Dis.* **2014**, *14*, 140–145. [CrossRef]

5. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. *Nature* **2020**, *579*, 270–273. [CrossRef] [PubMed]

6. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.M.; Lau, E.H.Y.; Wong, J.Y.; et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. *N. Engl. J. Med.* **2020**, *382*, 1199–1207. [CrossRef] [PubMed]

7. Worldometer. COVID-19 Coronavirus Pandemic. Available online: https://www.worldometers.info/coronavirus/ (accessed on 31 December 2021).

8. Walker, P.J.; Siddell, S.G.; Lefkowitz, E.J.; Mushegian, A.R.; Adriaenssens, E.M.; Alfenas-Zerbini, P.; Davison, A.J.; Dempsey, D.M.; Dutilh, B.E.; García, M.L.; et al. Changes to virus taxonomy and to the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses (2021). *Arch. Virol.* **2021**, *166*, 2633–2648. [CrossRef] [PubMed]

9. Duffy, S.; Shackelton, L.A.; Holmes, E.C. Rates of evolutionary change in viruses: Patterns and determinants. *Nat. Rev. Genet.* **2008**, *9*, 267–276. [CrossRef] [PubMed]

10. Bobay, L.-M.; O'Donnell, A.C.; Ochman, H. Recombination events are concentrated in the spike protein region of Betacoronaviruses. *PLoS Genet.* **2020**, *16*, e1009272. [CrossRef] [PubMed]

11. Huang, C.; Liu, W.J.; Xu, W.; Jin, T.; Zhao, Y.; Song, J.; Shi, Y.; Ji, W.; Jia, H.; Zhou, Y.; et al. A Bat-Derived Putative Cross-Family Recombinant Coronavirus with a Reovirus Gene. *PLoS Pathog.* **2016**, *12*, e1005883. [CrossRef] [PubMed]

12. Paskey, A.C.; Ng, J.H.J.; Rice, G.K.; Chia, W.N.; Philipson, C.W.; Foo, R.J.H.; Cer, R.Z.; Long, K.A.; Lueder, M.R.; Lim, X.F.; et al. Detection of Recombinant Rousettus Bat Coronavirus GCCDC1 in Lesser Dawn Bats (Eonycteris spelaea) in Singapore. *Viruses* **2020**, *12*, 539. [CrossRef]

13. Lim, X.F.; Lee, C.B.; Pascoe, S.M.; How, C.B.; Chan, S.; Tan, J.H.; Yang, X.; Zhou, P.; Shi, Z.; Sessions, O.M.; et al. Detection and characterization of a novel bat-borne coronavirus in Singapore using multiple molecular approaches. *J. Gen. Virol.* **2019**, *100*, 1363–1374. [CrossRef] [PubMed]

14. 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] [PubMed]

15. Delaune, D.; Hul, V.; Karlsson, E.A.; Hassanin, A.; Ou, T.P.; Baidaliuk, A.; Gambaro, 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]