MAP/ERK kinase kinase 1 (MEKK1) mediates transcriptional repression by interacting with polycystic kidney disease-1 (PKD1) promoter-bound p53 tumor suppressor protein.

Mitogen-activated protein kinase (MAPK) cascades regulate a wide variety of cellular processes that ultimately depend on changes in gene expression. We have found a novel mechanism whereby one of the key MAP3 kinases, Mekk1, regulates transcriptional activity through an interaction with p53. The tumor suppressor protein p53 down-regulates a number of genes, including the gene most frequently mutated in autosomal dominant polycystic kidney disease (PKD1). We have discovered that Mekk1 translocates to the nucleus and acts as a co-repressor with p53 to down-regulate PKD1 transcriptional activity. This repression does not require Mekk1 kinase activity, excluding the need for an Mekk1 phosphorylation cascade. However, this PKD1 repression can also be induced by the stress-pathway stimuli, including TNFα, suggesting that Mekk1 activation induces both JNK-dependent and JNK-independent pathways that target the PKD1 gene. An Mekk1-p53 interaction at the PKD1 promoter suggests a new mechanism by which abnormally elevated stress-pathway stimuli might directly down-regulate the PKD1 gene, possibly causing haploinsufficiency and cyst formation.

Functional dissection of the ash2 and ash1 transcriptomes provides insights into the transcriptional basis of wing phenotypes and reveals conserved protein interactions

Background: The trithorax group (trxG) genes absent, small or homeotic discs 1 (ash1) and 2 (ash2) were isolated in a screen for mutants with abnormal imaginal discs. Mutations in either gene cause homeotic transformations but Hox genes are not their only targets. Although analysis of double mutants revealed that ash2 and ash1 mutations enhance each other's phenotypes, suggesting they are functionally related, it was shown that these proteins are subunits of distinct complexes.Results: The analysis of wing imaginal disc transcriptomes from ash2 and ash1 mutants showed that they are highly similar. Functional annotation of regulated genes using Gene Ontology allowed identification of severely affected groups of genes that could be correlated to the wing phenotypes observed. Comparison of the differentially expressed genes with those from other genome-wide analyses revealed similarities between ASH2 and Sin3A, suggesting a putative functional relationship. Coimmunoprecipitation studies and immunolocalization on polytene chromosomes demonstrated that ASH2 and Sin3A interact with HCF (host-cell factor). The results of nucleosome western blots and clonal analysis indicated that ASH2 is necessary for trimethylation of the Lys4 on histone 3 (H3K4).Conclusion: The similarity between the transcriptomes of ash2 and ash1 mutants supports a model in which the two genes act together to maintain stable states of transcription. Like in humans, both ASH2 and Sin3A bind HCF. Finally, the reduction of H3K4 trimethylation in ash2 mutants is the first evidence in Drosophila regarding the molecular function of this trxG gene.

Retinal pigment epithelium protein of 65 kDA gene-linked retinal degeneration is not modulated by chicken acidic leucine-rich epidermal growth factor-like domain containing brain protein/Neuroglycan C/ chondroitin sulfate proteoglycan 5.

PURPOSE: To analyze in vivo the function of chicken acidic leucine-rich epidermal growth factor-like domain containing brain protein/Neuroglycan C (gene symbol: Cspg5) during retinal degeneration in the Rpe65⁻/⁻ mouse model of Leber congenital amaurosis. METHODS: We resorted to mice with targeted deletions in the Cspg5 and retinal pigment epithelium protein of 65 kDa (Rpe65) genes (Cspg5⁻/⁻/Rpe65⁻/⁻). Cone degeneration was assessed with cone-specific peanut agglutinin staining. Transcriptional expression of rhodopsin (Rho), S-opsin (Opn1sw), M-opsin (Opn1mw), rod transducin α subunit (Gnat1), and cone transducin α subunit (Gnat2) genes was assessed with quantitative PCR from 2 weeks to 12 months. The retinal pigment epithelium (RPE) was analyzed at P14 with immunodetection of the retinol-binding protein membrane receptor Stra6. RESULTS: No differences in the progression of retinal degeneration were observed between the Rpe65⁻/⁻ and Cspg5⁻/⁻/Rpe65⁻/⁻ mice. No retinal phenotype was detected in the late postnatal and adult Cspg5⁻/⁻ mice, when compared to the wild-type mice. CONCLUSIONS: Despite the previously reported upregulation of Cspg5 during retinal degeneration in Rpe65⁻/⁻ mice, no protective effect or any involvement of Cspg5 in disease progression was identified.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Up-regulation of macrophage migration inhibitory factor gene and protein expression in glial tumor cells during hypoxic and hypoglycemic stress indicates a critical role for angiogenesis in glioblastoma multiforme.

Glioblastoma multiforme (GBM) is the most malignant variant of human glial tumors. A prominent feature of this tumor is the occurrence of necrosis and vascular proliferation. The regulation of glial neovascularization is still poorly understood and the characterization of factors involved in this process is of major clinical interest. Macrophage migration inhibitory factor (MIF) is a pleiotropic cytokine released by leukocytes and by a variety of cells outside of the immune system. Recent work has shown that MIF may function to regulate cellular differentiation and proliferation in normal and tumor-derived cell lines, and may also contribute to the neovascularization of tumors. Our immunohistological analysis of MIF distribution in GBM tissues revealed the strong MIF protein accumulation in close association with necrotic areas and in tumor cells surrounding blood vessels. In addition, MIF expression was frequently associated with the presence of the tumor-suppressor gene p53. To substantiate the concept that MIF might be involved in the regulation of angiogenesis in GBM, we analyzed the MIF gene and protein expression under hypoxic and hypoglycemic stress conditions in vitro. Northern blot analysis showed a clear increase of MIF mRNA after hypoxia and hypoglycemia. We could also demonstrate that the increase of MIF transcripts on hypoxic stress can be explained by a profound transcriptional activation of the MIF gene. In parallel to the increase of MIF transcripts, we observed a significant rise in extracellular MIF protein on angiogenic stimulation. The data of our preliminary study suggest that the up-regulation of MIF expression during hypoxic and hypoglycemic stress might play a critical role for the neovascularization of glial tumors.

Transcriptional potentiation of the vitellogenin B1 promoter by a combination of both nucleosome assembly and transcription factors: an in vitro dissection.

The Xenopus laevis vitellogenin B1 promoter was assembled into nucleosomes in an oocyte extract. Subsequent RNA polymerase II-dependent transcription from these DNA templates fully reconstituted in chromatin in a HeLa nuclear extract was increased 50-fold compared with naked DNA. Remarkably, under specific conditions, production of a high level of transcripts occurred at very low DNA (1 ng/microliter) and HeLa nuclear protein (1.6 micrograms/microliters) concentrations. When partially reconstituted templates were used, transcription efficiency was intermediate between that of fully reconstituted and naked DNA. These results implicate chromatin in the process of the transcriptional activation observed. Depletion from the oocyte assembly extract of an NF-I-like factor which binds in the promoter region upstream of the TATA box (-114 to -101) or deletion from the promoter of the region interacting with this factor reduced the transcriptional efficiency of the assembled templates by a factor of 5, but transcription of these templates was still 10 times higher than that of naked DNA. Together, these results indicate that the NF-I-like factor participates in the very efficient transcriptional potentiation of the vitellogenin B1 promoter which occurs during nucleosome assembly.

Transcriptional dissection of pancreatic tumors engrafted in mice.

