Haplotype-resolved and integrated genome analysis of the cancer cell line HepG2

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

Discovery of Missing Proteins from an Aneuploidy Cell Line Using a Proteogenomic Approach

With the steadfast development of proteomic technology, the number of missing proteins (MPs) has been continuously shrinking, with approximately 1470 MPs that have not been explored yet. Due to this phenomenon, the discovery of MPs has been increasingly more difficult and elusive. In order to face this challenge, we have hypothesized that a stable aneuploid cell line with increased chromosomes serves as a useful material for assisting MP exploration. Ker-CT cell line with trisomy at chromosome 5 and 20 was selected for this purpose. With a combination strategy of RNA-Seq and LC-MS/MS, a total of 22 178 transcripts and 8846 proteins were identified in Ker-CT. Although the transcripts corresponding to 15 and 15 MP genes located at chromosome 5 and 20 were detected, none of the MPs were found in Ker-CT. Surprisingly, 3 MPs containing at least two unique non-nest peptides of length ≥9 amino acids were identified in Ker-CT, whose genes are located on chromosome 3 and 10, respectively. Furthermore, the 3 MPs were verified using the method of parallel reaction monitoring (PRM). These results suggest that the abnormal status of chromosomes may not only impact the expression of the corresponding genes in trisomy chromosomes, but also influence that of other chromosomes, which benefits MP discovery. The data obtained in this study are available via ProteomeXchange (PXD028647) and PeptideAtlas (PASS01700), respectively.

Antisense-Mediated Down-Regulation of Factor V-Short Splicing in a Liver Cell Line Model

Coagulation factor V (FV) is a liver-derived protein encoded by the F5 gene. Alternative splicing of F5 exon 13 produces a low-abundance splicing isoform, known as FV-short, which binds the anticoagulant protein tissue factor pathway inhibitor (TFPIα) with high affinity, stabilising it in the circulation and potently enhancing its anticoagulant activity. Accordingly, rare F5 gene mutations that up-regulate FV-short splicing are associated with bleeding. In this study we have explored the possibility of decreasing FV-short splicing by antisense-based splicing modulation. To this end, we have designed morpholino antisense oligonucleotides (MAOs) targeting the FV-short-specific donor and acceptor splice sites and tested their efficacy in a liver cell line (HepG2) that naturally expresses full-length FV and FV-short. Cells were treated with 0–20 µM MAO, and full-length FV and FV-short mRNA expression was analysed by RT-(q)PCR. Both MAOs, alone or in combination, decreased the FV-short/full-length FV mRNA ratio down to ~50% of its original value in a specific and dose-dependent manner. This pilot study provides proof-of-principle for the possibility to decrease FV-short expression by antisense-mediated splicing modulation. In turn, this may form the basis for novel therapeutic approaches to bleeding disorders caused by FV-short over-expression and/or elevated TFPIα (activity) levels.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.

New Insights on Molecular Mechanism of Hepatitis B Virus Covalently Closed Circular DNA Formation

The chronic factor of the Hepatitis B Virus (HBV), specifically the covalently closed circular DNA (cccDNA), is a highly stable and active viral episomal genome established in the livers of chronic hepatitis B patients as a constant source of disease. Being able to target and eliminate cccDNA is the end goal for a genuine cure for HBV. Yet how HBV cccDNA is formed from the viral genomic relaxed circular DNA (rcDNA) and by what host factors had been long-standing research questions. It is generally acknowledged that HBV hijacks cellular functions to turn the open circular DNA conformation of rcDNA into cccDNA through DNA repair mechanisms. With great efforts from the HBV research community, there have been several recent leaps in our understanding of cccDNA formation. It is our goal in this review to analyze the recent reports showing evidence of cellular factor's involvement in the molecular pathway of cccDNA biosynthesis.

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Background

DNA methylation is an important epigenetic modification that plays a critical role in most eukaryotic organisms. Parental alleles in haploid genomes may exhibit different methylation patterns, which can lead to different phenotypes and even different therapeutic and drug responses to diseases. However, to our knowledge, no software is available for the identification of DNA methylation haplotype regions with combined allele-specific DNA methylation, single nucleotide polymorphisms (SNPs) and high-throughput chromosome conformation capture (Hi-C) data.

Results

In this paper, we developed a new method, MethHaplo, that identify DNA methylation haplotype regions with allele-specific DNA methylation and SNPs from whole-genome bisulfite sequencing (WGBS) data. Our results showed that methylation haplotype regions were ten times longer than haplotypes with SNPs only. When we integrate WGBS and Hi-C data, MethHaplo could call even longer haplotypes.

Conclusions

This study illustrates the usefulness of methylation haplotypes. By constructing methylation haplotypes for various cell lines, we provide a clearer picture of the effect of DNA methylation on gene expression, histone modification and three-dimensional chromosome structure at the haplotype level. Our method could benefit the study of parental inheritance-related disease and hybrid vigor in agriculture.

A random forest-based framework for genotyping and accuracy assessment of copy number variations

Detection of copy number variations (CNVs) is essential for uncovering genetic factors underlying human diseases. However, CNV detection by current methods is prone to error, and precisely identifying CNVs from paired-end whole genome sequencing (WGS) data is still challenging. Here, we present a framework, CNV-JACG, for Judging the Accuracy of CNVs and Genotyping using paired-end WGS data. CNV-JACG is based on a random forest model trained on 21 distinctive features characterizing the CNV region and its breakpoints. Using the data from the 1000 Genomes Project, Genome in a Bottle Consortium, the Human Genome Structural Variation Consortium and in-house technical replicates, we show that CNV-JACG has superior sensitivity over the latest genotyping method, SV², particularly for the small CNVs (≤1 kb). We also demonstrate that CNV-JACG outperforms SV² in terms of Mendelian inconsistency in trios and concordance between technical replicates. Our study suggests that CNV-JACG would be a useful tool in assessing the accuracy of CNVs to meet the ever-growing needs for uncovering the missing heritability linked to CNVs.

Fully Phased Sequence of a Diploid Human Genome Determined de Novo from the DNA of a Single Individual

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

Structural variant analysis for linked-read sequencing data with gemtools

SUMMARY

Linked-read sequencing generates synthetic long reads which are useful for the detection and analysis of structural variants (SVs). The software associated with 10× Genomics linked-read sequencing, Long Ranger, generates the essential output files (BAM, VCF, SV BEDPE) necessary for downstream analyses. However, to perform downstream analyses requires the user to customize their own tools to handle the unique features of linked-read sequencing data. Here, we describe gemtools, a collection of tools for the downstream and in-depth analysis of SVs from linked-read data. Gemtools uses the barcoded aligned reads and the Megabase-scale phase blocks to determine haplotypes of SV breakpoints and delineate complex breakpoint configurations at the resolution of single DNA molecules. The gemtools package is a suite of tools that provides the user with the flexibility to perform basic functions on their linked-read sequencing output in order to address even more questions.

AVAILABILITY AND IMPLEMENTATION

The gemtools package is freely available for download at: https://github.com/sgreer77/gemtools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Allele-specific DNA methylation is increased in cancers and its dense mapping in normal plus neoplastic cells increases the yield of disease-associated regulatory SNPs

BACKGROUND

Mapping of allele-specific DNA methylation (ASM) can be a post-GWAS strategy for localizing regulatory sequence polymorphisms (rSNPs). The advantages of this approach, and the mechanisms underlying ASM in normal and neoplastic cells, remain to be clarified.

RESULTS

We perform whole genome methyl-seq on diverse normal cells and tissues and three cancer types. After excluding imprinting, the data pinpoint 15,112 high-confidence ASM differentially methylated regions, of which 1838 contain SNPs in strong linkage disequilibrium or coinciding with GWAS peaks. ASM frequencies are increased in cancers versus matched normal tissues, due to widespread allele-specific hypomethylation and focal allele-specific hypermethylation in poised chromatin. Cancer cells show increased allele switching at ASM loci, but disruptive SNPs in specific classes of CTCF and transcription factor binding motifs are similarly correlated with ASM in cancer and non-cancer. Rare somatic mutations affecting these same motif classes track with de novo ASM. Allele-specific transcription factor binding from ChIP-seq is enriched among ASM loci, but most ASM differentially methylated regions lack such annotations, and some are found in otherwise uninformative "chromatin deserts."

CONCLUSIONS

ASM is increased in cancers but occurs by a shared mechanism involving disruptive SNPs in CTCF and transcription factor binding sites in both normal and neoplastic cells. Dense ASM mapping in normal plus cancer samples reveals candidate rSNPs that are difficult to find by other approaches. Together with GWAS data, these rSNPs can nominate specific transcriptional pathways in susceptibility to autoimmune, cardiometabolic, neuropsychiatric, and neoplastic diseases.

Concentration-dependent toxicogenomic changes of silver nanoparticles in hepatocyte-like cells derived from human induced pluripotent stem cells

The application of silver nanoparticles (AgNPs) in consumer products has been increasing rapidly over the past decades. Therefore, in vitro models capable of accurately predicting the toxicity of AgNPs are much needed. Hepatocyte-like cells (HLCs) derived from human induced pluripotent stem cells (iPSCs) represent an attractive alternative in vitro hepatotoxicity model. Yet, the use of iPSC-derived HLCs (iPSC-HLCs) for the study of nanoparticle toxicity has not been reported so far. In the present study, transcriptomic changes induced by varying concentrations (5-25 μg/ml) of AgNPs were characterized in iPSC-HLCs after 24-h exposure. AgNPs caused concentration-dependent gene expression changes in iPSC-HLCs. At all the concentrations, members of the metallothionein (MT) and the heat shock protein (HSP) families were the dominating upregulated genes, suggesting that exposure to AgNPs induced oxidative stresses in iPSC-HLCs and as a result elicited cellular protective responses in the cells. Functional analysis showed that the differentially expressed genes (DEGs) were majorly involved in the biological processes of metabolism, response to stress, and cell organization and biogenesis. Ingenuity Pathway Analysis revealed that cancer was at the top of diseases and disorders associated with the DEGs at all concentrations. These results were in accordance with those reported previously on hepatoma cell lines and primary hepatocytes. Considering the advantages iPSC-HLCs have over other liver cell models in terms of unlimited supply, consistency in quality, sustainability of function in long-term culture, and, more importantly, affordability of donor specificity, the results of the current study suggest that iPSC-HLCs may serve as a better in vitro model for liver nanotoxicology.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Refined detection and phasing of structural aberrations in pediatric acute lymphoblastic leukemia by linked-read whole-genome sequencing

Structural chromosomal rearrangements that can lead to in-frame gene-fusions are a leading source of information for diagnosis, risk stratification, and prognosis in pediatric acute lymphoblastic leukemia (ALL). Traditional methods such as karyotyping and FISH struggle to accurately identify and phase such large-scale chromosomal aberrations in ALL genomes. We therefore evaluated linked-read WGS for detecting chromosomal rearrangements in primary samples of from 12 patients diagnosed with ALL. We assessed the effect of input DNA quality on phased haplotype block size and the detectability of copy number aberrations and structural variants in the ALL genomes. We found that biobanked DNA isolated by standard column-based extraction methods was sufficient to detect chromosomal rearrangements even at low 10x sequencing coverage. Linked-read WGS enabled precise, allele-specific, digital karyotyping at a base-pair resolution for a wide range of structural variants including complex rearrangements and aneuploidy assessment. With use of haplotype information from the linked-reads, we also identified previously unknown structural variants, such as a compound heterozygous deletion of ERG in a patient with the DUX4-IGH fusion gene. We conclude that linked-read WGS allows detection of important pathogenic variants in ALL genomes at a resolution beyond that of traditional karyotyping and FISH.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

rs953413 Regulates Polyunsaturated Fatty Acid Metabolism by Modulating ELOVL2 Expression

Long-chain polyunsaturated fatty acids (LC-PUFAs) influence human health in several areas, including cardiovascular disease, diabetes, fatty liver disease, and cancer. ELOVL2 encodes one of the key enzymes in the in vivo synthesis of LC-PUFAs from their precursors. Variants near ELOVL2 have repeatedly been associated with levels of LC-PUFA-derived metabolites in genome-wide association studies (GWAS), but the mechanisms behind these observations remain poorly defined. In this study, we found that rs953413, located in the first intron of ELOVL2, lies within a functional FOXA and HNF4α cooperative binding site. The G allele of rs953413 increases binding of FOXA1/FOXA2 and HNF4α to an evolutionarily conserved enhancer element, conferring allele-specific upregulation of the rs953413-associated gene ELOVL2. The expression of ELOVL2 was significantly downregulated by both FOXA1 and HNF4α knockdown and CRISPR/Cas9-mediated direct mutation to the enhancer element. Our results suggest that rs953413 regulates LC-PUFAs metabolism by altering ELOVL2 expression through FOXA1/FOXA2 and HNF4α cooperation.

•

rs953413 resides in an evolutionarily conserved enhancer region

•

rs953413 mediates the cooperative binding of FOXA and HNF4α to the enhancer region

•

The rs953413 locus plays a key role in regulating ELOVL2 expression

•

rs953413 is implicated in PUFA metabolism by regulating ELOVL2 expression

DNA Analysis by Restriction Enzyme (DARE) enables concurrent genomic and epigenomic characterization of single cells

Genome-wide profiling of copy number alterations and DNA methylation in single cells could enable detailed investigation into the genomic and epigenomic heterogeneity of complex cell populations. However, current methods to do this require complex sample processing and cleanup steps, lack consistency, or are biased in their genomic representation. Here, we describe a novel single-tube enzymatic method, DNA Analysis by Restriction Enzyme (DARE), to perform deterministic whole genome amplification while preserving DNA methylation information. This method was evaluated on low amounts of DNA and single cells, and provides accurate copy number aberration calling and representative DNA methylation measurement across the whole genome. Single-cell DARE is an attractive and scalable approach for concurrent genomic and epigenomic characterization of cells in a heterogeneous population.

Allele-specific binding of RNA-binding proteins reveals functional genetic variants in the RNA

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs.