Decorate: differential epigenetic correlation test.
Ontology highlight
ABSTRACT: Motivation
Identifying correlated epigenetic features and finding differences in correlation between individuals with disease compared to controls can give novel insight into disease biology. This framework has been successful in analysis of gene expression data, but application to epigenetic data has been limited by the computational cost, lack of scalable software and lack of robust statistical tests.Results
Decorate, differential epigenetic correlation test, identifies correlated epigenetic features and finds clusters of features that are differentially correlated between two or more subsets of the data. The software scales to genome-wide datasets of epigenetic assays on hundreds of individuals. We apply decorate to four large-scale datasets of DNA methylation, ATAC-seq and histone modification ChIP-seq.Availability and implementation
decorate R package is available from https://github.com/GabrielHoffman/decorate.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Hoffman GE
PROVIDER: S-EPMC7203736 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:When a statistical test is repeatedly applied to rows of a data matrix, correlations among data rows will give rise to correlations among corresponding test statistics. We investigate the relationship between test-statistic correlation and data-row correlation and discuss its implications.
| S-EPMC7723344 | biostudies-literature
Project description:Testing correlation structures has attracted extensive attention in the literature due to both its importance in real applications and several major theoretical challenges. The aim of this paper is to develop a general framework of testing correlation structures for the one-, two-, and multiple sample testing problems under a high-dimensional setting when both the sample size and data dimension go to infinity. Our test statistics are designed to deal with both the dense and sparse alternatives. We systematically investigate the asymptotic null distribution, power function, and unbiasedness of each test statistic. Theoretically, we make great efforts to deal with the non-independency of all random matrices of the sample correlation matrices. We use simulation studies and real data analysis to illustrate the versatility and practicability of our test statistics.
| S-EPMC6709985 | biostudies-literature
Project description:With the availability of an increasing number of 3D structures of bacteriophage components, combined with powerful in silico predictive tools, it has become possible to decipher the structural assembly and functionality of phage adhesion devices. In the current study, we examined 113 members of the 936 group of lactococcal siphophages, and identified a number of Carbohydrate Binding Modules (CBMs) in the neck passage structure and major tail protein, on top of evolved Dit proteins, as recently reported by us. The binding ability of such CBM-containing proteins was assessed through the construction of green fluorescent protein fusion proteins and subsequent binding assays. Two CBMs, one from the phage tail and another from the neck, demonstrated definite binding to their phage-specific host. Bioinformatic analysis of the structural proteins of 936 phages reveals that they incorporate binding modules which exhibit structural homology to those found in other lactococcal phage groups and beyond, indicating that phages utilize common structural "bricks" to enhance host binding capabilities. The omnipresence of CBMs in Siphophages supports their beneficial role in the infection process, as they can be combined in various ways to form appendages with different shapes and functionalities, ensuring their success in host detection in their respective ecological niches.
| S-EPMC6669499 | biostudies-literature
Project description:Distance correlation has gained much recent attention in the data science community: the sample statistic is straightforward to compute and asymptotically equals zero if and only if independence, making it an ideal choice to discover any type of dependency structure given sufficient sample size. One major bottleneck is the testing process: because the null distribution of distance correlation depends on the underlying random variables and metric choice, it typically requires a permutation test to estimate the null and compute the p-value, which is very costly for large amount of data. To overcome the difficulty, in this paper we propose a chi-square test for distance correlation. Method-wise, the chi-square test is non-parametric, extremely fast, and applicable to bias-corrected distance correlation using any strong negative type metric or characteristic kernel. The test exhibits a similar testing power as the standard permutation test, and can be utilized for K-sample and partial testing. Theory-wise, we show that the underlying chi-square distribution well approximates and dominates the limiting null distribution in upper tail, prove the chi-square test can be valid and universally consistent for testing independence, and establish a testing power inequality with respect to the permutation test.
| S-EPMC9191842 | biostudies-literature
Project description:In the context of multivariate multilevel data analysis, this paper focuses on the multivariate linear mixed-effects model, including all the correlations between the random effects when the dimensional residual terms are assumed uncorrelated. Using the EM algorithm, we suggest more general expressions of the model's parameters estimators. These estimators can be used in the framework of the multivariate longitudinal data analysis as well as in the more general context of the analysis of multivariate multilevel data. By using a likelihood ratio test, we test the significance of the correlations between the random effects of two dependent variables of the model, in order to investigate whether or not it is useful to model these dependent variables jointly. Simulation studies are done to assess both the parameter recovery performance of the EM estimators and the power of the test. Using two empirical data sets which are of longitudinal multivariate type and multivariate multilevel type, respectively, the usefulness of the test is illustrated.
| S-EPMC4990185 | biostudies-literature
Project description:BackgroundAn important emerging trend in the analysis of microarray data is to incorporate known pathway information a priori. Expression level "summaries" for pathways, obtained from the expression data for the genes constituting the pathway, permit the inclusion of pathway information, reduce the high dimensionality of microarray data, and have the power to elucidate gene-interaction dependencies which are not already accounted for through known pathway identification.ResultsWe present a novel method for the analysis of microarray data that identifies joint differential expression in gene-pathway pairs. This method takes advantage of known gene pathway memberships to compute a summary expression level for each pathway as a whole. Correlations between the pathway expression summary and the expression levels of genes not already known to be associated with the pathway provide clues to gene interaction dependencies that are not already accounted for through known pathway identification, and statistically significant differences between gene-pathway correlations in phenotypically different cells (e.g., where the expression level of a single gene and a given pathway summary correlate strongly in normal cells but weakly in tumor cells) may indicate biologically relevant gene-pathway interactions. Here, we detail the methodology and present the results of this method applied to two gene-expression datasets, identifying gene-pathway pairs which exhibit differential joint expression by phenotype.ConclusionThe method described herein provides a means by which interactions between large numbers of genes may be identified by incorporating known pathway information to reduce the dimensionality of gene interactions. The method is efficient and easily applied to data sets of ~102 arrays. Application of this method to two publicly-available cancer data sets yields suggestive and promising results. This method has the potential to complement gene-at-a-time analysis techniques for microarray analysis by indicating relationships between pathways and genes that have not previously been identified and which may play a role in disease.
| S-EPMC2613418 | biostudies-literature
Project description:Variation in cancer risk among somatic tissues has been attributed to variations in the underlying rate of stem cell division. For a given tissue type, variable cancer risk between individuals is thought to be influenced by extrinsic factors which modulate this rate of stem cell division. To date, no molecular mitotic clock has been developed to approximate the number of stem cell divisions in a tissue of an individual and which is correlated with cancer risk.Here, we integrate mathematical modeling with prior biological knowledge to construct a DNA methylation-based age-correlative model which approximates a mitotic clock in both normal and cancer tissue. By focusing on promoter CpG sites that localize to Polycomb group target genes that are unmethylated in 11 different fetal tissue types, we show that increases in DNA methylation at these sites defines a tick rate which correlates with the estimated rate of stem cell division in normal tissues. Using matched DNA methylation and RNA-seq data, we further show that it correlates with an expression-based mitotic index in cancer tissue. We demonstrate that this mitotic-like clock is universally accelerated in cancer, including pre-cancerous lesions, and that it is also accelerated in normal epithelial cells exposed to a major carcinogen.Unlike other epigenetic and mutational clocks or the telomere clock, the epigenetic clock proposed here provides a concrete example of a mitotic-like clock which is universally accelerated in cancer and precancerous lesions.
| S-EPMC5046977 | biostudies-literature
Project description:In this work, we developed a robust permutation test for the concordance correlation coefficient (ρc ) for testing the general hypothesis H0 : ρc = ρc(0) . The proposed test is based on an appropriately studentized statistic. Theoretically, the test is proven to be asymptotically valid in the general setting when two paired variables are uncorrelated but dependent. This desired property was demonstrated across a range of distributional assumptions and sample sizes in simulation studies, where the test exhibits robust type I error control in all settings tested, even when the sample size is small. We demonstrated the application of this test in two real world examples across cardiac output measurements and endocardiographic imaging.
| S-EPMC8359348 | biostudies-literature
Project description:Shotgun proteomics techniques infer the presence and quantity of proteins using peptide proxies produced by cleavage of the proteome with a protease. Most protein quantitation strategies assume that multiple peptides derived from a protein will behave quantitatively similar across treatment groups, but this assumption may be false due to (1) heterogeneous proteoforms and (2) technical artifacts. Here we describe a strategy called peptide correlation analysis (PeCorA) that detects quantitative disagreements between peptides mapped to the same protein. PeCorA fits linear models to assess whether a peptide's change across treatment groups differs from all other peptides assigned to the same protein. PeCorA revealed that ∼15% of proteins in a mouse microglia stress data set contain at least one discordant peptide. Inspection of the discordant peptides shows the utility of PeCorA for the direct and indirect detection of regulated post-translational modifications (PTMs) and also for the discovery of poorly quantified peptides. The exclusion of poorly quantified peptides before protein quantity summarization decreased false-positives in a benchmark data set. Finally, PeCorA suggests that the inactive isoform of prothrombin, a coagulation cascade protease, is more abundant in plasma from COVID-19 patients relative to non-COVID-19 controls. PeCorA is freely available as an R package that works with arbitrary tables of quantified peptides.
| S-EPMC8592057 | biostudies-literature
Project description:Adrenomedullin (AM) and proadrenomedullin N-terminal 20 peptide (PAMP) are secretory hormones, but it is not unusual to find them in intracellular compartments. Using yeast-2 hybrid technology, we found interactions between AM and several microtubule-associated proteins (MAPs), and between PAMP and tubulin. Expression of fluorescent-tagged AM and PAMP as well as immunofluorescence for the native peptides showed a complete decoration of the microtubules and colocalization with other MAPs. PAMP, but not AM, bound to tubulin in vitro and destabilized tubulin polymerization. Down-regulation of the gene coding for both AM and PAMP through small interfering RNA technology resulted in morphological changes, microtubule stabilization, increase in posttranslational modifications of tubulin such as acetylation and detyrosination, reduction in cell motility, and partial arrest at the G2 phase of the cell cycle, when compared with cells transfected with the same vector carrying a scrambled sequence. These results show that PAMP is a novel MAP, whereas AM may be exerting more subtle effects in regulating cytoskeleton function.
| S-EPMC2408810 | biostudies-literature