The noncoding RNA Mistral activates Hoxa6 and Hoxa7 expression and stem cell differentiation by recruiting MLL1 to chromatin.

The epigenetic activator Mixed lineage leukemia 1 (MLL1) is paramount for embryonic development and hematopoiesis. Here, we demonstrate that the long, noncoding RNA (lncRNA) Mistral (Mira) activates transcription of the homeotic genes Hoxa6 and Hoxa7 in mouse embryonic stem cells (mESC) by recruiting MLL1 to chromatin. The Mira gene is located in the spacer DNA region (SDR) separating Hoxa6 and Hoxa7, transcriptionally silent in mESCs, and activated by retinoic acid. Mira-mediated recruitment of MLL1 to the Mira gene triggers dynamic changes in chromosome conformation, culminating in activation of Hoxa6 and Hoxa7 transcription. Hoxa6 and Hoxa7 activate the expression of genes involved in germ layer specification during mESC differentiation in a cooperative and redundant fashion. Our results connect the lncRNA Mira with the recruitment of MLL1 to target genes and implicate lncRNAs in epigenetic activation of gene expression during vertebrate cell-fate determination.

Statin-induced changes in gene expression in EBV-transformed and native B-cells.

Human lymphoblastoid cell lines (LCLs), generated through Epstein-Barr Virus (EBV) transformation of B-lymphocytes (B-cells), are a commonly used model system for identifying genetic influences on human diseases and on drug responses. We have previously used LCLs to examine the cellular effects of genetic variants that modulate the efficacy of statins, the most prescribed class of cholesterol-lowering drugs used for the prevention and treatment of cardiovascular disease. However, statin-induced gene expression differences observed in LCLs may be influenced by their transformation, and thus differ from those observed in native B-cells. To assess this possibility, we prepared LCLs and purified B-cells from the same donors, and compared mRNA profiles after 24 h incubation with simvastatin (2 µm) or sham buffer. Genes involved in cholesterol metabolism were similarly regulated between the two cell types under both the statin and sham-treated conditions, and the statin-induced changes were significantly correlated. Genes whose expression differed between the native and transformed cells were primarily implicated in cell cycle, apoptosis and alternative splicing. We found that ChIP-seq signals for MYC and EBNA2 (an EBV transcriptional co-activator) were significantly enriched in the promoters of genes up-regulated in the LCLs compared with the B-cells, and could be involved in the regulation of cell cycle and alternative splicing. Taken together, the results support the use of LCLs for the study of statin effects on cholesterol metabolism, but suggest that drug effects on cell cycle, apoptosis and alternative splicing may be affected by EBV transformation. This dataset is now uploaded to GEO at the accession number GSE51444.

Identification of a binding motif specific to HNF4 by comparative analysis of multiple nuclear receptors.

Nuclear receptors (NRs) regulate gene expression by binding specific DNA sequences consisting of AG[G/T]TCA or AGAACA half site motifs in a variety of configurations. However, those motifs/configurations alone do not adequately explain the diversity of NR function in vivo. Here, a systematic examination of DNA binding specificity by protein-binding microarrays (PBMs) of three closely related human NRs--HNF4α, retinoid X receptor alpha (RXRα) and COUPTF2--reveals an HNF4-specific binding motif (H4-SBM), xxxxCAAAGTCCA, as well as a previously unrecognized polarity in the classical DR1 motif (AGGTCAxAGGTCA) for HNF4α, RXRα and COUPTF2 homodimers. ChIP-seq data indicate that the H4-SBM is uniquely bound by HNF4α but not 10 other NRs in vivo, while NRs PXR, FXRα, Rev-Erbα appear to bind adjacent to H4-SBMs. HNF4-specific DNA recognition and transactivation are mediated by residues Asp69 and Arg76 in the DNA-binding domain; this combination of amino acids is unique to HNF4 among all human NRs. Expression profiling and ChIP data predict ≈ 100 new human HNF4α target genes with an H4-SBM site, including several Co-enzyme A-related genes and genes with links to disease. These results provide important new insights into NR DNA binding.

Nuclear receptor HNF4α binding sequences are widespread in Alu repeats.

BACKGROUND: Alu repeats, which account for ~10% of the human genome, were originally considered to be junk DNA. Recent studies, however, suggest that they may contain transcription factor binding sites and hence possibly play a role in regulating gene expression. RESULTS: Here, we show that binding sites for a highly conserved member of the nuclear receptor superfamily of ligand-dependent transcription factors, hepatocyte nuclear factor 4alpha (HNF4α, NR2A1), are highly prevalent in Alu repeats. We employ high throughput protein binding microarrays (PBMs) to show that HNF4α binds > 66 unique sequences in Alu repeats that are present in ~1.2 million locations in the human genome. We use chromatin immunoprecipitation (ChIP) to demonstrate that HNF4α binds Alu elements in the promoters of target genes (ABCC3, APOA4, APOM, ATPIF1, CANX, FEMT1A, GSTM4, IL32, IP6K2, PRLR, PRODH2, SOCS2, TTR) and luciferase assays to show that at least some of those Alu elements can modulate HNF4α-mediated transactivation in vivo (APOM, PRODH2, TTR, APOA4). HNF4α-Alu elements are enriched in promoters of genes involved in RNA processing and a sizeable fraction are in regions of accessible chromatin. Comparative genomics analysis suggests that there may have been a gain in HNF4α binding sites in Alu elements during evolution and that non Alu repeats, such as Tiggers, also contain HNF4α sites. CONCLUSIONS: Our findings suggest that HNF4α, in addition to regulating gene expression via high affinity binding sites, may also modulate transcription via low affinity sites in Alu repeats.

Integrated approach for the identification of human hepatocyte nuclear factor 4alpha target genes using protein binding microarrays.

UNLABELLED: Hepatocyte nuclear factor 4 alpha (HNF4alpha), a member of the nuclear receptor superfamily, is essential for liver function and is linked to several diseases including diabetes, hemophilia, atherosclerosis, and hepatitis. Although many DNA response elements and target genes have been identified for HNF4alpha, the complete repertoire of binding sites and target genes in the human genome is unknown. Here, we adapt protein binding microarrays (PBMs) to examine the DNA-binding characteristics of two HNF4alpha species (rat and human) and isoforms (HNF4alpha2 and HNF4alpha8) in a high-throughput fashion. We identified approximately 1400 new binding sequences and used this dataset to successfully train a Support Vector Machine (SVM) model that predicts an additional approximately 10,000 unique HNF4alpha-binding sequences; we also identify new rules for HNF4alpha DNA binding. We performed expression profiling of an HNF4alpha RNA interference knockdown in HepG2 cells and compared the results to a search of the promoters of all human genes with the PBM and SVM models, as well as published genome-wide location analysis. Using this integrated approach, we identified approximately 240 new direct HNF4alpha human target genes, including new functional categories of genes not typically associated with HNF4alpha, such as cell cycle, immune function, apoptosis, stress response, and other cancer-related genes. CONCLUSION: We report the first use of PBMs with a full-length liver-enriched transcription factor and greatly expand the repertoire of HNF4alpha-binding sequences and target genes, thereby identifying new functions for HNF4alpha. We also establish a web-based tool, HNF4 Motif Finder, that can be used to identify potential HNF4alpha-binding sites in any sequence.

Prevalence of the initiator over the TATA box in human and yeast genes and identification of DNA motifs enriched in human TATA-less core promoters.

The core promoter of eukaryotic genes is the minimal DNA region that recruits the basal transcription machinery to direct efficient and accurate transcription initiation. The fraction of human and yeast genes that contain specific core promoter elements such as the TATA box and the initiator (INR) remains unclear and core promoter motifs specific for TATA-less genes remain to be identified. Here, we present genome-scale computational analyses indicating that approximately 76% of human core promoters lack TATA-like elements, have a high GC content, and are enriched in Sp1-binding sites. We further identify two motifs - M3 (SCGGAAGY) and M22 (TGCGCANK) - that occur preferentially in human TATA-less core promoters. About 24% of human genes have a TATA-like element and their promoters are generally AT-rich; however, only approximately 10% of these TATA-containing promoters have the canonical TATA box (TATAWAWR). In contrast, approximately 46% of human core promoters contain the consensus INR (YYANWYY) and approximately 30% are INR-containing TATA-less genes. Significantly, approximately 46% of human promoters lack both TATA-like and consensus INR elements. Surprisingly, mammalian-type INR sequences are present - and tend to cluster - in the transcription start site (TSS) region of approximately 40% of yeast core promoters and the frequency of specific core promoter types appears to be conserved in yeast and human genomes. Gene Ontology analyses reveal that TATA-less genes in humans, as in yeast, are frequently involved in basic "housekeeping" processes, while TATA-containing genes are more often highly regulated, such as by biotic or stress stimuli. These results reveal unexpected similarities in the occurrence of specific core promoter types and in their associated biological processes in yeast and humans and point to novel vertebrate-specific DNA motifs that might play a selective role in TATA-independent transcription.

