Transcriptional regulation of the stem cell leukemia gene (SCL)--comparative analysis of five vertebrate SCL loci

Alteration of CTCF-associated chromatin neighborhood inhibits TAL1-driven oncogenic transcription program and leukemogenesis

Aberrant activation of the TAL1 is associated with up to 60% of T-ALL cases and is involved in CTCF-mediated genome organization within the TAL1 locus, suggesting that CTCF boundary plays a pathogenic role in T-ALL. Here, we show that -31-Kb CTCF binding site (-31CBS) serves as chromatin boundary that defines topologically associating domain (TAD) and enhancer/promoter interaction required for TAL1 activation. Deleted or inverted -31CBS impairs TAL1 expression in a context-dependent manner. Deletion of -31CBS reduces chromatin accessibility and blocks long-range interaction between the +51 erythroid enhancer and TAL1 promoter-1 leading to inhibition of TAL1 expression in erythroid cells, but not T-ALL cells. However, in TAL1-expressing T-ALL cells, the leukemia-prone TAL1 promoter-IV specifically interacts with the +19 stem cell enhancer located 19 Kb downstream of TAL1 and this interaction is disrupted by the -31CBS inversion in T-ALL cells. Inversion of -31CBS in Jurkat cells alters chromatin accessibility, histone modifications and CTCF-mediated TAD leading to inhibition of TAL1 expression and TAL1-driven leukemogenesis. Thus, our data reveal that -31CBS acts as critical regulator to define +19-enhancer and the leukemic prone promoter IV interaction for TAL1 activation in T-ALL. Manipulation of CTCF boundary can alter TAL1 TAD and oncogenic transcription networks in leukemogenesis.

Overexpression of the immediate early response 5 gene increases the radiosensitivity of HeLa cells

The effects of the immediate early response 5 (IER5) gene on the sensitivity of HeLa cells to radiation remain unclear. In the present study, stably transfected HeLa cells resulting in the knockdown or overexpression of IER5 were investigated. In addition, xenografts of normal, IER5-silenced and -overexpressed HeLa cells were injected into nude mice and examined. The results demonstrated that the radiosensitivity of the IER5-overexpressed HeLa cells was significantly increased compared with that of the normal and IER5-silenced cells. The upregulation of IER5 effectively decreased cell proliferation and IER5 silencing promoted cell proliferation compared with that in the normal HeLa cells. Following irradiation of the cells with IER5 knockdown, cell cycle was arrested at the G₂/M phase and an increase in the proportion of S phase cells was observed. By contrast, the overexpression of IER5 led to an increase in the proportion of G₁ phase cells. Furthermore, the upregulation of IER5 inhibited tumor growth in vivo. The present findings demonstrate that the IER5 gene affects the radiosensitivity of HeLa cells and serves an important role in cell proliferation, suggesting that this gene may be a potential radiotherapeutic target in cervical cancer.

Robust Identification of Developmentally Active Endothelial Enhancers in Zebrafish Using FANS-Assisted ATAC-Seq

Identification of tissue-specific and developmentally active enhancers provides insights into mechanisms that control gene expression during embryogenesis. However, robust detection of these regulatory elements remains challenging, especially in vertebrate genomes. Here, we apply fluorescent-activated nuclei sorting (FANS) followed by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) to identify developmentally active endothelial enhancers in the zebrafish genome. ATAC-seq of nuclei from Tg(fli1a:egfp)^y1 transgenic embryos revealed expected patterns of nucleosomal positioning at transcriptional start sites throughout the genome and association with active histone modifications. Comparison of ATAC-seq from GFP-positive and -negative nuclei identified more than 5,000 open elements specific to endothelial cells. These elements flanked genes functionally important for vascular development and that displayed endothelial-specific gene expression. Importantly, a majority of tested elements drove endothelial gene expression in zebrafish embryos. Thus, FANS-assisted ATAC-seq using transgenic zebrafish embryos provides a robust approach for genome-wide identification of active tissue-specific enhancer elements.

Functional 5' UTR motif discovery with LESMoN: Local Enrichment of Sequence Motifs in biological Networks

Biological networks are rich representations of the relationships between entities such as genes or proteins and have become increasingly complete thanks to various high-throughput network mapping experimental approaches. Here, we propose a method to use such networks to guide the search for functional sequence motifs. Specifically, we introduce Local Enrichment of Sequence Motifs in biological Networks (LESMoN), an enumerative motif discovery algorithm that identifies 5' untranslated region (UTR) sequence motifs whose associated proteins form unexpectedly dense clusters in a given biological network. When applied to the human protein-protein interaction network from BioGRID, LESMoN identifies several highly significant 5' UTR sequence motifs, including both previously known motifs and uncharacterized ones. The vast majority of these motifs are evolutionary conserved and the genes containing them are significantly enriched for various gene ontology terms suggesting new associations between 5' UTR motifs and a number of biological processes. We validate in vivo the role in protein expression regulation of three motifs identified by LESMoN.

The hierarchical organization of natural protein interaction networks confers self-organization properties on pseudocells

Background

Cell organization is governed and maintained via specific interactions among its constituent macromolecules. Comparison of the experimentally determined protein interaction networks in different model organisms has revealed little conservation of the specific edges linking ortholog proteins. Nevertheless, some topological characteristics of the graphs representing the networks - namely non-random degree distribution and high clustering coefficient - are shared by networks of distantly related organisms. Here we investigate the role of the topological features of the protein interaction network in promoting cell organization.

Methods

We have used a stochastic model, dubbed ProtNet representing a computer stylized cell to answer questions about the dynamic consequences of the topological properties of the static graphs representing protein interaction networks.

Results

By using a novel metrics of cell organization, we show that natural networks, differently from random networks, can promote cell self-organization. Furthermore the ensemble of protein complexes that forms in pseudocells, which self-organize according to the interaction rules of natural networks, are more robust to perturbations.

Conclusions

The analysis of the dynamic properties of networks with a variety of topological characteristics lead us to conclude that self organization is a consequence of the high clustering coefficient, whereas the scale free degree distribution has little influence on this property.

Aberrant TAL1 activation is mediated by an interchromosomal interaction in human T-cell acute lymphoblastic leukemia

Long-range chromatin interactions control metazoan gene transcription. However, the involvement of intra- and interchromosomal interactions in development and oncogenesis remains unclear. TAL1/SCL is a critical transcription factor required for the development of all hematopoietic lineages; yet, aberrant TAL1 transcription often occurs in T-cell acute lymphoblastic leukemia (T-ALL). Here, we report that oncogenic TAL1 expression is regulated by different intra- and interchromosomal loops in normal hematopoietic and leukemic cells, respectively. These intra- and interchromosomal loops alter the cell-type-specific enhancers that interact with the TAL1 promoter. We show that human SET1 (hSET1)-mediated H3K4 methylations promote a long-range chromatin loop, which brings the +51 enhancer in close proximity to TAL1 promoter 1 in erythroid cells. The CCCTC-binding factor (CTCF) facilitates this long-range enhancer/promoter interaction of the TAL1 locus in erythroid cells while blocking the same enhancer/promoter interaction of the TAL1 locus in human T-cell leukemia. In human T-ALL, a T-cell-specific transcription factor c-Maf-mediated interchromosomal interaction brings the TAL1 promoter into close proximity with a T-cell-specific regulatory element located on chromosome 16, activating aberrant TAL1 oncogene expression. Thus, our study reveals a novel molecular mechanism involving changes in three-dimensional chromatin interactions that activate the TAL1 oncogene in human T-cell leukemia.

