Functional characterization of liver enhancers that regulate drug-associated transporters

Identification of drug responsive enhancers by predicting chromatin accessibility change from perturbed gene expression profiles

Individual may response to drug treatment differently due to their genetic variants located in enhancers. These variants can alter transcription factor's (TF) binding strength, affect enhancer's chromatin activity or interaction, and eventually change expression level of downstream gene. Here, we propose a computational framework, PERD, to Predict the Enhancers Responsive to Drug. A machine learning model was trained to predict the genome-wide chromatin accessibility from transcriptome data using the paired expression and chromatin accessibility data collected from ENCODE and ROADMAP. Then the model was applied to the perturbed gene expression data from Connectivity Map (CMAP) and Cancer Drug-induced gene expression Signature DataBase (CDS-DB) and identify drug responsive enhancers with significantly altered chromatin accessibility. Furthermore, the drug responsive enhancers were related to the pharmacogenomics genome-wide association studies (PGx GWAS). Stepping on the traditional drug-associated gene signatures, PERD holds the promise to enhance the causality of drug perturbation by providing candidate regulatory element of those drug associated genes.

Rare Variants in the ABCG2 Promoter Modulate In Vivo Activity

ABCG2 encodes the breast cancer resistance protein (BCRP), an efflux membrane transporter important in the detoxification of xenobiotics. In the present study, the basal activity of the ABCG2 promoter in liver, kidney, intestine, and breast cell lines was examined using luciferase reporter assays. The promoter activities of reference and variant ABCG2 sequences were compared in human hepatocellular carcinoma cell (HepG2), human embryonic kidney cell (HEK293T), human colorectal carcinoma cell (HCT116), and human breast adenocarcinoma cell (MCF-7) lines. The ABCG2 promoter activity was strongest in the kidney and intestine cell lines. Four variants in the basal ABCG2 promoter (rs76656413, rs66664036, rs139256004, and rs59370292) decreased the promoter activity by 25%-50% in at least three of the four cell lines. The activity of these four variants was also examined in vivo using the hydrodynamic tail vein assay, and two single nucleotide polymorphisms (rs76656413 and rs59370292) significantly decreased in vivo liver promoter activity by 50%-80%. Electrophoretic mobility shift assays confirmed a reduction in nuclear protein binding to the rs59370292 variant probe, whereas the rs76656413 probe had a shift in transcription factor binding specificity. Although both rs59370292 and rs76656413 are rare variants in all populations, they could contribute to patient-level variation in ABCG2 expression in the kidney, liver, and intestine.

ABCG2 regulatory single-nucleotide polymorphisms alter in vivo enhancer activity and expression

OBJECTIVES

The expression and activity of the breast cancer resistance protein (ABCG2) contributes toward the pharmacokinetics of endogenous and xenobiotic substrates. The effect of genetic variation on the activity of cis-regulatory elements and nuclear response elements in the ABCG2 locus and their contribution toward ABCG2 expression have not been investigated systematically. In this study, the effect of genetic variation on the in vitro and in vivo enhancer activity of six previously identified liver enhancers in the ABCG2 locus was examined.

METHODS

Reference and variant liver enhancers were tested for their ability to alter luciferase activity in vitro in HepG2 and HEK293T cell lines and in vivo using a hydrodynamic tail vein assay. Positive in vivo single-nucleotide polymorphisms (SNPs) were tested for association with gene expression and for altered protein binding in electrophoretic mobility shift assays.

RESULTS

Multiple SNPs were found to alter enhancer activity in vitro. Four of these variants (rs9999111, rs12508471, ABCG2RE1*2, and rs149713212) decreased and one (rs2725263) increased enhancer activity in vivo. In addition, rs9999111 and rs12508471 were associated with ABCG2 expression in lymphoblastoid cell lines, lymphocytes, and T cells, and showed increased HepG2 nuclear protein binding.

CONCLUSION

This study identifies SNPs within regulatory regions of the ABCG2 locus that alter enhancer activity in vitro and in vivo. Several of these SNPs correlate with tissue-specific ABCG2 expression and alter DNA/protein binding. These SNPs could contribute toward reported tissue-specific variability in ABCG2 expression and may influence the correlation between ABCG2 expression and disease risk or the pharmacokinetics and pharmacodynamics of breast cancer resistance protein substrates.

In Vivo Hepatic Enhancer Elements in the Human ABCG2 Locus

ABCG2 encodes the mitoxantrone resistance protein (MXR; breast cancer resistance protein), an ATP-binding cassette (ABC) efflux membrane transporter. Computational analysis of the ∼300 kb region of DNA surrounding ABCG2 (chr4:88911376-89220011, hg19) identified 30 regions with potential cis-regulatory capabilities. These putative regulatory regions were tested for their enhancer and suppressor activity in a human liver cell line using luciferase reporter assays. The in vitro enhancer and suppressor assays identified four regions that decreased gene expression and five regions that increased expression >1.6-fold. Four of five human hepatic in vitro enhancers were confirmed as in vivo liver enhancers using the mouse hydrodynamic tail vein injection assay. Two of the in vivo liver enhancers (ABCG2RE1 and ABCG2RE9) responded to 17β-estradiol or rifampin in human cell lines, and ABCG2RE9 had ChIP-seq evidence to support the binding of several transcription factors and the transcriptional coactivator p300 in human hepatocytes. This study identified genomic regions surrounding human ABCG2 that can function as regulatory elements, some with the capacity to alter gene expression upon environmental stimulus. The results from this research will drive future investigations of interindividual variation in ABCG2 expression and function that contribute to differences in drug response.

ABC transporter polymorphisms are associated with irinotecan pharmacokinetics and neutropenia

Neutropenia is a common dose-limiting toxicity associated with irinotecan treatment. Although UGT1A1 variants have been associated with neutropenia, a fraction of neutropenia risk remains unaccounted for. To identify additional genetic markers contributing to variability in irinotecan pharmacokinetics and neutropenia, a regression analysis was performed in 78 irinotecan-treated patients to analyze comprehensively three hepatic efflux transporter genes (ABCB1, ABCC1 and ABCG2). rs6498588 (ABCC1) and rs12720066 (ABCB1) were associated with increased SN-38 exposure, and rs17501331 (ABCC1) and rs12720066 were associated with lower absolute neutrophil count nadir. rs6498588 and a variant in high linkage disequilibrium are located in transcriptionally active regions or are predicted to alter transcription factor binding sites. While enhancer activity was not evident in vitro for genomic regions containing these single-nucleotide polymorphisms, rs6498588 was significantly associated with ABCC1 expression in human liver. These results suggest that genetic variation in ABCC1 and ABCB1 may contribute to irinotecan-induced neutropenia by altering expression of transporters involved in irinotecan metabolite disposition.

SUGP1 is a novel regulator of cholesterol metabolism

A large haplotype on chromosome 19p13.11 tagged by rs10401969 in intron 8 of SURP and G patch domain containing 1 (SUGP1) is associated with coronary artery disease (CAD), plasma LDL cholesterol levels, and other energy metabolism phenotypes. Recent studies have suggested that TM6SF2 is the causal gene within the locus, but we postulated that this locus could harbor additional CAD risk genes, including the putative splicing factor SUGP1. Indeed, we found that rs10401969 regulates SUGP1 exon 8 skipping, causing non-sense-mediated mRNA decay. Hepatic Sugp1 overexpression in CD1 male mice increased plasma cholesterol levels 20–50%. In human hepatoma cell lines, SUGP1 knockdown stimulated 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) alternative splicing and decreased HMGCR transcript stability, thus reducing cholesterol synthesis and increasing LDL uptake. Our results strongly support a role for SUGP1 as a novel regulator of cholesterol metabolism and suggest that it contributes to the relationship between rs10401969 and plasma cholesterol.

Multi-species, multi-transcription factor binding highlights conserved control of tissue-specific biological pathways

As exome sequencing gives way to genome sequencing, the need to interpret the function of regulatory DNA becomes increasingly important. To test whether evolutionary conservation of cis-regulatory modules (CRMs) gives insight into human gene regulation, we determined transcription factor (TF) binding locations of four liver-essential TFs in liver tissue from human, macaque, mouse, rat, and dog. Approximately, two thirds of the TF-bound regions fell into CRMs. Less than half of the human CRMs were found as a CRM in the orthologous region of a second species. Shared CRMs were associated with liver pathways and disease loci identified by genome-wide association studies. Recurrent rare human disease causing mutations at the promoters of several blood coagulation and lipid metabolism genes were also identified within CRMs shared in multiple species. This suggests that multi-species analyses of experimentally determined combinatorial TF binding will help identify genomic regions critical for tissue-specific gene control.

Pharmacogenetics of childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer in developed countries. Although treatment outcome has improved owing to advances in chemotherapy, there is still a group of patients for which therapy fails while some patients experience severe toxicity. In the last few years, several pharmacogenetic studies have been performed to search for markers of outcome and toxicity in pediatric ALL. However, to date, TPMT is the only pharmacogenetic marker in ALL with clinical guidelines for drug dosing. In this article, we will provide an overview of the most important findings carried out in pharmacogenetics for pediatric ALL, such as the interest drawn by methotrexate transporters in the context of methotrexate treatment. Even if most of the studies are centered on coding genes, we will also point to new approaches focusing on noncoding regions and epigenetic variation that could be interesting for consideration in the near future.

