Tumour DDR1 promotes collagen fibre alignment to instigate immune exclusion.

Immune exclusion predicts poor patient outcomes in multiple malignancies, including triple-negative breast cancer (TNBC)1. The extracellular matrix (ECM) contributes to immune exclusion2. However, strategies to reduce ECM abundance are largely ineffective or generate undesired outcomes3,4. Here we show that discoidin domain receptor 1 (DDR1), a collagen receptor with tyrosine kinase activity5, instigates immune exclusion by promoting collagen fibre alignment. Ablation of Ddr1 in tumours promotes the intratumoral penetration of T cells and obliterates tumour growth in mouse models of TNBC. Supporting this finding, in human TNBC the expression of DDR1 negatively correlates with the intratumoral abundance of anti-tumour T cells. The DDR1 extracellular domain (DDR1-ECD), but not its intracellular kinase domain, is required for immune exclusion. Membrane-untethered DDR1-ECD is sufficient to rescue the growth of Ddr1-knockout tumours in immunocompetent hosts. Mechanistically, the binding of DDR1-ECD to collagen enforces aligned collagen fibres and obstructs immune infiltration. ECD-neutralizing antibodies disrupt collagen fibre alignment, mitigate immune exclusion and inhibit tumour growth in immunocompetent hosts. Together, our findings identify a mechanism for immune exclusion and suggest an immunotherapeutic target for increasing immune accessibility through reconfiguration of the tumour ECM.

Phosphorylation of EZH2 by AMPK Suppresses PRC2 Methyltransferase Activity and Oncogenic Function.

Sustained energy starvation leads to activation of AMP-activated protein kinase (AMPK), which coordinates energy status with numerous cellular processes including metabolism, protein synthesis, and autophagy. Here, we report that AMPK phosphorylates the histone methyltransferase EZH2 at T311 to disrupt the interaction between EZH2 and SUZ12, another core component of the polycomb repressive complex 2 (PRC2), leading to attenuated PRC2-dependent methylation of histone H3 at Lys27. As such, PRC2 target genes, many of which are known tumor suppressors, were upregulated upon T311-EZH2 phosphorylation, which suppressed tumor cell growth both in cell culture and mouse xenografts. Pathologically, immunohistochemical analyses uncovered a positive correlation between AMPK activity and pT311-EZH2, and higher pT311-EZH2 correlates with better survival in both ovarian and breast cancer patients. Our finding suggests that AMPK agonists might be promising sensitizers for EZH2-targeting cancer therapies.

An essential signaling function of cytoplasmic NELFB is independent of RNA polymerase II pausing.

The four-subunit negative elongation factor (NELF) complex mediates RNA polymerase II (Pol II) pausing at promoter-proximal regions. Ablation of individual NELF subunits destabilizes the NELF complex and causes cell lethality, leading to the prevailing concept that NELF-mediated Pol II pausing is essential for cell proliferation. Using separation-of-function mutations, we show here that NELFB function in cell proliferation can be uncoupled from that in Pol II pausing. NELFB mutants sequestered in the cytoplasm and deprived of NELF nuclear function still support cell proliferation and part of the NELFB-dependent transcriptome. Mechanistically, cytoplasmic NELFB physically and functionally interacts with prosurvival signaling kinases, most notably phosphatidylinositol-3-kinase/AKT. Ectopic expression of membrane-tethered phosphatidylinositol-3-kinase/AKT partially bypasses the role of NELFB in cell proliferation, but not Pol II occupancy. Together, these data expand the current understanding of the physiological impact of Pol II pausing and underscore the multiplicity of the biological functions of individual NELF subunits.

Cas13d knockdown of lung protease Ctsl prevents and treats SARS-CoV-2 infection.

