The heads of annelids (earthworms, polychaetes, and others) and arthropods (insects, myriapods, spiders, and others) and the arthropod-related onychophorans (velvet worms) show similar brain architecture and for this reason have long been considered homologous. However, this view is challenged by the 'new phylogeny' placing arthropods and annelids into distinct superphyla, Ecdysozoa and Lophotrochozoa, together with many other phyla lacking elaborate heads or brains. To compare the organisation of annelid and arthropod heads and brains at the molecular level, we investigated head regionalisation genes in various groups. Regionalisation genes subdivide developing animals into molecular regions and can be used to align head regions between remote animal phyla.

RNA-seq based on short reads generated by next generation sequencing technologies has become the main approach to study differential gene expression. Until now, the main applications of this technique have been to study the variation of gene expression in a whole organism, tissue or cell type under different conditions or at different developmental stages. However, RNA-seq also has a great potential to be used in evolutionary studies to investigate gene expression divergence in closely related species.

A striking diversity of compound eye size and shape has evolved among insects. The number of ommatidia and their size are major determinants of the visual sensitivity and acuity of the compound eye. Each ommatidium is composed of eight photoreceptor cells that facilitate the discrimination of different colours via the expression of various light sensitive Rhodopsin proteins. It follows that variation in eye size, shape, and opsin composition is likely to directly influence vision. We analyzed variation in these three traits in D. melanogaster, D. simulans and D. mauritiana. We show that D. mauritiana generally has larger eyes than its sibling species, which is due to a combination of larger ommatidia and more ommatidia. In addition, intra- and inter-specific differences in eye size among D. simulans and D. melanogaster strains are mainly caused by variation in ommatidia number. By applying a geometric morphometrics approach to assess whether the formation of larger eyes influences other parts of the head capsule, we found that an increase in eye size is associated with a reduction in the adjacent face cuticle. Our shape analysis also demonstrates that D. mauritiana eyes are specifically enlarged in the dorsal region. Intriguingly, this dorsal enlargement is associated with enhanced expression of rhodopsin 3 in D. mauritiana. In summary, our data suggests that the morphology and functional properties of the compound eyes vary considerably within and among these closely related Drosophila species and may be part of coordinated morphological changes affecting the head capsule.

Several highly conserved genes play a role in anterior neural plate patterning of vertebrates and in head and brain patterning of insects. However, head involution in Drosophila has impeded a systematic identification of genes required for insect head formation. Therefore, we use the red flour beetle Tribolium castaneum in order to comprehensively test the function of orthologs of vertebrate neural plate patterning genes for a function in insect head development. RNAi analysis reveals that most of these genes are indeed required for insect head capsule patterning, and we also identified several genes that had not been implicated in this process before. Furthermore, we show that Tc-six3/optix acts upstream of Tc-wingless, Tc-orthodenticle1, and Tc-eyeless to control anterior median development. Finally, we demonstrate that Tc-six3/optix is the first gene known to be required for the embryonic formation of the central complex, a midline-spanning brain part connected to the neuroendocrine pars intercerebralis. These functions are very likely conserved among bilaterians since vertebrate six3 is required for neuroendocrine and median brain development with certain mutations leading to holoprosencephaly.

The head gap genes orthodenticle (otd), empty spiracles (ems) and buttonhead (btd) are required for metamerization and segment specification in Drosophila. We asked whether the function of their orthologs is conserved in the red flour beetle Tribolium castaneum which in contrast to Drosophila develops its larval head in a way typical for insects. We find that depending on dsRNA injection time, two functions of Tc-orthodenticle1 (Tc-otd1) can be identified. The early regionalization function affects all segments formed during the blastoderm stage while the later head patterning function is similar to Drosophila. In contrast, both expression and function of Tc-empty spiracles (Tc-ems) are restricted to the posterior part of the ocular and the anterior part of the antennal segment and Tc-buttonhead (Tc-btd) is not required for head cuticle formation at all. We conclude that the gap gene like roles of ems and btd are not conserved while at least the head patterning function of otd appears to be similar in fly and beetle. Hence, the ancestral mode of insect head segmentation remains to be discovered. With this work, we establish Tribolium as a model system for arthropod head development that does not suffer from the Drosophila specific problems like head involution and strongly reduced head structures.

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.

In Drosophila, large variations in rearrangement rate have been reported among different lineages and among Muller's elements. Nevertheless, the mechanisms that are involved in the generation of inversions, their increase in frequency, as well as their impact on the genome are not completely understood. This is in part due to the lack of comparative studies on species distantly related to Drosophila melanogaster. Therefore, we sequenced and assembled the genomes of two species of the virilis phylad (Drosophila novamexicana [15010-1031.00] and Drosophila americana [SF12]), which are diverging from D. melanogaster for more than 40 Myr. Based on these data, we identified the precise location of six novel inversion breakpoints. A molecular characterization provided clear evidence that DAIBAM (a miniature inverted-repeat transposable element) was involved in the generation of eight out of the nine inversions identified. In contrast to what has been previously reported for D. melanogaster and close relatives, ectopic recombination is thus the prevalent mechanism of generating inversions in species of the virilis phylad. Using pool-sequencing data for three populations of D. americana, we also show that common polymorphic inversions create a high degree of genetic differentiation between populations for chromosomes X, 4, and 5 over large physical distances. We did not find statistically significant differences in expression levels between D. americana (SF12) and D. novamexicana (15010-1031.00) strains for the three genes surveyed (CG9588, Fig 4, and fab1) flanking three inversion breakpoints.

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.

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.

The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality.

The size and shape of organs is tightly controlled to achieve optimal function. Natural morphological variations often represent functional adaptations to an ever-changing environment. For instance, variation in head morphology is pervasive in insects and the underlying molecular basis is starting to be revealed in the Drosophila genus for species of the melanogaster group. However, it remains unclear whether similar diversifications are governed by similar or different molecular mechanisms over longer timescales. To address this issue, we used species of the virilis phylad because they have been diverging from D. melanogaster for at least 40 million years. Our comprehensive morphological survey revealed remarkable differences in eye size and head shape among these species with D. novamexicana having the smallest eyes and southern D. americana populations having the largest eyes. We show that the genetic architecture underlying eye size variation is complex with multiple associated genetic variants located on most chromosomes. Our genome wide association study (GWAS) strongly suggests that some of the putative causative variants are associated with the presence of inversions. Indeed, northern populations of D. americana share derived inversions with D. novamexicana and they show smaller eyes compared to southern ones. Intriguingly, we observed a significant enrichment of genes involved in eye development on the 4th chromosome after intersecting chromosomal regions associated with phenotypic differences with those showing high differentiation among D. americana populations. We propose that variants associated with chromosomal inversions contribute to both intra- and interspecific variation in eye size among species of the virilis phylad.

Phenotypes on-demand generated by controlling activation and accumulation of proteins of interest are invaluable tools to analyse and engineer biological processes. While temperature-sensitive alleles are frequently used as conditional mutants in microorganisms, they are usually difficult to identify in multicellular species. Here we present a versatile and transferable, genetically stable system based on a low-temperature-controlled N-terminal degradation signal (lt-degron) that allows reversible and switch-like tuning of protein levels under physiological conditions in vivo. Thereby, developmental effects can be triggered and phenotypes on demand generated. The lt-degron was established to produce conditional and cell-type-specific phenotypes and is generally applicable in a wide range of organisms, from eukaryotic microorganisms to plants and poikilothermic animals. We have successfully applied this system to control the abundance and function of transcription factors and different enzymes by tunable protein accumulation.

The Onychophora are a probable sister group to Arthropoda, one of the most intensively studied animal phyla from a developmental perspective. Pioneering work on the fruit fly Drosophila melanogaster and subsequent investigation of other arthropods has revealed important roles for Wnt genes during many developmental processes in these animals.

Drosophila melanogaster head development represents a valuable process to study the developmental control of various organs, such as the antennae, the dorsal ocelli and the compound eyes from a common precursor, the eye-antennal imaginal disc. While the gene regulatory network underlying compound eye development has been extensively studied, the key transcription factors regulating the formation of other head structures from the same imaginal disc are largely unknown. We obtained the developmental transcriptome of the eye-antennal discs covering late patterning processes at the late 2nd larval instar stage to the onset and progression of differentiation at the end of larval development. We revealed the expression profiles of all genes expressed during eye-antennal disc development and we determined temporally co-expressed genes by hierarchical clustering. Since co-expressed genes may be regulated by common transcriptional regulators, we combined our transcriptome dataset with publicly available ChIP-seq data to identify central transcription factors that co-regulate genes during head development. Besides the identification of already known and well-described transcription factors, we show that the transcription factor Hunchback (Hb) regulates a significant number of genes that are expressed during late differentiation stages. We confirm that hb is expressed in two polyploid subperineurial glia cells (carpet cells) and a thorough functional analysis shows that loss of Hb function results in a loss of carpet cells in the eye-antennal disc. Additionally, we provide for the first time functional data indicating that carpet cells are an integral part of the blood-brain barrier. Eventually, we combined our expression data with a de novo Hb motif search to reveal stage specific putative target genes of which we find a significant number indeed expressed in carpet cells.

During the development of multicellular organisms the fate of individual cells is specified with great precision and reproducibility. Although classical genetic approaches led to the identification of many of the signaling pathways contributing to cell fate specification, they have provided little insight into the mechanisms that ensure robustness and reproducibility. We have used the specification of the R7 photoreceptor cells controlled by the Sevenless receptor tyrosine kinase (Sev) pathway to screen for modulators of pathway activity and to uncover the mechanisms underlying the robustness of cell fate decisions. Here we provide genetic evidence that the Drosophila SOCS36E adaptor protein containing an SH2 domain and a SOCS box acts as an attenuator of Sev signaling. Overexpression of Socs36E strongly suppresses the specification of extra R7 photoreceptor cells in response to constitutive activation of Sev, and loss of Socs36E function suppresses the loss of R7 cells when Sev activity is impaired. In a wild-type background, however, loss and gain of Socs36E function exhibits little effect on R7 specification. We also show that SH2 domain of SOCS36E is essential for this function in inhibiting Sev action and that Socs36E expression is suppressed by high Sev pathway activity. In our model, only the cell able to activate high levels of receptor tyrosine kinase signaling will repress SOCS36E expression, reduce the negative effect on Sev signaling and allow this cell to differentiate into R7. In contrast, the remaining cells fail to receive high signaling, and thus maintain high levels of SOCS36E. This represses residual Sev activity and blocks R7 development. Therefore, Socs36E constitutes a novel partially redundant feedback mechanism that contributes to the robustness of R7 specification. The SOCS family of adaptor proteins may have evolved as modulators of specific signaling pathways that contribute to the robustness and precision of cell fate specification.

Organizers play important roles during the embryonic development of many animals. The most famous example is the Spemann organizer that sets up embryonic axes in amphibian embryos. In spiders, a group of BMP secreting mesenchymal cells (the cumulus) functions as an organizer of the dorsoventral axis. Similar to experiments performed with the Spemann organizer, transplantation of the cumulus is able to induce a secondary axis in spiders. Despite the importance of this structure, it is unknown which factors are needed to activate cumulus specific gene expression. To address this question, we performed a transcriptomic analysis of early embryonic development in the spider Parasteatoda tepidariorum. Through this work, we found that the transcription factor Pt-Ets4 is needed for cumulus integrity, dorsoventral patterning and for the activation of Pt-hunchback and Pt-twist expression. Furthermore, ectopic expression of Pt-Ets4 is sufficient to induce cell delamination and migration by inducing a mesoderm-like cell fate.

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Two visual systems are present in most arthropod groups: median and lateral eyes. Most of our current knowledge about the developmental and molecular mechanisms involved in eye formation in arthropods comes from research in the model system Drosophila melanogaster. Here, a core set of retinal determination genes, namely, sine-oculis (so), eyes absent (eya), dachshund (dac), and the two pax6 orthologues eyeless (ey) and twin of eyeless (toy) govern early retinal development. By contrast, not much is known about the development of the up-to-eight eyes present in spiders. Therefore, we analyzed the embryonic expression of core retinal determination genes in the common house spider Parasteatoda tepidariorum.

Identifying the genetic mechanisms underlying phenotypic change is essential to understanding how gene regulatory networks and ultimately the genotype-to-phenotype map evolve. It is recognized that microRNAs (miRNAs) have the potential to facilitate evolutionary change [1-3]; however, there are no known examples of natural morphological variation caused by evolutionary changes in miRNA expression. Therefore, the contribution of miRNAs to evolutionary change remains unknown [1, 4]. Drosophila melanogaster subgroup species display a portion of trichome-free cuticle on the femur of the second leg called the "naked valley." It was previously shown that Ultrabithorax (Ubx) is involved in naked valley variation between D. melanogaster and D. simulans [5, 6]. However, naked valley size also varies among populations of D. melanogaster, ranging from 1,000 up to 30,000 μm(2). We investigated the genetic basis of intraspecific differences in the naked valley in D. melanogaster and found that neither Ubx nor shavenbaby (svb) [7, 8] contributes to this morphological difference. Instead, we show that changes in mir-92a expression underlie the evolution of naked valley size in D. melanogaster through repression of shavenoid (sha) [9]. Therefore, our results reveal a novel mechanism for morphological evolution and suggest that modulation of the expression of miRNAs potentially plays a prominent role in generating organismal diversity.