Dataset Information

Canonical correlation analysis for RNA-seq co-expression networks.

ABSTRACT: Digital transcriptome analysis by next-generation sequencing discovers substantial mRNA variants. Variation in gene expression underlies many biological processes and holds a key to unravelling mechanism of common diseases. However, the current methods for construction of co-expression networks using overall gene expression are originally designed for microarray expression data, and they overlook a large number of variations in gene expressions. To use information on exon, genomic positional level and allele-specific expressions, we develop novel component-based methods, single and bivariate canonical correlation analysis, for construction of co-expression networks with RNA-seq data. To evaluate the performance of our methods for co-expression network inference with RNA-seq data, they are applied to lung squamous cell cancer expression data from TCGA database and our bipolar disorder and schizophrenia RNA-seq study. The preliminary results demonstrate that the co-expression networks constructed by canonical correlation analysis and RNA-seq data provide rich genetic and molecular information to gain insight into biological processes and disease mechanism. Our new methods substantially outperform the current statistical methods for co-expression network construction with microarray expression data or RNA-seq data based on overall gene expression levels.

SUBMITTER: Hong S

PROVIDER: S-EPMC3632131 | biostudies-literature | 2013 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Canonical correlation analysis for RNA-seq co-expression networks.

Hong Shengjun S Chen Xiangning X Jin Li L Xiong Momiao M

Nucleic acids research 20130304 8

Digital transcriptome analysis by next-generation sequencing discovers substantial mRNA variants. Variation in gene expression underlies many biological processes and holds a key to unravelling mechanism of common diseases. However, the current methods for construction of co-expression networks using overall gene expression are originally designed for microarray expression data, and they overlook a large number of variations in gene expressions. To use information on exon, genomic positional lev ...[more]

PMID: 23460206

Similar Datasets

Project description:Simultaneous mapping of multiple behavioral domains into brain networks remains a major challenge. Here, we shed some light on this problem by employing a combination of machine learning, structural and functional brain networks at different spatial resolutions (also known as scales), together with performance scores across multiple neurobehavioral domains, including sensation, motor skills, and cognition. Provided by the Human Connectome Project, we make use of three cohorts: 640 participants for model training, 160 subjects for validation, and 200 subjects for model performance testing thus enhancing prediction generalization. Our modeling consists of two main stages, namely dimensionality reduction in brain network features at multiple scales, followed by canonical correlation analysis, which determines an optimal linear combination of connectivity features to predict multiple behavioral performance scores. To assess the differences in the predictive power of each modality, we separately applied three different strategies: structural unimodal, functional unimodal, and multimodal, that is, structural in combination with functional features of the brain network. Our results show that the multimodal association outperforms any of the unimodal analyses. Then, to answer which human brain structures were most involved in predicting multiple behavioral scores, we simulated different synthetic scenarios in which in each case we completely deleted a brain structure or a complete resting state network, and recalculated performance in its absence. In deletions, we found critical structures to affect performance when predicting single behavioral domains, but this occurred in a lesser manner for prediction of multi-domain behavior. Overall, our results confirm that although there are synergistic contributions between brain structure and function that enhance behavioral prediction, brain networks may also be mutually redundant in predicting multidomain behavior, such that even after deletion of a structure, the connectivity of the others can compensate for its lack in predicting behavior.

Project description:BackgroundGene co-expression networks (GCNs) describe relationships among expressed genes key to maintaining cellular identity and homeostasis. However, the small sample size of typical RNA-seq experiments which is several orders of magnitude fewer than the number of genes is too low to infer GCNs reliably. recount3, a publicly available dataset comprised of 316,443 uniformly processed human RNA-seq samples, provides an opportunity to improve power for accurate network reconstruction and obtain biological insight from the resulting networks.ResultsWe compared alternate aggregation strategies to identify an optimal workflow for GCN inference by data aggregation and inferred three consensus networks: a universal network, a non-cancer network, and a cancer network in addition to 27 tissue context-specific networks. Central network genes from our consensus networks were enriched for evolutionarily constrained genes and ubiquitous biological pathways, whereas central context-specific network genes included tissue-specific transcription factors and factorization based on the hubs led to clustering of related tissue contexts. We discovered that annotations corresponding to context-specific networks inferred from aggregated data were enriched for trait heritability beyond known functional genomic annotations and were significantly more enriched when we aggregated over a larger number of samples.ConclusionThis study outlines best practices for network GCN inference and evaluation by data aggregation. We recommend estimating and regressing confounders in each data set before aggregation and prioritizing large sample size studies for GCN reconstruction. Increased statistical power in inferring context-specific networks enabled the derivation of variant annotations that were enriched for concordant trait heritability independent of functional genomic annotations that are context-agnostic. While we observed strictly increasing held-out log-likelihood with data aggregation, we noted diminishing marginal improvements. Future directions aimed at alternate methods for estimating confounders and integrating orthogonal information from modalities such as Hi-C and ChIP-seq can further improve GCN inference.

Project description:BackgroundWith the rapid development of new genetic measurement methods, several types of genetic alterations can be quantified in a high-throughput manner. While the initial focus has been on investigating each data set separately, there is an increasing interest in studying the correlation structure between two or more data sets. Multivariate methods based on Canonical Correlation Analysis (CCA) have been proposed for integrating paired genetic data sets. The high dimensionality of microarray data imposes computational difficulties, which have been addressed for instance by studying the covariance structure of the data, or by reducing the number of variables prior to applying the CCA. In this work, we propose a new method for analyzing high-dimensional paired genetic data sets, which mainly emphasizes the correlation structure and still permits efficient application to very large data sets. The method is implemented by translating a regularized CCA to its dual form, where the computational complexity depends mainly on the number of samples instead of the number of variables. The optimal regularization parameters are chosen by cross-validation. We apply the regularized dual CCA, as well as a classical CCA preceded by a dimension-reducing Principal Components Analysis (PCA), to a paired data set of gene expression changes and copy number alterations in leukemia.ResultsUsing the correlation-maximizing methods, regularized dual CCA and PCA+CCA, we show that without pre-selection of known disease-relevant genes, and without using information about clinical class membership, an exploratory analysis singles out two patient groups, corresponding to well-known leukemia subtypes. Furthermore, the variables showing the highest relevance to the extracted features agree with previous biological knowledge concerning copy number alterations and gene expression changes in these subtypes. Finally, the correlation-maximizing methods are shown to yield results which are more biologically interpretable than those resulting from a covariance-maximizing method, and provide different insight compared to when each variable set is studied separately using PCA.ConclusionsWe conclude that regularized dual CCA as well as PCA+CCA are useful methods for exploratory analysis of paired genetic data sets, and can be efficiently implemented also when the number of variables is very large.

Dataset Information

Canonical correlation analysis for RNA-seq co-expression networks.

Publications

Canonical correlation analysis for RNA-seq co-expression networks.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets