Transcription factor expression landscape in Drosophila embryonic cell lines

BACKGROUND

Transcription factor (TF) proteins are a key component of the gene regulatory networks that control cellular fates and function. TFs bind DNA regulatory elements in a sequence-specific manner and modulate target gene expression through combinatorial interactions with each other, cofactors, and chromatin-modifying proteins. Large-scale studies over the last two decades have helped shed light on the complex network of TFs that regulate development in Drosophila melanogaster.

RESULTS

Here, we present a detailed characterization of expression of all known and predicted Drosophila TFs in two well-established embryonic cell lines, Kc167 and S2 cells. Using deep coverage RNA sequencing approaches we investigate the transcriptional profile of all 707 TF coding genes in both cell types. Only 103 TFs have no detectable expression in either cell line and 493 TFs have a read count of 5 or greater in at least one of the cell lines. The 493 TFs belong to 54 different DNA-binding domain families, with significant enrichment of those in the zf-C2H2 family. We identified 123 differentially expressed genes, with 57 expressed at significantly higher levels in Kc167 cells than S2 cells, and 66 expressed at significantly lower levels in Kc167 cells than S2 cells. Network mapping reveals that many of these TFs are crucial components of regulatory networks involved in cell proliferation, cell-cell signaling pathways, and eye development.

CONCLUSIONS

We produced a reference TF coding gene expression dataset in the extensively studied Drosophila Kc167 and S2 embryonic cell lines, and gained insight into the TF regulatory networks that control the activity of these cells.

Investigating the sequence landscape in the Drosophila initiator core promoter element using an enhanced MARZ algorithm

The core promoter elements are important DNA sequences for the regulation of RNA polymerase II transcription in eukaryotic cells. Despite the broad evolutionary conservation of these elements, there is extensive variation in the nucleotide composition of the actual sequences. In this study, we aim to improve our understanding of the complexity of this sequence variation in the TATA box and initiator core promoter elements in Drosophila melanogaster. Using computational approaches, including an enhanced version of our previously developed MARZ algorithm that utilizes gapped nucleotide matrices, several sequence landscape features are uncovered, including an interdependency between the nucleotides in position 2 and 5 in the initiator. Incorporating this information in an expanded MARZ algorithm improves predictive performance for the identification of the initiator element. Overall our results demonstrate the need to carefully consider detailed sequence composition features in core promoter elements in order to make more robust and accurate bioinformatic predictions.

The Dictyostelium discoideum genome lacks significant DNA methylation and uncovers palindromic sequences as a source of false positives in bisulfite sequencing

DNA methylation, the addition of a methyl (CH3) group to a cytosine residue, is an evolutionarily conserved epigenetic mark involved in a number of different biological functions in eukaryotes, including transcriptional regulation, chromatin structural organization, cellular differentiation and development. In the social amoeba Dictyostelium, previous studies have shown the existence of a DNA methyltransferase (DNMA) belonging to the DNMT2 family, but the extent and function of 5-methylcytosine in the genome are unclear. Here, we present the whole genome DNA methylation profile of Dictyostelium discoideum using deep coverage replicate sequencing of bisulfite-converted gDNA extracted from post-starvation cells. We find an overall very low number of sites with any detectable level of DNA methylation, occurring at significant levels in only 303-3432 cytosines out of the ∼7.5 million total cytosines in the genome depending on the replicate. Furthermore, a knockout of the DNMA enzyme leads to no overall decrease in DNA methylation. Of the identified sites, significant methylation is only detected at 11 sites in all four of the methylomes analyzed. Targeted bisulfite PCR sequencing and computational analysis demonstrate that the methylation profile does not change during development and that these 11 cytosines are most likely false positives generated by protection from bisulfite conversion due to their location in hairpin-forming palindromic DNA sequences. Our data therefore provide evidence that there is no significant DNA methylation in Dictyostelium before fruiting body formation and identify a reproducible experimental artifact from bisulfite sequencing.

Transcriptome profile in Drosophila Kc and S2 embryonic cell lines

Drosophila melanogaster cell lines are an important resource for a range of studies spanning genomics, molecular genetics, and cell biology. Amongst these valuable lines are Kc167 (Kc) and Schneider 2 (S2) cells, which were originally isolated in the late 1960s from embryonic sources and have been used extensively to investigate a broad spectrum of biological activities including cell-cell signaling and immune system function. Whole-genome tiling microarray analysis of total RNA from these two cell types was performed as part of the modENCODE project over a decade ago and revealed that they share a number of gene expression features. Here, we expand on these earlier studies by using deep-coverage RNA-sequencing approaches to investigate the transcriptional profile in Kc and S2 cells in detail. Comparison of the transcriptomes reveals that ∼75% of the 13,919 annotated genes are expressed at a detectable level in at least one of the cell lines, with the majority of these genes expressed at high levels in both cell lines. Despite the overall similarity of the transcriptional landscape in the two cell types, 2,588 differentially expressed genes are identified. Many of the genes with the largest fold change are known only by their "CG" designations, indicating that the molecular control of Kc and S2 cell identity may be regulated in part by a cohort of relatively uncharacterized genes. Our data also indicate that both cell lines have distinct hemocyte-like identities, but share active signaling pathways and express a number of genes in the network responsible for dorsal-ventral patterning of the early embryo.

Fitting thermodynamic-based models: Incorporating parameter sensitivity improves the performance of an evolutionary algorithm

