Modeling stochastic processes in disease spread across a heterogeneous social system.
Ontology highlight
ABSTRACT:
SUBMITTER: Kim M
PROVIDER: S-EPMC6329989 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We present results of a study of a simple, stochastic, agent-based model of influenza A infection, simulating its dynamics over the course of one flu season. Building on an early work of Bartlett, we define a model with a limited number of parameters and rates that have clear epidemiological interpretation and can be constrained by data. We demonstrate the occurrence of recurrent behavior in the infected number [more than one peak in a season], which is observed in data, in our simulations for populations consisting of cohorts with strong intra- and weak inter-cohort transmissibility. We examine the dependence of the results on epidemiological and population characteristics by investigating their dependence on a range of parameter values. Finally, we study infection with two strains of influenza, inspired by observations, and show a counter-intuitive result for the effect of inoculation against the strain that leads to the first wave of infection.
| S-EPMC7173783 | biostudies-literature
Project description:The use of RNA-seq as the preferred method for the discovery and validation of small RNA biomarkers is hindered by high variability and biased sequence counts. In this paper we develop a statistical model for sequence counts that accounts for ligase bias and stochastic variation in library amplification steps and sequencing depth variation. Our analytical contributions are the description of the Linear Quadratic (LQ) relation between the mean and variance of the sequence counts in an RNA-seq experiment and the derivation of the Poisson truncated mixture as the underlying probability distribution for RNA-seq data. Using a large number of sequencing datasets, we demonstrate here how one can use this modeling framework to calculate empirical correction factors for ligase bias, while accounting for random variation in sequence counts. Bias correction may remove the majority of bias in the absence of differential expression and more than 40% of the bias in the presence of variable expression of miRNAs. Empirical bias correction factors appear to be nearly constant over at least one and up to four orders of magnitude of total RNA input and independent of sample composition.
2017-03-15 | GSE93399 | GEO
Project description:Age-dependent tissue decline and increased cancer incidence are widely accepted to be rate-limited by the accumulation of somatic mutations over time. Current models of carcinogenesis are dominated by the assumption that oncogenic mutations have defined advantageous fitness effects on recipient stem and progenitor cells, promoting and rate-limiting somatic evolution. However, this assumption is markedly discrepant with evolutionary theory, whereby fitness is a dynamic property of a phenotype imposed upon and widely modulated by environment. We computationally modeled dynamic microenvironment-dependent fitness alterations in hematopoietic stem cells (HSC) within the Sprengel-Liebig system known to govern evolution at the population level. Our model for the first time integrates real data on age-dependent dynamics of HSC division rates, pool size, and accumulation of genetic changes and demonstrates that somatic evolution is not rate-limited by the occurrence of mutations, but instead results from aged microenvironment-driven alterations in the selective/fitness value of previously accumulated genetic changes. Our results are also consistent with evolutionary models of aging and thus oppose both somatic mutation-centric paradigms of carcinogenesis and tissue functional decline. In total, we demonstrate that aging directly promotes HSC fitness decline and somatic evolution via non-cell-autonomous mechanisms.
| S-EPMC4298364 | biostudies-literature
Project description:The use of RNA-seq as the preferred method for the discovery and validation of small RNA biomarkers has been hindered by high quantitative variability and biased sequence counts. In this paper we develop a statistical model for sequence counts that accounts for ligase bias and stochastic variation in sequence counts. This model implies a linear quadratic relation between the mean and variance of sequence counts. Using a large number of sequencing datasets, we demonstrate how one can use the generalized additive models for location, scale and shape (GAMLSS) distributional regression framework to calculate and apply empirical correction factors for ligase bias. Bias correction could remove more than 40% of the bias for miRNAs. Empirical bias correction factors appear to be nearly constant over at least one and up to four orders of magnitude of total RNA input and independent of sample composition. Using synthetic mixes of known composition, we show that the GAMLSS approach can analyze differential expression with greater accuracy, higher sensitivity and specificity than six existing algorithms (DESeq2, edgeR, EBSeq, limma, DSS, voom) for the analysis of small RNA-seq data.
| S-EPMC5499834 | biostudies-literature
Project description:Wolbachia are very common, maternally transmitted endosymbionts of insects. They often spread by a mechanism termed cytoplasmic incompatibility (CI) that involves reduced egg hatch when Wolbachia-free ova are fertilized by sperm from Wolbachia-infected males. Because the progeny of Wolbachia-infected females generally do not suffer CI-induced mortality, infected females are often at a reproductive advantage in polymorphic populations. Deterministic models show that Wolbachia that impose no costs on their hosts and have perfect maternal transmission will spread from arbitrarily low frequencies (though initially very slowly); otherwise, there will be a threshold frequency below which Wolbachia frequencies decline to extinction and above which they increase to fixation or a high stable equilibrium. Stochastic theory was used to calculate the probability of fixation in populations of different size for arbitrary current frequencies of Wolbachia, with special attention paid to the case of spread after the arrival of a single infected female. Exact results are given based on a Moran process that assumes a specific demographic model, and approximate results are obtained using the more general Wright-Fisher theory. A new analytical approximation for the probability of fixation is derived, which performs well for small population sizes. The significance of stochastic effects in the natural spread of Wolbachia and their relevance to the use of Wolbachia as a drive mechanism in vector and pest management are discussed.
| S-EPMC2605827 | biostudies-other
Project description:Natural and human-induced events are continuously altering the structure of our landscapes and as a result impacting the spatial relationships between individual landscape elements and the species living in the area. Yet, only recently has the influence of the surrounding landscape on invasive species spread started to be considered. The scientific community increasingly recognizes the need for broader modeling framework that focuses on cross-study comparisons at different spatiotemporal scales. Using two illustrative examples, we introduce a general modeling framework that allows for a systematic investigation of the effect of habitat change on invasive species establishment and spread. The essential parts of the framework are (i) a mechanistic spatially explicit model (a modular dispersal framework-MDIG) that allows population dynamics and dispersal to be modeled in a geographical information system (GIS), (ii) a landscape generator that allows replicated landscape patterns with partially controllable spatial properties to be generated, and (iii) landscape metrics that depict the essential aspects of landscape with which dispersal and demographic processes interact. The modeling framework provides functionality for a wide variety of applications ranging from predictions of the spatiotemporal spread of real species and comparison of potential management strategies, to theoretical investigation of the effect of habitat change on population dynamics. Such a framework allows to quantify how small-grain landscape characteristics, such as habitat size and habitat connectivity, interact with life-history traits to determine the dynamics of invasive species spread in fragmented landscape. As such, it will give deeper insights into species traits and landscape features that lead to establishment and spread success and may be key to preventing new incursions and the development of efficient monitoring, surveillance, control or eradication programs.
| S-EPMC5648669 | biostudies-literature
Project description:BackgroundThe COVID-19 (the disease caused by the SARS-CoV-2 virus) pandemic has underscored the need for additional data, tools, and methods that can be used to combat emerging and existing public health concerns. Since March 2020, there has been substantial interest in using social media data to both understand and intervene in the pandemic. Researchers from many disciplines have recently found a relationship between COVID-19 and a new data set from Facebook called the Social Connectedness Index (SCI).ObjectiveBuilding off this work, we seek to use the SCI to examine how social similarity of Missouri counties could explain similarities of COVID-19 cases over time. Additionally, we aim to add to the body of literature on the utility of the SCI by using a novel modeling technique.MethodsIn September 2020, we conducted this cross-sectional study using publicly available data to test the association between the SCI and COVID-19 spread in Missouri using exponential random graph models, which model relational data, and the outcome variable must be binary, representing the presence or absence of a relationship. In our model, this was the presence or absence of a highly correlated COVID-19 case count trajectory between two given counties in Missouri. Covariates included each county's total population, percent rurality, and distance between each county pair.ResultsWe found that all covariates were significantly associated with two counties having highly correlated COVID-19 case count trajectories. As the log of a county's total population increased, the odds of two counties having highly correlated COVID-19 case count trajectories increased by 66% (odds ratio [OR] 1.66, 95% CI 1.43-1.92). As the percent of a county classified as rural increased, the odds of two counties having highly correlated COVID-19 case count trajectories increased by 1% (OR 1.01, 95% CI 1.00-1.01). As the distance (in miles) between two counties increased, the odds of two counties having highly correlated COVID-19 case count trajectories decreased by 43% (OR 0.57, 95% CI 0.43-0.77). Lastly, as the log of the SCI between two Missouri counties increased, the odds of those two counties having highly correlated COVID-19 case count trajectories significantly increased by 17% (OR 1.17, 95% CI 1.09-1.26).ConclusionsThese results could suggest that two counties with a greater likelihood of sharing Facebook friendships means residents of those counties have a higher likelihood of sharing similar belief systems, in particular as they relate to COVID-19 and public health practices. Another possibility is that the SCI is picking up travel or movement data among county residents. This suggests the SCI is capturing a unique phenomenon relevant to COVID-19 and that it may be worth adding to other COVID-19 models. Additional research is needed to better understand what the SCI is capturing practically and what it means for public health policies and prevention practices.
| S-EPMC8675563 | biostudies-literature
Project description:Delay is an important and ubiquitous aspect of many biochemical processes. For example, delay plays a central role in the dynamics of genetic regulatory networks as it stems from the sequential assembly of first mRNA and then protein. Genetic regulatory networks are therefore frequently modeled as stochastic birth-death processes with delay. Here, we examine the relationship between delay birth-death processes and their appropriate approximating delay chemical Langevin equations. We prove a quantitative bound on the error between the pathwise realizations of these two processes. Our results hold for both fixed delay and distributed delay. Simulations demonstrate that the delay chemical Langevin approximation is accurate even at moderate system sizes. It captures dynamical features such as the oscillatory behavior in negative feedback circuits, cross-correlations between nodes in a network, and spatial and temporal information in two commonly studied motifs of metastability in biochemical systems. Overall, these results provide a foundation for using delay stochastic differential equations to approximate the dynamics of birth-death processes with delay.
| S-EPMC4039735 | biostudies-other
Project description:We introduce a multiscale dataset obtained from Heihe Watershed Allied Telemetry Experimental Research (HiWATER) in an oasis-desert area in 2012. Upscaling of eco-hydrological processes on a heterogeneous surface is a grand challenge. Progress in this field is hindered by the poor availability of multiscale observations. HiWATER is an experiment designed to address this challenge through instrumentation on hierarchically nested scales to obtain multiscale and multidisciplinary data. The HiWATER observation system consists of a flux observation matrix of eddy covariance towers, large aperture scintillometers, and automatic meteorological stations; an eco-hydrological sensor network of soil moisture and leaf area index; hyper-resolution airborne remote sensing using LiDAR, imaging spectrometer, multi-angle thermal imager, and L-band microwave radiometer; and synchronical ground measurements of vegetation dynamics, and photosynthesis processes. All observational data were carefully quality controlled throughout sensor calibration, data collection, data processing, and datasets generation. The data are freely available at figshare and the Cold and Arid Regions Science Data Centre. The data should be useful for elucidating multiscale eco-hydrological processes and developing upscaling methods.
| S-EPMC5501970 | biostudies-literature
Project description:In the framework of coding theory, under the assumption of a Markov process (Xt) on a finite alphabet A, the compressed representation of the data will be composed of a description of the model used to code the data and the encoded data. Given the model, the Huffman's algorithm is optimal for the number of bits needed to encode the data. On the other hand, modeling (Xt) through a Partition Markov Model (PMM) promotes a reduction in the number of transition probabilities needed to define the model. This paper shows how the use of Huffman code with a PMM reduces the number of bits needed in this process. We prove the estimation of a PMM allows for estimating the entropy of (Xt), providing an estimator of the minimum expected codeword length per symbol. We show the efficiency of the new methodology on a simulation study and, through a real problem of compression of DNA sequences of SARS-CoV-2, obtaining in the real data at least a reduction of 10.4%.
| S-EPMC8775014 | biostudies-literature