Protein quantitative trait loci identify novel candidates modulating cellular response to chemotherapy

Promises and Challenges of populational Proteomics in Health and Disease

Advances in proteomic assay technologies have significantly increased coverage and throughput, enabling recent increases in the number of large-scale population-based proteomic studies of human plasma and serum. Improvements in multiplexed protein assays have facilitated the quantification of thousands of proteins over a large dynamic range, a key requirement for detecting the lowest-ranging, and potentially the most disease-relevant, blood-circulating proteins. In this perspective, we examine how populational proteomic datasets in conjunction with other concurrent omic measures can be leveraged to better understand the genomic and non-genomic correlates of the soluble proteome, constructing biomarker panels for disease prediction, among others. Mass spectrometry workflows are discussed as they are becoming increasingly competitive with affinity-based array platforms in terms of speed, cost, and proteome coverage due to advances in both instrumentation and workflows. Despite much success, there remain considerable challenges such as orthogonal validation and absolute quantification. We also highlight emergent challenges associated with study design, analytical considerations, and data integration as population-scale studies are run in batches and may involve longitudinal samples collated over many years. Lastly, we take a look at the future of what the nascent next-generation proteomic technologies might provide to the analysis of large sets of blood samples, as well as the difficulties in designing large-scale studies that will likely require participation from multiple and complex funding sources and where data sharing, study designs, and financing must be solved.

Comprehensive information-based differential gene regulatory networks analysis (CIdrgn): Application to gastric cancer and chemotherapy-responsive gene network identification

Biological condition-responsive gene network analysis has attracted considerable research attention because of its ability to identify pathways or gene modules involved in the underlying mechanisms of diseases. Although many condition-specific gene network identification methods have been developed, they are based on partial or incomplete gene regulatory network information, with most studies only considering the differential expression levels or correlations among genes. However, a single gene-based analysis cannot effectively identify the molecular interactions involved in the mechanisms underlying diseases, which reflect perturbations in specific molecular network functions rather than disorders of a single gene. To comprehensively identify differentially regulated gene networks, we propose a novel computational strategy called comprehensive analysis of differential gene regulatory networks (CIdrgn). Our strategy incorporates comprehensive information on the networks between genes, including the expression levels, edge structures and regulatory effects, to measure the dissimilarity among networks. We extended the proposed CIdrgn to cell line characteristic-specific gene network analysis. Monte Carlo simulations showed the effectiveness of CIdrgn for identifying differentially regulated gene networks with different network structures and scales. Moreover, condition-responsive network identification in cell line characteristic-specific gene network analyses was verified. We applied CIdrgn to identify gastric cancer and itsf chemotherapy (capecitabine and oxaliplatin) -responsive network based on the Cancer Dependency Map. The CXC family of chemokines and cadherin gene family networks were identified as gastric cancer-specific gene regulatory networks, which was verified through a literature survey. The networks of the olfactory receptor family with the ASCL1/FOS family were identified as capecitabine- and oxaliplatin sensitive -specific gene networks. We expect that the proposed CIdrgn method will be a useful tool for identifying crucial molecular interactions involved in the specific biological conditions of cancer cell lines, such as the cancer stage or acquired anticancer drug resistance.

Genome-wide genotype-serum proteome mapping provides insights into the cross-ancestry differences in cardiometabolic disease susceptibility

Identification of protein quantitative trait loci (pQTL) helps understand the underlying mechanisms of diseases and discover promising targets for pharmacological intervention. For most important class of drug targets, genetic evidence needs to be generalizable to diverse populations. Given that the majority of the previous studies were conducted in European ancestry populations, little is known about the protein-associated genetic variants in East Asians. Based on data-independent acquisition mass spectrometry technique, we conduct genome-wide association analyses for 304 unique proteins in 2,958 Han Chinese participants. We identify 195 genetic variant-protein associations. Colocalization and Mendelian randomization analyses highlight 60 gene-protein-phenotype associations, 45 of which (75%) have not been prioritized in Europeans previously. Further cross-ancestry analyses uncover key proteins that contributed to the differences in the obesity-induced diabetes and coronary artery disease susceptibility. These findings provide novel druggable proteins as well as a unique resource for the trans-ancestry evaluation of protein-targeted drug discovery.

Proteome-wide Systems Genetics to Identify Functional Regulators of Complex Traits

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

Population-scale proteome variation in human induced pluripotent stem cells

Human disease phenotypes are driven primarily by alterations in protein expression and/or function. To date, relatively little is known about the variability of the human proteome in populations and how this relates to variability in mRNA expression and to disease loci. Here, we present the first comprehensive proteomic analysis of human induced pluripotent stem cells (iPSC), a key cell type for disease modelling, analysing 202 iPSC lines derived from 151 donors, with integrated transcriptome and genomic sequence data from the same lines. We characterised the major genetic and non-genetic determinants of proteome variation across iPSC lines and assessed key regulatory mechanisms affecting variation in protein abundance. We identified 654 protein quantitative trait loci (pQTLs) in iPSCs, including disease-linked variants in protein-coding sequences and variants with trans regulatory effects. These include pQTL linked to GWAS variants that cannot be detected at the mRNA level, highlighting the utility of dissecting pQTL at peptide level resolution.

Comprehensive Cis-Regulation Analysis of Genetic Variants in Human Lymphoblastoid Cell Lines

Genetic variants can influence the expression of mRNA and protein. Genetic regulatory loci such as expression quantitative trait loci (eQTLs) and protein quantitative trait loci (pQTLs) exist in several species. However, it remains unclear how human genetic variants regulate mRNA and protein expression. Here, we characterized six mechanistic models for the genetic regulatory patterns of single-nucleotide polymorphisms (SNPs) and their actions on post-transcriptional expression. Data from Yoruba HapMap lymphoblastoid cell lines were analyzed to identify human cis-eQTLs and pQTLs, as well as protein-specific QTLs (psQTLs). Our results indicated that genetic regulatory loci primarily affected mRNA and protein abundance in patterns where the two were well-correlated. While this finding was observed in both humans and mice (57.5% and 70.3%, respectively), the genetic regulatory patterns differed between species, implying evolutionary differences. Mouse SNPs generally targeted changes in transcript expression (51%), whereas in humans, they largely regulated protein abundance, independent of transcription levels (55.9%). The latter independent function can be explained by psQTLs. Our analysis suggests that local functional genetic variants in the human genome mainly modulate protein abundance independent of mRNA levels through post-transcriptional mechanisms. These findings clarify the impact of genetic variation on phenotype, which is of particular relevance to disease risk and treatment response.

Identification of novel methylated DNA marker ZNF569 for head and neck squamous cell carcinoma

Aberrant DNA methylation pattern plays an indispensable role in the initiation and development of head and neck squamous cell carcinoma (HNSCC). It is well recognized that lymph node metastasis is closely with unfavorable prognosis of HNSCC. Therefore, exploring the methylation events accounting for the lymph node metastasis of HNSCC is very important for improving the clinical outcome of HNSCC. Methylation data, RNA-seq data and clinical data were downloaded from The Cancer Genome Atlas (TCGA) and processed using the R package TCGA-Assembler. MethylMix was use for data analysis by integrating both methylation and gene expression data on HNSCC patients with lymph node metastasis and without lymph node metastasis. Pathway analysis was performed on significantly altered genes using ConsensusPathDB. The role of our interested gene zinc figure protein 569 (ZNF569) in HNSCC was further evaluated. Our results identified many novel hypermethylated/hypomethylated genes that might be closely associated with the lymph node metastasis of HNSCC. Pathway analysis revealed that increase in methylation of genes involved in generic transcription pathway including zinc figure proteins. ZNF569 was hypermethylated in HNSCC tissues especially those with lymph node metastasis. In addition, the expression levels of ZNF569 mRNA and protein were significantly lower in HNSCC tissues and cell lines compared to their respective controls. Moreover, overexpression of ZNF569 inhibited the proliferation, migration and invasion of HNSCC cells. HNSCC patients with lower ZNF569 expression suffered a significantly shorter overall survival than those with higher ZNF569 expression. In conclusion, we have identified many novel differentially methylated genes that might be important for the lymph node metastasis of HNSCC. In addition, ZNF569 might play a tumor suppressive role in carcinogenesis of HNSCC.

