A Sequence-based Method to Predict the Impact of Regulatory Variants Using Random Forest

Overview

Journal BMC Syst Biol

Publisher Biomed Central

Specialty Biology

Date 2017 Apr 1

PMID 28361702

Citations 7

Authors

Qiao Liu

Mingxin Gan

Rui Jiang

Affiliations

Soon will be listed here.

Abstract

Background: Most disease-associated variants identified by genome-wide association studies (GWAS) exist in noncoding regions. In spite of the common agreement that such variants may disrupt biological functions of their hosting regulatory elements, it remains a great challenge to characterize the risk of a genetic variant within the implicated genome sequence. Therefore, it is essential to develop an effective computational model that is not only capable of predicting the potential risk of a genetic variant but also valid in interpreting how the function of the genome is affected with the occurrence of the variant.

Results: We developed a method named kmerForest that used a random forest classifier with k-mer counts to predict accessible chromatin regions purely based on DNA sequences. We demonstrated that our method outperforms existing methods in distinguishing known accessible chromatin regions from random genomic sequences. Furthermore, the performance of our method can further be improved with the incorporation of sequence conservation features. Based on this model, we assessed importance of the k-mer features by a series of permutation experiments, and we characterized the risk of a single nucleotide polymorphism (SNP) on the function of the genome using the difference between the importance of the k-mer features affected by the occurrence of the SNP. We conducted a series of experiments and showed that our model can well discriminate between pathogenic and normal SNPs. Particularly, our model correctly prioritized SNPs that are proved to be enriched for the binding sites of FOXA1 in breast cancer cell lines from previous studies.

Conclusions: We presented a novel method to interpret functional genetic variants purely base on DNA sequences. The proposed k-mer based score offers an effective means of measuring the impact of SNPs on the function of the genome, and thus shedding light on the identification of genetic risk factors underlying complex traits and diseases.

Citing Articles

EpiGePT: a pretrained transformer-based language model for context-specific human epigenomics.

Gao Z, Liu Q, Zeng W, Jiang R, Wong W Genome Biol. 2024; 25(1):310.

PMID: 39696471 PMC: 11657395. DOI: 10.1186/s13059-024-03449-7.

Reproductomics: Exploring the Applications and Advancements of Computational Tools.

SENGUPTA P, Dutta S, Liew F, Samrot A, Dasgupta S, Rajput M Physiol Res. 2024; 73(5):687-702.

PMID: 39530905 PMC: 11629954. DOI: 10.33549/physiolres.935389.

Classification of group A rotavirus VP7 and VP4 genotypes using random forest.

Tran H, Friendship R, Poljak Z Front Genet. 2023; 14:1029185.

PMID: 37323680 PMC: 10267748. DOI: 10.3389/fgene.2023.1029185.

DeepCAGE: Incorporating Transcription Factors in Genome-wide Prediction of Chromatin Accessibility.

Liu Q, Hua K, Zhang X, Wong W, Jiang R Genomics Proteomics Bioinformatics. 2022; 20(3):496-507.

PMID: 35293310 PMC: 9801045. DOI: 10.1016/j.gpb.2021.08.015.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human .

Liu H, Duncan K, Helverson A, Kumari P, Mumm C, Xiao Y Elife. 2020; 9.

PMID: 32031521 PMC: 7039683. DOI: 10.7554/eLife.51325.

References

Barenboim M, Manke T . ChroMoS: an integrated web tool for SNP classification, prioritization and functional interpretation. Bioinformatics. 2013; 29(17):2197-8. PMC: 3740627. DOI: 10.1093/bioinformatics/btt356. View

Hindorff L, Sethupathy P, Junkins H, Ramos E, Mehta J, Collins F . Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009; 106(23):9362-7. PMC: 2687147. DOI: 10.1073/pnas.0903103106. View