A detailed comprehension of transcriptional regulation is critical to understanding the genetic control of development and disease across many different organisms. To more fully investigate the complex molecular interactions controlling the precise expression of genes, many groups have constructed mathematical models to complement their experimental approaches. A critical step in such studies is choosing the most appropriate parameter estimation algorithm to enable detailed analysis of the parameters that contribute to the models. In this study, we develop a novel set of evolutionary algorithms that use a pseudo-random Sobol Set to construct the initial population and incorporate parameter sensitivities into the adaptation of mutation rates, using local, global, and hybrid strategies. Comparison of the performance of these new algorithms to a number of current state-of-the-art global parameter estimation algorithms on a range of continuous test functions, as well as synthetic biological data representing models of gene regulatory systems, reveals improved performance of the new algorithms in terms of runtime, error and reproducibility. In addition, by analyzing the ability of these algorithms to fit datasets of varying quality, we provide the experimentalist with a guide to how the algorithms perform across a range of noisy data. These results demonstrate the improved performance of the new set of parameter estimation algorithms and facilitate meaningful integration of model parameters and predictions in our understanding of the molecular mechanisms of gene regulation.

Analyzing the stability of gene expression using a simple reaction-diffusion model in an early Drosophila embryo

In all complex organisms, the precise levels and timing of gene expression controls vital biological processes. In higher eukaryotes, including the fruit fly Drosophila melanogaster, the complex molecular control of transcription (the synthesis of RNA from DNA) and translation (the synthesis of proteins from RNA) events driving this gene expression are not fully understood. In particular, for Drosophila melanogaster, there is a plethora of experimental data, including quantitative measurements of both RNA and protein concentrations, but the precise mechanisms that control the dynamics of gene expression during early development and the processes which lead to steady-state levels of certain proteins remain elusive. This study analyzes a current mathematical modeling approach in an attempt to better understand the long-term behavior of gene regulation. The model is a modified reaction-diffusion equation which has been previously employed in predicting gene expression levels and studying the relative contributions of transcription and translation events to protein abundance [10,11,24]. Here, we use Matrix Algebra and Analysis techniques to study the stability of the gene expression system and analyze equilibria, using very general assumptions regarding the parameter values incorporated into the model. We prove that, given realistic biological parameter values, the system will result in a unique, stable equilibrium solution. Additionally, we give an example of this long-term behavior using the model alongside actual experimental data obtained from Drosophila embryos.

Rapid and efficient purification of Drosophila homeodomain transcription factors for biophysical characterization

Homeodomain transcription factors (HD TFs) are a large class of evolutionarily conserved DNA binding proteins that contain a basic 60-amino acid region required for binding to specific DNA sites. In Drosophila melanogaster, many of these HD TFs are expressed in the early embryo and control transcription of target genes in development through their interaction with cis-regulatory modules. Previous studies where some of the Drosophila HD TFs were purified required the use of strong denaturants (i.e. 6 M urea) and multiple chromatography columns, making the downstream biochemical examination of the isolated protein difficult. To circumvent these obstacles, we have developed a streamlined expression and purification protocol to produce large yields of Drosophila HD TFs. Using the HD TFs FUSHI-TARAZU (FTZ), ANTENNAPEDIA (ANTP), ABDOMINAL-A (ABD-A), ABDOMINAL-B (ABD-B), and ULTRABITHORAX (UBX) as examples, we demonstrate that our 3-day protocol involving the overexpression of His₆-SUMO fusion constructs in E. coli followed by a Ni²⁺-IMAC, SUMO-tag cleavage with the SUMO protease Ulp1, and a heparin column purification produces pure, soluble protein in biological buffers around pH 7 in the absence of denaturants. Electrophoretic mobility shift assays (EMSA) confirm that the purified HD proteins are functional and nuclear magnetic resonance (NMR) spectra confirm that the purified HDs are well-folded. These purified HD TFs can be used in future biophysical experiments to structurally and biochemically characterize how and why these HD TFs bind to different DNA sequences and further probe how nucleotide differences contribute to TF-DNA specificity in the HD family.

Nucleotide Interdependency in Transcription Factor Binding Sites in the Drosophila Genome

A long-standing objective in modern biology is to characterize the molecular components that drive the development of an organism. At the heart of eukaryotic development lies gene regulation. On the molecular level, much of the research in this field has focused on the binding of transcription factors (TFs) to regulatory regions in the genome known as cis-regulatory modules (CRMs). However, relatively little is known about the sequence-specific binding preferences of many TFs, especially with respect to the possible interdependencies between the nucleotides that make up binding sites. A particular limitation of many existing algorithms that aim to predict binding site sequences is that they do not allow for dependencies between nonadjacent nucleotides. In this study, we use a recently developed computational algorithm, MARZ, to compare binding site sequences using 32 distinct models in a systematic and unbiased approach to explore nucleotide dependencies within binding sites for 15 distinct TFs known to be critical to Drosophila development. Our results indicate that many of these proteins have varying levels of nucleotide interdependencies within their DNA recognition sequences, and that, in some cases, models that account for these dependencies greatly outperform traditional models that are used to predict binding sites. We also directly compare the ability of different models to identify the known KRUPPEL TF binding sites in CRMs and demonstrate that a more complex model that accounts for nucleotide interdependencies performs better when compared with simple models. This ability to identify TFs with critical nucleotide interdependencies in their binding sites will lead to a deeper understanding of how these molecular characteristics contribute to the architecture of CRMs and the precise regulation of transcription during organismal development.

Spatial distribution of predicted transcription factor binding sites in Drosophila ChIP peaks

