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.

A frame-based representation of genomic sequences for removing errors and rare variant detection in NGS data

Pro-Frame: similarity-based gene recognition in eukaryotic DNA sequences with errors

Multigems: Detection Of Snvs From Multiple Samples Using Model Selection On High-Throughput Sequencing Data

Error filtering, pair assembly and error correction for next-generation sequencing reads

BreaKmer: detection of structural variation in targeted massively parallel sequencing data using kmers

Empirical Bayes single nucleotide variant-calling for next-generation sequencing data

A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

GRAMEP: an alignment-free method based on the Maximum Entropy Principle for identifying SNPs

Too many needles in this haystack: algorithms for the analysis of next generation sequence data

Enhanced Detection and Genotyping of Disease-Associated Tandem Repeats Using HMMSTR and Targeted Long-Read Sequencing

A novel signal processing approach for the detection of copy number variations in the human genome

Accurate Typing of Short Tandem Repeats from Genome-Wide Sequencing Data and Its Applications.

RES-Scanner: a Software Package for Genome-Wide Identification of RNA-editing Sites

Interlacing Personal and Reference Genomes for Machine Learning Disease-Variant Detection

Turn ‘noise’ to signal: accurately rectify millions of erroneous short reads through graph learning on edit distances

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

SRBreak: A Read-Depth and Split-Read Framework to Identify Breakpoints of Different Events Inside Simple Copy-Number Variable Regions

RealSeq2: a Software Integrated with UMI Identification, Error Correction, and Methylation Modifications Storing

Comprehensive assessment of error correction methods for high-throughput sequencing data

HaplotagLR: An efficient and configurable utility for haplotagging long reads