Recent Zoonotic Spillover and Tropism Shift of a Canine Coronavirus Is Associated with Relaxed Selection and Putative Loss of Function in NTD Subdomain of Spike Protein

Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

Whole genome characteristics of hedgehog coronaviruses from Poland and analysis of the evolution of the Spike protein for its interspecies transmission potential

BACKGROUND

The hedgehogs have been recently identified as possible reservoir of Middle East respiratory syndrome coronavirus like (MERS-CoV-like). These viruses were classified as a distinct Betacoronavirus erinacei (BCoV-Eri) species within the MerBCoV-Eriirus subgenus. As coronaviruses are known for their ability to jump between different hosts, including humans, this can pose a particular threat to people in direct contact with hedgehogs, such as those working at animal asylums. Our previous studies have shown the presence of BCoV-Eri strains in animals collected in the wildlife rehabilitation centre. This study aimed to investigate the presence of CoV in subsequent hedgehogs collected from the urban area of Poland and their molecular characteristics.

RESULTS

Monitoring for the presence of coronavirus infection in hedgehogs revealed five positive individuals. The presence of BCoV-Eri was found in a total of 20% of animals tested. Our analyses revealed no correlation between CoVs positivity and animal health conditions but a higher probability of such infection in juveniles and females. The whole genome of two Polish Hedgehog coronavirus 1 strains were sequenced and compared with available counterparts from European and Asian countries. Phylogenetic analysis showed that both CoV strains formed common cluster with other similar MerBCoV-Eriirus, but they were also found to be genetically variable and most changes in the S protein were identified. Our analysis revealed that some S protein sites of the Hedgehog coronavirus 1 strains evolved under positive selection pressure and of five such sites, three are in the S1 region while the other two in the S2 region of the Spike.

CONCLUSIONS

BCoV-Eri is to some extent prevalent in wildlife asylums in Poland. Given that the S protein of BCoVs-Eri is highly variable and that some sites of this protein evolve under positive selection pressure, these strains could potentially acquire a favourable feature for cross-species transmission. Consequently, the threat to humans working in such asylums is particularly high. Adequate biosecurity safeguards, but also human awareness of such risks, are therefore essential.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, in saliva and feces of meso-carnivores in the northeastern United States

IMPORTANCE

Several coronaviruses (CoVs) have been detected in domesticated, farmed, and wild meso-carnivores, causing a wide range of diseases and infecting diverse species, highlighting their important but understudied role in the epidemiology of these viruses. Assessing the viral diversity hosted in wildlife species is essential to understand their significance in the cross-species transmission of CoVs. Our focus here was on CoV discovery in meso-carnivores in the Northeast United States as a potential "hotspot" area with high density of humans and urban wildlife. This study identifies novel alphacoronaviruses circulating in multiple free-ranging wild and domestic species in this area and explores their potential epidemiological importance based on regions of the Spike gene, which are relevant for virus-host interactions.

Characterization of CCoV-HuPn-2018 spike protein-mediated viral entry