In the development of the Drosophila embryo, gene expression is directed by the sequence-specific interactions of a large network of protein transcription factors (TFs) and DNA cis-regulatory binding sites. Once the identity of the typically 8-10bp binding sites for any given TF has been determined by one of several experimental procedures, the sequences can be represented in a position weight matrix (PWM) and used to predict the location of additional TF binding sites elsewhere in the genome. Often, alignments of large (>200bp) genomic fragments that have been experimentally determined to bind the TF of interest in Chromatin Immunoprecipitation (ChIP) studies are trimmed under the assumption that the majority of the binding sites are located near the center of all the aligned fragments. In this study, ChIP/chip datasets are analyzed using the corresponding PWMs for the well-studied TFs; CAUDAL, HUNCHBACK, KNIRPS and KRUPPEL, to determine the distribution of predicted binding sites. All four TFs are critical regulators of gene expression along the anterio-posterior axis in early Drosophila development. For all four TFs, the ChIP peaks contain multiple binding sites that are broadly distributed across the genomic region represented by the peak, regardless of the prediction stringency criteria used. This result suggests that ChIP peak trimming may exclude functional binding sites from subsequent analyses.

Parent-of-origin effects on genome-wide DNA methylation in the Cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation

BACKGROUND

Intersexual genomic conflict sometimes leads to unequal expression of paternal and maternal alleles in offspring, resulting in parent-of-origin effects. In honey bees reciprocal crosses can show strong parent-of-origin effects, supporting theoretical predictions that genomic imprinting occurs in this species. Mechanisms behind imprinting in honey bees are unclear but differential DNA methylation in eggs and sperm suggests that DNA methylation could be involved. Nonetheless, because DNA methylation is multifunctional, it is difficult to separate imprinting from other roles of methylation. Here we use a novel approach to investigate parent-of-origin DNA methylation in honey bees. In the subspecies Apis mellifera capensis, reproduction of females occurs either sexually by fertilization of eggs with sperm, or via thelytokous parthenogenesis, producing female embryos derived from two maternal genomes.

RESULTS

We compared genome-wide methylation patterns of sexually-produced, diploid embryos laid by a queen, with parthenogenetically-produced diploid embryos laid by her daughters. Thelytokous embryos inheriting two maternal genomes had fewer hypermethylated genes compared to fertilized embryos, supporting the prediction that fertilized embryos have increased methylation due to inheritance of a paternal genome. However, bisulfite PCR and sequencing of a differentially methylated gene, Stan (GB18207) showed strong allele-specific methylation that was maintained in both fertilized and thelytokous embryos. For this gene, methylation was associated with haplotype, not parent of origin.

CONCLUSIONS

The results of our study are consistent with predictions from the kin theory of genomic imprinting. However, our demonstration of allele-specific methylation based on sequence shows that genome-wide differential methylation studies can potentially confound imprinting and allele-specific methylation. It further suggests that methylation patterns are heritable or that specific sequence motifs are targets for methylation in some genes.

Global sensitivity analysis of a dynamic model for gene expression in Drosophila embryos

It is well known that gene regulation is a tightly controlled process in early organismal development. However, the roles of key processes involved in this regulation, such as transcription and translation, are less well understood, and mathematical modeling approaches in this field are still in their infancy. In recent studies, biologists have taken precise measurements of protein and mRNA abundance to determine the relative contributions of key factors involved in regulating protein levels in mammalian cells. We now approach this question from a mathematical modeling perspective. In this study, we use a simple dynamic mathematical model that incorporates terms representing transcription, translation, mRNA and protein decay, and diffusion in an early Drosophila embryo. We perform global sensitivity analyses on this model using various different initial conditions and spatial and temporal outputs. Our results indicate that transcription and translation are often the key parameters to determine protein abundance. This observation is in close agreement with the experimental results from mammalian cells for various initial conditions at particular time points, suggesting that a simple dynamic model can capture the qualitative behavior of a gene. Additionally, we find that parameter sensitivites are temporally dynamic, illustrating the importance of conducting a thorough global sensitivity analysis across multiple time points when analyzing mathematical models of gene regulation.

MARZ: an algorithm to combinatorially analyze gapped n-mer models of transcription factor binding

Background

A key challenge in understanding the molecular mechanisms that control gene regulation is the characterization of the specificity with which transcription factor proteins bind to specific DNA sequences. A number of computational approaches have been developed to examine these interactions, including simple mononucleotide and dinucleotide position weight matrix models.

Results

Here we develop a novel, unbiased computational algorithm, MARZ, that systematically analyzes all possible gapped matrices across a fixed number of nucleotides. In addition, to evaluate the ability of these matrix models to predict in vivo binding sites, we utilize a new scoring system and, in combination with established scoring methods and statistical analysis, test the performance of 32 different gapped matrices on the well characterized HUNCHBACK transcription factor in Drosophila.

Conclusions

Our results indicate that in many cases gapped matrix models can outperform traditional models, but that the relative strength of the binding sites considered in the analysis can profoundly influence the predictive ability of specific models.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0446-3) contains supplementary material, which is available to authorized users.

Deciphering the combinatorial architecture of a Drosophila homeotic gene enhancer

In Drosophila, the 330 kb bithorax complex regulates cellular differentiation along the anterior–posterior axis during development in the thorax and abdomen and is comprised of three homeotic genes: Ultrabithorax, abdominal-A, and Abdominal-B. The expression of each of these genes is in turn controlled through interactions between transcription factors and a number of cis-regulatory modules in the neighboring intergenic regions. In this study, we examine how the sequence architecture of transcription factor binding sites mediates the functional activity of one of these cis-regulatory modules. Using computational, mathematical modeling and experimental molecular genetic approaches we investigate the IAB7b enhancer, which regulates Abdominal-B expression specifically in the presumptive seventh and ninth abdominal segments of the early embryo. A cross-species comparison of the IAB7b enhancer reveals an evolutionarily conserved signature motif containing two FUSHI-TARAZU activator transcription factor binding sites. We find that the transcriptional repressors KNIRPS, KRUPPEL and GIANT are able to restrict reporter gene expression to the posterior abdominal segments, using different molecular mechanisms including short-range repression and competitive binding. Additionally, we show the functional importance of the spacing between the two FUSHI-TARAZU binding sites and discuss the potential importance of cooperativity for transcriptional activation. Our results demonstrate that the transcriptional output of the IAB7b cis-regulatory module relies on a complex set of combinatorial inputs mediated by specific transcription factor binding and that the sequence architecture at this enhancer is critical to maintain robust regulatory function.

