A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets

Close to 50 genetic loci have been associated with type 2 diabetes (T2D), but they explain only 15% of the heritability. In an attempt to identify additional T2D genes, we analyzed global gene expression in human islets from 63 donors. Using 48 genes located near T2D risk variants, we identified gene coexpression and protein-protein interaction networks that were strongly associated with islet insulin secretion and HbA(1c). We integrated our data to form a rank list of putative T2D genes, of which CHL1, LRFN2, RASGRP1, and PPM1K were validated in INS-1 cells to influence insulin secretion, whereas GPR120 affected apoptosis in islets. Expression variation of the top 20 genes explained 24% of the variance in HbA(1c) with no claim of the direction. The data present a global map of genes associated with islet dysfunction and demonstrate the value of systems genetics for the identification of genes potentially involved in T2D.     

Detecting shared pathways linked to rheumatoid arthritis with other autoimmune diseases in a in silico analysis

Pathway-based analysis approach has exploded in use during the last several years. It is successful in recognizing additional biological insight of disease and finding groupings of risk genes that represent disease developing processes. Therefore, shared pathways, with pleiotropic effects, are important for understanding similar pathogenesis and indicating the common genetic origin of certain diseases. Here, we present a pathway analysis to reveal the potential disease associations between RA and three potential RA-related autoimmune diseases: psoriasis, diabetes mellitus, type 1 (T1D) and systemic lupus erythematosus (SLE). First, a comprehensive knowledge mining of public databases is performed to discover risk genes associated with RA, T1D, SLE and psoriasis; then by enrichment test of these genes, disease-related risk pathways are detected to recognize the pathways common for RA and three other diseases. Finally, the underlying disease associations are evaluated with the association rules mining method. In total, we identify multiple RA risk pathways with significant pleiotropic effects, the most unsurprising of which are the immunology related pathways. Meanwhile for the first time we highlight the involvement of the viral myocarditis pathway related to cardiovascular disease (CVD) in autoimmune diseases such as RA, psoriasis, T1D and SLE. Further Association rule mining results validate the strong association between RA and T1D and RA and SLE. It is clear that pleiotropy is a common property of pathways associated with disease traits. We provide novel pathway associations among RA and three autoimmune diseases. These results ascertain that there are shared genetic risk profiles that predispose individuals to autoimmune diseases.     

Community Structure Analysis of Gene Interaction Networks in Duchenne Muscular Dystrophy

Duchenne Muscular Dystrophy (DMD) is an important pathology associated with the human skeletal muscle and has been studied extensively. Gene expression measurements on skeletal muscle of patients afflicted with DMD provides the opportunity to understand the underlying mechanisms that lead to the pathology. Community structure analysis is a useful computational technique for understanding and modeling genetic interaction networks. In this paper, we leverage this technique in combination with gene expression measurements from normal and DMD patient skeletal muscle tissue to study the structure of genetic interactions in the context of DMD. We define a novel framework for transforming a raw dataset of gene expression measurements into an interaction network, and subsequently apply algorithms for community structure analysis for the extraction of topological communities. The emergent communities are analyzed from a biological standpoint in terms of their constituent biological pathways, and an interpretation that draws correlations between functional and structural organization of the genetic interactions is presented. We also compare these communities and associated functions in pathology against those in normal human skeletal muscle. In particular, differential enhancements are observed in the following pathways between pathological and normal cases: Metabolic, Focal adhesion, Regulation of actin cytoskeleton and Cell adhesion, and implication of these mechanisms are supported by prior work. Furthermore, our study also includes a gene-level analysis to identify genes that are involved in the coupling between the pathways of interest. We believe that our results serve to highlight important distinguishing features in the structural/functional organization of constituent biological pathways, as it relates to normal and DMD cases, and provide the mechanistic basis for further biological investigations into specific pathways differently regulated between normal and DMD patients. These findings have the potential to serve as fertile ground for therapeutic applications involving targeted drug development for DMD.     

Integrating gene expression and metabolic profiles

