Ian J. Donaldson

ORegAnno: an open-access community-driven resource for regulatory annotation

ORegAnno is an open-source, open-access database and literature curation system for community-based annotation of experimentally identified DNA regulatory regions, transcription factor binding sites and regulatory variants. The current release comprises 30 145 records curated from 922 publications and describing regulatory sequences for over 3853 genes and 465 transcription factors from 19 species. A new feature called the ‘publication queue’ allows users to input relevant papers from scientific literature as targets for annotation. The queue contains 4438 gene regulation papers entered by experts and another 54 351 identified by text-mining methods. Users can enter or ‘check out’ papers from the queue for manual curation using a series of user-friendly annotation pages. A typical record entry consists of species, sequence type, sequence, target gene, binding factor, experimental outcome and one or more lines of experimental evidence. An evidence ontology was developed to describe and categorize these experiments. Records are cross-referenced to Ensembl or Entrez gene identifiers, PubMed and dbSNP and can be visualized in the Ensembl or UCSC genome browsers. All data are freely available through search pages, XML data dumps or web services at: http://www.oreganno.org.

The transcriptional program controlled by the stem cell leukemia gene Scl/Tal1 during early embryonic hematopoietic development

The basic helix-loop-helix transcription factor Scl/Tal1 controls the development and subsequent differentiation of hematopoietic stem cells (HSCs). However, because few Scl target genes have been validated to date, the underlying mechanisms have remained largely unknown. In this study, we have used ChIP-Seq technology (coupling chromatin immunoprecipitation with deep sequencing) to generate a genome-wide catalog of Scl-binding events in a stem/progenitor cell line, followed by validation using primary fetal liver cells and comprehensive transgenic mouse assays. Transgenic analysis provided in vivo validation of multiple new direct Scl target genes and allowed us to reconstruct an in vivo validated network consisting of 17 factors and their respective regulatory elements. By coupling ChIP-Seq in model cell lines with in vivo transgenic validation and sophisticated bioinformatic analysis, we have identified a widely applicable strategy for the reconstruction of stem cell regulatory networks in which biologic material is otherwise limiting. Moreover, in addition to revealing multiple previously unrecognized links to known HSC regulators, as well as novel links to genes not previously implicated in HSC function, comprehensive transgenic analysis of regulatory elements provided substantial new insights into the transcriptional control of several important hematopoietic regulators, including Cbfa2t3h/Eto2, Cebpe, Nfe2, Zfpm1/Fog1, Erg, Mafk, Gfi1b, and Myb.

Identification of the REST regulon reveals extensive transposable element-mediated binding site duplication

The genome-wide mapping of gene-regulatory motifs remains a major goal that will facilitate the modelling of gene-regulatory networks and their evolution. The repressor element 1 is a long, conserved transcription factor-binding site which recruits the transcriptional repressor REST to numerous neuron-specific target genes. REST plays important roles in multiple biological processes and disease states. To map RE1 sites and target genes, we created a position specific scoring matrix representing the RE1 and used it to search the human and mouse genomes. We identified 1301 and 997 RE1s inhuman and mouse genomes, respectively, of which >40% are novel. By employing an ontological analysis we show that REST target genes are significantly enriched in a number of functional classes. Taking the novel REST target gene CACNA1A as an experimental model, we show that it can be regulated by multiple RE1s of different binding affinities, which are only partially conserved between human and mouse. A novel BLAST methodology indicated that many RE1s belong to closely related families. Most of these sequences are associated with transposable elements, leading us to propose that transposon-mediated duplication and insertion of RE1s has led to the acquisition of novel target genes by REST during evolution.

The SCL transcriptional network and BMP signaling pathway interact to regulate RUNX1 activity

Hematopoietic stem cell (HSC) development is regulated by several signaling pathways and a number of key transcription factors, which include Scl/Tal1, Runx1, and members of the Smad family. However, it remains unclear how these various determinants interact. Using a genome-wide computational screen based on the well characterized Scl +19 HSC enhancer, we have identified a related Smad6 enhancer that also targets expression to blood and endothelial cells in transgenic mice. Smad6, Bmp4, and Runx1 transcripts are concentrated along the ventral aspect of the E10.5 dorsal aorta in the aorta–gonad–mesonephros region from which HSCs originate. Moreover, Smad6, an inhibitor of Bmp4 signaling, binds and inhibits Runx1 activity, whereas Smad1, a positive mediator of Bmp4 signaling, transactivates the Runx1 promoter. Taken together, our results integrate three key determinants of HSC development; the Scl transcriptional network, Runx1 activity, and the Bmp4/Smad signaling pathway.