The dynamic DNA methylation cycle from egg to sperm in the honey bee Apis mellifera

In honey bees (Apis mellifera), the epigenetic mark of DNA methylation is central to the developmental regulation of caste differentiation, but may also be involved in additional biological functions. In this study, we examine the whole genome methylation profiles of three stages of the haploid honey bee genome: unfertilised eggs, the adult drones that develop from these eggs and the sperm produced by these drones. These methylomes reveal distinct patterns of methylation. Eggs and sperm show 381 genes with significantly different CpG methylation patterns, with the vast majority being more methylated in eggs. Adult drones show greatly reduced levels of methylation across the genome when compared with both gamete samples. This suggests a dynamic cycle of methylation loss and gain through the development of the drone and during spermatogenesis. Although fluxes in methylation during embryogenesis may account for some of the differentially methylated sites, the distinct methylation patterns at some genes suggest parent-specific epigenetic marking in the gametes. Extensive germ line methylation of some genes possibly explains the lower-than-expected frequency of CpG sites in these genes. We discuss the potential developmental and evolutionary implications of methylation in eggs and sperm in this eusocial insect species.

A parent-of-origin effect on honeybee worker ovary size

Apis mellifera capensis is unique among honeybees in that unmated workers can produce pseudo-clonal female offspring via thelytokous parthenogenesis. Workers use this ability to compete among themselves and with their queen to be the mother of new queens. Males could therefore enhance their reproductive success by imprinting genes that enhance fertility in their daughter workers. This possibility sets the scene for intragenomic conflict between queens and drones over worker reproductive traits. Here, we show a strong parent-of-origin effect for ovary size (number of ovarioles) in reciprocal crosses between two honeybee subspecies, A. m. capensis and Apis mellifera scutellata. In this cross, workers with an A. m. capensis father had 30% more ovarioles than genotypically matched workers with an A. m. scutellata father. Other traits we measured (worker weight at emergence and the presence/absence of a spermatheca) are influenced more by rearing conditions than by parent-of-origin effects. Our study is the first to show a strong epigenetic (or, less likely, cytoplasmic maternal) effect for a reproductive trait in the honeybee and suggests that a search for parent-of-origin effects in other social insects may be fruitful.

Flanking sequence context-dependent transcription factor binding in early Drosophila development

BACKGROUND

Gene expression in the Drosophila embryo is controlled by functional interactions between a large network of protein transcription factors (TFs) and specific sequences in DNA cis-regulatory modules (CRMs). The binding site sequences for any TF can be experimentally determined and represented in a position weight matrix (PWM). PWMs can then be used to predict the location of TF binding sites in other regions of the genome, although there are limitations to this approach as currently implemented.

RESULTS

In this proof-of-principle study, we analyze 127 CRMs and focus on four TFs that control transcription of target genes along the anterio-posterior axis of the embryo early in development. For all four of these TFs, there is some degree of conserved flanking sequence that extends beyond the predicted binding regions. A potential role for these conserved flanking sequences may be to enhance the specificity of TF binding, as the abundance of these sequences is greatly diminished when we examine only predicted high-affinity binding sites.

CONCLUSIONS

Expanding PWMs to include sequence context-dependence will increase the information content in PWMs and facilitate a more efficient functional identification and dissection of CRMs.

Molecular dissection of cis-regulatory modules at the Drosophila bithorax complex reveals critical transcription factor signature motifs

At the Drosophila melanogaster bithorax complex (BX-C) over 330kb of intergenic DNA is responsible for directing the transcription of just three homeotic (Hox) genes during embryonic development. A number of distinct enhancer cis-regulatory modules (CRMs) are responsible for controlling the specific expression patterns of the Hox genes in the BX-C. While it has proven possible to identify orthologs of known BX-C CRMs in different Drosophila species using overall sequence conservation, this approach has not proven sufficiently effective for identifying novel CRMs or defining the key functional sequences within enhancer CRMs. Here we demonstrate that the specific spatial clustering of transcription factor (TF) binding sites is important for BX-C enhancer activity. A bioinformatic search for combinations of putative TF binding sites in the BX-C suggests that simple clustering of binding sites is frequently not indicative of enhancer activity. However, through molecular dissection and evolutionary comparison across the Drosophila genus we discovered that specific TF binding site clustering patterns are an important feature of three known BX-C enhancers. Sub-regions of the defined IAB5 and IAB7b enhancers were both found to contain an evolutionarily conserved signature motif of clustered TF binding sites which is critical for the functional activity of the enhancers. Together, these results indicate that the spatial organization of specific activator and repressor binding sites within BX-C enhancers is of greater importance than overall sequence conservation and is indicative of enhancer functional activity.

Transcriptional repression by the proximal exonic region at the human TERT gene

