Impacts of Vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 Variants Alpha and Delta on Coronavirus Disease 2019 Transmission Dynamics in Four Metropolitan Areas of the United States.

To characterize Coronavirus Disease 2019 (COVID-19) transmission dynamics in each of the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix in 2020 and 2021, we extended a previously reported compartmental model accounting for effects of multiple distinct periods of non-pharmaceutical interventions by adding consideration of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2). For each MSA, we found region-specific parameterizations of the model using daily reports of new COVID-19 cases available from January 21, 2020 to October 31, 2021. In the process, we obtained estimates of the relative infectiousness of Alpha and Delta as well as their takeoff times in each MSA (the times at which sustained transmission began). The estimated infectiousness of Alpha ranged from 1.1x to 1.4x that of viral strains circulating in 2020 and early 2021. The estimated relative infectiousness of Delta was higher in all cases, ranging from 1.6x to 2.1x. The estimated Alpha takeoff times ranged from February 1 to February 28, 2021. The estimated Delta takeoff times ranged from June 2 to June 26, 2021. Estimated takeoff times are consistent with genomic surveillance data.

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit.

SUMMARY: Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. AVAILABILITY AND IMPLEMENTATION: PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Optimal management of stochastic invasion in a metapopulation with Allee effects.

Invasive species account for incalculable damages worldwide, in both ecological and bioeconomic terms. The question of how a network of invasive populations can be optimally managed is one that deserves further exploration. A study accounting for partial observability and imperfect detection, in particular, could yield useful insights into species eradication efforts. Here, we generalized a simple model system that we developed in previous work. This model consists of three interacting populations with underlying strong Allee effects and stochastic dynamics, inhabiting distinct locations connected by dispersal, which can generate bistability. To explore the stochastic dynamics, we formulated an individual-based modeling approach. Next, using the theory of continuous-time Markov chains, we approximated the original high-dimensional model by a Markov chain with eight states, with each state corresponding to a combination of population thresholds. We then used the reduced model as the core for a powerful decision-making tool, referred to as a Partially Observable Markov Decision Process (POMDP). Analysis of this POMDP indicates when the system results in optimal management outcomes.

Robustness and the evolution of length control strategies in the T3SS and flagellar hook.

Bacterial cells construct many structures, such as the flagellar hook and the type III secretion system (T3SS) injectisome, that aid in crucial physiological processes such as locomotion and pathogenesis. Both of these structures involve long extracellular channels, and the length of these channels must be highly regulated in order for these structures to perform their intended functions. There are two leading models for how length control is achieved in the flagellar hook and T3SS needle: the substrate switching model, in which the length is controlled by assembly of an inner rod, and the ruler model, in which a molecular ruler controls the length. Although there is qualitative experimental evidence to support both models, comparatively little has been done to quantitatively characterize these mechanisms or make detailed predictions that could be used to unambiguously test these mechanisms experimentally. In this work, we constructed a mathematical model of length control based on the ruler mechanism and found that the predictions of this model are consistent with experimental data-not just for the scaling of the average length with the ruler protein length, but also for the variance. Interestingly, we found that the ruler mechanism allows for the evolution of needles with large average lengths without the concomitant large increase in variance that occurs in the substrate switching mechanism. In addition to making further predictions that can be tested experimentally, these findings shed new light on the trade-offs that may have led to the evolution of different length control mechanisms in different bacterial species.

Daily Forecasting of Regional Epidemics of Coronavirus Disease with Bayesian Uncertainty Quantification, United States.

To increase situational awareness and support evidence-based policymaking, we formulated a mathematical model for coronavirus disease transmission within a regional population. This compartmental model accounts for quarantine, self-isolation, social distancing, a nonexponentially distributed incubation period, asymptomatic persons, and mild and severe forms of symptomatic disease. We used Bayesian inference to calibrate region-specific models for consistency with daily reports of confirmed cases in the 15 most populous metropolitan statistical areas in the United States. We also quantified uncertainty in parameter estimates and forecasts. This online learning approach enables early identification of new trends despite considerable variability in case reporting.

Crosstalk and ultrasensitivity in protein degradation pathways.

