Inferring rates of metastatic dissemination using stochastic network models.
Ontology highlight
ABSTRACT:
SUBMITTER: Gerlee P
PROVIDER: S-EPMC6459558 | biostudies-literature | 2019 Apr
REPOSITORIES: biostudies-literature
Publications
Inferring rates of metastatic dissemination using stochastic network models.
Gerlee Philip P Johansson Mia M
PLoS computational biology 20190401 4
The formation of metastases is driven by the ability of cancer cells to disseminate from the site of the primary tumour to target organs. The process of dissemination is constrained by anatomical features such as the flow of blood and lymph in the circulatory system. We exploit this fact in a stochastic network model of metastasis formation, in which only anatomically feasible routes of dissemination are considered. By fitting this model to two different clinical datasets (tongue & ovarian cance ...[more]
Similar Datasets
Project description:Estimation of epidemiological and population parameters from molecular sequence data has become central to the understanding of infectious disease dynamics. Various models have been proposed to infer details of the dynamics that describe epidemic progression. These include inference approaches derived from Kingman's coalescent theory. Here, we use recently described coalescent theory for epidemic dynamics to develop stochastic and deterministic coalescent susceptible-infected-removed (SIR) tree priors. We implement these in a Bayesian phylogenetic inference framework to permit joint estimation of SIR epidemic parameters and the sample genealogy. We assess the performance of the two coalescent models and also juxtapose results obtained with a recently published birth-death-sampling model for epidemic inference. Comparisons are made by analyzing sets of genealogies simulated under precisely known epidemiological parameters. Additionally, we analyze influenza A (H1N1) sequence data sampled in the Canterbury region of New Zealand and HIV-1 sequence data obtained from known United Kingdom infection clusters. We show that both coalescent SIR models are effective at estimating epidemiological parameters from data with large fundamental reproductive number [Formula: see text] and large population size [Formula: see text]. Furthermore, we find that the stochastic variant generally outperforms its deterministic counterpart in terms of error, bias, and highest posterior density coverage, particularly for smaller [Formula: see text] and [Formula: see text]. However, each of these inference models is shown to have undesirable properties in certain circumstances, especially for epidemic outbreaks with [Formula: see text] close to one or with small effective susceptible populations.
| S-EPMC4317665 | biostudies-literature
Project description:Adaptive introgression-the flow of adaptive genetic variation between species or populations-has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.
| S-EPMC8097282 | biostudies-literature
Project description:Often the time derivative of a measured variable is of as much interest as the variable itself. For a growing population of biological cells, for example, the population's growth rate is typically more important than its size. Here we introduce a non-parametric method to infer first and second time derivatives as a function of time from time-series data. Our approach is based on Gaussian processes and applies to a wide range of data. In tests, the method is at least as accurate as others, but has several advantages: it estimates errors both in the inference and in any summary statistics, such as lag times, and allows interpolation with the corresponding error estimation. As illustrations, we infer growth rates of microbial cells, the rate of assembly of an amyloid fibril and both the speed and acceleration of two separating spindle pole bodies. Our algorithm should thus be broadly applicable.
| S-EPMC5159892 | biostudies-literature
Project description:As population-level patterns of interest in forests emerge from individual vital rates, modelling forest dynamics requires making the link between the scales at which data are collected (individual stems) and the scales at which questions are asked (e.g. populations and communities). Structured population models (e.g. integral projection models (IPMs)) are useful tools for linking vital rates to population dynamics. However, the application of such models to forest trees remains challenging owing to features of tree life cycles, such as slow growth, long lifespan and lack of data on crucial ontogenic stages. We developed a survival model that accounts for size-dependent mortality and a growth model that characterizes individual heterogeneity. We integrated vital rate models into two types of population model; an analytically tractable form of IPM and an individual-based model (IBM) that is applied with stochastic simulations. We calculated longevities, passage times to, and occupancy time in, different life cycle stages, important metrics for understanding how demographic rates translate into patterns of forest turnover and carbon residence times. Here, we illustrate the methods for three tropical forest species with varying life-forms. Population dynamics from IPMs and IBMs matched a 34 year time series of data (albeit a snapshot of the life cycle for canopy trees) and highlight differences in life-history strategies between species. Specifically, the greater variation in growth rates within the two canopy species suggests an ability to respond to available resources, which in turn manifests as faster passage times and greater occupancy times in larger size classes. The framework presented here offers a novel and accessible approach to modelling the population dynamics of forest trees.
| S-EPMC5879618 | biostudies-literature
Project description:BackgroundStochastic effects can be important for the behavior of processes involving small population numbers, so the study of stochastic models has become an important topic in the burgeoning field of computational systems biology. However analysis techniques for stochastic models have tended to lag behind their deterministic cousins due to the heavier computational demands of the statistical approaches for fitting the models to experimental data. There is a continuing need for more effective and efficient algorithms. In this article we focus on the parameter inference problem for stochastic kinetic models of biochemical reactions given discrete time-course observations of either some or all of the molecular species.ResultsWe propose an algorithm for inference of kinetic rate parameters based upon maximum likelihood using stochastic gradient descent (SGD). We derive a general formula for the gradient of the likelihood function given discrete time-course observations. The formula applies to any explicit functional form of the kinetic rate laws such as mass-action, Michaelis-Menten, etc. Our algorithm estimates the gradient of the likelihood function by reversible jump Markov chain Monte Carlo sampling (RJMCMC), and then gradient descent method is employed to obtain the maximum likelihood estimation of parameter values. Furthermore, we utilize flux balance analysis and show how to automatically construct reversible jump samplers for arbitrary biochemical reaction models. We provide RJMCMC sampling algorithms for both fully observed and partially observed time-course observation data. Our methods are illustrated with two examples: a birth-death model and an auto-regulatory gene network. We find good agreement of the inferred parameters with the actual parameters in both models.ConclusionsThe SGD method proposed in the paper presents a general framework of inferring parameters for stochastic kinetic models. The method is computationally efficient and is effective for both partially and fully observed systems. Automatic construction of reversible jump samplers and general formulation of the likelihood gradient function makes our method applicable to a wide range of stochastic models. Furthermore our derivations can be useful for other purposes such as using the gradient information for parametric sensitivity analysis or using the reversible jump samplers for full Bayesian inference. The software implementing the algorithms is publicly available at http://cbcl.ics.uci.edu/sgd.
| S-EPMC2914651 | biostudies-literature
Project description:Stochastic population models are widely used to model phenomena in different areas such as cyber-physical systems, chemical kinetics, collective animal behaviour, and beyond. Quantitative analysis of stochastic population models easily becomes challenging due to the combinatorial number of possible states of the population. Moreover, while the modeller easily hypothesises the mechanistic aspects of the model, the quantitative parameters associated to these mechanistic transitions are difficult or impossible to measure directly. In this paper, we investigate how formal verification methods can aid parameter inference for population discrete-time Markov chains in a scenario where only a limited sample of population-level data measurements-sample distributions among terminal states-are available. We first discuss the parameter identifiability and uncertainty quantification in this setup, as well as how the existing techniques of formal parameter synthesis and Bayesian inference apply. Then, we propose and implement four different methods, three of which incorporate formal parameter synthesis as a pre-computation step. We empirically evaluate the performance of the proposed methods over four representative case studies. We find that our proposed methods incorporating formal parameter synthesis as a pre-computation step allow us to significantly enhance the accuracy, precision, and scalability of inference. Specifically, in the case of unidentifiable parameters, we accurately capture the subspace of parameters which is data-compliant at a desired confidence level.
| S-EPMC10642793 | biostudies-literature
Project description:The study of the initial phase of the adaptive immune response after first antigen encounter provides essential information on the magnitude and quality of the immune response. This phase is characterized by proliferation and dissemination of T cells in the lymphoid organs. Modeling and identifying the key features of this phenomenon may provide a useful tool for the analysis and prediction of the effects of immunization. This knowledge can be effectively exploited in vaccinology, where it is of interest to evaluate and compare the responses to different vaccine formulations. The objective of this paper is to construct a stochastic model based on branching process theory, for the dissemination network of antigen-specific CD4+ T cells. The devised model is validated on in vivo animal experimental data. The model presented has been applied to the vaccine immunization context making references to simple proliferation laws that take into account division, death and quiescence, but it can also be applied to any context where it is of interest to study the dynamic evolution of a population.
| S-EPMC4547705 | biostudies-literature
Project description:Mouse colorectal cancer (CRC) models generated by orthotopic microinjection of human CRC cell lines reproduce the pattern of lymphatic, haematological and transcoelomic spread but generate low metastatic efficiency. Our aim was to develop a new strategy that could increase the metastatic efficiency of these models. We used subcutaneous implantation of the human CRC cell lines HCT116 or SW48 prior to their orthotopic microinjection in the cecum of nude mice (SC+ORT). This subcutaneous preconditioning significantly enhanced metastatic dissemination. In the HCT116 model it increased the number and size of metastatic foci in lymph nodes, lung, liver and peritoneum, whereas, in the SW48 model, it induced a shift from non-metastatic to metastatic. In both models the number of apoptotic bodies in the primary tumour in the SC+ORT group was significantly reduced compared with that in the direct orthotopic injection (ORT) group. Moreover, in HCT116 tumours the number of keratin-positive tumour buddings and single epithelial cells increased at the invasion front in SC+ORT mice. In the SW48 tumour model, we observed a trend towards a higher number of tumour buds and single cells in the SC+ORT group but this did not reach statistical significance. At a molecular level, the enhanced metastatic efficiency observed in the HCT116 SC+ORT model was associated with an increase in AKT activation, VEGF-A overexpression and downregulation of ?1 integrin in primary tumour tissue, whereas, in SW48 SC+ORT mice, the level of expression of these proteins remained unchanged. In summary, subcutaneous preconditioning increased the metastatic dissemination of both orthotopic CRC models by increasing tumour cell survival and invasion at the tumour invasion front. This approach could be useful to simultaneously study the mechanisms of metastases and to evaluate anti-metastatic drugs against CRC.
| S-EPMC3944498 | biostudies-literature
Project description:BackgroundMonkeypox virus is gaining attention due to its severity and spread among people. This study sheds light on the modeling and forecasting of new monkeypox cases. Knowledge about the future situation of the virus using a more accurate time series and stochastic models is required for future actions and plans to cope with the challenge.MethodsWe conduct a side-by-side comparison of the machine learning approach with the traditional time series model. The multilayer perceptron model (MLP), a machine learning technique, and the Box-Jenkins methodology, also known as the ARIMA model, are used for classical modeling. Both methods are applied to the Monkeypox cumulative data set and compared using different model selection criteria such as root mean square error, mean square error, mean absolute error, and mean absolute percentage error.ResultsWith a root mean square error of 150.78, the monkeypox series follows the ARIMA (7,1,7) model among the other potential models. Comparatively, we use the multilayer perceptron (MLP) model, which employs the sigmoid activation function and has a different number of hidden neurons in a single hidden layer. The root mean square error of the MLP model, which uses a single input and ten hidden neurons, is 54.40, significantly lower than that of the ARIMA model. The actual confirmed cases versus estimated or fitted plots also demonstrate that the multilayer perceptron model has a better fit for the monkeypox data than the ARIMA model.Conclusions and recommendationWhen it comes to predicting monkeypox, the machine learning method outperforms the traditional time series. A better match can be achieved in future studies by applying the extreme learning machine model (ELM), support vector machine (SVM), and some other methods with various activation functions. It is thus concluded that the selected data provide a real picture of the virus. If the situations remain the same, governments and other stockholders should ensure the follow-up of Standard Operating Procedures (SOPs) among the masses, as the trends will continue rising in the upcoming 10 days. However, governments should take some serious interventions to cope with the virus.LimitationIn the ARIMA models selected for forecasting, we did not incorporate the effect of covariates such as the effect of net migration of monkeypox virus patients, government interventions, etc.
| S-EPMC9659136 | biostudies-literature
Project description:The most common evolutionary events at the molecular level are single-base substitutions, as well as insertions and deletions (indels) of short DNA segments. A large body of research has been devoted to develop probabilistic substitution models and to infer their parameters using likelihood and Bayesian approaches. In contrast, relatively little has been done to model indel dynamics, probably due to the difficulty in writing explicit likelihood functions. Here, we contribute to the effort of modeling indel dynamics by presenting SpartaABC, an approximate Bayesian computation (ABC) approach to infer indel parameters from sequence data (either aligned or unaligned). SpartaABC circumvents the need to use an explicit likelihood function by extracting summary statistics from simulated sequences. First, summary statistics are extracted from the input sequence data. Second, SpartaABC samples indel parameters from a prior distribution and uses them to simulate sequences. Third, it computes summary statistics from the simulated sets of sequences. By computing a distance between the summary statistics extracted from the input and each simulation, SpartaABC can provide an approximation to the posterior distribution of indel parameters as well as point estimates. We study the performance of our methodology and show that it provides accurate estimates of indel parameters in simulations. We next demonstrate the utility of SpartaABC by studying the impact of alignment errors on the inference of positive selection. A C?++?program implementing SpartaABC is freely available in http://spartaabc.tau.ac.il.
| S-EPMC5438127 | biostudies-literature