Identification of Evolutionary Trajectories Shared across Human Betacoronaviruses

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, although the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), although a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1, and SARS-CoV-2), we developed a methodological pipeline to classify shared nonsynonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection and draw upon protein structure data to identify potential biological implications. We find 30 candidate mutations, from which 4 (codon sites 18121 [nsp14/residue 28], 21623 [spike/21], 21635 [spike/25], and 23948 [spike/796]; SARS-CoV-2 genome numbering) further display evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness

Coronavirus receptor switch explained from the stereochemistry of protein-carbohydrate interactions and a single mutation

Hemagglutinin-esterases (HEs) are bimodular envelope proteins of orthomyxoviruses, toroviruses, and coronaviruses with a carbohydrate-binding "lectin" domain appended to a receptor-destroying sialate-O-acetylesterase ("esterase"). In concert, these domains facilitate dynamic virion attachment to cell-surface sialoglycans. Most HEs (type I) target 9-O-acetylated sialic acids (9-O-Ac-Sias), but one group of coronaviruses switched to using 4-O-Ac-Sias instead (type II). This specificity shift required quasisynchronous adaptations in the Sia-binding sites of both lectin and esterase domains. Previously, a partially disordered crystal structure of a type II HE revealed how the shift in lectin ligand specificity was achieved. How the switch in esterase substrate specificity was realized remained unresolved, however. Here, we present a complete structure of a type II HE with a receptor analog in the catalytic site and identify the mutations underlying the 9-O- to 4-O-Ac-Sia substrate switch. We show that (i) common principles pertaining to the stereochemistry of protein-carbohydrate interactions were at the core of the transition in lectin ligand and esterase substrate specificity; (ii) in consequence, the switch in O-Ac-Sia specificity could be readily accomplished via convergent intramolecular coevolution with only modest architectural changes in lectin and esterase domains; and (iii) a single, inconspicuous Ala-to-Ser substitution in the catalytic site was key to the emergence of the type II HEs. Our findings provide fundamental insights into how proteins "see" sugars and how this affects protein and virus evolution

Human coronaviruses OC43 and HKU1 bind to 9-O-acetylated sialic acids via a conserved receptor-binding site in spike protein domain A

Human betacoronaviruses OC43 and HKU1 are endemic respiratory pathogens and, while related, originated from independent zoonotic introductions. OC43 is in fact a host-range variant of the species Betacoronavirus-1, and more closely related to bovine coronavirus (BCoV)-its presumptive ancestor-and porcine hemagglutinating encephalomyelitis virus (PHEV). The β1-coronaviruses (β1CoVs) and HKU1 employ glycan-based receptors carrying 9-O-acetylated sialic acid (9-O-Ac-Sia). Receptor binding is mediated by spike protein S, the main determinant of coronavirus host specificity. For BCoV, a crystal structure for the receptor-binding domain S1A is available and for HKU1 a cryoelectron microscopy structure of the complete S ectodomain. However, the location of the receptor-binding site (RBS), arguably the single-most important piece of information, is unknown. Here we solved the 3.0-Å crystal structure of PHEV S1A We then took a comparative structural analysis approach to map the β1CoV S RBS, using the general design of 9-O-Ac-Sia-binding sites as blueprint, backed-up by automated ligand docking, structure-guided mutagenesis of OC43, BCoV, and PHEV S1A, and infectivity assays with BCoV-S-pseudotyped vesicular stomatitis viruses. The RBS is not exclusive to OC43 and related animal viruses, but is apparently conserved and functional also in HKU1 S1A The binding affinity of the HKU1 S RBS toward short sialoglycans is significantly lower than that of OC43, which we attribute to differences in local architecture and accessibility, and which may be indicative for differences between the two viruses in receptor fine-specificity. Our findings challenge reports that would map the OC43 RBS elsewhere in S1A and that of HKU1 in domain S1B

Human coronaviruses OC43 and HKU1 bind to 9-O-acetylated sialic acids via a conserved receptor-binding site in spike protein domain A

