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.

Pre-mRNA splicing is a critical step of gene expression in eukaryotes. Transcriptome-wide splicing patterns are complex and primarily regulated by a diverse set of recognition elements and associated RNA-binding proteins. The retention and splicing (RES) complex is formed by three different proteins (Bud13p, Pml1p and Snu17p) and is involved in splicing in yeast. However, the importance of the RES complex for vertebrate splicing, the intronic features associated with its activity, and its role in development are unknown. In this study, we have generated loss-of-function mutants for the three components of the RES complex in zebrafish and showed that they are required during early development. The mutants showed a marked neural phenotype with increased cell death in the brain and a decrease in differentiated neurons. Transcriptomic analysis of bud13, snip1 (pml1) and rbmx2 (snu17) mutants revealed a global defect in intron splicing, with strong mis-splicing of a subset of introns. We found these RES-dependent introns were short, rich in GC and flanked by GC depleted exons, all of which are features associated with intron definition. Using these features, we developed and validated a predictive model that classifies RES dependent introns. Altogether, our study uncovers the essential role of the RES complex during vertebrate development and provides new insights into its function during splicing.

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.

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.

Recursive splicing (RS) starts by defining an "RS-exon," which is then spliced to the preceding exon, thus creating a recursive 5' splice site (RS-5ss). Previous studies focused on cryptic RS-exons, and now we find that the exon junction complex (EJC) represses RS of hundreds of annotated, mainly constitutive RS-exons. The core EJC factors, and the peripheral factors PNN and RNPS1, maintain RS-exon inclusion by repressing spliceosomal assembly on RS-5ss. The EJC also blocks 5ss located near exon-exon junctions, thus repressing inclusion of cryptic microexons. The prevalence of annotated RS-exons is high in deuterostomes, while the cryptic RS-exons are more prevalent in Drosophila, where EJC appears less capable of repressing RS. Notably, incomplete repression of RS also contributes to physiological alternative splicing of several human RS-exons. Finally, haploinsufficiency of the EJC factor Magoh in mice is associated with skipping of RS-exons in the brain, with relevance to the microcephaly phenotype and human diseases.

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.

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.

Epithelial-mesenchymal interactions are crucial for the development of numerous animal structures. Thus, unraveling how molecular tools are recruited in different lineages to control interplays between these tissues is key to understanding morphogenetic evolution. Here, we study Esrp genes, which regulate extensive splicing programs and are essential for mammalian organogenesis. We find that Esrp homologs have been independently recruited for the development of multiple structures across deuterostomes. Although Esrp is involved in a wide variety of ontogenetic processes, our results suggest ancient roles in non-neural ectoderm and regulating specific mesenchymal-to-epithelial transitions in deuterostome ancestors. However, consistent with the extensive rewiring of Esrp-dependent splicing programs between phyla, most developmental defects observed in vertebrate mutants are related to other types of morphogenetic processes. This is likely connected to the origin of an event in Fgfr, which was recruited as an Esrp target in stem chordates and subsequently co-opted into the development of many novel traits in vertebrates.

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.

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.

The specification and growth of organs is controlled simultaneously by networks of transcription factors. While the connection between these transcription factors with fate determinants is increasingly clear, how they establish the link with the cell cycle is far less understood. Here we investigate this link in the developing Drosophila eye, where two transcription factors, the MEIS1 homologue hth and the Zn-finger tsh, synergize to stimulate the proliferation of naïve eye progenitors. Experiments combining transcriptomics, open-chromatin profiling, motif analysis and functional assays indicate that these progenitor transcription factors exert a global regulation of the proliferation program. Rather than directly regulating cell cycle genes, they control proliferation through an intermediary layer of nuclear receptors of the ecdysone/estrogen-signaling pathway. This regulatory subnetwork between hth, tsh and nuclear receptors might be conserved from Drosophila to mammals, as we find a significant co-overexpression of their human homologues in specific cancer types.

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.

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.