RESUMEN
There are outstanding evolutionary questions on the recent emergence of coronavirus SARS-CoV-2/hCoV-19 in Hubei province that caused the COVID-19 pandemic, including (1) the relationship of the new virus to the SARS-related coronaviruses, (2) the role of bats as a reservoir species, (3) the potential role of other mammals in the emergence event, and (4) the role of recombination in viral emergence. Here, we address these questions and find that the sarbecoviruses - the viral subgenus responsible for the emergence of SARS-CoV and SARS-CoV-2 - exhibit frequent recombination, but the SARS-CoV-2 lineage itself is not a recombinant of any viruses detected to date. In order to employ phylogenetic methods to date the divergence events between SARS-CoV-2 and the bat sarbecovirus reservoir, recombinant regions of a 68-genome sarbecovirus alignment were removed with three independent methods. Bayesian evolutionary rate and divergence date estimates were consistent for all three recombination-free alignments and robust to two different prior specifications based on HCoV-OC43 and MERS-CoV evolutionary rates. Divergence dates between SARS-CoV-2 and the bat sarbecovirus reservoir were estimated as 1948 (95% HPD: 1879-1999), 1969 (95% HPD: 1930-2000), and 1982 (95% HPD: 1948-2009). Despite intensified characterization of sarbecoviruses since SARS, the lineage giving rise to SARS-CoV-2 has been circulating unnoticed for decades in bats and been transmitted to other hosts such as pangolins. The occurrence of a third significant coronavirus emergence in 17 years together with the high prevalence and virus diversity in bats implies that these viruses are likely to cross species boundaries again. In BriefThe Betacoronavirus SARS-CoV-2 is a member of the sarbecovirus subgenus which shows frequent recombination in its evolutionary history. We characterize the extent of this genetic exchange and identify non-recombining regions of the sarbecovirus genome using three independent methods to remove the effects of recombination. Using these non-recombining genome regions and prior information on coronavirus evolutionary rates, we obtain estimates from three approaches that the most likely divergence date of SARS-CoV-2 from its most closely related available bat sequences ranges from 1948 to 1982. Key PointsO_LIRaTG13 is the closest available bat virus to SARS-CoV-2; a sub-lineage of these bat viruses is able to infect humans. Two sister lineages of the RaTG13/SARS-CoV-2 lineage infect Malayan pangolins. C_LIO_LIThe sarbecoviruses show a pattern of deep recombination events, indicating that there are high levels of co-infection in horseshoe bats and that the viral pool can generate novel allele combinations and substantial genetic diversity; the sarbecoviruses are efficient explorers of phenotype space. C_LIO_LIThe SARS-CoV-2 lineage is not a recent recombinant, at least not involving any of the bat or pangolin viruses sampled to date. C_LIO_LINon-recombinant regions of the sarbecoviruses can be identified, allowing for phylogenetic inference and dating to be performed. We constructed three such regions using different methods. C_LIO_LIWe estimate that RaTG13 and SARS-CoV-2 diverged 40 to 70 years ago. There is a diverse unsampled reservoir of generalist viruses established in horseshoe bats. C_LIO_LIWhile an intermediate host responsible for the zoonotic event cannot be ruled out, the relevant evolution for spillover to humans very likely occurred in horseshoe bats. C_LI
RESUMEN
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.
RESUMEN
As three SARS-CoV-2 vaccines come to market in Europe and North America in the winter of 2020-2021, distribution networks will be in a race against a major epidemiological wave of SARS-CoV-2 that began in autumn 2020. Rapid and optimized vaccine allocation is critical during this time. With 95% efficacy reported for two of the vaccines, near-term public health needs require that distribution is prioritized to the elderly, health-care workers, teachers, essential workers, and individuals with co-morbidities putting them at risk of severe clinical progression. Here, we evaluate various age-based vaccine distributions using a validated mathematical model based on current epidemic trends in Rhode Island and Massachusetts. We allow for varying waning efficacy of vaccine-induced immunity, as this has not yet been measured. We account for the fact that known COVID-positive cases may not be included in the first round of vaccination. And, we account for current age-specific immune patterns in both states. We find that allocating a substantial proportion (> 75%) of vaccine supply to individuals over the age of 70 is optimal in terms of reducing total cumulative deaths through mid-2021. As we do not explicitly model other high mortality groups, this result on vaccine allocation applies to all groups at high risk of mortality if infected. Our analysis confirms that for an easily transmissible respiratory virus, allocating a large majority of vaccinations to groups with the highest mortality risk is optimal. Our analysis assumes that health systems during winter 2020-2021 have equal staffing and capacity to previous phases of the SARS-CoV-2 epidemic; we do not consider the effects of understaffed hospitals or unvaccinated medical staff. Vaccinating only seronegative individuals avoids redundancy in vaccine use on individuals that may already be immune, and will result in 1% to 2% reductions in cumulative hospitalizations and deaths by mid-2021. Assuming high vaccination coverage (> 28%) and no major relaxations in distancing, masking, gathering size, or hygiene guidelines between now and spring 2021, our model predicts that a combination of vaccination and population immunity will lead to low or near-zero transmission levels by the second quarter of 2021.
RESUMEN
Estimating an infectious disease attack rate requires inference on the number of reported symptomatic cases of a disease, the number of unreported symptomatic cases, and the number of asymptomatic infections. Population-level immunity can then be estimated as the attack rate plus the number of vaccine recipients who had not been previously infected; this requires an estimate of the fraction of vaccines that were distributed to seropositive individuals. To estimate attack rates and population immunity in southern New England, we fit a validated dynamic epidemiological model to case, clinical, and death data streams reported by Rhode Island, Massachusetts, and Connecticut for the first 15 months of the COVID-19 pandemic, from March 1 2020 to May 31 2021. This period includes the initial spring 2020 wave, the major winter wave of 2020-2021, and the lagging wave of lineage B.1.1.7(Alpha) infections during March-April 2021. In autumn 2020, SARS-CoV-2 population immunity (equal to the attack rate at that point) in southern New England was still below 15%, setting the stage for a large winter wave. After the roll-out of vaccines in early 2021, population immunity in many states was expected to approach 70% by spring 2021, with more than half of this immune population coming from vaccinations. Our population immunity estimates for May 31 2021 are 73.4% (95% CrI: 72.9% - 74.1%) for Rhode Island, 64.1% (95% CrI: 64.0% - 64.4%) for Connecticut, and 66.3% (95% CrI: 65.9% - 66.9%) for Massachusetts, indicating that >33% of southern Englanders were still susceptible to infection when the Delta variant began spreading in July 2021. Despite high vaccine coverage in these states, population immunity in summer 2021 was lower than planned due to 34% (Rhode Island), 25% (Connecticut), and 28% (Massachusetts) of vaccine distribution going to seropositive individuals. Future emergency-setting vaccination planning will likely have to consider over-vaccination as a strategy to ensure that high levels of population immunity are reached during the course of an ongoing epidemic.
RESUMEN
In the United States, state-level re-openings in spring 2020 presented an opportunity for the resurgence of SARS-CoV-2 transmission. One important question during this time was whether human contact and mixing patterns could increase gradually without increasing viral transmission, the rationale being that new mixing patterns would likely be associated with improved distancing, masking, and hygiene practices. A second key question to follow during this time was whether clinical characteristics of the epidemic would improve after the initial surge of cases. Here, we analyze age-structured case, hospitalization, and death time series from three states - Rhode Island, Massachusetts, and Pennsylvania - that had successful re-openings in May 2020 without summer waves of infection. Using a Bayesian inference framework on eleven daily data streams and flexible daily population contact parameters, we show that population-average mixing rates dropped by >50% during the lockdown period in March/April, and that the correlation between overall population mobility and transmission-capable mobility was broken in May as these states partially re-opened. We estimate the reporting rates (fraction of symptomatic cases reporting to health system) at 96.0% (RI), 72.1% (MA), and 75.5% (PA); in Rhode Island, when accounting for cases caught through general-population screening programs, the reporting rate estimate is 94.5%. We show that elderly individuals were less able to reduce contacts during the lockdown period when compared to younger individuals. Attack rate estimates through August 31 2020 are 6.4% (95% CI: 5.8% - 7.3%) of the total population infected for Rhode Island, 5.7% (95% CI: 5.0% - 6.8%) in Massachusetts, and 3.7% (95% CI: 3.1% - 4.5%) in Pennsylvania, with some validation available through published seroprevalence studies. Infection fatality rates (IFR) estimates for the spring epidemic are higher in our analysis (>2%) than previously reported values, likely resulting from the epidemics in these three states affecting the most vulnerable sub-populations, especially the most vulnerable of the [≥]80 age group.
RESUMEN
We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern, but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of two months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7s set of mutations.
RESUMEN
Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.
RESUMEN
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic in southern Africa has been characterised by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, whilst the second and third waves were driven by the Beta and Delta variants respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng Province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, predicted to influence antibody neutralization and spike function4. Here, we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.