Sequence-specific DNA binding by MYC/MAX to low-affinity non-E-box motifs.

The MYC oncoprotein regulates transcription of a large fraction of the genome as an obligatory heterodimer with the transcription factor MAX. The MYC:MAX heterodimer and MAX:MAX homodimer (hereafter MYC/MAX) bind Enhancer box (E-box) DNA elements (CANNTG) and have the greatest affinity for the canonical MYC E-box (CME) CACGTG. However, MYC:MAX also recognizes E-box variants and was reported to bind DNA in a "non-specific" fashion in vitro and in vivo. Here, in order to identify potential additional non-canonical binding sites for MYC/MAX, we employed high throughput in vitro protein-binding microarrays, along with electrophoretic mobility-shift assays and bioinformatic analyses of MYC-bound genomic loci in vivo. We identified all hexameric motifs preferentially bound by MYC/MAX in vitro, which include the low-affinity non-E-box sequence AACGTT, and found that the vast majority (87%) of MYC-bound genomic sites in a human B cell line contain at least one of the top 21 motifs bound by MYC:MAX in vitro. We further show that high MYC/MAX concentrations are needed for specific binding to the low-affinity sequence AACGTT in vitro and that elevated MYC levels in vivo more markedly increase the occupancy of AACGTT sites relative to CME sites, especially at distal intergenic and intragenic loci. Hence, MYC binds diverse DNA motifs with a broad range of affinities in a sequence-specific and dose-dependent manner, suggesting that MYC overexpression has more selective effects on the tumor transcriptome than previously thought.

Opposing roles of nuclear receptor HNF4α isoforms in colitis and colitis-associated colon cancer.

HNF4α has been implicated in colitis and colon cancer in humans but the role of the different HNF4α isoforms expressed from the two different promoters (P1 and P2) active in the colon is not clear. Here, we show that P1-HNF4α is expressed primarily in the differentiated compartment of the mouse colonic crypt and P2-HNF4α in the proliferative compartment. Exon swap mice that express only P1- or only P2-HNF4α have different colonic gene expression profiles, interacting proteins, cellular migration, ion transport and epithelial barrier function. The mice also exhibit altered susceptibilities to experimental colitis (DSS) and colitis-associated colon cancer (AOM+DSS). When P2-HNF4α-only mice (which have elevated levels of the cytokine resistin-like ß, RELMß, and are extremely sensitive to DSS) are crossed with Retnlb(-/-) mice, they are rescued from mortality. Furthermore, P2-HNF4α binds and preferentially activates the RELMß promoter. In summary, HNF4α isoforms perform non-redundant functions in the colon under conditions of stress, underscoring the importance of tracking them both in colitis and colon cancer.

Differential Effects of Hepatocyte Nuclear Factor 4α Isoforms on Tumor Growth and T-Cell Factor 4/AP-1 Interactions in Human Colorectal Cancer Cells.

The nuclear receptor hepatocyte nuclear factor 4α (HNF4α) is tumor suppressive in the liver but amplified in colon cancer, suggesting that it also might be oncogenic. To investigate whether this discrepancy is due to different HNF4α isoforms derived from its two promoters (P1 and P2), we generated Tet-On-inducible human colon cancer (HCT116) cell lines that express either the P1-driven (HNF4α2) or P2-driven (HNF4α8) isoform and analyzed them for tumor growth and global changes in gene expression (transcriptome sequencing [RNA-seq] and chromatin immunoprecipitation sequencing [ChIP-seq]). The results show that while HNF4α2 acts as a tumor suppressor in the HCT116 tumor xenograft model, HNF4α8 does not. Each isoform regulates the expression of distinct sets of genes and recruits, colocalizes, and competes in a distinct fashion with the Wnt/ß-catenin mediator T-cell factor 4 (TCF4) at CTTTG motifs as well as at AP-1 motifs (TGAXTCA). Protein binding microarrays (PBMs) show that HNF4α and TCF4 share some but not all binding motifs and that single nucleotide polymorphisms (SNPs) in sites bound by both HNF4α and TCF4 can alter binding affinity in vitro, suggesting that they could play a role in cancer susceptibility in vivo. Thus, the HNF4α isoforms play distinct roles in colon cancer, which could be due to differential interactions with the Wnt/ß-catenin/TCF4 and AP-1 pathways.

