Modeling identifies variability in SARS-CoV-2 uptake and eclipse phase by infected cells as principal drivers of extreme variability in nasal viral load in the 48 h post infection

SARS-CoV-2 variant of concern fitness and adaptation in primary human airway epithelia

The severe acute respiratory syndrome coronavirus 2 pandemic is characterized by the emergence of novel variants of concern (VOCs) that replace ancestral strains. Here, we dissect the complex selective pressures by evaluating variant fitness and adaptation in human respiratory tissues. We evaluate viral properties and host responses to reconstruct forces behind D614G through Omicron (BA.1) emergence. We observe differential replication in airway epithelia, differences in cellular tropism, and virus-induced cytotoxicity. D614G accumulates the most mutations after infection, supporting zoonosis and adaptation to the human airway. We perform head-to-head competitions and observe the highest fitness for Gamma and Delta. Under these conditions, RNA recombination favors variants encoding the B.1.617.1 lineage 3' end. Based on viral growth kinetics, Alpha, Gamma, and Delta exhibit increased fitness compared to D614G. In contrast, the global success of Omicron likely derives from increased transmission and antigenic variation. Our data provide molecular evidence to support epidemiological observations of VOC emergence.

Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Computational Modeling Insights into Extreme Heterogeneity in COVID-19 Nasal Swab Data

Throughout the COVID-19 pandemic, an unprecedented level of clinical nasal swab data from around the globe has been collected and shared. Positive tests have consistently revealed viral titers spanning six orders of magnitude! An open question is whether such extreme population heterogeneity is unique to SARS-CoV-2 or possibly generic to viral respiratory infections. To probe this question, we turn to the computational modeling of nasal tract infections. Employing a physiologically faithful, spatially resolved, stochastic model of respiratory tract infection, we explore the statistical distribution of human nasal infections in the immediate 48 h of infection. The spread, or heterogeneity, of the distribution derives from variations in factors within the model that are unique to the infected host, infectious variant, and timing of the test. Hypothetical factors include: (1) reported physiological differences between infected individuals (nasal mucus thickness and clearance velocity); (2) differences in the kinetics of infection, replication, and shedding of viral RNA copies arising from the unique interactions between the host and viral variant; and (3) differences in the time between initial cell infection and the clinical test. Since positive clinical tests are often pre-symptomatic and independent of prior infection or vaccination status, in the model we assume immune evasion throughout the immediate 48 h of infection. Model simulations generate the mean statistical outcomes of total shed viral load and infected cells throughout 48 h for each "virtual individual", which we define as each fixed set of model parameters (1) and (2) above. The "virtual population" and the statistical distribution of outcomes over the population are defined by collecting clinically and experimentally guided ranges for the full set of model parameters (1) and (2). This establishes a model-generated "virtual population database" of nasal viral titers throughout the initial 48 h of infection of every individual, which we then compare with clinical swab test data. Support for model efficacy comes from the sampling of infection dynamics over the virtual population database, which reproduces the six-order-of-magnitude clinical population heterogeneity. However, the goal of this study is to answer a deeper biological and clinical question. What is the impact on the dynamics of early nasal infection due to each individual physiological feature or virus-cell kinetic mechanism? To answer this question, global data analysis methods are applied to the virtual population database that sample across the entire database and de-correlate (i.e., isolate) the dynamic infection outcome sensitivities of each model parameter. These methods predict the dominant, indeed exponential, driver of population heterogeneity in dynamic infection outcomes is the latency time of infected cells (from the moment of infection until onset of viral RNA shedding). The shedding rate of the viral RNA of infected cells in the shedding phase is a strong, but not exponential, driver of infection. Furthermore, the unknown timing of the nasal swab test relative to the onset of infection is an equally dominant contributor to extreme population heterogeneity in clinical test data since infectious viral loads grow from undetectable levels to more than six orders of magnitude within 48 h.

The power of weak, transient interactions across biology: A paradigm of emergent behavior

A growing list of diverse biological systems and their equally diverse functionalities provides realizations of a paradigm of emergent behavior. In each of these biological systems, pervasive ensembles of weak, short-lived, spatially local interactions act autonomously to convey functionalities at larger spatial and temporal scales. In this article, a range of diverse systems and functionalities are presented in a cursory manner with literature citations for further details. Then two systems and their properties are discussed in more detail: yeast chromosome biology and human respiratory mucus.