Comprehensive identification of Drosophila dorsal-ventral patterning genes using a whole-genome tiling array

Prediction and analysis of cis-regulatory elements in Dorsal and Ventral patterning genes of Tribolium castaneum and its comparison with Drosophila melanogaster

Insect embryonic development and morphology are characterized by their anterior-posterior and dorsal-ventral (DV) patterning. In Drosophila embryos, DV patterning is mediated by a dorsal protein gradient which activates twist and snail proteins, the important regulators of DV patterning. To activate or repress gene expression, some regulatory proteins bind in clusters to their target gene at sites known as cis-regulatory elements or enhancers. To understand how variations in gene expression in different lineages might lead to different phenotypes, it is necessary to understand enhancers and their evolution. Drosophila melanogaster has been widely studied to understand the interactions between transcription factors and the transcription factor binding sites. Tribolium castaneum is an upcoming model animal which is catching the interest of biologists and the research on the enhancer mechanisms in the insect's axes patterning is still in infancy. Therefore, the current study was designed to compare the enhancers of DV patterning in the two insect species. The sequences of ten proteins involved in DV patterning of D. melanogaster were obtained from Flybase. The protein sequences of T. castaneum orthologous to those obtained from D. melanogaster were acquired from NCBI BLAST, and these were then converted to DNA sequences which were modified by adding 20 kb sequences both upstream and downstream to the gene. These modified sequences were used for further analysis. Bioinformatics tools (Cluster-Buster and MCAST) were used to search for clusters of binding sites (enhancers) in the modified DV genes. The results obtained showed that the transcription factors in Drosophila melanogaster and Tribolium castaneum are nearly identical; however, the number of binding sites varies between the two species, indicating transcription factor binding site evolution, as predicted by two different computational tools. It was observed that dorsal, twist, snail, zelda, and Supressor of Hairless are the transcription factors responsible for the regulation of DV patterning in the two insect species.

miR-9a regulates levels of both rhomboid mRNA and protein in the early Drosophila melanogaster embryo

MicroRNAs can have subtle and combinatorial effects on the levels of the targets and pathways they act on. Studying the consequences of a single microRNA knockout often proves difficult as many such knockouts exhibit phenotypes only under stress conditions. This has often led to the hypothesis that microRNAs buffer the effects of intrinsic and environmental stochasticity on gene expression. Observing and understanding this buffering effect entails quantitative analysis of microRNA and target expression in single cells. To this end, we have employed single-molecule fluorescence in situ hybridization, immunofluorescence, and high-resolution confocal microscopy to investigate the effects of miR-9a loss on the expression of the serine-protease Rhomboid in Drosophila melanogaster early embryos. Our single-cell quantitative approach shows that spatially, the rhomboid mRNA pattern is identical in WT and miR-9a knockout embryos. However, we find that the number of mRNA molecules per cell is higher when miR-9a is absent, and the level and temporal accumulation of rhomboid protein shows a more dramatic increase in the miR-9a knockout. Specifically, we see accumulation of rhomboid protein in miR-9a mutants by stage 5, much earlier than in WT. The data, therefore, show that miR-9a functions in the regulation of rhomboid mRNA and protein levels. While further work is required to establish whether rhomboid is a direct target of miR-9 in Drosophila, our results further establish the miR-9 family microRNAs as conserved regulators of timing in neurogenic processes. This study shows the power of single-cell quantification as an experimental tool to study phenotypic consequences of microRNA mis-regulation.

Dynamic patterning by morphogens illuminated by cis-regulatory studies

Morphogen concentration changes in space as well as over time during development. However, how these dynamics are interpreted by cells to specify fate is not well understood. Here, we focus on two morphogens: the maternal transcription factors Bicoid and Dorsal, which directly regulate target genes to pattern Drosophila embryos. The actions of these factors at enhancers has been thoroughly dissected and provides a rich platform for understanding direct input by morphogens and their changing roles over time. Importantly, Bicoid and Dorsal do not work alone; we also discuss additional inputs that work with morphogens to control spatiotemporal gene expression in embryos.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

Transcriptome analysis to characterize the genes related to gonad growth and fatty acid metabolism in the sea urchin Strongylocentrotus intermedius

BACKGROUND

Sea urchin gonads of both sexes, commonly termed "roe", are highly valued seafood delicacies, and Strongylocentrotus intermedius is considered one of the tastiest sea urchins. In order to produce high-quality gonads for consumption and clarify the mechanism of gonad growth and development of the sea urchin, more genetic information, especially at the transcriptome level, is needed.

OBJECTIVE

A more thorough understanding of sea urchin gonad growth and development in both sexes could enable regulation of these processes at several stages with the aim of suppressing gametogenesis in order to produce high-quality gonads for consumption.

METHODS

The adult sea urchins S. intermedius were cultured for 3 months, and were sampled for the gonadal transcriptome analysis which has been performed on the RNAs of three male and female adults of S. intermedius in each gonad development stage.

RESULTS

Illumina sequencing raw sequence data was deposited in the NCBI Sequence Read Archive (SRA) database (PRJNA532998). It generated 560,196,356 raw reads and 548,956,944 clean reads were acquired, which were assembled into 107,850 transcripts with 44,124 genes. Comparative analysis showed the differentially expressed genes (DEGs) from 114 to 2566. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to determine the functional significance of these DEGs. We have selected 9 genes related to growth and 12 genes related to fatty acid biosynthesis and metabolism in sea urchin gonads.

CONCLUSION

These data for sea urchins were intended to provide markers for gonad growth and development that can be accumulated for use in aquaculture applications.

Mechanisms of Groucho-mediated repression revealed by genome-wide analysis of Groucho binding and activity

Background

The transcriptional corepressor Groucho (Gro) is required for the function of many developmentally regulated DNA binding repressors, thus helping to define the gene expression profile of each cell during development. The ability of Gro to repress transcription at a distance together with its ability to oligomerize and bind to histones has led to the suggestion that Gro may spread along chromatin. However, much is unknown about the mechanism of Gro-mediated repression and about the dynamics of Gro targeting.

Results