Association and linkage of leprosy phenotypes with HLA class II and tumour necrosis factor genes.

Previous analyses indicate major gene control of susceptibility to leprosy per se and the HLA class II region has been implicated in determining susceptibility and control of clinical phenotype. Segregation analysis using data from 76 Brazilian leprosy multi-case pedigrees (1166 individuals) supported a two locus model as the best fit: a recessive major gene and a recessive modifier gene(s) (single locus vs two locus model, P = 0.0007). Combined segregation and linkage analysis to the major locus, showed strong linkage to HLA class II (HLA-DQB1 P = 0.000002, HLA-DQA1 P = 0.000002, HLA-DRB1 P = 0.0000003) and tumour necrosis factor genes (TNF P = 0.00002, LTA P = 0.003). Extended transmission disequilibrium testing, using multiple affected family members, demonstrated that the common allele TNF*1 of the −308 promoter region polymorphism showed linkage and/or association with disease per se, at a high level of significance (P < 0.0001). Two locus transmission disequilibrium testing suggested susceptibility (TNF*1/LTA*2) and protective (TNF*2/LTA*2) haplotypes in the class iii region. Taken together the segregation and HLA analyses suggest the possibility of more than one susceptibility locus in the MHC.

Elucidation of the ELK1 target gene network reveals a role in the coordinate regulation of core components of the gene regulation machinery

Transcription factors play an important role in orchestrating the activation of specific networks of genes through targeting their proximal promoter and distal enhancer regions. However, it is unclear how the specificity of downstream responses is maintained by individual members of transcription-factor families and, in most cases, what their target repertoire is. We have used ChIP-chip analysis to identify the target genes of the ETS-domain transcription factor ELK1. Two distinct modes of ELK1 target gene selection are identified; the first involves redundant promoter binding with other ETS-domain family members; the second occurs through combinatorial binding with a second transcription factor SRF, which specifies a unique group of target genes. One of the most prominent groups of genes forming the ELK1 target network includes classes involved in core gene expression control, namely, components of the basal transcriptional machinery, the spliceosome and the ribosome. Amongst the set of genes encoding the basal transcription machinery components, are a functionally linked subset of GTFs and TAFs. Our study, therefore, reveals an unsuspected level of coordinate regulation of components of the core gene expression control machinery and also identifies two different modes of promoter targeting through binding with a second transcription factor or redundant binding with other ETS-domain family members.

Expression of the leukemia oncogene Lmo2 is controlled by an array of tissue-specific elements dispersed over 100 kb and bound by Tal1/Lmo2, Ets, and Gata factors.

The Lmo2 gene encodes a transcriptional cofactor critical for the development of hematopoietic stem cells. Ectopic LMO2 expression causes leukemia in T-cell acute lymphoblastic leukemia (T-ALL) patients and severe combined immunodeficiency patients undergoing retroviral gene therapy. Tightly controlled Lmo2 expression is therefore essential, yet no comprehensive analysis of Lmo2 regulation has been published so far. By comparative genomics, we identified 17 highly conserved noncoding elements, 9 of which revealed specific acetylation marks in chromatin-immunoprecipitation and microarray (ChIP-chip) assays performed across 250 kb of the Lmo2 locus in 11 cell types covering different stages of hematopoietic differentiation. All candidate regulatory regions were tested in transgenic mice. An extended LMO2 proximal promoter fragment displayed strong endothelial activity, while the distal promoter showed weak forebrain activity. Eight of the 15 distal candidate elements functioned as enhancers, which together recapitulated the full expression pattern of Lmo2, directing expression to endothelium, hematopoietic cells, tail, and forebrain. Interestingly, distinct combinations of specific distal regulatory elements were required to extend endothelial activity of the LMO2 promoter to yolk sac or fetal liver hematopoietic cells. Finally, Sfpi1/Pu.1, Fli1, Gata2, Tal1/Scl, and Lmo2 were shown to bind to and transactivate Lmo2 hematopoietic enhancers, thus identifying key upstream regulators and positioning Lmo2 within hematopoietic regulatory networks.

