Global estimates of the fitness advantage of SARS-CoV-2 variant Omicron

New variants of SARS-CoV-2 show remarkable heterogeneity in their relative fitness both over time and space. In this paper we extend a previously published model for estimating the selection strength for new SARS-CoV-2 variants to a hierarchical, mixed-effects, renewal equation model. This formulation allows us to globally estimate selection effects at different spatial levels while controlling for complex patterns of transmission and jointly inferring the effects of unit-level covariates in the spatial heterogeneity of SARS-CoV-2 selection effects. Applying this model to the spread of Omicron in 40 counties finding evidence for very strong (64%) but very heterogeneous selection effects at the country level. We further considered different measures of vaccination levels and measures of recent population-level infection as possible explanations. However, none of those variables were found to explain a significant proportion of the heterogeneity in country-level selection effects. We did find a significant positive correlation between the selection advantage of Delta and Omicron at the country level, suggesting that region-specific explanatory variables of fitness differences do exist. Our method is implemented in the Stan programming language, can be run on standard commercial-grade computing resources, and should be straightforward to apply to future variants.

A simple model of COVID-19 explains disease severity and the effect of treatments

Considerable effort was made to better understand why some people suffer from severe COVID-19 while others remain asymptomatic. This has led to important clinical findings; people with severe COVID-19 generally experience persistently high levels of inflammation, slower viral load decay, display a dysregulated type-I interferon response, have less active natural killer cells and increased levels of neutrophil extracellular traps. How these findings are connected to the pathogenesis of COVID-19 remains unclear. We propose a mathematical model that sheds light on this issue. The model focuses on cells that trigger inflammation through molecular patterns: infected cells carrying pathogen-associated molecular patterns (PAMPs) and damaged cells producing damage-associated molecular patterns (DAMPs). The former signals the presence of pathogens while the latter signals danger such as hypoxia or the lack of nutrients. Analyses show that SARS-CoV-2 infections can lead to a self-perpetuating feedback loop between DAMP expressing cells and inflammation. It identifies the inability to quickly clear PAMPs and DAMPs as the main contributor to hyperinflammation. The model explains clinical findings and the conditional impact of treatments on disease severity. The simplicity of the model and its high level of consistency with clinical findings motivate its use for the formulation of new treatment strategies.

Emergence of SARS-CoV-2 Resistance with Monoclonal Antibody Therapy

Resistance mutations to monoclonal antibody (mAb) therapy has been reported, but in the non-immunosuppressed population, it is unclear if in vivo emergence of SARS-CoV-2 resistance mutations alters either viral replication dynamics or therapeutic efficacy. In ACTIV-2/A5401, non-hospitalized participants with symptomatic SARS-CoV-2 infection were randomized to bamlanivimab (700mg or 7000mg) or placebo. Treatment-emergent resistance mutations were significantly more likely detected after bamlanivimab 700mg treatment than placebo (7% of 111 vs 0% of 112 participants, P=0.003). There were no treatment-emergent resistance mutations among the 48 participants who received bamlanivimab 7000mg. Participants with emerging mAb resistant virus had significantly higher pre-treatment nasopharyngeal and anterior nasal viral load. Intensive respiratory tract viral sampling revealed the dynamic nature of SARS-CoV-2 evolution, with evidence of rapid and sustained viral rebound after emergence of resistance mutations, and worsened symptom severity. Participants with emerging bamlanivimab resistance often accumulated additional polymorphisms found in current variants of concern/interest and associated with immune escape. These results highlight the potential for rapid emergence of resistance during mAb monotherapy treatment, resulting in prolonged high level respiratory tract viral loads and clinical worsening. Careful virologic assessment should be prioritized during the development and clinical implementation of antiviral treatments for COVID-19.

Longitudinal analysis of SARS-CoV-2 vaccine breakthrough infections reveal limited infectious virus shedding and restricted tissue distribution

The global effort to vaccinate people against SARS-CoV-2 in the midst of an ongoing pandemic has raised questions about the nature of vaccine breakthrough infections and the potential for vaccinated individuals to transmit the virus. These questions have become even more urgent as new variants of concern with enhanced transmissibility, such as Delta, continue to emerge. To shed light on how vaccine breakthrough infections compare with infections in immunologically naive individuals, we examined viral dynamics and infectious virus shedding through daily longitudinal sampling in a small cohort of adults infected with SARS-CoV-2 at varying stages of vaccination. The durations of both infectious virus shedding and symptoms were significantly reduced in vaccinated individuals compared with unvaccinated individuals. We also observed that breakthrough infections are associated with strong tissue compartmentalization and are only detectable in saliva in some cases. These data indicate that vaccination shortens the duration of time of high transmission potential, minimizes symptom duration, and may restrict tissue dissemination.

Daily sampling of early SARS-CoV-2 infection reveals substantial heterogeneity in infectiousness

The dynamics of SARS-CoV-2 replication and shedding in humans remain poorly understood. We captured the dynamics of infectious virus and viral RNA shedding during acute infection through daily longitudinal sampling of 60 individuals for up to 14 days. By fitting mechanistic models, we directly estimate viral reproduction and clearance rates, and overall infectiousness for each individual. Significant person-to-person variation in infectious virus shedding suggests that individual-level heterogeneity in viral dynamics contributes to superspreading. Viral genome load often peaked days earlier in saliva than in nasal swabs, indicating strong compartmentalization and suggesting that saliva may serve as a superior sampling site for early detection of infection. Viral loads and clearance kinetics of B.1.1.7 and non-B.1.1.7 viruses in nasal swabs were indistinguishable, however B.1.1.7 exhibited a significantly slower pre-peak growth rate in saliva. These results provide a high-resolution portrait of SARS-CoV-2 infection dynamics and implicate individual-level heterogeneity in infectiousness in superspreading.

In vivo kinetics of SARS-CoV-2 infection and its relationship with a person's infectiousness

The within-host viral kinetics of SARS-CoV-2 infection and how they relate to a persons infectiousness are not well understood. This limits our ability to quantify the impact of interventions on viral transmission. Here, we develop data-driven viral dynamic models of SARS-CoV-2 infection and estimate key within-host parameters such as the infected cell half-life and the within-host reproductive number. We then develop a model linking VL to infectiousness, showing that a persons infectiousness increases sub-linearly with VL. We show that the logarithm of the VL in the upper respiratory tract (URT) is a better surrogate of infectiousness than the VL itself. Using data on VL and the predicted infectiousness, we further incorporated data on antigen and reverse transcription polymerase chain reaction (RT-PCR) tests and compared their usefulness in detecting infection and preventing transmission. We found that RT-PCR tests perform better than antigen tests assuming equal testing frequency; however, more frequent antigen testing may perform equally well with RT-PCR tests at a lower cost, but with many more false-negative tests. Overall, our models provide a quantitative framework for inferring the impact of therapeutics and vaccines that lower VL on the infectiousness of individuals and for evaluating rapid testing strategies. SignificanceQuantifying the kinetics of SARS-CoV-2 infection and individual infectiousness is key to quantitatively understanding SARS-CoV-2 transmission and evaluating intervention strategies. Here we developed data-driven within-host models of SARS-CoV-2 infection and by fitting them to clinical data we estimated key within-host viral dynamic parameters. We also developed a mechanistic model for viral transmission and show that the logarithm of the viral load in the upper respiratory tract serves an appropriate surrogate for a persons infectiousness. Using data on how viral load changes during infection, we further evaluated the effectiveness of PCR and antigen-based testing strategies for averting transmission and identifying infected individuals.

Assessing the Impact of Human Mobility to Predict Regional Excess Death in Ecuador

COVID-19 outbreaks have had high mortality in low- and medium-income countries such as Ecuador. Human mobility is an important factor influencing the spread of diseases possibly leading to a high burden of disease at the country level. Drastic control measures, such as complete lockdown are effective epidemic controls, yet in practice, one hopes that a partial shutdown would suffice. It is an open problem to determine how much mobility can be allowed while controlling an outbreak. In this paper, we use statistical models to relate human mobility to the excess death in Ecuador while controlling for demographic factors. The mobility index provided by GRANDATA, based on mobile phone users, represents the change of number of out-of-home events with respect to a benchmark date (March 2nd, the first date the data is available). The study confirms the global trend that more men are dying than expected compared to women, and that people under 30 show less deaths than expected. Specifically, individuals in the age groups younger than 20, we found have their death rate reduced during the pandemic between 22% and 27% of the expected deaths in the absence of COVID-19. The weekly median mobility time series shows a sharp decrease in human mobility immediately after a national lockdown was declared on March 17, 2020 and a progressive increase towards the pre-lockdown level within two months. Relating median mobility to excess death shows a lag in its effect: first, a decrease in mobility in the previous two to three weeks decreases the excess death and more novel, we found that an increase of mobility variability four weeks prior, increases the number of excess deaths.

Estimating the strength of selection for new SARS-CoV-2 variants

Controlling the SARS-CoV-2 pandemic becomes increasingly challenging as the virus adapts to human hosts through the continual emergence of more transmissible variants. Simply observing that a variant is increasing in frequency is relatively straightforward, but more sophisticated methodology is needed to determine whether a new variant is a global threat and the magnitude of its selective advantage. We present three methods for quantifying the strength of selection for new and emerging variants of SARS-CoV-2 relative to the background of contemporaneous variants. These methods range from a detailed model of dynamics within one country to a broad analysis across all countries, and they include alternative explanations such as migration and drift. We find evidence for strong selection favoring the D614G spike mutation and B.1.1.7 (Alpha), weaker selection favoring B.1.351 (Beta), and no advantage of R.1 after it spreads beyond Japan. Cutting back data to earlier time horizons reveals large uncertainty very soon after emergence, but that estimates of selection stabilize after several weeks. Our results also show substantial heterogeneity among countries, demonstrating the need for a truly global perspective on the molecular epidemiology of SARS-CoV-2.

Longitudinal assessment of diagnostic test performance over the course of acute SARS-CoV-2 infection

What is already known about this topic?Diagnostic tests and sample types for SARS-CoV-2 vary in sensitivity across the infection period. What is added by this report?We show that both RTqPCR (from nasal swab and saliva) and the Quidel SARS Sofia FIA rapid antigen tests peak in sensitivity during the period in which live virus can be detected in nasal swabs, but that the sensitivity of RTqPCR tests rises more rapidly in the pre-infectious period. We also use empirical data to estimate the sensitivities of RTqPCR and antigen tests as a function of testing frequency. What are the implications for public health practice?RTqPCR tests will be more effective than rapid antigen tests at identifying infected individuals prior to or early during the infectious period and thus for minimizing forward transmission (provided results reporting is timely). All modalities, including rapid antigen tests, showed >94% sensitivity to detect infection if used at least twice per week. Regular surveillance/screening using rapid antigen tests 2-3 times per week can be an effective strategy to achieve high sensitivity (>95%) for identifying infected individuals.

Excess deaths reveal the true spatial, temporal, and demographic impact of COVID-19 on mortality in Ecuador

BackgroundIn early 2020, Ecuador reported one of the highest surges of per capita deaths across the globe. MethodsWe collected a comprehensive dataset containing individual death records between 2015 and 2020 from the Ecuadorian National Institute of Statistics and Census and the Ecuadorian Ministry of Government. We computed the number of excess deaths across time, geographical locations and demographic groups using Poisson regression methods. ResultsBetween January 1st and September 23rd, 2020, the number of excess deaths in Ecuador is 36,402 (95% CI: 35,762-36,827) or 208 per 105 population, which is 171% of the expected deaths in that period in a typical year. Only 20% of the excess deaths are attributable to confirmed COVID-19 deaths. Strikingly, in provinces that were most affected by COVID-19, such as Guayas and Santa Elena, the all-cause deaths are more than double the expected number of deaths that would have occurred in a normal year. The extent of excess deaths in men is higher than in women, and the number of excess deaths increases with age. Indigenous populations had the highest level of excess deaths among all ethnic groups. ConclusionsOverall, the exceptionally high level of excess deaths in Ecuador highlights the enormous burden and heterogeneous impact of COVID-19 on mortality especially in older age groups and indigenous populations in Ecuador that was not fully revealed by COVID-19 death counts. Together with the limited testing in Ecuador, our results suggest that the majority of the excess deaths were likely to be undocumented COVID-19 deaths.

Kinetics of SARS-CoV-2 infection in the human upper and lower respiratory tracts and their relationship with infectiousness

SARS-CoV-2 is a human pathogen that causes infection in both the upper respiratory tract (URT) and the lower respiratory tract (LRT). The viral kinetics of SARS-CoV-2 infection and how they relate to infectiousness and disease progression are not well understood. Here, we develop data-driven viral dynamic models of SARS-CoV-2 infection in both the URT and LRT. We fit the models to viral load data from patients with likely infection dates known, we estimated that infected individuals with a longer incubation period had lower rates of viral growth, took longer to reach peak viremia in the URT, and had higher chances of presymptomatic transmission. We then developed a model linking viral load to infectiousness. We found that to explain the substantial fraction of transmissions occurring presymptomatically, the infectiousness of a person should depend on a saturating function of the viral load, making the logarithm of the URT viral load a better surrogate of infectiousness than the viral load itself. Comparing the roles of target-cell limitation, the innate immune response, proliferation of target cells and spatial infection in the LRT, we found that spatial dissemination in the lungs is likely to be an important process in sustaining the prolonged high viral loads. Overall, our models provide a quantitative framework for predicting how SARS-CoV-2 within-host dynamics determine infectiousness and represent a step towards quantifying how viral load dynamics and the immune responses determine disease severity.

Estimating the reproductive number R0 of SARS-CoV-2 in the United States and eight European countries and implications for vaccination