Our chromatin immunoprecipitation sequencing analysis of temporally staged Drosophila embryos shows that Gro binds in a highly dynamic manner primarily to clusters of discrete (<1 kb) segments. Consistent with the idea that Gro may facilitate communication between silencers and promoters, Gro binding is enriched at both cis-regulatory modules, as well as within the promotors of potential target genes. While this Gro-recruitment is required for repression, our data show that it is not sufficient for repression. Integration of Gro binding data with transcriptomic analysis suggests that, contrary to what has been observed for another Gro family member, Drosophila Gro is probably a dedicated repressor. This analysis also allows us to define a set of high confidence Gro repression targets. Using publically available data regarding the physical and genetic interactions between these targets, we are able to place them in the regulatory network controlling development. Through analysis of chromatin associated pre-mRNA levels at these targets, we find that genes regulated by Gro in the embryo are enriched for characteristics of promoter proximal paused RNA polymerase II.

Conclusions

Our findings are inconsistent with a one-dimensional spreading model for long-range repression and suggest that Gro-mediated repression must be regulated at a post-recruitment step. They also show that Gro is likely a dedicated repressor that sits at a prominent highly interconnected regulatory hub in the developmental network. Furthermore, our findings suggest a role for RNA polymerase II pausing in Gro-mediated repression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3589-6) contains supplementary material, which is available to authorized users.

Genome-wide identification of Drosophila dorso-ventral enhancers by differential histone acetylation analysis

BACKGROUND

Drosophila dorso-ventral (DV) patterning is one of the best-understood regulatory networks to date, and illustrates the fundamental role of enhancers in controlling patterning, cell fate specification, and morphogenesis during development. Histone acetylation such as H3K27ac is an excellent marker for active enhancers, but it is challenging to obtain precise locations for enhancers as the highest levels of this modification flank the enhancer regions. How to best identify tissue-specific enhancers in a developmental system de novo with a minimal set of data is still unclear.

RESULTS

Using DV patterning as a test system, we develop a simple and effective method to identify tissue-specific enhancers de novo. We sample a broad set of candidate enhancer regions using data on CREB-binding protein co-factor binding or ATAC-seq chromatin accessibility, and then identify those regions with significant differences in histone acetylation between tissues. This method identifies hundreds of novel DV enhancers and outperforms ChIP-seq data of relevant transcription factors when benchmarked with mRNA expression data and transgenic reporter assays. These DV enhancers allow the de novo discovery of the relevant transcription factor motifs involved in DV patterning and contain additional motifs that are evolutionarily conserved and for which the corresponding transcription factors are expressed in a DV-biased fashion. Finally, we identify novel target genes of the regulatory network, implicating morphogenesis genes as early targets of DV patterning.

CONCLUSIONS

Taken together, our approach has expanded our knowledge of the DV patterning network even further and is a general method to identify enhancers in any developmental system, including mammalian development.

Global analysis of dorsoventral patterning in the wasp Nasonia reveals extensive incorporation of novelty in a regulatory network

Background

Gene regulatory networks (GRNs) underlie developmental patterning and morphogenetic processes, and changes in the interactions within the underlying GRNs are a major driver of evolutionary processes. In order to make meaningful comparisons that can provide significant insights into the evolution of regulatory networks, homologous networks from multiple taxa must be deeply characterized. One of the most thoroughly characterized GRNs is the dorsoventral (DV) patterning system of the Drosophila melanogaster embryo. We have developed the wasp Nasonia as a comparative DV patterning model because it has shown the convergent evolution of a mode of early embryonic patterning very similar to that of the fly, and it is of interest to know whether the similarity at the gross level also extends to the molecular level.

Results

We used RNAi to dorsalize and ventralize Nasonia embryos, RNAseq to quantify transcriptome-wide expression levels, and differential expression analysis to identify genes whose expression levels change in either RNAi case. This led to the identification of >100 genes differentially expressed and regulated along the DV axis. Only a handful of these genes are shared DV components in both fly and wasp. Many of those unique to Nasonia are cytoskeletal and adhesion molecules, which may be related to the divergent cell and tissue behavior observed at gastrulation. In addition, many transcription factors and signaling components are only DV regulated in Nasonia, likely reflecting the divergent upstream patterning mechanisms involved in producing the conserved pattern of cell fates observed at gastrulation. Finally, several genes that lack Drosophila orthologs show robust and distinct expression patterns. These include genes with vertebrate homologs that have been lost in the fly lineage, genes that are found only among Hymenoptera, and several genes that entered the Nasonia genome through lateral transfer from endosymbiotic bacteria.

Conclusions

Altogether, our results provide insights into how GRNs respond to new functional demands and how they can incorporate novel components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0285-y) contains supplementary material, which is available to authorized users.

Stepwise Progression of Embryonic Patterning

It is long established that the graded distribution of Dorsal transcription factor influences spatial domains of gene expression along the dorsoventral (DV) axis of Drosophila melanogaster embryos. However, the more recent realization that Dorsal levels also change with time raises the question of whether these dynamics are instructive. An overview of DV axis patterning is provided, focusing on new insights identified through quantitative analysis of temporal changes in Dorsal target gene expression from one nuclear cycle to the next ('steps'). Possible roles for the stepwise progression of this patterning program are discussed including (i) tight temporal regulation of signaling pathway activation, (ii) control of gene expression cohorts, and (iii) ensuring the irreversibility of the patterning and cell fate specification process.

From blastocyst to gastrula: gene regulatory networks of embryonic stem cells and early mouse embryogenesis

To date, many regulatory genes and signalling events coordinating mammalian development from blastocyst to gastrulation stages have been identified by mutational analyses and reverse-genetic approaches, typically on a gene-by-gene basis. More recent studies have applied bioinformatic approaches to generate regulatory network models of gene interactions on a genome-wide scale. Such models have provided insights into the gene networks regulating pluripotency in embryonic and epiblast stem cells, as well as cell-lineage determination in vivo. Here, we review how regulatory networks constructed for different stem cell types relate to corresponding networks in vivo and provide insights into understanding the molecular regulation of the blastocyst-gastrula transition.

Rice stripe virus affects the viability of its vector offspring by changing developmental gene expression in embryos

Plant viruses may affect the viability and development process of their herbivore vectors. Small brown planthopper (SBPH) is main vector of Rice stripe virus (RSV), which causes serious rice stripe disease. Here, we reported the effects of RSV on SBPH offspring by crossing experiments between viruliferous and non-viruliferous strains. The life parameters of offspring from different cross combinations were compared. The hatchability of F1 progeny from viruliferous parents decreased significantly, and viruliferous rate was completely controlled by viruliferous maternal parent. To better elucidate the underlying biological mechanisms, the morphology of eggs, viral propagation and distribution in the eggs and expression profile of embryonic development genes were investigated. The results indicated that RSV replicated and accumulated in SBPH eggs resulting in developmental stunt or delay of partial eggs; in addition, RSV was only able to infect ovum but not sperm. According to the expression profile, expression of 13 developmental genes was regulated in the eggs from viruliferous parents, in which two important regulatory genes (Ls-Dorsal and Ls-CPO) were most significantly down-regulated. In general, RSV exerts an adverse effect on SBPH, which is unfavourable for the expansion of viruliferous populations. The viewpoint is also supported by systematic monitoring of SBPH viruliferous rate.

Toll signals regulate dorsal-ventral patterning and anterior-posterior placement of the embryo in the hemipteran Rhodnius prolixus

BACKGROUND

Insect embryonic dorso-ventral patterning depends greatly on two pathways: the Toll pathway and the Bone Morphogenetic Protein pathway. While the relative contribution of each pathway has been investigated in holometabolous insects, their role has not been explored in insects with a hemimetabolous type of development. The hemimetabolous insect Rhodnius prolixus, an important vector of Chagas disease in the Americas, develops from an intermediate germ band and displays complex movements during katatrepsis that are not observed in other orders. However, little is known about the molecular events that regulate its embryogenesis. Here we investigate the expression and function of genes potentially involved in the initial patterning events that establish the embryonic dorso-ventral axis in this hemipteran.

RESULTS

We establish a staging system for early embryogenesis that allows us to correlate embryo morphology with gene expression profiles. Using this system, we investigate the role of Toll pathway genes during embryogenesis. Detailed analyses of gene expression throughout development, coupled with functional analyses using parental RNA interference, revealed that maternal Toll is required to establish germ layers along the dorso-ventral axis and for embryo placement along the anterior-posterior axis. Interestingly, knockdown of the Toll pathway effector Rp-dorsal appears to regulate the expression of the Bone Morphogenetic Protein antagonist Rp-short-gastrulation.

CONCLUSIONS

Our results indicate that Toll signals are the initiating event in dorso-ventral patterning during Rhodnius embryogenesis, and this is the first report of a conserved role for Toll in a hemipteran. Furthermore, as Rp-dorsal RNA interference generates anteriorly misplaced embryos, our results indicate a novel role for Toll signals in establishment of the anterior-posterior axis in Rhodnius.

Maternal control of the Drosophila dorsal-ventral body axis

The pathway that generates the dorsal-ventral (DV) axis of the Drosophila embryo has been the subject of intense investigation over the previous three decades. The initial asymmetric signal originates during oogenesis by the movement of the oocyte nucleus to an anterior corner of the oocyte, which establishes DV polarity within the follicle through signaling between Gurken, the Drosophila Transforming Growth Factor (TGF)-α homologue secreted from the oocyte, and the Drosophila Epidermal Growth Factor Receptor (EGFR) that is expressed by the follicular epithelium cells that envelop the oocyte. Follicle cells that are not exposed to Gurken follow a ventral fate and express Pipe, a sulfotransferase that enzymatically modifies components of the inner vitelline membrane layer of the eggshell, thereby transferring DV spatial information from the follicle to the egg. These ventrally sulfated eggshell proteins comprise a localized cue that directs the ventrally restricted formation of the active Spätzle ligand within the perivitelline space between the eggshell and the embryonic membrane. Spätzle activates Toll, a transmembrane receptor in the embryonic membrane. Transmission of the Toll signal into the embryo leads to the formation of a ventral-to-dorsal gradient of the transcription factor Dorsal within the nuclei of the syncytial blastoderm stage embryo. Dorsal controls the spatially specific expression of a large constellation of zygotic target genes, the Dorsal gene regulatory network, along the embryonic DV circumference. This article reviews classic studies and integrates them with the details of more recent work that has advanced our understanding of the complex pathway that establishes Drosophila embryo DV polarity. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

A graphical model method for integrating multiple sources of genome-scale data

Making effective use of multiple data sources is a major challenge in modern bioinformatics. Genome-wide data such as measures of transcription factor binding, gene expression, and sequence conservation, which are used to identify binding regions and genes that are important to major biological processes such as development and disease, can be difficult to use together due to the different biological meanings and statistical distributions of the heterogeneous data types, but each can provide valuable information for understanding the processes under study. Here we present methods for integrating multiple data sources to gain a more complete picture of gene regulation and expression. Our goal is to identify genes and cis-regulatory regions which play specific biological roles. We describe a graphical mixture model approach for data integration, examine the effect of using different model topologies, and discuss methods for evaluating the effectiveness of the models. Model fitting is computationally efficient and produces results which have clear biological and statistical interpretations. The Hedgehog and Dorsal signaling pathways in Drosophila, which are critical in embryonic development, are used as examples.

Inference of the Xenopus tropicalis embryonic regulatory network and spatial gene expression patterns

BACKGROUND

During embryogenesis, signaling molecules produced by one cell population direct gene regulatory changes in neighboring cells and influence their developmental fates and spatial organization. One of the earliest events in the development of the vertebrate embryo is the establishment of three germ layers, consisting of the ectoderm, mesoderm and endoderm. Attempts to measure gene expression in vivo in different germ layers and cell types are typically complicated by the heterogeneity of cell types within biological samples (i.e., embryos), as the responses of individual cell types are intermingled into an aggregate observation of heterogeneous cell types. Here, we propose a novel method to elucidate gene regulatory circuits from these aggregate measurements in embryos of the frog Xenopus tropicalis using gene network inference algorithms and then test the ability of the inferred networks to predict spatial gene expression patterns.

RESULTS

We use two inference models with different underlying assumptions that incorporate existing network information, an ODE model for steady-state data and a Markov model for time series data, and contrast the performance of the two models. We apply our method to both control and knockdown embryos at multiple time points to reconstruct the core mesoderm and endoderm regulatory circuits. Those inferred networks are then used in combination with known dorsal-ventral spatial expression patterns of a subset of genes to predict spatial expression patterns for other genes. Both models are able to predict spatial expression patterns for some of the core mesoderm and endoderm genes, but interestingly of different gene subsets, suggesting that neither model is sufficient to recapitulate all of the spatial patterns, yet they are complementary for the patterns that they do capture.

CONCLUSION

The presented methodology of gene network inference combined with spatial pattern prediction provides an additional layer of validation to elucidate the regulatory circuits controlling the spatial-temporal dynamics in embryonic development.

A genetic network conferring canalization to a bistable patterning system in Drosophila

To achieve the "constancy of the wild-type," the developing organism must be buffered against stochastic fluctuations and environmental perturbations. This phenotypic buffering has been theorized to arise from a variety of genetic mechanisms and is widely thought to be adaptive and essential for viability. In the Drosophila blastoderm embryo, staining with antibodies against the active, phosphorylated form of the bone morphogenetic protein (BMP) signal transducer Mad, pMad, or visualization of the spatial pattern of BMP-receptor interactions reveals a spatially bistable pattern of BMP signaling centered on the dorsal midline. This signaling event is essential for the specification of dorsal cell fates, including the extraembryonic amnioserosa. BMP signaling is initiated by facilitated extracellular diffusion that localizes BMP ligands dorsally. BMP signaling then activates an intracellular positive feedback circuit that promotes future BMP-receptor interactions. Here, we identify a genetic network comprising three genes that canalizes this BMP signaling event. The BMP target eiger (egr) acts in the positive feedback circuit to promote signaling, while the BMP binding protein encoded by crossveinless-2 (cv-2) antagonizes signaling. Expression of both genes requires the early activity of the homeobox gene zerknüllt (zen). Two Drosophila species lacking early zen expression have high variability in BMP signaling. These data both detail a new mechanism that generates developmental canalization and identify an example of a species with noncanalized axial patterning.

Temporal regulation of single-minded target genes in the ventral midline of the Drosophila central nervous system

Differentiation of a specific organ or tissue requires sequential activation of regulatory genes. However, little is known about how serial gene expression is temporally regulated. Here, we present evidence that differential expression of single-minded (sim) target genes can be attributed, in part, to the number of Sim and Tango (Tgo) heterodimer binding sites within their enhancer regions. The Sim, termed a master regulator, directs ventral midline differentiation of Drosophila central nervous system (CNS). According to data on the onset timing of ventral midline gene expression, sim target genes are classified into at least 2 groups (early and late). The sim and rhomboid (rho) genes are activated during early midline differentiation whereas orthodenticle (otd), CG10249, and slit (sli) genes undergo activation during later stages of midline differentiation. Germline transformation and in situ hybridization with transgenic embryos demonstrate that enhancers activating sim and rho expression contain 4 Sim-Tgo binding sites whereas only 1 Sim-Tgo binding site is found in an enhancer of sli. A mutagenized version of the rho enhancer lacking either 1, 2, or 3 Sim-Tgo binding sites mediated progressively more delayed expression of a lacZ reporter gene in the ventral midline. In contrast, a modified sli enhancer displayed progressively earlier onset of lacZ expression when 1, 2, or 3 more Sim-Tgo binding sites were added. Taken together, these results suggest that the number of Sim-Tgo-binding sites is decisive in determining the timing of gene expression in the developing ventral midline. We also discuss a combinatorial model accounting for the sequential expression of sim target genes.

RNAi reveals the key role of Nervana 1 in cockroach oogenesis and embryo development

Na(+), K(+)-ATPases is a heterodimer protein consisting of α- and β-subunits that control the ion transport through cell membranes. In insects the β-subunit of the Na(+), K(+)-ATPase, known as Nervana, was characterized as a nervous system-specific glycoprotein antigen from adult Drosophila melanogaster heads. Nervana is expressed ubiquitously in all insect tissues, and in epithelial cells appeared located in a basolateral position as part of the septate junctions. Herein we study two Nervana isoforms from Blattella germanica, a cockroach species with panoistic ovaries. The sequencing and the phylogenetic analysis results suggest that these two isoforms are orthologs of D. melanogaster Nervana 1 and Nervana 2, respectively. Nervana 1 is highly expressed in the ovary of B. germanica, and depleting its expression results in changes in oocyte shape that do not impair oviposition. However, the resulting embryos show different defects and never hatch. These findings highlight the importance of this type of membrane pump in insect oogenesis as well as in embryo development, and its possible regulation by juvenile hormone.

Deep sequencing and expression of microRNAs from early honeybee (Apis mellifera) embryos reveals a role in regulating early embryonic patterning

Background

Recent evidence supports the proposal that the observed diversity of animal body plans has been produced through alterations to the complexity of the regulatory genome rather than increases in the protein-coding content of a genome. One significant form of gene regulation is the contribution made by the non-coding content of the genome. Non-coding RNAs play roles in embryonic development of animals and these functions might be expected to evolve rapidly. Using next-generation sequencing and in situ hybridization, we have examined the miRNA content of early honeybee embryos.

Results

Through small RNA sequencing we found that 28% of known miRNAs are expressed in the early embryo. We also identified developmentally expressed microRNAs that are unique to the Apoidea clade. Examination of expression patterns implied these miRNAs have roles in patterning the anterior-posterior and dorso-ventral axes as well as the extraembryonic membranes. Knockdown of Dicer, a key component of miRNA processing, confirmed that miRNAs are likely to have a role in patterning these tissues.

Conclusions

Examination of the expression patterns of novel miRNAs, some unique to the Apis group, indicated that they are likely to play a role in early honeybee development. Known miRNAs that are deeply conserved in animal phyla display differences in expression pattern between honeybee and Drosophila, particularly at early stages of development. This may indicate miRNAs play a rapidly evolving role in regulating developmental pathways, most likely through changes to the way their expression is regulated.

NF-κB: where did it come from and why?

The vast majority of research on nuclear factor κB (NF-κB) signaling in the past 25 years has focused on its roles in normal and disease-related processes in vertebrates, especially mice and humans. Recent genome and transcriptome sequencing efforts have shown that homologs of NF-κB transcription factors, inhibitor of NF-κB (IκB) proteins, and IκB kinases are present in a variety of invertebrates, including several in phyla simpler than Arthropoda, the phylum containing insects such Drosophila. Moreover, many invertebrates also contain genes encoding homologs of upstream signaling proteins in the Toll-like receptor signaling pathway, which is well-known for its downstream activation of NF-κB for innate immunity. This review describes what we now know or can infer and speculate about the evolution of the core elements of NF-κB signaling as well as the biological processes controlled by NF-κB in invertebrates. Further research on NF-κB in invertebrates is likely to uncover information about the evolutionary origins of this key human signaling pathway and may have relevance to our management of the responses of ecologically and economically important organisms to environmental and adaptive pressures.

[Wnt signaling pathway and the Evo-Devo of deuterostome axis]