The human brainome: network analysis identifies HSPA2 as a novel Alzheimer’s disease target

Our hypothesis is that changes in gene and protein expression are crucial to the development of late-onset Alzheimer’s disease. Previously we examined how DNA alleles control downstream expression of RNA transcripts and how those relationships are changed in late-onset Alzheimer’s disease. We have now examined how proteins are incorporated into networks in two separate series and evaluated our outputs in two different cell lines. Our pipeline included the following steps: (i) predicting expression quantitative trait loci; (ii) determining differential expression; (iii) analysing networks of transcript and peptide relationships; and (iv) validating effects in two separate cell lines. We performed all our analysis in two separate brain series to validate effects. Our two series included 345 samples in the first set (177 controls, 168 cases; age range 65–105; 58% female; KRONOSII cohort) and 409 samples in the replicate set (153 controls, 141 cases, 115 mild cognitive impairment; age range 66–107; 63% female; RUSH cohort). Our top target is heat shock protein family A member 2 (HSPA2), which was identified as a key driver in our two datasets. HSPA2 was validated in two cell lines, with overexpression driving further elevation of amyloid-β40 and amyloid-β42 levels in APP mutant cells, as well as significant elevation of microtubule associated protein tau and phosphorylated-tau in a modified neuroglioma line. This work further demonstrates that studying changes in gene and protein expression is crucial to understanding late onset disease and further nominates HSPA2 as a specific key regulator of late-onset Alzheimer’s disease processes.10.1093/brain/awy215_video1awy215media15824729224001.

CHO-Omics Review: The Impact of Current and Emerging Technologies on Chinese Hamster Ovary Based Bioproduction

CHO cells are the most prevalent platform for modern bio-therapeutic production. Currently, there are several CHO cell lines used in bioproduction with distinct characteristics and unique genotypes and phenotypes. These differences limit advances in productivity and quality that can be achieved by the most common approaches to bioprocess optimization and cell line engineering. Incorporating omics-based approaches into current bioproduction processes will complement traditional methodologies to maximize gains from CHO engineering and bioprocess improvements. In order to highlight the utility of omics technologies in CHO bioproduction, the authors discuss current applications as well as limitations of genomics, transcriptomics, proteomics, metabolomics, lipidomics, fluxomics, glycomics, and multi-omics approaches and the potential they hold for the future of bioproduction. Multiple omics approaches are currently being used to improve CHO bioprocesses; however, the application of these technologies is still limited. As more CHO-omic datasets become available and integrated into systems models, the authors expect significant gains in product yield and quality. While individual omics technologies provide incremental improvements in bioproduction, the authors will likely see the most significant gains by applying multi-omics and systems biology approaches to individual CHO cell lines.

Exploring genetic associations with ceRNA regulation in the human genome

Competing endogenous RNAs (ceRNAs) are RNA molecules that sequester shared microRNAs (miRNAs) thereby affecting the expression of other targets of the miRNAs. Whether genetic variants in ceRNA can affect its biological function and disease development is still an open question. Here we identified a large number of genetic variants that are associated with ceRNA's function using Geuvaids RNA-seq data for 462 individuals from the 1000 Genomes Project. We call these loci competing endogenous RNA expression quantitative trait loci or 'cerQTL', and found that a large number of them were unexplored in conventional eQTL mapping. We identified many cerQTLs that have undergone recent positive selection in different human populations, and showed that single nucleotide polymorphisms in gene 3΄UTRs at the miRNA seed binding regions can simultaneously regulate gene expression changes in both cis and trans by the ceRNA mechanism. We also discovered that cerQTLs are significantly enriched in traits/diseases associated variants reported from genome-wide association studies in the miRNA binding sites, suggesting that disease susceptibilities could be attributed to ceRNA regulation. Further in vitro functional experiments demonstrated that a cerQTL rs11540855 can regulate ceRNA function. These results provide a comprehensive catalog of functional non-coding regulatory variants that may be responsible for ceRNA crosstalk at the post-transcriptional level.

Pervasive coexpression of spatially proximal genes is buffered at the protein level

Genes are not randomly distributed in the genome. In humans, 10% of protein-coding genes are transcribed from bidirectional promoters and many more are organised in larger clusters. Intriguingly, neighbouring genes are frequently coexpressed but rarely functionally related. Here we show that coexpression of bidirectional gene pairs, and closeby genes in general, is buffered at the protein level. Taking into account the 3D architecture of the genome, we find that co-regulation of spatially close, functionally unrelated genes is pervasive at the transcriptome level, but does not extend to the proteome. We present evidence that non-functional mRNA coexpression in human cells arises from stochastic chromatin fluctuations and direct regulatory interference between spatially close genes. Protein-level buffering likely reflects a lack of coordination of post-transcriptional regulation of functionally unrelated genes. Grouping human genes together along the genome sequence, or through long-range chromosome folding, is associated with reduced expression noise. Our results support the hypothesis that the selection for noise reduction is a major driver of the evolution of genome organisation.

Methods to analyze big data in pharmacogenomics research

The scale and scope of pharmacogenomics research continues to expand as the cost and efficiency of molecular data generation techniques advance. These new technologies give rise to enormous opportunity for the identification of important genetic and genomic factors important for drug treatment response. With this opportunity come significant challenges. Most of these can be categorized as 'big data' issues, facing not only pharmacogenomics, but other fields in the life sciences as well. In this review, we describe some of the analysis techniques and tools being implemented for genetic/genomic discovery in pharmacogenomics.

Genetic Variants Contributing to Colistin Cytotoxicity: Identification of TGIF1 and HOXD10 Using a Population Genomics Approach

Colistin sulfate (polymixin E) is an antibiotic prescribed with increasing frequency for severe Gram-negative bacterial infections. As nephrotoxicity is a common side effect, the discovery of pharmacogenomic markers associated with toxicity would benefit the utility of this drug. Our objective was to identify genetic markers of colistin cytotoxicity that were also associated with expression of key proteins using an unbiased, whole genome approach and further evaluate the functional significance in renal cell lines. To this end, we employed International HapMap lymphoblastoid cell lines (LCLs) of Yoruban ancestry with known genetic information to perform a genome-wide association study (GWAS) with cellular sensitivity to colistin. Further association studies revealed that single nucleotide polymorphisms (SNPs) associated with gene expression and protein expression were significantly enriched in SNPs associated with cytotoxicity (p ≤ 0.001 for gene and p = 0.015 for protein expression). The most highly associated SNP, chr18:3417240 (p = 6.49 × 10⁻⁸), was nominally a cis-expression quantitative trait locus (eQTL) of the gene TGIF1 (transforming growth factor β (TGFβ)-induced factor-1; p = 0.021) and was associated with expression of the protein HOXD10 (homeobox protein D10; p = 7.17 × 10⁻⁵). To demonstrate functional relevance in a murine colistin nephrotoxicity model, HOXD10 immunohistochemistry revealed upregulated protein expression independent of mRNA expression in response to colistin administration. Knockdown of TGIF1 resulted in decreased protein expression of HOXD10 and increased resistance to colistin cytotoxicity. Furthermore, knockdown of HOXD10 in renal cells also resulted in increased resistance to colistin cytotoxicity, supporting the physiological relevance of the initial genomic associations.

Proteomic Profiling of the Microsomal Root Fraction: Discrimination of Pisum sativum L. Cultivars and Identification of Putative Root Growth Markers

Legumes are a large and economically important family, containing a variety of crop plants. Alongside different cereals, some fruits, and tropical roots, a number of leguminosae evolved for millennia as crops with human society. One of these legumes is Pisum sativum L., the common garden pea. In the past, breeding has been largely selective on improved above-ground organs. However, parameters, such as root-growth, which determines acquisition of nutrients and water, have largely been underestimated. Although the genome of P. sativum is still not fully sequenced, multiple proteomic studies have been published on a variety of physiological aspects in the last years. The presented work focused on the connection between root length and the influence of the microsomal root proteome of four different pea cultivars after five days of germination (cultivar Vroege, Girl from the Rhineland, Kelvedon Wonder, and Blauwschokker). In total, 60 proteins were identified to have significantly differential abundances in the four cultivars. Root growth of five-days old seedlings and their microsomal proteome revealed a similar separation pattern, suggesting that cultivar-specific root growth performance is explained by differential membrane and ribosomal protein levels. Hence, we reveal and discuss several putative root growth protein markers possibly playing a key role for improved primary root growth breeding strategies.

Multidimensional Integrative Genomics Approaches to Dissecting Cardiovascular Disease

Elucidating the mechanisms of complex diseases such as cardiovascular disease (CVD) remains a significant challenge due to multidimensional alterations at molecular, cellular, tissue, and organ levels. To better understand CVD and offer insights into the underlying mechanisms and potential therapeutic strategies, data from multiple omics types (genomics, epigenomics, transcriptomics, metabolomics, proteomics, microbiomics) from both humans and model organisms have become available. However, individual omics data types capture only a fraction of the molecular mechanisms. To address this challenge, there have been numerous efforts to develop integrative genomics methods that can leverage multidimensional information from diverse data types to derive comprehensive molecular insights. In this review, we summarize recent methodological advances in multidimensional omics integration, exemplify their applications in cardiovascular research, and pinpoint challenges and future directions in this incipient field.

Computational discovery of transcription factors associated with drug response

This study integrates gene expression, genotype and drug response data in lymphoblastoid cell lines with transcription factor (TF)-binding sites from ENCODE (Encyclopedia of Genomic Elements) in a novel methodology that elucidates regulatory contexts associated with cytotoxicity. The method, GENMi (Gene Expression iN the Middle), postulates that single-nucleotide polymorphisms within TF-binding sites putatively modulate its regulatory activity, and the resulting variation in gene expression leads to variation in drug response. Analysis of 161 TFs and 24 treatments revealed 334 significantly associated TF-treatment pairs. Investigation of 20 selected pairs yielded literature support for 13 of these associations, often from studies where perturbation of the TF expression changes drug response. Experimental validation of significant GENMi associations in taxanes and anthracyclines across two triple-negative breast cancer cell lines corroborates our findings. The method is shown to be more sensitive than an alternative, genome-wide association study-based approach that does not use gene expression. These results demonstrate the utility of GENMi in identifying TFs that influence drug response and provide a number of candidates for further testing.

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

Linking short tandem repeat polymorphisms with cytosine modifications in human lymphoblastoid cell lines

Inter-individual variation in cytosine modifications has been linked to complex traits in humans. Cytosine modification variation is partially controlled by single nucleotide polymorphisms (SNPs), known as modified cytosine quantitative trait loci (mQTL). However, little is known about the role of short tandem repeat polymorphisms (STRPs), a class of structural genetic variants, in regulating cytosine modifications. Utilizing the published data on the International HapMap Project lymphoblastoid cell lines (LCLs), we assessed the relationships between 721 STRPs and the modification levels of 283,540 autosomal CpG sites. Our findings suggest that, in contrast to the predominant cis-acting mode for SNP-based mQTL, STRPs are associated with cytosine modification levels in both cis-acting (local) and trans-acting (distant) modes. In local scans within the ±1 Mb windows of target CpGs, 21, 9, and 21 cis-acting STRP-based mQTL were detected in CEU (Caucasian residents from Utah, USA), YRI (Yoruba people from Ibadan, Nigeria), and the combined samples, respectively. In contrast, 139,420, 76,817, and 121,866 trans-acting STRP-based mQTL were identified in CEU, YRI, and the combined samples, respectively. A substantial proportion of CpG sites detected with local STRP-based mQTL were not associated with SNP-based mQTL, suggesting that STRPs represent an independent class of mQTL. Functionally, genetic variants neighboring CpG-associated STRPs are enriched with genome-wide association study (GWAS) loci for a variety of complex traits and diseases, including cancers, based on the National Human Genome Research Institute (NHGRI) GWAS Catalog. Therefore, elucidating these STRP-based mQTL in addition to SNP-based mQTL can provide novel insights into the genetic architectures of complex traits.

