Accounting for cellular-level variation in lysis: implications for virus-host dynamics.

Viral impacts on microbial populations depend on interaction phenotypes-including viral traits spanning the adsorption rate, latent period, and burst size. The latent period is a key viral trait in lytic infections. Defined as the time from viral adsorption to viral progeny release, the latent period of bacteriophage is conventionally inferred via one-step growth curves in which the accumulation of free virus is measured over time in a population of infected cells. Developed more than 80 years ago, one-step growth curves do not account for cellular-level variability in the timing of lysis, potentially biasing inference of viral traits. Here, we use nonlinear dynamical models to understand how individual-level variation of the latent period impacts virus-host dynamics. Our modeling approach shows that inference of the latent period via one-step growth curves is systematically biased-generating estimates of shorter latent periods than the underlying population-level mean. The bias arises because variability in lysis timing at the cellular level leads to a fraction of early burst events, which are interpreted, artefactually, as an earlier mean time of viral release. We develop a computational framework to estimate latent period variability from joint measurements of host and free virus populations. Our computational framework recovers both the mean and variance of the latent period within simulated infections including realistic measurement noise. This work suggests that reframing the latent period as a distribution to account for variability in the population will improve the study of viral traits and their role in shaping microbial populations.IMPORTANCEQuantifying viral traits-including the adsorption rate, burst size, and latent period-is critical to characterize viral infection dynamics and develop predictive models of viral impacts across scales from cells to ecosystems. Here, we revisit the gold standard of viral trait estimation-the one-step growth curve-to assess the extent to which assumptions at the core of viral infection dynamics lead to ongoing and systematic biases in inferences of viral traits. We show that latent period estimates obtained via one-step growth curves systematically underestimate the mean latent period and, in turn, overestimate the rate of viral killing at population scales. By explicitly incorporating trait variability into a dynamical inference framework that leverages both virus and host time series, we provide a practical route to improve estimates of the mean and variance of viral traits across diverse virus-microbe systems.

Transmission Bottleneck Size Estimation from De Novo Viral Genetic Variation.

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, these approaches have the potential to substantially underestimate true transmission bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arise de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these 2 respiratory viruses.

Transmission bottleneck size estimation from de novo viral genetic variation.

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, allele frequencies can change dramatically over the course of an individual's infection, such that sites that are polymorphic in the donor at the time of transmission may not be polymorphic in the donor at the time of sampling and allele frequencies at donor-polymorphic sites may change dramatically over the course of a recipient's infection. Because of this, transmission bottleneck sizes estimated using allele frequencies observed at a donor's polymorphic sites may be considerable underestimates of true bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arose de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these two respiratory viruses, using an approach that does not tend to underestimate transmission bottleneck sizes.

How time-scale differences in asymptomatic and symptomatic transmission shape SARS-CoV-2 outbreak dynamics.

Asymptomatic and symptomatic SARS-CoV-2 infections can have different characteristic time scales of transmission. These time-scale differences can shape outbreak dynamics as well as bias population-level estimates of epidemic strength, speed, and controllability. For example, prior work focusing on the initial exponential growth phase of an outbreak found that larger time scales for asymptomatic vs. symptomatic transmission can lead to under-estimates of the basic reproduction number as inferred from epidemic case data. Building upon this work, we use a series of nonlinear epidemic models to explore how differences in asymptomatic and symptomatic transmission time scales can lead to changes in the realized proportion of asymptomatic transmission throughout an epidemic. First, we find that when asymptomatic transmission time scales are longer than symptomatic transmission time scales, then the effective proportion of asymptomatic transmission increases as total incidence decreases. Moreover, these time-scale-driven impacts on epidemic dynamics are enhanced when infection status is correlated between infector and infectee pairs (e.g., due to dose-dependent impacts on symptoms). Next we apply these findings to understand the impact of time-scale differences on populations with age-dependent assortative mixing and in which the probability of having a symptomatic infection increases with age. We show that if asymptomatic generation intervals are longer than corresponding symptomatic generation intervals, then correlations between age and symptoms lead to a decrease in the age of infection during periods of epidemic decline (whether due to susceptible depletion or intervention). Altogether, these results demonstrate the need to explore the role of time-scale differences in transmission dynamics alongside behavioral changes to explain outbreak features both at early stages (e.g., in estimating the basic reproduction number) and throughout an epidemic (e.g., in connecting shifts in the age of infection to periods of changing incidence).

Cellular co-infection can modulate the efficiency of influenza A virus production and shape the interferon response.

During viral infection, the numbers of virions infecting individual cells can vary significantly over time and space. The functional consequences of this variation in cellular multiplicity of infection (MOI) remain poorly understood. Here, we rigorously quantify the phenotypic consequences of cellular MOI during influenza A virus (IAV) infection over a single round of replication in terms of cell death rates, viral output kinetics, interferon and antiviral effector gene transcription, and superinfection potential. By statistically fitting mathematical models to our data, we precisely define specific functional forms that quantitatively describe the modulation of these phenotypes by MOI at the single cell level. To determine the generality of these functional forms, we compare two distinct cell lines (MDCK cells and A549 cells), both infected with the H1N1 strain A/Puerto Rico/8/1934 (PR8). We find that a model assuming that infected cell death rates are independent of cellular MOI best fits the experimental data in both cell lines. We further observe that a model in which the rate and efficiency of virus production increase with cellular co-infection best fits our observations in MDCK cells, but not in A549 cells. In A549 cells, we also find that induction of type III interferon, but not type I interferon, is highly dependent on cellular MOI, especially at early timepoints. This finding identifies a role for cellular co-infection in shaping the innate immune response to IAV infection. Finally, we show that higher cellular MOI is associated with more potent superinfection exclusion, thus limiting the total number of virions capable of infecting a cell. Overall, this study suggests that the extent of cellular co-infection by influenza viruses may be a critical determinant of both viral production kinetics and cellular infection outcomes in a host cell type-dependent manner.