A PDX1 cistrome and single-cell transcriptome resource of the developing pancreas.

Pancreatic and duodenal homeobox 1 (PDX1) is crucial for pancreas organogenesis, yet the dynamic changes in PDX1 binding in human or mouse developing pancreas have not been examined. To address this knowledge gap, we performed PDX1 ChIP-seq and single-cell RNA-seq using fetal human pancreata. We integrated our datasets with published datasets and revealed the dynamics of PDX1 binding and potential cell lineage-specific PDX1-bound genes in the pancreas from fetal to adult stages. We identified a core set of developmentally conserved PDX1-bound genes that reveal the broad multifaceted role of PDX1 in pancreas development. Despite the well-known dramatic changes in PDX1 function and expression, we found that PDX1-bound genes are largely conserved from embryonic to adult stages. This points towards a dual role of PDX1 in regulating the expression of its targets at different ages, dependent on other functionally congruent or directly interacting partners. We also showed that PDX1 binding is largely conserved in mouse pancreas. Together, our study reveals PDX1 targets in the developing pancreas in vivo and provides an essential resource for future studies on pancreas development.

A comparative analysis of ENCODE and Cistrome in the context of TF binding signal.

BACKGROUND: With the rise of publicly available genomic data repositories, it is now common for scientists to rely on computational models and preprocessed data, either as control or to discover new knowledge. However, different repositories adhere to the different principles and guidelines, and data processing plays a significant role in the quality of the resulting datasets. Two popular repositories for transcription factor binding sites data - ENCODE and Cistrome - process the same biological samples in alternative ways, and their results are not always consistent. Moreover, the output format of the processing (BED narrowPeak) exposes a feature, the signalValue, which is seldom used in consistency checks, but can offer valuable insight on the quality of the data. RESULTS: We provide evidence that data points with high signalValue(s) (top 25% of values) are more likely to be consistent between ENCODE and Cistrome in human cell lines K562, GM12878, and HepG2. In addition, we show that filtering according to said high values improves the quality of predictions for a machine learning algorithm that detects transcription factor interactions based only on positional information. Finally, we provide a set of practices and guidelines, based on the signalValue feature, for scientists who wish to compare and merge narrowPeaks from ENCODE and Cistrome. CONCLUSIONS: The signalValue feature is an informative feature that can be effectively used to highlight consistent areas of overlap between different sources of TF binding sites that expose it. Its applicability extends to downstream to positional machine learning algorithms, making it a powerful tool for performance tweaking and data aggregation.

SMAD4 target genes are part of a transcriptional network that integrates the response to BMP and SHH signaling during early limb bud patterning.

SMAD4 regulates gene expression in response to BMP and TGFß signal transduction, and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we identify the SMAD4 target genes in mouse limb buds using an epitope-tagged Smad4 allele for ChIP-seq analysis in combination with transcription profiling. This analysis shows that SMAD4 predominantly mediates BMP signal transduction during early limb bud development. Unexpectedly, the expression of cholesterol biosynthesis enzymes is precociously downregulated and intracellular cholesterol levels are reduced in Smad4-deficient limb bud mesenchymal progenitors. Most importantly, our analysis reveals a predominant function of SMAD4 in upregulating target genes in the anterior limb bud mesenchyme. Analysis of differentially expressed genes shared between Smad4- and Shh-deficient limb buds corroborates this function of SMAD4 and also reveals the repressive effect of SMAD4 on posterior genes that are upregulated in response to SHH signaling. This analysis uncovers opposing trans-regulatory inputs from SHH- and SMAD4-mediated BMP signal transduction on anterior and posterior gene expression during the digit patterning and outgrowth in early limb buds.

Modulation of early gene expression responses to water deprivation stress by the E3 ubiquitin ligase ATL80: implications for retrograde signaling interplay.

BACKGROUND: Primary response genes play a pivotal role in translating short-lived stress signals into sustained adaptive responses. In this study, we investigated the involvement of ATL80, an E3 ubiquitin ligase, in the dynamics of gene expression following water deprivation stress. We observed that ATL80 is rapidly activated within minutes of water deprivation stress perception, reaching peak expression around 60 min before gradually declining. ATL80, despite its post-translational regulation role, emerged as a key player in modulating early gene expression responses to water deprivation stress. RESULTS: The impact of ATL80 on gene expression was assessed using a time-course microarray analysis (0, 15, 30, 60, and 120 min), revealing a burst of differentially expressed genes, many of which were associated with various stress responses. In addition, the diversity of early modulation of gene expression in response to water deprivation stress was significantly abolished in the atl80 mutant compared to wild-type plants. A subset of 73 genes that exhibited a similar expression pattern to ATL80 was identified. Among them, several are linked to stress responses, including ERF/AP2 and WRKY transcription factors, calcium signaling genes, MAP kinases, and signaling peptides. Promoter analysis predicts enrichment of binding sites for CAMTA1 and CAMTA5, which are known regulators of rapid stress responses. Furthermore, we have identified a group of differentially expressed ERF/AP2 transcription factors, proteins associated with folding and refolding, as well as pinpointed core module genes which are known to play roles in retrograde signaling pathways that cross-referenced with the early ATL80 transcriptome. CONCLUSIONS: Based on these findings, we propose that ATL80 may target one or more components within the retrograde signaling pathways for degradation. In essence, ATL80 serves as a bridge connecting these signaling pathways and effectively functions as an alarm signal.

Modulation of androgen receptor DNA binding activity through direct interaction with the ETS transcription factor ERG.

The androgen receptor (AR) is a type I nuclear hormone receptor and the primary drug target in prostate cancer due to its role as a lineage survival factor in prostate luminal epithelium. In prostate cancer, the AR cistrome is reprogrammed relative to normal prostate epithelium and particularly in cancers driven by oncogenic ETS fusion genes. The molecular basis for this change has remained elusive. Using purified proteins, we report a minimal cell-free system that demonstrates interdomain cooperativity between the ligand (LBD) and DNA binding domains (DBD) of AR, and its autoinhibition by the N terminus of AR. Furthermore, we identify ERG as a cofactor that activates AR's ability to bind DNA in both high and lower affinity contexts through direct interaction within a newly identified AR-interacting motif (AIM) in the ETS domain, independent of ERG's own DNA binding ability. Finally, we present evidence that this interaction is conserved among ETS factors whose expression is altered in prostate cancer. Our work highlights, at a biochemical level, how tumor-initiating ETS translocations result in reprogramming of the AR cistrome.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR.

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Genomic Insights into Non-steroidal Nuclear Receptors in Prostate and Breast Cancer.

Alterations in transcriptional programs are a fundamental feature of prostate (PCa) and breast cancer (BrCa), and frequently target the actions of the principal steroidal nuclear receptors (NRs), namely the androgen receptor (AR) and the estrogen receptor alpha (ERα), respectively. Indeed, the functions of AR and ERα are central to both prostate and mammary gland biology. The genomic interactions of these NRs become highly distorted in part by changing how they functionally interact with a cohort of non-steroidal Type II NRs, which are by contrast relatively understudied compared to their steroidal cousins. For example, the AR cistrome overlaps with cistromes of different Type II NRs, which suggests a high potential for integrated NR functions to tailor transcriptional signals. Over recent years the cistromes of these Type II NRs, including HNF4s, RARs, PPARs and VDR, have been studied in PCa and BrCa revealing convergence and functional consequences, and are reviewed in the current chapter.

Prostate Cancer Epigenetic Plasticity and Enhancer Heterogeneity: Molecular Causes, Consequences and Clinical Implications.

