Conserved recombination patterns across coronavirus subgenera

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination breakpoints at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination breakpoints across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination breakpoint hot-spot locations. We find that while the locations of recombination breakpoints are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination breakpoints most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination breakpoint distributions in coronavirus genomes sampled from nature.

A Prenylated dsRNA Sensor Protects Against Severe COVID-19 and is Absent in Horseshoe Bats

Cell autonomous antiviral defenses can inhibit the replication of viruses and reduce transmission and disease severity. To better understand the antiviral response to SARS-CoV-2, we used interferon-stimulated gene (ISG) expression screening to reveal that OAS1, through RNase L, potently inhibits SARS-CoV-2. We show that while some people can express a prenylated OAS1 variant, that is membrane-associated and blocks SARS-CoV-2 infection, other people express a cytosolic, nonprenylated OAS1 variant which does not detect SARS-CoV-2 (determined by the splice-acceptor SNP Rs10774671). Alleles encoding nonprenylated OAS1 predominate except in people of African descent. Importantly, in hospitalized patients, expression of prenylated OAS1 was associated with protection from severe COVID-19, suggesting this antiviral defense is a major component of a protective antiviral response. Remarkably, approximately 55 million years ago, retrotransposition ablated the OAS1 prenylation signal in horseshoe bats (the presumed source of SARS-CoV-2). Thus, SARS-CoV-2 never had to adapt to evade this defense. As prenylated OAS1 is widespread in animals, the billions of people that lack a prenylated OAS1 could make humans particularly vulnerable to the spillover of coronaviruses from horseshoe bats.

Exploring the natural origins of SARS-CoV-2

The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (Sarbecoviruses) have been identified, is fuelling speculation on the natural origins of SARS-CoV-2. We have comprehensively analysed phylogenetic relations between SARS-CoV-2, and the related bat and pangolin Sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Clustering of the inferred recombination events is non-random with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift in the ancestry of bat Sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm horseshoe bats, Rhinolophus, are the likely SARS-CoV-2 progenitor reservoir species. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and co-circulation of these viruses ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years or so. We confirm that a direct proximal ancestor to SARS-CoV-2 is yet to be sampled, since the closest relative shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for more wildlife sampling to (i) pinpoint the exact origins of SARS-CoV-2s animal progenitor, and (ii) survey the extent of the diversity in the related Sarbecoviruses phylogeny that present high risk for future spillover. HighlightsO_LIThe origin of SARS-CoV-2 can be traced to horseshoe bats, genus Rhinolophus, with ranges in both China and Southeast Asia. C_LIO_LIThe closest known relatives of SARS-CoV-2 exhibit frequent transmission among their Rhinolophus host species. C_LIO_LISarbecoviruses have undergone extensive recombination throughout their evolutionary history. C_LIO_LIAccounting for the mosaic patterns of these recombinants is important when inferring relatedness to SARS-CoV-2. C_LIO_LIBreakpoint patterns are consistent with recombination hotspots in the coronavirus genome, particularly upstream of the pike open reading frame with a coldspot in S1. C_LI

Recurrent emergence and transmission of a SARS-CoV-2 Spike deletion ΔH69/V70

SARS-CoV-2 amino acid replacements in the receptor binding domain (RBD) occur relatively frequently and some have a consequence for immune recognition. Here we report recurrent emergence and significant onward transmission of a six-nucleotide out of frame deletion in the S gene, which results in loss of two amino acids: H69 and V70. We report that in human infections {Delta}H69/V70 often co-occurs with the receptor binding motif amino acid replacements N501Y, N439K and Y453F, and in the latter two cases has followed the RBD mutation. One of the {Delta}H69/V70+ N501Y lineages, now known as B.1.1.7, has undergone rapid expansion and includes eight S gene mutations: RBD (N501Y and A570D), S1 ({Delta}H69/V70 and {Delta}144) and S2 (P681H, T716I, S982A and D1118H). In vitro, we show that {Delta}H69/V70 does not reduce serum neutralisation across multiple convalescent sera. However, {Delta}H69/V70 increases infectivity and is associated with increased incorporation of cleaved spike into virions. {Delta}H69/V70 is able to compensate for small infectivity defects induced by RBD mutations N501Y, N439K and Y453F. In addition, replacement of H69 and V70 residues in the B.1.1.7 spike reduces its infectivity and spike mediated cell-cell fusion. Based on our data {Delta}H69/V70 likely acts as a permissive mutation that allows acquisition of otherwise deleterious immune escape mutations. Enhanced surveillance for the {Delta}H69/V70 deletion with and without RBD mutations should be considered as a global priority not only as a marker for the B.1.1.7 variant, but potentially also for other emerging variants of concern. Vaccines designed to target the deleted spike protein could mitigate against its emergence as increased selective forces from immunity and vaccines increase globally. HighlightsO_LI{Delta}H69/V70 is present in at least 28 SARS-CoV-2 lineages C_LIO_LI{Delta}H69/V70 does not confer escape from convalescent sera C_LIO_LI{Delta}H69/V70 increases spike infectivity and compensates for RBD mutations C_LIO_LI{Delta}H69/V70 is associated with greater spike cleavage C_LIO_LIB.1.1.7 requires {Delta}H69/V70 for optimal spike cleavage and infectivity C_LI

The circulating SARS-CoV-2 spike variant N439K maintains fitness while evading antibody-mediated immunity

SARS-CoV-2 can mutate to evade immunity, with consequences for the efficacy of emerging vaccines and antibody therapeutics. Herein we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is the most divergent region of S, and provide epidemiological, clinical, and molecular characterization of a prevalent RBM variant, N439K. We demonstrate that N439K S protein has enhanced binding affinity to the hACE2 receptor, and that N439K virus has similar clinical outcomes and in vitro replication fitness as compared to wild- type. We observed that the N439K mutation resulted in immune escape from a panel of neutralizing monoclonal antibodies, including one in clinical trials, as well as from polyclonal sera from a sizeable fraction of persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Evidence of significant natural selection in the evolution of SARS-CoV-2 in bats, not humans

RNA viruses are proficient at switching host species, and evolving adaptations to exploit the new hosts cells efficiently. Surprisingly, SARS-CoV-2 has apparently required no significant adaptation to humans since the start of the COVID-19 pandemic, with no observed selective sweeps since genome sampling began. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects significant positive episodic diversifying selection acting on the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in ancestral hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor [~]1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. Collectively our results demonstrate the progenitor of SARS-CoV-2 was capable of near immediate human-human transmission as a consequence of its adaptive evolutionary history in bats, not humans.