Mark D. Long

Integration of VDR genome wide binding and GWAS genetic variation data reveals co-occurrence of VDR and NF-κB binding that is linked to immune phenotypes

BackgroundThe nuclear hormone receptor superfamily acts as a genomic sensor of diverse signals. Their actions are often intertwined with other transcription factors. Nuclear hormone receptors are targets for many therapeutic drugs, and include the vitamin D receptor (VDR). VDR signaling is pleotropic, being implicated in calcaemic function, antibacterial actions, growth control, immunomodulation and anti-cancer actions. Specifically, we hypothesized that the biologically significant relationships between the VDR transcriptome and phenotype-associated biology could be discovered by integrating the known VDR transcription factor binding sites and all published trait- and disease-associated SNPs. By integrating VDR genome-wide binding data (ChIP-seq) with the National Human Genome Research Institute (NHGRI) GWAS catalog of SNPs we would see where and which target gene interactions and pathways are impacted by inherited genetic variation in VDR binding sites, indicating which of VDR’s multiple functions are most biologically significant.ResultsTo examine how genetic variation impacts VDR function we overlapped 23,409 VDR genomic binding peaks from six VDR ChIP-seq datasets with 191,482 SNPs, derived from GWAS-significant SNPs (Lead SNPs) and their correlated variants (r2 > 0.8) from HapMap3 and the 1000 genomes project. In total, 574 SNPs (71 Lead and 503 SNPs in linkage disequilibrium with Lead SNPs) were present at VDR binding loci and associated with 211 phenotypes. For each phenotype a hypergeometric test was used to determine if SNPs were enriched at VDR binding sites. Bonferroni correction for multiple testing across the 211 phenotypes yielded 42 SNPs that were either disease- or phenotype-associated with seven predominately immune related including self-reported allergy; esophageal cancer was the only cancer phenotype. Motif analyses revealed that only two of these 42 SNPs reside within a canonical VDR binding site (DR3 motif), and that 1/3 of the 42 SNPs significantly impacted binding and gene regulation by other transcription factors, including NF-κB. This suggests a plausible link for the potential cross-talk between VDR and NF-κB.ConclusionsThese analyses showed that VDR peaks are enriched for SNPs associated with immune phenotypes suggesting that VDR immunomodulatory functions are amongst its most important actions. The enrichment of genetic variation in non-DR3 motifs suggests a significant role for the VDR to bind in multimeric complexes containing other transcription factors that are the primary DNA binding component. Our work provides a framework for the combination of ChIP-seq and GWAS findings to provide insight into the underlying phenotype-associated biology of a given transcription factor.

Vitamin D Receptor and RXR in the Post-Genomic Era

Following the elucidation of the human genome and components of the epigenome, it is timely to revisit what is known of vitamin D receptor (VDR) function. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. Investigators also considered VDR regulation of non‐protein coding RNA and again, cell and time dependency was observed. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co‐regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. Alternative splicing in the trasncriptome has emerged as a critical process in transcriptional control and there is evidence of the VDR interacting with components of the splicesome. ChIP‐Seq approaches have proved to be pivotal to reveal the diversity of the VDR binding choices across cell types and following treatment, and have revealed that the majority of these are non‐canonical in nature. The underlying causes driving the diversity of VDR binding choices remain enigmatic. Finally, genetic variation has emerged as important to impact the transcription factor affinity towards genomic binding sites, and recently the impact of this on VDR function has begun to be considered. J. Cell. Physiol. 230: 758–766, 2015.

Cooperative behavior of the nuclear receptor superfamily and its deregulation in prostate cancer

The current study aimed to assess the topology of the nuclear receptor (NR) superfamily in normal prostate epithelial cells and its distortion in prostate cancer. Both in vitro and in silico approaches were utilized to profile NRs expressed in non-malignant RWPE-1 cells, which were subsequently investigated by treating cells with 132 binary NR ligand combinations. Nine significant cooperative interactions emerged including both superadditive [22(R)-hydroxycholesterol and eicosatetraenoic acid] and subadditive [1α,25(OH)2D3 and chenodeoxycholic acid] cellular responses, which could be explained in part by cooperative control of cell-cycle progression and candidate gene expression. In addition, publicly available data were employed to assess NR expression in human prostate tissue. Common and significant loss of NR superfamily expression was established in publicly available data from prostate tumors, in part predicting parallel distortion of targeting microRNA. These findings suggest that the NR superfamily in the prostate cooperatively integrates signals from dietary, hormonal and metabolic cues, and is significantly distorted in prostate cancer.

VDR regulation of microRNA differs across prostate cell models suggesting extremely flexible control of transcription

The Vitamin D Receptor (VDR) is a member of the nuclear receptor superfamily and is of therapeutic interest in cancer and other settings. Regulation of microRNA (miRNA) by the VDR appears to be important to mediate its actions, for example, to control cell growth. To identify if and to what extent VDR-regulated miRNA patterns change in prostate cancer progression, we undertook miRNA microarray analyses in 7 cell models representing non-malignant and malignant prostate cells (RWPE-1, RWPE-2, HPr1, HPr1AR, LNCaP, LNCaP-C4–2, and PC-3). To focus on primary VDR regulatory events, we undertook expression analyses after 30 minutes treatment with 1α,25(OH)2D3. Across all models, 111 miRNAs were significantly modulated by 1α,25(OH)2D3 treatment. Of these, only 5 miRNAs were modulated in more than one cell model, and of these, only 3 miRNAs were modulated in the same direction. The patterns of miRNA regulation, and the networks they targeted, significantly distinguished the different cell types. Integration of 1α,25(OH)2D3-regulated miRNAs with published VDR ChIP-seq data showed significant enrichment of VDR peaks in flanking regions of miRNAs. Furthermore, mRNA and miRNA expression analyses in non-malignant RWPE-1 cells revealed patterns of miRNA and mRNA co-regulation; specifically, 13 significant reciprocal patterns were identified and these patterns were also observed in TCGA prostate cancer data. Lastly, motif search analysis revealed differential motif enrichment within VDR peaks flanking mRNA compared to miRNA genes. Together, this study revealed that miRNAs are rapidly regulated in a highly cell-type specific manner, and are significantly co-integrated with mRNA regulation.

Integrative genomic analysis in K562 chronic myelogenous leukemia cells reveals that proximal NCOR1 binding positively regulates genes that govern erythroid differentiation and Imatinib sensitivity

To define the functions of NCOR1 we developed an integrative analysis that combined ENCODE and NCI-60 data, followed by in vitro validation. NCOR1 and H3K9me3 ChIP-Seq, FAIRE-seq and DNA CpG methylation interactions were related to gene expression using bootstrapping approaches. Most NCOR1 combinations (24/44) were associated with significantly elevated level expression of protein coding genes and only very few combinations related to gene repression. DAVIDs biological process annotation revealed that elevated gene expression was uniquely associated with acetylation and ETS binding. A matrix of gene and drug interactions built on NCI-60 data identified that Imatinib significantly targeted the NCOR1 governed transcriptome. Stable knockdown of NCOR1 in K562 cells slowed growth and significantly repressed genes associated with NCOR1 cistrome, again, with the GO terms acetylation and ETS binding, and significantly dampened sensitivity to Imatinib-induced erythroid differentiation. Mining public microarray data revealed that NCOR1-targeted genes were significantly enriched in Imatinib response gene signatures in cell lines and chronic myelogenous leukemia (CML) patients. These approaches integrated cistrome, transcriptome and drug sensitivity relationships to reveal that NCOR1 function is surprisingly most associated with elevated gene expression, and that these targets, both in CML cell lines and patients, associate with sensitivity to Imatinib.

The Genomic Impact of DNA CpG Methylation on Gene Expression; Relationships in Prostate Cancer