Prostate cancer (PCa) proliferation is dictated by androgen receptor (AR) signaling, which regulates gene expression through cis-regulatory regions including proximal and distal enhancers. The repertoire of AR interactions at enhancers is dependent on tissue and cellular contexts and thus shape a spectrum of phenotypes through such epigenetic heterogeneity. Moreover, PCa is a multifocal disease and displays a high degree of intra- and inter-tumor heterogeneity, adding to the phenotypic complexity. It is increasingly becoming clear that PCa may be considered an epigenetic disease caused by various molecular causes with profound consequences and clinical implications which are underpinned by enhancer interaction heterogeneity.In this review, we provide a detailed overview of molecular interactors that affect prostate cancer epigenetic heterogeneity, such as coding and non-coding somatic variants, large scale structural variations, pioneer factor binding at enhancers and various contexts that influence enhancer engagement heterogeneity in PCa development and progression. Finally, we explore how the vast heterogeneity in epigenetic profiles identified in recent omics studies results in distinct genomic subtypes which predict disease progression and thus offer opportunities in biomarker discovery and further personalizing cancer treatment. As such, heterogeneous enhancer interactions take center stage in elucidating mechanisms of prostate cancer progression, patient prognostication, therapy discovery and overcoming acquired treatment resistance.

Deciphering essential cistromes using genome-wide CRISPR screens.

Although millions of transcription factor binding sites, or cistromes, have been identified across the human genome, defining which of these sites is functional in a given condition remains challenging. Using CRISPR/Cas9 knockout screens and gene essentiality or fitness as the readout, we systematically investigated the essentiality of over 10,000 FOXA1 and CTCF binding sites in breast and prostate cancer cells. We found that essential FOXA1 binding sites act as enhancers to orchestrate the expression of nearby essential genes through the binding of lineage-specific transcription factors. In contrast, CRISPR screens of the CTCF cistrome revealed 2 classes of essential binding sites. The first class of essential CTCF binding sites act like FOXA1 sites as enhancers to regulate the expression of nearby essential genes, while a second class of essential CTCF binding sites was identified at topologically associated domain (TAD) boundaries and display distinct characteristics. Using regression methods trained on our screening data and public epigenetic profiles, we developed a model to predict essential cis-elements with high accuracy. The model for FOXA1 essentiality correctly predicts noncoding variants associated with cancer risk and progression. Taken together, CRISPR screens of cis-regulatory elements can define the essential cistrome of a given factor and can inform the development of predictive models of cistrome function.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines.

BACKGROUND & AIMS: We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins. METHODS: We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured. RESULTS: We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone. CONCLUSIONS: In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

Dynamic transcriptome, accessible genome, and PGR cistrome profiles in the human myometrium.

The myometrium undergoes structural and functional remodeling during pregnancy. We hypothesize that myometrial genomic elements alter correspondingly in preparation for parturition. Human myometrial tissues from nonpregnant (NP) and term pregnant (TP) human subjects were examined by RNAseq, ATACseq, and PGR ChIPseq assays to profile transcriptome, assessible genome, and PGR occupancy. NP and TP specimens exhibit 2890 differentially expressed genes, reflecting an increase of metabolic, inflammatory, and PDGF signaling, among others, in adaptation to pregnancy. At the epigenome level, patterns of accessible genome change between NP and TP myometrium, leading to the altered enrichment of binding motifs for hormone and muscle regulators such as the progesterone receptor (PGR), Krüppel-like factors, and MEF2A transcription factors. PGR genome occupancy exhibits a significant difference between the two stages of the myometrium, concomitant with distinct transcriptomic profiles including genes such as ENO1, LHDA, and PLCL1 in the glycolytic and calcium signaling pathways. Over-representation of SRF, MYOD, and STAT binding motifs in PGR occupying sites further suggests interactions between PGR and major muscle regulators for myometrial gene expression. In conclusion, changes in accessible genome and PGR occupancy are part of the myometrial remodeling process and may serve as mechanisms to formulate the state-specific transcriptome profiles.

Genes Associated with Calcium Signaling are Involved in Alcohol-Induced Breast Cancer Growth.

BACKGROUND: Alcohol consumption is a risk factor for breast cancer, contributing to up to nearly 23,000 new cases each year. Mechanistic studies show that alcohol increases tumor aggressiveness and metastatic potential, promotes angiogenesis, induces chronic inflammation, and dysregulates RNA polymerase III-related genes. Alcohol has also been shown to affect estrogen signaling in breast cancer, including in our study of the transcriptomic effects of alcohol in breast cancer cells. METHODS: To elucidate mechanisms of action of alcohol in breast cancer, we carried out secondary analyses of our alcohol-responsive transcriptome data using gene ontology and pathway databases and analysis tools and cistromic data analysis of candidate transcription factors which may mediate the transcriptomic alterations. Predicted alcohol-responsive pathways and mechanisms were perturbed and examined experimentally in breast cancer cells. The clinical relevance of identified genes was determined by expression profiles in patient samples and correlation with disease outcomes and alcohol consumption in previously published study cohorts. RESULTS: Gene ontology analysis showed that alcohol alters the expression of many metabolism-related genes, and cistromic data of differentially expressed genes revealed the potential involvement of nuclear factor of activated T cells 3 (NFATC3) in mediating the transcriptomic effects of alcohol. Pathway analysis also predicted regulation of calcium signaling by alcohol in breast cancer cells. Chemical perturbation of this pathway reversed the effect of alcohol on breast cancer cell growth and reduced the elevated cytosolic Ca2+ levels induced by alcohol. Expression levels of alcohol-responsive genes in tumor samples from breast cancer patients are associated with poor disease outcomes. Moreover, expression of some of these genes was altered in breast cancer patients who consumed alcohol previously as compared to those who did not drink. CONCLUSION: Alcohol alters expression of genes that regulate intracellular calcium levels and downstream signaling pathways which drive breast cancer cell proliferation and disease progression.

Gene networks orchestrated by MeGI: a single-factor mechanism underlying sex determination in persimmon.

Separating male and female sex organs is one of the main strategies used to maintain genetic diversity within a species. However, the genetic determinants and their regulatory mechanisms have been identified in only a few species. In dioecious persimmons, the homeodomain transcription factor, MeGI, which is the target of a Y chromosome-encoded small-RNA, OGI, can determine floral sexuality. The basic features of this system are conserved in the monoecious hexaploid Oriental persimmon, in which an additional epigenetic regulation of MeGI determines floral sexuality. The downstream regulatory pathways of MeGI remain uncharacterized. In this study, we examined transcriptomic data for male and female flowers from monoecious persimmon cultivars to unveil the gene networks orchestrated by MeGI. A network visualization and cistrome assessment suggested that class-1 KNOTTED-like homeobox (KNOX)/ovate family protein (OFP)/growth regulating factors (GRFs) and short vegetative phase (SVP) genes mediate the differences in gynoecium and androecium development between male and female flowers, respectively. The expression of these genes is directly controlled by MeGI. The gene networks also suggested that some cytokinin, auxin, and gibberellin signaling genes function cooperatively in the KNOX/OFP/GRF pathway during gynoecium differentiation. Meanwhile, SVP may repress PI expression in developing androecia. Overall, our results suggest that MeGI evolved the ability to promote gynoecium development and suppress androecium development by regulating KNOX/OFP/GRF and SVP expression levels, respectively. These insights may help to clarify the molecular mechanism underlying the production of unisexual flowers, while also elucidating the physiological background enabling a single-factor system to establish dioecy in plants.

Embryonic transcription factor SOX9 drives breast cancer endocrine resistance.

The estrogen receptor (ER) drives the growth of most luminal breast cancers and is the primary target of endocrine therapy. Although ER blockade with drugs such as tamoxifen is very effective, a major clinical limitation is the development of endocrine resistance especially in the setting of metastatic disease. Preclinical and clinical observations suggest that even following the development of endocrine resistance, ER signaling continues to exert a pivotal role in tumor progression in the majority of cases. Through the analysis of the ER cistrome in tamoxifen-resistant breast cancer cells, we have uncovered a role for an RUNX2-ER complex that stimulates the transcription of a set of genes, including most notably the stem cell factor SOX9, that promote proliferation and a metastatic phenotype. We show that up-regulation of SOX9 is sufficient to cause relative endocrine resistance. The gain of SOX9 as an ER-regulated gene associated with tamoxifen resistance was validated in a unique set of clinical samples supporting the need for the development of improved ER antagonists.

NR2F2 Orphan Nuclear Receptor is Involved in Estrogen Receptor Alpha-Mediated Transcriptional Regulation in Luminal A Breast Cancer Cells.

Nuclear Receptor Subfamily 2 Group F Member 2 (NR2F2) is a member of the steroid/thyroid hormone receptor superfamily with a crucial role in organogenesis, angiogenesis, cardiovascular development and tumorigenesis. However, there is limited knowledge about the cistrome and transcriptome of NR2F2 in breast cancer. In this study, we mapped the regulatory mechanism by NR2F2 using functional genomic methods. To investigate the clinical significance of NR2F2 in breast cancer, The Cancer Genome Atlas (TCGA) data were used. These results show that a high NR2F2 is associated with better survival of a specific subset of patients, namely those with luminal A breast cancer. Therefore, genome-wide NR2F2 and estrogen receptor alpha (ERα) binding sites were mapped in luminal A breast cancer cells using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), revealing that most NR2F2 overlap with ERα that are co-occupied by forkhead box A1 (FOXA1) and GATA binding protein 3 (GATA3) in active enhancer regions. NR2F2 overlaps with highly frequent ERα chromatin interactions, which are essential for the formation of ERα-bound super-enhancers. In the process of the transcriptome profiling of NR2F2-depleted breast cancer cells such differentially expressed genes have been identified that are involved in endocrine therapy resistance and are also ERα target genes. Overall, these findings demonstrate that the NR2F2 nuclear receptor has a key role in ERα-mediated transcription and it can offer a potential therapeutic target in patients with luminal A breast cancer.

Exogenous ERα Expression in the Mammary Epithelium Decreases Over Time and Does Not Contribute to p53-Deficient Mammary Tumor Formation in Mice.

Approximately 75% of all breast cancers express the nuclear hormone receptor estrogen receptor α (ERα). However, the majority of mammary tumors from genetically engineered mouse models (GEMMs) are ERα-negative. To model ERα-positive breast cancer in mice, we exogenously introduced expression of mouse and human ERα in an existing GEMM of p53-deficient breast cancer. After initial ERα expression during mammary gland development, expression was reduced or lost in adult glands and p53-deficient mammary tumors. Chromatin immunoprecipitation (ChIP)-sequencing analysis of primary mouse mammary epithelial cells (MMECs) derived from these models, in which expression of the ERα constructs was induced in vitro, confirmed interaction of ERα with the DNA. In human breast and endometrial cancer, and also in healthy breast tissue, DNA binding of ERα is facilitated by the pioneer factor FOXA1. Surprisingly, the ERα binding sites identified in primary MMECs, but also in mouse mammary gland and uterus, showed an high enrichment of ERE motifs, but were devoid of Forkhead motifs. Furthermore, exogenous introduction of FOXA1 and GATA3 in ERα-expressing MMECs was not sufficient to promote ERα-responsiveness of these cells. Together, this suggests that species-specific differences in pioneer factor usage between mouse and human are dictated by the DNA sequence, resulting in ERα-dependencies in mice that are not FOXA1 driven. These species-specific differences in ERα-biology may limit the utility of mice for in vivo modeling of ERα-positive breast cancer.

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.