BACKGROUND: Engraftment of primary pancreas ductal adenocarcinomas (PDAC) in mice to generate patient-derived xenograft (PDX) models is a promising platform for biological and therapeutic studies in this disease. However, these models are still incompletely characterized. Here, we measured the impact of the murine tumor environment on the gene expression of the engrafted human tumoral cells. METHODS: We have analyzed gene expression profiles from 35 new PDX models and compared them with previously published microarray data of 18 PDX models, 53 primary tumors and 41 cell lines from PDAC. The results obtained in the PDAC system were further compared with public available microarray data from 42 PDX models, 108 primary tumors and 32 cell lines from hepatocellular carcinoma (HCC). We developed a robust analysis protocol to explore the gene expression space. In addition, we completed the analysis with a functional characterization of PDX models, including if changes were caused by murine environment or by serial passing. RESULTS: Our results showed that PDX models derived from PDAC, or HCC, were clearly different to the cell lines derived from the same cancer tissues. Indeed, PDAC- and HCC-derived cell lines are indistinguishable from each other based on their gene expression profiles. In contrast, the transcriptomes of PDAC and HCC PDX models can be separated into two different groups that share some partial similarity with their corresponding original primary tumors. Our results point to the lack of human stromal involvement in PDXs as a major factor contributing to their differences from the original primary tumors. The main functional differences between pancreatic PDX models and human PDAC are the lower expression of genes involved in pathways related to extracellular matrix and hemostasis and the up- regulation of cell cycle genes. Importantly, most of these differences are detected in the first passages after the tumor engraftment. CONCLUSIONS: Our results suggest that PDX models of PDAC and HCC retain, to some extent, a gene expression memory of the original primary tumors, while this pattern is not detected in conventional cancer cell lines. Expression changes in PDXs are mainly related to pathways reflecting the lack of human infiltrating cells and the adaptation to a new environment. We also provide evidence of the stability of gene expression patterns over subsequent passages, indicating early phases of the adaptation process.

Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

Relationship between estrogen receptor alpha location and gene induction reveals the importance of downstream sites and cofactors.

BACKGROUND: To understand cancer-related modifications to transcriptional programs requires detailed knowledge about the activation of signal-transduction pathways and gene expression programs. To investigate the mechanisms of target gene regulation by human estrogen receptor alpha (hERalpha), we combine extensive location and expression datasets with genomic sequence analysis. In particular, we study the influence of patterns of DNA occupancy by hERalpha on expression phenotypes. RESULTS: We find that strong ChIP-chip sites co-localize with strong hERalpha consensus sites and detect nucleotide bias near hERalpha sites. The localization of ChIP-chip sites relative to annotated genes shows that weak sites are enriched near transcription start sites, while stronger sites show no positional bias. Assessing the relationship between binding configurations and expression phenotypes, we find binding sites downstream of the transcription start site (TSS) to be equally good or better predictors of hERalpha-mediated expression as upstream sites. The study of FOX and SP1 cofactor sites near hERalpha ChIP sites shows that induced genes frequently have FOX or SP1 sites. Finally we integrate these multiple datasets to define a high confidence set of primary hERalpha target genes. CONCLUSION: Our results support the model of long-range interactions of hERalpha with the promoter-bound cofactor SP1 residing at the promoter of hERalpha target genes. FOX motifs co-occur with hERalpha motifs along responsive genes. Importantly we show that the spatial arrangement of sites near the start sites and within the full transcript is important in determining response to estrogen signaling.

Adipose reduction in the activity of the repressor ICER is responsible for insulin resistance elicited by the transcriptional factors CREB in obesity

BACKGROUND AND AIMS: Sustained adipose activation of the transcriptional activators cAMP response binding proteins (CREB) in obesity leads to impaired expression of the glucose transporter GLUT4 and adiponectin (adipoq) in mice model of obesity. Diminution of GLUT4 and adipoq caused by CREB is indirect and relies on the increased repressive activity of the CREB target gene activating transcription factor 3 (ATF3). Specific inactivation of CREB in adipocytes decreases ATF3 production and improves whole-body insulin sensitivity of mice in the context of diet-induced obesity. Thus, elevation of CREB activity is a key mechanism responsible for adipocyte dysfunction and systemic insulin resistance. The inducible cAMP early repressor (ICER) is a negative regulator of the CREB activity. In fact, ICER antagonizes the CREB factor by competing for the regulation of similar target genes. The goal of the study was to investigate whether loss of ICER expression in adipocytes could be responsible for increased CREB activity in obesity. MATERIALS AND METHODS: Mice C57bl6 were fed with a high fat diet (HFD) for 12 weeks to increase body weight and generate insulin resistance. Biopsies of visceral adipose tissues (VAT) were prepared from human lean (BMI=24}0.5 Kg/m2) or obese subjects (BMI>35 Kg/m2). Total RNA and protein were prepared from white adipose tissues (WAT) of chow- or HFD-fed mice and VAT of lean and obese subjects. Activities of CREBs and ICER were monitored by electromobility shift assays (EMSA). The role of ICER on CREB activity was confirmed in 3T3-L1 adipocytes cells. Briefly after differentiation, the cells were electroporated with the plasmid coding for ICER cDNA. Gene expression was quantified by quantitative real-time PCR and western Blotting experiments. RESULTS: The expression of ICER is reduced in WAT of HFD-induced obese mice when compared to chow mice as measured by real-time PCR and EMSA. Similar result was found in human tissues. Reduction in ICER expression was associated with increased ATF3 expression and decreased adipoq and GLUT4 contents. Diminution in ICER levels was observed in adipocytes fraction whereas its expression was unchanged in stroma vascular fraction of WAT. Overexpression of ICER in 3T3-L1 adipocytes silenced the expression of ATF3, confirming the regulation of the factor by ICER. The expression of ICER is regulated by histone deacetylases activity (HDAC). Inhibition of HDACs in 3T3-L1 adipocytes cells using trichostatin inhibited the production of ICER. The whole activity of HDAC was reduced in WAT and VAT of obese mice and human obese subjects. CONCLUSION: Impaired adipose expression of ICER is responsible of increased CREB activity in adipocytes in obesity. This mechanism relies on reduction of the HDAC activity.

NLRC5 exclusively transactivates MHC class I and related genes through a distinctive SXY module.

MHC class II (MHCII) genes are transactivated by the NOD-like receptor (NLR) family member CIITA, which is recruited to SXY enhancers of MHCII promoters via a DNA-binding "enhanceosome" complex. NLRC5, another NLR protein, was recently found to control transcription of MHC class I (MHCI) genes. However, detailed understanding of NLRC5's target gene specificity and mechanism of action remained lacking. We performed ChIP-sequencing experiments to gain comprehensive information on NLRC5-regulated genes. In addition to classical MHCI genes, we exclusively identified novel targets encoding non-classical MHCI molecules having important functions in immunity and tolerance. ChIP-sequencing performed with Rfx5(-/-) cells, which lack the pivotal enhanceosome factor RFX5, demonstrated its strict requirement for NLRC5 recruitment. Accordingly, Rfx5-knockout mice phenocopy Nlrc5 deficiency with respect to defective MHCI expression. Analysis of B cell lines lacking RFX5, RFXAP, or RFXANK further corroborated the importance of the enhanceosome for MHCI expression. Although recruited by common DNA-binding factors, CIITA and NLRC5 exhibit non-redundant functions, shown here using double-deficient Nlrc5(-/-)CIIta(-/-) mice. These paradoxical findings were resolved by using a "de novo" motif-discovery approach showing that the SXY consensus sequence occupied by NLRC5 in vivo diverges significantly from that occupied by CIITA. These sequence differences were sufficient to determine preferential occupation and transactivation by NLRC5 or CIITA, respectively, and the S box was found to be the essential feature conferring NLRC5 specificity. These results broaden our knowledge on the transcriptional activities of NLRC5 and CIITA, revealing their dependence on shared enhanceosome factors but their recruitment to distinct enhancer motifs in vivo. Furthermore, we demonstrated selectivity of NLRC5 for genes encoding MHCI or related proteins, rendering it an attractive target for therapeutic intervention. NLRC5 and CIITA thus emerge as paradigms for a novel class of transcriptional regulators dedicated for transactivating extremely few, phylogenetically related genes.