Canine coronavirus-human pneumonia-2018 (CCoV-HuPn-2018) was recently isolated from a child with pneumonia. This novel human pathogen resulted from cross-species transmission of a canine coronavirus. It has been known that CCoV-HuPn-2018 uses aminopeptidase N (APN) from canines, felines, and porcines, but not humans, as functional receptors for cell entry. The molecular mechanism of cell entry in CCoV-HuPn-2018 remains poorly understood. In this study, we demonstrated that among the nine APN orthologs tested, the APN of the Mexican free-tailed bat could also efficiently support CCoV-HuPn-2018 spike (S) protein-mediated entry, raising the possibility that bats may also be an alternative host epidemiologically important for the transmission of this virus. The glycosylation at residue N747 of canine APN is critical for its receptor activity. The gain of glycosylation at the corresponding residues in human and rabbit APNs converted them to functional receptors for CCoV-HuPn-2018. Interestingly, the CCoV-HuPn-2018 spike protein pseudotyped virus infected multiple human cancer cell lines in a human APN-independent manner, whereas sialic acid appeared to facilitate the entry of the pseudotyped virus into human cancer cells. Moreover, while host cell surface proteases trypsin and TMPRSS2 did not promote the entry of CCoV-HuPn-2018, endosomal proteases cathepsin L and B are required for the entry of CCoV-HuPn-2018 in a pH-dependent manner. IFITMs and LY6E are host restriction factors for the CCoV-HuPn-2018 entry. Our results thus suggest that CCoV-HuPn-2018 has not yet evolved to be an efficient human pathogen. Collectively, this study helps us understand the cell tropism, receptor usage, cross-species transmission, natural reservoir, and pathogenesis of this potential human coronavirus. IMPORTANCE Viral entry is driven by the interaction between the viral spike protein and its specific cellular receptor, which determines cell tropism and host range and is the major constraint to interspecies transmission of coronaviruses. Aminopeptidase N (APN; also called CD13) is a cellular receptor for HCoV-229E, the newly discovered canine coronavirus-human pneumonia-2018 (CCoV-HuPn-2018), and many other animal alphacoronaviruses. We examined the receptor activity of nine APN orthologs and found that CCoV-HuPn-2018 utilizes APN from a broad range of animal species, including bats but not humans, to enter host cells. To our surprise, we found that CCoV-HuPn-2018 spike protein pseudotyped viral particles successfully infected multiple human hepatoma-derived cell lines and a lung cancer cell line, which is independent of the expression of human APN. Our findings thus provide mechanistic insight into the natural hosts and interspecies transmission of CCoV-HuPn-2018-like coronaviruses.

Canine respiratory coronavirus in Thailand undergoes mutation and evidences a potential putative parent for genetic recombination

Canine respiratory coronavirus (CRCoV) is associated with canine infectious respiratory disease complex. Although its detection has been reported worldwide, the genomic characteristics and evolutionary patterns of this virus remain poorly defined. In this study, 21 CRCoV sequences obtained from dogs in Thailand during two episodes (2013-2015, group A; 2021-2022, group B) were characterized and analyzed. The genomic characteristics of Thai CRCoVs changed from 2013 to 2022 and showed a distinct phylogenetic cluster. Phylogenetic analysis of the spike (S) genes divided the analyzed CRCoV strains into five clades. The full-length genome characterization revealed that all Thai CRCoVs possessed a nonsense mutation within the nonstructural gene located between the S and envelope genes, leading to a truncated putative nonstructural protein. Group B Thai CRCoV strains represented the signature nonsynonymous mutations in the S gene that was not identified in group A Thai CRCoVs, suggesting the ongoing evolutionary process of Thai CRCoVs. Although no evidence of recombination of Thai CRCoV strains was found, our analysis identified one Thai CRCoV strain as a potential parent virus for a CRCoV strain found in the United States. Selective pressure analysis of the hypervariable S region indicated that the CRCoV had undergone purifying selection during evolution. Evolutionary analysis suggested that the CRCoV was emerged in 1992 and was first introduced in Thailand in 2004, sharing a common ancestor with Korean CRCoV strains. These findings regarding the genetic characterization and evolutionary analysis of CRCoVs add to the understanding of CRCoVs. IMPORTANCE Knowledge of genomic characterization of the CRCoV is still limited and its evolution remains poorly investigated. We, therefore, investigated the full-length genome of CRCoV in Thailand for the first time and analyzed the evolutionary dynamic of CRCoV. Genomic characterization of Thai CRCoV strains revealed that they possess unique genome structures and have undergone nonsynonymous mutations, which have not been reported in previously described CRCoV strains. Our work suggests that the Thai CRCoVs were not undergone mutation through genetic recombination for their evolution. However, one Thai CRCoV strain PP158_THA_2015 was found to be a potential parent virus for the CRCoV strains found in the United States. This study provides an understanding of the genomic characterization and highlights the signature mutations and ongoing evolutionary process of CRCoV that could be crucial for monitoring in the future.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, shed in saliva and feces of meso-carnivores in the northeastern United States