The process of DNA CpG methylation has been extensively investigated for over 50 years and revealed associations between changing methylation status of CpG islands and gene expression. As a result, DNA CpG methylation is implicated in the control of gene expression in developmental and homeostasis processes, as well as being a cancer-driver mechanism. The development of genome-wide technologies and sophisticated statistical analytical approaches has ushered in an era of widespread analyses, for example in the cancer arena, of the relationships between altered DNA CpG methylation, gene expression, and tumor status. The remarkable increase in the volume of such genomic data, for example, through investigators from the Cancer Genome Atlas (TCGA), has allowed dissection of the relationships between DNA CpG methylation density and distribution, gene expression, and tumor outcome. In this manner, it is now possible to test that the genome-wide correlations are measurable between changes in DNA CpG methylation and gene expression. Perhaps surprisingly is that these associations can only be detected for hundreds, but not thousands, of genes, and the direction of the correlations are both positive and negative. This, perhaps, suggests that CpG methylation events in cancer systems can act as disease drivers but the effects are possibly more restricted than suspected. Additionally, the positive and negative correlations suggest direct and indirect events and an incomplete understanding. Within the prostate cancer TCGA cohort, we examined the relationships between expression of genes that control DNA methylation, known targets of DNA methylation and tumor status. This revealed that genes that control the synthesis of S-adenosyl-l-methionine (SAM) associate with altered expression of DNA methylation targets in a subset of aggressive tumors.

Pan-Cancer Analyses of the Nuclear Receptor Superfamily

Nuclear receptors (NR) act as an integrated conduit for environmental and hormonal signals to govern genomic responses, which relate to cell fate decisions. We review how their integrated actions with each other, shared co-factors and other transcription factors are disrupted in cancer. Steroid hormone nuclear receptors are oncogenic drivers in breast and prostate cancer and blockade of signaling is a major therapeutic goal. By contrast to blockade of receptors, in other cancers enhanced receptor function is attractive, as illustrated initially with targeting of retinoic acid receptors in leukemia. In the post-genomic era large consortia, such as The Cancer Genome Atlas, have developed a remarkable volume of genomic data with which to examine multiple aspects of nuclear receptor status in a pan-cancer manner. Therefore to extend the review of NR function we have also undertaken bioinformatics analyses of NR expression in over 3000 tumors, spread across six different tumor types (bladder, breast, colon, head and neck, liver and prostate). Specifically, to ask how the NR expression was distorted (altered expression, mutation and CNV) we have applied bootstrapping approaches to simulate data for comparison, and also compared these NR findings to 12 other transcription factor families. Nuclear receptors were uniquely and uniformly downregulated across all six tumor types, more than predicted by chance. These approaches also revealed that each tumor type had a specific NR expression profile but these were most similar between breast and prostate cancer. Some NRs were down-regulated in at least five tumor types (e.g. NR3C2/MR and NR5A2/LRH-1)) whereas others were uniquely down-regulated in one tumor (e.g. NR1B3/RARG). The downregulation was not driven by copy number variation or mutation and epigenetic mechanisms maybe responsible for the altered nuclear receptor expression.

Integrative genomic approaches to dissect clinically-significant relationships between the VDR cistrome and gene expression in primary colon cancer

Recently, we undertook a pan-cancer analyses of the nuclear hormone receptor (NR) superfamily in The Cancer Genome Atlas (TCGA), and revealed that the vitamin D receptor (NR1I1/VDR) was commonly and significantly down-regulated specifically in colon adenocarcinoma cohort (COAD). To examine the consequence of down-regulated VDR expression we re-analyzed VDR chromatin immunoprecipitation sequencing (ChIP-Seq) data from LS180 colon cancer cells (GSE31939). This analysis identified 1809 loci that displayed significant (p.adj<0.01) differential binding of the VDR in response 1,25(OH)2D3 treatment; 947 peaks annotated to 672 genes. We examined expression patterns in the COAD cohort of 286 tumors compared to 41 normal samples and revealed that VDR bound genes were significantly positively correlated to VDR expression compared to the background transcriptome, suggesting direct regulation by VDR. Gene set enrichment analyses revealed significant enrichment for genes known to be regulated by a number of other transcription factors including SMADs and JUN. Filtering VDR associated genes for those that were commonly and significantly altered in COAD revealed a cohort of 27 differentially expressed genes. The expression patterns of these genes clustered tumors and significantly associated with disease free survival. For instance, males with low expression of Lectin, Galactoside Binding Soluble 4 (LGALS4, encodes the colon tumor suppressor, Galactin 4) had significantly shorted disease free survival. These analyses suggest that reduced expression of VDR in colon cancer (but neither loss nor mutation) changes the actions of the VDR by both dampening the expression of tumor suppressors (e.g. LGALS4) whilst either stabilizing or not down-regulating expression of oncogenes (e.g. Carbonic Anhydrase 9 (CA9)). These integrative genomic approaches are relatively generic and applicable to the study of any transcription factor.

