A Bayesian Semi-parametric Approach for the Differential Analysis of Sequence Counts Data.
Ontology highlight
ABSTRACT:
SUBMITTER: Guindani M
PROVIDER: S-EPMC4017673 | biostudies-literature | 2014 Apr
REPOSITORIES: biostudies-literature
Publications
A Bayesian Semi-parametric Approach for the Differential Analysis of Sequence Counts Data.
Guindani Michele M Sepúlveda Nuno N Paulino Carlos Daniel CD Müller Peter P
Journal of the Royal Statistical Society. Series C, Applied statistics 20140401 3
Data obtained using modern sequencing technologies are often summarized by recording the frequencies of observed sequences. Examples include the analysis of T cell counts in immunological research and studies of gene expression based on counts of RNA fragments. In both cases the items being counted are sequences, of proteins and base pairs, respectively. The resulting sequence-abundance distribution is usually characterized by overdispersion. We propose a Bayesian semi-parametric approach to imp ...[more]
Similar Datasets
Project description:Differential co-expression analysis is an emerging strategy for characterizing disease related dysregulation of gene expression regulatory networks. Given pre-defined sets of biological samples, such analysis aims at identifying genes that are co-expressed in one, but not in the other set of samples.We developed a novel probabilistic framework for jointly uncovering contexts (i.e. groups of samples) with specific co-expression patterns, and groups of genes with different co-expression patterns across such contexts. In contrast to current clustering and bi-clustering procedures, the implicit similarity measure in this model used for grouping biological samples is based on the clustering structure of genes within each sample and not on traditional measures of gene expression level similarities. Within this framework, biological samples with widely discordant expression patterns can be placed in the same context as long as the co-clustering structure of genes is concordant within these samples. To the best of our knowledge, this is the first method to date for unsupervised differential co-expression analysis in this generality. When applied to the problem of identifying molecular subtypes of breast cancer, our method identified reproducible patterns of differential co-expression across several independent expression datasets. Sample groupings induced by these patterns were highly informative of the disease outcome. Expression patterns of differentially co-expressed genes provided new insights into the complex nature of the ERalpha regulatory network.We demonstrated that the use of the co-clustering structure as the similarity measure in the unsupervised analysis of sample gene expression profiles provides valuable information about expression regulatory networks.
| S-EPMC2876132 | biostudies-literature
Project description:In the U.S., the Centers for Medicare and Medicaid Services uses 30-day readmission, following hospitalization, as a proxy outcome to monitor quality of care. These efforts generally focus on treatable health conditions, such as pneumonia and heart failure. Expanding quality of care systems to monitor conditions for which treatment options are limited or non-existent, such as pancreatic cancer, is challenging because of the non-trivial force of mortality; 30-day mortality for pancreatic cancer is approximately 30%. In the statistical literature, data that arise when the observation of the time to some non-terminal event is subject to some terminal event are referred to as 'semi-competing risks data'. Given such data, scientific interest may lie in at least one of three areas: (i) estimation/inference for regression parameters, (ii) characterization of dependence between the two events, and (iii) prediction given a covariate profile. Existing statistical methods focus almost exclusively on the first of these; methods are sparse or non-existent, however, when interest lies with understanding dependence and performing prediction. In this paper we propose a Bayesian semi-parametric regression framework for analyzing semi-competing risks data that permits the simultaneous investigation of all three of the aforementioned scientific goals. Characterization of the induced posterior and posterior predictive distributions is achieved via an efficient Metropolis-Hastings-Green algorithm, which has been implemented in an R package. The proposed framework is applied to data on 16,051 individuals diagnosed with pancreatic cancer between 2005-2008, obtained from Medicare Part A. We found that increased risk for readmission is associated with a high comorbidity index, a long hospital stay at initial hospitalization, non-white race, male, and discharge to home care.
| S-EPMC4427057 | biostudies-literature
Project description:BACKGROUND:Identifying differentially abundant features between different experimental groups is a common goal for many metabolomics and proteomics studies. However, analyzing data from mass spectrometry (MS) is difficult because the data may not be normally distributed and there is often a large fraction of zero values. Although several statistical methods have been proposed, they either require the data normality assumption or are inefficient. RESULTS:We propose a new semi-parametric differential abundance analysis (SDA) method for metabolomics and proteomics data from MS. The method considers a two-part model, a logistic regression for the zero proportion and a semi-parametric log-linear model for the possibly non-normally distributed non-zero values, to characterize data from each feature. A kernel-smoothed likelihood method is developed to estimate model coefficients and a likelihood ratio test is constructed for differential abundant analysis. The method has been implemented into an R package, SDAMS, which is available at https://www.bioconductor.org/packages/release/bioc/html/SDAMS.html . CONCLUSION:By introducing the two-part semi-parametric model, SDA is able to handle both non-normally distributed data and large fraction of zero values in a MS dataset. It also allows for adjustment of covariates. Simulations and real data analyses demonstrate that SDA outperforms existing methods.
| S-EPMC6798423 | biostudies-literature
Project description:Impairment caused by Parkinson's disease (PD) is multidimensional (e.g., sensoria, functions, and cognition) and progressive. Its multidimensional nature precludes a single outcome to measure disease progression. Clinical trials of PD use multiple categorical and continuous longitudinal outcomes to assess the treatment effects on overall improvement. A terminal event such as death or dropout can stop the follow-up process. Moreover, the time to the terminal event may be dependent on the multivariate longitudinal measurements. In this article, we consider a joint random-effects model for the correlated outcomes. A multilevel item response theory model is used for the multivariate longitudinal outcomes and a parametric accelerated failure time model is used for the failure time because of the violation of proportional hazard assumption. These two models are linked via random effects. The Bayesian inference via MCMC is implemented in 'BUGS' language. Our proposed method is evaluated by a simulation study and is applied to DATATOP study, a motivating clinical trial to determine if deprenyl slows the progression of PD.
| S-EPMC3947121 | biostudies-literature
Project description:The thermal stability of proteins can be altered when they interact with small molecules, other biomolecules or are subject to post-translation modifications. Thus monitoring the thermal stability of proteins under various cellular perturbations can provide insights into protein function, as well as potentially determine drug targets and off-targets. Thermal proteome profiling is a highly multiplexed mass-spectrommetry method for monitoring the melting behaviour of thousands of proteins in a single experiment. In essence, thermal proteome profiling assumes that proteins denature upon heating and hence become insoluble. Thus, by tracking the relative solubility of proteins at sequentially increasing temperatures, one can report on the thermal stability of a protein. Standard thermodynamics predicts a sigmoidal relationship between temperature and relative solubility and this is the basis of current robust statistical procedures. However, current methods do not model deviations from this behaviour and they do not quantify uncertainty in the melting profiles. To overcome these challenges, we propose the application of Bayesian functional data analysis tools which allow complex temperature-solubility behaviours. Our methods have improved sensitivity over the state-of-the art, identify new drug-protein associations and have less restrictive assumptions than current approaches. Our methods allows for comprehensive analysis of proteins that deviate from the predicted sigmoid behaviour and we uncover potentially biphasic phenomena with a series of published datasets.
| S-EPMC8241860 | biostudies-literature
Project description:Understanding sub-cellular protein localisation is an essential component in the analysis of context specific protein function. Recent advances in quantitative mass-spectrometry (MS) have led to high resolution mapping of thousands of proteins to sub-cellular locations within the cell. Novel modelling considerations to capture the complex nature of these data are thus necessary. We approach analysis of spatial proteomics data in a non-parametric Bayesian framework, using K-component mixtures of Gaussian process regression models. The Gaussian process regression model accounts for correlation structure within a sub-cellular niche, with each mixture component capturing the distinct correlation structure observed within each niche. The availability of marker proteins (i.e. proteins with a priori known labelled locations) motivates a semi-supervised learning approach to inform the Gaussian process hyperparameters. We moreover provide an efficient Hamiltonian-within-Gibbs sampler for our model. Furthermore, we reduce the computational burden associated with inversion of covariance matrices by exploiting the structure in the covariance matrix. A tensor decomposition of our covariance matrices allows extended Trench and Durbin algorithms to be applied to reduce the computational complexity of inversion and hence accelerate computation. We provide detailed case-studies on Drosophila embryos and mouse pluripotent embryonic stem cells to illustrate the benefit of semi-supervised functional Bayesian modelling of the data.
| S-EPMC7613899 | biostudies-literature
Project description:Spatial statistics has seen rapid application in many fields, especially epidemiology and public health. Many studies, nonetheless, make limited use of the geographical location information and also usually assume that the covariates, which are related to the response variable, have linear effects. We develop a Bayesian semi-parametric regression model for HIV prevalence data. Model estimation and inference is based on fully Bayesian approach via Markov Chain Monte Carlo (McMC). The model is applied to HIV prevalence data among men in Kenya, derived from the Kenya AIDS indicator survey, with n?=?3,662. Past studies have concluded that HIV infection has a nonlinear association with age. In this study a smooth function based on penalized regression splines is used to estimate this nonlinear effect. Other covariates were assumed to have a linear effect. Spatial references to the counties were modeled as both structured and unstructured spatial effects. We observe that circumcision reduces the risk of HIV infection. The results also indicate that men in the urban areas were more likely to be infected by HIV as compared to their rural counterpart. Men with higher education had the lowest risk of HIV infection. A nonlinear relationship between HIV infection and age was established. Risk of HIV infection increases with age up to the age of 40 then declines with increase in age. Men who had STI in the last 12 months were more likely to be infected with HIV. Also men who had ever used a condom were found to have higher likelihood to be infected by HIV. A significant spatial variation of HIV infection in Kenya was also established. The study shows the practicality and flexibility of Bayesian semi-parametric regression model in analyzing epidemiological data.
| S-EPMC4111556 | biostudies-literature
Project description:BackgroundIntimate partner violence (IPV) is an important public health problem with health and socioeconomic consequences and is endemic in Namibia. Studies assessing risk factors for IPV often use logistic and Poisson regression without geographical location information and spatial effects. We used a Bayesian spatial semi-parametric regression model to determine the risk factors for IPV in Namibia; assess the non-linear effects of age difference between partners and determine spatial effects in the different regions on IPV prevalence.MethodsWe used the couples' dataset of the 2013-2014 Namibia Demographic and Health Survey (DHS) obtained on request from Measure DHS. The DHS domestic violence module included 2226 women. We generated a binary variable measuring IPV from the questions "ever experienced physical, sexual or emotional violence?" Covariates included respondent's educational level, age, couples' age difference, place of residence and partner's educational level. All estimation was done with the full Bayesian approach using R version 3.5.2 implementing the R2BayesX package.ResultsIPV country prevalence was 33.3% (95% CI = 30.1-36.5%); Kavango had the highest [50.6% (95% CI = 41.2-60.1%)] and Oshana the lowest [11.5% (95% CI = 3.2-19.9%)] regional prevalence. IPV prevalence was highest among teenagers [60.8% (95% CI = 36.9-84.7%)]). The spatial semi-parametric model used for adjusted results controlled for regional spatial effects, respondent's age, age difference, respondent's years of education, residence, wealth, and education levels. Women with higher education were 50% less likely to experience IPV [aOR: 0.46, 95% CI = 0.23-0.87]. For non-linear effects, the risk of IPV was high for women ≥ 5 years older or ≥ 25 years younger than their partners. Younger and older women had higher risks of IPV than those between 25 and 45 years. For spatial variation of IPV prevalence, northern regions had low spatial effects while western regions had very high spatial effects.ConclusionThe prevalence of IPV among Namibia women was high especially among teenagers, with higher educational levels being protective. The risk of IPV was lower in rural than urban areas and higher with wide partner age differences. Interventions and policies for IPV prevention in Namibia are needed for couples with wide age differences as well as for younger women, women with lower educational attainment and in urban and western regions.
| S-EPMC8340378 | biostudies-literature
Project description:Early evaluation of new drug entities for their potential to cause mitochondrial dysfunction is becoming an important task for drug development. Multi-parametric high-content screening (mp-HCS) of mitochondrial toxicity holds promise as a lead in-vitro strategy for drug testing and safety evaluations. In this study, we have developed a mp-HCS and multi-parametric data analysis scheme for assessing cell responses to induced mitochondrial perturbation. The mp-HCS measurements are shown to be robust enough to allow for quantitative comparison of biological systems with different metabolic pathways simulated by alteration of growth media. Substitution of medium glucose for galactose sensitized cells to drug action and revealed novel response parameters. Each compound was quantitatively characterized according to induced phenotypic changes of cell morphology and functionality measured by fluorescent biomarkers for mitochondrial activity, plasma membrane permeability, and nuclear morphology. Descriptors of drug effects were established by generation of a SCRIT (Specialized-Cell-Response-to-Induced-Toxicity) vector, consisting of normalized statistical measures of each parameter at each dose and growth condition. The dimensionality of SCRIT vectors depends on the number of parameters chosen, which in turn depends on the hypothesis being tested. Specifically, incorporation of three parameters of response into SCRIT vectors enabled clustering of 84 training compounds with known pharmacological and toxicological activities according to the degree of toxicity and mitochondrial involvement. Inclusion of 6 parameters enabled the resolution of more subtle differences between compounds within a common therapeutic class; scoring enabled a ranking of statins in direct agreement with clinical outcomes. Comparison of drug-induced changes required variations in glucose for separation of mitochondrial dysfunction from other types of cytotoxicity. These results also demonstrate that the number of drugs in a training set, the choice of parameters used in analysis, and statistical measures are fundamental for specific hypothesis testing and assessment of quantitative phenotypic differences.
| S-EPMC3471932 | biostudies-literature
Project description:In the analysis of semi-competing risks data interest lies in estimation and inference with respect to a so-called non-terminal event, the observation of which is subject to a terminal event. Multi-state models are commonly used to analyse such data, with covariate effects on the transition/intensity functions typically specified via the Cox model and dependence between the non-terminal and terminal events specified, in part, by a unit-specific shared frailty term. To ensure identifiability, the frailties are typically assumed to arise from a parametric distribution, specifically a Gamma distribution with mean 1.0 and variance, say, σ2. When the frailty distribution is misspecified, however, the resulting estimator is not guaranteed to be consistent, with the extent of asymptotic bias depending on the discrepancy between the assumed and true frailty distributions. In this paper, we propose a novel class of transformation models for semi-competing risks analysis that permit the non-parametric specification of the frailty distribution. To ensure identifiability, the class restricts to parametric specifications of the transformation and the error distribution; the latter are flexible, however, and cover a broad range of possible specifications. We also derive the semi-parametric efficient score under the complete data setting and propose a non-parametric score imputation method to handle right censoring; consistency and asymptotic normality of the resulting estimators is derived and small-sample operating characteristics evaluated via simulation. Although the proposed semi-parametric transformation model and non-parametric score imputation method are motivated by the analysis of semi-competing risks data, they are broadly applicable to any analysis of multivariate time-to-event outcomes in which a unit-specific shared frailty is used to account for correlation. Finally, the proposed model and estimation procedures are applied to a study of hospital readmission among patients diagnosed with pancreatic cancer.
| S-EPMC5400113 | biostudies-literature