A novel test for selection on cis-regulatory elements reveals positive and negative selection acting on mammalian transcriptional enhancers

Measuring natural selection on genomic elements involved in the cis-regulation of gene expression--such as transcriptional enhancers and promoters--is critical for understanding the evolution of genomes, yet it remains a major challenge. Many studies have attempted to detect positive or negative selection in these noncoding elements by searching for those with the fastest or slowest rates of evolution, but this can be problematic. Here, we introduce a new approach to this issue, and demonstrate its utility on three mammalian transcriptional enhancers. Using results from saturation mutagenesis studies of these enhancers, we classified all possible point mutations as upregulating, downregulating, or silent, and determined which of these mutations have occurred on each branch of a phylogeny. Applying a framework analogous to Ka/Ks in protein-coding genes, we measured the strength of selection on upregulating and downregulating mutations, in specific branches as well as entire phylogenies. We discovered distinct modes of selection acting on different enhancers: although all three have experienced negative selection against downregulating mutations, the selection pressures on upregulating mutations vary. In one case, we detected positive selection for upregulation, whereas the other two had no detectable selection on upregulating mutations. Our methodology is applicable to the growing number of saturation mutagenesis data sets, and provides a detailed picture of the mode and strength of natural selection acting on cis-regulatory elements.

Massively parallel decoding of mammalian regulatory sequences supports a flexible organizational model

Despite continual progress in the cataloging of vertebrate regulatory elements, little is known about their organization and regulatory architecture. Here we describe a massively parallel experiment to systematically test the impact of copy number, spacing, combination and order of transcription factor binding sites on gene expression. A complex library of ∼5,000 synthetic regulatory elements containing patterns from 12 liver-specific transcription factor binding sites was assayed in mice and in HepG2 cells. We find that certain transcription factors act as direct drivers of gene expression in homotypic clusters of binding sites, independent of spacing between sites, whereas others function only synergistically. Heterotypic enhancers are stronger than their homotypic analogs and favor specific transcription factor binding site combinations, mimicking putative native enhancers. Exhaustive testing of binding site permutations suggests that there is flexibility in binding site order. Our findings provide quantitative support for a flexible model of regulatory element activity and suggest a framework for the design of synthetic tissue-specific enhancers.

GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer

Genome-wide association studies (GWAS) have identified four susceptibility loci for epithelial ovarian cancer (EOC), with another two suggestive loci reaching near genome-wide significance. We pooled data from a GWAS conducted in North America with another GWAS from the UK. We selected the top 24,551 SNPs for inclusion on the iCOGS custom genotyping array. We performed follow-up genotyping in 18,174 individuals with EOC (cases) and 26,134 controls from 43 studies from the Ovarian Cancer Association Consortium. We validated the two loci at 3q25 and 17q21 that were previously found to have associations close to genome-wide significance and identified three loci newly associated with risk: two loci associated with all EOC subtypes at 8q21 (rs11782652, P = 5.5 × 10(-9)) and 10p12 (rs1243180, P = 1.8 × 10(-8)) and another locus specific to the serous subtype at 17q12 (rs757210, P = 8.1 × 10(-10)). An integrated molecular analysis of genes and regulatory regions at these loci provided evidence for functional mechanisms underlying susceptibility and implicated CHMP4C in the pathogenesis of ovarian cancer.

Drug-induced liver injury: the role of drug metabolism and transport

Many studies have pinpointed the significant contribution of liver-mediated drug metabolism and transport to the complexity of drug-induced liver injury (DILI). Phase I cytochrome P450 (CYP450) enzymes can lead to altered drug metabolism and formation of toxic metabolites, whilst Phase II enzymes are also associated with DILI. The emerging role of hepatic transporters in regulating the movement of endogenous and exogenous chemicals (e.g., bile acids and drugs) across cellular and tissue membranes is critical in determining the pathophysiology of liver disease as well as drug toxicity and efficacy. Genetic and environmental factors can have a significant impact on drug metabolism and transporter proteins, consequently increasing the risk of DILI in susceptible individuals. The assessment of these factors therefore represents an important approach for predicting and preventing DILI, by better understanding the pharmacological profile of a specific drug. This review focuses on the mechanisms of DILI associated with drug metabolism and hepatic transport, and how they can be influenced by underlying factors.

Sequence signatures extracted from proximal promoters can be used to predict distal enhancers

BACKGROUND

Gene expression is controlled by proximal promoters and distal regulatory elements such as enhancers. While the activity of some promoters can be invariant across tissues, enhancers tend to be highly tissue-specific.

RESULTS

