The Transcription Factor Bhlhe40 Programs Mitochondrial Regulation of Resident CD8+ T Cell Fitness and Functionality.

Tissue-resident memory CD8+ T (Trm) cells share core residency gene programs with tumor-infiltrating lymphocytes (TILs). However, the transcriptional, metabolic, and epigenetic regulation of Trm cell and TIL development and function is largely undefined. Here, we found that the transcription factor Bhlhe40 was specifically required for Trm cell and TIL development and polyfunctionality. Local PD-1 signaling inhibited TIL Bhlhe40 expression, and Bhlhe40 was critical for TIL reinvigoration following anti-PD-L1 blockade. Mechanistically, Bhlhe40 sustained Trm cell and TIL mitochondrial fitness and a functional epigenetic state. Building on these findings, we identified an epigenetic and metabolic regimen that promoted Trm cell and TIL gene signatures associated with tissue residency and polyfunctionality. This regimen empowered the anti-tumor activity of CD8+ T cells and possessed therapeutic potential even at an advanced tumor stage in mouse models. Our results provide mechanistic insights into the local regulation of Trm cell and TIL function.

Phosphorylated RB Promotes Cancer Immunity by Inhibiting NF-κB Activation and PD-L1 Expression.

Aberrant expression of programmed death ligand-1 (PD-L1) in tumor cells promotes cancer progression by suppressing cancer immunity. The retinoblastoma protein RB is a tumor suppressor known to regulate the cell cycle, DNA damage response, and differentiation. Here, we demonstrate that RB interacts with nuclear factor κB (NF-κB) protein p65 and that their interaction is primarily dependent on CDK4/6-mediated serine-249/threonine-252 (S249/T252) phosphorylation of RB. RNA-seq analysis shows a subset of NF-κB pathway genes including PD-L1 are selectively upregulated by RB knockdown or CDK4/6 inhibitor. S249/T252-phosphorylated RB inversely correlates with PD-L1 expression in patient samples. Expression of a RB-derived S249/T252 phosphorylation-mimetic peptide suppresses radiotherapy-induced upregulation of PD-L1 and augments therapeutic efficacy of radiation in vivo. Our findings reveal a previously unrecognized tumor suppressor function of hyperphosphorylated RB in suppressing NF-κB activity and PD-L1 expression and suggest that the RB-NF-κB axis can be exploited to overcome cancer immune evasion triggered by conventional or targeted therapies.

DUB3 Promotes BET Inhibitor Resistance and Cancer Progression by Deubiquitinating BRD4.

The bromodomain and extra-terminal domain (BET) protein BRD4 is emerging as a promising anticancer therapeutic target. However, resistance to BET inhibitors often occurs, and it has been linked to aberrant degradation of BRD4 protein in cancer. Here, we demonstrate that the deubiquitinase DUB3 binds to BRD4 and promotes its deubiquitination and stabilization. Expression of DUB3 is transcriptionally repressed by the NCOR2-HDAC10 complex. The NCOR2 gene is frequently deleted in castration-resistant prostate cancer patient specimens, and loss of NCOR2 induces elevation of DUB3 and BRD4 proteins in cancer cells. DUB3-proficient prostate cancer cells are resistant to the BET inhibitor JQ1 in vitro and in mice, but this effect is diminished by DUB3 inhibitory agents such as CDK4/6 inhibitor in a RB-independent manner. Our findings identify a previously unrecognized mechanism causing BRD4 upregulation and drug resistance, suggesting that DUB3 is a viable therapeutic target to overcome BET inhibitor resistance in cancer.

Targeting Src SH3 domain-mediated glycolysis of HSC suppresses transcriptome, myofibroblastic activation, and colorectal liver metastasis.

BACKGROUND AND AIMS: Transforming growth factor-beta 1 (TGFß1) induces HSC activation into metastasis-promoting cancer-associated fibroblasts (CAFs), but how the process is fueled remains incompletely understood. We studied metabolic reprogramming induced by TGFß1 in HSCs. APPROACHES AND RESULTS: Activation of cultured primary human HSCs was assessed by the expression of myofibroblast markers. Glucose transporter 1 (Glut1) of murine HSC was disrupted by Cre recombinase/LoxP sequence derived from bacteriophage P1 recombination (Cre/LoxP). Plasma membrane (PM) Glut1 and glycolysis were studied by biotinylation assay and the Angilent Seahorse XFe96 Analyzer. S.c. HSC/tumor co-implantation and portal vein injection of MC38 colorectal cancer cells into HSC-specific Glut1 knockout mice were performed to determine in vivo relevance. Transcriptome was obtained by RNA sequencing of HSCs and spatialomics with MC38 liver metastases. TGFß1-induced CAF activation of HSCs was accompanied by elevation of PM Glut1, glucose uptake, and glycolysis. Targeting Glut1 or Src by short hairpin RNA, pharmacologic inhibition, or a Src SH3 domain deletion mutant abrogated TGFß1-stimulated PM accumulation of Glut1, glycolysis, and CAF activation. Mechanistically, binding of the Src SH3 domain to SH3 domain-binding protein 5 led to a Src/SH3 domain-binding protein 5/Rab11/Glut1 complex that activated Rab11-dependent Glut1 PM transport under TGFß1 stimulation. Deleting the Src SH3 domain or targeting Glut1 of HSCs by short hairpin RNA or Cre recombinase/LoxP sequence derived from bacteriophage P1 recombination suppressed CAF activation in mice and MC38 colorectal liver metastasis. Multi-omics revealed that Glut1 deficiency in HSCs/CAFs suppressed HSC expression of tumor-promoting factors and altered MC38 transcriptome, contributing to reduced MC38 liver metastases. CONCLUSION: The Src SH3 domain-facilitated metabolic reprogramming induced by TGFß1 represents a target to inhibit CAF activation and the pro-metastatic liver microenvironment.

Loss of BRD4 induces cell senescence in HSC/HPCs by deregulating histone H3 clipping.

Bromodomain-containing protein 4 (BRD4) is overexpressed and functionally implicated in various myeloid malignancies. However, the role of BRD4 in normal hematopoiesis remains largely unknown. Here, utilizing an inducible Brd4 knockout mouse model, we find that deletion of Brd4 (Brd4Δ/Δ ) in the hematopoietic system impairs hematopoietic stem cell (HSC) self-renewal and differentiation, which associates with cell cycle arrest and senescence. ATAC-seq analysis shows increased chromatin accessibility in Brd4Δ/Δ hematopoietic stem/progenitor cells (HSC/HPCs). Genome-wide mapping with cleavage under target and release using nuclease (CUT&RUN) assays demonstrate that increased global enrichment of H3K122ac and H3K4me3 in Brd4Δ/Δ HSC/HPCs is associated with the upregulation of senescence-specific genes. Interestingly, Brd4 deletion increases clipped H3 (cH3) which correlates with the upregulation of senescence-specific genes and results in a higher frequency of senescent HSC/HPCs. Re-expression of BRD4 reduces cH3 levels and rescues the senescence rate in Brd4Δ/Δ HSC/HPCs. This study unveils an important role of BRD4 in HSC/HPC function by preventing H3 clipping and suppressing senescence gene expression.

Hypoxia-Inducible Factor 1α Stabilization Restores Epigenetic Control of Nitric Oxide Synthase 1 Expression and Reverses Gastroparesis in Female Diabetic Mice.

BACKGROUND & AIMS: Although depletion of neuronal nitric oxide synthase (NOS1)-expressing neurons contributes to gastroparesis, stimulating nitrergic signaling is not an effective therapy. We investigated whether hypoxia-inducible factor 1α (HIF1A), which is activated by high O2 consumption in central neurons, is a Nos1 transcription factor in enteric neurons and whether stabilizing HIF1A reverses gastroparesis. METHODS: Mice with streptozotocin-induced diabetes, human and mouse tissues, NOS1+ mouse neuroblastoma cells, and isolated nitrergic neurons were studied. Gastric emptying of solids and volumes were determined by breath test and single-photon emission computed tomography, respectively. Gene expression was analyzed by RNA-sequencing, microarrays, immunoblotting, and immunofluorescence. Epigenetic assays included chromatin immunoprecipitation sequencing (13 targets), chromosome conformation capture sequencing, and reporter assays. Mechanistic studies used Cre-mediated recombination, RNA interference, and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated epigenome editing. RESULTS: HIF1A signaling from physiological intracellular hypoxia was active in mouse and human NOS1+ myenteric neurons but reduced in diabetes. Deleting Hif1a in Nos1-expressing neurons reduced NOS1 protein by 50% to 92% and delayed gastric emptying of solids in female but not male mice. Stabilizing HIF1A with roxadustat (FG-4592), which is approved for human use, restored NOS1 and reversed gastroparesis in female diabetic mice. In nitrergic neurons, HIF1A up-regulated Nos1 transcription by binding and activating proximal and distal cis-regulatory elements, including newly discovered super-enhancers, facilitating RNA polymerase loading and pause-release, and by recruiting cohesin to loop anchors to alter chromosome topology. CONCLUSIONS: Pharmacologic HIF1A stabilization is a novel, translatable approach to restoring nitrergic signaling and treating diabetic gastroparesis. The newly recognized effects of HIF1A on chromosome topology may provide insights into physioxia- and ischemia-related organ function.

