One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Insulin and insulin-like growth factor I receptors utilize different G protein signaling components.

We examined the role of heterotrimeric G protein signaling components in insulin and insulin-like growth factor I (IGF-I) action. In HIRcB cells and in 3T3L1 adipocytes, treatment with the Galpha(i) inhibitor (pertussis toxin) or microinjection of the Gbetagamma inhibitor (glutathione S-transferase-betaARK) inhibited IGF-I and lysophosphatidic acid-stimulated mitogenesis but had no effect on epidermal growth factor (EGF) or insulin action. In basal state, Galpha(i) and Gbeta were associated with the IGF-I receptor (IGF-IR), and after ligand stimulation the association of IGF-IR with Galpha(i) increased concomitantly with a decrease in Gbeta association. No association of Galpha(i) was found with either the insulin or EGF receptor. Microinjection of anti-beta-arrestin-1 antibody specifically inhibited IGF-I mitogenic action but had no effect on EGF or insulin action. beta-Arrestin-1 was associated with the receptors for IGF-I, insulin, and EGF in a ligand-dependent manner. We demonstrated that Galpha(i), betagamma subunits, and beta-arrestin-1 all play a critical role in IGF-I mitogenic signaling. In contrast, neither metabolic, such as GLUT4 translocation, nor mitogenic signaling by insulin is dependent on these protein components. These results suggest that insulin receptors and IGF-IRs can function as G protein-coupled receptors and engage different G protein partners for downstream signaling.

The transcription factor c-Maf controls touch receptor development and function.

The sense of touch relies on detection of mechanical stimuli by specialized mechanosensory neurons. The scarcity of molecular data has made it difficult to analyze development of mechanoreceptors and to define the basis of their diversity and function. We show that the transcription factor c-Maf/c-MAF is crucial for mechanosensory function in mice and humans. The development and function of several rapidly adapting mechanoreceptor types are disrupted in c-Maf mutant mice. In particular, Pacinian corpuscles, a type of mechanoreceptor specialized to detect high-frequency vibrations, are severely atrophied. In line with this, sensitivity to high-frequency vibration is reduced in humans carrying a dominant mutation in the c-MAF gene. Thus, our work identifies a key transcription factor specifying development and function of mechanoreceptors and their end organs.

The transcription factor NFATc2 controls IL-6-dependent T cell activation in experimental colitis.

The nuclear factor of activated T cells (NFAT) family of transcription factors controls calcium signaling in T lymphocytes. In this study, we have identified a crucial regulatory role of the transcription factor NFATc2 in T cell-dependent experimental colitis. Similar to ulcerative colitis in humans, the expression of NFATc2 was up-regulated in oxazolone-induced chronic intestinal inflammation. Furthermore, NFATc2 deficiency suppressed colitis induced by oxazolone administration. This finding was associated with enhanced T cell apoptosis in the lamina propria and strikingly reduced production of IL-6, -13, and -17 by mucosal T lymphocytes. Further studies using knockout mice showed that IL-6, rather than IL-23 and -17, are essential for oxazolone colitis induction. Administration of hyper-IL-6 blocked the protective effects of NFATc2 deficiency in experimental colitis, suggesting that IL-6 signal transduction plays a major pathogenic role in vivo. Finally, adoptive transfer of IL-6 and wild-type T cells demonstrated that oxazolone colitis is critically dependent on IL-6 production by T cells. Collectively, these results define a unique regulatory role for NFATc2 in colitis by controlling mucosal T cell activation in an IL-6-dependent manner. NFATc2 in T cells thus emerges as a potentially new therapeutic target for inflammatory bowel diseases.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

The δ-Opioid Receptor Affects Epidermal Homeostasis via ERK-Dependent Inhibition of Transcription Factor POU2F3.

Neuropeptides and their receptors are present in human skin, and their importance for cutaneous homeostasis and during wound healing is increasingly appreciated. However, there is currently a lack of understanding of the molecular mechanisms by which their signaling modulates keratinocyte function. Here, we show that δ-opioid receptor (DOPr) activation inhibits proliferation of human keratinocytes, resulting in decreased epidermal thickness in an organotypic skin model. DOPr signaling markedly delayed induction of keratin intermediate filament (KRT10) during in vitro differentiation and abolished its induction in the organotypic skin model. This was accompanied by deregulation of involucrin (IVL), loricrin, and filaggrin. Analysis of the transcription factor POU2F3, which is involved in regulation of KRT10, IVL, and profilaggrin expression, revealed a DOPr-mediated extracellular signal-regulated kinase (ERK)-dependent downregulation of this factor. We propose that DOPr signaling specifically activates the ERK 1/2 mitogen-activated protein kinase pathway to regulate keratinocyte functions. Complementing our earlier studies in DOPr-deficient mice, these data suggest that DOPr activation in human keratinocytes profoundly influences epidermal morphogenesis and homeostasis.

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.

The Ciliogenic Transcription Factor RFX3 Regulates Early Midline Distribution of Guidepost Neurons Required for Corpus Callosum Development.

The corpus callosum (CC) is the major commissure that bridges the cerebral hemispheres. Agenesis of the CC is associated with human ciliopathies, but the origin of this default is unclear. Regulatory Factor X3 (RFX3) is a transcription factor involved in the control of ciliogenesis, and Rfx3-deficient mice show several hallmarks of ciliopathies including left-right asymmetry defects and hydrocephalus. Here we show that Rfx3-deficient mice suffer from CC agenesis associated with a marked disorganisation of guidepost neurons required for axon pathfinding across the midline. Using transplantation assays, we demonstrate that abnormalities of the mutant midline region are primarily responsible for the CC malformation. Conditional genetic inactivation shows that RFX3 is not required in guidepost cells for proper CC formation, but is required before E12.5 for proper patterning of the cortical septal boundary and hence accurate distribution of guidepost neurons at later stages. We observe focused but consistent ectopic expression of Fibroblast growth factor 8 (Fgf8) at the rostro commissural plate associated with a reduced ratio of GLIoma-associated oncogene family zinc finger 3 (GLI3) repressor to activator forms. We demonstrate on brain explant cultures that ectopic FGF8 reproduces the guidepost neuronal defects observed in Rfx3 mutants. This study unravels a crucial role of RFX3 during early brain development by indirectly regulating GLI3 activity, which leads to FGF8 upregulation and ultimately to disturbed distribution of guidepost neurons required for CC morphogenesis. Hence, the RFX3 mutant mouse model brings novel understandings of the mechanisms that underlie CC agenesis in ciliopathies.

Coordination between zinc and phosphate homeostasis involves the transcription factor PHR1, the phosphate exporter PHO1, and its homologue PHO1;H3 in Arabidopsis.

Interactions between zinc (Zn) and phosphate (Pi) nutrition in plants have long been recognized, but little information is available on their molecular bases and biological significance. This work aimed at examining the effects of Zn deficiency on Pi accumulation in Arabidopsis thaliana and uncovering genes involved in the Zn-Pi synergy. Wild-type plants as well as mutants affected in Pi signalling and transport genes, namely the transcription factor PHR1, the E2-conjugase PHO2, and the Pi exporter PHO1, were examined. Zn deficiency caused an increase in shoot Pi content in the wild type as well as in the pho2 mutant, but not in the phr1 or pho1 mutants. This indicated that PHR1 and PHO1 participate in the coregulation of Zn and Pi homeostasis. Zn deprivation had a very limited effect on transcript levels of Pi-starvation-responsive genes such as AT4, IPS1, and microRNA399, or on of members of the high-affinity Pi transporter family PHT1. Interestingly, one of the PHO1 homologues, PHO1;H3, was upregulated in response to Zn deficiency. The expression pattern of PHO1 and PHO1;H3 were similar, both being expressed in cells of the root vascular cylinder and both localized to the Golgi when expressed transiently in tobacco cells. When grown in Zn-free medium, pho1;h3 mutant plants displayed higher Pi contents in the shoots than wild-type plants. This was, however, not observed in a pho1 pho1;h3 double mutant, suggesting that PHO1;H3 restricts root-to-shoot Pi transfer requiring PHO1 function for Pi homeostasis in response to Zn deficiency.

Identification of the minimal promoter region of the mouse NKX5-3, a transcription factor implicated in eye development.

Early ocular development is controlled by a complex network of transcription factors, cell cycle regulators, and diffusible signalling molecules. Together, these molecules regulate cell proliferation and apoptosis, and specify retinal fate. NKX5-3 is a homeobox transcription factor implicated in eye development. The analysis of the 5'-flanking region of the mouse Nkx5-3 gene revealed a predicted TATA-less promoter sequence between -416 and -166 of the translation start site. To functionally characterise Nkx5-3 promoter activity, serial deletions of the promoter sequence were introduced in pGL-3 basic vector and promoter activity of these 5'- and 3'-deleted constructions was tested in HeLa and CHO cells. Transactivation assays identified a region between -350 and -296 exhibiting promoter-like activity. Combined analysis by deletions and point mutations showed that this sequence, containing multiple Sp1 binding sites was necessary to promote transcriptional activity. Binding of Sp1 to this region was confirmed by electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation, using an antibody specific for Sp1. Altogether, these results demonstrated that the immediate upstream region of Nkx5-3 gene possessed a strong intrinsic promoter activity in vitro, suggesting a potential role in Nkx5-3 transcription in vivo.

The complementary strand of the human T-cell leukemia virus type 1 RNA genome encodes a bZIP transcription factor that down-regulates viral transcription.

The RNA genome of the human T-cell leukemia virus type 1 (HTLV-1) codes for proteins involved in infectivity, replication, and transformation. We report in this study the characterization of a novel viral protein encoded by the complementary strand of the HTLV-1 RNA genome. This protein, designated HBZ (for HTLV-1 bZIP factor), contains a N-terminal transcriptional activation domain and a leucine zipper motif in its C terminus. We show here that HBZ is able to interact with the bZIP transcription factor CREB-2 (also called ATF-4), known to activate the HTLV-1 transcription by recruiting the viral trans-activator Tax on the Tax-responsive elements (TxREs). However, we demonstrate that the HBZ/CREB-2 heterodimers are no more able to bind to the TxRE and cyclic AMP response element sites. Taking these findings together, the functional inactivation of CREB-2 by HBZ is suggested to contribute to regulation of the HTLV-1 transcription. Moreover, the characterization of a minus-strand gene protein encoded by HTLV-1 has never been reported until now.

Role of the transcription factor sox4 in schwann cell development

Previous studies demonstrated that both Schwann cell differentiation and de-differentiation (in the situation of a nerve injury or demyelinating disease) are regulated by cell-intrinsic regulators including several transcription factors. In particular, the de-differentiation of mature Schwann cells is driven by the activation of multiple negative regulators of myelination including c-Jun, Notch, Sox-2 and Pax-3, all usually expressed in the immature Schwann cells and suppressed at the onset of myelination. In order to identify new negative regulators of myelination involved in the development of the peripheral nervous system (PNS) we analyzed the data from a previously performed transcriptional analysis of myelinating Schwann cells. Based on its transcriptional expression profile during myelination, Sox4, a member of the Sox gene family, was identified as a potential candidate. Previous studies demonstrated that prolonged Sox4 expression in oligodendrocytes maintains these cells in a premyelinating state, further suggesting its role as a negative regulator of myelination. Concomitantly, we observed upregulation of Sox4 mRNA and protein expression levels in the PNS of three different models of demyelinating neuropathies (Pmp22, Lpin1, and Scap KOs). To better characterize the molecular function of Sox4, we used a viral vector allowing Sox4 overexpression in cultured Schwann cells and in neuron-Schwann cell co-cultures. In parallel, we generated two transgenic lines of mice in which the overexpression of Sox4 is driven specifically in Schwann cells by the Myelin Protein Zero gene promoter. The preliminary data from these in vitro and in vivo experiments show that overexpression of Sox4 in PNS causes a delay in progression of myelination thus indicating that Sox4 acts as a negative regulator of Schwann cell myelination.

Prospective study of insulin-like growth factor-I, insulin-like growth factor-binding protein 3, genetic variants in the IGF1 and IGFBP3 genes and risk of coronary artery disease.

Although experimental studies have suggested that insulin-like growth factor I (IGF-I) and its binding protein IGFBP-3 might have a role in the aetiology of coronary artery disease (CAD), the relevance of circulating IGFs and their binding proteins in the development of CAD in human populations is unclear. We conducted a nested case-control study, with a mean follow-up of six years, within the EPIC-Norfolk cohort to assess the association between circulating levels of IGF-I and IGFBP-3 and risk of CAD in up to 1,013 cases and 2,055 controls matched for age, sex and study enrolment date. After adjustment for cardiovascular risk factors, we found no association between circulating levels of IGF-I or IGFBP-3 and risk of CAD (odds ratio: 0.98 (95% Cl 0.90-1.06) per 1 SD increase in circulating IGF-I; odds ratio: 1.02 (95% Cl 0.94-1.12) for IGFBP-3). We examined associations between tagging single nucleotide polymorphisms (tSNPs) at the IGF1 and IGFBP3 loci and circulating IGF-I and IGFBP-3 levels in up to 1,133 cases and 2,223 controls and identified three tSNPs (rs1520220, rs3730204, rs2132571) that showed independent association with either circulating IGF-I or IGFBP-3 levels. In an assessment of 31 SNPs spanning the IGF1 or IGFBP3 loci, none were associated with risk of CAD in a meta-analysis that included EPIC-Norfolk and eight additional studies comprising up to 9,319 cases and 19,964 controls. Our results indicate that IGF-I and IGFBP-3 are unlikely to be importantly involved in the aetiology of CAD in human populations.

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

Role of forkhead transcription factor FOXC2 in lymphatic vascular developement

Le système vasculaire lymphatique est le second réseau de vaisseaux du corps humain. Sa fonction principale est de retourner le fluide interstitiel excédentaire au système cardiovasculaire. Il est également impliqué dans la défense immunitaire de l'organisme, ainsi que dans le transport initial des graisses alimentaires. De multiples pathologies sont associées au dysfonctionnement du développement vasculaire lymphatique, dont les lymphoedèmes. Un des gènes clés dans le contrôle de l'étape de maturation du système lymphatique est le facteur de transcription FOXC2. De précédentes études utilisant des modèles génétiques mutins déficients en Foxc2 ont montré son rôle dans la régulation du processus de spécification des vaisseaux lymphatiques en capillaires versus vaisseaux collecteurs, ainsi que dans la formation des valves lymphatiques. Chez l'homme, les mutations dans le gène FOXC2 causent le syndrome lymphoedème- distichiasis. Dans ce travail, nous avons étudié les mécanismes moléculaires qui régulent l'expression et l'activité de FOXC2 dans les vaisseaux lymphatiques. Nous avons découvert que la fonction de FOXC2 est régulée par phosphorylation de la protéine, qui détermine son activité transcriptionnelle au niveau génomique, jouant ainsi un rôle important dans le développement vasculaire in vivo. Les vaisseaux lymphatiques sont soumis à des forces de stress générées par le flux de la lymphe (FSS). Nous avons donc testé l'hypothèse que ces forces contribuent à la morphogenèse et à l'organisation des vaisseaux lymphatiques. In vitro, les cellules endothéliales lymphatiques répondent aux forces mécaniques, qui induisent l'expression de FOXC2, activent la voie de signalisation Ca2+/calcineurin/NFATcl et régulent l'expression de la protéine de jonction gap connexin37. Nous avons également montré que le stress de flux mécanique, FOXC2, calcineurin/NFATcl et connexin37 coopèrent dans le contrôle de la maturation des vaisseaux lymphatiques in vivo. En dernier lieu, nous avons cherché à identifier les récepteurs de surface cellulaires permettant le transfert du signal de stress mécanique qui induit l'expression de FOXC2. Nous présentons ici des données préliminaires, qui suggèrent le rôle de la voie de signalisation TGFß ainsi que l'implication des jonctions adhérentes dans ce processus. En conclusion, la présente étude met en lumière les mécanismes de l'activité de FOXC2 dans les cellules endothéliales lymphatiques et l'importance du rôle des forces mécaniques de flux dans le contrôle de son l'expression, ainsi que dans le développement et la fonction du système vasculaire lymphatique. - The lymphatic vascular system is a second vascular system of human body. Its main fonction is to transfer excess interstitial fluid back to cardiovascular system. In addition, it is involved in immune defense and responsible for the uptake of dietary fat. A number of pathologies called lymphedemas are associated with lymphatic vascular system dysfunction. Hereditary lymphedemas are caused by mutations in genes controlling lymphatic vascular development. One of the key genes responsible for lymphatic vascular maturation is forkhead transcription factor FOXC2. Previous studies of Foxc2 knockout mice showed that Foxc2 controls the process of lymphatic capillary versus collecting vessel fate specification and formation of lymphatic valves. Importantly, mutations in FOXC2 cause human lymphedema-distichiasis syndrome. In this work we investigated the molecular mechanisms regulating the expression and activity of FOXC2 in lymphatic vasculature. We discovered that FOXC2 function is regulated by phosphorylation. We describe how phosphorylation controls FOXC2 transcriptional activity on a genome-wide level and show that FOXC2 phosphorylation plays an important role in vascular development in vivo. Lymphatic vessels are subjected to fluid shear stress (FSS). Therefore we investigated whether mechanical forces contribute to lymphatic vascular patterning and morphogenesis. We found that FSS induces the expression of FOXC2, activates Ca2+/calcineurin/NFATcl signaling and induces the expression of gap junction protein connexin37 in lymphatic endothelial cells in vitro. Importantly, we were able to show that shear stress, FOXC2, calcineurin/NFATcl and connexin37, control maturation of lymphatic vessels in vivo. Finally, we searched for cell surface receptors that mediate the induction of FOXC2 by shear stress, and we present some preliminary data, suggesting the role of TGF-beta signaling and adherens junctions in this process. In conclusion, the present study sheds light on the mechanisms of FOXC2 activity and suggests an important role of mechanical forces in controlling FOXC2 expression as well as lymphatic system development and function.