Dataset Information

Multi-experiment nonlinear mixed effect modeling of single-cell translation kinetics after transfection.

ABSTRACT: Single-cell time-lapse studies have advanced the quantitative understanding of cellular pathways and their inherent cell-to-cell variability. However, parameters retrieved from individual experiments are model dependent and their estimation is limited, if based on solely one kind of experiment. Hence, methods to integrate data collected under different conditions are expected to improve model validation and information content. Here we present a multi-experiment nonlinear mixed effect modeling approach for mechanistic pathway models, which allows the integration of multiple single-cell perturbation experiments. We apply this approach to the translation of green fluorescent protein after transfection using a massively parallel read-out of micropatterned single-cell arrays. We demonstrate that the integration of data from perturbation experiments allows the robust reconstruction of cell-to-cell variability, i.e., parameter densities, while each individual experiment provides insufficient information. Indeed, we show that the integration of the datasets on the population level also improves the estimates for individual cells by breaking symmetries, although each of them is only measured in one experiment. Moreover, we confirmed that the suggested approach is robust with respect to batch effects across experimental replicates and can provide mechanistic insights into the nature of batch effects. We anticipate that the proposed multi-experiment nonlinear mixed effect modeling approach will serve as a basis for the analysis of cellular heterogeneity in single-cell dynamics.

SUBMITTER: Frohlich F

PROVIDER: S-EPMC6288153 | biostudies-literature | 2019

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Ensuring future food security for a growing population while climate change and urban sprawl put pressure on agricultural land will require sustainable intensification of current farming practices. For the crop breeder this means producing higher crop yields with less resources due to greater environmental stresses. While easy gains in crop yield have been made mostly "above ground," little progress has been made "below ground"; and yet it is these root system traits that can improve productivity and resistance to drought stress. Wheat pre-breeders use soil coring and core-break counts to phenotype root architecture traits, with data collected on rooting density for hundreds of genotypes in small increments of depth. The measured densities are both large datasets and highly variable even within the same genotype, hence, any rigorous, comprehensive statistical analysis of such complex field data would be technically challenging. Traditionally, most attributes of the field data are therefore discarded in favor of simple numerical summary descriptors which retain much of the high variability exhibited by the raw data. This poses practical challenges: although plant scientists have established that root traits do drive resource capture in crops, traits that are more randomly (rather than genetically) determined are difficult to breed for. In this paper we develop a hierarchical nonlinear mixed modeling approach that utilizes the complete field data for wheat genotypes to fit, under the Bayesian paradigm, an "idealized" relative intensity function for the root distribution over depth. Our approach was used to determine heritability: how much of the variation between field samples was purely random vs. being mechanistically driven by the plant genetics? Based on the genotypic intensity functions, the overall heritability estimate was 0.62 (95% Bayesian confidence interval was 0.52 to 0.71). Despite root count profiles that were statistically very noisy, our approach led to denoised profiles which exhibited rigorously discernible phenotypic traits. Profile-specific traits could be representative of a genotype, and thus, used as a quantitative tool to associate phenotypic traits with specific genotypes. This would allow breeders to select for whole root system distributions appropriate for sustainable intensification, and inform policy for mitigating crop yield risk and food insecurity.

Project description:BackgroundAlthough in vivo studies of internal exposure to hazardous substances have been carried out for many years, there is room for progress to improve their informative value while adhering to the four R's: replacement, reduction, refinement, and responsibility rule.ObjectivesThe objective of the study was to illustrate how toxicokinetic (TK) study design and data analysis can be implemented under the 4R rule to plan a chronic dosage regimen for investigating TK/toxicodynamic (TD) relationships.MethodsThe intravenous (IV) and oral serum concentrations of eight hazardous environmental contaminants including 1,1-Dichloro-2,2-bis(p-chlorophenyl)ethylene (pp'DDE), ß-Hexachlorocyclohexane (β-HCH), hexachlorobenzene (HCB), 2,2'4,4'-tetrabromodiphenyl ether (BDE-47), perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), di(2ethylhexyl)phthalate (DEHP), and bisphenol S (BPS) were obtained after mixture dosing in rabbits using a sparse sampling design. Data were comprehensively analyzed using nonlinear mixed effect (NLME) modeling.ResultsThe short persistence of BPS and of the DEHP metabolite (mono-2-ethylhexyl phthalate), reflected by their mean residence times (MRT) of a few hours, was due to their efficient clearance (CL, 3.2 and 0.47L/kg/h). The longer MRT of the other compounds (1-48 d) resulted either from their extremely low clearance (lower than 0.01L/kg/h for PFOA and PFOS) or from their very large volume of distribution (VSS) ranging from 33 to 45L/kg. Estimates of CL, VSS, and bioavailability were used to compute the oral loading and daily maintenance doses required to attain a nominal steady-state serum concentration of 1 ng/mL. Simulations with the NLME model were applied to predict the serum concentration profile and to contrast the differential rates of accumulation in the central vs. peripheral compartments.ConclusionNLME modeling of the IV and oral TK of hazardous environmental contaminants, in rabbits while fulfilling the 4R rule, was able to provide the physiological basis for interspecies extrapolation of exposure rates in a TK/TD approach to risk assessment. https://doi.org/10.1289/EHP8957.

Project description:In synapses, calcium is required for modulating synaptic transmission, plasticity, synaptogenesis, and synaptic pruning. The regulation of calcium dynamics within neurons involves cellular mechanisms such as synaptically activated channels and pumps, calcium buffers, and calcium sequestrating organelles. Many experimental studies tend to focus on only one or a small number of these mechanisms, as technical limitations make it difficult to observe all features at once. Computational modeling enables incorporation of many of these properties together, allowing for more complete and integrated studies. However, the scale of existing detailed models is often limited to synaptic and dendritic compartments as the computational burden rapidly increases when these models are integrated in cellular or network level simulations. In this article we present a computational model of calcium dynamics at the postsynaptic spine of a CA1 pyramidal neuron, as well as a methodology that enables its implementation in multi-scale, large-scale simulations. We first present a mechanistic model that includes individually validated models of various components involved in the regulation of calcium at the spine. We validated our mechanistic model by comparing simulated calcium levels to experimental data found in the literature. We performed additional simulations with the mechanistic model to determine how the simulated calcium activity varies with respect to presynaptic-postsynaptic stimulation intervals and spine distance from the soma. We then developed an input-output (IO) model that complements the mechanistic calcium model and provide a computationally efficient representation for use in larger scale modeling studies; we show the performance of the IO model compared to the mechanistic model in terms of accuracy and speed. The models presented here help achieve two objectives. First, the mechanistic model provides a comprehensive platform to describe spine calcium dynamics based on individual contributing factors. Second, the IO model is trained on the main dynamical features of the mechanistic model and enables nonlinear spine calcium modeling on the cell and network level simulation scales. Utilizing both model representations provide a multi-level perspective on calcium dynamics, originating from the molecular interactions at spines and propagating the effects to higher levels of activity involved in network behavior.

Dataset Information

Multi-experiment nonlinear mixed effect modeling of single-cell translation kinetics after transfection.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets