Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a ...systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.
Fibroblast growth factor 21 (FGF21) is a hormone that has insulin-sensitizing properties. Some trials of FGF21 analogs show weight loss and lipid-lowering effects. Recent studies have shown that a ...common allele in the FGF21 gene alters the balance of macronutrients consumed, but there was little evidence of an effect on metabolic traits. We studied a common FGF21 allele (A:rs838133) in 451,099 people from the UK Biobank study, aiming to use the human allele to inform potential adverse and beneficial effects of targeting FGF21. We replicated the association between the A allele and higher percentage carbohydrate intake. We then showed that this allele is more strongly associated with higher blood pressure and waist-hip ratio, despite an association with lower total body-fat percentage, than it is with BMI or type 2 diabetes. These human phenotypes of variation in the FGF21 gene will inform research into FGF21’s mechanisms and therapeutic potential.
Display omitted
•A human allele in FGF21 is associated with higher alcohol as well as sugar intake•Phenotypes associated with this allele likely mimic those of lower FGF21 function•The human allele in FGF21 is associated with lower total body-fat percentage•The human allele in FGF21 is associated with higher blood pressure and waist-hip ratio
Drugs targeting the hormone FGF21 may have beneficial health effects. Variations in human DNA in the FGF21 gene provide an indication of what those effects may be. Here, we show that variation in the FGF21 gene is associated with higher blood pressure and altered body shape, despite lower total body-fat percentage.
Abstract
In this study, we generated whole-transcriptome RNA-Seq from
n
= 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL ...(microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.
Hypertension is a common toxicity induced by bevacizumab and other antiangiogenic drugs. There are no biomarkers to predict the risk of bevacizumab-induced hypertension. This study aimed to identify ...plasma proteins related to the function of the vasculature to predict the risk of severe bevacizumab-induced hypertension. Using pretreated plasma samples from 398 bevacizumab-treated patients in two clinical trials (CALGB 80303 and 90401), the levels of 17 proteins were measured via ELISA. The association between proteins and grade 3 bevacizumab-induced hypertension was performed by calculating the odds ratio (OR) from logistic regression adjusting for age, sex, and clinical trial. Using the optimal cut-point of each protein, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) for hypertension were estimated. Five proteins showed no difference in levels between clinical trials and were used for analyses. Lower levels of angiopoietin-2 (
p
= 0.0013, OR 3.41, 95% CI 1.67–7.55), VEGF-A (
p
= 0.0008, OR 4.25, 95% CI 1.93–10.72), and VCAM-1 (
p
= 0.0067, OR 2.68, 95% CI 1.34–5.63) were associated with an increased risk of grade 3 hypertension. The multivariable model suggests independent effects of angiopoietin-2 (
p
= 0.0111, OR 2.71, 95% CI 1.29–6.10), VEGF-A (
p
= 0.0051, OR 3.66, 95% CI 1.54–9.73), and VCAM-1 (
p
= 0.0308, OR 2.27, 95% CI 1.10–4.92). The presence of low levels of 2–3 proteins had an OR of 10.06 (95% CI 3.92–34.18,
p
= 1.80 × 10
–5
) for the risk of hypertension, with sensitivity of 89.7%, specificity of 53.5%, PPV of 17.3%, and NPV of 97.9%. This is the first study providing evidence of plasma proteins with potential value to predict patients at risk of developing bevacizumab-induced hypertension.
Clinical trial registration:
ClinicalTrials.gov Identifier: NCT00088894 (CALGB 80303); and NCT00110214 (CALGB 90401).
Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including ...therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD,
gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and
that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls
transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and
promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.
Understanding the molecular mechanisms of complex traits is essential for developing targeted interventions. We analyzed liver expression quantitative-trait locus (eQTL) meta-analysis data on 1,183 ...participants to identify conditionally distinct signals. We found 9,013 eQTL signals for 6,564 genes; 23% of eGenes had two signals, and 6% had three or more signals. We then integrated the eQTL results with data from 29 cardiometabolic genome-wide association study (GWAS) traits and identified 1,582 GWAS-eQTL colocalizations for 747 eGenes. Non-primary eQTL signals accounted for 17% of all colocalizations. Isolating signals by conditional analysis prior to coloc resulted in 37% more colocalizations than using marginal eQTL and GWAS data, highlighting the importance of signal isolation. Isolating signals also led to stronger evidence of colocalization: among 343 eQTL-GWAS signal pairs in multi-signal regions, analyses that isolated the signals of interest resulted in higher posterior probability of colocalization for 41% of tests. Leveraging allelic heterogeneity, we predicted causal effects of gene expression on liver traits for four genes. To predict functional variants and regulatory elements, we colocalized eQTL with liver chromatin accessibility QTL (caQTL) and found 391 colocalizations, including 73 with non-primary eQTL signals and 60 eQTL signals that colocalized with both a caQTL and a GWAS signal. Finally, we used publicly available massively parallel reporter assays in HepG2 to highlight 14 eQTL signals that include at least one expression-modulating variant. This multi-faceted approach to unraveling the genetic underpinnings of liver-related traits could lead to therapeutic development.
Understanding molecular mechanisms of complex traits is essential for developing targeted interventions. We performed a liver eQTL meta-analysis, identifying 9,013 signals for 6,564 genes. We highlighted the value of analyzing distinct signals. Integrating eQTL with GWAS, chromatin QTL, and transcriptional reporter assays provided insights into genetic underpinnings of liver-related traits.
Drug-induced liver injury (DILI) due to amoxicillin–clavulanate (AC) has been associated with HLA-A∗02:01, HLA-DRB1∗15:01, and rs2476601, a missense variant in PTPN22. The aim of this study was to ...identify novel risk factors for AC-DILI and to construct a genetic risk score (GRS). METHODS: Transcriptome-wide association study and genome-wide association study analyses were performed on 444 AC-DILI cases and 10,397 population-based controls of European descent. Associations were confirmed in a validation cohort (n = 133 cases and 17,836 population-based controls). Discovery and validation AC-DILI cases were also compared with 1358 and 403 non–AC-DILI cases.
Transcriptome-wide association study revealed a significant association of AC-DILI risk with reduced liver expression of ERAP2 (P = 3.7 × 10–7), coding for an aminopeptidase involved in antigen presentation. The lead eQTL single nucleotide polymorphism, rs1363907 (G), was associated with AC-DILI risk in the discovery (odds ratio OR, 1.68; 95% CI, 1.23–1.66; P = 1.7 × 10–7) and validation cohorts (OR, 1.2; 95% CI, 1.04–2.05; P = .03), following a recessive model. We also identified HLA-B∗15:18 as a novel AC-DILI risk factor in both discovery (OR, 4.19; 95% CI, 2.09–8.36; P = 4.9 × 10–5) and validation (OR, 7.78; 95% CI, 2.75–21.99; P = .0001) cohorts. GRS, incorporating rs1363907, rs2476601, HLA-B∗15:18, HLA-A∗02:01, and HLA-DRB1∗15:01, was highly predictive of AC-DILI risk when cases were analyzed against both general population and non–AC-DILI control cohorts. GRS was the most significant predictor in a regression model containing known AC-DILI clinical risk characteristics and significantly improved the predictive model.
We identified novel associations of AC-DILI risk with ERAP2 low expression and with HLA-B∗15:18. GRS based on the 5 risk variants may assist AC-DILI causality assessment and risk management.
▪
Advanced pancreatic cancer is a highly refractory disease almost always associated with survival of little more than a year. New interventions based on novel targets are needed. We aim to identify ...new genetic determinants of overall survival (OS) in patients after treatment with gemcitabine using genome-wide screens of germline DNA. We aim also to support these findings with in vitro functional analysis.
Genome-wide screens of germline DNA in two independent cohorts of pancreatic cancer patients (from the Cancer and Leukemia Group B (CALGB) 80303 and the Mayo Clinic) were used to select new genes associated with OS. The vitamin D receptor gene (VDR) was selected, and the interactions of genetic variation in VDR with circulating vitamin D levels and gemcitabine treatment were evaluated. Functional effects of common VDR variants were also evaluated in experimental assays in human cell lines.
The rs2853564 variant in VDR was associated with OS in patients from both the Mayo Clinic (HR 0.81, 95% CI 0.70-0.94, p = 0.0059) and CALGB 80303 (HR 0.74, 0.63-0.87, p = 0.0002). rs2853564 interacted with high pre-treatment levels of 25-hydroxyvitamin D (25(OH)D, a measure of endogenous vitamin D) (p = 0.0079 for interaction) and with gemcitabine treatment (p = 0.024 for interaction) to confer increased OS. rs2853564 increased transcriptional activity in luciferase assays and reduced the binding of the IRF4 transcription factor.
Our findings propose VDR as a novel determinant of survival in advanced pancreatic cancer patients. Common functional variation in this gene might interact with endogenous vitamin D and gemcitabine treatment to determine improved patient survival. These results support evidence for a modulatory role of the vitamin D pathway for the survival of advanced pancreatic cancer patients.
Irinotecan (CPT-11) is an anticancer agent widely used to treat adult solid tumours. Large interindividual variability in the clearance of irinotecan and SN-38, its active and toxic metabolite, ...results in highly unpredictable toxicity.
In 217 cancer patients treated with intravenous irinotecan single agent or in combination, germline DNA was used to interrogate the variation in 84 genes by next-generation sequencing. A stepwise analytical framework including a population pharmacokinetic model with SNP- and gene-based testing was used to identify demographic/clinical/genetic factors that influence the clearance of irinotecan and SN-38.
Irinotecan clearance was influenced by rs4149057 in SLCO1B1, body surface area, and co-administration of 5-fluorouracil/leucovorin/bevacizumab. SN-38 clearance was influenced by rs887829 in UGT1A1, pre-treatment total bilirubin, and EGFR rare variant burden. Within each UGT1A1 genotype group, elevated pre-treatment total bilirubin and/or presence of at least one rare variant in EGFR resulted in significantly lower SN-38 clearance. The model reduced the interindividual variability in irinotecan clearance from 38 to 34% and SN-38 clearance from 49 to 32%.
This new model significantly reduced the interindividual variability in the clearance of irinotecan and SN-38. New genetic factors of variability in clearance have been identified.
Expression quantitative trait locus (eQTL) studies in human liver are crucial for elucidating how genetic variation influences variability in disease risk and therapeutic outcomes and may help guide ...strategies to obtain maximal efficacy and safety of clinical interventions. Associations between expression microarray and genome‐wide genotype data from four human liver eQTL studies (n = 1,183) were analyzed. More than 2.3 million cis‐eQTLs for 15,668 genes were identified. When eQTLs were filtered against a list of 1,496 drug response genes, 187,829 cis‐eQTLs for 1,191 genes were identified. Additionally, 1,683 sex‐biased cis‐eQTLs were identified, as well as 49 and 73 cis‐eQTLs that colocalized with genome‐wide association study signals for blood metabolite or lipid levels, respectively. Translational relevance of these results is evidenced by linking DPYD eQTLs to differences in safety of chemotherapy, linking the sex‐biased regulation of PCSK9 expression to anti‐lipid therapy, and identifying the G‐protein coupled receptor GPR180 as a novel drug target for hypertriglyceridemia.