Identifying and mapping cell-type-specific chromatin programming of gene expression

A problem of substantial interest is to systematically map variation in chromatin structure to gene-expression regulation across conditions, environments, or differentiated cell types. We developed and applied a quantitative framework for determining the existence, strength, and type of relationship between high-resolution chromatin structure in terms of DNaseI hypersensitivity and genome-wide gene-expression levels in 20 diverse human cell types. We show that ∼25% of genes show cell-type-specific expression explained by alterations in chromatin structure. We find that distal regions of chromatin structure (e.g., ±200 kb) capture more genes with this relationship than local regions (e.g., ±2.5 kb), yet the local regions show a more pronounced effect. By exploiting variation across cell types, we were capable of pinpointing the most likely hypersensitive sites related to cell-type-specific expression, which we show have a range of contextual uses. This quantitative framework is likely applicable to other settings aimed at relating continuous genomic measurements to gene-expression variation.

Transcriptional hierarchies regulating early blood cell development

Hematopoiesis represents one of the paradigmatic systems for studying stem cell biology, but our understanding of how the hematopoietic system develops during embryogenesis is still incomplete. While many lessons have been learned from studying the mouse embryo, embryonic stem cells have come to the fore as an alternative and more tractable model to recapitulate hematopoietic development. Here we review what is known about the embryonic origin of blood from these complementary systems and how transcription factor networks regulate the emergence of hematopoietic tissue from the mesoderm. Furthermore, we have performed an integrated analysis of genome-wide microarray and ChIP-seq data sets from mouse embryos and embryonic stem (ES) cell lines deficient in key regulators and demonstrate how this type of analysis can be used to reconstruct regulatory hierarchies that both confirm existing regulatory linkages and suggest additional interactions.

MAP17, a ROS-dependent oncogene

MAP17 is a small 17 kDa non-glycosylated membrane protein previously identified as being overexpressed in carcinomas. Breast tumor cells that overexpress MAP17 show an increased tumoral phenotype with enhanced proliferative capabilities both in the presence or the absence of contact inhibition, decreased apoptotic sensitivity, and increased migration. MAP17-expressing clones also grow better in nude mice. The increased malignant cell behavior induced by MAP17 is associated with an increase in reactive oxygen species (ROS) production, and the treatment of MAP17-expressing cells with antioxidants results in a reduction in the tumorigenic properties of these cells. The MAP17-dependent increase in ROS and tumorigenesis relies on its PDZ-binding domain because disruption of this sequence by point mutations abolishes the ability of MAP17 to enhance ROS production and tumorigenesis. MAP17 is overexpressed in a great variety of human carcinomas, including breast tumors. Immunohistochemical analysis of MAP17 during cancer progression demonstrates that overexpression of the protein strongly correlates with tumoral progression. Generalized MAP17 overexpression in human carcinomas indicates that MAP17 can be a good marker for tumorigenesis and, especially, for malignant progression.

MAP17 and the double-edged sword of ROS

Reactive oxygen species, ROS, are beneficially involved in many signaling pathways that control development and maintain cellular homeostasis. In physiological conditions, a tightly regulated redox balance protects cells from injurious ROS activity, but if the balance is altered, it promotes various pathological conditions including cancer. Understanding the duality of ROS as cytotoxic molecules and key mediators in signaling cascades may provide novel opportunities for improved cancer therapy. MAP17 is a small 17-kDa non-glycosylated membrane protein that is overexpressed in many tumors of different origins, including carcinomas. Immunohistochemical analysis of MAP17 during cancer progression demonstrates that overexpression of the protein strongly correlates with the progression of most types of tumor. Tumor cells that overexpress MAP17 show an increased tumoral phenotype associated with an increase in ROS. However, in non-tumor cells MAP17 increases ROS, resulting in senescence or apoptosis. Therefore, in tumor cells, MAP17 could be a marker for increased oxidative stress and could define new therapeutic approaches. Here, we review the role of MAP17 as a putative oncogene, as well as its role in cancer and anticancer therapies.

Using zebrafish (Danio rerio) to assess gene function in thrombus formation

Cardiovascular and cerebrovascular disease is the major cause of death in the developed world, with a high burden of disease and substantial pharmaceutical investment to manage it (WHO, Global Burden of Disease, 2004 Update, W.H. Organisation, Editor. 2008). Platelets, as the principal mediators of thrombus formation, are a primary pharmaceutical target, with attenuation of platelet function and thrombus formation significantly reducing the incidence of myocardial infarction and stroke. Haemostasis, however, may also be affected by antithrombotics, leading to spontaneous and/or prolonged bleeding as a potentially severe side effect. Developing a comprehensive understanding of the mechanisms involved in platelet function and thrombus formation is anticipated to identify drug targets that may effectively manage vascular disease without an impact on haemostasis. Despite the progress in characterising individual genes in platelet function and thrombosis, using gene knockout and transgenic mice over the past decade or so, there is still much to be uncovered. Investigating gene function using mouse models is a substantial investment and a considerable amount of work, with a relevant phenotype not guaranteed. As such, a new model is needed for the effective screening of novel genes that have been identified as having potential roles in platelet function or cardiovascular disease by genomic association and comparative expression studies (Nature, 447(7145): 661-678, 2007; Nat Genet, 41(11): 1182-1190, 2009; N Engl J Med, 357(5): 443-453, 2007; Blood, 109(8): 3260-3269, 2007). Here, we highlight and discuss the relevance of the zebrafish (Danio rerio) as a model for studying thrombosis, the current techniques that are employed to assess gene function in a zebrafish model of thrombosis, and how an effective genetic screen may be constructed.

Protein kinase C alpha and beta are positive regulators of thrombus formation in vivo in a zebrafish (Danio rerio) model of thrombosis

BACKGROUND

The zebrafish (Danio rerio) is becoming an attractive model organism for the assessment of gene function in thrombosis in vivo. Zebrafish, as a thrombosis model, have several advantages, with the capacity to follow thrombus formation at high resolution in real time using intravital microscopy, without the need for complex surgical techniques, and the capability to rapidly knockdown gene expression using morpholino antisense approaches.

OBJECTIVES

We have recently shown, in mouse models, that protein kinase C alpha (PKCα) plays a critical role in regulating thrombus formation in vivo. PKC beta (β) plays a non-redundant role also in platelet function in vitro, but the function of this gene had not yet been assessed in vivo.

METHODS

In the present study, we analyzed the function of both PKCα and PKCβ in the zebrafish model in vivo, by live imaging using a laser-induced injury of the main caudal artery in 3-day-old larvae.

RESULTS

We showed that D. rerio express orthologs of both the PKCα and PKCβ genes, with high sequence identity. Translation blocking and splice-blocking morpholinos effectively and specifically knockdown expression of these genes and knockdown with either morpholino leads to attenuated thrombus formation, as assessed by several quantitative parameters including time to initial adhesion and peak thrombus surface area.

