Dataset Information

Integrative analysis of GWASs, human protein interaction, and gene expression identified gene modules associated with BMDs.

ABSTRACT: To date, few systems genetics studies in the bone field have been performed. We designed our study from a systems-level perspective by integrating genome-wide association studies (GWASs), human protein-protein interaction (PPI) network, and gene expression to identify gene modules contributing to osteoporosis risk.First we searched for modules significantly enriched with bone mineral density (BMD)-associated genes in human PPI network by using 2 large meta-analysis GWAS datasets through a dense module search algorithm. One included 7 individual GWAS samples (Meta7). The other was from the Genetic Factors for Osteoporosis Consortium (GEFOS2). One was assigned as a discovery dataset and the other as an evaluation dataset, and vice versa.In total, 42 modules and 129 modules were identified significantly in both Meta7 and GEFOS2 datasets for femoral neck and spine BMD, respectively. There were 3340 modules identified for hip BMD only in Meta7. As candidate modules, they were assessed for the biological relevance to BMD by gene set enrichment analysis in 2 expression profiles generated from circulating monocytes in subjects with low versus high BMD values. Interestingly, there were 2 modules significantly enriched in monocytes from the low BMD group in both gene expression datasets (nominal P value <.05). Two modules had 16 nonredundant genes. Functional enrichment analysis revealed that both modules were enriched for genes involved in Wnt receptor signaling and osteoblast differentiation.We highlighted 2 modules and novel genes playing important roles in the regulation of bone mass, providing important clues for therapeutic approaches for osteoporosis.

SUBMITTER: He H

PROVIDER: S-EPMC4223444 | biostudies-literature | 2014 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Integrative analysis of GWASs, human protein interaction, and gene expression identified gene modules associated with BMDs.

He Hao H Zhang Lei L Li Jian J Wang Yu-Ping YP Zhang Ji-Gang JG Shen Jie J Guo Yan-Fang YF Deng Hong-Wen HW

The Journal of clinical endocrinology and metabolism 20140813 11

<h4>Context</h4>To date, few systems genetics studies in the bone field have been performed. We designed our study from a systems-level perspective by integrating genome-wide association studies (GWASs), human protein-protein interaction (PPI) network, and gene expression to identify gene modules contributing to osteoporosis risk.<h4>Methods</h4>First we searched for modules significantly enriched with bone mineral density (BMD)-associated genes in human PPI network by using 2 large meta-analysi ...[more]

PMID: 25119315

Similar Datasets

Project description:BackgroundMore than 90% of neuroblastoma patients are cured in the low-risk group while only less than 50% for those with high-risk disease can be cured. Since the high-risk patients still have poor outcomes, we need more accurate stratification to establish an individualized precise treatment plan for the patients to improve the long-term survival rate.ResultsWe focus on extracting features and providing a workflow to improve survival prediction for neuroblastoma patients. With a workflow for gene co-expression network (GCN) mining in microarray and RNA-Seq datasets, we extracted molecular features from each co-expressed module and summarized them into eigengenes. Then we adopted the lasso-regularized Cox proportional hazards model to select the most informative eigengene features regarding association to the risk of metastasis. Nine eigengenes were selected which show strong association with patient survival prognosis. All of the nine corresponding gene modules also have highly enriched biological functions or cytoband locations. Three of them are unique modules to RNA-Seq data, which complement the modules from microarray data in terms of survival prognosis. We then merged all eigengenes from these unique modules and used an integrative method called Similarity Network Fusion to test the prognostic power of these eigengenes for prognosis. The prognostic accuracies are significantly improved as compared to using all eigengenes, and a subgroup of patients with very poor survival rate was identified.ConclusionsWe first compared GCNs mined from microarray and RNA-seq data. We discovered that each data modality yields unique GCNs, which are enriched with clear biological functions. Then we do module unique analysis and use lasso-cox model to select survival-associated eigengenes. Integration of unique and survival-associated eigengenes from both data types provides complementary information that leads to more accurate survival prognosis.ReviewersReviewed by Susmita Datta, Marco Chierici and Dimitar Vassilev.