SARS-CoV-2 rapidly spread from a regional outbreak to a global pandemic in just a few months. Global research efforts have focused on developing effective vaccines against SARS-CoV-2 and the disease it causes, COVID-19. However, some of the basic epidemiological parameters, such as the exponential epidemic growth rate and the basic reproductive number, R0, across geographic areas are still not well quantified. Here, we developed and fit a mathematical model to case and death count data collected from the United States and eight European countries during the early epidemic period before broad control measures were implemented. Results show that the early epidemic grew exponentially at rates between 0.19-0.29/day (epidemic doubling times between 2.4-3.6 days). We discuss the current estimates of the mean serial interval, and argue that existing evidence suggests that the interval is between 6-8 days in the absence of active isolation efforts. Using parameters consistent with this range, we estimated the median R0 value to be 5.8 (confidence interval: 4.7-7.3) in the United States and between 3.6 and 6.1 in the eight European countries. This translates to herd immunity thresholds needed to stop transmission to be between 73% and 84%. We further analyze how vaccination schedules depends on R0, the duration of vaccine-induced immunity to SARS-CoV-2, and show that individual-level heterogeneity in vaccine induced immunity can significantly affect vaccination schedules. SignificanceWith the global efforts to develop vaccines for COVID-19, it is important to understand the contagiousness of the virus to design regional vaccination policy. To that end, we fit a mathematical model to data collected from the early epidemic period in the United States and eight European countries, estimating that the early epidemic doubles between 2.4-3.6 days. This suggests that SARS-CoV-2 is highly transmissible in the absence of strong control measures irrespective of heterogeneity in geographic and social settings. We estimated the median basic reproduction number, R0 to be 5.8 (confidence interval: 4.7-7.3) in the United States and between 3.6 and 6.1 in the eight European countries. The herd immunity needed to stop transmission is high (between 73% and 84%).

Decline in global COVID-19 transmission

AO_SCPLOWBSTRACTC_SCPLOWWe analyzed COVID-19 data through May 6th, 2020 using a partially observed Markov process. Our method uses a hybrid deterministic and stochastic formalism that allows for time variable transmission rates and detection probabilities. The model was fit using iterated particle filtering to case count and death count time series from 55 countries. We found evidence for a shrinking epidemic in 30 of the 55 examined countries. Of those 30 countries, 27 have significant evidence for subcritical transmission rates, although the decline in new cases is relatively slow compared to the initial growth rates. Generally, the transmission rates in Europe were lower than in the Americas and Asia. This suggests that global scale social distancing efforts to slow the spread of COVID-19 are effective although they need to be strengthened in many regions and maintained in others to avoid further resurgence of COVID-19. The slow decline also suggests alternative strategies to control the virus are needed before social distancing efforts are partially relaxed.

Timing of antiviral treatment initiation is critical to reduce SARS-Cov-2 viral load

We modeled the viral dynamics of 13 untreated patients infected with SARS-CoV-2 to infer viral growth parameters and predict the effects of antiviral treatments. In order to reduce peak viral load by more than 2 logs, drug efficacy needs to be greater than 90% if treatment is administered after symptom onset; an efficacy of 60% could be sufficient if treatment is initiated before symptom onset. Given their pharmacokinetic/pharmacodynamic properties, current investigated drugs may be in a range of 6-87% efficacy. They may help control virus if administered very early, but may not have a major effect in severe patients.

Fast spread of COVID-19 in Europe and the US and its implications: even modest public health goals require comprehensive intervention

AbstratThe COVID-19 pandemic caused more than 800,000 infections and 40,000 deaths by the end of March 2020. However, some of the basic epidemiological parameters, such as the exponential epidemic growth rate and R0 are debated. We developed an inference approach to control for confounding factors in data collection, such as underreporting and changes in surveillance intensities, and fitted a mathematical model to infection and death count data collected from eight European countries and the US. In all countries, the early epidemic grew exponentially at rates between 0.19-0.29/day (epidemic doubling times between 2.4-3.7 days). This suggests a highly infectious virus with an R0 likely between 4.0 and 7.1. We show that similar levels of intervention efforts are needed, no matter the goal is mitigation or containment. Early, strong and comprehensive intervention efforts to achieve greater than 74-86% reduction in transmission are necessary. One-sentence SummaryWe estimated that COVID-19 spreads rapidly in 8 European countries and the US, suggesting that early, strong and comprehensive interventions are necessary.

The Novel Coronavirus, 2019-nCoV, is Highly Contagious and More Infectious Than Initially Estimated

The novel coronavirus (2019-nCoV) is a recently emerged human pathogen that has spread widely since January 2020. Initially, the basic reproductive number, R0, was estimated to be 2.2 to 2.7. Here we provide a new estimate of this quantity. We collected extensive individual case reports and estimated key epidemiology parameters, including the incubation period. Integrating these estimates and high-resolution real-time human travel and infection data with mathematical models, we estimated that the number of infected individuals during early epidemic double every 2.4 days, and the R0 value is likely to be between 4.7 and 6.6. We further show that quarantine and contact tracing of symptomatic individuals alone may not be effective and early, strong control measures are needed to stop transmission of the virus. One-sentence summaryBy collecting and analyzing spatiotemporal data, we estimated the transmission potential for 2019-nCoV.