Dietary folate levels alter the kinetics and molecular mechanism of prostate cancer recurrence in the CWR22 model

Folate impacts the genome and epigenome by feeding into one-carbon metabolism to produce critical metabolites, deoxythymidine monophosphate and s-adenosylmethionine. The impact of folate exposure and intervention timing on cancer progression remains controversial. Due to polyamine metabolism’s extraordinary biosynthetic flux in prostate cancer (CaP) we demonstrated androgen stimulated CaP is susceptible to dietary folate deficiency. We hypothesized dietary folate levels may also affect castration recurrent CaP. We used the CWR22 human xenograft model which recurs following androgen withdrawal. Engrafted mice were fed a folate depleted or supplemented diet beginning at androgen withdrawal, or prior to xenograft implantation. Both folate depletion and supplementation at the time of withdrawal significantly decreased recurrence incidence. Folate supplementation prior to xenograft implantation increased time to recurrence, suggesting a protective role. By contrast, folate depleted recurrent tumors exhibited transcriptional adaptive responses that maintained high polyamine levels at the expense of increased DNA damage and DNA methylation alterations. Mining of publically available data demonstrated folate related pathways are exceptionally dysregulated in human CaP, which correlated with decreased time to biochemical recurrence. These findings highlight the potential for novel therapeutic interventions that target these metabolic pathways in CaP and provide a rationale to apply such strategies alongside androgen withdrawal.

The miR-96 and RARγ signaling axis governs androgen signaling and prostate cancer progression

Expression levels of retinoic acid receptor gamma (NR1B3/RARG, encodes RARγ) are commonly reduced in prostate cancer (PCa). Therefore, we sought to establish the cellular and gene regulatory consequences of reduced RARγ expression, and determine RARγ regulatory mechanisms. RARG shRNA approaches in non-malignant (RWPE-1 and HPr1-AR) and malignant (LNCaP) prostate models revealed that reducing RARγ levels, rather than adding exogenous retinoid ligand, had the greatest impact on prostate cell viability and gene expression. ChIP-Seq defined the RARγ cistrome, which was significantly enriched at active enhancers associated with AR binding sites. Reflecting a significant genomic role for RARγ to regulate androgen signaling, RARγ knockdown in HPr1-AR cells significantly regulated the magnitude of the AR transcriptome. RARγ downregulation was explained by increased miR-96 in PCa cell and mouse models, and TCGA PCa cohorts. Biochemical approaches confirmed that miR-96 directly regulated RARγ expression and function. Capture of the miR-96 targetome by biotin-miR-96 identified that RARγ and a number of RARγ interacting co-factors including TACC1 were all targeted by miR-96, and expression of these genes were prominently altered, positively and negatively, in the TCGA-PRAD cohort. Differential gene expression analyses between tumors in the TCGA-PRAD cohort with lower quartile expression levels of RARG and TACC1 and upper quartile miR-96, compared to the reverse, identified a gene network including several RARγ target genes (e.g., SOX15) that significantly associated with worse disease-free survival (hazard ratio 2.23, 95% CI 1.58 to 2.88, p = 0.015). In summary, miR-96 targets a RARγ network to govern AR signaling, PCa progression and disease outcome.