Search

In the last decades, thousands of common genetic variants have been associated with human diseases by genome-wide association studies (GWAS). However, the functional interpretation of GWAS hits is usually nontrivial, especially because most of them lay outside the coding genome. These noncoding variants presumably exert their effect by altering gene expression levels; therefore, expression quantitative trait loci (eQTL) mapping analyses represent an important step in understanding their functional relevance and identifying the target genes. Here we describe an alternative strategy for the detection of eQTL that takes into account the combined effect of genetic variants within regulatory regions and leverages the idea that changes in gene expression often are the consequence of the alteration of transcription factor (TF) binding.

Alternative polyadenylation (APA) is a process that changes the posttranscriptional regulation and translation potential of mRNAs via addition or deletion of 3' untranslated region (3' UTR) sequences. To identify posttranscriptional-regulatory events affected by APA in breast tumors, tumor datasets were analyzed for recurrent APA events. Motif mapping of the changed 3' UTR regions found that APA-mediated removal of Pumilio regulatory elements (PRE) was unusually common. Breast tumor subtype-specific APA profiling identified triple-negative breast tumors as having the highest levels of APA. To determine the frequency of these events, an independent cohort of triple-negative breast tumors and normal breast tissue was analyzed for APA. APA-mediated shortening of NRAS and c-JUN was seen frequently, and this correlated with changes in the expression of downstream targets. mRNA stability and luciferase assays demonstrated APA-dependent alterations in RNA and protein levels of affected candidate genes. Examination of clinical parameters of these tumors found those with APA of NRAS and c-JUN to be smaller and less proliferative, but more invasive than non-APA tumors. RT-PCR profiling identified elevated levels of polyadenylation factor CSTF3 in tumors with APA. Overexpression of CSTF3 was common in triple-negative breast cancer cell lines, and elevated CSTF3 levels were sufficient to induce APA of NRAS and c-JUN. Our results support the hypothesis that PRE-containing mRNAs are disproportionately affected by APA, primarily due to high sequence similarity in the motifs utilized by polyadenylation machinery and the PUM complex. Cancer Res; 76(24); 7231-41. ©2016 AACR., (©2016 American Association for Cancer Research.)

Expression of the SS18/SYT-SSX fusion protein is believed to underlie the pathogenesis of synovial sarcoma (SS). Recent evidence suggests that deregulation of the Wnt pathway may play an important role in SS but the mechanisms whereby SS18-SSX might affect Wnt signaling remain to be elucidated. Here, we show that SS18/SSX tightly regulates the elevated expression of the key Wnt target AXIN2 in primary SS. SS18-SSX is shown to interact with TCF/LEF, TLE and HDAC but not β-catenin in vivo and to induce Wnt target gene expression by forming a complex containing promoter-bound TCF/LEF and HDAC but lacking β-catenin. Our observations provide a tumor-specific mechanistic basis for Wnt target gene induction in SS that can occur in the absence of Wnt ligand stimulation.

Transcription factors regulate gene expression by binding regulatory DNA. Understanding the rules governing such binding is an essential step in describing the network of regulatory interactions, and its pathological alterations. We show that describing regulatory regions in terms of their profile of total binding affinities for transcription factors leads to increased predictive power compared to methods based on the identification of discrete binding sites. This applies both to the prediction of transcription factor binding as revealed by ChIP-seq experiments and to the prediction of gene expression through RNA-seq. Further significant improvements in predictive power are obtained when regulatory regions are defined based on chromatin states inferred from histone modification data.

