key: cord-1006795-cd6qcntn
authors: Seyran, Murat; Hassan, Sk. Sarif; Uversky, Vladimir N.; Pal Choudhury, Pabitra; Uhal, Bruce D.; Lundstrom, Kenneth; Attrish, Diksha; Rezaei, Nima; Aljabali, Alaa A. A.; Ghosh, Shinjini; Pizzol, Damiano; Adadi, Parise; El-Aziz, Tarek Mohamed Abd; Kandimalla, Ramesh; Tambuwala, Murtaza M.; Lal, Amos; Azad, Gajendra Kumar; Sherchan, Samendra P.; Baetas-da-Cruz, Wagner; Palù, Giorgio; Brufsky, Adam M.
title: Urgent Need for Field Surveys of Coronaviruses in Southeast Asia to Understand the SARS-CoV-2 Phylogeny and Risk Assessment for Future Outbreaks †
date: 2021-03-09
journal: Biomolecules
DOI: 10.3390/biom11030398
sha: 715be8c504edc8f325c0ab1546e36c6564f9c3a8
doc_id: 1006795
cord_uid: cd6qcntn

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not adequately represent genetically related coronaviruses (CoVs) [...].

Biomolecules 2021, 11, 398 2 of 7 adequately represent genetically related coronaviruses (CoVs). The unique bat CoV RaTG13 is the only identified sequence genetically associated with SARS-CoV-2. Data scarcity of bat CoV sequences raises concerns over several fundamental experimental and biostatistical aspects, e.g., repeatability of sequences and intraspecies variations in critical gene regions, such as the receptor-binding domain of the spike protein. The Sunda pangolin has been proposed as the intermediate host and source of SARS-CoV-2, but no pangolin CoV isolates have been reported in its habitats in Southeast Asia. Most pangolin CoVs were isolated from pangolins captured during illegal animal trafficking, raising questions about such isolates' reliability and quality. Problems with pangolin CoV sampling are also related to the substandard quality of deposited sequences. There is an urgent need for field surveys of bat CoVs and possible intermediate hosts, such as pangolins, ferrets, and civets, in Southeast Asia to investigate the genomic source of SARS-CoV-2 and assess possible future risks for new outbreaks.