A series of signal transduction pathways have been found to regulate the polarity establishment and formation of animal primary body axis. Among them, Wnt signaling pathway is extremely conserved and several key components in the pathway have been identified in the demosponge lineage. This implies that it is one of the earliest pathways involved in the ancestral metazoan axis development and might play an important role in specification and development of posterior and ventral fate of animal axis. Recently, with the establishment of functional experiments in vitro, the body plan formation has been found to be affected, in varying degrees, by many genes in the Wnt signaling pathway, such as members of wnt gene family, maternal gene beta-catenin and some transcription factor encoding genes. In this review, we analyzed the evolutionary origin of the wnt gene family involved in development of metazoan body plans, and then made a brief review on the roles of canonical Wnt/beta-catenin signaling in the polarity establishment and formation of primary body axis in diverse deuterostomes including sea urchin, amphioxus, zebrafish, frog, and mouse.

Temporal coordination of gene networks by Zelda in the early Drosophila embryo

In past years, much attention has focused on the gene networks that regulate early developmental processes, but less attention has been paid to how multiple networks and processes are temporally coordinated. Recently the discovery of the transcriptional activator Zelda (Zld), which binds to CAGGTAG and related sequences present in the enhancers of many early-activated genes in Drosophila, hinted at a mechanism for how batteries of genes could be simultaneously activated. Here we use genome-wide binding and expression assays to identify Zld target genes in the early embryo with the goal of unraveling the gene circuitry regulated by Zld. We found that Zld binds to genes involved in early developmental processes such as cellularization, sex determination, neurogenesis, and pattern formation. In the absence of Zld, many target genes failed to be activated, while others, particularly the patterning genes, exhibited delayed transcriptional activation, some of which also showed weak and/or sporadic expression. These effects disrupted the normal sequence of patterning-gene interactions and resulted in highly altered spatial expression patterns, demonstrating the significance of a timing mechanism in early development. In addition, we observed prevalent overlap between Zld-bound regions and genomic "hotspot" regions, which are bound by many developmental transcription factors, especially the patterning factors. This, along with the finding that the most over-represented motif in hotspots, CAGGTA, is the Zld binding site, implicates Zld in promoting hotspot formation. We propose that Zld promotes timely and robust transcriptional activation of early-gene networks so that developmental events are coordinated and cell fates are established properly in the cellular blastoderm embryo.

Role of sequence encoded κB DNA geometry in gene regulation by Dorsal

Many proteins of the Rel family can act as both transcriptional activators and repressors. However, mechanism that discerns the ‘activator/repressor’ functions of Rel-proteins such as Dorsal (Drosophila homologue of mammalian NFκB) is not understood. Using genomic, biophysical and biochemical approaches, we demonstrate that the underlying principle of this functional specificity lies in the ‘sequence-encoded structure’ of the κB-DNA. We show that Dorsal-binding motifs exist in distinct activator and repressor conformations. Molecular dynamics of DNA-Dorsal complexes revealed that repressor κB-motifs typically have A-tract and flexible conformation that facilitates interaction with co-repressors. Deformable structure of repressor motifs, is due to changes in the hydrogen bonding in A:T pair in the ‘A-tract’ core. The sixth nucleotide in the nonameric κB-motif, ‘A’ (A₆) in the repressor motifs and ‘T’ (T₆) in the activator motifs, is critical to confer this functional specificity as A₆ → T₆ mutation transformed flexible repressor conformation into a rigid activator conformation. These results highlight that ‘sequence encoded κB DNA-geometry’ regulates gene expression by exerting allosteric effect on binding of Rel proteins which in turn regulates interaction with co-regulators. Further, we identified and characterized putative repressor motifs in Dl-target genes, which can potentially aid in functional annotation of Dorsal gene regulatory network.

Graded dorsal and differential gene regulation in the Drosophila embryo

A gradient of Dorsal activity patterns the dorsoventral (DV) axis of the early Drosophila melanogaster embryo by controlling the expression of genes that delineate presumptive mesoderm, neuroectoderm, and dorsal ectoderm. The availability of the Drosophila melanogaster genome sequence has accelerated the study of embryonic DV patterning, enabling the use of systems-level approaches. As a result, our understanding of Dorsal-dependent gene regulation has expanded to encompass a collection of more than 50 genes and 30 cis-regulatory sequences. This information, which has been integrated into a spatiotemporal atlas of gene regulatory interactions, comprises one of the best-understood networks controlling any developmental process to date. In this article, we focus on how Dorsal controls differential gene expression and how recent studies have expanded our understanding of Drosophila embryonic development from the cis-regulatory level to that controlling morphogenesis of the embryo.

Dynamic repositioning of dorsal to two different kappaB motifs controls its autoregulation during immune response in Drosophila

Autoregulation is one of the mechanisms of imparting feedback control on gene expression. Positive autoregulatory feedback results in induction of a gene, and negative feedback leads to its suppression. Here, we report an interesting mechanism of autoregulation operating on Drosophila Rel gene dorsal that can activate as well as repress its expression. Using biochemical and genetic approaches, we show that upon immune challenge Dorsal regulates its activation as well as repression by dynamically binding to two different kappaB motifs, kappaB(I) (intronic kappaB) and kappaB(P) (promoter kappaB), present in the dorsal gene. Although the kappaB(I) motif functions as an enhancer, the kappaB(P) motif acts as a transcriptional repressor. Interestingly, Dorsal binding to these two motifs is dynamic; immediately upon immune challenge, Dorsal binds to the kappaB(I) leading to auto-activation, whereas at the terminal phase of the immune response, it is removed from the kappaB(I) and repositioned at the kappaB(P), resulting in its repression. Furthermore, we show that repression of Dorsal as well as its binding to the kappaB(P) depends on the transcription factor AP1. Depletion of AP1 by RNA interference resulted in constitutive expression of Dorsal. In conclusion, this study suggests that during acute phase response dorsal is regulated by following two subcircuits: (i) Dl-kappaB(I) for activation and (ii) Dl-AP1-kappaB(P) for repression. These two subcircuits are temporally delineated and bring about overall regulation of dorsal during immune response. These results suggest the presence of a previously unknown mechanism of Dorsal autoregulation in immune-challenged Drosophila.

Analysis of antisense expression by whole genome tiling microarrays and siRNAs suggests mis-annotation of Arabidopsis orphan protein-coding genes

Background

MicroRNAs (miRNAs) and trans-acting small-interfering RNAs (tasi-RNAs) are small (20–22 nt long) RNAs (smRNAs) generated from hairpin secondary structures or antisense transcripts, respectively, that regulate gene expression by Watson-Crick pairing to a target mRNA and altering expression by mechanisms related to RNA interference. The high sequence homology of plant miRNAs to their targets has been the mainstay of miRNA prediction algorithms, which are limited in their predictive power for other kingdoms because miRNA complementarity is less conserved yet transitive processes (production of antisense smRNAs) are active in eukaryotes. We hypothesize that antisense transcription and associated smRNAs are biomarkers which can be computationally modeled for gene discovery.

Principal Findings

We explored rice (Oryza sativa) sense and antisense gene expression in publicly available whole genome tiling array transcriptome data and sequenced smRNA libraries (as well as C. elegans) and found evidence of transitivity of MIRNA genes similar to that found in Arabidopsis. Statistical analysis of antisense transcript abundances, presence of antisense ESTs, and association with smRNAs suggests several hundred Arabidopsis ‘orphan’ hypothetical genes are non-coding RNAs. Consistent with this hypothesis, we found novel Arabidopsis homologues of some MIRNA genes on the antisense strand of previously annotated protein-coding genes. A Support Vector Machine (SVM) was applied using thermodynamic energy of binding plus novel expression features of sense/antisense transcription topology and siRNA abundances to build a prediction model of miRNA targets. The SVM when trained on targets could predict the “ancient” (deeply conserved) class of validated Arabidopsis MIRNA genes with an accuracy of 84%, and 76% for “new” rapidly-evolving MIRNA genes.

Conclusions

Antisense and smRNA expression features and computational methods may identify novel MIRNA genes and other non-coding RNAs in plants and potentially other kingdoms, which can provide insight into antisense transcription, miRNA evolution, and post-transcriptional gene regulation.

Intergenic, gene terminal, and intragenic CpG islands in the human genome

Background

Recently, it has been discovered that the human genome contains many transcription start sites for non-coding RNA. Regulatory regions related to transcription of this non-coding RNAs are poorly studied. Some of these regulatory regions may be associated with CpG islands located far from transcription start-sites of any protein coding gene. The human genome contains many such CpG islands; however, until now their properties were not systematically studied.

Results

We studied CpG islands located in different regions of the human genome using methods of bioinformatics and comparative genomics. We have observed that CpG islands have a preference to overlap with exons, including exons located far from transcription start site, but usually extend well into introns. Synonymous substitution rate of CpG-containing codons becomes substantially reduced in regions where CpG islands overlap with protein-coding exons, even if they are located far downstream from transcription start site. CAGE tag analysis displayed frequent transcription start sites in all CpG islands, including those found far from transcription start sites of protein coding genes. Computational prediction and analysis of published ChIP-chip data revealed that CpG islands contain an increased number of sites recognized by Sp1 protein. CpG islands containing more CAGE tags usually also contain more Sp1 binding sites. This is especially relevant for CpG islands located in 3' gene regions. Various examples of transcription, confirmed by mRNAs or ESTs, but with no evidence of protein coding genes, were found in CAGE-enriched CpG islands located far from transcription start site of any known protein coding gene.

Conclusions

CpG islands located far from transcription start sites of protein coding genes have transcription initiation activity and display Sp1 binding properties. In exons, overlapping with these islands, the synonymous substitution rate of CpG containing codons is decreased. This suggests that these CpG islands are involved in transcription initiation, possibly of some non-coding RNAs.

Genes encoding novel secreted and transmembrane proteins are temporally and spatially regulated during Drosophila melanogaster embryogenesis

BACKGROUND

Morphogenetic events that shape the Drosophila melanogaster embryo are tightly controlled by a genetic program in which specific sets of genes are up-regulated. We used a suppressive subtractive hybridization procedure to identify a group of developmentally regulated genes during early stages of D. melanogaster embryogenesis. We studied the spatiotemporal activity of these genes in five different intervals covering 12 stages of embryogenesis.

RESULTS

Microarrays were constructed to confirm induction of expression and to determine the temporal profile of isolated subtracted cDNAs during embryo development. We identified a set of 118 genes whose expression levels increased significantly in at least one developmental interval compared with a reference interval. Of these genes, 53% had a phenotype and/or molecular function reported in the literature, whereas 47% were essentially uncharacterized. Clustering analysis revealed demarcated transcript groups with maximum gene activity at distinct developmental intervals. In situ hybridization assays were carried out on 23 uncharacterized genes, 15 of which proved to have spatiotemporally restricted expression patterns. Among these 15 uncharacterized genes, 13 were found to encode putative secreted and transmembrane proteins. For three of them we validated our protein sequence predictions by expressing their cDNAs in Drosophila S2R+ cells and analyzed the subcellular distribution of recombinant proteins. We then focused on the functional characterization of the gene CG6234. Inhibition of CG6234 by RNA interference resulted in morphological defects in embryos, suggesting the involvement of this gene in germ band retraction.

CONCLUSION

Our data have yielded a list of developmentally regulated D. melanogaster genes and their expression profiles during embryogenesis and provide new information on the spatiotemporal expression patterns of several uncharacterized genes. In particular, we recovered a substantial number of unknown genes encoding putative secreted and transmembrane proteins, suggesting new components of signaling pathways that might be incorporated within the existing regulatory networks controlling D. melanogaster embryogenesis. These genes are also good candidates for additional targeted functional analyses similar to those we conducted for CG6234.See related minireview by Vichas and Zallen: http://www.jbiol.com/content/8/8/76.

Combinatorial patterning mechanisms in the Drosophila embryo

The classical concept of the morphogen gradient proposes that small differences in the levels of a signalling molecule or transcription factor are responsible for producing a continuous spectrum of distinctive cellular identities across a naïve field of cells. In this review, we discuss how the Dorsal gradient controls the dorsal-ventral patterning of the early Drosophila embryo. This gradient extends from the ventral midline of the embryo into dorso-lateral regions, encompassing a cross-sectional field of approximately 20 cells. There is no evidence that these cells acquire distinctive identities due to subtle changes in the nuclear concentrations of the Dorsal protein. Rather, a variety of evidence suggests that the Dorsal gradient generates just three primary thresholds of gene activity. High levels activate gene expression in the presumptive mesoderm, while intermediate and low levels activate gene expression in the ventral and dorsal neurogenic ectoderm, respectively. We discuss how these primary readouts of the gradient establish localized domains of cell signalling, which work in a combinatorial manner with transcriptional networks to produce complex patterns of gene expression and tissue differentiation.

