A novel SARS-CoV-2 related coronavirus with complex recombination isolated from bats in Yunnan province, China

ABSTRACT

At the end of 2019, A new type of beta-CoV, SARS-CoV-2 emerged and triggered the COVID-19 pandemic, which spread overwhelmingly around the world in less than a year. However, the origin and direct ancestral viruses of SARS-CoV-2 remain unknown. RaTG13, a novel coronavirus found in bats in China’s Yunnan Province, is the closest relative virus of the SARS-CoV-2 identified so far. In this study, a new SARS-CoV-2 related virus, provisionally named PrC31, was discovered in Yunnan province by retrospectively analyse bat next generation sequencing (NGS) data of intestinal samples collected in 2018. PrC31 shared 90.7% and 92.0% nucleotide identities to the genomes of SARS-CoV-2 and the bat SARSr-CoV ZC45, respectively. Sequence alignment of PrC31 showed that several genomic regions, especially orf1a and orf8 had the highest homology with those corresponding genomic regions of SARS-CoV-2 than any other related viruses. Phylogenetic analysis indicated that PrC31 shared a common ancestor with SARS-CoV-2 in evolutionary history. The differences between the PrC31 and SARS-CoV-2 genomes were mainly manifested in the spike genes. The amino acid homology between the receptor binding domains of PrC31 and SARS-CoV-2 was only 64.2%. Importantly, recombination analysis revealed that PrC31 underwent multiple complex recombination events (including three recombination breakpoints) involving the SARS-CoV and SARS-CoV-2 sub-lineages, indicating that PrC31 evolved from yet-to-be-identified intermediate recombination strains. Combined with previous studies, it is revealed that the beta-CoVs may possess a more complex recombination mechanism than we thought.

KEYWORDS:

• SARS-CoV-2
• bats
• coronaviruses
• recombination
• evolution

Disclosure statement

No potential conflict of interest was reported by the author(s).

Data availability statement

The sequence of PrC31 generated in this study was deposited in the GISAID and GenBank databases with the accession numbers EPI_ISL_1098866 and MW703458, respectively.

Additional information

Funding

This work was supported by National Science and Technology Major Project of China (2018ZX10305409-004-002), Basic Research Projects of Yunnan Province (2019FA015) and Yunnan Science and Technology Talents and Platform Plan (2018HB046).