Zhenqin Wu1,2, Nilah M Ioannidis2, James Zou2,3. 1. Department of Chemistry, Stanford University, CA 94305, USA. 2. Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, 94305 CA, USA. 3. Chan-Zuckerberg Biohub, San Francisco, 94158 CA, USA.
Abstract
SUMMARY: Interpreting genetic variants of unknown significance (VUS) is essential in clinical applications of genome sequencing for diagnosis and personalized care. Non-coding variants remain particularly difficult to interpret, despite making up a large majority of trait associations identified in genome-wide association studies (GWAS) analyses. Predicting the regulatory effects of non-coding variants on candidate genes is a key step in evaluating their clinical significance. Here, we develop a machine-learning algorithm, Inference of Connected expression quantitative trait loci (eQTLs) (IRT), to predict the regulatory targets of non-coding variants identified in studies of eQTLs. We assemble datasets using eQTL results from the Genotype-Tissue Expression (GTEx) project and learn to separate positive and negative pairs based on annotations characterizing the variant, gene and the intermediate sequence. IRT achieves an area under the receiver operating characteristic curve (ROC-AUC) of 0.799 using random cross-validation, and 0.700 for a more stringent position-based cross-validation. Further evaluation on rare variants and experimentally validated regulatory variants shows a significant enrichment in IRT identifying the true target genes versus negative controls. In gene-ranking experiments, IRT achieves a top-1 accuracy of 50% and top-3 accuracy of 90%. Salient features, including GC-content, histone modifications and Hi-C interactions are further analyzed and visualized to illustrate their influences on predictions. IRT can be applied to any VUS of interest and each candidate nearby gene to output a score reflecting the likelihood of regulatory effect on the expression level. These scores can be used to prioritize variants and genes to assist in patient diagnosis and GWAS follow-up studies. AVAILABILITY AND IMPLEMENTATION: Codes and data used in this work are available at https://github.com/miaecle/eQTL_Trees. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
SUMMARY: Interpreting genetic variants of unknown significance (VUS) is essential in clinical applications of genome sequencing for diagnosis and personalized care. Non-coding variants remain particularly difficult to interpret, despite making up a large majority of trait associations identified in genome-wide association studies (GWAS) analyses. Predicting the regulatory effects of non-coding variants on candidate genes is a key step in evaluating their clinical significance. Here, we develop a machine-learning algorithm, Inference of Connected expression quantitative trait loci (eQTLs) (IRT), to predict the regulatory targets of non-coding variants identified in studies of eQTLs. We assemble datasets using eQTL results from the Genotype-Tissue Expression (GTEx) project and learn to separate positive and negative pairs based on annotations characterizing the variant, gene and the intermediate sequence. IRT achieves an area under the receiver operating characteristic curve (ROC-AUC) of 0.799 using random cross-validation, and 0.700 for a more stringent position-based cross-validation. Further evaluation on rare variants and experimentally validated regulatory variants shows a significant enrichment in IRT identifying the true target genes versus negative controls. In gene-ranking experiments, IRT achieves a top-1 accuracy of 50% and top-3 accuracy of 90%. Salient features, including GC-content, histone modifications and Hi-C interactions are further analyzed and visualized to illustrate their influences on predictions. IRT can be applied to any VUS of interest and each candidate nearby gene to output a score reflecting the likelihood of regulatory effect on the expression level. These scores can be used to prioritize variants and genes to assist in patient diagnosis and GWAS follow-up studies. AVAILABILITY AND IMPLEMENTATION: Codes and data used in this work are available at https://github.com/miaecle/eQTL_Trees. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Authors: Maryam M Asgari; Wei Wang; Nilah M Ioannidis; Jacqueline Itnyre; Thomas Hoffmann; Eric Jorgenson; Alice S Whittemore Journal: J Invest Dermatol Date: 2016-01-29 Impact factor: 8.551
Authors: Nilah M Ioannidis; Joe R Davis; Marianne K DeGorter; Nicholas B Larson; Shannon K McDonnell; Amy J French; Alexis J Battle; Trevor J Hastie; Stephen N Thibodeau; Stephen B Montgomery; Carlos D Bustamante; Weiva Sieh; Alice S Whittemore Journal: Bioinformatics Date: 2017-12-15 Impact factor: 6.937
Authors: Xin Li; Alexis Battle; Konrad J Karczewski; Zach Zappala; David A Knowles; Kevin S Smith; Kim R Kukurba; Eric Wu; Noah Simon; Stephen B Montgomery Journal: Am J Hum Genet Date: 2014-09-04 Impact factor: 11.025
Authors: Robert Lesurf; Kelsy C Cotto; Grace Wang; Malachi Griffith; Katayoon Kasaian; Steven J M Jones; Stephen B Montgomery; Obi L Griffith Journal: Nucleic Acids Res Date: 2015-11-17 Impact factor: 16.971
Authors: Harvind S Chahal; Yuan Lin; Katherine J Ransohoff; David A Hinds; Wenting Wu; Hong-Ji Dai; Abrar A Qureshi; Wen-Qing Li; Peter Kraft; Jean Y Tang; Jiali Han; Kavita Y Sarin Journal: Nat Commun Date: 2016-07-18 Impact factor: 14.919
Authors: Susanna Pagni; James D Mills; Adam Frankish; Jonathan M Mudge; Sanjay M Sisodiya Journal: Neuropathol Appl Neurobiol Date: 2021-12-16 Impact factor: 6.250