CONCLUSIONS

Our data indicate that these two highly related genes play non-redundant roles in regulating thrombosis, an observation that supports our previous in vitro murine data, and suggests unique roles, and possibly unique regulation, for PKCα and PKCβ in controlling platelet function in vivo.

Genome-wide analysis of simultaneous GATA1/2, RUNX1, FLI1, and SCL binding in megakaryocytes identifies hematopoietic regulators

Hematopoietic differentiation critically depends on combinations of transcriptional regulators controlling the development of individual lineages. Here, we report the genome-wide binding sites for the five key hematopoietic transcription factors--GATA1, GATA2, RUNX1, FLI1, and TAL1/SCL--in primary human megakaryocytes. Statistical analysis of the 17,263 regions bound by at least one factor demonstrated that simultaneous binding by all five factors was the most enriched pattern and often occurred near known hematopoietic regulators. Eight genes not previously appreciated to function in hematopoiesis that were bound by all five factors were shown to be essential for thrombocyte and/or erythroid development in zebrafish. Moreover, one of these genes encoding the PDZK1IP1 protein shared transcriptional enhancer elements with the blood stem cell regulator TAL1/SCL. Multifactor ChIP-Seq analysis in primary human cells coupled with a high-throughput in vivo perturbation screen therefore offers a powerful strategy to identify essential regulators of complex mammalian differentiation processes.

cis-Regulatory remodeling of the SCL locus during vertebrate evolution

Development progresses through a sequence of cellular identities which are determined by the activities of networks of transcription factor genes. Alterations in cis-regulatory elements of these genes play a major role in evolutionary change, but little is known about the mechanisms responsible for maintaining conserved patterns of gene expression. We have studied the evolution of cis-regulatory mechanisms controlling the SCL gene, which encodes a key transcriptional regulator of blood, vasculature, and brain development and exhibits conserved function and pattern of expression throughout vertebrate evolution. SCL cis-regulatory elements are conserved between frog and chicken but accrued alterations at an accelerated rate between 310 and 200 million years ago, with subsequent fixation of a new cis-regulatory pattern at the beginning of the mammalian radiation. As a consequence, orthologous elements shared by mammals and lower vertebrates exhibit functional differences and binding site turnover between widely separated cis-regulatory modules. However, the net effect of these alterations is constancy of overall regulatory inputs and of expression pattern. Our data demonstrate remarkable cis-regulatory remodelling across the SCL locus and indicate that stable patterns of expression can mask extensive regulatory change. These insights illuminate our understanding of vertebrate evolution.

Genomic approaches uncover increasing complexities in the regulatory landscape at the human SCL (TAL1) locus

The SCL (TAL1) transcription factor is a critical regulator of haematopoiesis and its expression is tightly controlled by multiple cis-acting regulatory elements. To elaborate further the DNA elements which control its regulation, we used genomic tiling microarrays covering 256 kb of the human SCL locus to perform a concerted analysis of chromatin structure and binding of regulatory proteins in human haematopoietic cell lines. This approach allowed us to characterise further or redefine known human SCL regulatory elements and led to the identification of six novel elements with putative regulatory function both up and downstream of the SCL gene. They bind a number of haematopoietic transcription factors (GATA1, E2A LMO2, SCL, LDB1), CTCF or components of the transcriptional machinery and are associated with relevant histone modifications, accessible chromatin and low nucleosomal density. Functional characterisation shows that these novel elements are able to enhance or repress SCL promoter activity, have endogenous promoter function or enhancer-blocking insulator function. Our analysis opens up several areas for further investigation and adds new layers of complexity to our understanding of the regulation of SCL expression.

Ontogeny of haematopoiesis: recent advances and open questions

Unravelling the embryonic origins of the haematopoietic system has been the subject of sustained research for more than a century. Nevertheless, many important questions are still either unanswered or remain a matter of intense debate. Recent progress in mouse and embryonic stem cell model systems as well as imaging and post-genomic technologies has provided new insights into many of these open questions. Here we place into context recent reports on the anatomical site of blood stem cell emergence and, using red blood cells as an example, illustrate how the development of stem cells and the other blood lineages is both temporally and spatially decoupled. In addition, we outline how embryonic stem cell assays are increasingly used as a powerful surrogate for studying lineage relationships and developmental potential of early embryonic blood progenitors. Finally, we review how recent progress in the reconstruction of transcriptional regulatory networks is beginning to define the connectivity between key regulators that control early blood development. In light of these rapid recent advances, research into the embryonic origins of the haematopoietic system should remain one of the most vibrant disciplines within the wider field of haematology for the foreseeable future.

Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins

In this study, we demonstrate the suitability of the vertebrate Danio rerio (zebrafish) for functional screening of novel platelet genes in vivo by reverse genetics. Comparative transcript analysis of platelets and their precursor cell, the megakaryocyte, together with nucleated blood cell elements, endothelial cells, and erythroblasts, identified novel platelet membrane proteins with hitherto unknown roles in thrombus formation. We determined the phenotype induced by antisense morpholino oligonucleotide (MO)-based knockdown of 5 of these genes in a laser-induced arterial thrombosis model. To validate the model, the genes for platelet glycoprotein (GP) IIb and the coagulation protein factor VIII were targeted. MO-injected fish showed normal thrombus initiation but severely impaired thrombus growth, consistent with the mouse knockout phenotypes, and concomitant knockdown of both resulted in spontaneous bleeding. Knockdown of 4 of the 5 novel platelet proteins altered arterial thrombosis, as demonstrated by modified kinetics of thrombus initiation and/or development. We identified a putative role for BAMBI and LRRC32 in promotion and DCBLD2 and ESAM in inhibition of thrombus formation. We conclude that phenotypic analysis of MO-injected zebrafish is a fast and powerful method for initial screening of novel platelet proteins for function in thrombosis.

A novel mode of enhancer evolution: the Tal1 stem cell enhancer recruited a MIR element to specifically boost its activity

Altered cis-regulation is thought to underpin much of metazoan evolution, yet the underlying mechanisms remain largely obscure. The stem cell leukemia TAL1 (also known as SCL) transcription factor is essential for the normal development of blood stem cells and we have previously shown that the Tal1 +19 enhancer directs expression to hematopoietic stem cells, hematopoietic progenitors, and to endothelium. Here we demonstrate that an adjacent region 1 kb upstream (+18 element) is in an open chromatin configuration and carries active histone marks but does not function as an enhancer in transgenic mice. Instead, it boosts activity of the +19 enhancer both in stable transfection assays and during differentiation of embryonic stem (ES) cells carrying single-copy reporter constructs targeted to the Hprt locus. The +18 element contains a mammalian interspersed repeat (MIR) which is essential for the +18 function and which was transposed to the Tal1 locus approximately 160 million years ago at the time of the mammalian/marsupial branchpoint. Our data demonstrate a previously unrecognized mechanism whereby enhancer activity is modulated by a transposon exerting a "booster" function which would go undetected by conventional transgenic approaches.

Qualifying the relationship between sequence conservation and molecular function

Quantification of evolutionary constraints via sequence conservation can be leveraged to annotate genomic functional sequences. Recent efforts addressing the converse of this relationship have identified many sites in metazoan genomes with molecular function but without detectable conservation between related species. Here, we discuss explanations and implications for these results considering both practical and theoretical issues. In particular, phylogenetic scope influences the relationship between sequence conservation and function. Comparisons of distantly related species can detect constraint with high specificity due to the loss of conserved neutral sequence, but sensitivity is sacrificed as a result of functional changes related to lineage-specific biology. The strength of natural selection operating on functional sequence is also important. Mutations to functional sequences that result in small fitness effects are subject to weaker constraints. Therefore, particularly when comparing highly divergent species, functional sequences that are degenerate or biologically redundant will be prone to turnover, wherein functional sequences are replaced by effectively equivalent, but nonorthologous counterparts. Finally, considering the size and complexity of metazoan genomes and the fact that many nonconserved sequences are associated with sequence-degenerate, low-level molecular functions, we find it likely that there exist many biochemically functional sequences that are not under constraint. This hypothesis does not lead to the conclusion that huge amounts of vertebrate genomes are functionally important, but rather that such "functionality" represents molecular noise that has weak or no effect on organismal phenotypes.

Comparative genomics for detecting human disease genes

Originally, comparative genomics was geared toward defining the synteny of genes between species. As the human genome project accelerated, there was an increase in the number of tools and means to make comparisons culminating in having the genomic sequence for a large number of organisms spanning the evolutionary tree. With this level of resolution and a long history of comparative biology and comparative genetics, it is now possible to use comparative genomics to build or select better animal models and to facilitate gene discovery. Comparative genomics takes advantage of the functional genetic information from other organisms, (vertebrates and invertebrates), to apply it to the study of human physiology and disease. It allows for the identification of genes and regulatory regions, and for acquiring knowledge about gene function. In this chapter, the current state of comparative genomics and the available tools are discussed in the context of developing animal model systems that reflect the clinical picture.

Computational approaches to finding and analyzing cis-regulatory elements

The transcription of almost all developmental genes is driven by tissue- and time-specific regulatory elements. These transcriptional regulatory elements lie in the genomic DNA proximal to the gene, and hence are cis-regulatory (as opposed to trans-regulatory elements like transcription factor genes). Over the past three decades, a number of techniques have been applied to the problem of finding and characterizing these regulatory elements. In this chapter, I discuss some computational approaches that have been particularly useful in identifying developmental cis-regulatory regions, and provide a tutorial on how to apply these approaches to the study of chick development.

Distinct functions for different scl isoforms in zebrafish primitive and definitive hematopoiesis

The stem-cell leukemia (SCL, also known as TAL1) gene encodes a basic helix-loop-helix transcription factor that is essential for the initiation of primitive and definitive hematopoiesis, erythrocyte and megakarocyte differentiation, angiogenesis, and astrocyte development. Here we report that the zebrafish produces, through an alternative promoter site, a novel truncated scl (tal1) isoform, scl-beta, which manifests a temporal and spatial expression distinct from the previously described full-length scl-alpha. Functional analysis reveals that while scl-alpha and -beta are redundant for the initiation of primitive hematopoiesis, these two isoforms exert distinct functions in the regulation of primitive erythroid differentiation and definitive hematopoietic stem cell specification. We further demonstrate that differences in the protein expression levels of scl-alpha and -beta, by regulating their protein stability, are likely to give rise to their distinct functions. Our findings suggest that hematopoietic cells at different levels of hierarchy are likely governed by a gradient of the Scl protein established through temporal and spatial patterns of expression of the different isoforms.

MAP17 overexpression is a common characteristic of carcinomas

We undertook a large-scale genetic screen to identify genes able to alter the cellular response to physiological signals and provide selective advantage once tumorigenesis has begun. We identified MAP17, a small 17 kDa non-glycosylated membrane protein previously identified, being overexpressed in carcinomas. We found that MAP17 is overexpressed in a great variety of human carcinomas. Immunohistochemical analysis of MAP17 during cancer progression shows, at least in prostate and ovarian carcinomas, that overexpression of the protein strongly correlates with tumoral progression (P < 0.0001). Many tumor cells also express MAP17 and its expression does not correlate with expression of SCL, a neighbor gene reported to be co-expressed in some hematopoietic cell lines. SCL neither is expressed in most MAP17-positive tumors, indicating the independent transcription of MAP17, at least in carcinomas. We cloned 5' genomic region to MAP17 and described the minimal promoter necessary to produce independent activation of MAP17. Moreover, we have found that MAP17 promoter is activated by oncogenes. Taken together, our data show an independent activation of MAP17 promoter that can be driven by oncogenes and that might explain the common overexpression of MAP17 in human carcinomas.

SCL-GFP transgenic zebrafish: in vivo imaging of blood and endothelial development and identification of the initial site of definitive hematopoiesis

The bHLH transcription factor SCL plays a central role in the generation of hematopoietic cells in vertebrates. We modified a PAC containing the whole zebrafish scl locus, inserting GFP into the first coding exon of scl. In germline-transgenic zebrafish generated using this construct, GFP expression completely recapitulates the endogenous expression of scl in blood, endothelium and CNS. We performed in vivo timelapse imaging of blood and endothelial precursor migration at the single-cell level and show that these cells migrate from the posterior lateral plate mesoderm to their site of differentiation in the intermediate cell mass. The anterior lateral plate domain of GFP expression gives rise to primitive macrophages and the blood vessels of the head. In later embryos, GFP expression identifies clusters of hematopoietic cells that develop between the dorsal aorta and posterior cardinal veins after primitive erythrocytes have entered circulation. Two treatments that block definitive hematopoiesis (treatment with dioxin (TCDD), and injection of an antisense morpholino oligonucleotide targeted to runx1) ablate these hematopoietic clusters. This indicates that these clusters represent the first site of definitive hematopoiesis in zebrafish. This site is anatomically homologous to the proposed source of hematopoietic stem cells in amniotes, the aorta-gonad-mesonephros (AGM) region. A second transgenic line, containing the promoter of scl driving GFP, lacks expression in the definitive clusters.

Multiple sequence alignment: in pursuit of homologous DNA positions

DNA sequence alignment is a prerequisite to virtually all comparative genomic analyses, including the identification of conserved sequence motifs, estimation of evolutionary divergence between sequences, and inference of historical relationships among genes and species. While it is mere common sense that inaccuracies in multiple sequence alignments can have detrimental effects on downstream analyses, it is important to know the extent to which the inferences drawn from these alignments are robust to errors and biases inherent in all sequence alignments. A survey of investigations into strengths and weaknesses of sequence alignments reveals, as expected, that alignment quality is generally poor for two distantly related sequences and can often be improved by adding additional sequences as stepping stones between distantly related species. Errors in sequence alignment are also found to have a significant negative effect on subsequent inference of sequence divergence, phylogenetic trees, and conserved motifs. However, our understanding of alignment biases remains rudimentary, and sequence alignment procedures continue to be used somewhat like benign formatting operations to make sequences equal in length. Because of the central role these alignments now play in our endeavors to establish the tree of life and to identify important parts of genomes through evolutionary functional genomics, we see a need for increased community effort to investigate influences of alignment bias on the accuracy of large-scale comparative genomics.