SARS-CoV-2 is the causative agent of the coronavirus disease 2019 (COVID-19) pandemic in which the first cases were reported in Wuhan, Hubei Province, China. SARS-CoV-2, a member of Betacoronavirus and subgenus Sarbecovirus, is phylogenetically related to bat coronaviruses (bat CoVs) RaTG13 (detected in Pu'er City in 2013) and RmYN02 (detected in Xishuangbanna City in 2019) and were detected approximately 2000 km from Wuhan [1, 2] . Betacoronavirus and its genetic reservoir bat species, such as the intermediate horseshoe bat Rhinolophus affinis, inhabit Southeast Asia [1, 3, 4] . It has been hypothesized that pangolin CoVs may have originated from cross-species transmission in bats [5] . To date, no additional bat CoVs RaTG13 or RmYN02 isolates and no pangolin CoV isolates originating from healthy or sick pangolins in their natural habitats in Southeast Asian countries have been sequenced ( Figure 1 ). Therefore, here, we seek to draw the attention of relevant institutions toward the need to obtain CoV isolates from potential hosts, such as healthy or sick pangolins, in their natural habitats to match the existing sequences.

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not ade quately represent genetically related coronaviruses (CoVs). The unique bat CoV RaTG1 is the only identified sequence genetically associated with SARS-CoV-2. Data scarcity o bat CoV sequences raises concerns over several fundamental experimental and biostatis tical aspects, e.g., repeatability of sequences and intraspecies variations in critical gen regions, such as the receptor-binding domain of the spike protein. The Sunda pangolin has been proposed as the intermediate host and source of SARS-CoV-2, but no pangolin CoV isolates have been reported in its habitats in Southeast Asia. Most pangolin CoV were isolated from pangolins captured during illegal animal trafficking, raising question about such isolates' reliability and quality. Problems with pangolin CoV sampling are als related to the substandard quality of deposited sequences. There is an urgent need fo field surveys of bat CoVs and possible intermediate hosts, such as pangolins, ferrets, and civets, in Southeast Asia to investigate the genomic source of SARS-CoV-2 and assess pos sible future risks for new outbreaks.

SARS-CoV-2 is the causative agent of the coronavirus disease 2019 (COVID-19) pan demic in which the first cases were reported in Wuhan, Hubei Province, China. SARS CoV-2, a member of Betacoronavirus and subgenus Sarbecovirus, is phylogenetically related to bat coronaviruses (bat CoVs) RaTG13 (detected in Pu'er City in 2013) and RmYN0 (detected in Xishuangbanna City in 2019) and were detected approximately 2000 km from Wuhan [1, 2] . Betacoronavirus and its genetic reservoir bat species, such as the intermediat horseshoe bat Rhinolophus affinis, inhabit Southeast Asia [1, 3, 4] . It has been hypothesized that pangolin CoVs may have originated from cross-species transmission in bats [5] . T date, no additional bat CoVs RaTG13 or RmYN02 isolates and no pangolin CoV isolate originating from healthy or sick pangolins in their natural habitats in Southeast Asian countries have been sequenced ( Figure 1 ). Therefore, here, we seek to draw the attention of relevant institutions toward the need to obtain CoV isolates from potential hosts, such as healthy or sick pangolins, in their natural habitats to match the existing sequences. Figure 1 . Zone 1 is the area where Sunda pangolins are distributed in Southeast Asia and the potential source area of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that has never been screened for genetically related clades of coronaviruses (CoVs). Zone 2 is the location where bat coronaviruses (bat CoVs) RaTG13 and RmYN02 were isolated in Yunnan, China. Zone is where Sunda pangolin CoVs were isolated from dead pangolins in Guangxi and Guangdong, China. Zone 4 is the location where SARS-CoV-2 was first reported from its epicenter in Wuhan [3] .

Based on genomic and sequence data submitted on pangolin CoVs in early 2020, pan golins were proposed as the intermediate host for SARS-CoV-2 [5] [6] [7] [8] . Recent genomi Figure 1 . Zone 1 is the area where Sunda pangolins are distributed in Southeast Asia and the potential source area of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that has never been screened for genetically related clades of coronaviruses (CoVs). Zone 2 is the location where bat coronaviruses (bat CoVs) RaTG13 and RmYN02 were isolated in Yunnan, China. Zone 3 is where Sunda pangolin CoVs were isolated from dead pangolins in Guangxi and Guangdong, China. Zone 4 is the location where SARS-CoV-2 was first reported from its epicenter in Wuhan [3] .

Based on genomic and sequence data submitted on pangolin CoVs in early 2020, pangolins were proposed as the intermediate host for SARS-CoV-2 [5] [6] [7] [8] . Recent genomic analysis suggests the possible recombination of bat CoVs and pangolin CoVs, which occurred at least twice, leading to the creation of SARS-CoV-2 [9] . Furthermore, pangolin cells lacking the interferon-induced helicase C domain 1 (IFIH1) and Z-DNA-binding protein (ZBP1) have been postulated to contribute to the switch from resistance to tolerance of CoV infections [10] . Pangolins such as the giant pangolin, Smutsia gigantea, inhabit the same caves as different bat species such as Hipposideridae sp., Emballonuridae sp., and Miniopterus sp. in Gabon, which could also be the case in Southeast Asia [11] .

Illegal animal trafficking of infected pangolins has been postulated as the possible transmission route of SARS-CoV-2 to its epicenter at Wuhan's wet seafood animal market [6] [7] [8] . The Sunda pangolin (Manis javanica) is native to Southeast Asian countries, such as Cambodia, Myanmar, Thailand, Laos, Vietnam, and a small area in Yunnan, China, where the potential source of the SARS-CoV-2 intermediate horseshoe bat is also present [1, 3, 12] . Under these conditions, the Sunda pangolin would serve as the intermediate host, and SARS-CoV-2 would have emerged through spillover from pangolins to humans once (scenario I) or several times (scenario II) [5] . It is also possible that SARS-CoV-2 and pangolin CoVs would have originated independently through the crosstransmission of bat species (scenario III) [5] . Interestingly, pangolin CoVs have also been detected approximately 1000 km from Wuhan (Guangxi Province, 2017, and Guangdong Province, 2019) [6] [7] [8] .

SARS-CoV-2 has been suggested to diverge from the lineage of bat CoV Ratg13 in 1969, with the highest posterior density interval of 95% for the years 1930 to 2000 [13] . To date, no isolates other than bat CoVs RaTG13 and RmYN02 have been discovered, resulting in incomplete phylogenetic analyses, such as the selection of midpoint rooting, limited taxon sampling, and inappropriate emphasis on one single element such as ACE2 [14] . The intraspecies phylogeny of the SARS-CoV-2 clinical isolate has been criticized, and it has been requested to utilize larger datasets and sample variability, which are also required for the determination of the interspecies SARS-CoV-2 phylogeny with bat CoVs, e.g., RaTG13 [14, 15] . Most phylogenetic studies ignore the necessity of sufficient intraspecies sample sizes and the assessment of intraspecies genomic variations [16] .

The findings from comparative phylogenetic analyses can be influenced by intraspecies sample sizes [16] , which might be the case for SARS-CoV-2. If the intraspecific genomic variation is high, a large dataset is necessary due to potential heterogeneity [16] . For instance, the pangolin CoVs and SARS-CoV-2 spike (S) protein receptor-binding domains (RBDs) are almost identical [8] . However, if the intraspecies variability in the S protein RBD of pangolin CoVs is high, the assessment of the genetic relatedness to SARS-CoV-2 is unreliable. This also applies to the phylogenetic analysis of SARS-CoV-2 with bat CoV species.

However, if the S protein RBD in the SARS-CoV-2 clade can be considered as conserved, the genomic relatedness to RmYN02 is irrelevant since their RBD homology is low [2] . Similarly, based on a single sequence of RmYN02, the furin protease cleavage insert of SARS-CoV-2 was proposed [2, 17] . If S1/S2 in the S protein of RmYN02 is a variable and a discrete trait of this viral species, these genetic conclusions are irrelevant [2] .

Reproducibility, seasonality, population differences, and measurement errors of sequences of bat CoVs and pangolin CoVs have not been validated [2, [6] [7] [8] 13] . Additionally, there could be sequencing errors due to several factors such as poor sample quality, improper handling, secondary PCR enrichment, and low-quality measurements [16, 18] . The RaTG13 strain was isolated in 2013, but its complete genomic sequence (GenBank ID MN996532) was submitted after the emergence of SARS-CoV-2 in 2020 [19] . Additionally, the RNA-dependent RNA polymerase (RdRp) gene of RaTG13 is identical to that of another bat CoV sequence BtCoV/4991 (GenBank ID KP876546) submitted in 2015 [19] . Strikingly, NCBI KRONA analysis of the RaTG13 sequence suggested the possibility of DNA contamination, and the sequence was considered a fossil record [7, 20] . For example, the analysis relied on the pangolin CoV sequence MP789 on GISAID (https: //www.gisaid.org/EPI_ISL_412860hCoV-19/pangolin/China/MP789/2019, accessed on 14 September 2020). However, the database commented on sequences with long stretches of unreadable bases in nucleic acid sequences annotated as NNNs (about 7%) and missing NSP14. Pangolin CoV MP789 was collected from the lungs of dead pangolins [7] . The sequencing was completed with gap-filling PCR, and the filled version was deposited in GenBank (MT121216.1) [7] . Although gap filling is a standard protocol, this indicates the need for better sample and sequence quality for pangolin CoVs. The gap filling had quality issues, missing data, and unexpected reads possibly due to contamination by other viruses, including SARS-CoV-2-related viruses, and mitochondrial genes from pangolins (NC_026781), humans (NC_012920), tigers (NC_010642), and mice (NC_005089) were identified [21] . Similarly, pangolin CoV isolate P3B, collected from Sunda pangolin blood, possesses long sequence stretches of unreadable NNNs (8% of the sequence), and NSP3, NSP2, NSP6, NSP4, NSP15, and NSP8 were missing in this sample. The sequence of P3B was characterized using a gap-filling protocol, but the completed sequence has not been deposited in GenBank [6] .

Based on three single isolates of bat CoVs (RaTG13, ZC45, and ZXC21) and gap-filled pangolin CoV (2017 and 2019 isolates) sequences, it was postulated that SARS-CoV-2 diversified as a species 70 years ago and remained undetected [12] . These analyzed sequences are an insufficient representation of all genetically related sarbecoviruses to SARS-CoV-2 in Southeast Asia [1, 3, 13] . This raises the question of how the wild-type SARS-CoV-2 has existed for 70 years without infecting humans. Many people involved in the hunting and trade of pangolins should have been exposed to the SARS-CoV-2 ancestor before the Wuhan outbreak [12] . For instance, approximately 83 bat species are consumed in 33 different countries, including 13 species in New Guinea and 14 species in the Philippines [22] . Considering the high mobility of bat species consumed by humans as bushmeats, there is a high possibility of direct contact of humans with SARS-CoV-2 [1] . Additionally, the higher risk of CoV transmission to other animals and humans is not due to meat consumption but contamination of water sources, since CoVs exist in bat feces [23] . For example, bats contaminated wells and ponds with MERS-CoV, which infected camels, leading to humans [23] .

However, the reason why SARS-CoV-2 remained dormant for 70 years could be related to its different pathology in bats. CoVs like RaTG13 infect the bat gastrointestinal system, since they have been detected in the feces and intestines of bats [23] . Bats cope well with CoVs, eliciting robust amounts of antibody responses [23] . In a study conducted in Yunnan Province on feces from the Chinese rufous horseshoe bat Rhinolophus sinicus, CoV was detected only in 12 of the 164 samples [4] . Therefore, bat-mediated host tropism events could have lesser proximity than expected due to the low frequency of bat CoV infection. Moreover, in a study conducted in Malaysia from 2009 to 2019 on 334 Sunda pangolins, none of the samples were positive for CoVs as verified by PCR [24] .

The SARS-CoV-2 host range has been evaluated based on the interaction with the ACE2 entry receptor [25] [26] [27] . Despite residues 24, 30, 34, 38, 82 , and 354 being common in both human and pangolin ACE2 protein, the binding affinity is much lower for pangolin ACE2, indicating that pangolins are not an intermediate host for the COVID-19 pandemic [26, 27] . Moreover, SARS-CoV-2 can use other entry pathways such as C-type lectin receptors (CLRs) and neuropilin [28] . Therefore, SARS-CoV-2 might be capable of infecting host species with an incompatible ACE2 protein structure, such as mice and chickens [25, 27] . Moreover, the presence of SARS-CoV-2 has been serologically confirmed in minks (Neovison vison), cats (Felis catus), ferrets (Mustela putorius furo), Chinese tree shrews (Tupaia belangeri chinensis), rhesus macaques (Macaca mulatta), domestic pigs (Sus domesticus), cynomolgus or crabeating macaques (Macaca fascicularis), racoon dogs (Nyctereutes procyonoides), common marmosets (Callithrix jacchus), hamsters (Mesocricetus auratus), African green or vervet monkeys (Chlorocebus aethiops), dogs (Canis familiaris), fruit bats (Rousettus aegyptiacus), tigers (Panthera tigris), and lions (Panthera leo), along with mild or moderate infections in some species [29] [30] [31] [32] [33] [34] [35] [36] [37] . Most cases could be traced back to human-to-animal transmission, but in some cases, fatal infections with an animal-to-animal transmission in minks and cats were detected, which indicates the highly contagious nature of SARS-CoV-2 [38] [39] [40] [41] . Even mink-to-human back-infections have been described [39] . Given the pathological potential of SARS-CoV-2, a field survey to acquire large complete-genome samples of CoVs for various inter-and intraspecies analyses is essential to investigate the zoonotic origin of SARS-CoV-2 or to discover genetic diversity/unity among various CoVs and especially sarbecoviruses in Southeast Asia [1, 2] . Furthermore, it will enrich information in the existing bat CoV and pangolin CoV species pool, leading to the discovery of CoVs in other potential hosts and expanding intraspecies sampling.

The foundation of humankind is vulnerable in the context of the pathological capacity of CoVs, especially SARS-CoV-2, a member of the Sarbecovirus subgenus that is responsible for the ongoing COVID-19 pandemic [42] . Expanded intra-and interspecies sampling for bat CoV RaTG13 and pangolin CoVs in Southeast Asia seems to be a reasonable scientific step critical for reliably elucidating the phylogenetic composition of the Sarbecovirus subgenus and its members, which are genetically related to SARS-CoV-2 [1, 2] . Multiple isolates from different locations, times, and hosts of bats, pangolins, and other potential populations of CoVs are essential for improving understanding of the adaptations, mutations, and recombination patterns of genetically related clades of SARS-CoV-2. Furthermore, if SARS-CoV-2 is unrelated to RaTG13, its true bat CoV ancestor is still awaiting to be discovered. Therefore, there is an urgent need for CoV surveys on bats and possible intermediate hosts in Southeast Asia to reliably investigate the SARS-CoV-2 ancestry and identify genetically related sarbecoviruses with a similar pathological capacity to prevent future CoV outbreaks [42] .

Possibility for reverse zoonotic transmission of SARS-CoV-2 to free-ranging wildlife: A case study of bats

A Novel Bat Coronavirus Closely Related to SARS-CoV-2 Contains Natural Insertions at the S1/S2 Cleavage Site of the Spike Protein

Global Epidemiology of Bat Coronaviruses

Identification of Diverse Bat Alphacoronaviruses and Betacoronaviruses in China Provides New Insights into the Evolution and Origin of Coronavirus-Related Diseases

Pangolins Harbor SARS-CoV-2-Related Coronaviruses

Identifying SARS-CoV-2-related coronaviruses in Malayan pangolins

PLoS Pathog. 2020

Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins

Insight into the origin of 5'UTR and source of CpG reduction in SARS-CoV-2 genome

Pangolins Lack IFIH1/MDA5, a Cytoplasmic RNA Sensor That Initiates Innate Immune Defense Upon Coronavirus Infection

Pangolins and bats living together in underground burrows in Lopé National Park

Covid-19: Natural or anthropic origin? Mammalia

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

Origins of SARS-CoV-1 and SARS-CoV-2 are often poorly explored in leading publications

Sampling bias and incorrect rooting make phylogenetic network tracing of SARS-COV-2 infections unreliable

Biomolecules 2021

Effects of sample size and intraspecific variation in phylogenetic comparative studies: A metaanalytic review

Questions concerning the proximal origin of SARS-CoV-2

Analysis of error profiles in deep next-generation sequencing data

The genetic structure of SARS-CoV-2 does not rule out a laboratory origin: SARS-COV-2 chimeric structure and furin cleavage site might be the result of genetic manipulation

De-novo assembly of RaTG13 Genome Reveals Inconsistencies further Obscuring SARS-CoV-2 Origins

The SARS-CoV-2-like virus found in captive pangolins from Guangdong should be better sequenced

Bats as bushmeat: A global review

Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS

No Evidence of Coronaviruses or Other Potentially Zoonotic Viruses in Sunda pangolins (Manis javanica) Entering the Wildlife Trade via Malaysia

Predicting susceptibility to SARS-CoV-2 infection based on structural differences in ACE2 across species

COVID-19: Time to exonerate the pangolin from the transmission of SARS-CoV-2 to humans

Broad host range of SARS-CoV-2 and the molecular basis for SARS-CoV-2 binding to cat ACE2

The Structural Basis of Accelerated Host Cell Entry by SARS-CoV-2

Susceptibility of raccoon dogs for experimental SARS-CoV-2 infection

High prevalence of SARS-CoV-2 antibodies in pets from COVID-19+ households

Animals and SARS-CoV-2: Species susceptibility and viral transmission in experimental and natural conditions, and the potential implications for community transmission

Susceptibility of swine cells and domestic pigs to SARS-CoV-2

From People to Panthera: Natural SARS-CoV-2 Infection in Tigers and Lions at the Bronx Zoo

Viral CpG Deficiency Provides No Evidence That Dogs Were Intermediate Hosts for SARS-CoV-2

SARS-CoV-2 in fruit bats, ferrets, pigs, and chickens: An experimental transmission study

Household pets and SARS-CoV2 transmissibility in the light of the ACE2 intrinsic disorder status

COVID-19-like symptoms observed in Chinese tree shrews infected with SARS-CoV-2

SARS-CoV-2 infection, disease and transmission in domestic cats

Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans

New SARS-CoV-2 Infection Detected in an Italian Pet Cat by RT-qPCR from Deep Pharyngeal Swab

3D reconstruction of SARS-CoV-2 infection in ferrets emphasizes focal infection pattern in the upper respiratory tract

Opinion: To stop the next pandemic, we need to unravel the origins of COVID-19

Funding: This research received no external funding.

The authors declare no conflict of interest.