Human betacoronaviruses OC43 and HKU1 are endemic respiratory pathogens and, while related, originated from independent zoonotic introductions. OC43 is in fact a host-range variant of the species Betacoronavirus-1, and more closely related to bovine coronavirus (BCoV)-its presumptive ancestor-and porcine hemagglutinating encephalomyelitis virus (PHEV). The β1-coronaviruses (β1CoVs) and HKU1 employ glycan-based receptors carrying 9-O-acetylated sialic acid (9-O-Ac-Sia). Receptor binding is mediated by spike protein S, the main determinant of coronavirus host specificity. For BCoV, a crystal structure for the receptor-binding domain S1A is available and for HKU1 a cryoelectron microscopy structure of the complete S ectodomain. However, the location of the receptor-binding site (RBS), arguably the single-most important piece of information, is unknown. Here we solved the 3.0-Å crystal structure of PHEV S1A We then took a comparative structural analysis approach to map the β1CoV S RBS, using the general design of 9-O-Ac-Sia-binding sites as blueprint, backed-up by automated ligand docking, structure-guided mutagenesis of OC43, BCoV, and PHEV S1A, and infectivity assays with BCoV-S-pseudotyped vesicular stomatitis viruses. The RBS is not exclusive to OC43 and related animal viruses, but is apparently conserved and functional also in HKU1 S1A The binding affinity of the HKU1 S RBS toward short sialoglycans is significantly lower than that of OC43, which we attribute to differences in local architecture and accessibility, and which may be indicative for differences between the two viruses in receptor fine-specificity. Our findings challenge reports that would map the OC43 RBS elsewhere in S1A and that of HKU1 in domain S1B

Broad receptor engagement of an emerging global coronavirus may potentiate its diverse cross-species transmissibility

Porcine deltacoronavirus (PDCoV), identified in 2012, is a common enteropathogen of swine with worldwide distribution. The source and evolutionary history of this virus is, however, unknown. PDCoV belongs to the Deltacoronavirus genus that comprises predominantly avian CoV. Phylogenetic analysis suggests that PDCoV originated relatively recently from a host-switching event between birds and mammals. Insight into receptor engagement by PDCoV may shed light into such an exceptional phenomenon. Here we report that PDCoV employs host aminopeptidase N (APN) as an entry receptor and interacts with APN via domain B of its spike (S) protein. Infection of porcine cells with PDCoV was drastically reduced by APN knockout and rescued after reconstitution of APN expression. In addition, we observed that PDCoV efficiently infects cells of unusual broad species range, including human and chicken. Accordingly, PDCoV S was found to target the phylogenetically conserved catalytic domain of APN. Moreover, transient expression of porcine, feline, human, and chicken APN renders cells susceptible to PDCoV infection. Binding of PDCoV to an interspecies conserved site on APN may facilitate direct transmission of PDCoV to nonreservoir species, including humans, potentially reflecting the mechanism that enabled a virus, ancestral to PDCoV, to breach the species barrier between birds and mammals. The APN cell surface protein is also used by several members of the Alphacoronavirus genus. Hence, our data constitute the second identification of CoVs from different genera that use the same receptor, implying that CoV receptor selection is subjected to specific restrictions that are still poorly understood

Betacoronavirus Adaptation to Humans Involved Progressive Loss of Hemagglutinin-Esterase Lectin Activity

Human beta1-coronavirus (β1CoV) OC43 emerged relatively recently through a single zoonotic introduction. Like related animal β1CoVs, OC43 uses 9-O-acetylated sialic acid as receptor determinant. β1CoV receptor binding is typically controlled by attachment/fusion spike protein S and receptor-binding/receptor-destroying hemagglutinin-esterase protein HE. We show that following OC43's introduction into humans, HE-mediated receptor binding was selected against and ultimately lost through progressive accumulation of mutations in the HE lectin domain. Consequently, virion-associated receptor-destroying activity toward multivalent glycoconjugates was reduced and altered such that some clustered receptor populations are no longer cleaved. Loss of HE lectin function was also observed for another respiratory human coronavirus, HKU1. This thus appears to be an adaptation to the sialoglycome of the human respiratory tract and for replication in human airways. The findings suggest that the dynamics of virion-glycan interactions contribute to host tropism. Our observations are relevant also to other human respiratory viruses of zoonotic origin, particularly influenza A virus