Real-time PCR mapping of DNaseI-hypersensitive sites using a novel ligation-mediated amplification technique

Mapping sites within the genome that are hypersensitive to digestion with DNaseI is an important method for identifying DNA elements that regulate transcription. The standard approach to locating these DNaseI-hypersensitive sites (DHSs) has been to use Southern blotting techniques, although we, and others, have recently published alternative methods using a range of technologies including high-throughput sequencing and genomic array tiling paths. In this article, we describe a novel protocol to use real-time PCR to map DHS. Advantages of the technique reported here include the small cell numbers required for each analysis, rapid, relatively low-cost experiments with minimal need for specialist equipment. Presented examples include comparative DHS mapping of known TAL1/SCL regulatory elements between human embryonic stem cells and K562 cells.

Transcription factors in myeloid development: balancing differentiation with transformation

In recent years, great progress has been made in elucidating the progenitor-cell hierarchy of the myeloid lineage. Transcription factors have been shown to be key determinants in the orchestration of myeloid identity and differentiation fates. Most transcription factors show cell-lineage-restricted and stage-restricted expression patterns, indicating the requirement for tight regulation of their activities. Moreover, if dysregulated or mutated, these transcription factors cause the differentiation block observed in many myeloid leukaemias. Consequently, therapies designed to restore defective transcription factor functions are an attractive option in the treatment of myeloid and other human cancers.

WeederH: an algorithm for finding conserved regulatory motifs and regions in homologous sequences

Background

This work addresses the problem of detecting conserved transcription factor binding sites and in general regulatory regions through the analysis of sequences from homologous genes, an approach that is becoming more and more widely used given the ever increasing amount of genomic data available.

Results

We present an algorithm that identifies conserved transcription factor binding sites in a given sequence by comparing it to one or more homologs, adapting a framework we previously introduced for the discovery of sites in sequences from co-regulated genes. Differently from the most commonly used methods, the approach we present does not need or compute an alignment of the sequences investigated, nor resorts to descriptors of the binding specificity of known transcription factors. The main novel idea we introduce is a relative measure of conservation, assuming that true functional elements should present a higher level of conservation with respect to the rest of the sequence surrounding them. We present tests where we applied the algorithm to the identification of conserved annotated sites in homologous promoters, as well as in distal regions like enhancers.

Conclusion

Results of the tests show how the algorithm can provide fast and reliable predictions of conserved transcription factor binding sites regulating the transcription of a gene, with better performances than other available methods for the same task. We also show examples on how the algorithm can be successfully employed when promoter annotations of the genes investigated are missing, or when regulatory sites and regions are located far away from the genes.

Comparative genomics: a tool to functionally annotate human DNA

The availability of an increasing number of vertebrate genomes has enabled comparative methods to infer functional sequences based on evolutionary constraint. Although this has proved powerful for gene identification, significant progress has also been made in uncovering gene regulatory sequences such as distant acting transcriptional enhancers. These pursuits have led to the development of a variety of valuable databases and resources that should serve as a routine toolbox for biological discovery.

Alignment of genomic sequences using DIALIGN

DIALIGN is a software program for multiple alignment of DNA or protein sequences that combines global and local alignment features. During the last years, the program has been used extensively to compare syntenic regions in genomic sequences. An anchoring option speeds up the alignment procedure and makes it possible to use user-defined constraints to improve the quality of the program output. This chapter explains features of DIALIGN that are useful if genomic sequences are to be aligned. The program is online available through Göttingen Bioinformatics Compute Server at http://dialign.gobics.de/.

Identifying gene regulatory elements by genomic microarray mapping of DNaseI hypersensitive sites

The identification of cis-regulatory elements is central to understanding gene transcription. Hypersensitivity of cis-regulatory elements to digestion with DNaseI remains the gold-standard approach to locating such elements. Traditional methods used to identify DNaseI hypersensitive sites are cumbersome and can only be applied to short stretches of DNA at defined locations. Here we report the development of a novel genomic array-based approach to DNaseI hypersensitive site mapping (ADHM) that permits precise, large-scale identification of such sites from as few as 5 million cells. Using ADHM we identified all previously recognized hematopoietic regulatory elements across 200 kb of the mouse T-cell acute lymphocytic leukemia-1 (Tal1) locus, and, in addition, identified two novel elements within the locus, which show transcriptional regulatory activity. We further validated the ADHM protocol by mapping the DNaseI hypersensitive sites across 250 kb of the human TAL1 locus in CD34+ primary stem/progenitor cells and K562 cells and by mapping the previously known DNaseI hypersensitive sites across 240 kb of the human alpha-globin locus in K562 cells. ADHM provides a powerful approach to identifying DNaseI hypersensitive sites across large genomic regions.

SinicView: a visualization environment for comparisons of multiple nucleotide sequence alignment tools

Background

Deluged by the rate and complexity of completed genomic sequences, the need to align longer sequences becomes more urgent, and many more tools have thus been developed. In the initial stage of genomic sequence analysis, a biologist is usually faced with the questions of how to choose the best tool to align sequences of interest and how to analyze and visualize the alignment results, and then with the question of whether poorly aligned regions produced by the tool are indeed not homologous or are just results due to inappropriate alignment tools or scoring systems used. Although several systematic evaluations of multiple sequence alignment (MSA) programs have been proposed, they may not provide a standard-bearer for most biologists because those poorly aligned regions in these evaluations are never discussed. Thus, a tool that allows cross comparison of the alignment results obtained by different tools simultaneously could help a biologist evaluate their correctness and accuracy.

Results

In this paper, we present a versatile alignment visualization system, called SinicView, (for Sequence-aligning INnovative and Interactive Comparison VIEWer), which allows the user to efficiently compare and evaluate assorted nucleotide alignment results obtained by different tools. SinicView calculates similarity of the alignment outputs under a fixed window using the sum-of-pairs method and provides scoring profiles of each set of aligned sequences. The user can visually compare alignment results either in graphic scoring profiles or in plain text format of the aligned nucleotides along with the annotations information. We illustrate the capabilities of our visualization system by comparing alignment results obtained by MLAGAN, MAVID, and MULTIZ, respectively.

Conclusion

With SinicView, users can use their own data sequences to compare various alignment tools or scoring systems and select the most suitable one to perform alignment in the initial stage of sequence analysis.

Many genes in fish have species-specific asymmetric rates of molecular evolution

BACKGROUND

Gene and genome duplication events increase the amount of genetic material that might then contribute to an increase in the genomic and phenotypic complexity of organisms during evolution. Thus, it has been argued that there is a relationship between gene copy number and morphological complexity and/or species diversity. This hypothesis implies that duplicated genes have subdivided or evolved novel functions compared to their pre-duplication proto-orthologs. Such a functional divergence might be caused by an increase in evolutionary rates in one ortholog, by changes in expression, regulatory evolution, insertion of repetitive elements, or due to positive Darwinian selection in one copy. We studied a set of 2466 genes that were present in Danio rerio, Takifugu rubripes, Tetraodon nigroviridis and Oryzias latipes to test (i) for forces of positive Darwinian selection; (ii) how frequently duplicated genes are retained, and (iii) whether novel gene functions might have evolved.

RESULTS

25% (610) of all investigated genes show significantly smaller or higher genetic distances in the genomes of particular fish species compared to their human ortholog than their orthologs in other fish according to relative rate tests. We identified 49 new paralogous pairs of duplicated genes in fish, in which one of the paralogs is under positive Darwinian selection and shows a significantly higher rate of molecular evolution in one of the four fish species, whereas the other copy apparently did not undergo adaptive changes since it retained the original rate of evolution. Among the genes under positive Darwinian selection, we found a surprisingly high number of ATP binding proteins and transcription factors.

CONCLUSION

The significant rate difference suggests that the function of these rate-changed genes might be essential for the respective fish species. We demonstrate that the measurement of positive selection is a powerful tool to identify divergence rates of duplicated genes and that this method has the capacity to identify potentially interesting candidates for adaptive gene evolution.

A novel approach to identifying regulatory motifs in distantly related genomes

A two-step procedure for identifying regulatory motifs in distantly related organisms is described that combines the advantages of sequence alignment and motif detection approaches.

Although proven successful in the identification of regulatory motifs, phylogenetic footprinting methods still show some shortcomings. To assess these difficulties, most apparent when applying phylogenetic footprinting to distantly related organisms, we developed a two-step procedure that combines the advantages of sequence alignment and motif detection approaches. The results on well-studied benchmark datasets indicate that the presented method outperforms other methods when the sequences become either too long or too heterogeneous in size.

Trade-offs in detecting evolutionarily constrained sequence by comparative genomics

As whole-genome sequencing efforts extend beyond more traditional model organisms to include a deep diversity of species, comparative genomic analyses will be further empowered to reveal insights into the human genome and its evolution. The discovery and annotation of functional genomic elements is a necessary step toward a detailed understanding of our biology, and sequence comparisons have proven to be an integral tool for that task. This review is structured to broadly reflect the statistical challenges in discriminating these functional elements from the bulk of the genome that has evolved neutrally. Specifically, we review the comparative genomics literature in terms of specificity, sensitivity, and phylogenetic scope, as well as the trade-offs that relate these factors in standard analyses. We consider the impact of an expanding diversity of orthologous sequences on our ability to resolve functional elements. This impact is assessed through both recent comparative analyses of deep alignments and mathematical modeling.

Distinct Roles for SCL in Erythroid Specification and Maturation in Zebrafish

The stem cell leukemia (SCL) transcription factor is essential for vertebrate hematopoiesis. Using the powerful zebrafish model for embryonic analysis, we compared the effects of either reducing or ablating Scl using morpholino-modified antisense RNAs. Ablation of Scl resulted in the loss of primitive and definitive hematopoiesis, consistent with its essential role in these processes. Interestingly, in embryos with severely reduced Scl levels, erythroid progenitors expressing gata1 and embryonic globin developed. Erythroid maturation was deficient in these Scl hypomorphs, supporting that Scl was required both for the erythroid specification and for the maturation steps, with maturation requiring higher Scl levels than specification. Although all hematopoietic functions were rescued by wild-type Scl mRNA, an Scl DNA binding mutant rescued primitive and definitive hematopoiesis but did not rescue primitive erythroid maturation. Together, we showed that there is a distinct Scl hypomorphic phenotype and demonstrated that distinct functions are required for the roles of Scl in the specification and differentiation of primitive and definitive hematopoietic lineages. Our results revealed that Scl participates in multiple processes requiring different levels and functions. Further, we identified an Scl hypomorphic phenotype distinct from the null state.

Large multiple organism gene finding by collapsed Gibbs sampling

The Gibbs sampling method has been widely used for sequence analysis after it was successfully applied to the problem of identifying regulatory motif sequences upstream of genes. Since then, numerous variants of the original idea have emerged: however, in all cases the application has been to finding short motifs in collections of short sequences (typically less than 100 nucleotides long). In this paper, we introduce a Gibbs sampling approach for identifying genes in multiple large genomic sequences up to hundreds of kilobases long. This approach leverages the evolutionary relationships between the sequences to improve the gene predictions, without explicitly aligning the sequences. We have applied our method to the analysis of genomic sequence from 14 genomic regions, totaling roughly 1.8 Mb of sequence in each organism. We show that our approach compares favorably with existing ab initio approaches to gene finding, including pairwise comparison based gene prediction methods which make explicit use of alignments. Furthermore, excellent performance can be obtained with as little as four organisms, and the method overcomes a number of difficulties of previous comparison based gene finding approaches: it is robust with respect to genomic rearrangements, can work with draft sequence, and is fast (linear in the number and length of the sequences). It can also be seamlessly integrated with Gibbs sampling motif detection methods.

Selecting transpositions using phage P1 headful packaging: new markerless transposons for functionally mapping long-range regulatory sequences in bacterial artificial chromosomes and P1-derived artificial chromosomes

New Tn10 minitransposons were constructed to functionally map long-range transcription regulatory sequences in bacterial artificial chromosomes (BACs) and P1-derived artificial chromosomes (PACs). Each contained a wild-type loxP site but, significantly, contained no mammalian or bacterial genes and/or promoter elements within the transposed portion of DNA. In contrast to loxP transposons described previously, the new ones do not introduce transcription regulatory elements capable of interfering with those endogenous to the BAC clone in functional mapping studies. Progressive deletions from the loxP end of genomic DNA were efficiently generated using these transposons, and a series of truncations generated in a green fluorescence protein (GFP)-BAC fusion clone unambiguously identified three new long-range enhancer sequences functionally in the Nkx2-5 gene in transgenic mice. Insertions of these new transposons lacking antibiotic resistance genes into a BAC or PAC were indirectly selected by their ability to delete enough DNA from the clone so as to enable its packaging within a P1 phage head with both loxP sites intact for subsequent recovery of the large plasmid. The outcome of such an indirect mode of selection is both desirable and undesirable. First, because the screen is not antibiotic resistance marker dependent, the same transposon can be used to generate nested deletions efficiently in both BACs and PACs. Second, deletions through intrainsert recombinations unrelated to loxP/Cre also get packaged and recovered, and size analyses of the BAC/PAC vector band after NotI digestion is indispensable to identify authentic loxP/Cre deletions. The procedure nevertheless offers a potential approach to map recombinogenic sequences in BACs and PACs.

Identification of genetic polymorphisms through comparative DNA sequence analysis on the K-ras gene: implications for lung tumor susceptibility

In the present study, the authors performed a comparative sequence analysis of the K-ras gene. By comparing sequences from different mouse inbred strains, the authors have identified new nucleotide polymorphisms in the noncoding regions of mouse K-ras gene. They have also identified noncoding DNA segments evolutionarily conserved among the human, mouse, and rat. Computational analysis for transcription factor binding sites suggests that these polymorphic and conserved DNA sequences harbor potential cis-regulatory elements, which may contribute to the transcriptional regulation of the K-ras gene. Further studies on these potential regulatory sites may help to elucidate the fundamental mechanism underlying allele-specific activation and expression of K-ras gene in hybrid mouse lung tumors, which determines lung tumor susceptibility in mice.