In humans, the enzyme telomerase (hTERT) is responsible for the synthesis of new repeat sequences at the telomeres of chromosomes. Although active in early embryogenesis, the hTERT gene is transcriptionally silenced in almost all somatic cells in the adult, but is aberrantly re-activated in over 90% of human cancers. The molecular mechanisms responsible for repression of this gene are thought to involve the transcription factor CTCF. In this study, we bioinformatically identify putative CTCF binding sites in the hTERT proximal exonic region (PER) and determine their functional relevance in mediating transcriptional silencing at this gene. Tests using a reporter gene assay in HeLa cancer cells demonstrate that a sub-region of the PER exhibits strong transcriptional repressive activity. This repression is independent of the previously identified CTCF binding site near the transcriptional start site of the hTERT gene. In addition, site directed mutagenesis of three predicted CTCF binding sites, including a previously characterized in vivo site in exon 2, does not result in a loss of the repression mediated by the PER. The results from this study indicate that expression of the hTERT gene in HeLa cells is regulated by sequences in the PER. This transcriptional control is mediated through additional regulatory molecular mechanisms, independent of CTCF binding.

Disruption of the abdominal-B promoter tethering element results in a loss of long-range enhancer-directed Hox gene expression in Drosophila

There are many examples within gene complexes of transcriptional enhancers interacting with only a subset of target promoters. A number of molecular mechanisms including promoter competition, insulators and chromatin looping are thought to play a role in regulating these interactions. At the Drosophila bithorax complex (BX-C), the IAB5 enhancer specifically drives gene expression only from the Abdominal-B (Abd-B) promoter, even though the enhancer and promoter are 55 kb apart and are separated by at least three insulators. In previous studies, we discovered that a 255 bp cis-regulatory module, the promoter tethering element (PTE), located 5′ of the Abd-B transcriptional start site is able to tether IAB5 to the Abd-B promoter in transgenic embryo assays. In this study we examine the functional role of the PTE at the endogenous BX-C using transposon-mediated mutagenesis. Disruption of the PTE by P element insertion results in a loss of enhancer-directed Abd-B expression during embryonic development and a homeotic transformation of abdominal segments. A partial deletion of the PTE and neighboring upstream genomic sequences by imprecise excision of the P element also results in a similar loss of Abd-B expression in embryos. These results demonstrate that the PTE is an essential component of the regulatory network at the BX-C and is required in vivo to mediate specific long-range enhancer-promoter interactions.

Dissecting the regulatory switches of development: lessons from enhancer evolution in Drosophila

Cis-regulatory modules are non-protein-coding regions of DNA essential for the control of gene expression. One class of regulatory modules is embryonic enhancers, which drive gene expression during development as a result of transcription factor protein binding at the enhancer sequences. Recent comparative studies have begun to investigate the evolution of the sequence architecture within enhancers. These analyses are illuminating the way that developmental biologists think about enhancers by revealing their molecular mechanism of function.

Functional evolution of cis-regulatory modules at a homeotic gene in Drosophila

It is a long-held belief in evolutionary biology that the rate of molecular evolution for a given DNA sequence is inversely related to the level of functional constraint. This belief holds true for the protein-coding homeotic (Hox) genes originally discovered in Drosophila melanogaster. Expression of the Hox genes in Drosophila embryos is essential for body patterning and is controlled by an extensive array of cis-regulatory modules (CRMs). How the regulatory modules functionally evolve in different species is not clear. A comparison of the CRMs for the Abdominal-B gene from different Drosophila species reveals relatively low levels of overall sequence conservation. However, embryonic enhancer CRMs from other Drosophila species direct transgenic reporter gene expression in the same spatial and temporal patterns during development as their D. melanogaster orthologs. Bioinformatic analysis reveals the presence of short conserved sequences within defined CRMs, representing gap and pair-rule transcription factor binding sites. One predicted binding site for the gap transcription factor KRUPPEL in the IAB5 CRM was found to be altered in Superabdominal (Sab) mutations. In Sab mutant flies, the third abdominal segment is transformed into a copy of the fifth abdominal segment. A model for KRUPPEL-mediated repression at this binding site is presented. These findings challenge our current understanding of the relationship between sequence evolution at the molecular level and functional activity of a CRM. While the overall sequence conservation at Drosophila CRMs is not distinctive from neighboring genomic regions, functionally critical transcription factor binding sites within embryonic enhancer CRMs are highly conserved. These results have implications for understanding mechanisms of gene expression during embryonic development, enhancer function, and the molecular evolution of eukaryotic regulatory modules.

Non-genic transcription at the Drosophila bithorax complex functional activity of the dark matter of the genome

Drosophila melanogaster is a powerful model system for the study of gene regulation due to its short generation time, high fertility and the availability of various genetic tools to manipulate the genome. Investigation into the regulation of homeotic genes and their role in embryonic patterning during development was pioneered in Drosophila. Recently, the molecular mechanisms responsible for regulating gene expression in the bithorax complex have been the focus of active study. Many of these studies have pointed to the importance of cis-regulatory modules, genetic sequences that direct the temporal and spatial patterns of gene expression over large genomic distances. Additional components of the regulatory code have emerged beyond the primary DNA sequence. In particular, non-genic transcription is an important mechanism for controlling gene expression either through direct transcriptional mechanisms that mediate dynamic epigenetic control of the chromatin environment or through functional activity of the RNA products.

Between transcription and translation: Re-defining RNA and regulation

The diverse functional roles for RNA molecules in cells of the developing embryo have been an area of intense study in the last few years. Progress reported at the 49(th) Annual Drosophila Research Conference in San Diego, California highlighted many of the varied mechanistic activities for RNAs. In particular, talks at the 'RNA Biology' platform session provided a great deal of insight into the function of RNA transcripts and their associated protein complexes. The topics covered included: (1) a large-scale screen examining the localization of mRNAs during embryonic development, (2) mechanisms of mRNA transport in different cell types, (3) localization-dependent repression of mRNA translation and (4) the activity of the RNAi machinery in insulator-mediated chromatin structures. Our journey through the modern RNA world clearly indicates that we should be considering a much more expansive role for RNAs in molecular biology.

A novel promoter-tethering element regulates enhancer-driven gene expression at the bithorax complex in the Drosophila embryo

A key question in our understanding of the cis-regulation of gene expression during embryonic development has been the molecular mechanism that directs enhancers to specific promoters within a gene complex. Promoter competition and insulators are thought to play a role in regulating these interactions. In the bithorax complex of Drosophila, the IAB5 enhancer is located 55 kb 3' of the Abdominal-B (Abd-B) promoter and 48 kb 5' of the abdominal-A (abd-A) promoter. Although roughly equidistant from the two promoters, IAB5 specifically interacts only with the Abdominal-B promoter, even though the enhancer and promoter are separated by at least two insulators. Here we demonstrate that a 255 bp element, located 40 bp 5' of the Abd-B transcriptional start site, has a novel cis-regulatory activity as it is able to tether IAB5 to the Abd-B promoter in transgenic embryos. The tethering element is sufficient to direct IAB5 to an ectopic promoter in competition assays. Deletion of the promoter-tethering element results in the redirection of enhancer-driven gene expression on transgenes. Taken together, these results provide evidence that specific long-range enhancer-promoter interactions in the bithorax complex are regulated by a tethering element 5' of the Abd-B promoter. We discuss a bioinformatic analysis of the tethering element across different Drosophila species and a possible molecular mechanism by which this element functions. We also examine existing evidence that this novel class of cis-regulatory elements might regulate enhancer-promoter specificity at other gene complexes.

The abdominal-B promoter tethering element mediates promoter-enhancer specificity at the Drosophila bithorax complex

At the Drosophila bithorax complex many distinct classes of cis-regulatory modules work collectively during development to control gene expression. Abdominal-B (Abd-B) is one of three homeotic genes in the BX-C and is expressed in specific presumptive abdominal segments in the embryo. The transcription of Abd-B is tightly controlled by an array of cis-regulatory modules that direct its expression over extended genomic distances. These regulatory modules include promoters, insulators, silencers, enhancers, promoter targeting sequences and the recently identified promoter tethering element (PTE). To activate gene expression at the endogenous complex, enhancers located >50 kb away must bypass intervening insulators to interact with the Abd-B promoter. The molecular mechanisms that allow enhancers to bypass insulators are not currently well understood. In this short article, we report on a novel mechanism for insulator bypass involving the PTE. In addition, we use bioinformatic analysis across twelve Drosophila genomes to identify putative cis-regulatory sequences that may be capable of facilitating specific promoter-enhancer interactions at the bithorax complex and propose a model for their molecular function during development.

Chromatin looping mediates boundary element promoter interactions

One facet of the control of gene expression is long-range promoter regulation by distant enhancers. It is an important component of the regulation of genes that control metazoan development and has been appreciated for some time but the molecular mechanisms underlying this regulation have remained poorly understood. A recent study by Cleard and colleagues1 reports the first in vivo evidence of chromatin looping and boundary element promoter interaction. Specifically, they studied the function of a boundary element within the cis-regulatory region of the Abdominal-B (Abd-B) gene of Drosophila melanogaster.

The human and mouse H19 imprinting control regions harbor an evolutionarily conserved silencer element that functions on transgenes in Drosophila

Differentially methylated regions have been characterized at a number of imprinted gene complexes with important roles in the regulation of monoallelic expression of one or more genes. The differentially methylated imprinting control region (ICR) located upstream of the murine H19 gene has been shown to control the imprinted expression of H19 and the coordinately regulated Igf2 gene by acting as a transcriptional silencer. In this study, we show that the murine ICR maintains this function when tested in an in vivo transgenic Drosophila assay in the absence of DNA methylation. Furthermore, the H19 ICR interacts distinctively with Drosophila promoters of different regulatory strengths. We also demonstrate that the comparable region upstream of the human H19 gene is a multipartite cis-regulatory element, demonstrating silencing function when tested in mammalian and Drosophila systems. These results indicate a conservation of the H19/Igf2 imprinting mechanism between humans and mice and further elucidate the functional activities of the H19 ICR. They demonstrate the value of Drosophila as an in vivo system for testing function and interaction of eukaryotic regulatory elements and that mechanisms of transcriptional cis-regulation in mammals and Drosophila are conserved.

Unraveling cis-regulatory mechanisms at the abdominal-A and Abdominal-B genes in the Drosophila bithorax complex

Genome sequencing has revealed that in metazoans, only a small percentage of DNA actually codes for functional proteins. Research efforts have focused on elucidating the purpose of the rest of the genome, which was initially largely thought of as mere 'junk' DNA. One genomic region that is proving to be a rich source of new information is the Drosophila bithorax complex (BX-C). At this homeotic gene complex, many different classes of cis-regulatory elements, such as insulators, silencers, enhancers, and promoters, work together to tightly control gene expression during development. Recent studies have begun to unravel the intricate nature of these regulatory interactions. The BX-C was first discovered and characterized by Ed Lewis over three decades ago. In his seminal 1978 Nature paper, Lewis speculated that "substances" originating from the nongenic regions of the BX-C may regulate expression of the neighboring abdominal-A and Abdominal-B homeotic genes. A number of discoveries in the last few years suggest that he was right. The activation of some of the cis-sequences at the complex appears to be controlled by nongenic transcription, providing a further level of regulatory complexity to regions of nonprotein coding DNA. The hope is that these studies of gene regulation at the BX-C in the humble fruit fly will provide clues as to how vast intergenic regions contribute to the incredible complexity of gene regulation in other species, including humans.