We compiled sets of tissue-specific promoters based on gene expression profiles of 79 human tissues and cell types. Putative transcription factor binding sites within each set of sequences were used to train a support vector machine classifier capable of distinguishing tissue-specific promoters from control sequences. We obtained reliable classifiers for 92% of the tissues, with an area under the receiver operating characteristic curve between 60% (for subthalamic nucleus promoters) and 98% (for heart promoters). We next used these classifiers to identify tissue-specific enhancers, scanning distal non-coding sequences in the loci of the 200 most highly and lowly expressed genes. Thirty percent of reliable classifiers produced consistent enhancer predictions, with significantly higher densities in the loci of the most highly expressed compared to lowly expressed genes. Liver enhancer predictions were assessed in vivo using the hydrodynamic tail vein injection assay. Fifty-eight percent of the predictions yielded significant enhancer activity in the mouse liver, whereas a control set of five sequences was completely negative.

CONCLUSIONS

We conclude that promoters of tissue-specific genes often contain unambiguous tissue-specific signatures that can be learned and used for the de novo prediction of enhancers.

The hydrodynamic tail vein assay as a tool for the study of liver promoters and enhancers

The hydrodynamic tail vein injection is a technique that is used to deliver nucleic acids into live mice. Delivery through this method results in the in vivo transfection of foreign DNA primarily in the liver. Here, we describe the use of this technique to test for regulatory activity of liver promoters and enhancers, using a dual luciferase reporter system as the readable/measureable output and how this application can be used for pharmacogenomic studies.

Genome-wide study of methotrexate clearance replicates SLCO1B1

Methotrexate clearance can influence the cure of and toxicity in children with acute lymphoblastic leukemia (ALL). We estimated methotrexate plasma clearance for 1279 patients with ALL treated with methotrexate (24-hour infusion of a 1 g/m2 dose or 4-hour infusion of a 2 g/m2 dose) on the Children’s Oncology Group P9904 and P9905 protocols. Methotrexate clearance was lower in older children (P = 7 x 10(-7)), girls (P = 2.7 x 10(-4)), and those who received a delayed-intensification phase (P = .0022). A genome-wide analysis showed that methotrexate clearance was associated with polymorphisms in the organic anion transporter gene SLCO1B1 (P = 2.1 x 10(-11)). This replicates findings using different schedules of high-dose methotrexate in St Jude ALL treatment protocols; a combined meta-analysis yields a P value of 5.7 x 10(-19) for the association of methotrexate clearance with SLCO1B1 SNP rs4149056. Validation of this variant with 5 different treatment regimens of methotrexate solidifies the robustness of this pharmacogenomic determinant of methotrexate clearance. This study is registered at http://www.clinicaltrials.gov as NCT00005585 and NCT00005596.

Pharmacogene regulatory elements: from discovery to applications

Regulatory elements play an important role in the variability of individual responses to drug treatment. This has been established through studies on three classes of elements that regulate RNA and protein abundance: promoters, enhancers and microRNAs. Each of these elements, and genetic variants within them, are being characterized at an exponential pace by next-generation sequencing (NGS) technologies. In this review, we outline examples of how each class of element affects drug response via regulation of drug targets, transporters and enzymes. We also discuss the impact of NGS technologies such as chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq), and the ramifications of new techniques such as high-throughput chromosome capture (Hi-C), chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and massively parallel reporter assays (MPRA). NGS approaches are generating data faster than they can be analyzed, and new methods will be required to prioritize laboratory results before they are ready for the clinic. However, there is no doubt that these approaches will bring about a systems-level understanding of the interplay between genetic variants and drug response. An understanding of the importance of regulatory variants in pharmacogenomics will facilitate the identification of responders versus non-responders, the prevention of adverse effects and the optimization of therapies for individual patients.

Massively parallel functional dissection of mammalian enhancers in vivo

The functional consequences of genetic variation in mammalian regulatory elements are poorly understood. We report the in vivo dissection of three mammalian enhancers at single-nucleotide resolution through a massively parallel reporter assay. For each enhancer, we synthesized a library of >100,000 mutant haplotypes with 2-3% divergence from the wild-type sequence. Each haplotype was linked to a unique sequence tag embedded within a transcriptional cassette. We introduced each enhancer library into mouse liver and measured the relative activities of individual haplotypes en masse by sequencing the transcribed tags. Linear regression analysis yielded highly reproducible estimates of the effect of every possible single-nucleotide change on enhancer activity. The functional consequence of most mutations was modest, with ∼22% affecting activity by >1.2-fold and ∼3% by >2-fold. Several, but not all, positions with higher effects showed evidence for purifying selection, or co-localized with known liver-associated transcription factor binding sites, demonstrating the value of empirical high-resolution functional analysis.