Deep learning tackles single-cell analysis-a survey of deep learning for scRNA-seq analysis.

Since its selection as the method of the year in 2013, single-cell technologies have become mature enough to provide answers to complex research questions. With the growth of single-cell profiling technologies, there has also been a significant increase in data collected from single-cell profilings, resulting in computational challenges to process these massive and complicated datasets. To address these challenges, deep learning (DL) is positioned as a competitive alternative for single-cell analyses besides the traditional machine learning approaches. Here, we survey a total of 25 DL algorithms and their applicability for a specific step in the single cell RNA-seq processing pipeline. Specifically, we establish a unified mathematical representation of variational autoencoder, autoencoder, generative adversarial network and supervised DL models, compare the training strategies and loss functions for these models, and relate the loss functions of these models to specific objectives of the data processing step. Such a presentation will allow readers to choose suitable algorithms for their particular objective at each step in the pipeline. We envision that this survey will serve as an important information portal for learning the application of DL for scRNA-seq analysis and inspire innovative uses of DL to address a broader range of new challenges in emerging multi-omics and spatial single-cell sequencing.

Phosphorylated MED1 links transcription recycling and cancer growth.

Mediator activates RNA polymerase II (Pol II) function during transcription, but it remains unclear whether Mediator is able to travel with Pol II and regulate Pol II transcription beyond the initiation and early elongation steps. By using in vitro and in vivo transcription recycling assays, we find that human Mediator 1 (MED1), when phosphorylated at the mammal-specific threonine 1032 by cyclin-dependent kinase 9 (CDK9), dynamically moves along with Pol II throughout the transcribed genes to drive Pol II recycling after the initial round of transcription. Mechanistically, MED31 mediates the recycling of phosphorylated MED1 and Pol II, enhancing mRNA output during the transcription recycling process. Importantly, MED1 phosphorylation increases during prostate cancer progression to the lethal phase, and pharmacological inhibition of CDK9 decreases prostate tumor growth by decreasing MED1 phosphorylation and Pol II recycling. Our results reveal a novel role of MED1 in Pol II transcription and identify phosphorylated MED1 as a targetable driver of dysregulated Pol II recycling in cancer.

Glucocorticoids unmask silent non-coding genetic risk variants for common diseases.

Understanding the function of non-coding genomic sequence variants represents a challenge for biomedicine. Many diseases are products of gene-by-environment interactions with complex mechanisms. This study addresses these themes by mechanistic characterization of non-coding variants that influence gene expression only after drug or hormone exposure. Using glucocorticoid signaling as a model system, we integrated genomic, transcriptomic, and epigenomic approaches to unravel mechanisms by which variant function could be revealed by hormones or drugs. Specifically, we identified cis-regulatory elements and 3D interactions underlying ligand-dependent associations between variants and gene expression. One-quarter of the glucocorticoid-modulated variants that we identified had already been associated with clinical phenotypes. However, their affected genes were 'unmasked' only after glucocorticoid exposure and often with function relevant to the disease phenotypes. These diseases involved glucocorticoids as risk factors or therapeutic agents and included autoimmunity, metabolic and mood disorders, osteoporosis and cancer. For example, we identified a novel breast cancer risk gene, MAST4, with expression that was repressed by glucocorticoids in cells carrying the risk genotype, repression that correlated with MAST4 expression in breast cancer and treatment outcomes. These observations provide a mechanistic framework for understanding non-coding genetic variant-chemical environment interactions and their role in disease risk and drug response.

EZH2 cooperates with gain-of-function p53 mutants to promote cancer growth and metastasis.