Protein turnover is vital to cellular homeostasis. Many proteins are degraded efficiently only after they have been post-translationally "tagged" with a polyubiquitin chain. Ubiquitylation is a form of Post-Translational Modification (PTM): addition of a ubiquitin to the chain is catalyzed by E3 ligases, and removal of ubiquitin is catalyzed by a De-UBiquitylating enzyme (DUB). Nearly four decades ago, Goldbeter and Koshland discovered that reversible PTM cycles function like on-off switches when the substrates are at saturating concentrations. Although this finding has had profound implications for the understanding of switch-like behavior in biochemical networks, the general behavior of PTM cycles subject to synthesis and degradation has not been studied. Using a mathematical modeling approach, we found that simply introducing protein turnover to a standard modification cycle has profound effects, including significantly reducing the switch-like nature of the response. Our findings suggest that many classic results on PTM cycles may not hold in vivo where protein turnover is ubiquitous. We also found that proteins sharing an E3 ligase can have closely related changes in their expression levels. These results imply that it may be difficult to interpret experimental results obtained from either overexpressing or knocking down protein levels, since changes in protein expression can be coupled via E3 ligase crosstalk. Understanding crosstalk and competition for E3 ligases will be key in ultimately developing a global picture of protein homeostasis.

Tipping Cascades in a Multi-patch System with Noise and Spatial Coupling.

Forecasting tipping points in spatially extended systems is a key area of interest to ecologists. A slowly declining spatially distributed population is an important example of an ecological system that could exhibit a cascade of tipping points. Here, we develop a spatial two-patch model with environmental stochasticity that is slowly forced through population collapse, in the presence of changing environmental conditions. We begin with a basic spatial model, then introduce a fast-slow version of the model using geometric singular perturbation theory, followed by the inclusion of stochasticity. Using the spectral density of the fluctuating subpopulation in each patch, we derive analytic expressions for candidate indicators of population extinction and evaluate their performance through a simulation study. We find that coupling and spatial heterogeneity decrease the magnitude of the proposed indicators in coupled populations relative to isolated populations. Moreover, the degree of coupling dictates the trends in summary statistics. We conclude that this theory may be applied to other contexts, including the control of invasive species.

The Role of Stochasticity in Noise-Induced Tipping Point Cascades: A Master Equation Approach.

Tipping points have been shown to be ubiquitous, both in models and empirically in a range of physical and biological systems. The question of how tipping points cascade through systems has been less explored and is an important one. A study of noise-induced tipping, in particular, could provide key insights into tipping cascades. Here, we consider a specific example of a simple model system that could have cascading tipping points. This model consists of two interacting populations with underlying Allee effects and stochastic dynamics, in separate patches connected by dispersal, which can generate bistability. From an ecological standpoint, we look for rescue effects whereby one population can prevent the collapse of a second population. As a way to investigate the stochastic dynamics, we use an individual-based modeling approach rooted in chemical reaction network theory. Then, using continuous-time Markov chains and the theory of first passage times, we essentially approximate, or emulate, the original high-dimensional model by a Markov chain with just three states, where each state corresponds to a combination of population thresholds. Analysis of this reduced model shows when the system is likely to recover, as well as when tipping cascades through the whole system.

Impacts of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 variants Alpha and Delta on Coronavirus Disease 2019 transmission dynamics in four metropolitan areas of the United States.

To characterize Coronavirus Disease 2019 (COVID-19) transmission dynamics in each of the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix in 2020 and 2021, we extended a previously reported compartmental model accounting for effects of multiple distinct periods of non-pharmaceutical interventions by adding consideration of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2). For each MSA, we found region-specific parameterizations of the model using daily reports of new COVID-19 cases available from January 21, 2020 to October 31, 2021. In the process, we obtained estimates of the relative infectiousness of Alpha and Delta as well as their takeoff times in each MSA (the times at which sustained transmission began). The estimated infectiousness of Alpha ranged from 1.1x to 1.4x that of viral strains circulating in 2020 and early 2021. The estimated relative infectiousness of Delta was higher in all cases, ranging from 1.6x to 2.1x. The estimated Alpha takeoff times ranged from February 1 to February 28, 2021. The estimated Delta takeoff times ranged from June 2 to June 26, 2021. Estimated takeoff times are consistent with genomic surveillance data. One-Sentence Summary: Using a compartmental model parameterized to reproduce available reports of new Coronavirus Disease 2019 (COVID-19) cases, we quantified the impacts of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2) on regional epidemics in the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix.

