CoMoDis: composite motif discovery in mammalian genomes

Evaluation of the NOD/SCID xenograft model for glucocorticoid-regulated gene expression in childhood B-cell precursor acute lymphoblastic leukemia

Background

Glucocorticoids such as prednisolone and dexamethasone are critical drugs used in multi-agent chemotherapy protocols used to treat acute lymphoblastic leukemia (ALL), and response to glucocorticoids is highly predictive of outcome. The NOD/SCID xenograft mouse model of ALL is a clinically relevant model in which the mice develop a systemic leukemia which retains the fundamental biological characteristics of the original disease. Here we report a study evaluating the NOD/SCID xenograft mouse model to investigate glucocorticoid-induced gene expression. Cells from a glucocorticoid-sensitive xenograft derived from a child with B-cell precursor ALL were inoculated into NOD/SCID mice. When highly engrafted the mice were randomized into groups of 4 to receive dexamethasone 15 mg/kg by intraperitoneal injection or vehicle control. Leukemia cells were harvested from mice spleens at 0, 8, 24 or 48 hours thereafter, and gene expression analyzed on Illumina WG-6_V3 chips, comparing all groups to time 0 hours.

Results

The 8 hour dexamethasone-treated timepoint had the highest number of significantly differentially expressed genes, with fewer observed at the 24 and 48 hour timepoints, and with minimal changes seen across the time-matched controls. When compared to publicly available datasets of glucocorticoid-induced gene expression from an in vitro cell line study and from an in vivo study of patients with ALL, at the level of pathways, expression changes in the 8 hour xenograft samples showed a similar response to patients treated with glucocorticoids. Replicate analysis revealed that at the 8 hour timepoint, a dataset with high signal and differential expression, using data from 3 replicates instead of 4 resulted in excellent recovery scores of > 0.9. However at other timepoints with less signal very poor recovery scores were obtained with 3 replicates.

Conclusions

The NOD/SCID xenograft mouse model provides a reproducible experimental system in which to investigate clinically-relevant mechanisms of drug-induced gene regulation in ALL; the 8 hour timepoint provides the highest number of significantly differentially expressed genes; time-matched controls are redundant and excellent recovery scores can be obtained with 3 replicates.

A Runx1-Smad6 rheostat controls Runx1 activity during embryonic hematopoiesis

The oncogenic transcription factor Runx1 is required for the specification of definitive hematopoietic stem cells (HSC) in the developing embryo. The activity of this master regulator is tightly controlled during development. The transcription factors that upregulate the expression of Runx1 also upregulate the expression of Smad6, the inhibitory Smad, which controls Runx1 activity by targeting it to the proteasome. Here we show that Runx1, in conjunction with Fli1, Gata2, and Scl, directly regulates the expression of Smad6 in the aorta-gonad-mesonephros (AGM) region in the developing embryo, where HSCs originate. Runx1 regulates Smad6 activity via a novel upstream enhancer, and Runx1 null embryos show reduced Smad6 transcripts in the yolk-sac and c-Kit-positive fetal liver cells. By directly regulating the expression of Smad6, Runx1 sets up a functional rheostat to control its own activity. The perturbation of this rheostat, using a proteasomal inhibitor, results in an increase in Runx1 and Smad6 levels that can be directly attributed to increased Runx1 binding to tissue-specific regulatory elements of these genes. Taken together, we describe a scenario in which a key hematopoietic transcription factor controls its own expression levels by transcriptionally controlling its controller.

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.

Transcription factor site dependencies in human, mouse and rat genomes

BACKGROUND

It is known that transcription factors frequently act together to regulate gene expression in eukaryotes. In this paper we describe a computational analysis of transcription factor site dependencies in human, mouse and rat genomes.

RESULTS

Our approach for quantifying tendencies of transcription factor binding sites to co-occur is based on a binding site scoring function which incorporates dependencies between positions, the use of information about the structural class of each transcription factor (major/minor groove binder), and also considered the possible implications of varying GC content of the sequences. Significant tendencies (dependencies) have been detected by non-parametric statistical methodology (permutation tests). Evaluation of obtained results has been performed in several ways: reports from literature (many of the significant dependencies between transcription factors have previously been confirmed experimentally); dependencies between transcription factors are not biased due to similarities in their DNA-binding sites; the number of dependent transcription factors that belong to the same functional and structural class is significantly higher than would be expected by chance; supporting evidence from GO clustering of targeting genes. Based on dependencies between two transcription factor binding sites (second-order dependencies), it is possible to construct higher-order dependencies (networks). Moreover results about transcription factor binding sites dependencies can be used for prediction of groups of dependent transcription factors on a given promoter sequence. Our results, as well as a scanning tool for predicting groups of dependent transcription factors binding sites are available on the Internet.