The 3' portion of the mouse H19 Imprinting-Control Region is required for proper tissue-specific expression of the Igf2 gene

Genomic imprinting at the H19/Igf2 locus is governed by a cis-acting Imprinting-Control Region (ICR), located 2 kb upstream of the H19 gene. This region possesses an insulator function which is activated on the unmethylated maternal allele through the binding of the CTCF factor. It has been previously reported that paternal transmission of the H19(SilK) deletion, which removes the 3' portion of H19 ICR, leads to the loss of H19 imprinting. Here we show that, in the liver, this reactivation of the paternal H19 gene is concomitant to a dramatic decrease in Igf2 mRNA levels. This deletion alters higher-order chromatin architecture, Igf2 promoter usage and tissue-specific expression. Therefore, when methylated, the 3' portion of the H19 ICR is a bi-functional regulatory element involved not only in H19 imprinting but also in 'formatting' the higher-order chromatin structure for proper tissue-specific expression of both H19 and Igf2 genes.

Transcription defines the embryonic domains of cis-regulatory activity at the Drosophila bithorax complex

The extensive infraabdominal (iab) region contains a number of cis-regulatory elements, including enhancers, silencers, and insulators responsible for directing the developmental expression of the abdominal-A and Abdominal-B homeotic genes at the Drosophila bithorax complex. It is unclear how these regulatory elements are primed for activity early in embryogenesis, but the 100-kb intergenic region is subject to a complex transcriptional program. Here, we use molecular and genetic methods to examine the functional activity of the RNAs produced from this region and their role in cis regulation. We show that a subset of these transcripts demonstrates a distinct pattern of cellular localization. Furthermore, the transcripts from each iab region are discrete and the transcripts do not spread across the insulator elements that delineate the iab regions. In embryos carrying a Mcp deletion, the intergenic transcription pattern is disrupted in the iab4 region and the fourth abdominal segment is transformed into the fifth. We propose that intergenic transcription is required early in embryogenesis to initiate the activation of the Drosophila bithorax complex and define the domains of activity for the iab cis-regulatory elements. We also discuss a possible mechanism by which this may occur.

Characterization of the intergenic RNA profile at abdominal-A and Abdominal-B in the Drosophila bithorax complex

The correct spatial expression of two Drosophila bithorax complex (BX-C) genes, abdominal-A (abdA) and Abdominal-B (AbdB), is dependent on the 100-kb intergenic infraabdominal (iab) region. The iab region is known to contain a number of different domains (iab2 through iab8) that harbor cis-regulatory elements responsible for directing expression of abdA and AbdB in the second through eighth abdominal segments. Here, we use in situ hybridization to perform high-resolution mapping of the transcriptional activity in the iab control regions. We show that transcription of the control regions themselves is abundant and precedes activation of the abdA and AbdB genes. As with the homeotic genes of the BX-C, the transcription patterns of the RNAs from the iab control regions demonstrate colinearity with the sequence of the iab regions along the chromosome and the domains in the embryo under the control of the specific iab regions. These observations suggest that the intergenic RNAs may play a role in initiating cis regulation at the BX-C early in development.

Methylation-dependent silencing at the H19 imprinting control region by MeCP2

Methylation of CpG dinucleotides is correlated with transcriptional repression of genes, including imprinted genes. In the case of the imprinted H19 gene, a 2 kb imprinting control region (ICR) is subject to differential methylation, as it is methylated only on the silenced paternal allele. This region has previously been shown to act as a silencer element at the endogenous locus. The proteins that bind at the H19 differentially methylated domain (DMD) and mediate transcriptional silencing have yet to be identified, although a family of proteins containing a methyl-CpG-binding domain (MBD), of which MeCP2 is the best characterised, are obvious candidates. MeCP2 can bind to a single methylated CpG dinucleotide through its MBD and also contains a transcriptional repression domain (TRD). The TRD interacts with Sin3a and histone deacetylases (HDACs) in vivo, forming a repressive complex. Here we show that MeCP2 is recruited to the H19 DMD in vivo and can silence a reporter gene regulated by the H19 DMD in a methylation-dependent manner. This repression can be alleviated by deletion of the TRD from MeCP2 or by inhibition of HDAC activity. These data indicate that transcriptional silencing from the H19 ICR involves recruitment of MeCP2 and presumably an associated protein complex with deacetylase activity. This complex may also be recruited to the ICR in vivo, resulting in a compact, repressive chromatin structure capable of silencing the paternal H19 allele.

Novel conserved elements upstream of theH19gene are transcribed and act as mesodermal enhancers

The reciprocally imprinted H19 and Igf2 genes form a co-ordinately regulated 130 kb unit in the mouse controlled by widely dispersed enhancers, epigenetically modified silencers and an imprinting control region (ICR). Comparative human and mouse genomic sequencing between H19 and Igf2 revealed two novel regions of strong homology upstream of the ICR termed H19 upstream conserved regions (HUCs). Mouse HUC1 and HUC2 act as potent enhancers capable of driving expression of an H19 reporter gene in a range of mesodermal tissues. Intriguingly, the HUC sequences are also transcribed bi-allelically in mouse and human, but their expression pattern in neural and endodermal tissues in day 13.5 embryos is distinct from their enhancer function. The location of the HUC mesodermal enhancers upstream of the ICR and H19, and their capacity for interaction with both H19 and Igf2 requires critical re-evaluation of the cis-regulation of imprinted gene expression of H19 and Igf2 in a range of mesodermal tissues. We propose that these novel sequences interact with the ICR at H19 and the epigenetically regulated silencer at differentially methylated region 1 (DMR1) of Igf2.