Prediction of LDL cholesterol response to statin using transcriptomic and genetic variation.

BACKGROUND: Statins are widely prescribed for lowering LDL-cholesterol (LDLC) levels and risk of cardiovascular disease. There is, however, substantial inter-individual variation in the magnitude of statin-induced LDLC reduction. To date, analysis of individual DNA sequence variants has explained only a small proportion of this variability. The present study was aimed at assessing whether transcriptomic analyses could be used to identify additional genetic contributions to inter-individual differences in statin efficacy. RESULTS: Using expression array data from immortalized lymphoblastoid cell lines derived from 372 participants of the Cholesterol and Pharmacogenetics clinical trial, we identify 100 signature genes differentiating high versus low statin responders. A radial-basis support vector machine prediction model of these signature genes explains 12.3% of the variance in statin-mediated LDLC change. Addition of SNPs either associated with expression levels of the signature genes (eQTLs) or previously reported to be associated with statin response in genome-wide association studies results in a combined model that predicts 15.0% of the variance. Notably, a model of the signature gene associated eQTLs alone explains up to 17.2% of the variance in the tails of a separate subset of the Cholesterol and Pharmacogenetics population. Furthermore, using a support vector machine classification model, we classify the most extreme 15% of high and low responders with high accuracy. CONCLUSIONS: These results demonstrate that transcriptomic information can explain a substantial proportion of the variance in LDLC response to statin treatment, and suggest that this may provide a framework for identifying novel pathways that influence cholesterol metabolism.

The transcription factor encyclopedia.

Here we present the Transcription Factor Encyclopedia (TFe), a new web-based compendium of mini review articles on transcription factors (TFs) that is founded on the principles of open access and collaboration. Our consortium of over 100 researchers has collectively contributed over 130 mini review articles on pertinent human, mouse and rat TFs. Notable features of the TFe website include a high-quality PDF generator and web API for programmatic data retrieval. TFe aims to rapidly educate scientists about the TFs they encounter through the delivery of succinct summaries written and vetted by experts in the field. TFe is available at http://www.cisreg.ca/tfe.

Identification of an endogenous ligand bound to a native orphan nuclear receptor.

Orphan nuclear receptors have been instrumental in identifying novel signaling pathways and therapeutic targets. However, identification of ligands for these receptors has often been based on random compound screens or other biased approaches. As a result, it remains unclear in many cases if the reported ligands are the true endogenous ligands,--i.e., the ligand that is bound to the receptor in an unperturbed in vivo setting. Technical limitations have limited our ability to identify ligands based on this rigorous definition. The orphan receptor hepatocyte nuclear factor 4 alpha (HNF4alpha) is a key regulator of many metabolic pathways and linked to several diseases including diabetes, atherosclerosis, hemophilia and cancer. Here we utilize an affinity isolation/mass-spectrometry (AIMS) approach to demonstrate that HNF4alpha is selectively occupied by linoleic acid (LA, C18:2omega6) in mammalian cells and in the liver of fed mice. Receptor occupancy is dramatically reduced in the fasted state and in a receptor carrying a mutation derived from patients with Maturity Onset Diabetes of the Young 1 (MODY1). Interestingly, however, ligand occupancy does not appear to have a significant effect on HNF4alpha transcriptional activity, as evidenced by genome-wide expression profiling in cells derived from human colon. We also use AIMS to show that LA binding is reversible in intact cells, indicating that HNF4alpha could be a viable drug target. This study establishes a general method to identify true endogenous ligands for nuclear receptors (and other lipid binding proteins), independent of transcriptional function, and to track in vivo receptor occupancy under physiologically relevant conditions.

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data.

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.