Dynamic nucleosome-depleted regions at androgen receptor enhancers in the absence of ligand in prostate cancer cells

Nucleosome positioning at transcription start sites is known to regulate gene expression by altering DNA accessibility to transcription factors; however, its role at enhancers is poorly understood. We investigated nucleosome positioning at the androgen receptor (AR) enhancers of TMPRSS2, KLK2, and KLK3/PSA in prostate cancer cells. Surprisingly, a population of enhancer modules in androgen-deprived cultures showed nucleosome-depleted regions (NDRs) in all three loci. Under androgen-deprived conditions, NDRs at the TMPRSS2 enhancer were maintained by the pioneer AR transcriptional collaborator GATA-2. Androgen treatment resulted in AR occupancy, an increased number of enhancer modules with NDRs without changes in footprint width, increased levels of histone H3 acetylation (AcH3), and dimethylation (H3K4me2) at nucleosomes flanking the NDRs. Our data suggest that, in the absence of ligand, AR enhancers exist in an equilibrium in which a percentage of modules are occupied by nucleosomes while others display NDRs. We propose that androgen treatment leads to the disruption of the equilibrium toward a nucleosome-depleted state, rather than to enhancer de novo “remodeling.” This allows the recruitment of histone modifiers, chromatin remodelers, and ultimately gene activation. The “receptive” state described here could help explain AR signaling activation under very low ligand concentrations

Familial defective apolipoprotein-B is rare in hypercholesterolaemic South African Afrikaners, coloureds and Indians

No Abstract

Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells

Abstract\ud \ud Background\ud The precise nature of how cell type specific chromatin structures at enhancer sites affect gene expression is largely unknown. Here we identified cell type specific enhancers coupled with gene expression in two different types of breast epithelial cells, HMEC (normal breast epithelial cells) and MDAMB231 (triple negative breast cancer cell line).\ud \ud \ud Results\ud Enhancers were defined by modified neighboring histones [using chromatin immunoprecipitation followed by sequencing (ChIP-seq)] and nucleosome depletion [using formaldehyde-assisted isolation of regulatory elements followed by sequencing (FAIRE-seq)]. Histone modifications at enhancers were related to the expression levels of nearby genes up to 750 kb away. These expression levels were correlated with enhancer status (poised or active), defined by surrounding histone marks. Furthermore, about fifty percent of poised and active enhancers contained nucleosome-depleted regions. We also identified response element motifs enriched at these enhancer sites that revealed key transcription factors (e.g. TP63) likely involved in regulating breast epithelial enhancer-mediated gene expression. By utilizing expression data, potential target genes of more than 600 active enhancers were identified. These genes were involved in proteolysis, epidermis development, cell adhesion, mitosis, cell cycle, and DNA replication.\ud \ud \ud Conclusions\ud These findings facilitate the understanding of epigenetic regulation specifically, such as the relationships between regulatory elements and gene expression and generally, how breast epithelial cellular phenotypes are determined by cell type specific enhancers.National Institutes of Health [R01 CA136924 to GAC, T32CA009320 to HN]Genetic Associations and Mechanisms in Oncology (GAME-ON

Principles for the post-GWAS functional characterisation of risk loci

Several challenges lie ahead in assigning functionality to susceptibility SNPs. For example, most effect sizes are small relative to effects seen in monogenic diseases, with per allele odds ratios usually ranging from 1.15 to 1.3. It is unclear whether current molecular biology methods have enough resolution to differentiate such small effects. Our objective here is therefore to provide a set of recommendations to optimize the allocation of effort and resources in order maximize the chances of elucidating the functional contribution of specific loci to the disease phenotype. It has been estimated that 88% of currently identified disease-associated SNP are intronic or intergenic. Thus, in this paper we will focus our attention on the analysis of non-coding variants and outline a hierarchical approach for post-GWAS functional studies

Genomic Androgen Receptor-Occupied Regions with Different Functions, Defined by Histone Acetylation, Coregulators and Transcriptional Capacity

Background: The androgen receptor (AR) is a steroid-activated transcription factor that binds at specific DNA locations and plays a key role in the etiology of prostate cancer. While numerous studies have identified a clear connection between AR binding and expression of target genes for a limited number of loci, high-throughput elucidation of these sites allows for a deeper understanding of the complexities of this process. Methodology/Principal Findings: We have mapped 189 AR occupied regions (ARORs) and 1,388 histone H3 acetylation (AcH3) loci to a 3 % continuous stretch of human genomic DNA using chromatin immunoprecipitation (ChIP) microarray analysis. Of 62 highly reproducible ARORs, 32 (52%) were also marked by AcH3. While the number of ARORs detected in prostate cancer cells exceeded the number of nearby DHT-responsive genes, the AcH3 mark defined a subclass of ARORs much more highly associated with such genes – 12 % of the genes flanking AcH3+ARORs were DHT-responsive, compared to only 1 % of genes flanking AcH32ARORs. Most ARORs contained enhancer activities as detected in luciferase reporte

Runx2 transcriptome of prostate cancer cells: insights into invasiveness and bone metastasis

<p>Abstract</p> <p>Background</p> <p>Prostate cancer (PCa) cells preferentially metastasize to bone at least in part by acquiring osteomimetic properties. Runx2, an osteoblast master transcription factor, is aberrantly expressed in PCa cells, and promotes their metastatic phenotype. The transcriptional programs regulated by Runx2 have been extensively studied during osteoblastogenesis, where it activates or represses target genes in a context-dependent manner. However, little is known about the gene regulatory networks influenced by Runx2 in PCa cells. We therefore investigated genome wide mRNA expression changes in PCa cells in response to Runx2.</p> <p>Results</p> <p>We engineered a C4-2B PCa sub-line called C4-2B/Rx2^dox, in which Doxycycline (Dox) treatment stimulates Runx2 expression from very low to levels observed in other PCa cells. Transcriptome profiling using whole genome expression array followed by <it>in silico </it>analysis indicated that Runx2 upregulated a multitude of genes with prominent cancer associated functions. They included secreted factors (CSF2, SDF-1), proteolytic enzymes (MMP9, CST7), cytoskeleton modulators (SDC2, Twinfilin, SH3PXD2A), intracellular signaling molecules (DUSP1, SPHK1, RASD1) and transcription factors (Sox9, SNAI2, SMAD3) functioning in epithelium to mesenchyme transition (EMT), tissue invasion, as well as homing and attachment to bone. Consistent with the gene expression data, induction of Runx2 in C4-2B cells enhanced their invasiveness. It also promoted cellular quiescence by blocking the G1/S phase transition during cell cycle progression. Furthermore, the cell cycle block was reversed as Runx2 levels declined after Dox withdrawal.</p> <p>Conclusions</p> <p>The effects of Runx2 in C4-2B/Rx2^doxcells, as well as similar observations made by employing LNCaP, 22RV1 and PC3 cells, highlight multiple mechanisms by which Runx2 promotes the metastatic phenotype of PCa cells, including tissue invasion, homing to bone and induction of high bone turnover. Runx2 is therefore an attractive target for the development of novel diagnostic, prognostic and therapeutic approaches to PCa management. Targeting Runx2 may prove more effective than focusing on its individual downstream genes and pathways.</p

Comprehensive functional annotation of 77 prostate cancer risk loci.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process

Functional Enhancers at the Gene-Poor 8q24 Cancer-Linked Locus

Multiple discrete regions at 8q24 were recently shown to contain alleles that predispose to many cancers including prostate, breast, and colon. These regions are far from any annotated gene and their biological activities have been unknown. Here we profiled a 5-megabase chromatin segment encompassing all the risk regions for RNA expression, histone modifications, and locations occupied by RNA polymerase II and androgen receptor (AR). This led to the identification of several transcriptional enhancers, which were verified using reporter assays. Two enhancers in one risk region were occupied by AR and responded to androgen treatment; one contained a single nucleotide polymorphism (rs11986220) that resides within a FoxA1 binding site, with the prostate cancer risk allele facilitating both stronger FoxA1 binding and stronger androgen responsiveness. The study reported here exemplifies an approach that may be applied to any risk-associated allele in non-protein coding regions as it emerges from genome-wide association studies to better understand the genetic predisposition of complex diseases

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus.

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439

Chromatin accessibility reveals insights into androgen receptor activation and transcriptional specificity

Abstract Background Epigenetic mechanisms such as chromatin accessibility impact transcription factor binding to DNA and transcriptional specificity. The androgen receptor (AR), a master regulator of the male phenotype and prostate cancer pathogenesis, acts primarily through ligand-activated transcription of target genes. Although several determinants of AR transcriptional specificity have been elucidated, our understanding of the interplay between chromatin accessibility and AR function remains incomplete. Results We used deep sequencing to assess chromatin structure via DNase I hypersensitivity and mRNA abundance, and paired these datasets with three independent AR ChIP-seq datasets. Our analysis revealed qualitative and quantitative differences in chromatin accessibility that corresponded to both AR binding and an enrichment of motifs for potential collaborating factors, one of which was identified as SP1. These quantitative differences were significantly associated with AR-regulated mRNA transcription across the genome. Base-pair resolution of the DNase I cleavage profile revealed three distinct footprinting patterns associated with the AR-DNA interaction, suggesting multiple modes of AR interaction with the genome. Conclusions In contrast with other DNA-binding factors, AR binding to the genome does not only target regions that are accessible to DNase I cleavage prior to hormone induction. AR binding is invariably associated with an increase in chromatin accessibility and, consequently, changes in gene expression. Furthermore, we present the first in vivo evidence that a significant fraction of AR binds only to half of the full AR DNA motif. These findings indicate a dynamic quantitative relationship between chromatin structure and AR-DNA binding that impacts AR transcriptional specificity