Epigenetic reprogramming of the genome--from the germ line to the embryo and back again

Mammalian parental genomes are not functionally equivalent, and both a maternal and paternal contribution is required for normal development. The differences between the parental genomes are the result of genomic imprinting--a form of gene regulation that results in monoallelic expression of imprinted genes. Cis-regulatory elements at imprinted loci are responsible for directing allele-specific epigenetic marks required for correct gene expression. This cis information must be interpreted at various points in development, including in the germline where existing imprints are erased and reset. Imprints must also be maintained during preimplantation development, when the genome undergoes dramatic global epigenetic changes.

A conserved imprinting control region at the HYMAI/ZAC domain is implicated in transient neonatal diabetes mellitus

Transient neonatal diabetes mellitus (TNDM) is associated with intra-uterine growth retardation, dehydration and a lack of insulin. Some TNDM patients exhibit paternal uniparental disomy (UPD) of chromosome 6q24, where at least two imprinted genes, HYMAI and ZAC, have so far been characterized. Here we show that the differentially methylated CpG island that partially overlaps mZac1 and mHymai at the syntenic mouse locus is a likely imprinting control region (ICR) for the approximately 120--200 kb domain. The region is unmethylated in sperm but probably methylated in oocytes, a difference that persists between parental alleles throughout pre- and post-implantation development. We also show that within this ICR, there is a region that exhibits a high degree of homology between mouse and human. Using a reporter expression assay, we demonstrate that this conserved region acts as a strong transcriptional repressor when methylated. Finally, we provide in vivo evidence that in the majority of TNDM patients with a normal karyotype, there is a loss of methylation within the highly homologous region. We propose that this ICR regulates expression of imprinted genes within the domain; epigenetic or genetic mutations of this region probably result in TNDM, possibly by affecting expression of ZAC in the pancreas and/or the pituitary.

Deletion of a silencer element disrupts H19 imprinting independently of a DNA methylation epigenetic switch

The H19 imprinted gene is silenced when paternally inherited and active only when inherited maternally. This is thought to involve a cis-acting control region upstream of H19 that is responsible for regulating a number of functions including DNA methylation, asynchronous replication of parental chromosomes and an insulator. Here we report on the function of a 1.2 kb upstream element in the mouse, which was previously shown to function as a bi-directional silencer in Drosophila. The cre-loxP-mediated targeted deletion of the 1.2 kb region had no effect on the maternal allele. However, there was loss of silencing of the paternal allele in many endodermal and other tissues. The pattern of expression was very similar to the expression pattern conferred by the enhancer elements downstream of H19. We could not detect an effect on the expression of the neighbouring imprinted Igf2 gene, suggesting that the proposed boundary element insulating this gene from the downstream enhancers was unaffected. Despite derepression of the paternal H19 allele, the deletion surprisingly did not affect the differential DNA methylation of the locus, which displayed an appropriate epigenetic switch in the parental germlines. Furthermore, the characteristic asynchronous pattern of DNA replication at H19 was also not disrupted by the deletion, suggesting that the sequences that mediate this were also intact. The silencer is therefore part of a complex cis-regulatory region upstream of the H19 gene and acts specifically to ensure the repression of the paternal allele, without a predominant effect on the epigenetic switch in the germline.

A novel imprinted gene, HYMAI, is located within an imprinted domain on human chromosome 6 containing ZAC

Transient neonatal diabetes mellitus (TNDM) is a rare disease characterized by intrauterine growth retardation, dehydration, and failure to thrive due to a lack of normal insulin secretion. This disease is associated with paternal uniparental disomy or paternal duplication of chromosome 6, suggesting that the causative gene(s) for TNDM is imprinted. Recently, Gardner et al. (1999, J. Med. Genet. 36: 192-196) proposed that a candidate gene for TNDM lies within chromosome 6q24.1-q24.3. To find human imprinted genes, we performed a database search for EST sequences that mapped to this region, followed by RT-PCR analysis using monochromosomal hybrid cells with a human chromosome 6 of defined parental origin. Here we report the identification of a novel imprinted gene, HYMAI. This gene exhibits differential DNA methylation between the two parental alleles at an adjacent CpG island and is expressed only from the paternal chromosome. A previously characterized imprinted gene, ZAC/LOT1, is located 70 kb downstream of HYMAI and is also expressed only from the paternal allele. In the pancreas, both genes are moderately expressed. HYMAI and ZAC/LOT1 are therefore candidate genes involved in TNDM. Furthermore, the human chromosome 6q24 region is syntenic to mouse chromosome 10 and represents a novel imprinted domain.

A silencer element identified in Drosophila is required for imprinting of H19 reporter transgenes in mice

The H19 gene is subject to genomic imprinting because it is methylated and repressed after paternal inheritance and is unmethylated and expressed after maternal inheritance. We recently identified a 1.1-kb control element in the upstream region of the H19 gene that functions as a cis-acting silencer element in Drosophila. Here we investigate the function of this element in mice. We demonstrate that both H19-lacZ and H19-PLAP reporter transgenes can undergo imprinting with repression and hypermethylation after paternal transmission at many integration sites. However, transgenes that were deleted for the 1.1-kb silencer element showed loss of paternal repression, but they did not show marked changes in the paternal methylation of the remaining upstream region. This study demonstrates that the 1.1-kb control element identified in Drosophila is required to silence paternally transmitted H19 minitransgenes in mice.