Efficient and accurate causal inference with hidden confounders from genome-transcriptome variation data
Open Access
- 18 August 2017
- journal article
- research article
- Published by Public Library of Science (PLoS) in PLOS Computational Biology
- Vol. 13 (8) , e1005703
Abstract
Mapping gene expression as a quantitative trait using whole genome-sequencing and transcriptome analysis allows to discover the functional consequences of genetic variation. We developed a novel method and ultra-fast software Findr for higly accurate causal inference between gene expression traits using cis-regulatory DNA variations as causal anchors, which improves current methods by taking into consideration hidden confounders and weak regulations. Findr outperformed existing methods on the DREAM5 Systems Genetics challenge and on the prediction of microRNA and transcription factor targets in human lymphoblastoid cells, while being nearly a million times faster. Findr is publicly available at https://github.com/lingfeiwang/findr. Understanding how genetic variation between individuals determines variation in observable traits or disease risk is one of the core aims of genetics. It is known that genetic variation often affects gene regulatory DNA elements and directly causes variation in expression of nearby genes. This effect in turn cascades down to other genes via the complex pathways and gene interaction networks that ultimately govern how cells operate in an ever changing environment. In theory, when genetic variation and gene expression levels are measured simultaneously in a large number of individuals, the causal effects of genes on each other can be inferred using statistical models similar to those used in randomized controlled trials. We developed a novel method and ultra-fast software Findr which, unlike existing methods, takes into account the complex but unknown network context when predicting causality between specific gene pairs. Findr’s predictions have a significantly higher overlap with known gene networks compared to existing methods, using both simulated and real data. Findr is also nearly a million times faster, and hence the only software in its class that can handle modern datasets where the expression levels of ten-thousands of genes are simultaneously measured in hundreds to thousands of individuals.Keywords
All Related Versions
Funding Information
- Biotechnology and Biological Sciences Research Council (GB) (BB/J004235/1)
- Biotechnology and Biological Sciences Research Council (GB) (BB/M020053/1)
This publication has 39 references indexed in Scilit:
- Architecture of the human regulatory network derived from ENCODE dataNature, 2012
- Matrix eQTL: ultra fast eQTL analysis via large matrix operationsBioinformatics, 2012
- Simulating systems genetics data with SysGenSIMBioinformatics, 2011
- Characterizing the role of miRNAs within gene regulatory networks using integrative genomics techniquesMolecular Systems Biology, 2011
- Critical reasoning on causal inference in genome-wide linkage and association studiesTrends in Genetics, 2010
- Illustrating bias due to conditioning on a colliderInternational Journal of Epidemiology, 2009
- Backup in gene regulatory networks explains differences between binding and knockout resultsMolecular Systems Biology, 2009
- miRecords: an integrated resource for microRNA-target interactionsNucleic Acids Research, 2008
- An integrative genomics approach to infer causal associations between gene expression and diseaseNature Genetics, 2005
- Statistical significance for genomewide studiesProceedings of the National Academy of Sciences, 2003