Eeckhoute J, Carroll J, Geistlinger T, Torres-Arzayus M, Brown M . A cell-type-specific transcriptional network required for estrogen regulation of cyclin D1 and cell cycle progression in breast cancer. Genes Dev. 2006; 20(18):2513-26. PMC: 1578675. DOI: 10.1101/gad.1446006. View

McGaughey D, Vinton R, Huynh J, Al-Saif A, Beer M, McCallion A . Metrics of sequence constraint overlook regulatory sequences in an exhaustive analysis at phox2b. Genome Res. 2007; 18(2):252-60. PMC: 2203623. DOI: 10.1101/gr.6929408. View

Hirschhorn J, Daly M . Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005; 6(2):95-108. DOI: 10.1038/nrg1521. View

Stenson P, Mort M, Ball E, Howells K, Phillips A, Thomas N . The Human Gene Mutation Database: 2008 update. Genome Med. 2009; 1(1):13. PMC: 2651586. DOI: 10.1186/gm13. View

Abecasis G, Auton A, Brooks L, DePristo M, Durbin R, Handsaker R . An integrated map of genetic variation from 1,092 human genomes. Nature. 2012; 491(7422):56-65. PMC: 3498066. DOI: 10.1038/nature11632. View

Robertson G, Hirst M, Bainbridge M, Bilenky M, Zhao Y, Zeng T . Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing. Nat Methods. 2007; 4(8):651-7. DOI: 10.1038/nmeth1068. View

Compeau P, Pevzner P, Tesler G . How to apply de Bruijn graphs to genome assembly. Nat Biotechnol. 2011; 29(11):987-91. PMC: 5531759. DOI: 10.1038/nbt.2023. View

10.

Frazer K, Murray S, Schork N, Topol E . Human genetic variation and its contribution to complex traits. Nat Rev Genet. 2009; 10(4):241-51. DOI: 10.1038/nrg2554. View

11.

Ward L, Kellis M . HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2011; 40(Database issue):D930-4. PMC: 3245002. DOI: 10.1093/nar/gkr917. View

12.

Kent W, Sugnet C, Furey T, Roskin K, Pringle T, Zahler A . The human genome browser at UCSC. Genome Res. 2002; 12(6):996-1006. PMC: 186604. DOI: 10.1101/gr.229102. View

13.

Jiang R, Yang H, Zhou L, Kuo C, Sun F, Chen T . Sequence-based prioritization of nonsynonymous single-nucleotide polymorphisms for the study of disease mutations. Am J Hum Genet. 2007; 81(2):346-60. PMC: 1950793. DOI: 10.1086/519747. View

14.

Pennacchio L, Ahituv N, Moses A, Prabhakar S, Nobrega M, Shoukry M . In vivo enhancer analysis of human conserved non-coding sequences. Nature. 2006; 444(7118):499-502. DOI: 10.1038/nature05295. View

15.

Marcais G, Kingsford C . A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011; 27(6):764-70. PMC: 3051319. DOI: 10.1093/bioinformatics/btr011. View

16.

. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012; 489(7414):57-74. PMC: 3439153. DOI: 10.1038/nature11247. View

17.

McCarthy M, Abecasis G, Cardon L, Goldstein D, Little J, Ioannidis J . Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008; 9(5):356-69. DOI: 10.1038/nrg2344. View

18.

Visel A, Prabhakar S, Akiyama J, Shoukry M, Lewis K, Holt A . Ultraconservation identifies a small subset of extremely constrained developmental enhancers. Nat Genet. 2008; 40(2):158-60. PMC: 2647775. DOI: 10.1038/ng.2007.55. View

19.

Diaz-Uriarte R, de Andres S . Gene selection and classification of microarray data using random forest. BMC Bioinformatics. 2006; 7:3. PMC: 1363357. DOI: 10.1186/1471-2105-7-3. View

20.

Stranger B, Stahl E, Raj T . Progress and promise of genome-wide association studies for human complex trait genetics. Genetics. 2010; 187(2):367-83. PMC: 3030483. DOI: 10.1534/genetics.110.120907. View