Betacoronavirus Adaptation to Humans Involved Progressive Loss of Hemagglutinin-Esterase Lectin Activity

Human beta1-coronavirus (β1CoV) OC43 emerged relatively recently through a single zoonotic introduction. Like related animal β1CoVs, OC43 uses 9-O-acetylated sialic acid as receptor determinant. β1CoV receptor binding is typically controlled by attachment/fusion spike protein S and receptor-binding/receptor-destroying hemagglutinin-esterase protein HE. We show that following OC43's introduction into humans, HE-mediated receptor binding was selected against and ultimately lost through progressive accumulation of mutations in the HE lectin domain. Consequently, virion-associated receptor-destroying activity toward multivalent glycoconjugates was reduced and altered such that some clustered receptor populations are no longer cleaved. Loss of HE lectin function was also observed for another respiratory human coronavirus, HKU1. This thus appears to be an adaptation to the sialoglycome of the human respiratory tract and for replication in human airways. The findings suggest that the dynamics of virion-glycan interactions contribute to host tropism. Our observations are relevant also to other human respiratory viruses of zoonotic origin, particularly influenza A virus

Broad receptor engagement of an emerging global coronavirus may potentiate its diverse cross-species transmissibility

Porcine deltacoronavirus (PDCoV), identified in 2012, is a common enteropathogen of swine with worldwide distribution. The source and evolutionary history of this virus is, however, unknown. PDCoV belongs to the Deltacoronavirus genus that comprises predominantly avian CoV. Phylogenetic analysis suggests that PDCoV originated relatively recently from a host-switching event between birds and mammals. Insight into receptor engagement by PDCoV may shed light into such an exceptional phenomenon. Here we report that PDCoV employs host aminopeptidase N (APN) as an entry receptor and interacts with APN via domain B of its spike (S) protein. Infection of porcine cells with PDCoV was drastically reduced by APN knockout and rescued after reconstitution of APN expression. In addition, we observed that PDCoV efficiently infects cells of unusual broad species range, including human and chicken. Accordingly, PDCoV S was found to target the phylogenetically conserved catalytic domain of APN. Moreover, transient expression of porcine, feline, human, and chicken APN renders cells susceptible to PDCoV infection. Binding of PDCoV to an interspecies conserved site on APN may facilitate direct transmission of PDCoV to nonreservoir species, including humans, potentially reflecting the mechanism that enabled a virus, ancestral to PDCoV, to breach the species barrier between birds and mammals. The APN cell surface protein is also used by several members of the Alphacoronavirus genus. Hence, our data constitute the second identification of CoVs from different genera that use the same receptor, implying that CoV receptor selection is subjected to specific restrictions that are still poorly understood

Multifunctional human monoclonal antibody combination mediates protection against Rift Valley fever virus at low doses

Abstract The zoonotic Rift Valley fever virus (RVFV) can cause severe disease in humans and has pandemic potential, yet no approved vaccine or therapy exists. Here we describe a dual-mechanism human monoclonal antibody (mAb) combination against RVFV that is effective at minimal doses in a lethal mouse model of infection. We structurally analyze and characterize the binding mode of a prototypical potent Gn domain-A-binding antibody that blocks attachment and of an antibody that inhibits infection by abrogating the fusion process as previously determined. Surprisingly, the Gn domain-A antibody does not directly block RVFV Gn interaction with the host receptor low density lipoprotein receptor-related protein 1 (LRP1) as determined by a competitive assay. This study identifies a rationally designed combination of human mAbs deserving of future investigation for use in humans against RVFV infection. Using a two-pronged mechanistic approach, we demonstrate the potent efficacy of a rationally designed combination mAb therapeutic