In light of the increasing number of identified cancer-driven gain-of-function (GOF) mutants of p53, it is important to define a common mechanism to systematically target several mutants, rather than developing strategies tailored to inhibit each mutant individually. Here, using RNA immunoprecipitation-sequencing (RIP-seq), we identified the Polycomb-group histone methyltransferase EZH2 as a p53 mRNA-binding protein. EZH2 bound to an internal ribosome entry site (IRES) in the 5'UTR of p53 mRNA and enhanced p53 protein translation in a methyltransferase-independent manner. EZH2 augmented p53 GOF mutant-mediated cancer growth and metastasis by increasing protein levels of mutant p53. EZH2 overexpression was associated with worsened outcome selectively in patients with p53-mutated cancer. Depletion of EZH2 by antisense oligonucleotides inhibited p53 GOF mutant-mediated cancer growth. Our findings reveal a non-methyltransferase function of EZH2 that controls protein translation of p53 GOF mutants, inhibition of which causes synthetic lethality in cancer cells expressing p53 GOF mutants.

CpGtools: a python package for DNA methylation analysis.

MOTIVATION: DNA methylation can be measured at the single CpG level using sodium bisulfite conversion of genomic DNA followed by sequencing or array hybridization. Many analytic tools have been developed, yet there is still a high demand for a comprehensive and multifaceted tool suite to analyze, annotate, QC and visualize the DNA methylation data. RESULTS: We developed the CpGtools package to analyze DNA methylation data generated from bisulfite sequencing or Illumina methylation arrays. The CpGtools package consists of three types of modules: (i) 'CpG position modules' focus on analyzing the genomic positions of CpGs, including associating other genomic and epigenomic features to a given list of CpGs and generating the DNA motif logo enriched in the genomic contexts of a given list of CpGs; (ii) 'CpG signal modules' are designed to analyze DNA methylation values, such as performing the PCA or t-SNE analyses, using Bayesian Gaussian mixture modeling to classify CpG sites into fully methylated, partially methylated and unmethylated groups, profiling the average DNA methylation level over user-specified genomics regions and generating the bean/violin plots and (iii) 'differential CpG analysis modules' focus on identifying differentially methylated CpGs between groups using different statistical methods including Fisher's Exact Test, Student's t-test, ANOVA, non-parametric tests, linear regression, logistic regression, beta-binomial regression and Bayesian estimation. AVAILABILITY AND IMPLEMENTATION: CpGtools is written in Python under the open-source GPL license. The source code and documentation are freely available at https://github.com/liguowang/cpgtools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

TCF7L2 lncRNA: a link between bipolar disorder and body mass index through glucocorticoid signaling.

Bipolar disorder (BD) and obesity are highly comorbid. We previously performed a genome-wide association study (GWAS) for BD risk accounting for the effect of body mass index (BMI), which identified a genome-wide significant single-nucleotide polymorphism (SNP) in the gene encoding the transcription factor 7 like 2 (TCF7L2). However, the molecular function of TCF7L2 in the central nervous system (CNS) and its possible role in the BD and BMI interaction remained unclear. In the present study, we demonstrated by studying human induced pluripotent stem cell (hiPSC)-derived astrocytes, cells that highly express TCF7L2 in the CNS, that the BD-BMI GWAS risk SNP is associated with glucocorticoid-dependent repression of the expression of a previously uncharacterized TCF7L2 transcript variant. That transcript is a long non-coding RNA (lncRNA-TCF7L2) that is highly expressed in the CNS but not in peripheral tissues such as the liver and pancreas that are involved in metabolism. In astrocytes, knockdown of the lncRNA-TCF7L2 resulted in decreased expression of the parent gene, TCF7L2, as well as alterations in the expression of a series of genes involved in insulin signaling and diabetes. We also studied the function of TCF7L2 in hiPSC-derived astrocytes by integrating RNA sequencing data after TCF7L2 knockdown with TCF7L2 chromatin-immunoprecipitation sequencing (ChIP-seq) data. Those studies showed that TCF7L2 directly regulated a series of BD risk genes. In summary, these results support the existence of a CNS-based mechanism underlying BD-BMI genetic risk, a mechanism based on a glucocorticoid-dependent expression quantitative trait locus that regulates the expression of a novel TCF7L2 non-coding transcript.

BCL6 modulates tissue neutrophil survival and exacerbates pulmonary inflammation following influenza virus infection.

