A novel joint sparse partial correlation method for estimating group functional networks.
Ontology highlight
ABSTRACT:
SUBMITTER: Liang X
PROVIDER: S-EPMC8022278 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In this paper, we propose a computationally efficient approach -space(Sparse PArtial Correlation Estimation)- for selecting non-zero partial correlations under the high-dimension-low-sample-size setting. This method assumes the overall sparsity of the partial correlation matrix and employs sparse regression techniques for model fitting. We illustrate the performance of space by extensive simulation studies. It is shown that space performs well in both non-zero partial correlation selection and the identification of hub variables, and also outperforms two existing methods. We then apply space to a microarray breast cancer data set and identify a set of hub genes which may provide important insights on genetic regulatory networks. Finally, we prove that, under a set of suitable assumptions, the proposed procedure is asymptotically consistent in terms of model selection and parameter estimation.
| S-EPMC2770199 | biostudies-literature
Project description:Histone modifications are known to play an important role in the regulation of transcription. While individual modifications have received much attention in genome-wide analyses, little is known about their relationships. Some authors have built Bayesian networks of modifications, however most often they have used discretized data, and relied on unrealistic assumptions such as the absence of feedback mechanisms or hidden confounding factors. Here, we propose to infer undirected networks based on partial correlations between histone modifications. Within the partial correlation framework, correlations among two variables are controlled for associations induced by the other variables. Partial correlation networks thus focus on direct associations of histone modifications. We apply this methodology to data in CD4+ cells. The resulting network is well supported by common knowledge. When pairs of modifications show a large difference between their correlation and their partial correlation, a potential confounding factor is identified and provided as explanation. Data from different cell types (IMR90, H1) is also exploited in the analysis to assess the stability of the networks. The results are remarkably similar across cell types. Based on this observation, the networks from the three cell types are integrated into a consensus network to increase robustness. The data and the results discussed in the manuscript can be found, together with code, on http://spcn.molgen.mpg.de/index.html.
| S-EPMC3764007 | biostudies-literature
Project description:MOTIVATION:MicroRNAs (miRNAs) are important non-coding post-transcriptional regulators that are involved in many biological processes and human diseases. Individual miRNAs may regulate hundreds of genes, giving rise to a complex gene regulatory network in which transcripts carrying miRNA binding sites act as competing endogenous RNAs (ceRNAs). Several methods for the analysis of ceRNA interactions exist, but these do often not adjust for statistical confounders or address the problem that more than one miRNA interacts with a target transcript. RESULTS:We present SPONGE, a method for the fast construction of ceRNA networks. SPONGE uses 'multiple sensitivity correlation', a newly defined measure for which we can estimate a distribution under a null hypothesis. SPONGE can accurately quantify the contribution of multiple miRNAs to a ceRNA interaction with a probabilistic model that addresses previously neglected confounding factors and allows fast P-value calculation, thus outperforming existing approaches. We applied SPONGE to paired miRNA and gene expression data from The Cancer Genome Atlas for studying global effects of miRNA-mediated cross-talk. Our results highlight already established and novel protein-coding and non-coding ceRNAs which could serve as biomarkers in cancer. AVAILABILITY AND IMPLEMENTATION:SPONGE is available as an R/Bioconductor package (doi: 10.18129/B9.bioc.SPONGE). SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
| S-EPMC6612827 | biostudies-literature
Project description:Coordinated variations in brain morphology (e.g., cortical thickness) across individuals have been widely used to infer large-scale population brain networks. These structural correlation networks (SCNs) have been shown to reflect synchronized maturational changes in connected brain regions. Further, evidence suggests that SCNs, to some extent, reflect both anatomical and functional connectivity and hence provide a complementary measure of brain connectivity in addition to diffusion weighted networks and resting-state functional networks. Although widely used to study between-group differences in network properties, SCNs are inferred only at the group-level using brain morphology data from a set of participants, thereby not providing any knowledge regarding how the observed differences in SCNs are associated with individual behavioral, cognitive and disorder states. In the present study, we introduce two novel distance-based approaches to extract information regarding individual differences from the group-level SCNs. We applied the proposed approaches to a moderately large dataset (n=100) consisting of individuals with fragile X syndrome (FXS; n=50) and age-matched typically developing individuals (TD; n=50). We tested the stability of proposed approaches using permutation analysis. Lastly, to test the efficacy of our method, individual contributions extracted from the group-level SCNs were examined for associations with intelligence scores and genetic data. The extracted individual contributions were stable and were significantly related to both genetic and intelligence estimates, in both typically developing individuals and participants with FXS. We anticipate that the approaches developed in this work could be used as a putative biomarker for altered connectivity in individuals with neurodevelopmental disorders.
| S-EPMC4589508 | biostudies-literature
Project description:Characterizing interactions between multiple brain regions is important for understanding brain function. Functional connectivity measures based on partial correlation provide an estimate of the linear conditional dependence between brain regions after removing the linear influence of other regions. Estimation of partial correlations is, however, difficult when the number of regions is large, as is now increasingly the case with a growing number of large-scale brain connectivity studies. To address this problem, we develop novel methods for estimating sparse partial correlations between multiple regions in fMRI data using elastic net penalty (SPC-EN), which combines L1- and L2-norm regularization We show that L1-norm regularization in SPC-EN provides sparse interpretable solutions while L2-norm regularization improves the sensitivity of the method when the number of possible connections between regions is larger than the number of time points, and when pair-wise correlations between brain regions are high. An issue with regularization-based methods is choosing the regularization parameters which in turn determine the selection of connections between brain regions. To address this problem, we deploy novel stability selection methods to infer significant connections between brain regions. We also compare the performance of SPC-EN with existing methods which use only L1-norm regularization (SPC-L1) on simulated and experimental datasets. Detailed simulations show that the performance of SPC-EN, measured in terms of sensitivity and accuracy is superior to SPC-L1, especially at higher rates of feature prevalence. Application of our methods to resting-state fMRI data obtained from 22 healthy adults shows that SPC-EN reveals a modular architecture characterized by strong inter-hemispheric links, distinct ventral and dorsal stream pathways, and a major hub in the posterior medial cortex - features that were missed by conventional methods. Taken together, our findings suggest that SPC-EN provides a powerful tool for characterizing connectivity involving a large number of correlated regions that span the entire brain.
| S-EPMC3288428 | biostudies-literature
Project description:The emergence of high-throughput genomic datasets from different sources and platforms (e.g., gene expression, single nucleotide polymorphisms (SNP), and copy number variation (CNV)) has greatly enhanced our understandings of the interplay of these genomic factors as well as their influences on the complex diseases. It is challenging to explore the relationship between these different types of genomic data sets. In this paper, we focus on a multivariate statistical method, canonical correlation analysis (CCA) method for this problem. Conventional CCA method does not work effectively if the number of data samples is significantly less than that of biomarkers, which is a typical case for genomic data (e.g., SNPs). Sparse CCA (sCCA) methods were introduced to overcome such difficulty, mostly using penalizations with l-1 norm (CCA-l1) or the combination of l-1and l-2 norm (CCA-elastic net). However, they overlook the structural or group effect within genomic data in the analysis, which often exist and are important (e.g., SNPs spanning a gene interact and work together as a group).We propose a new group sparse CCA method (CCA-sparse group) along with an effective numerical algorithm to study the mutual relationship between two different types of genomic data (i.e., SNP and gene expression). We then extend the model to a more general formulation that can include the existing sCCA models. We apply the model to feature/variable selection from two data sets and compare our group sparse CCA method with existing sCCA methods on both simulation and two real datasets (human gliomas data and NCI60 data). We use a graphical representation of the samples with a pair of canonical variates to demonstrate the discriminating characteristic of the selected features. Pathway analysis is further performed for biological interpretation of those features.The CCA-sparse group method incorporates group effects of features into the correlation analysis while performs individual feature selection simultaneously. It outperforms the two sCCA methods (CCA-l1 and CCA-group) by identifying the correlated features with more true positives while controlling total discordance at a lower level on the simulated data, even if the group effect does not exist or there are irrelevant features grouped with true correlated features. Compared with our proposed CCA-group sparse models, CCA-l1 tends to select less true correlated features while CCA-group inclines to select more redundant features.
| S-EPMC3751310 | biostudies-literature
Project description:Intrinsic functional connectivity analysis using resting-state functional magnetic resonance imaging (rsfMRI) has become a powerful tool for examining brain functional organization. Global artifacts such as physiological noise pose a significant problem in estimation of intrinsic functional connectivity. Here we develop and test a novel random subspace method for functional connectivity (RSMFC) that effectively removes global artifacts in rsfMRI data. RSMFC estimates the partial correlation between a seed region and each target brain voxel using multiple subsets of voxels sampled randomly across the whole brain. We evaluated RSMFC on both simulated and experimental rsfMRI data and compared its performance with standard methods that rely on global mean regression (GSReg) which are widely used to remove global artifacts. Using extensive simulations we demonstrate that RSMFC is effective in removing global artifacts in rsfMRI data. Critically, using a novel simulated dataset we demonstrate that, unlike GSReg, RSMFC does not artificially introduce anti-correlations between inherently uncorrelated networks, a result of paramount importance for reliably estimating functional connectivity. Furthermore, we show that the overall sensitivity, specificity and accuracy of RSMFC are superior to GSReg. Analysis of posterior cingulate cortex connectivity in experimental rsfMRI data from 22 healthy adults revealed strong functional connectivity in the default mode network, including more reliable identification of connectivity with left and right medial temporal lobe regions that were missed by GSReg. Notably, compared to GSReg, negative correlations with lateral fronto-parietal regions were significantly weaker in RSMFC. Our results suggest that RSMFC is an effective method for minimizing the effects of global artifacts and artificial negative correlations, while accurately recovering intrinsic functional brain networks.
| S-EPMC3759623 | biostudies-literature
Project description:Imaging genetics combines brain imaging and genetic information to identify the relationships between genetic variants and brain activities. When the data samples belong to different classes (e.g. disease status), the relationships may exhibit class-specific patterns that can be used to facilitate the understanding of a disease. Conventional approaches often perform separate analysis on each class and report the differences, but ignore important shared patterns.In this paper, we develop a multivariate method to analyze the differential dependency across multiple classes. We propose a joint sparse canonical correlation analysis method, which uses a generalized fused lasso penalty to jointly estimate multiple pairs of canonical vectors with both shared and class-specific patterns. Using a data fusion approach, the method is able to detect differentially correlated modules effectively and efficiently. The results from simulation studies demonstrate its higher accuracy in discovering both common and differential canonical correlations compared to conventional sparse CCA. Using a schizophrenia dataset with 92 cases and 116 controls including a single nucleotide polymorphism (SNP) array and functional magnetic resonance imaging data, the proposed method reveals a set of distinct SNP-voxel interaction modules for the schizophrenia patients, which are verified to be both statistically and biologically significant.The Matlab code is available at https://sites.google.com/site/jianfang86/JSCCA CONTACT: wyp@tulane.eduSupplementary information: Supplementary data are available at Bioinformatics online.
| S-EPMC5181564 | biostudies-literature
Project description:BackgroundAlternative pre-mRNA splicing is an important gene regulation mechanism for expanding proteomic diversity in higher eukaryotes. Each splicing regulator can potentially influence a large group of alternative exons. Meanwhile, each alternative exon is controlled by multiple splicing regulators. The rapid accumulation of high-throughput data provides us with a unique opportunity to study the complicated alternative splicing regulatory network.ResultsWe propose the use of partial correlation analysis to identify association links between exons and their upstream regulators or their downstream target genes (exon-gene links) and links between co-spliced exons (exon-exon links). The partial correlation analysis avoids taking the ratio of two noisy random variables, exon expression level and gene expression level, so that it achieves a higher statistical power. We named this analysis procedure pCastNet (partial Correlation analysis of splicing transcriptome Network). Through studies of known alternative exons, conservation patterns, relative positions, functional annotations, and RT-PCR experiments, we concluded that pCastNet can effectively identify exon-gene or exon-exon links. We further found that gene pairs with exon-gene or exon-exon links tend to have similar functions or are present in the same pathways. More interestingly, gene pairs with exon-gene or exon-exon links tend to share cis-elements in promoter regions and microRNA binding elements in 3' untranslated regions, which suggests the coupling of co-alternative-splicing, co-transcription-factor-binding, and co-microRNA-binding.ConclusionsAlternative splicing regulatory networks reconstructed by pCastNet can help us better understand the coordinate and combinatorial nature of alternative splicing regulation. The proposed tool can be readily applied to other high-throughput data such as transcriptome sequencing data.
| S-EPMC2687791 | biostudies-literature
Project description:Both genetic variants and brain region abnormalities are recognized as important factors for complex diseases (e.g., schizophrenia). In this paper, we investigated the correspondence between single nucleotide polymorphism (SNP) and brain activity measured by functional magnetic resonance imaging (fMRI) to understand how genetic variation influences the brain activity. A group sparse canonical correlation analysis method (group sparse CCA) was developed to explore the correlation between these two datasets which are high dimensional-the number of SNPs/voxels is far greater than the number of samples. Different from the existing sparse CCA methods (sCCA), our approach can exploit structural information in the correlation analysis by introducing group constraints. A simulation study demonstrates that it outperforms the existing sCCA. We applied this method to the real data analysis and identified two pairs of significant canonical variates with average correlations of 0.4527 and 0.4292 respectively, which were used to identify genes and voxels associated with schizophrenia. The selected genes are mostly from 5 schizophrenia (SZ)-related signalling pathways. The brain mappings of the selected voxles also indicate the abnormal brain regions susceptible to schizophrenia. A gene and brain region of interest (ROI) correlation analysis was further performed to confirm the significant correlations between genes and ROIs.
| S-EPMC4007390 | biostudies-literature