Differential contagiousness of respiratory disease across the United States.

The initial contagiousness of a communicable disease within a given population is quantified by the basic reproduction number, R0. This number depends on both pathogen and population properties. On the basis of compartmental models that reproduce Coronavirus Disease 2019 (COVID-19) surveillance data, we used Bayesian inference and the next-generation matrix approach to estimate region-specific R0 values for 280 of 384 metropolitan statistical areas (MSAs) in the United States (US), which account for 95% of the US population living in urban areas and 82% of the total population. We focused on MSA populations after finding that these populations were more uniformly impacted by COVID-19 than state populations. Our maximum a posteriori (MAP) estimates for R0 range from 1.9 to 7.7 and quantify the relative susceptibilities of regional populations to spread of respiratory diseases. ONE-SENTENCE SUMMARY: Initial contagiousness of Coronavirus Disease 2019 varied over a 4-fold range across urban areas of the United States.

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah).

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. To investigate the causes of this difference, we used a compartmental model accounting for distinct periods of non-pharmaceutical interventions (NPIs ) ( e.g., behaviors that limit disease transmission) to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah).

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States.

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1-1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States.

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ 0 , the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ 0 relates to a herd immunity threshold (HIT), which is given by 1 - 1/ℛ 0 . When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ 0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ 0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21-January-2020 to 21-June-2020. Our ℛ 0 estimates characterize infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we find that no state has achieved herd immunity as of 20-September-2021. SIGNIFICANCE STATEMENT: COVID-19 will continue to threaten non-immune persons in the presence of ongoing disease transmission. We can estimate when sustained disease transmission will end by calculating the population-specific basic reproduction number ℛ 0 , which relates to a herd immunity threshold (HIT), given by 1 - 1/ℛ 0 . When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely. Here, we report state-level ℛ 0 estimates indicating that disease transmission varies considerably across states. Our ℛ 0 estimates can also be used to determine HITs for the Delta variant of COVID-19. On the basis of Delta-adjusted HITs, vaccination data, and serological survey results, we find that no state has yet achieved herd immunity.

Daily Forecasting of New Cases for Regional Epidemics of Coronavirus Disease 2019 with Bayesian Uncertainty Quantification.

To increase situational awareness and support evidence-based policy-making, we formulated a mathematical model for COVID-19 transmission within a regional population. This compartmental model accounts for quarantine, self-isolation, social distancing, a non-exponentially distributed incubation period, asymptomatic individuals, and mild and severe forms of symptomatic disease. Using Bayesian inference, we have been calibrating region-specific models daily for consistency with new reports of confirmed cases from the 15 most populous metropolitan statistical areas in the United States and quantifying uncertainty in parameter estimates and predictions of future case reports. This online learning approach allows for early identification of new trends despite considerable variability in case reporting. ARTICLE SUMMARY LINE: We report models for regional COVID-19 epidemics and use of Bayesian inference to quantify uncertainty in daily predictions of expected reporting of new cases, enabling identification of new trends in surveillance data.

Daily Forecasting of New Cases for Regional Epidemics of Coronavirus Disease 2019 with Bayesian Uncertainty Quantification.

To increase situational awareness and support evidence-based policy-making, we formulated two types of mathematical models for COVID-19 transmission within a regional population. One is a fitting function that can be calibrated to reproduce an epidemic curve with two timescales (e.g., fast growth and slow decay). The other is a compartmental model that accounts for quarantine, self-isolation, social distancing, a non-exponentially distributed incubation period, asymptomatic individuals, and mild and severe forms of symptomatic disease. Using Bayesian inference, we have been calibrating our models daily for consistency with new reports of confirmed cases from the 15 most populous metropolitan statistical areas in the United States and quantifying uncertainty in parameter estimates and predictions of future case reports. This online learning approach allows for early identification of new trends despite considerable variability in case reporting. We infer new significant upward trends for five of the metropolitan areas starting between 19-April-2020 and 12-June-2020.