Project description:Schizophrenia (SCZ) is a severe mental disorder that may result in hallucinations, delusions, and extremely disordered thinking. How each cell type in the brain contributes to SCZ occurrence is still unclear. Here, we leveraged the human dorsolateral prefrontal cortex bulk RNA-seq data, then used the RNA-seq deconvolution algorithm CIBERSORTx to generate SCZ brain single-cell RNA-seq data for a comprehensive analysis to understand SCZ-associated brain cell types and gene expression changes. Firstly, we observed that the proportions of brain cell types in SCZ differed from normal samples. Among these cell types, astrocyte, pericyte, and PAX6 cells were found to have a higher proportion in SCZ patients (astrocyte: SCZ = 0.163, control = 0.145, P.adj = 4.9 × 10-4, effect size = 0.478; pericyte: SCZ = 0.057, control = 0.066, P.adj = 1.1 × 10-4, effect size = 0.519; PAX6: SCZ = 0.014, control = 0.011, P.adj = 0.014, effect size = 0.377), while the L5/6_IT_CAR3 cells and LAMP5 cells are the exact opposite (L5/6_IT_Car3: SCZ = 0.102, control = 0.108, P.adj = 0.016, effect size = 0.369; LAMP5: SCZ = 0.057, control = 0.066, P.adj = 2.2 × 10-6, effect size = 0.617). Next, we investigated gene expression in cell types and functional pathways in SCZ. We observed chemical synaptic transmission dysregulation in two types of GABAergic neurons (PVALB and LAMP5), and immune reaction involvement in GABAergic neurons (SST) and non-neuronal cell types (endothelial and oligodendrocyte). Furthermore, we observed that some differential expression genes from bulk RNA-seq displayed cell-type-specific abnormalities in the expression of molecules in SCZ. Finally, the cell types with the SCZ-related transcriptomic changes could be considered to belong to the same module since we observed two major similar coordinated transcriptomic changes across these cell types. Together, our results offer novel insights into cellular heterogeneity and the molecular mechanisms underlying SCZ.

Project description:Gene expression profiles have been frequently integrated with the human protein interactome to uncover functional modules under specific conditions like disease state. Beyond traditional differential expression analysis, differential co-expression analysis has emerged as a robust approach to reveal condition-specific network modules, with successful applications in a few human disease studies. Hepatocellular carcinoma (HCC), which is often interrelated with the Hepatitis C virus, typically develops through multiple stages. A comprehensive investigation of HCC progression-specific differential co-expression modules may advance our understanding of HCC's pathophysiological mechanisms.Compared with differentially expressed genes, differentially co-expressed genes were found more likely enriched with Hepatitis C virus binding proteins and cancer-mutated genes, and they were clustered more densely in the human reference protein interaction network. These observations indicated that a differential co-expression approach could outperform the standard differential expression network analysis in searching for disease-related modules. We then proposed a differential co-expression network approach to uncover network modules involved in HCC development. Specifically, we discovered subnetworks that enriched differentially co-expressed gene pairs in each HCC transition stage, and further resolved modules with coherent co-expression change patterns over all HCC developmental stages. Our identified network modules were enriched with HCC-related genes and implicated in cancer-related biological functions. In particular, APC and YWHAZ were highlighted as two most remarkable genes in the network modules, and their dynamic interaction partnership was resolved in HCC development.We demonstrated that integration of differential co-expression with the protein interactome could outperform the traditional differential expression approach in discovering network modules of human diseases. In our application of this approach to HCC's gene expression data, we successfully identified subnetworks with marked differential co-expression in individual HCC stage transitions and network modules with coherent co-expression change patterns over all HCC developmental stages. Our results shed light on subtle HCC mechanisms, including temporal activation and dismissal of pivotal functions and dynamic interaction partnerships of key genes.

Dataset Information

Integrative analysis of GWASs, human protein interaction, and gene expression identified gene modules associated with BMDs.

Publications

Integrative analysis of GWASs, human protein interaction, and gene expression identified gene modules associated with BMDs.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets