Assessing computational tools for the discovery of transcription factor binding sites

The prediction of regulatory elements is a problem where computational methods offer great hope. Over the past few years, numerous tools have become available for this task. The purpose of the current assessment is twofold: to provide some guidance to users regarding the accuracy of currently available tools in various settings, and to provide a benchmark of data sets for assessing future tools.

Evolutionary turnover of mammalian transcription start sites

Alignments of homologous genomic sequences are widely used to identify functional genetic elements and study their evolution. Most studies tacitly equate homology of functional elements with sequence homology. This assumption is violated by the phenomenon of turnover, in which functionally equivalent elements reside at locations that are nonorthologous at the sequence level. Turnover has been demonstrated previously for transcription-factor-binding sites. Here, we show that transcription start sites of equivalent genes do not always reside at equivalent locations in the human and mouse genomes. We also identify two types of partial turnover, illustrating evolutionary pathways that could lead to complete turnover. These findings suggest that the signals encoding transcription start sites are highly flexible and evolvable, and have cautionary implications for the use of sequence-level conservation to detect gene regulatory elements.

The human sulfotransferase SULT1A1 gene is regulated in a synergistic manner by Sp1 and GA binding protein

Human sulfotransferase SULT1A1 is an important phase II xenobiotic metabolizing enzyme that is highly expressed in the liver and mediates the sulfonation of drugs, carcinogens, and steroids. Until this study, the transcriptional regulation of the SULT1A subfamily had been largely unexplored. Preliminary experiments in primary human hepatocytes showed that SULT1A mRNA levels were not changed in response to nuclear receptor activators, such as dexamethasone and 3-methylcolanthrene, unlike other metabolizing enzymes. Using HepG2 cells, the high activity of the TATA-less SULT1A1 promoter was shown to be dependent on the presence of Sp1 and Ets transcription factor binding sites (EBS), located within - 112 nucleotides from the transcriptional start site. The homologous promoter of the closely related SULT1A3 catecholamine sulfotransferase, which is expressed at negligible levels in the adult liver, displayed 70% less activity than SULT1A1. This was shown to be caused by a two-base pair difference in the EBS. The Ets transcription factor GA binding protein (GABP) was shown to bind the SULT1A1 EBS and could transactivate the SULT1A1 promoter in Drosophila melanogaster S2 cells. Cotransfection of Sp1 could synergistically enhance GABP-mediated activation by 10-fold. Although Sp1 and GABP alone could induce SULT1A3 promoter activity, the lack of the EBS on this promoter prevented a synergistic interaction between the two factors. This study reports the first insight into the transcriptional regulation of the SULT1A1 gene and identifies a crucial difference in regulation of the closely related SULT1A3 gene, which accounts for the two enzymes' differential expression patterns observed in the adult liver.

Identification of functional sequences in the pregenomic RNA promoter of the Banana streak virus Cavendish strain (BSV-Cav)

The promoter regions of plant pararetroviruses direct transcription of the full-length viral genome into a pregenomic RNA that is an intermediate in the replication of the virus. It serves as template for reverse transcription and as polycistronic mRNA for translation to viral proteins. We have identified functional promoter elements in the intergenic region of the Cavendish isolate of Banana streak virus (BSV-Cav), a member of the genus Badnavirus. Potential binding sites for plant transcription factors were found both upstream and downstream of the transcription start site by homology search in the PLACE database of plant cis-acting elements. The functionality of these putative cis-acting elements was tested by constructing loss-of-function and regain-of-function mutant promoters whose activity was quantified in embryogenic sugarcane suspension cells. Four regions that are important for activity of the BSV-Cav promoter were identified: the region containing an as-l-like element, the region around-141 and down to -77, containing several putative transcription factor binding sites, the region including the CAAT-box, and the leader region. The results could help explain the high BSV-Cav promoter activity that was observed previously in transgenic sugarcane plants and give more insight into the plant cell-mediated replication of the viral genome in banana streak disease. (C) 2004 Elsevier B.V. All rights reserved.

The rat Na+-sulfate cotransporter rNaS2: functional characterization, tissue distribution, and gene (slc13a4) structure

