Predicting genotype-specific gene regulatory networks

Understanding how each person’s unique genotype influences their individual patterns of gene regulation has the potential to improve our understanding of human health and development and to refine genotype-specific disease risk assessments and treatments. However, the effects of genetic variants are not typically considered when constructing gene regulatory networks, despite the fact that many disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding. We developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network (GRN) for each individual in a study population. EGRET begins by constructing a genotype-informed TF-gene prior network derived using TF motif predictions, eQTL data, individual genotypes, and the predicted effects of genetic variants on TF binding. It then uses message passing to integrate this prior network with gene expression and TF protein-protein interaction data to produce a refined, genotype-specific regulatory network. We used EGRET to infer GRNs for two blood-derived cell lines and identified genotype-associated, cell-line specific regulatory differences that we subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential ChIP-seq TF binding. We also inferred EGRET GRNs for three cell types from each of 119 individuals and identified cell type-specific regulatory differences associated with diseases related to those cell types. EGRET is, to our knowledge, the first method that infers networks that reflect individual genetic variation in a way that provides insight into genetic regulatory associations that drive complex phenotypes. EGRET is available through the Network Zoo R package (netZooR v0.9; netzoo.github.io ).