SARS-CoV-2 entry into cells requires specific host proteases; however, no successful in vivo applications of host protease inhibitors have yet been reported for treatment of SARS-CoV-2 pathogenesis. Here we describe a chemically engineered nanosystem encapsulating CRISPR-Cas13d, developed to specifically target lung protease cathepsin L (Ctsl) messenger RNA to block SARS-CoV-2 infection in mice. We show that this nanosystem decreases lung Ctsl expression in normal mice efficiently, specifically and safely. We further show that this approach extends survival of mice lethally infected with SARS-CoV-2, correlating with decreased lung virus burden, reduced expression of proinflammatory cytokines/chemokines and diminished severity of pulmonary interstitial inflammation. Postinfection treatment by this nanosystem dramatically lowers the lung virus burden and alleviates virus-induced pathological changes. Our results indicate that targeting lung protease mRNA by Cas13d nanosystem represents a unique strategy for controlling SARS-CoV-2 infection and demonstrate that CRISPR can be used as a potential treatment for SARS-CoV-2 infection.

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.

Autistic-like behavior and cerebellar dysfunction in Bmal1 mutant mice ameliorated by mTORC1 inhibition.

Although circadian and sleep disorders are frequently associated with autism spectrum disorders (ASD), it remains elusive whether clock gene disruption can lead to autistic-like phenotypes in animals. The essential clock gene Bmal1 has been associated with human sociability and its missense mutations are identified in ASD. Here we report that global Bmal1 deletion led to significant social impairments, excessive stereotyped and repetitive behaviors, as well as motor learning disabilities in mice, all of which resemble core behavioral deficits in ASD. Furthermore, aberrant cell density and immature morphology of dendritic spines were identified in the cerebellar Purkinje cells (PCs) of Bmal1 knockout (KO) mice. Electrophysiological recordings uncovered enhanced excitatory and inhibitory synaptic transmission and reduced firing rates in the PCs of Bmal1 KO mice. Differential expression of ASD- and ataxia-associated genes (Ntng2, Mfrp, Nr4a2, Thbs1, Atxn1, and Atxn3) and dysregulated pathways of translational control, including hyperactivated mammalian target of rapamycin complex 1 (mTORC1) signaling, were identified in the cerebellum of Bmal1 KO mice. Interestingly, the antidiabetic drug metformin reversed mTORC1 hyperactivation and alleviated major behavioral and PC deficits in Bmal1 KO mice. Importantly, conditional Bmal1 deletion only in cerebellar PCs was sufficient to recapitulate autistic-like behavioral and cellular changes akin to those identified in Bmal1 KO mice. Together, these results unveil a previously unidentified role for Bmal1 disruption in cerebellar dysfunction and autistic-like behaviors. Our findings provide experimental evidence supporting a putative role for dysregulation of circadian clock gene expression in the pathogenesis of ASD.

abc4pwm: affinity based clustering for position weight matrices in applications of DNA sequence analysis.

BACKGROUND: Transcription factor (TF) binding motifs are identified by high throughput sequencing technologies as means to capture Protein-DNA interactions. These motifs are often represented by consensus sequences in form of position weight matrices (PWMs). With ever-increasing pool of TF binding motifs from multiple sources, redundancy issues are difficult to avoid, especially when every source maintains its own database for collection. One solution can be to cluster biologically relevant or similar PWMs, whether coming from experimental detection or in silico predictions. However, there is a lack of efficient tools to cluster PWMs. Assessing quality of PWM clusters is yet another challenge. Therefore, new methods and tools are required to efficiently cluster PWMs and assess quality of clusters. RESULTS: A new Python package Affinity Based Clustering for Position Weight Matrices (abc4pwm) was developed. It efficiently clustered PWMs from multiple sources with or without using DNA-Binding Domain (DBD) information, generated a representative motif for each cluster, evaluated the clustering quality automatically, and filtered out incorrectly clustered PWMs. Additionally, it was able to update human DBD family database automatically, classified known human TF PWMs to the respective DBD family, and performed TF motif searching and motif discovery by a new ensemble learning approach. CONCLUSION: This work demonstrates applications of abc4pwm in the DNA sequence analysis for various high throughput sequencing data using ~ 1770 human TF PWMs. It recovered known TF motifs at gene promoters based on gene expression profiles (RNA-seq) and identified true TF binding targets for motifs predicted from ChIP-seq experiments. Abc4pwm is a useful tool for TF motif searching, clustering, quality assessment and integration in multiple types of sequence data analysis including RNA-seq, ChIP-seq and ATAC-seq.

Haploinsufficiency of a Circadian Clock Gene Bmal1 (Arntl or Mop3) Causes Brain-Wide mTOR Hyperactivation and Autism-like Behavioral Phenotypes in Mice.