Small to mid-sized carnivores, or meso-carnivores, comprise a group of diverse mammals, many of which can adapt to anthropogenically disturbed environments. Wild meso-carnivores living in urban areas may get exposed to or spread pathogens to other species, including stray/feral domestic animals. Several coronaviruses (CoVs) have been detected in domesticated and farmed meso-carnivores, but knowledge of CoVs circulating in free-ranging wild meso-carnivores remains limited. In this study, we analyzed 321 samples collected between 2016 and 2022 from 9 species of free-ranging wild meso-carnivores and stray/feral domestic cats in the northeastern United States. Using a pan-CoV PCR, we screened tissues, feces, and saliva, nasal, and rectal swabs. We detected CoV RNA in fecal and saliva samples of animals in four species: fisher (Pekania pennanti), bobcat (Lynx rufus), red fox (Vulpes vulpes), and domestic cat (Felis catus). Next-generation sequencing revealed that all these viruses belonged to the Luchacovirus subgenus (Alphacoronavirus genus), previously reported only in rodents and lagomorphs (i.e., rabbits). Genetic comparison of the 3'-end of the genome (~12,000bp) revealed that although the viruses detected group with, and have a genetic organization similar to other luchacoviruses, they are genetically distinct from those from rodents and lagomorphs. Genetic characterization of the spike protein revealed that the meso-carnivore luchacoviruses do not have an S1/S2 cleavage motif but do have highly variable structural loops containing cleavage motifs similar to those identified in certain pathogenic CoVs. This study highlights the importance of characterizing the spike protein of CoVs in wild species for further targeted epidemiologic monitoring.

Human neutralizing antibodies to cold linear epitopes and subdomain 1 of the SARS-CoV-2 spike glycoprotein

Emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants diminishes the efficacy of vaccines and antiviral monoclonal antibodies. Continued development of immunotherapies and vaccine immunogens resilient to viral evolution is therefore necessary. Using coldspot-guided antibody discovery, a screening approach that focuses on portions of the virus spike glycoprotein that are both functionally relevant and averse to change, we identified human neutralizing antibodies to highly conserved viral epitopes. Antibody fp.006 binds the fusion peptide and cross-reacts against coronaviruses of the four genera, including the nine human coronaviruses, through recognition of a conserved motif that includes the S2' site of proteolytic cleavage. Antibody hr2.016 targets the stem helix and neutralizes SARS-CoV-2 variants. Antibody sd1.040 binds to subdomain 1, synergizes with antibody rbd.042 for neutralization, and, similar to fp.006 and hr2.016, protects mice expressing human angiotensin-converting enzyme 2 against infection when present as a bispecific antibody. Thus, coldspot-guided antibody discovery reveals donor-derived neutralizing antibodies that are cross-reactive with Orthocoronavirinae, including SARS-CoV-2 variants.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Human neutralizing antibodies to cold linear epitopes and to subdomain 1 of SARS-CoV-2

Emergence of SARS-CoV-2 variants diminishes the efficacy of vaccines and antiviral monoclonal antibodies. Continued development of immunotherapies and vaccine immunogens resilient to viral evolution is therefore necessary. Using coldspot-guided antibody discovery, a screening approach that focuses on portions of the virus spike that are both functionally relevant and averse to change, we identified human neutralizing antibodies to highly conserved viral epitopes. Antibody fp.006 binds the fusion peptide and cross-reacts against coronaviruses of the four genera , including the nine human coronaviruses, through recognition of a conserved motif that includes the S2' site of proteolytic cleavage. Antibody hr2.016 targets the stem helix and neutralizes SARS-CoV-2 variants. Antibody sd1.040 binds to subdomain 1, synergizes with antibody rbd.042 for neutralization and, like fp.006 and hr2.016, protects mice when present as bispecific antibody. Thus, coldspot-guided antibody discovery reveals donor-derived neutralizing antibodies that are cross-reactive with Orthocoronavirinae , including SARS-CoV-2 variants.

One sentence summary

Broadly cross-reactive antibodies that protect from SARS-CoV-2 variants are revealed by virus coldspot-driven discovery.