Recent advances in high throughput technologies have generated an abundance of biological information, such as gene expression, protein-protein interaction, and metabolic data. These various types of data capture different aspects of the cellular response to environmental factors. Integrating data from different measurements enhances the ability of modeling frameworks to predict cellular function more accurately and can lead to a more coherent reconstruction of the underlying regulatory network structure. Different techniques, newly developed and borrowed, have been applied for the purpose of extracting this information from experimental data. In this study, we developed a framework to integrate metabolic and gene expression profiles for a hepatocellular system. Specifically, we applied genetic algorithm and partial least square analysis to identify important genes relevant to a specific cellular function. We identified genes 1) whose expression levels quantitatively predict a metabolic function and 2) that play a part in regulating a hepatocellular function and reconstructed their role in the metabolic network. The framework 1) preprocesses the gene expression data using statistical techniques, 2) selects genes using a genetic algorithm and couples them to a partial least squares analysis to predict cellular function, and 3) reconstructs, with the assistance of a literature search, the pathways that regulate cellular function, namely intracellular triglyceride and urea synthesis. This provides a framework for identifying cellular pathways that are active as a function of the environment and in turn helps to uncover the interplay between gene and metabolic networks.     

The Type 1 Diabetes - HLA Susceptibility Interactome - Identification of HLA Genotype-Specific Disease Genes for Type 1 Diabetes

Background

The individual contribution of genes in the HLA region to the risk of developing type 1 diabetes (T1D) is confounded by the high linkage disequilibrium (LD) in this region. Using a novel approach we have combined genetic association data with information on functional protein-protein interactions to elucidate risk independent of LD and to place the genetic association into a functional context.

Methodology/Principal Findings

Genetic association data from 2300 single nucleotide polymorphisms (SNPs) in the HLA region was analysed in 2200 T1D family trios divided into six risk groups based on HLA-DRB1 genotypes. The best SNP signal in each gene was mapped to proteins in a human protein interaction network and their significance of clustering in functional network modules was evaluated. The significant network modules identified through this approach differed between the six HLA risk groups, which could be divided into two groups based on carrying the DRB1*0301 or the DRB1*0401 allele. Proteins identified in networks specific for DRB1*0301 carriers were involved in stress response and inflammation whereas in DRB1*0401 carriers the proteins were involved in antigen processing and presentation.

Conclusions/Significance

In this study we were able to hypothesise functional differences between individuals with T1D carrying specific DRB1 alleles. The results point at candidate proteins involved in distinct cellular processes that could not only help the understanding of the pathogenesis of T1D, but also the distinction between individuals at different genetic risk for developing T1D.     

Integrative Genomics Reveals Novel Molecular Pathways and Gene Networks for Coronary Artery Disease

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions.     

RASSF1A对黑色素瘤A375细胞基因网络的影响

RASSF1A(Ras association domain family 1 isoform A)是定位于染色体3p21.3区域的抑瘤基因,编码一个由340个氨基酸残基构成的微管相关蛋白.该基因在包括恶性黑色素瘤在内的多种肿瘤中因启动子高甲基化而表达沉默.本研究建立了RASSF1A稳定表达的恶性黑色素瘤A375细胞系,通过全基因组表达谱基因芯片分析RASSF1A过表达对A375细胞基因表达谱的影响,发现RASSF1A引起184个基因表达上调,26个基因表达下调.通过Realtime RT-PCR对部分差异表达基因进行验证,结果表明与芯片筛选结果一致.RASSF1A影响的差异表达基因功能上归属于细胞生长与增殖、细胞周期、细胞凋亡、细胞间黏附、信号传导等生物过程.采用STRING软件构建了RASSF1A影响的差异表达基因调控网络,结果表明RASSF1A调控的差异表达基因构成一个高连接度的基因网络.其中,炎症细胞因子、转录因子位于网络中央.RASSF1A通过影响炎症细胞因子与转录因子之间的表达,影响A375细胞基因网络,调节黑色素瘤恶性生物学行为.     

Correlation between Oncogenic Mutations and Parameter Sensitivity of the Apoptosis Pathway Model

One of the major breakthroughs in oncogenesis research in recent years is the discovery that, in most patients, oncogenic mutations are concentrated in a few core biological functional pathways. This discovery indicates that oncogenic mechanisms are highly related to the dynamics of biologic regulatory networks, which govern the behaviour of functional pathways. Here, we propose that oncogenic mutations found in different biological functional pathways are closely related to parameter sensitivity of the corresponding networks. To test this hypothesis, we focus on the DNA damage-induced apoptotic pathway—the most important safeguard against oncogenesis. We first built the regulatory network that governs the apoptosis pathway, and then translated the network into dynamics equations. Using sensitivity analysis of the network parameters and comparing the results with cancer gene mutation spectra, we found that parameters that significantly affect the bifurcation point correspond to high-frequency oncogenic mutations. This result shows that the position of the bifurcation point is a better measure of the functionality of a biological network than gene expression levels of certain key proteins. It further demonstrates the suitability of applying systems-level analysis to biological networks as opposed to studying genes or proteins in isolation.     

Network motif analysis of a multi-mode genetic-interaction network

Different modes of genetic interaction indicate different functional relationships between genes. The extraction of biological information from dense multi-mode genetic-interaction networks demands appropriate statistical and computational methods. We developed such methods and implemented them in open-source software. Motifs extracted from multi-mode genetic-interaction networks form functional subnetworks, highlight genes dominating these subnetworks, and reveal genetic reflections of the underlying biochemical system.     

Identification of key candidate genes and pathways in multiple myeloma by integrated bioinformatics analysis

Multiple myeloma (MM) is a common hematologic malignancy for which the underlying molecular mechanisms remain largely unclear. This study aimed to elucidate key candidate genes and pathways in MM by integrated bioinformatics analysis. Expression profiles GSE6477 and GSE47552 were obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) with p < .05 and [logFC] > 1 were identified. Functional enrichment, protein–protein interaction network construction and survival analyses were then performed. First, 51 upregulated and 78 downregulated DEGs shared between the two GSE datasets were identified. Second, functional enrichment analysis showed that these DEGs are mainly involved in the B cell receptor signaling pathway, hematopoietic cell lineage, and NF-kappa B pathway. Moreover, interrelation analysis of immune system processes showed enrichment of the downregulated DEGs mainly in B cell differentiation, positive regulation of monocyte chemotaxis and positive regulation of T cell proliferation. Finally, the correlation between DEG expression and survival in MM was evaluated using the PrognoScan database. In conclusion, we identified key candidate genes that affect the outcomes of patients with MM, and these genes might serve as potential therapeutic targets.     

Can Genetic Pleiotropy Replicate Common Clinical Constellations of Cardiovascular Disease and Risk?

The relationship between obesity, diabetes, hyperlipidemia, hypertension, kidney disease and cardiovascular disease (CVD) is established when looked at from a clinical, epidemiological or pathophysiological perspective. Yet, when viewed from a genetic perspective, there is comparatively little data synthesis that these conditions have an underlying relationship. We sought to investigate the overlap of genetic variants independently associated with each of these commonly co-existing conditions from the NHGRI genome-wide association study (GWAS) catalog, in an attempt to replicate the established notion of shared pathophysiology and risk. We used pathway-based analyses to detect subsets of pleiotropic genes involved in similar biological processes. We identified 107 eligible GWAS studies related to CVD and its established comorbidities and risk factors and assigned genes that correspond to the associated signals based on their position. We found 44 positional genes shared across at least two CVD-related phenotypes that independently recreated the established relationship between the six phenotypes, but only if studies representing non-European populations were included. Seven genes revealed pleiotropy across three or more phenotypes, mostly related to lipid transport and metabolism. Yet, many genes had no relationship to each other or to genes with established functional connection. Whilst we successfully reproduced established relationships between CVD risk factors using GWAS findings, interpretation of biological pathways involved in the observed pleiotropy was limited. Further studies linking genetic variation to gene expression, as well as describing novel biological pathways will be needed to take full advantage of GWAS results.