Approximately 50-80% of children with autism spectrum disorders (ASDs) exhibit sleep problems, but the contribution of circadian clock dysfunction to the development of ASDs remains largely unknown. The essential clock gene Bmal1 (Arntl or Mop3) has been associated with human sociability, and its missense mutation is found in ASD. Our recent study found that Bmal1-null mice exhibit a variety of autism-like phenotypes. Here, we further investigated whether an incomplete loss of Bmal1 function could cause significant autism-like behavioral changes in mice. Our results demonstrated that heterozygous Bmal1 deletion (Bmal1+/-) reduced the Bmal1 protein levels by ~50-75%. Reduced Bmal1 expression led to decreased levels of clock proteins, including Per1, Per2, Cry 1, and Clock but increased mTOR activities in the brain. Accordingly, Bmal1+/- mice exhibited aberrant ultrasonic vocalizations during maternal separation, deficits in sociability and social novelty, excessive repetitive behaviors, impairments in motor coordination, as well as increased anxiety-like behavior. The novel object recognition memory remained intact. Together, these results demonstrate that haploinsufficiency of Bmal1 can cause autism-like behavioral changes in mice, akin to those identified in Bmal1-null mice. This study provides further experimental evidence supporting a potential role for disrupted clock gene expression in the development of ASD.

Progesterone Receptor Attenuates STAT1-Mediated IFN Signaling in Breast Cancer.

Why some tumors remain indolent and others progress to clinical relevance remains a major unanswered question in cancer biology. IFN signaling in nascent tumors, mediated by STAT1, is a critical step through which the surveilling immune system can recognize and destroy developing tumors. In this study, we have identified an interaction between the progesterone receptor (PR) and STAT1 in breast cancer cells. This interaction inhibited efficient IFN-induced STAT1 phosphorylation, as we observed a decrease in phospho-STAT1 in response to IFN treatment in PR-positive breast cancer cell lines. This phenotype was further potentiated in the presence of PR ligand. In human breast cancer samples, PR-positive tumors exhibited lower levels of phospho-STAT1 as compared with their PR-negative counterparts, indicating that this phenotype translates to human tumors. Breast cancer cells lacking PR exhibited higher levels of IFN-stimulated gene (ISG) RNA, the transcriptional end point of IFN activation, indicating that unliganded PR alone could decrease transcription of ISGs. Moreover, the absence of PR led to increased recruitment of STAT1, STAT2, and IRF9 (key transcription factors necessary for ISG transcription) to ISG promoters. These data indicate that PR, both in the presence and absence of ligand, attenuates IFN-induced STAT1 signaling, culminating in significantly abrogated activation of genes transcribed in response to IFNs. PR-positive tumors may use downregulation of STAT1-mediated IFN signaling to escape immune surveillance, leading to the development of clinically relevant tumors. Selective immune evasion of PR-positive tumors may be one explanation as to why over 65% of breast cancers are PR positive at the time of diagnosis.

Diverse AR-V7 cistromes in castration-resistant prostate cancer are governed by HoxB13.

The constitutively active androgen receptor (AR) splice variant 7 (AR-V7) plays an important role in the progression of castration-resistant prostate cancer (CRPC). Although biomarker studies established the role of AR-V7 in resistance to AR-targeting therapies, how AR-V7 mediates genomic functions in CRPC remains largely unknown. Using a ChIP-exo approach, we show AR-V7 binds to distinct genomic regions and recognizes a full-length androgen-responsive element in CRPC cells and patient tissues. Remarkably, we find dramatic differences in AR-V7 cistromes across diverse CRPC cells and patient tissues, regulating different target gene sets involved in CRPC progression. Surprisingly, we discover that HoxB13 is universally required for and colocalizes with AR-V7 binding to open chromatin across CRPC genomes. HoxB13 pioneers AR-V7 binding through direct physical interaction, and collaborates with AR-V7 to up-regulate target oncogenes. Transcriptional coregulation by HoxB13 and AR-V7 was further supported by their coexpression in tumors and circulating tumor cells from CRPC patients. Importantly, HoxB13 silencing significantly decreases CRPC growth through inhibition of AR-V7 oncogenic function. These results identify HoxB13 as a pivotal upstream regulator of AR-V7-driven transcriptomes that are often cell context-dependent in CRPC, suggesting that HoxB13 may serve as a therapeutic target for AR-V7-driven prostate tumors.

Epigenomics-based identification of oestrogen-regulated long noncoding RNAs in ER+ breast cancer.

Breast cancer is one of the most prevalent cancers in women worldwide. Through the regulation of many coding and non-coding target genes, oestrogen (E2 or 17ß-oestradiol) and its nuclear receptor ERα play important roles in breast cancer development and progression. Despite the astounding advances in our understanding of oestrogen-regulated coding genes over the past decades, our knowledge on oestrogen-regulated non-coding targets has just begun to expand. Here we leverage epigenomic approaches to systematically analyse oestrogen-regulated long non-coding RNAs (lncRNAs). Similar to the coding targets of ERα, the transcription of oestrogen-regulated lncRNAs correlates with the activation status of ERα enhancers, measured by eRNA production, chromatin accessibility, and the occupancy of the enhancer regulatory components including P300, MED1, and ARID1B. Our 3D chromatin architecture analyses suggest that lncRNAs and their neighbouring E2-resonsive coding genes, exemplified by LINC00160 and RUNX1, might be regulated as a 3D structural unit resulted from enhancer-promoter interactions. Finally, we evaluated the expression levels of LINC00160 and RUNX1 in various types of breast cancer and found that their expression positively correlated with the survival rate in ER+ breast cancer patients, implying that the oestrogen-regulated LINC00160 and its neighbouring RUNX1 might represent potential biomarkers for ER+ breast cancers.

Three-dimensional analysis reveals altered chromatin interaction by enhancer inhibitors harbors TCF7L2-regulated cancer gene signature.

Distal regulatory elements influence the activity of gene promoters through chromatin looping. Chromosome conformation capture (3C) methods permit identification of chromatin contacts across different regions of the genome. However, due to limitations in the resolution of these methods, the detection of functional chromatin interactions remains a challenge. In the current study, we employ an integrated approach to define and characterize the functional chromatin contacts of human pancreatic cancer cells. We applied tethered chromatin capture to define classes of chromatin domains on a genome-wide scale. We identified three types of structural domains (topologically associated, boundary, and gap) and investigated the functional relationships of these domains with respect to chromatin state and gene expression. We uncovered six distinct sub-domains associated with epigenetic states. Interestingly, specific epigenetically active domains are sensitive to treatment with histone acetyltransferase (HAT) inhibitors and decrease in H3K27 acetylation levels. To examine whether the subdomains that change upon drug treatment are functionally linked to transcription factor regulation, we compared TCF7L2 chromatin binding and gene regulation to HAT inhibition. We identified a subset of coding RNA genes that together can stratify pancreatic cancer patients into distinct survival groups. Overall, this study describes a process to evaluate the functional features of chromosome architecture and reveals the impact of epigenetic inhibitors on chromosome architecture and identifies genes that may provide insight into disease outcome.

BRCA1 mutations attenuate super-enhancer function and chromatin looping in haploinsufficient human breast epithelial cells.

BACKGROUND: BRCA1-associated breast cancer originates from luminal progenitor cells. BRCA1 functions in multiple biological processes, including double-strand break repair, replication stress suppression, transcriptional regulation, and chromatin reorganization. While non-malignant cells carrying cancer-predisposing BRCA1 mutations exhibit increased genomic instability, it remains unclear whether BRCA1 haploinsufficiency affects transcription and chromatin dynamics in breast epithelial cells. METHODS: H3K27ac-associated super-enhancers were compared in primary breast epithelial cells from BRCA1 mutation carriers (BRCA1mut/+) and non-carriers (BRCA1+/+). Non-tumorigenic MCF10A breast epithelial cells with engineered BRCA1 haploinsufficiency were used to confirm the H3K27ac changes. The impact of BRCA1 mutations on enhancer function and enhancer-promoter looping was assessed in MCF10A cells. RESULTS: Here, we show that primary mammary epithelial cells from women with BRCA1 mutations display significant loss of H3K27ac-associated super-enhancers. These BRCA1-dependent super-enhancers are enriched with binding motifs for the GATA family. Non-tumorigenic BRCA1mut/+ MCF10A cells recapitulate the H3K27ac loss. Attenuated histone mark and enhancer activity in these BRCA1mut/+ MCF10A cells can be partially restored with wild-type BRCA1. Furthermore, chromatin conformation analysis demonstrates impaired enhancer-promoter looping in BRCA1mut/+ MCF10A cells. CONCLUSIONS: H3K27ac-associated super-enhancer loss is a previously unappreciated functional deficiency in ostensibly normal BRCA1 mutation-carrying breast epithelium. Our findings offer new mechanistic insights into BRCA1 mutation-associated transcriptional and epigenetic abnormality in breast epithelial cells and tissue/cell lineage-specific tumorigenesis.

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.

Genome-wide analysis reveals positional-nucleosome-oriented binding pattern of pioneer factor FOXA1.

The compaction of nucleosomal structures creates a barrier for DNA-binding transcription factors (TFs) to access their cognate cis-regulatory elements. Pioneer factors (PFs) such as FOXA1 are able to directly access these cis-targets within compact chromatin. However, how these PFs interplay with nucleosomes remains to be elucidated, and is critical for us to understand the underlying mechanism of gene regulation. Here, we have conducted a computational analysis on a strand-specific paired-end ChIP-exo (termed as ChIP-ePENS) data of FOXA1 in LNCaP cells by our novel algorithm ePEST. We find that FOXA1 chromatin binding occurs via four distinct border modes (or footprint boundary patterns), with a preferential footprint boundary patterns relative to FOXA1 motif orientation. In addition, from this analysis three fundamental nucleotide positions (oG, oS and oH) emerged as major determinants for blocking exo-digestion and forming these four distinct border modes. By integrating histone MNase-seq data, we found an astonishingly consistent, 'well-positioned' configuration occurs between FOXA1 motifs and dyads of nucleosomes genome-wide. We further performed ChIP-seq of eight chromatin remodelers and found an increased occupancy of these remodelers on FOXA1 motifs for all four border modes (or footprint boundary patterns), indicating the full occupancy of FOXA1 complex on the three blocking sites (oG, oS and oH) likely produces an active regulatory status with well-positioned phasing for protein binding events. Together, our results suggest a positional-nucleosome-oriented accessing model for PFs seeking target motifs, in which FOXA1 can examine each underlying DNA nucleotide and is able to sense all potential motifs regardless of whether they face inward or outward from histone octamers along the DNA helix axis.

Integrative analysis identifies targetable CREB1/FoxA1 transcriptional co-regulation as a predictor of prostate cancer recurrence.

Identifying prostate cancer-driving transcription factors (TFs) in addition to the androgen receptor promises to improve our ability to effectively diagnose and treat this disease. We employed an integrative genomics analysis of master TFs CREB1 and FoxA1 in androgen-dependent prostate cancer (ADPC) and castration-resistant prostate cancer (CRPC) cell lines, primary prostate cancer tissues and circulating tumor cells (CTCs) to investigate their role in defining prostate cancer gene expression profiles. Combining genome-wide binding site and gene expression profiles we define CREB1 as a critical driver of pro-survival, cell cycle and metabolic transcription programs. We show that CREB1 and FoxA1 co-localize and mutually influence each other's binding to define disease-driving transcription profiles associated with advanced prostate cancer. Gene expression analysis in human prostate cancer samples found that CREB1/FoxA1 target gene panels predict prostate cancer recurrence. Finally, we showed that this signaling pathway is sensitive to compounds that inhibit the transcription co-regulatory factor MED1. These findings not only reveal a novel, global transcriptional co-regulatory function of CREB1 and FoxA1, but also suggest CREB1/FoxA1 signaling is a targetable driver of prostate cancer progression and serves as a biomarker of poor clinical outcomes.

Transcription factor-associated combinatorial epigenetic pattern reveals higher transcriptional activity of TCF7L2-regulated intragenic enhancers.

BACKGROUND: Recent studies have suggested that combinations of multiple epigenetic modifications are essential for controlling gene expression. Despite numerous computational approaches have been developed to decipher the combinatorial epigenetic patterns or "epigenetic code", none of them has explicitly addressed the relationship between a specific transcription factor (TF) and the patterns. METHODS: Here, we developed a novel computational method, T-cep, for annotating chromatin states associated with a specific TF. T-cep is composed of three key consecutive modules: (i) Data preprocessing, (ii) HMM training, and (iii) Potential TF-states calling. RESULTS: We evaluated T-cep on a TCF7L2-omics data. Unexpectedly, our method has uncovered a novel set of TCF7L2-regulated intragenic enhancers missed by other software tools, where the associated genes exert the highest gene expression. We further used siRNA knockdown, Co-transfection, RT-qPCR and Luciferase Reporter Assay not only to validate the accuracy and efficiency of prediction by T-cep, but also to confirm the functionality of TCF7L2-regulated enhancers in both MCF7 and PANC1 cells respectively. CONCLUSIONS: Our study for the first time at a genome-wide scale reveals the enhanced transcriptional activity of cell-type-specific TCF7L2 intragenic enhancers in regulating gene expression.

Computational analysis reveals a correlation of exon-skipping events with splicing, transcription and epigenetic factors.

Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, 'skipping', exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon-intron boundary establishment leading to skipping events.

Three-tiered role of the pioneer factor GATA2 in promoting androgen-dependent gene expression in prostate cancer.

In prostate cancer, androgen receptor (AR) binding and androgen-responsive gene expression are defined by hormone-independent binding patterns of the pioneer factors FoxA1 and GATA2. Insufficient evidence of the mechanisms by which GATA2 contributes to this process precludes complete understanding of a key determinant of tissue-specific AR activity. Our observations suggest that GATA2 facilitates androgen-responsive gene expression by three distinct modes of action. By occupying novel binding sites within the AR gene locus, GATA2 positively regulates AR expression before and after androgen stimulation. Additionally, GATA2 engages AR target gene enhancers prior to hormone stimulation, producing an active and accessible chromatin environment via recruitment of the histone acetyltransferase p300. Finally, GATA2 functions in establishing and/or sustaining basal locus looping by recruiting the Mediator subunit MED1 in the absence of androgen. These mechanisms may contribute to the generally positive role of GATA2 in defining AR genome-wide binding patterns that determine androgen-responsive gene expression profiles. We also find that GATA2 and FoxA1 exhibit both independent and codependent co-occupancy of AR target gene enhancers. Identifying these determinants of AR transcriptional activity may provide a foundation for the development of future prostate cancer therapeutics that target pioneer factor function.

DNA methylation screening of primary prostate tumors identifies SRD5A2 and CYP11A1 as candidate markers for assessing risk of biochemical recurrence.

BACKGROUND: Altered DNA methylation in CpG islands of gene promoters has been implicated in prostate cancer (PCa) progression and can be used to predict disease outcome. In this study, we determine whether methylation changes of androgen biosynthesis pathway (ABP)-related genes in patients' plasma cell-free DNA (cfDNA) can serve as prognostic markers for biochemical recurrence (BCR). METHODS: Methyl-binding domain capture sequencing (MBDCap-seq) was used to identify differentially methylated regions (DMRs) in primary tumors of patients who subsequently developed BCR or not, respectively. Methylation pyrosequencing of candidate loci was validated in cfDNA samples of 86 PCa patients taken at and/or post-radical prostatectomy (RP) using univariate and multivariate prediction analyses. RESULTS: Putative DMRs in 13 of 30 ABP-related genes were found between tumors of BCR (n = 12) versus no evidence of disease (NED) (n = 15). In silico analysis of The Cancer Genome Atlas data confirmed increased DNA methylation of two loci-SRD5A2 and CYP11A1, which also correlated with their decreased expression, in tumors with subsequent BCR development. Their aberrant cfDNA methylation was also associated with detectable levels of PSA taken after patients' post-RP. Multivariate analysis of the change in cfDNA methylation at all of CpG sites measured along with patient's treatment history predicted if a patient will develop BCR with 77.5% overall accuracy. CONCLUSIONS: Overall, increased DNA methylation of SRD5A2 and CYP11A1 related to androgen biosynthesis functions may play a role in BCR after patients' RP. The correlation between aberrant cfDNA methylation and detectable PSA in post-RP further suggests their utility as predictive markers for PCa recurrence. .