Transcriptional landscape estimation from tiling array data using a model of signal shift and drift

Motivation: High-density oligonucleotide tiling array technology holds the promise of a better description of the complexity and the dynamics of transcriptional landscapes. In organisms such as bacteria and yeasts, transcription can be measured on a genome-wide scale with a resolution >25 bp. The statistical models currently used to handle these data remain however very simple, the most popular being the piecewise constant Gaussian model with a fixed number of breakpoints.

Results: This article describes a new methodology based on a hidden Markov model that embeds the segmentation of a continuous-valued signal in a probabilistic setting. For a computationally affordable cost, this framework (i) alleviates the difficulty of choosing a fixed number of breakpoints, and (ii) permits retrieving more information than a unique segmentation by giving access to the whole probability distribution of the transcription profile. Importantly, the model is also enriched and accounts for subtle effects such as signal ‘drift’ and covariates. Relevance of this framework is demonstrated on a Bacillus subtilis dataset.

Availability: A software is distributed under the GPL.

Contact:pierre.nicolas@jouy.inra.fr

Supplementary information:Supplementary data is available at Bioinformatics online.

A distinct class of small RNAs arises from pre-miRNA-proximal regions in a simple chordate

MicroRNAs (miRNAs) have been implicated in various cellular processes. They are thought to function primarily as inhibitors of gene activity by attenuating translation or promoting mRNA degradation. A typical miRNA gene produces a predominant approximately 21-nucleotide (nt) RNA (the miRNA) along with a less abundant miRNA(*) product. We sought to identify miRNAs from the simple chordate Ciona intestinalis through comprehensive sequencing of small RNA libraries created from different developmental stages. Unexpectedly, half of the identified miRNA loci encode up to four distinct, stable small RNAs. The additional RNAs, miRNA-offset RNAs (moRs), are generated from sequences immediately adjacent to the predicted approximately 60-nt pre-miRNA. moRs seem to be produced by RNAse III-like processing, are approximately 20 nt long and, like miRNAs, are observed at specific developmental stages. We present evidence suggesting that the biogenesis of moRs results from an intrinsic property of the miRNA processing machinery in C. intestinalis.

How the Dorsal gradient works: insights from postgenome technologies

Gradients of extracellular signaling molecules and transcription factors are used in a variety of developmental processes, including the patterning of the Drosophila embryo, the establishment of diverse neuronal cell types in the vertebrate neural tube, and the anterior-posterior patterning of vertebrate limbs. Here, we discuss how a gradient of the maternal transcription factor Dorsal produces complex patterns of gene expression across the dorsal-ventral (DV) axis of the early Drosophila embryo. The identification of 60-70 Dorsal target genes, along with the characterization of approximately 35 associated regulatory DNAs, suggests that there are at least six different regulatory codes driving diverse DV expression profiles.

Toll-like receptor 4 restricts retinal progenitor cell proliferation

Retinal neurogenesis ceases by the early postnatal period, although retinal progenitor cells (RPCs) persist throughout life. In this study, we show that in the mammalian eye, the function of Toll-like receptor 4 (TLR4) extends beyond regulation of the innate immune response; it restricts RPC proliferation. In TLR4-deficient mice, enhanced proliferation of cells reminiscent of RPCs is evident during the early postnatal period. In vitro experiments demonstrate that TLR4 acts as an intrinsic regulator of RPC fate decision. Increased TLR4 expression in the eye correlates with the postnatal cessation of cell proliferation. However, deficient TLR4 expression is not sufficient to extend the proliferative period but rather contributes to resumption of proliferation in combination with growth factors. Proliferation in vivo is inhibited by both MyD88-dependent and -independent pathways, similar to the mechanisms activated by TLR4 in immune cells. Thus, our study attributes a novel role to TLR4 as a negative regulator of RPC proliferation.

The BMP-binding protein Crossveinless 2 is a short-range, concentration-dependent, biphasic modulator of BMP signaling in Drosophila

In Drosophila, the secreted BMP-binding protein Short gastrulation (Sog) inhibits signaling by sequestering BMPs from receptors, but enhances signaling by transporting BMPs through tissues. We show that Crossveinless 2 (Cv-2) is also a secreted BMP-binding protein that enhances or inhibits BMP signaling. Unlike Sog, however, Cv-2 does not promote signaling by transporting BMPs. Rather, Cv-2 binds cell surfaces and heparan sulfate proteoglygans and acts over a short range. Cv-2 binds the type I BMP receptor Thickveins (Tkv), and we demonstrate how the exchange of BMPs between Cv-2 and receptor can produce the observed biphasic response to Cv-2 concentration, where low levels promote and high levels inhibit signaling. Importantly, we show also how the concentration or type of BMP present can determine whether Cv-2 promotes or inhibits signaling. We also find that Cv-2 expression is controlled by BMP signaling, and these combined properties enable Cv-2 to exquisitely tune BMP signaling.

Promoter elements associated with RNA Pol II stalling in the Drosophila embryo

RNA Polymerase II (Pol II) is bound to the promoter regions of many or most developmental control genes before their activation during Drosophila embryogenesis. It has been suggested that Pol II stalling is used to produce dynamic and rapid responses of developmental patterning genes to transient cues such as extracellular signaling molecules. Here, we present a combined computational and experimental analysis of stalled promoters to determine how they come to bind Pol II in the early Drosophila embryo. At least one-fourth of the stalled promoters contain a shared sequence motif, the "pause button" (PB): KCGRWCG. The PB motif is sometimes located in the position of the DPE, and over one-fifth of the stalled promoters contain the following arrangement of core elements: GAGA, Inr, PB, and/or DPE. This arrangement was used to identify additional stalled promoters in the Drosophila genome, and permanganate footprint assays were used to confirm that the segmentation gene engrailed contains paused Pol II as seen for heat-shock genes. We discuss different models for Pol II binding and gene activation in the early embryo.

Self-regulatory circuits in dorsoventral axis formation of the short-germ beetle Tribolium castaneum

The rel/NF-kappaB transcription factor Dorsal controls dorsoventral (DV) axis formation in Drosophila. A stable nuclear gradient of Dorsal directly regulates approximately 50 target genes. In Tribolium castaneum (Tc), a beetle with an ancestral type of embryogenesis, the Dorsal nuclear gradient is not stable, but rapidly shrinks and disappears. We find that negative feedback accounts for this dynamic behavior: Tc-Dorsal and one of its target genes activate transcription of the IkB homolog Tc-cactus, terminating Dorsal function. Despite its transient role, Tc-Dorsal is strictly required to initiate DV polarity, as in Drosophila. However, unlike in Drosophila, embryos lacking Tc-Dorsal display a periodic pattern of DV cell fates along the AP axis, indicating that a self-organizing ectodermal patterning system operates independently of mesoderm or maternal DV polarity cues. Our results also elucidate how extraembryonic tissues are organized in short-germ embryos, and how patterning information is transmitted from the early embryo to the growth zone.

Wnt signalling and its impact on development and cancer

The Wnt signalling pathway is an ancient system that has been highly conserved during evolution. It has a crucial role in the embryonic development of all animal species, in the regeneration of tissues in adult organisms and in many other processes. Mutations or deregulated expression of components of the Wnt pathway can induce disease, most importantly cancer. The first gene to be identified that encodes a Wnt signalling component, Int1 (integration 1), was molecularly characterized from mouse tumour cells 25 years ago. In parallel, the homologous gene Wingless in Drosophila melanogaster, which produces developmental defects in embryos, was characterized. Since then, further components of the Wnt pathway have been identified and their epistatic relationships have been defined. This article is a Timeline of crucial discoveries about the components and functions of this essential pathway.

RNA polymerase stalling at developmental control genes in the Drosophila melanogaster embryo

It is widely assumed that the key rate-limiting step in gene activation is the recruitment of RNA polymerase II (Pol II) to the core promoter. Although there are well-documented examples in which Pol II is recruited to a gene but stalls, a general role for Pol II stalling in development has not been established. We have carried out comprehensive Pol II chromatin immunoprecipitation microarray (ChIP-chip) assays in Drosophila embryos and identified three distinct Pol II binding behaviors: active (uniform binding across the entire transcription unit), no binding, and stalled (binding at the transcription start site). The notable feature of the approximately 10% genes that are stalled is that they are highly enriched for developmental control genes, which are either repressed or poised for activation during later stages of embryogenesis. We propose that Pol II stalling facilitates rapid temporal and spatial changes in gene activity during development.

Functional conservation of Rel binding sites in drosophilid genomes

Evolutionary constraints on gene regulatory elements are poorly understood: Little is known about how the strength of transcription factor binding correlates with DNA sequence conservation, and whether transcription factor binding sites can evolve rapidly while retaining their function. Here we use the model of the NFKB/Rel-dependent gene regulation in divergent Drosophila species to examine the hypothesis that the functional properties of authentic transcription factor binding sites are under stronger evolutionary constraints than the genomic background. Using molecular modeling we compare tertiary structures of the Drosophila Rel family proteins Dorsal, Dif, and Relish and demonstrate that their DNA-binding and protein dimerization domains undergo distinct rates of evolution. The accumulated amino acid changes, however, are unlikely to affect DNA sequence recognition and affinity. We employ our recently developed microarray-based experimental platform and principal coordinates statistical analysis to quantitatively and systematically profile DNA binding affinities of three Drosophila Rel proteins to 10,368 variants of the NFKB recognition sequences. We then correlate the evolutionary divergence of gene regulatory regions with differences in DNA binding affinities. Genome-wide analyses reveal a significant increase in the number of conserved Rel binding sites in promoters of developmental and immune genes. Significantly, the affinity of Rel proteins to these sites was higher than to less conserved sites and was maintained by the conservation of the DNA binding site sequence (static conservation) or in some cases despite significantly diverged sequences (dynamic conservation). We discuss how two types of conservation may contribute to the stabilization and optimization of a functional gene regulatory code in evolution.

Whole-genome ChIP-chip analysis of Dorsal, Twist, and Snail suggests integration of diverse patterning processes in the Drosophila embryo

Genetic studies have identified numerous sequence-specific transcription factors that control development, yet little is known about their in vivo distribution across animal genomes. We determined the genome-wide occupancy of the dorsoventral (DV) determinants Dorsal, Twist, and Snail in the Drosophila embryo using chromatin immunoprecipitation coupled with microarray analysis (ChIP-chip). The in vivo binding of these proteins correlate tightly with the limits of known enhancers. Our analysis predicts substantially more target genes than previous estimates, and includes Dpp signaling components and anteroposterior (AP) segmentation determinants. Thus, the ChIP-chip data uncover a much larger than expected regulatory network, which integrates diverse patterning processes during development.

Threshold response of C15 to the Dpp gradient in Drosophila is established by the cumulative effect of Smad and Zen activators and negative cues

Morphogen gradients determine a range of cell fates by specifying multiple transcriptional threshold responses. In the dorsal ectoderm of the Drosophila embryo, a BMP gradient is translated into an activated Smad transcription factor gradient, which elicits at least three threshold responses - high, intermediate and low. However, the mechanism underlying differential response to Dpp is poorly understood, due in part to the insufficient number of well-studied target genes. We analyzed the regulation of the C15 gene, which can be activated in cells containing intermediate levels of Dpp. We show that C15 expression requires both dpp and zen, thus forming a genetic feed-forward loop. The C15 regulatory element contains clusters of Smad- and Zen-binding sites in close proximity. Mutational analysis shows that the number of intact Smad- and Zen-binding sites is essential for the C15 transcriptional response, and that the spatial limits of C15 expression are established through a repression mechanism in the dorsolateral cells of the embryo. Thus, the combinatorial action of Smad and Zen activators bound to a number of adjacent sites, and competing negative cues allows for proper gene response to lower than peak levels of the Dpp morphogen.

The Transcriptome of the Sea Urchin Embryo

The sea urchin Strongylocentrotus purpuratus is a model organism for study of the genomic control circuitry underlying embryonic development. We examined the complete repertoire of genes expressed in the S. purpuratus embryo, up to late gastrula stage, by means of high-resolution custom tiling arrays covering the whole genome. We detected complete spliced structures even for genes known to be expressed at low levels in only a few cells. At least 11,000 to 12,000 genes are used in embryogenesis. These include most of the genes encoding transcription factors and signaling proteins, as well as some classes of general cytoskeletal and metabolic proteins, but only a minor fraction of genes encoding immune functions and sensory receptors. Thousands of small asymmetric transcripts of unknown function were also detected in intergenic regions throughout the genome. The tiling array data were used to correct and authenticate several thousand gene models during the genome annotation process.