Functional and topological characterization of transcriptional cooperativity in yeast

This work was supported by grants from Spanish Ministry of Science andInnovation (MICINN) BIO2011-22568 & BIO2008-205.

Positional bias of general and tissue-specific regulatory motifs in mouse gene promoters

Background: The arrangement of regulatory motifs in gene promoters, or promoterarchitecture, is the result of mutation and selection processes that have operated over manymillions of years. In mammals, tissue-specific transcriptional regulation is related to the presence ofspecific protein-interacting DNA motifs in gene promoters. However, little is known about therelative location and spacing of these motifs. To fill this gap, we have performed a systematic searchfor motifs that show significant bias at specific promoter locations in a large collection ofhousekeeping and tissue-specific genes.Results: We observe that promoters driving housekeeping gene expression are enriched inparticular motifs with strong positional bias, such as YY1, which are of little relevance in promotersdriving tissue-specific expression. We also identify a large number of motifs that show positionalbias in genes expressed in a highly tissue-specific manner. They include well-known tissue-specificmotifs, such as HNF1 and HNF4 motifs in liver, kidney and small intestine, or RFX motifs in testis,as well as many potentially novel regulatory motifs. Based on this analysis, we provide predictionsfor 559 tissue-specific motifs in mouse gene promoters.Conclusion: The study shows that motif positional bias is an important feature of mammalianproximal promoters and that it affects both general and tissue-specific motifs. Motif positionalconstraints define very distinct promoter architectures depending on breadth of expression andtype of tissue.

Structure-function analysis of the human TFIIB-related factor II protein reveals an essential role for the C-terminal domain in RNA polymerase III transcription.

The transcription factors TFIIB, Brf1, and Brf2 share related N-terminal zinc ribbon and core domains. TFIIB bridges RNA polymerase II (Pol II) with the promoter-bound preinitiation complex, whereas Brf1 and Brf2 are involved, as part of activities also containing TBP and Bdp1 and referred to here as Brf1-TFIIIB and Brf2-TFIIIB, in the recruitment of Pol III. Brf1-TFIIIB recruits Pol III to type 1 and 2 promoters and Brf2-TFIIIB to type 3 promoters such as the human U6 promoter. Brf1 and Brf2 both have a C-terminal extension absent in TFIIB, but their C-terminal extensions are unrelated. In yeast Brf1, the C-terminal extension interacts with the TBP/TATA box complex and contributes to the recruitment of Bdp1. Here we have tested truncated Brf2, as well as Brf2/TFIIB chimeric proteins for U6 transcription and for assembly of U6 preinitiation complexes. Our results characterize functions of various human Brf2 domains and reveal that the C-terminal domain is required for efficient association of the protein with U6 promoter-bound TBP and SNAP(c), a type 3 promoter-specific transcription factor, and for efficient recruitment of Bdp1. This in turn suggests that the C-terminal extensions in Brf1 and Brf2 are crucial to specific recruitment of Pol III over Pol II.

Integrative analysis of the Regulatory Region of the FGFR3 Oncogene

The study of transcriptional regulation often needs the integration of diverse yet independent data. In the present work, sequence conservation, predic-tion of transcription factor binding sites (TFBS) and gene expression analysis have been applied to the detection of putative transcription factor (TF) modules in the regulatory region of the FGFR3 oncogene. Several TFs with conserved binding sites in the FGFR3 regulatory region have shown high positive or negative corre-lation with FGFR3 expression both in urothelial carcinoma and in benign nevi. By means of conserved TF cluster analysis, two different TF modules have been iden-tified in the promoter and first intron of FGFR3 gene. These modules contain acti-vating AP2, E2F, E47 and SP1 binding sites plus motifs for EGR with possible repressor function.