Neutrophils are vital for antimicrobial defense; however, their role during viral infection is less clear. Furthermore, the molecular regulation of neutrophil fate and function at the viral infected sites is largely elusive. Here we report that BCL6 deficiency in myeloid cells exhibited drastically enhanced host resistance to severe influenza A virus (IAV) infection. In contrast to the notion that BCL6 functions to suppress innate inflammation, we find that myeloid BCL6 deficiency diminished lung inflammation without affecting viral loads. Using a series of Cre-transgenic, reporter, and knockout mouse lines, we demonstrate that BCL6 deficiency in neutrophils, but not in monocytes or lung macrophages, attenuated host inflammation and morbidity following IAV infection. Mechanistically, BCL6 bound to the neutrophil gene loci involved in cellular apoptosis in cells specifically at the site of infection. As such, BCL6 disruption resulted in increased expression of apoptotic genes in neutrophils in the respiratory tract, but not in the circulation or bone marrow. Consequently, BCL6 deficiency promoted tissue neutrophil apoptosis. Partial neutrophil depletion led to diminished pulmonary inflammation and decreased host morbidity. Our results reveal a previously unappreciated role of BCL6 in modulating neutrophil apoptosis at the site of infection for the regulation of host disease development following viral infection. Furthermore, our studies indicate that tissue-specific regulation of neutrophil survival modulates host inflammation and tissue immunopathology during acute respiratory viral infection.

TGFß-induced fibroblast activation requires persistent and targeted HDAC-mediated gene repression.

Tissue fibrosis is a chronic disease driven by persistent fibroblast activation that has recently been linked to epigenetic modifications. Here, we screened a small library of epigenetic small-molecule modulators to identify compounds capable of inhibiting or reversing TGFß-mediated fibroblast activation. We identified pracinostat, an HDAC inhibitor, as a potent attenuator of lung fibroblast activation and confirmed its efficacy in patient-derived fibroblasts isolated from fibrotic lung tissue. Mechanistically, we found that HDAC-dependent transcriptional repression was an early and essential event in TGFß-mediated fibroblast activation. Treatment of lung fibroblasts with pracinostat broadly attenuated TGFß-mediated epigenetic repression and promoted fibroblast quiescence. We confirmed a specific role for HDAC-dependent histone deacetylation in the promoter region of the anti-fibrotic gene PPARGC1A (PGC1α) in response to TGFß stimulation. Finally, we identified HDAC7 as a key factor whose siRNA-mediated knockdown attenuates fibroblast activation without altering global histone acetylation. Together, these results provide novel mechanistic insight into the essential role HDACs play in TGFß-mediated fibroblast activation via targeted gene repression.

Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

Androgen receptor splice variants bind to constitutively open chromatin and promote abiraterone-resistant growth of prostate cancer.

Androgen receptor (AR) splice variants (ARVs) are implicated in development of castration-resistant prostate cancer (CRPC). Upregulation of ARVs often correlates with persistent AR activity after androgen deprivation therapy (ADT). However, the genomic and epigenomic characteristics of ARV-dependent cistrome and the disease relevance of ARV-mediated transcriptome remain elusive. Through integrated chromatin immunoprecipitation coupled sequencing (ChIP-seq) and RNA sequencing (RNA-seq) analysis, we identified ARV-preferential-binding sites (ARV-PBS) and a set of genes preferentially transactivated by ARVs in CRPC cells. ARVs preferentially bind to enhancers located in nucleosome-depleted regions harboring the full AR-response element (AREfull), while full-length AR (ARFL)-PBS are enhancers resided in closed chromatin regions containing the composite FOXA1-nnnn-AREhalf motif. ARV-PBS exclusively overlapped with AR binding sites in castration-resistant (CR) tumors in patients and ARV-preferentially activated genes were up-regulated in abiraterone-resistant patient specimens. Expression of ARV-PBS target genes, such as oncogene RAP2A and cell cycle gene E2F7, were significantly associated with castration resistance, poor survival and tumor progression. We uncover distinct genomic and epigenomic features of ARV-PBS, highlighting that ARVs are useful tools to depict AR-regulated oncogenic genome and epigenome landscapes in prostate cancer. Our data also suggest that the ARV-preferentially activated transcriptional program could be targeted for effective treatment of CRPC.

Enhanced and controlled chromatin extraction from FFPE tissues and the application to ChIP-seq.

