During development, extracellular signaling molecules interact with intracellular gene networks to control the specification, pattern and size of organs. One such signaling molecule is Hedgehog (Hh). Hh is known to act as a morphogen, instructing different fates depending on the distance to its source. However, how Hh, when signaling across a cell field, impacts organ-specific transcriptional networks is still poorly understood. Here, we investigate this issue during the development of the Drosophila ocellar complex. The development of this sensory structure, which is composed of three simple eyes (or ocelli) located at the vertices of a triangular patch of cuticle on the dorsal head, depends on Hh signaling and on the definition of three domains: two areas of eya and so expression--the prospective anterior and posterior ocelli--and the intervening interocellar domain. Our results highlight the role of the homeodomain transcription factor engrailed (en) both as a target and as a transcriptional repressor of hh signaling in the prospective interocellar region. Furthermore, we identify a requirement for the Notch pathway in the establishment of en maintenance in a Hh-independent manner. Therefore, hh signals transiently during the specification of the interocellar domain, with en being required here for hh signaling attenuation. Computational analysis further suggests that this network design confers robustness to signaling noise and constrains phenotypic variation. In summary, using genetics and modeling we have expanded the ocellar gene network to explain how the interaction between the Hh gradient and this gene network results in the generation of stable mutually exclusive gene expression domains. In addition, we discuss some general implications our model may have in some Hh-driven gene networks.

During organ development, the progenitor state is transient, and depends on specific combinations of transcription factors and extracellular signals. Not surprisingly, abnormal maintenance of progenitor transcription factors may lead to tissue overgrowth, and the concurrence of signals from the local environment is often critical to trigger this overgrowth. Therefore, identifying specific combinations of transcription factors/signals promoting -or opposing- proliferation in progenitors is essential to understand normal development and disease. We have investigated this issue using the Drosophila eye as model. Transcription factors hth and tsh are transiently expressed in eye progenitors causing the expansion of the progenitor pool. However, if their co-expression is maintained experimentally, cell proliferation continues and differentiation is halted. Here we show that Hth+Tsh-induced tissue overgrowth requires the BMP2 Dpp and the abnormal hyperactivation of its pathway. Rather than using autocrine Dpp expression, Hth+Tsh cells increase their avidity for Dpp, produced locally, by upregulating extracellular matrix components. During normal development, Dpp represses hth and tsh ensuring that the progenitor state is transient. However, cells in which Hth+Tsh expression is forcibly maintained use Dpp to enhance their proliferation.

The morphology and function of organs depend on coordinated changes in gene expression during development. These changes are controlled by transcription factors, signaling pathways, and their regulatory interactions, which are represented by gene regulatory networks (GRNs). Therefore, the structure of an organ GRN restricts the morphological and functional variations that the organ can experience-its potential morphospace. Therefore, two important questions arise when studying any GRN: what is the predicted available morphospace and what are the regulatory linkages that contribute the most to control morphological variation within this space. Here, we explore these questions by analyzing a small "three-node" GRN model that captures the Hh-driven regulatory interactions controlling a simple visual structure: the ocellar region of Drosophila. Analysis of the model predicts that random variation of model parameters results in a specific non-random distribution of morphological variants. Study of a limited sample of drosophilids and other dipterans finds a correspondence between the predicted phenotypic range and that found in nature. As an alternative to simulations, we apply Bayesian networks methods in order to identify the set of parameters with the largest contribution to morphological variation. Our results predict the potential morphological space of the ocellar complex and identify likely candidate processes to be responsible for ocellar morphological evolution using Bayesian networks. We further discuss the assumptions that the approach we have taken entails and their validity.

Odd-skipped family of proteins (Odd in Drosophila and Osr in vertebrates) are evolutionarily conserved zinc finger transcription factors. Two Osr genes are present in mammalian genomes, and it was recently reported that Osr1, but not Osr2, is required for murine kidney development. Here, we show that in Xenopus and zebrafish both Osr1 and Osr2 are necessary and sufficient for the development of the pronephros. Osr genes are expressed in early prospective pronephric territories, and morphants for either of the two genes show severely impaired kidney development. Conversely, overexpression of Osr genes promotes formation of ectopic kidney tissue. Molecularly, Osr proteins function as transcriptional repressors during kidney formation. We also show that Drosophila Odd induces kidney tissue in Xenopus. This might be accomplished through recruitment of Groucho-like co-repressors. Odd genes may also be required for proper development of the Malpighian tubules, the Drosophila renal organs. Our results highlight the evolutionary conserved involvement of Odd-skipped transcription factors in the development of kidneys.

Signalling by TGFβ superfamily factors plays an important role in tissue growth and cell proliferation. In Drosophila, the activity of the TGFβ/Activin signalling branch has been linked to the regulation of cell growth and proliferation, but the cellular and molecular basis for these functions are not fully understood. In this study, we show that both the RII receptor Punt (Put) and the R-Smad Smad2 are strongly required for cell and tissue growth. Knocking down the expression of Put or Smad2 in salivary glands causes alterations in nucleolar structure and functions. Cells with decreased TGFβ/Activin signalling accumulate intermediate pre-rRNA transcripts containing internal transcribed spacer 1 regions accompanied by the nucleolar retention of ribosomal proteins. Thus, our results show that TGFβ/Activin signalling is required for ribosomal biogenesis, a key aspect of cellular growth control. Importantly, overexpression of Put enhanced cell growth induced by Drosophila Myc, a well-characterized inducer of nucleolar hypertrophy and ribosome biogenesis.

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Alternative splicing increases neuronal transcriptomic complexity throughout animal phylogeny. To delve into the mechanisms controlling the assembly and evolution of this regulatory layer, we characterized the neuronal microexon program in Drosophila and compared it with that of mammals. In nonvertebrate bilaterians, this splicing program is restricted to neurons by the posttranscriptional processing of the enhancer of microexons (eMIC) domain in Srrm234. In Drosophila, this processing is dependent on regulation by Elav/Fne. eMIC deficiency or misexpression leads to widespread neurological alterations largely emerging from impaired neuronal activity, as revealed by a combination of neuronal imaging experiments and cell type-specific rescues. These defects are associated with the genome-wide skipping of short neural exons, which are strongly enriched in ion channels. We found no overlap of eMIC-regulated exons between flies and mice, illustrating how ancient posttranscriptional programs can evolve independently in different phyla to affect distinct cellular modules while maintaining cell-type specificity.

The homeodomain transcription factor Prrxl1/DRG11 has emerged as a crucial molecule in the establishment of the pain circuitry, in particular spinal cord targeting of dorsal root ganglia (DRG) axons and differentiation of nociceptive glutamatergic spinal cord neurons. Despite Prrxl1 importance in the establishment of the DRG-spinal nociceptive circuit, the molecular mechanisms that regulate its expression along development remain largely unknown. Here, we show that Prrxl1 transcription is regulated by three alternative promoters (named P1, P2, and P3), which control the expression of three distinct Prrxl1 5'-UTR variants, named 5'-UTR-A, 5'-UTR-B, and 5'-UTR-C. These 5'-UTR sequences confer distinct mRNA stability and translation efficiency to the Prrxl1 transcript. The most conserved promoter (P3) contains a TATA-box and displays in vivo enhancer activity in a pattern that overlaps with the zebrafish Prrxl1 homologue, drgx. Regulatory modules present in this sequence were identified and characterized, including a binding site for Phox2b. Concomitantly, we demonstrate that zebrafish Phox2b is required for the expression of drgx in the facial, glossopharyngeal, and vagal cranial ganglia.

Planarians regenerate a whole animal from a small body piece within a few days. Recent studies have shown that the bone morphogenetic protein (BMP) pathway is required to reestablish the dorsoventral (DV) axis. In vertebrates, the specification of the DV axis depends on the coordinated action of a dual organizer defined by BMP and antidorsalizing morphogenetic protein (ADMP) under the control of several factors, including the inhibitors chordin and noggin. Planarians have an expanded noggin family (up to ten members), which have been classified as canonical noggin (nog) and noggin-like (nlg) genes, the latter carrying an insertion within the noggin domain. Here we show that a BMP/ADMP organizer governs DV axis reestablishment during planarian regeneration, highlighting a greater-than-thought conservation of the mechanisms that establish this axis in protostomes and deuterostomes. Also, we report that whereas noggin genes function as canonical BMP inhibitors, the silencing of planarian nlg8 induces ectopic neurogenesis and enhances ventralizing bmp(RNAi) phenotypes. Finally, we show that noggin-like genes are conserved from cnidarian to vertebrates and that both planarian nlg8 and Xenopus nlg ventralize Xenopus embryos when overexpressed. Remarkably, this ventralization is not associated with an increase in SMAD1/5/8 phosphorylation.

During Drosophila eye development, recruitment of retinal precursors from a pool of progenitor cells is tightly coupled to proliferation control. However, how this coupling operates is still unclear. Here we show that the transcription factor hth, together with eyeless, is required to stimulate proliferation of progenitor cells. Accordingly, knocking down hth expression results in severely reduced eyes. Our experiments reveal three additional functions for hth: the cell cycle of progenitors is characterized by a relatively long G2 phase, which makes them prone to enter mitosis; hth represses the burst of string/cdc25 expression that precedes G1 arrest, and also the early expression of the proneural gene atonal. Thereby, hth maintains the proliferative and undifferentiated state of eye progenitors. Furthermore, we show that the G1 synchronization that characterizes retinal precursors is the result of the spatially controlled repression of hth by Dpp and Hh, and not of an actively induced cell cycle arrest. We integrate these results in a model of the early steps of eye development that links proliferation control and differential gene expression with patterning signals.

Meis1, a conserved transcription factor of the TALE-homeodomain class, is expressed in a wide variety of tissues during development. Its complex expression pattern is likely to be controlled by an equally complex regulatory landscape. Here we have scanned the Meis1 locus for regulatory elements and found 13 non-coding regions, highly conserved between humans and teleost fishes, that have enhancer activity in stable transgenic zebrafish lines. All these regions are syntenic in most vertebrates. The composite expression of all these enhancer elements recapitulate most of Meis1 expression during early embryogenesis, indicating they comprise a basic set of regulatory elements of the Meis1 gene. Using bioinformatic tools, we identify a number of potential binding sites for transcription factors that are compatible with the regulation of these enhancers. Specifically, HHc2:066650, which is expressed in the developing retina and optic tectum, harbors several predicted Pax6 sites. Biochemical, functional and transgenic assays indicate that pax6 genes directly regulate HHc2:066650 activity.

The evolution of winged insects revolutionized terrestrial ecosystems and led to the largest animal radiation on Earth. However, we still have an incomplete picture of the genomic changes that underlay this diversification. Mayflies, as one of the sister groups of all other winged insects, are key to understanding this radiation. Here, we describe the genome of the mayfly Cloeon dipterum and its gene expression throughout its aquatic and aerial life cycle and specific organs. We discover an expansion of odorant-binding-protein genes, some expressed specifically in breathing gills of aquatic nymphs, suggesting a novel sensory role for this organ. In contrast, flying adults use an enlarged opsin set in a sexually dimorphic manner, with some expressed only in males. Finally, we identify a set of wing-associated genes deeply conserved in the pterygote insects and find transcriptomic similarities between gills and wings, suggesting a common genetic program. Globally, this comprehensive genomic and transcriptomic study uncovers the genetic basis of key evolutionary adaptations in mayflies and winged insects.

ETCHbox genes are eutherian-specific homeobox genes expressed during preimplantation development at a time when the first cell lineage decisions are being made. The mouse has an unusual repertoire of ETCHbox genes with several gene families lost in evolution and the remaining two, Crxos and Obox, greatly divergent in sequence and number. Each has undergone duplication to give a double homeodomain Crxos locus and a large cluster of over 60 Obox loci. The gene content differences between species raise important questions about how evolution can tolerate loss of genes implicated in key developmental events.

Organ-selector transcription factors control simultaneously cell differentiation and proliferation, ensuring the development of functional organs and their homeostasis. How this is achieved at the molecular level is still unclear. Here we have investigated how the transcriptional pulse of string/cdc25 (stg), the universal mitotic trigger, is regulated during Drosophila retina development as an example of coordinated deployment of differentiation and proliferation programs. We identify the eye specific stg enhancer, stg-FMW, and show that Pax6 selector genes, in cooperation with Eya and So, two members of the retinal determination network, activate stg-FMW, establishing a positive feed-forward loop. This loop is negatively modulated by the Meis1 protein, Hth. This regulatory logic is reminiscent of that controlling the expression of differentiation transcription factors. Our work shows that subjecting transcription factors and key cell cycle regulators to the same regulatory logic ensures the coupling between differentiation and proliferation programs during organ development.

The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians.

Non-coding DNA conservation across species has been often used as a predictor for transcriptional enhancer activity. However, only a few systematic analyses of the function of these highly conserved non-coding regions (HCNRs) have been performed. Here we use zebrafish transgenic assays to perform a systematic study of 113 HCNRs from human chromosome 16. By comparing transient and stable transgenesis, we show that the first method is highly inefficient, leading to 40% of false positives and 20% of false negatives. When analyzed in stable transgenic lines, a great majority of HCNRs were active in the central nervous system, although some of them drove expression in other organs such as the eye and the excretory system. Finally, by testing a fraction of the HCNRs lacking enhancer activity for in vivo insulator activity, we find that 20% of them may contain enhancer-blocking function. Altogether our data indicate that HCNRs may contain different types of cis-regulatory activity, including enhancer, insulators as well as other not yet discovered functions.

Genome-wide association studies (GWAS) identified the MEIS1 locus for Restless Legs Syndrome (RLS), but causal single nucleotide polymorphisms (SNPs) and their functional relevance remain unknown. This locus contains a large number of highly conserved noncoding regions (HCNRs) potentially functioning as cis-regulatory modules. We analyzed these HCNRs for allele-dependent enhancer activity in zebrafish and mice and found that the risk allele of the lead SNP rs12469063 reduces enhancer activity in the Meis1 expression domain of the murine embryonic ganglionic eminences (GE). CREB1 binds this enhancer and rs12469063 affects its binding in vitro. In addition, MEIS1 target genes suggest a role in the specification of neuronal progenitors in the GE, and heterozygous Meis1-deficient mice exhibit hyperactivity, resembling the RLS phenotype. Thus, in vivo and in vitro analysis of a common SNP with small effect size showed allele-dependent function in the prospective basal ganglia representing the first neurodevelopmental region implicated in RLS.

Proper organ patterning depends on a tight coordination between cell proliferation and differentiation. The patterning of Drosophila retina occurs both very fast and with high precision. This process is driven by the dynamic changes in signaling activity of the conserved Hedgehog (Hh) pathway, which coordinates cell fate determination, cell cycle and tissue morphogenesis. Here we show that during Drosophila retinogenesis, the retinal determination gene dachshund (dac) is not only a target of the Hh signaling pathway, but is also a modulator of its activity. Using developmental genetics techniques, we demonstrate that dac enhances Hh signaling by promoting the accumulation of the Gli transcription factor Cubitus interruptus (Ci) parallel to or downstream of fused. In the absence of dac, all Hh-mediated events associated to the morphogenetic furrow are delayed. One of the consequences is that, posterior to the furrow, dac- cells cannot activate a Roadkill-Cullin3 negative feedback loop that attenuates Hh signaling and which is necessary for retinal cells to continue normal differentiation. Therefore, dac is part of an essential positive feedback loop in the Hh pathway, guaranteeing the speed and the accuracy of Drosophila retinogenesis.

The functional consequences of whole genome duplications in vertebrate evolution are not fully understood. It remains unclear, for instance, why paralogues were retained in some gene families but extensively lost in others. Cdx homeobox genes encode conserved transcription factors controlling posterior development across diverse bilaterians. These genes are part of the ParaHox gene cluster. Multiple Cdx copies were retained after genome duplication, raising questions about how functional divergence, overlap, and redundancy respectively contributed to their retention and evolutionary fate.

All vertebrate brains develop following a common Bauplan defined by anteroposterior (AP) and dorsoventral (DV) subdivisions, characterized by largely conserved differential expression of gene markers. However, it is still unclear how this Bauplan originated during evolution. We studied the relative expression of 48 genes with key roles in vertebrate neural patterning in a representative amphioxus embryonic stage. Unlike nonchordates, amphioxus develops its central nervous system (CNS) from a neural plate that is homologous to that of vertebrates, allowing direct topological comparisons. The resulting genoarchitectonic model revealed that the amphioxus incipient neural tube is unexpectedly complex, consisting of several AP and DV molecular partitions. Strikingly, comparison with vertebrates indicates that the vertebrate thalamus, pretectum, and midbrain domains jointly correspond to a single amphioxus region, which we termed Di-Mesencephalic primordium (DiMes). This suggests that these domains have a common developmental and evolutionary origin, as supported by functional experiments manipulating secondary organizers in zebrafish and mice.