CONCLUSION

We show that the computational analysis of transcription factor site dependencies is a valuable complement to experimental approaches for discovering transcription regulatory interactions and networks. Scanning promoter sequences with dependent groups of transcription factor binding sites improve the quality of transcription factor predictions.

Integrating sequence, evolution and functional genomics in regulatory genomics

With genome analysis expanding from the study of genes to the study of gene regulation, 'regulatory genomics' utilizes sequence information, evolution and functional genomics measurements to unravel how regulatory information is encoded in the genome.

Chalcone synthase gene lineage diversification confirms allopolyploid evolutionary relationships of European rostrate violets

Phylogenetic relationships among and within the subsections of the genus Viola are still far from resolved. We present the first organismal phylogeny of predominantly western European species of subsection Rostratae based on the plastid trnS-trnG intron and intergenic spacer and the nuclear low-copy gene chalcone synthase (CHS) sequences. CHS is a key enzyme in the synthesis of flavonoids, which are important for flower pigmentation. Genes encoding for CHS are members of a multigene family. In Viola, 3 different CHS copies are present. CHS gene lineages obtained confirmed earlier hypotheses about reticulate relationships between species of Viola subsection Rostratae based on karyotype data. Comparison of the CHS gene lineage tree and the plastid species phylogeny of Viola reconstructed in this study indicates that the different CHS copies present in Viola are the products of both recent and more ancient duplications.

MOPAT: a graph-based method to predict recurrent cis-regulatory modules from known motifs

The identification of cis-regulatory modules (CRMs) can greatly advance our understanding of eukaryotic regulatory mechanism. Current methods to predict CRMs from known motifs either depend on multiple alignments or can only deal with a small number of known motifs provided by users. These methods are problematic when binding sites are not well aligned in multiple alignments or when the number of input known motifs is large. We thus developed a new CRM identification method MOPAT (motif pair tree), which identifies CRMs through the identification of motif modules, groups of motifs co-ccurring in multiple CRMs. It can identify ‘orthologous’ CRMs without multiple alignments. It can also find CRMs given a large number of known motifs. We have applied this method to mouse developmental genes, and have evaluated the predicted CRMs and motif modules by microarray expression data and known interacting motif pairs. We show that the expression profiles of the genes containing CRMs of the same motif module correlate significantly better than those of a random set of genes do. We also show that the known interacting motif pairs are significantly included in our predictions. Compared with several current methods, our method shows better performance in identifying meaningful CRMs.

DiRE: identifying distant regulatory elements of co-expressed genes

Regulation of gene expression in eukaryotic genomes is established through a complex cooperative activity of proximal promoters and distant regulatory elements (REs) such as enhancers, repressors and silencers. We have developed a web server named DiRE, based on the Enhancer Identification (EI) method, for predicting distant regulatory elements in higher eukaryotic genomes, namely for determining their chromosomal location and functional characteristics. The server uses gene co-expression data, comparative genomics and profiles of transcription factor binding sites (TFBSs) to determine TFBS-association signatures that can be used for discriminating specific regulatory functions. DiRE's unique feature is its ability to detect REs outside of proximal promoter regions, as it takes advantage of the full gene locus to conduct the search. DiRE can predict common REs for any set of input genes for which the user has prior knowledge of co-expression, co-function or other biologically meaningful grouping. The server predicts function-specific REs consisting of clusters of specifically-associated TFBSs and it also scores the association of individual transcription factors (TFs) with the biological function shared by the group of input genes. Its integration with the Array2BIO server allows users to start their analysis with raw microarray expression data. The DiRE web server is freely available at http://dire.dcode.org.

Runx genes are direct targets of Scl/Tal1 in the yolk sac and fetal liver

Transcription factors such as Scl/Tal1, Lmo2, and Runx1 are essential for the development of hematopoietic stem cells (HSCs). However, the precise mechanisms by which these factors interact to form transcriptional networks, as well as the identity of the genes downstream of these regulatory cascades, remain largely unknown. To this end, we generated an Scl(-/-) yolk sac cell line to identify candidate Scl target genes by global expression profiling after reintroduction of a TAT-Scl fusion protein. Bioinformatics analysis resulted in the identification of 9 candidate Scl target transcription factor genes, including Runx1 and Runx3. Chromatin immunoprecipitation confirmed that both Runx genes are direct targets of Scl in the fetal liver and that Runx1 is also occupied by Scl in the yolk sac. Furthermore, binding of an Scl-Lmo2-Gata2 complex was demonstrated to occur on the regions flanking the conserved E-boxes of the Runx1 loci and was shown to transactivate the Runx1 element. Together, our data provide a key component of the transcriptional network of early hematopoiesis by identifying downstream targets of Scl that can explain key aspects of the early Scl(-/-) phenotype.