BACKGROUND: Epigenetic dysregulation is involved in the etiology and progression of various human diseases. Formalin-fixed paraffin-embedded (FFPE) samples represent the gold standard for archiving pathology samples, and thus FFPE samples are a major resource of samples in clinical research. However, chromatin-based epigenetic assays in the clinical settings are limited to fresh or frozen samples, and are hampered by low chromatin yield in FFPE samples due to the lack of a reliable and efficient chromatin preparation method. Here, we introduce a new chromatin extraction method from FFPE tissues (Chrom-EX PE) for chromatin-based epigenetic assays. RESULTS: During rehydration of FFPE tissues, applying a tissue-level cross-link reversal into the deparaffinized tissue at 65 °C dramatically increased chromatin yield in the soluble fraction. The resulting chromatin is compatible with targeted ChIP-qPCR and genome-wide ChIP-seq approaches. The chromatin prepared by Chrom-EX PE showed a gradual fragmentation pattern with varying incubation temperature. At temperatures below 37 °C, the majority of soluble chromatin is over 1 kb. The soluble chromatin prepared in the range of 45-60 °C showed a typical nucleosomal pattern. And the majority of chromatin prepared at 65 °C is close to mononucleosomal size. These observations indicate that chromatin preparation from FFPE samples can be controlled for downstream chromatin-based epigenetic assays. CONCLUSIONS: This study provided a new method that achieves efficient extraction of high-quality chromatin suitable for chromatin-based epigenetic assays with less damage on chromatin. This approach may provide a way to circumvent the over-fixed nature of FFPE tissues for future technology development.

Agonist and antagonist switch DNA motifs recognized by human androgen receptor in prostate cancer.

Human transcription factors recognize specific DNA sequence motifs to regulate transcription. It is unknown whether a single transcription factor is able to bind to distinctly different motifs on chromatin, and if so, what determines the usage of specific motifs. By using a motif-resolution chromatin immunoprecipitation-exonuclease (ChIP-exo) approach, we find that agonist-liganded human androgen receptor (AR) and antagonist-liganded AR bind to two distinctly different motifs, leading to distinct transcriptional outcomes in prostate cancer cells. Further analysis on clinical prostate tissues reveals that the binding of AR to these two distinct motifs is involved in prostate carcinogenesis. Together, these results suggest that unique ligands may switch DNA motifs recognized by ligand-dependent transcription factors in vivo. Our findings also provide a broad mechanistic foundation for understanding ligand-specific induction of gene expression profiles.

Single Nucleotide Polymorphisms at a Distance from Aryl Hydrocarbon Receptor (AHR) Binding Sites Influence AHR Ligand-Dependent Gene Expression.

Greater than 90% of significant genome-wide association study (GWAS) single-nucleotide polymorphisms (SNPs) are in noncoding regions of the genome, but only 25.6% are known expression quantitative trait loci (eQTLs). Therefore, the function of many significant GWAS SNPs remains unclear. We have identified a novel type of eQTL for which SNPs distant from ligand-activated transcription factor (TF) binding sites can alter target gene expression in a SNP genotype-by-ligand-dependent fashion that we refer to as pharmacogenomic eQTLs (PGx-eQTLs)-loci that may have important pharmacotherapeutic implications. In the present study, we integrated chromatin immunoprecipitation-seq with RNA-seq and SNP genotype data for a panel of lymphoblastoid cell lines to identify 10 novel cis PGx-eQTLs dependent on the ligand-activated TF aryl hydrocarbon receptor (AHR)-a critical environmental sensor for xenobiotic (drug) and immune response. Those 10 cis PGx-eQTLs were eQTLs only after AHR ligand treatment, even though the SNPs did not create/destroy an AHR response element-the DNA sequence motif recognized and bound by AHR. Additional functional studies in multiple cell lines demonstrated that some cis PGx-eQTLs are functional in multiple cell types, whereas others displayed SNP-by-ligand-dependent effects in just one cell type. Furthermore, four of those cis PGx-eQTLs had previously been associated with clinical phenotypes, indicating that those loci might have the potential to inform clinical decisions. Therefore, SNPs across the genome that are distant from TF binding sites for ligand-activated TFs might function as PGx-eQTLs and, as a result, might have important clinical implications for interindividual variation in drug response. SIGNIFICANCE STATEMENT: More than 90% of single-nucleotide polymorphisms (SNPs) that are associated with clinical phenotypes are located in noncoding regions of the genome. However, the mechanisms of action of many of those SNPs have not been elucidated, and drugs may unmask functional expression quantitative trail loci (eQTLs). In the current study, we used drugs that bind to the ligand-activated transcription factor aryl hydrocarbon receptor (AHR) and identified SNPs that were associated with interindividual variation in gene expression following drug exposure-termed pharmacogenomic (PGx)-eQTLs. Possibly of greater significance, those PGx-eQTL SNPs were outside of AHR binding sites, indicating that they do not interrupt AHR DNA recognition. PGx-eQTLs such as those described in this work may have crucial implications for interindividual variation in drug.