Hematopoietic differentiation of human embryonic stem cells progresses through sequential hematoendothelial, primitive, and definitive stages resembling human yolk sac development

We elucidate the cellular and molecular kinetics of the stepwise differentiation of human embryonic stem cells (hESCs) to primitive and definitive erythromyelopoiesis from human embryoid bodies (hEBs) in serum-free clonogenic assays. Hematopoiesis initiates from CD45 hEB cells with emergence of semiadherent mesodermal-hematoendothelial (MHE) colonies that can generate endothelium and form organized, yolk sac-like structures that secondarily generate multipotent primitive hematopoietic stem progenitor cells (HSPCs), erythroblasts, and CD13+CD45+ macrophages. A first wave of hematopoiesis follows MHE colony emergence and is predominated by primitive erythropoiesis characterized by a brilliant red hemoglobinization, CD71/CD325a (glycophorin A) expression, and exclusively embryonic/fetal hemoglobin expression. A second wave of definitive-type erythroid burst-forming units (BFU-e's), erythroid colony-forming units (CFU-e's), granulocyte-macrophage colony-forming cells (GM-CFCs), and multilineage CFCs follows next from hEB progenitors. These stages of hematopoiesis proceed spontaneously from hEB-derived cells without requirement for supplemental growth factors during hEB differentiation. Gene expression analysis of differentiating hEBs revealed that initiation of hematopoiesis correlated with increased levels of SCL/TAL1, GATA1, GATA2, CD34, CD31, and the homeobox gene-regulating factor CDX4 These data indicate that hematopoietic differentiation of hESCs models the earliest events of embryonic and definitive hematopoiesis in a manner resembling human yolk sac development, thus providing a valuable tool for dissecting the earliest events in human HSPC genesis.

Fli1, Elf1, and Ets1 regulate the proximal promoter of the LMO2 gene in endothelial cells

Transcriptional control has been identified as a key mechanism regulating the formation and subsequent behavior of hematopoietic stem cells. We have used a comparative genomics approach to identify transcriptional regulatory elements of the LMO2 gene, a transcriptional cofactor originally identified through its involvement in T-cell leukemia and subsequently shown to be critical for normal hematopoietic and endothelial development. Of the 2 previously characterized LMO2 promoters, the second (proximal) promoter was highly conserved in vertebrates ranging from mammals to fish. Real-time reverse transcriptase-polymerase chain reaction (RT-PCR) expression analysis identified this promoter as the predominant source of transcription in hematopoietic tissue. Transient and stable transfections indicated that the proximal promoter was active in hematopoietic progenitor and endothelial cell lines and this activity was shown to depend on 3 conserved Ets sites that were bound in vivo by E74-like factor 1 (Elf1), Friend leukemia integration 1 (Fli1), and erythroblastosis virus oncogene homolog E twenty-six-1 (Ets1). Finally, transgenic analysis demonstrated that the LMO2 proximal promoter is sufficient for expression in endothelial cells in vivo. No hematopoietic expression was observed, indicating that additional enhancers are required to mediate transcription from the proximal promoter in hematopoietic cells. Together, these results suggest that the conserved proximal promoter is central to LMO2 transcription in hematopoietic and endothelial cells, where it is regulated by Ets factors.

Distribution and intensity of constraint in mammalian genomic sequence

Comparisons of orthologous genomic DNA sequences can be used to characterize regions that have been subject to purifying selection and are enriched for functional elements. We here present the results of such an analysis on an alignment of sequences from 29 mammalian species. The alignment captures approximately 3.9 neutral substitutions per site and spans approximately 1.9 Mbp of the human genome. We identify constrained elements from 3 bp to over 1 kbp in length, covering approximately 5.5% of the human locus. Our estimate for the total amount of nonexonic constraint experienced by this locus is roughly twice that for exonic constraint. Constrained elements tend to cluster, and we identify large constrained regions that correspond well with known functional elements. While constraint density inversely correlates with mobile element density, we also show the presence of unambiguously constrained elements overlapping mammalian ancestral repeats. In addition, we describe a number of elements in this region that have undergone intense purifying selection throughout mammalian evolution, and we show that these important elements are more numerous than previously thought. These results were obtained with Genomic Evolutionary Rate Profiling (GERP), a statistically rigorous and biologically transparent framework for constrained element identification. GERP identifies regions at high resolution that exhibit nucleotide substitution deficits, and measures these deficits as "rejected substitutions". Rejected substitutions reflect the intensity of past purifying selection and are used to rank and characterize constrained elements. We anticipate that GERP and the types of analyses it facilitates will provide further insights and improved annotation for the human genome as mammalian genome sequence data become richer.

Transcriptional regulation of the SCL locus: identification of an enhancer that targets the primitive erythroid lineage in vivo

The stem cell leukemia (SCL) gene, also known as TAL-1, encodes a basic helix-loop-helix protein that is essential for the formation of all hematopoietic lineages, including primitive erythropoiesis. Appropriate transcriptional regulation is essential for the biological functions of SCL, and we have previously identified five distinct enhancers which target different subdomains of the normal SCL expression pattern. However, it is not known whether these SCL enhancers also regulate neighboring genes within the SCL locus, and the erythroid expression of SCL remains unexplained. Here, we have quantitated transcripts from SCL and neighboring genes in multiple hematopoietic cell types. Our results show striking coexpression of SCL and its immediate downstream neighbor, MAP17, suggesting that they share regulatory elements. A systematic survey of histone H3 and H4 acetylation throughout the SCL locus in different hematopoietic cell types identified several peaks of histone acetylation between SIL and MAP17, all of which corresponded to previously characterized SCL enhancers or to the MAP17 promoter. Downstream of MAP17 (and 40 kb downstream of SCL exon 1a), an additional peak of acetylation was identified in hematopoietic cells and was found to correlate with expression of SCL but not other neighboring genes. This +40 region is conserved in human-dog-mouse-rat sequence comparisons, functions as an erythroid cell-restricted enhancer in vitro, and directs beta-galactosidase expression to primitive, but not definitive, erythroblasts in transgenic mice. The SCL +40 enhancer provides a powerful tool for studying the molecular and cellular biology of the primitive erythroid lineage.

Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements

With the imminent completion of the whole genome sequence of humans, increasing attention is being focused on the annotation of cis-regulatory elements in the human genome. Comparative genomics approaches based on evolutionary conservation have proved useful in the detection of conserved cis-regulatory elements. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparative genomics, by virtue of its compact genome and maximal phylogenetic distance from mammals. Fugu has lost a large proportion of nonessential DNA, and retained single orthologs for many duplicate genes that arose in the fish lineage. Non-coding sequences conserved between fugu and mammals have been shown to be functional cis-regulatory elements. Thus, fugu is a model fish genome of choice for discovering evolutionarily conserved regulatory elements in the human genome. Such evolutionarily conserved elements are likely to be shared by all vertebrates, and related to regulatory interactions fundamental to all vertebrates. The functions of these conserved vertebrate elements can be rapidly assayed in mammalian cell lines or in transgenic systems such as zebrafish/medaka and Xenopus, followed by validation of crucial elements in transgenic rodents.

Long-range comparison of human and mouse Sprr loci to identify conserved noncoding sequences involved in coordinate regulation

Mammalian epidermis provides a permeability barrier between an organism and its environment. Under homeostatic conditions, epidermal cells produce structural proteins, which are cross-linked in an orderly fashion to form a cornified envelope (CE). However, under genetic or environmental stress, specific genes are induced to rapidly build a temporary barrier. Small proline-rich (SPRR) proteins are the primary constituents of the CE. Under stress the entire family of 14 Sprr genes is upregulated. The Sprr genes are clustered within the larger epidermal differentiation complex on mouse chromosome 3, human chromosome 1q21. The clustering of the Sprr genes and their upregulation under stress suggest that these genes may be coordinately regulated. To identify enhancer elements that regulate this stress response activation of the Sprr locus, we utilized bioinformatic tools and classical biochemical dissection. Long-range comparative sequence analysis identified conserved noncoding sequences (CNSs). Clusters of epidermal-specific DNaseI-hypersensitive sites (HSs) mapped to specific CNSs. Increased prevalence of these HSs in barrier-deficient epidermis provides in vivo evidence of the regulation of the Sprr locus by these conserved sequences. Individual components of these HSs were cloned, and one was shown to have strong enhancer activity specific to conditions when the Sprr genes are coordinately upregulated.

Towards quantitative biology: integration of biological information to elucidate disease pathways and to guide drug discovery

Developing a new drug is a tedious and expensive undertaking. The recently developed high-throughput experimental technologies, summarised by the terms genomics, transcriptomics, proteomics and metabolomics provide for the first time ever the means to comprehensively monitor the molecular level of disease processes. The "-omics" technologies facilitate the systematic characterisation of a drug target's physiology, thereby helping to reduce the typically high attrition rates in discovery projects, and improving the overall efficiency of pharmaceutical research processes. Currently, the bottleneck for taking full advantage of the new experimental technologies are the rapidly growing volumes of automatically produced biological data. A lack of scalable database systems and computational tools for target discovery has been recognised as a major hurdle. In this review, an overview will be given on recent progress in computational biology that has an impact on drug discovery applications. The focus will be on novel in silico methods to reconstruct regulatory networks, signalling cascades, and metabolic pathways, with an emphasis on comparative genomics and microarray-based approaches. Promising methods, such as the mathematical simulation of pathway dynamics are discussed in the context of applications in discovery projects. The review concludes by exemplifying concrete data-driven studies in pharmaceutical research that demonstrate the value of integrated computational systems for drug target identification and validation, screening assay development, as well as drug candidate efficacy and toxicity evaluations.

The 5′ zebrafishscl promoter targets transcription to the brain, spinal cord, and hematopoietic and endothelial progenitors

The stem cell leukemia (SCL) gene encodes a basic helix-loop-helix transcription factor and is essential for embryonic angiogenesis, hematopoietic stem cell specification, and erythrocyte maturation. Here, we report the isolation and characterization of the zebrafish scl promoter. We show that a 5-kilobase (kb) genomic fragment immediately upstream of the translation start site is capable of targeting the enhanced green fluorescence protein (EGFP) expression to the anterior and posterior lateral mesoderm where the endogenous scl normally expresses. Detailed analysis of the stable transgenic fish reveals that this 5-kb upstream sequence is sufficient to direct the EGFP transcription to the brain, spinal cord, and hematopoietic-endothelial progenitors, possibly the hemangioblast, but not primitive erythrocyte, suggesting that the zebrafish scl transcription in hematopoietic-endothelial progenitors and erythrocyte is regulated by distinct cis element(s). Our study has defined the cis regulatory element(s) for zebrafish scl expression in the brain, spinal cord, and hematopoietic-endothelial progenitors and established a valuable transgenic line Tg(5'5kbscl:EGFP) for studying hematopoietic lineage development.

ABC: software for interactive browsing of genomic multiple sequence alignment data

Background

Alignment and comparison of related genome sequences is a powerful method to identify regions likely to contain functional elements. Such analyses are data intensive, requiring the inclusion of genomic multiple sequence alignments, sequence annotations, and scores describing regional attributes of columns in the alignment. Visualization and browsing of results can be difficult, and there are currently limited software options for performing this task.

Results

The Application for Browsing Constraints (ABC) is interactive Java software for intuitive and efficient exploration of multiple sequence alignments and data typically associated with alignments. It is used to move quickly from a summary view of the entire alignment via arbitrary levels of resolution to individual alignment columns. It allows for the simultaneous display of quantitative data, (e.g., sequence similarity or evolutionary rates) and annotation data (e.g. the locations of genes, repeats, and constrained elements). It can be used to facilitate basic comparative sequence tasks, such as export of data in plain-text formats, visualization of phylogenetic trees, and generation of alignment summary graphics.

Conclusions

The ABC is a lightweight, stand-alone, and flexible graphical user interface for browsing genomic multiple sequence alignments of specific loci, up to hundreds of kilobases or a few megabases in length. It is coded in Java for cross-platform use and the program and source code are freely available under the General Public License. Documentation and a sample data set are also available .

Noncoding regulatory sequences of Ciona exhibit strong correspondence between evolutionary constraint and functional importance

We show that sequence comparisons at different levels of resolution can efficiently guide functional analyses of regulatory regions in the ascidians Ciona savignyi and Ciona intestinalis. Sequence alignments of several tissue-specific genes guided discovery of minimal regulatory regions that are active in whole-embryo reporter assays. Using the Troponin I (TnI) locus as a case study, we show that more refined local sequence analyses can then be used to reveal functional substructure within a regulatory region. A high-resolution saturation mutagenesis in conjunction with comparative sequence analyses defined essential sequence elements within the TnI regulatory region. Finally, we found a significant, quantitative relationship between function and sequence divergence of noncoding functional elements. This work demonstrates the power of comparative sequence analysis between the two Ciona species for guiding gene regulatory experiments.

Striking differences between the mouse and the human alpha-fetoprotein enhancers

The alpha-fetoprotein (AFP) gene is expressed abundantly in the fetal liver and transcriptionally repressed in the adult liver, but can be reactivated during liver regeneration and in liver tumors. Previous studies identified three enhancers, E1, E2, and E3, upstream of the mouse and rat Afp genes and a single enhancer upstream of the human gene. We have compared the sequences upstream of the rodent and primate AFP genes. Our analysis demonstrates that the previously identified human enhancer is the counterpart to mouse E2. This comparison also reveals that a functional primate counterpart to the rodent E1 is absent due to a deletion that removes the core region of this enhancer. Furthermore, our studies identify a novel human enhancer corresponding to rodent E3. Despite the overall similarity of E3 between human and mouse, we found differences in transcription factor binding sites between these species. A C/EBP binding site is conserved but two other motifs in rodent E3, one that binds orphan nuclear receptors and a second that binds FoxA proteins, are not conserved in humans. The human counterpart to the rodent FoxA site can bind COUP-TF factors. Despite the overall sequence similarity in E3 between mice and humans, the difference in factor binding sites in E3, as well as the absence of E1 in primates, indicates that different mechanisms regulate AFP transcription in these different species.