OTX proteins, homologs of the Drosophila orthodenticle (Otd), are important for the morphogenesis of the neuroectoderm, and for the central nervous system formation. OTX1 and OTX2 are important for the cochlea and macula development, indeed when OTX1 is knocked down, these organs undergo developmental failure. Moreover OTX2 transfection revert this effect in OTX1(-/-) mice. The TA isoform of TP63, involved in Notch regulation pathway, has a critical function in the cochlear neuroepithelium differentiation. TAp63 positively regulates Hes5 and Atoh1 transcription. This pathway has been also demonstrated in p63(-/-) mice, and in patients p63 mutated, affected by Ectodermal Dysplasia (ED, OMIM 129810). These patients are affected by mild sensorineural deafness, most likely related to the mutation in p63 gene impairing the Notch pathway. We demonstrated the role of OTX2 on TAp63 regulation necessary for the correct formation of macular neuroepithelium and we confirmed the impairment of vestibular function caused by p63 mutations. Although the abnormalities found in our patient were still at a subclinical extent, aging could exacerbate this impairment and cause a decrease in quality of life.

Changes in gene regulation are believed to play an important role in the evolution of animals. It has been suggested that changes in cis-regulatory regions are responsible for many or most of the anatomical and behavioral differences between humans and apes. However, the study of the evolution of cis-regulatory regions is made problematic by the degeneracy of transcription factor (TF) binding sites and the shuffling of their positions. In this work, we use the predicted total affinity of a promoter for a large collection of TFs as the basis for studying the evolution of cis-regulatory regions in mammals. We introduce the human specificity of a promoter, measuring the divergence between the affinity profile of a human promoter and its orthologous promoters in other mammals. The promoters of genes involved in functional categories such as neural processes and signal transduction, among others, have higher human specificity compared with the rest of the genome. Clustering of the human-specific affinities (HSAs) of neural genes reveals patterns of promoter evolution associated with functional categories such as synaptic transmission and brain development and to diseases such as bipolar disorder and autism.

Overexpression of the polycomb group protein enhancer of zeste homologue 2 (EZH2) occurs in diverse malignancies, including prostate cancer, breast cancer, and glioblastoma multiforme (GBM). Based on its ability to modulate transcription of key genes implicated in cell cycle control, DNA repair, and cell differentiation, EZH2 is believed to play a crucial role in tissue-specific stem cell maintenance and tumor development. Here, we show that targeted pharmacologic disruption of EZH2 by the S-adenosylhomocysteine hydrolase inhibitor 3-deazaneplanocin A (DZNep), or its specific downregulation by short hairpin RNA (shRNA), strongly impairs GBM cancer stem cell (CSC) self-renewal in vitro and tumor-initiating capacity in vivo. Using genome-wide expression analysis of DZNep-treated GBM CSCs, we found the expression of c-myc, recently reported to be essential for GBM CSCs, to be strongly repressed upon EZH2 depletion. Specific shRNA-mediated downregulation of EZH2 in combination with chromatin immunoprecipitation experiments revealed that c-myc is a direct target of EZH2 in GBM CSCs. Taken together, our observations provide evidence that direct transcriptional regulation of c-myc by EZH2 may constitute a novel mechanism underlying GBM CSC maintenance and suggest that EZH2 may be a valuable new therapeutic target for GBM management.

The identification of direct targets of transcription factors is a key problem in the study of gene regulatory networks. However, the use of high throughput experimental methods, such as ChIP-chip and ChIP-sequencing, is limited by their high cost and strong dependence on cellular type and context. We developed a computational method for the genome-wide identification of functional transcription factor binding sites based on positional weight matrices, comparative genomics, and gene expression profiling. The method was applied to Stat3, a transcription factor playing crucial roles in inflammation, immunity and oncogenesis, and able to induce distinct subsets of target genes in different cell types or conditions. A newly generated positional weight matrix enabled us to assign affinity scores of high specificity, as measured by EMSA competition assays. Phylogenetic conservation with 7 vertebrate species was used to select the binding sites most likely to be functional. Validation was carried out on predicted sites within genes identified as differentially expressed in the presence or absence of Stat3 by microarray analysis. Twelve of the fourteen sites tested were bound by Stat3 in vivo, as assessed by Chromatin Immunoprecipitation, allowing us to identify 9 Stat3 transcriptional targets. Given its high validation rate, and the availability of large transcription factor-dependent gene expression datasets obtained under diverse experimental conditions, our approach appears to be a valid alternative to high-throughput experimental assays for the discovery of novel direct targets of transcription factors.

Background: The EWS-FLI-1 fusion protein is associated with 85-90% of Ewing's sarcoma family tumors (ESFT), the remaining 10-15% of cases expressing chimeric genes encoding EWS or FUS fused to one of several ets transcription factor family members, including ERG-1, FEV, ETV1 and ETV6. ESFT are dependent on insulin-like growth factor-1 (IGF-1) for growth and survival and recent evidence suggests that mesenchymal progenitor/stem cells constitute a candidate ESFT origin., Methodology/principal Findings: To address the functional relatedness between ESFT-associated fusion proteins, we compared mouse progenitor cell (MPC) permissiveness for EWS-FLI-1, EWS-ERG and FUS-ERG expression and assessed the corresponding expression profile changes. Whereas all MPC isolates tested could stably express EWS-FLI-1, only some sustained stable EWS-ERG expression and none could express FUS-ERG for more than 3-5 days. Only 14% and 4% of the total number of genes that were respectively induced and repressed in MPCs by the three fusion proteins were shared. However, all three fusion proteins, but neither FLI-1 nor ERG-1 alone, activated the IGF1 promoter and induced IGF1 expression., Conclusion/significance: Whereas expression of different ESFT-associated fusion proteins may require distinct cellular microenvironments and induce transcriptome changes of limited similarity, IGF1 induction may provide one common mechanism for their implication in ESFT pathogenesis.

Background: Transcriptional regulation is a key mechanism in the functioning of the cell, and is mostly effected through transcription factors binding to specific recognition motifs located upstream of the coding region of the regulated gene. The computational identification of such motifs is made easier by the fact that they often appear several times in the upstream region of the regulated genes, so that the number of occurrences of relevant motifs is often significantly larger than expected by pure chance., Results: To exploit this fact, we construct sets of genes characterized by the statistical overrepresentation of a certain motif in their upstream regions. Then we study the functional characterization of these sets by analyzing their annotation to Gene Ontology terms. For the sets showing a statistically significant specific functional characterization, we conjecture that the upstream motif characterizing the set is a binding site for a transcription factor involved in the regulation of the genes in the set., Conclusions: The method we propose is able to identify many known binding sites in S. cerevisiae and new candidate targets of regulation by known transcription factors. Its application to less well studied organisms is likely to be valuable in the exploration of their regulatory interaction network.

In pluripotent stem cells, there is increasing evidence for crosstalk between post-transcriptional and transcriptional networks, offering multifold steps at which pluripotency can be controlled. In addition to well-studied transcription factors, chromatin modifiers and miRNAs, RNA-binding proteins are emerging as fundamental players in pluripotency regulation. Here, we report a new role for the RNA-binding protein ESRP1 in the control of pluripotency. Knockdown of Esrp1 in mouse embryonic stem cells induces, other than the well-documented epithelial to mesenchymal-like state, also an increase in expression of the core transcription factors Oct4, Nanog and Sox2, thereby enhancing self-renewal of these cells. Esrp1-depleted embryonic stem cells displayed impaired early differentiation in vitro and formed larger teratomas in vivo when compared to control embryonic stem cells. We also show that ESRP1 binds to Oct4 and Sox2 mRNAs and decreases their polysomal loading. ESRP1 thus acts as a physiological regulator of the finely-tuned balance between self-renewal and commitment to a restricted developmental fate. Importantly, both mouse and human epithelial stem cells highly express ESRP1, pinpointing the importance of this RNA-binding protein in stem cell biology. [ABSTRACT FROM AUTHOR]

ABSTRACT Polymorphisms in regulatory DNA regions are believed to play an important role in determining phenotype, including disease, and in providing raw material for evolution. We devised a new pipeline for the systematic identification of functional variation in human regulatory sequences. The algorithm is based on the identification of SNPs leading to significant changes in both the affinity of a regulatory region for transcription factors ( TFs) and the expression in vivo of the regulated gene. We tested the algorithm by identifying SNPs leading to altered regulation by STAT3 in human promoters and introns, and experimentally validated the top-scoring ones, showing that most of the SNPs identified by the algorithm indeed correspond to differential binding of STAT3 and differential induction of the target gene upon stimulation with IL6. Using the same computational approach, we compiled a database of thousands of predicted functional regulatory SNPs for hundreds of human TFs, which we provide as online Supporting Information. We discuss possible applications to the interpretation of noncoding SNPs associated with human diseases. The method we propose and the database of predicted functional cis-regulatory polymorphisms will be useful in future studies of regulatory variation and in particular to interpret the results of past and future genome-wide association studies. [ABSTRACT FROM AUTHOR]

The Dlx and Msx homeodomain transcription factors play important roles in the control of limb development. The combined disruption of Msx1 and Msx2, as well as that of D lx5 and Dlx6, lead to limb patterning defects with anomalies in digit number and shape. Msx1;Msx2 double mutants are characterized by the loss of derivatives of the anterior limb mesoderm which is not observed in either of the simple mutants. Dlx5;Dlx6 double mutants exhibit hindlimb ectrodactyly. While the morphogenetic action of Msx genes seems to involve the BMP molecules, the mode of action of Dlx genes still remains elusive. Here, examining the limb phenotypes of combined Dlx and Msx mutants we reveal a new Dlx-Msx regulatory loop directly involving BMPs. In Msx1;Dlx5;Dlx6 triple mutant mice (TKO), beside the expected ectrodactyly, we also observe the hallmark morphological anomalies of Msx1;Msx2 double mutants suggesting an epistatic role of Dlx5 and Dlx6 over Msx2. In Msx2;Dlx5;Dlx6 TKO mice we only observe an aggravation of the ectrodactyly defect without changes in the number of the individual components of the limb. Using a combination of qPCR, ChIP and bioinformatic analyses, we identify two Dlx/Msx regulatory pathways: 1) in the anterior limb mesoderm a non-cell autonomous Msx-Dlx regulatory loop involves BMP molecules through the AER and 2) in AER cells and, at later stages, in the limb mesoderm the regulation of Msx2 by Dlx5 and Dlx6 occurs also cell autonomously. These data bring new elements to decipher the complex AER-mesoderm dialogue that takes place during limb development and provide clues to understanding the etiology of congenital limb malformations. [ABSTRACT FROM AUTHOR]

Ewing's sarcoma family tumors (ESFT) are the second most common bone malignancy in children and young adults, characterized by unique chromosomal translocations that in 85% of cases lead to expression of the EWS-FLI-1 fusion protein. EWS-FLI-1 functions as an aberrant transcription factor that can both induce and suppress members of its target gene repertoire. We have recently demonstrated that EWS-FLI-1 can alter microRNA (miRNA) expression and that miRNA145 is a direct EWS-FLI-1 target whose suppression is implicated in ESFT development. Here, we use miRNA arrays to compare the global miRNA expression profile of human mesenchymal stem cells (MSC) and ESFT cell lines, and show that ESFT display a distinct miRNA signature that includes induction of the oncogenic miRNA 17-92 cluster and repression of the tumor suppressor let-7 family. We demonstrate that direct repression of let-7a by EWS-FLI-1 participates in the tumorigenic potential of ESFT cells in vivo. The mechanism whereby let-7a expression regulates ESFT growth is shown to be mediated by its target gene HMGA2, as let-7a overexpression and HMGA2 repression both block ESFT cell tumorigenicity. Consistent with these observations, systemic delivery of synthetic let-7a into ESFT-bearing mice restored its expression in tumor cells, decreased HMGA2 expression levels and resulted in ESFT growth inhibition in vivo. Our observations provide evidence that deregulation of let-7a target gene expression participates in ESFT development and identify let-7a as promising new therapeutic target for one of the most aggressive pediatric malignancies. [ABSTRACT FROM AUTHOR]

Background: Different approaches have been developed to dissect the interplay between transcription factors (TFs) and their cis-acting sequences on DNA in order to identify TF target genes. Here we used a combination of computational and experimental approaches to identify novel direct targets of TFAP2A, a key TF for a variety of physiological and pathological cellular processes. Gene expression profiles of HeLa cells either silenced for TFAP2A by RNA interference or not were previously compared and a set of differentially expressed genes was revealed. Results: The regulatory regions of 494 TFAP2A-modulated genes were analyzed for the presence of TFAP2A binding sites, employing the canonical TFAP2A Positional Weight Matrix (PWM) reported in Jaspar http://jaspar.genereg.net/. 264 genes containing at least 2 high score TFAP2A binding sites were identified, showing a central role in "Cellular Movement" and "Cellular Development". In an attempt to identify TFs that could cooperate with TFAP2A, a statistically significant enrichment for SP1 binding sites was found for TFAP2A-activated but not repressed genes. The direct binding of TFAP2A or SP1 to a random subset of TFAP2A-modulated genes was demonstrated by Chromatin ImmunoPrecipitation (ChIP) assay and the TFAP2A-driven regulation of DCBLD2/ESDN/CLCP1 gene studied in details. Conclusions: We proved that our computational approaches applied to microarray selected genes are valid tools to identify functional TF binding sites in gene regulatory regions as confirmed by experimental validations. In addition, we demonstrated a fine-tuned regulation of DCBLD2/ESDN transcription by TFAP2A. [ABSTRACT FROM AUTHOR]

The identification of direct targets of transcription factors is a key problem in the study of gene regulatory networks. However, the use of high throughput experimental methods, such as ChIP-chip and ChIP-sequencing, is limited by their high cost and strong dependence on cellular type and context. We developed a computational method for the genome-wide identification of functional transcription factor binding sites based on positional weight matrices, comparative genomics, and gene expression profiling. The method was applied to Stat3, a transcription factor playing crucial roles in inflammation, immunity and oncogenesis, and able to induce distinct subsets of target genes in different cell types or conditions. A newly generated positional weight matrix enabled us to assign affinity scores of high specificity, as measured by EMSA competition assays. Phylogenetic conservation with 7 vertebrate species was used to select the binding sites most likely to be functional. Validation was carried out on predicted sites within genes identified as differentially expressed in the presence or absence of Stat3 by microarray analysis. Twelve of the fourteen sites tested were bound by Stat3 in vivo, as assessed by Chromatin Immunoprecipitation, allowing us to identify 9 Stat3 transcriptional targets. Given its high validation rate, and the availability of large transcription factor-dependent gene expression datasets obtained under diverse experimental conditions, our approach appears to be a valid alternative to high-throughput experimental assays for the discovery of novel direct targets of transcription factors. [ABSTRACT FROM AUTHOR]

Background: Transcriptional regulation is a key mechanism in the functioning of the cell, and is mostly effected through transcription factors binding to specific recognition motifs located upstream of the coding region of the regulated gene. The computational identification of such motifs is made easier by the fact that they often appear several times in the upstream region of the regulated genes, so that the number of occurrences of relevant motifs is often significantly larger than expected by pure chance. Results: To exploit this fact, we construct sets of genes characterized by the statistical overrepresentation of a certain motif in their upstream regions. Then we study the functional characterization of these sets by analyzing their annotation to Gene Ontology terms. For the sets showing a statistically significant specific functional characterization, we conjecture that the upstream motif characterizing the set is a binding site for a transcription factor involved in the regulation of the genes in the set. Conclusions: The method we propose is able to identify many known binding sites in S. cerevisiae and new candidate targets of regulation by known transcritpion factors. Its application to less well studied organisms is likely to be valuable in the exploration of their regulatory interaction network. [ABSTRACT FROM AUTHOR]