Inorganic sulfate is essential for numerous functions in mammalian physiology. In the present study, we characterized the functional properties of the rat Na+-sulfate cotransporter NaS2 (rNaS2), determined its tissue distribution, and identified its gene (slc13a4) structure. Expression of rNaS2 protein in Xenopus oocytes led to a Na+-dependent transport of sulfate that was inhibited by phosphate, thiosulfate, tungstate, selenate, oxalate, and molybdate, but not by citrate, succinate, or DIDS. Transport kinetics of rNaS2 determined a K-M for sulfate of 1.26 mM. Na+ kinetics determined a Hill coefficient of n=3.0 +/- 0.7, suggesting a Na+:SO42- stoichiometry of 3:1. rNaS2 mRNA was highly expressed in placenta, with lower levels found in the brain and liver. slc13a4 maps to rat chromosome 4 and contains 17 exons, spanning over 46 kb in length. This gene produces two alternatively spliced transcripts, of which the transcript lacking exon 2 is the most abundant form. Its 5' flanking region contains CAAT- and GC-box motifs and a number of putative transcription factor binding sites, including GATA-1, SP1, and AP-2 consensus sequences. This is the first study to characterize rNaS2 transport kinetics, define its tissue distribution, and resolve its gene (slc13a4) structure and 5' flanking region.

Genomic organization of human arylamine N-acetyltransferase Type I reveals alternative promoters that generate different 5'-UTR splice variants with altered translational activities

In humans, a polymorphic gene encodes the drug-metabolizing enzyme NATI (arylamine N-acetyltransferase Type 1), which is widely expressed throughout the body. While the protein-coding region of NATI is contained within a single exon, examination of the human EST (expressed sequence tag) database at the NCBI revealed the presence of nine separate exons, eight of which were located in the 5'non-coding region of NATI. Differential splicing produced at least eight unique mRNA isoforms that could be grouped according to the location of the first exon, which suggested that NATI expression occurs from three alternative promoters. Using RT (reverse transcriptase)-PCR, we identified one major transcript in various epithelial cells derived from different tissues. In contrast, multiple transcripts were observed in blood-derived cell lines (CEM, THP-1 and Jurkat), with a novel variant, not identified in the EST database, found in CEM cells only. The major splice variant increased gene expression 9-11-fold in a luciferase reporter assay, while the other isoforrns were similar or slightly greater than the control. We examined the upstream region of the most active splice variant in a promoter-reporter assay, and isolated a 257 bp sequence that produced maximal promoter activity. This sequence lacked a TATA box, but contained a consensus Sp1 site and a CAAT box, as well as several other putative transcription-factor-binding sites. Cell-specific expression of the different NATI transcripts may contribute to the variation in NATI activity in vivo.

Functional characterization and genomic organization of the human Na+-sulfate cotransporter hNaS2 gene (SLC13A4)

Sulfate plays an essential role in human growth and development. Here, we characterized the functional properties of the human Na+-sulfate cotransporter (hNaS2), determined its tissue distribution, and identified its gene (SLC13A4) structure. Expression of hNaS2 protein in Xenopus oocytes led to a Na+-dependent transport of sulfate that was inhibited by thiosulfate, phosphate, molybdate. selenate and tungstate, but not by oxalate, citrate, succinate, phenol red or DIDS. Transport kinetics of hNaS2 determined a K, for sulfate of 0.38 mM, suggestive of a high affinity sulfate transporter. Na+ kinetics determined a Hill coefficient of 1.6 +/- 0.6, suggesting a Na: SO42- stoichiometry of 2:1. hNaS2 mRNA was highly expressed in placenta and testis, with intermediate levels in brain and lower levels found in the heart, thymus, and liver. The SLC13A4 gene contains 16 exons, spanning over 47 kb in length. Its 5'-flanking region contains CAAT- and GC-box motifs, and a number of putative transcription factor binding sites, including GATA-1, AP-1, and AP-2 consensus sequences. This is the first study to characterize hNaS2 transport kinetics, define its tissue distribution, and resolve its gene (SLC13A4) structure and 5' flanking region. (C) 2004 Elsevier Inc. All rights reserved.

Transcriptional network dynamics in macrophage activation

Transcriptional regulatory networks govern cell differentiation and the cellular response to external stimuli. However, mammalian model systems have not yet been accessible for network analysis. Here, we present a genome-wide network analysis of the transcriptional regulation underlying the mouse macrophage response to bacterial lipopolysaccharide (LPS). Key to uncovering the network structure is our combination of time-series cap analysis of gene expression with in silico prediction of transcription factor binding sites. By integrating microarray and qPCR time-series expression data with a promoter analysis, we find dynamic subnetworks that describe how signaling pathways change dynamically during the progress of the macrophage LPS response, thus defining regulatory modules characteristic of the inflammatory response. In particular, our integrative analysis enabled us to suggest novel roles for the transcription factors ATF-3 and NRF-2 during the inflammatory response. We believe that our system approach presented here is applicable to understanding cellular differentiation in higher eukaryotes. (c) 2006 Elsevier Inc. All rights reserved.

Molecular cloning and characterization of the mouse Na+ sulfate cotransporter gene (Slc13a4): structure and expression

Sulfate is an essential ion required for numerous functions in mammalian physiology. Due to its hydrophilic nature, cells require sulfate transporters on their plasma membranes to allow entry of sulfate into cells. In this study, we identified a new mouse Na+-sulfate cotransporter (mNaS2), characterized its tissue distribution and determined its cDNA and gene (Slc13a4) structures. mNaS2 mRNA was expressed in placenta, brain, lung, eye, heart, testis, thymus and liver. The mouse NaS2 cDNA spans 3384 nucleotides and its open frame encodes a protein of 624 amino acids. Slc13a4 maps to mouse chromosome 6131 and contains 16 exons, spanning over 40 kb in length. Its 5'-flanking region contains CART- and GC-box motifs and a number of putative transcription factor binding sites, including GATA-1, MTF-1, STAT6 and HNF4 consensus sequences. This is the first study to define the tissue distribution of mNaS2 and resolve its cDNA and gene structures, which will allow us to investigate mNaS2 gene expression in vivo and determine its role in mammalian physiology.

neptune, a Kruppel-like transcription factor that participates in primitive erythropoiesis in Xenopus

The specification of the erythroid lineage from hematopoietic stem cells requires the expression and activity of lineage-specific transcription factors. One transcription factor family that has several members involved in hematopoiesis is the Kruppel-like factor (KLF) family [1]. For example, erythroid KLF (EKLF) regulates beta -globin expression during erythroid differentiation [2-6]. KLFs share a highly conserved zinc finger-based DNA binding domain (DBD) that mediates binding to CACCC-box and GC-rich sites, both of which are frequently found in the promoters of hematopoietic genes. Here, we identified a novel Xenopus KLF gene, neptune, which is highly expressed in the ventral blood island (VBI), cranial ganglia, and hatching and cement glands. neptune expression is induced in response to components of the BMP-4 signaling pathway in injected animal cap explants. Similar to its family member, EKLF, Neptune can bind CACCC-box and GC-rich DNA elements. We show that Neptune cooperates with the hematopoietic transcription factor XGATA-1 to enhance globin induction in animal cap explants. A fusion protein comprised of Neptune's DBD and the Drosophila engrailed repressor domain suppresses the induction of globin in ventral marginal zones and in animal caps. These studies demonstrate that Neptune is a positive regulator of primitive erythropoiesis in Xenopus.

Genetic and physical interactions between Microphthalmia transcription factor and PU.1 are necessary for osteoclast gene expression and differentiation

The microphthalmia transcription factor (MITF), a basic-helix-loop-helix zipper factor, regulates distinct target genes in several cell types. We hypothesized that interaction with the Ets family factor PU.1, whose expression is limited to hematopoietic cells, might be necessary for activation of target genes like tartrate-resistant acid phosphatase (TRAP) in osteoclasts. Several lines of evidence were consistent with this model. The combination of MITF and PU.1 synergistically activated the TRAP promoter in transient assays. This activation was dependent on intact binding sites for both factors in the TRAP promoter. MITF and PU.1 physically interacted when coexpressed in COS cells or in vitro when purified recombinant proteins were studied. The minimal regions of MITF and PU.1 required for the interaction were the basic-helix-loop-helix zipper domain and the Ets DNA binding domain, respectively. Significantly, mice heterozygous for both the mutant mi allele and a PU.1 null allele developed osteopetrosis early in life which resolved with age. The size and number of osteoclasts were not altered in the double heterozygous mutant mice, indicating that the defect lies in mature osteoclast function. Taken in total, the results afford an example of how lineage-specific gene regulation can be achieved by the combinatorial action of two broadly expressed transcription factors.

Transcription factor Tfec contributes to the IL-4-inducible expression of a small group of genes in mouse macrophages including the granulocyte colony-stimulating factor receptor

Expression of the mouse transcription factor EC (Tfec) is restricted to the myeloid compartment, suggesting a function for Tfec in the development or function of these cells. However, mice lacking Tfec develop normally, indicating a redundant role for Tfec in myeloid cell development. We now report that Tfec is specifically induced in bone marrow-derived macrophages upon stimulation with the Th2 cytokines, IL-4 and IL-13, or LPS. LPS induced a rapid and transient up-regulation of Tfec mRNA expression and promoter activity, which was dependent on a functional NF-kappa B site. IL-4, however, induced a rapid, but long-lasting, increase in Tfec mRNA, which, in contrast to LPS stimulation, also resulted in detectable levels of Tfec protein. IL-4-induced transcription of Tfec was absent in macrophages lacking Stat6, and its promoter depended on two functional Stat6-binding sites. A global comparison of IL-4-induced genes in both wild-type and Tfec mutant macrophages revealed a surprisingly mild phenotype with only a few genes affected by Tfec deficiency. These included the G-CSFR (Csf3r) gene that was strongly up-regulated by IL-4 in wild-type macrophages and, to a lesser extent, in Tfec mutant macrophages. Our study also provides a general definition of the transcriptome in alternatively activated mouse macrophages and identifies a large number of novel genes characterizing this cell type.

The Runx1 transcription factor controls CSF-1-dependent and -independent growth and survival of macrophages

Gene translocations that repress the function of the Runx1 transcription factor play a critical role in the development of myeloid leukemia. In this report, we demonstrate that Runx1 precisely regulates c-fms (CSF-1 receptor) gene expression. Runx1 controlled expression by binding to multiple sites within the mouse c-fms gene, allowing interaction between promoter and downstream enhancer elements. The runx1 and c-fms genes showed an identical pattern of expression in mature macrophages. Runx1 expression was repressed in CSF-1 stimulated, proliferating bone marrow-derived macrophages (BMM) and significantly increased in quiescent, CSF-1 starved cells. The RAW264.7 and Mono-Mac-6, macrophage-like cell lines expressed low levels of Runx1 and both showed growth arrest and cell death with ectopic expression of Runx1. The EM-3 cell line, which represents an early myeloid progenitor cell line, showed growth arrest with Runx1 expression in the absence of any detectable changes in cell differentiation. These findings suggest that Runx1 regulates growth and survival of myeloid cells and provide a novel insight into the role of Runx family gene translocations in leukemogenesis.

Differentiation of the mononuclear phagocyte system during mouse embryogenesis: The role of transcription factor PU.1

During mouse embryogenesis, macrophage-like cells arise first in the yolk sac and are produced subsequently in the liver. The onset of liver hematopoiesis is associated with the transition from primitive to definitive erythrocyte production. This report addresses the hypothesis that a similar transition in phenotype occurs in myelopoiesis. We have used whole mount in situ hybridization to detect macrophage-specific genes expressed during mouse development. The mouse c-fms mRNA, encoding the receptor for macrophage colony-stimulating factor (CSF-1), was expressed on phagocytic cells in the yolk sac and throughout the embryo before the onset of liver hematopoiesis, Similar cells were detected using the mannose receptor, the complement receptor (CR3), or the Microphthalmia transcription factor (MITF) as mRNA markers. By contrast, other markers including the F4/80 antigen, the macrophage scavenger receptor, the S-100 proteins, S100A8 and S100A9, and the secretory product lysozyme appeared later in development and appeared restricted to only a subset of c-fms-positive cells. Two-color immunolabeling on disaggregated cells confirmed that CR3 and c-fms proteins are expressed on the same cells. Among the genes appearing later in development was the macrophage-restricted transcription factor, PU.1, which has been shown to be required for normal adult myelopoiesis. Mice with null mutations in PU.1 had normal numbers of c-fms-positive phagocytes at 11.5dpc. PU.1(-/-) embryonic stem cells were able to give rise to macrophagelike cells after cultivation in vitro. The results support previous evidence that yolk sac-derived fetal phagocytes are functionally distinct from those arising in the liver and develop via a different pathway. (C) 1999 by The American Society of Hematology.