BACKGROUND: The 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. RESULTS: To 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 (r 2 > 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. CONCLUSIONS: These 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.

1,25-Dihydroxyvitamin D3 Controls a Cohort of Vitamin D Receptor Target Genes in the Proximal Intestine That Is Enriched for Calcium-regulating Components.

1,25-Dihydroxyvitamin D3 (1,25(OH)2D3) plays an integral role in calcium homeostasis in higher organisms through its actions in the intestine, kidney, and skeleton. Interestingly, although several intestinal genes are known to play a contributory role in calcium homeostasis, the entire caste of key components remains to be identified. To examine this issue, Cyp27b1 null mice on either a normal or a high calcium/phosphate-containing rescue diet were treated with vehicle or 1,25(OH)2D3 and evaluated 6 h later. RNA samples from the duodena were then subjected to RNA sequence analysis, and the data were analyzed bioinformatically. 1,25(OH)2D3 altered expression of large collections of genes in animals under either dietary condition. 45 genes were found common to both 1,25(OH)2D3-treated groups and were composed of genes previously linked to intestinal calcium uptake, including S100g, Trpv6, Atp2b1, and Cldn2 as well as others. An additional distinct network of 56 genes was regulated exclusively by diet. We then conducted a ChIP sequence analysis of binding sites for the vitamin D receptor (VDR) across the proximal intestine in vitamin D-sufficient normal mice treated with vehicle or 1,25(OH)2D3. The residual VDR cistrome was composed of 4617 sites, which was increased almost 4-fold following hormone treatment. Interestingly, the majority of the genes regulated by 1,25(OH)2D3 in each diet group as well as those found in common in both groups contained frequent VDR sites that likely regulated their expression. This study revealed a global network of genes in the intestine that both represent direct targets of vitamin D action in mice and are involved in calcium absorption.

A temporary compendium of thyroid hormone target genes in brain.

BACKGROUND: Thyroid hormone controls a number of developmental and physiological processes in the brain by directly acting on gene expression. Transcriptome analyses in rodent identified a number of thyroid hormone regulated genes in several brain areas at different stages. Genome wide analysis of chromatin occupancy in a neural cell line also identified a subset of genes which transcription is likely to be directly regulated by thyroid hormone receptors in neurons. However, the abundance of these data and apparent discrepancies between studies brought some confusion. RESULTS: We present here a meta-analysis of available data to identify recurrent themes in thyroid hormone action in brain cells. This provides a curated list of 734 regulated genes in rodent brain, and highlights a small number of likely direct target genes. Some of these genes are also regulated in amphibians during metamorphosis. This article is part of a Special Issue entitled: Nuclear receptors in animal development.

The RUNX2 cistrome in osteoblasts: characterization, down-regulation following differentiation, and relationship to gene expression.

RUNX2 is a transcription factor that is first expressed in early osteoblast-lineage cells and represents a primary determinant of osteoblastogenesis. While numerous target genes are regulated by RUNX2, little is known of sites on the genome occupied by RUNX2 or of the gene networks that are controlled by these sites. To explore this, we conducted a genome-wide analysis of the RUNX2 cistrome in both pre-osteoblastic MC3T3-E1 cells (POB) and their mature osteoblast progeny (OB), characterized the two cistromes and assessed their relationship to changes in gene expression. We found that although RUNX2 was widely bound to the genome in POB cells, this binding profile was reduced upon differentiation to OBs. Numerous sites were lost upon differentiation, new sites were also gained; many sites remained common to both cell states. Additional features were identified as well including location relative to potential target genes, abundance with respect to single genes, the frequent presence of a consensus TGTGGT RUNX2 binding motif, co-occupancy by C/EBPß and the presence of a typical epigenetic histone enhancer signature. This signature was changed quantitatively following differentiation. While RUNX2 binding sites were associated extensively with adjacent genes, the distal nature of the majority of these sites prevented assessment of whether they represented direct targets of RUNX2 action. Changes in gene expression, however, revealed an abundance of genes that contained RUNX2 binding sites and were regulated in concert. These studies establish a basis for further analysis of the role of RUNX2 activity and its function during osteoblast lineage maturation.