In vitro human cell line models to predict clinical response to anticancer drugs

In vitro human cell line models have been widely used for cancer pharmacogenomic studies to predict clinical response, to help generate pharmacogenomic hypothesis for further testing, and to help identify novel mechanisms associated with variation in drug response. Among cell line model systems, immortalized cell lines such as Epstein-Barr virus (EBV)-transformed lymphoblastoid cell lines (LCLs) have been used most often to test the effect of germline genetic variation on drug efficacy and toxicity. Another model, especially in cancer research, uses cancer cell lines such as the NCI-60 panel. These models have been used mainly to determine the effect of somatic alterations on response to anticancer therapy. Even though these cell line model systems are very useful for initial screening, results from integrated analyses of multiple omics data and drug response phenotypes using cell line model systems still need to be confirmed by functional validation and mechanistic studies, as well as validation studies using clinical samples. Future models might include the use of patient-specific inducible pluripotent stem cells and the incorporation of 3D culture which could further optimize in vitro cell line models to improve their predictive validity.

Pharmacoethnicity in Paclitaxel-Induced Sensory Peripheral Neuropathy

PURPOSE

Paclitaxel is used worldwide in the treatment of breast, lung, ovarian, and other cancers. Sensory peripheral neuropathy is an associated adverse effect that cannot be predicted, prevented, or mitigated. To better understand the contribution of germline genetic variation to paclitaxel-induced peripheral neuropathy, we undertook an integrative approach that combines genome-wide association study (GWAS) data generated from HapMap lymphoblastoid cell lines (LCL) and Asian patients.

METHODS

GWAS was performed with paclitaxel-induced cytotoxicity generated in 363 LCLs and with paclitaxel-induced neuropathy from 145 Asian patients. A gene-based approach was used to identify overlapping genes and compare with a European clinical cohort of paclitaxel-induced neuropathy. Neurons derived from human-induced pluripotent stem cells were used for functional validation of candidate genes.

RESULTS

SNPs near AIPL1 were significantly associated with paclitaxel-induced cytotoxicity in Asian LCLs (P < 10(-6)). Decreased expression of AIPL1 resulted in decreased sensitivity of neurons to paclitaxel by inducing neurite morphologic changes as measured by increased relative total outgrowth, number of processes and mean process length. Using a gene-based analysis, there were 32 genes that overlapped between Asian LCL cytotoxicity and Asian patient neuropathy (P < 0.05), including BCR. Upon BCR knockdown, there was an increase in neuronal sensitivity to paclitaxel as measured by neurite morphologic characteristics.

CONCLUSIONS

We identified genetic variants associated with Asian paclitaxel-induced cytotoxicity and functionally validated the AIPL1 and BCR in a neuronal cell model. Furthermore, the integrative pharmacogenomics approach of LCL/patient GWAS may help prioritize target genes associated with chemotherapeutic-induced peripheral neuropathy.

Approaches for establishing the function of regulatory genetic variants involved in disease

The diversity of regulatory genetic variants and their mechanisms of action reflect the complexity and context-specificity of gene regulation. Regulatory variants are important in human disease and defining such variants and establishing mechanism is crucial to the interpretation of disease-association studies. This review describes approaches for identifying and functionally characterizing regulatory variants, illustrated using examples from common diseases. Insights from recent advances in resolving the functional epigenomic regulatory landscape in which variants act are highlighted, showing how this has enabled functional annotation of variants and the generation of hypotheses about mechanism of action. The utility of quantitative trait mapping at the transcript, protein and metabolite level to define association of specific genes with particular variants and further inform disease associations are reviewed. Establishing mechanism of action is an essential step in resolving functional regulatory variants, and this review describes how this is being facilitated by new methods for analyzing allele-specific expression, mapping chromatin interactions and advances in genome editing. Finally, integrative approaches are discussed together with examples highlighting how defining the mechanism of action of regulatory variants and identifying specific modulated genes can maximize the translational utility of genome-wide association studies to understand the pathogenesis of diseases and discover new drug targets or opportunities to repurpose existing drugs to treat them.