HIPred: an integrative approach to predicting haploinsufficient genes.
Ontology highlight
ABSTRACT: Motivation
A major cause of autosomal dominant disease is haploinsufficiency, whereby a single copy of a gene is not sufficient to maintain the normal function of the gene. A large proportion of existing methods for predicting haploinsufficiency incorporate biological networks, e.g. protein-protein interaction networks that have recently been shown to introduce study bias. As a result, these methods tend to perform best on well-studied genes, but underperform on less studied genes. The advent of large genome sequencing consortia, such as the 1000 genomes project, NHLBI Exome Sequencing Project and the Exome Aggregation Consortium creates an urgent need for unbiased haploinsufficiency prediction methods.Results
Here, we describe a machine learning approach, called HIPred, that integrates genomic and evolutionary information from ENSEMBL, with functional annotations from the Encyclopaedia of DNA Elements consortium and the NIH Roadmap Epigenomics Project to predict haploinsufficiency, without the study bias described earlier. We benchmark HIPred using several datasets and show that our unbiased method performs as well as, and in most cases, outperforms existing biased algorithms.Availability and implementation
HIPred scores for all gene identifiers are available at: https://github.com/HAShihab/HIPred .Contact
h.shihab@bristol.ac.uk or tom.gaunt@bristol.ac.uk.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Shihab HA
PROVIDER: S-EPMC5581952 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.
| S-EPMC8354477 | biostudies-literature
Project description:Variation in gene expression has been observed in natural populations and associated with complex traits or phenotypes such as disease susceptibility and drug response. Gene expression itself is controlled by various genetic and non-genetic factors. The binding of a class of small RNA molecules, microRNAs (miRNAs), to mRNA transcript targets has recently been demonstrated to be an important mechanism of gene regulation. Because individual miRNAs may regulate the expression of multiple gene targets, a comprehensive and reliable catalogue of miRNA-regulated targets is critical to understanding gene regulatory networks. Though experimental approaches have been used to identify many miRNA targets, due to cost and efficiency, current miRNA target identification still relies largely on computational algorithms that aim to take advantage of different biochemical/thermodynamic properties of the sequences of miRNAs and their gene targets. A novel approach, ExprTarget, therefore, is proposed here to integrate some of the most frequently invoked methods (miRanda, PicTar, TargetScan) as well as the genome-wide HapMap miRNA and mRNA expression datasets generated in our laboratory. To our knowledge, this dataset constitutes the first miRNA expression profiling in the HapMap lymphoblastoid cell lines. We conducted diagnostic tests of the existing computational solutions using the experimentally supported targets in TarBase as gold standard. To gain insight into the biases that arise from such an analysis, we investigated the effect of the choice of gold standard on the evaluation of the various computational tools. We analyzed the performance of ExprTarget using both ROC curve analysis and cross-validation. We show that ExprTarget greatly improves miRNA target prediction relative to the individual prediction algorithms in terms of sensitivity and specificity. We also developed an online database, ExprTargetDB, of human miRNA targets predicted by our approach that integrates gene expression profiling into a broader framework involving important features of miRNA target site predictions.
| S-EPMC2958831 | biostudies-literature
Project description:BACKGROUND: Elucidation of the direct/indirect protein interactions and gene associations is required to fully understand the workings of the cell. This can be achieved through the use of both low- and high-throughput biological experiments and in silico methods. We present GAP (Gene functional Association Predictor), an integrative method for predicting and characterizing gene functional associations. GAP integrates different biological features using a novel taxonomy-based semantic similarity measure in predicting and prioritizing high-quality putative gene associations. The proposed similarity measure increases information gain from the available gene annotations. The annotation information is incorporated from several public pathway databases, Gene Ontology annotations as well as drug and disease associations from the scientific literature. RESULTS: We evaluated GAP by comparing its prediction performance with several other well-known functional interaction prediction tools over a comprehensive dataset of known direct and indirect interactions, and observed significantly better prediction performance. We also selected a small set of GAP's highly-scored novel predicted pairs (i.e., currently not found in any known database or dataset), and by manually searching the literature for experimental evidence accessible in the public domain, we confirmed different categories of predicted functional associations with available evidence of interaction. We also provided extra supporting evidence for subset of the predicted functionally-associated pairs using an expert curated database of genes associated to autism spectrum disorders. CONCLUSIONS: GAP's predicted "functional interactome" contains ?1M highly-scored predicted functional associations out of which about 90% are novel (i.e., not experimentally validated). GAP's novel predictions connect disconnected components and singletons to the main connected component of the known interactome. It can, therefore, be a valuable resource for biologists by providing corroborating evidence for and facilitating the prioritization of potential direct or indirect interactions for experimental validation. GAP is freely accessible through a web portal: http://ophid.utoronto.ca/gap.
| S-EPMC3663825 | biostudies-other
Project description:The discovery of novel disease-associated variations in genes is often a daunting task in highly heterogeneous disease classes. We seek a generalizable algorithm that integrates multiple publicly available genomic data sources in a machine-learning model for the prioritization of candidates identified in patients with retinal disease. To approach this problem, we generate a set of feature vectors from publicly available microarray, RNA-seq, and ChIP-seq datasets of biological relevance to retinal disease, to observe patterns in gene expression specificity among tissues of the body and the eye, in addition to photoreceptor-specific signals by the CRX transcription factor. Using these features, we describe a novel algorithm, positive and unlabeled learning for prioritization (PULP). This article compares several popular supervised learning techniques as the regression function for PULP. The results demonstrate a highly significant enrichment for previously characterized disease genes using a logistic regression method. Finally, a comparison of PULP with the popular gene prioritization tool ENDEAVOUR shows superior prioritization of retinal disease genes from previous studies. The java source code, compiled binary, assembled feature vectors, and instructions are available online at https://github.com/ahwagner/PULP.
| S-EPMC4509594 | biostudies-literature
Project description:BackgroundOne of the main current challenges in computational biology is to make sense of the huge amounts of multidimensional experimental data that are being produced. For instance, large cohorts of patients are often screened using different high-throughput technologies, effectively producing multiple patient-specific molecular profiles for hundreds or thousands of patients.ResultsWe propose and implement a network-based method that integrates such patient omics data into Patient Similarity Networks. Topological features derived from these networks were then used to predict relevant clinical features. As part of the 2017 CAMDA challenge, we have successfully applied this strategy to a neuroblastoma dataset, consisting of genomic and transcriptomic data. In particular, we observe that models built on our network-based approach perform at least as well as state of the art models. We furthermore explore the effectiveness of various topological features and observe, for instance, that redundant centrality metrics can be combined to build more powerful models.ConclusionWe demonstrate that the networks inferred from omics data contain clinically relevant information and that patient clinical outcomes can be predicted using only network topological data.ReviewersThis article was reviewed by Yang-Yu Liu, Tomislav Smuc and Isabel Nepomuceno.
| S-EPMC5992838 | biostudies-literature
Project description:Technological advances have enabled the identification of an increasingly large spectrum of single nucleotide variants within the human genome, many of which may be associated with monogenic disease or complex traits. Here, we propose an integrative approach, named FATHMM-MKL, to predict the functional consequences of both coding and non-coding sequence variants. Our method utilizes various genomic annotations, which have recently become available, and learns to weight the significance of each component annotation source.We show that our method outperforms current state-of-the-art algorithms, CADD and GWAVA, when predicting the functional consequences of non-coding variants. In addition, FATHMM-MKL is comparable to the best of these algorithms when predicting the impact of coding variants. The method includes a confidence measure to rank order predictions.
| S-EPMC4426838 | biostudies-other
Project description:Rationale: GWAS (Genome-Wide Association Studies) have identified hundreds of genetic loci associated with atrial fibrillation (AF). However, these loci explain only a small proportion of AF heritability. Objective: To develop an approach to identify additional AF-related genes by integrating multiple omics data. Methods and Results: Three types of omics data were integrated: (1) summary statistics from the AFGen 2017 GWAS; (2) a whole blood EWAS (Epigenome-Wide Association Study) of AF; and (3) a whole blood TWAS (Transcriptome-Wide Association Study) of AF. The variant-level GWAS results were collapsed into gene-level associations using fast set-based association analysis. The CpG-level EWAS results were also collapsed into gene-level associations by an adapted SNP-set Kernel Association Test approach. Both GWAS and EWAS gene-based associations were then meta-analyzed with TWAS using a fixed-effects model weighted by the sample size of each data set. A tissue-specific network was subsequently constructed using the NetWAS (Network-Wide Association Study). The identified genes were then compared with the AFGen 2018 GWAS that contained more than triple the number of AF cases compared with AFGen 2017 GWAS. We observed that the multiomics approach identified many more relevant AF-related genes than using AFGen 2018 GWAS alone (1931 versus 206 genes). Many of these genes are involved in the development and regulation of heart- and muscle-related biological processes. Moreover, the gene set identified by multiomics approach explained much more AF variance than those identified by GWAS alone (10.4% versus 3.5%). Conclusions: We developed a strategy to integrate multiple omics data to identify AF-related genes. Our integrative approach may be useful to improve the power of traditional GWAS, which might be particularly useful for rare traits and diseases with limited sample size.
| S-EPMC7004281 | biostudies-literature
Project description:We propose a statistical boosting method, termed I-Boost, to integrate multiple types of high-dimensional genomics data with clinical data for predicting survival time. I-Boost provides substantially higher prediction accuracy than existing methods. By applying I-Boost to The Cancer Genome Atlas, we show that the integration of multiple genomics platforms with clinical variables improves the prediction of survival time over the use of clinical variables alone; gene expression values are typically more prognostic of survival time than other genomics data types; and gene modules/signatures are at least as prognostic as the collection of individual gene expression data.
| S-EPMC6404283 | biostudies-literature
Project description:The current multistep carcinogenesis models of colon cancer do not fully capture the genetic heterogeneity of the disease, which is additionally complicated by the presence of passenger and driver genetic alterations. The aim of the present study was to select in the context of this significant heterogeneity additional genes functionally related to colon cancer development. Methods: High-throughput copy number and gene expression data of 36 microsatellite stable sporadic colon cancers resected from patients of a single Institution characterized for mutations in APC, KRAS, TP53 and loss of 18q were analyzed. Genes whose expression correlated with the underlying copy number pattern were selected and their association with the above listed mutations and overall survival was evaluated. Results: Gain of 20q was strongly associated with TP53 mutation, and overall survival with alterations on 7p, 8p, 13q, 18q and 20q. An association with 18q loss and gain of 8q24 was also observed. New candidate genes with a potential role in colon cancer are PLCG1 on 20q, DBC1 on 8q21 and NDGR1 on 8p24. In addition an unexpected pattern of loss and mutability was found in the region upstream of KRAS gene. Conclusions: By integrating copy number alterations with gene expression and mutations in colon cancer associated genes we have developed a strategy that identifies previously known molecular features and additional players in the molecular landscape of colon cancer.
2009-08-27 | GSE16125 | GEO
Project description:The current multistep carcinogenesis models of colon cancer do not fully capture the genetic heterogeneity of the disease, which is additionally complicated by the presence of passenger and driver genetic alterations. The aim of the present study was to select in the context of this significant heterogeneity additional genes functionally related to colon cancer development. Methods: High-throughput copy number and gene expression data of 36 microsatellite stable sporadic colon cancers resected from patients of a single Institution characterized for mutations in APC, KRAS, TP53 and loss of 18q were analyzed. Genes whose expression correlated with the underlying copy number pattern were selected and their association with the above listed mutations and overall survival was evaluated. Results: Gain of 20q was strongly associated with TP53 mutation, and overall survival with alterations on 7p, 8p, 13q, 18q and 20q. An association with 18q loss and gain of 8q24 was also observed. New candidate genes with a potential role in colon cancer are PLCG1 on 20q, DBC1 on 8q21 and NDGR1 on 8p24. In addition an unexpected pattern of loss and mutability was found in the region upstream of KRAS gene. Conclusions: By integrating copy number alterations with gene expression and mutations in colon cancer associated genes we have developed a strategy that identifies previously known molecular features and additional players in the molecular landscape of colon cancer. A total of 48 sporadic colon cancer samples were analyzed by Affymetrix Mapping 250K Nsp SNP Arrays and 36 of them were also analyzed by Affymetrix Human Exon 1.0 ST Array [transcript (gene) version]. Short summary: Expression data was correlated to copy number data to identify genes whose expression was induced by copy number changes. Gene dosage candidates were then evaluated for their association with gene mutation status of APC, KRAS and TP53, loss of heterozigosity of 18q and overall survival. Long summary: Raw intensity .CEL files of the SNP arrays were processed with Chromosome Copy Number Analysis Tool (v.1.5.6 Affymetrix, Santa Clara, CA) to identify chromosomal gains and losses. Forty eight normal samples from the HapMap project supplied by Affymetrix were used as an un-paired reference set [http://www.affymetrix.com/support/technical/sample_data/500k_data.affx]. All genomic coordinates of the SNP array probes were mapped to the Human Mar. 2006 assembly the UCSC genome browser. Raw intensity .CEL files of the exon arrays were processed with the Robust Multi-Array implementation of Affymetrix Power Tools (v.1.8.6) using the core set of features (22011 probesets). All plots and analysis steps following this processing were done using the R programming language version 2.9.0Bioconductor packages. The identification of statistically significant segments of aberration in the copy number data was performed using the default parameters of the KC-SMART algorithm and 1000 permutations (KCsmart v.2.2.0). The identification and annotation of genes within each aberrant segment was performed using biomaRt v.2.0.0. The list of genes was further annotated by using the cancerGenes resource, Cancer Gene Census and the list of breast and colon CAN genes listed in Wood et al. The genomic landscapes of human breast and colorectal cancers. Science 2007;318:1108-1113. Gene dosage effects across 36 samples for which gene expression and copy number data was available was assessed by evaluating the Spearman correlation of the raw continuous copy number (log-ratios) expression (log-intensity) for each genomic region surrounding each gene falling within the segments identified by KC-SMART. Prior to performing this task, to reduce the total number of correlation tests to perform, we filtered the gene expression dataset by removing all entries having no Gene Symbol annotation 17291 probesets and half the remaining dataset exhibiting the lowest variance resulting in 8645 probesets. Category package version 2.10.0 was used to apply a linear model-based test to detect enrichment of systematic high correlation in specific chromosomal bands taking into consideration the hierarchical structure of the bands. Gene expression class comparisons and survival analysis of the selected gene dosage candidates were performed using two sample t-tests and Cox proportional hazards regression. P-value adjustment of the correlation tests and the gene expression associations was performed using the step-up false discovery rate (FDR) controlling procedure of Benjamini and Hochberg. Raw data, supplementary methods and figures supplied as reproducible documentation (Sweave file) available from http://pierotti.group.ifom-ieo-campus.it/suppl/crc.html.
2010-04-10 | E-GEOD-16125 | biostudies-arrayexpress