ID1 promotes expansion and survival of primary erythroid cells and is a target of JAK2V617F-STAT5 signaling

The discovery of JAK2V617F as an acquired mutation in the majority of patients with myeloproliferative disorders (MPDs) and the key role of the JAK2-STAT5 signaling cascade in normal hematopoiesis has focused attention on the downstream transcriptional targets of STAT5. Despite evidence of its vital role in normal erythropoiesis and its ability to recapitulate many of the features of myeloid malignancies, including the MPDs, few functionally validated targets of STAT5 have been described. Here we used a combination of comparative genomics and chromatin immunoprecipitation assays to identify ID1 as a novel target of the JAK2-STAT5 signaling axis in erythroid cells. STAT5 binds and transactivates a downstream enhancer of ID1, and ID1 expression levels correlate with the JAK2V617F mutation in both retrovirally transfected fetal liver cells and polycythemia vera patients. Knockdown and overexpression studies in a well-characterized erythroid differentiation assay from primary murine fetal liver cells demonstrated a survival-promoting action of ID1. This hitherto unrecognized function implicates ID1 in the expansion of erythroblasts during terminal differentiation and suggests that ID1 plays an important role in the pathogenesis of polycythemia vera. Furthermore, our findings contribute to an increasing body of evidence implicating ID proteins in a wider range of cellular functions than initially appreciated.

Genome-wide occupancy links Hoxa2 to Wnt–β-catenin signaling in mouse embryonic development

The regulation of gene expression is central to developmental programs and largely depends on the binding of sequence-specific transcription factors with cis-regulatory elements in the genome. Hox transcription factors specify the spatial coordinates of the body axis in all animals with bilateral symmetry, but a detailed knowledge of their molecular function in instructing cell fates is lacking. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) to identify Hoxa2 genomic locations in a time and space when it is actively instructing embryonic development in mouse. Our data reveals that Hoxa2 has large genome coverage and potentially regulates thousands of genes. Sequence analysis of Hoxa2-bound regions identifies high occurrence of two main classes of motifs, corresponding to Hox and Pbx–Hox recognition sequences. Examination of the binding targets of Hoxa2 faithfully captures the processes regulated by Hoxa2 during embryonic development; in addition, it uncovers a large cluster of potential targets involved in the Wnt-signaling pathway. In vivo examination of canonical Wnt–β-catenin signaling reveals activity specifically in Hoxa2 domain of expression, and this is undetectable in Hoxa2 mutant embryos. The comprehensive mapping of Hoxa2-binding sites provides a framework to study Hox regulatory networks in vertebrate developmental processes.

Sequencing and Characterisation of Rearrangements in Three S. pastorianus Strains Reveals the Presence of Chimeric Genes and Gives Evidence of Breakpoint Reuse

Gross chromosomal rearrangements have the potential to be evolutionarily advantageous to an adapting organism. The generation of a hybrid species increases opportunity for recombination by bringing together two homologous genomes. We sought to define the location of genomic rearrangements in three strains of Saccharomyces pastorianus, a natural lager-brewing yeast hybrid of Saccharomyces cerevisiae and Saccharomyces eubayanus, using whole genome shotgun sequencing. Each strain of S. pastorianus has lost species-specific portions of its genome and has undergone extensive recombination, producing chimeric chromosomes. We predicted 30 breakpoints that we confirmed at the single nucleotide level by designing species-specific primers that flank each breakpoint, and then sequencing the PCR product. These rearrangements are the result of recombination between areas of homology between the two subgenomes, rather than repetitive elements such as transposons or tRNAs. Interestingly, 28/30 S. cerevisiae- S. eubayanus recombination breakpoints are located within genic regions, generating chimeric genes. Furthermore we show evidence for the reuse of two breakpoints, located in HSP82 and KEM1, in strains of proposed independent origin.