Breast cancer is a heterogeneous disease whose clinical management is very challenging. Although specific molecular features characterize breast cancer subtypes with different prognosis, the identification of specific markers predicting disease outcome within the single subtypes still lags behind. Both the non-canonical Wingless-type MMTV Integration site (WNT) and the Signal Transducer and Activator of Transcription (STAT)3 pathways are often constitutively activated in breast tumors, and both can induce the small GTPase Ras Homolog Family Member U RhoU. Here we show that RhoU transcription can be triggered by both canonical and non-canonical WNT ligands via the activation of c-JUN N-terminal kinase (JNK) and the recruitment of the Specificity Protein 1 (SP1) transcription factor to the RhoU promoter, identifying for the first time SP1 as a JNK-dependent mediator of WNT signaling. RhoU down-regulation by silencing or treatment with JNK, SP1 or STAT3 inhibitors leads to impaired migration and invasion in basal-like MDA-MB-231 and BT-549 cells, suggesting that STAT3 and SP1 can cooperate to induce high RhoU expression and enhance breast cancer cells migration. Moreover, in vivo concomitant binding of STAT3 and SP1 defines a subclass of genes belonging to the non-canonical WNT and the Interleukin (IL)-6/STAT3 pathways and contributing to breast cancer aggressiveness, suggesting the relevance of developing novel targeted therapies combining inhibitors of the STAT3 and WNT pathways or of their downstream mediators. [ABSTRACT FROM AUTHOR]

During neuronal development and maturation, microRNAs (miRs) play diverse functions ranging from early patterning, proliferation and commitment to differentiation, survival, homeostasis, activity and plasticity of more mature and adult neurons. The role of miRs in the differentiation of olfactory receptor neurons (ORNs) is emerging from the conditional inactivation of Dicer in immature ORN, and the depletion of all mature miRs in this system. Here, we identify specific miRs involved in olfactory development, by focusing on mice null for Dlx5 , a homeogene essential for both ORN differentiation and axon guidance and connectivity. Analysis of miR expression in Dlx5 −/− olfactory epithelium pointed to reduced levels of miR-9 , miR-376a and four miRs of the - 200 class in the absence of Dlx5 . To functionally examine the role of these miRs, we depleted miR-9 and miR-200 class in reporter zebrafish embryos and observed delayed ORN differentiation, altered axonal trajectory/targeting, and altered genesis and position of olfactory-associated GnRH neurons, i.e. a phenotype known as Kallmann syndrome in humans. miR-9 and miR-200 -class negatively control Foxg1 mRNA, a fork-head transcription factor essential for development of the olfactory epithelium and of the forebrain, known to maintain progenitors in a stem state. Increased levels of z- foxg1 mRNA resulted in delayed ORN differentiation and altered axon trajectory, in zebrafish embryos. This work describes for the first time the role of specific miR (-9 and -200) in olfactory/GnRH development, and uncovers a Dlx5–Foxg1 regulation whose alteration affects receptor neuron differentiation, axonal targeting, GnRH neuron development, the hallmarks of the Kallmann syndrome. [ABSTRACT FROM AUTHOR]

Genome expansion is believed to be an important driver of the evolution of gene regulation. To investigate the role of a newly arising sequence in rewiring regulatory networks, we estimated the age of each region of the human genome by applying maximum parsimony to genome-wide alignments with 100 vertebrates. We then studied the age distribution of several types of functional regions, with a focus on regulatory elements. The age distribution of regulatory elements reveals the extensive use of newly formed genomic sequence in the evolution of regulatory interactions. Many transcription factors have expanded their repertoire of targets through waves of genomic expansions that can be traced to specific evolutionary times. Repeated elements contributed a major part of such expansion: many classes of such elements are enriched in binding sites of one or a few specific transcription factors, whose binding sites are localized in specific portions of the element and characterized by distinctive motif words. These features suggest that the binding sites were available as soon as the new sequence entered the genome, rather than being created later by accumulation of point mutations. By comparing the age of regulatory regions to the evolutionary shift in expression of nearby genes, we show that rewiring through genome expansion played an important role in shaping human regulatory networks. [ABSTRACT FROM AUTHOR]