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.

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 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.

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.

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.

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.

The re-use of genes in new organs forms the base of many evolutionary novelties. A well-characterised case is the recruitment of the posterior spiracle gene network to the Drosophila male genitalia. Here we find that this network has also been co-opted to the testis mesoderm where is required for sperm liberation, providing an example of sequentially repeated developmental co-options. Associated to this co-option event, an evolutionary expression novelty appeared, the activation of the posterior segment determinant Engrailed to the anterior A8 segment controlled by common testis and spiracle regulatory elements. Enhancer deletion shows that A8 anterior Engrailed activation is not required for spiracle development but only necessary in the testis. Our study presents an example of pre-adaptive developmental novelty: the activation of the Engrailed transcription factor in the anterior compartment of the A8 segment where, despite having no specific function, opens the possibility of this developmental factor acquiring one. We propose that recently co-opted networks become interlocked, so that any change to the network because of its function in one organ, will be mirrored by other organs even if it provides no selective advantage to them.

Insect compound eyes are composed of ommatidia, which contain photoreceptor cells that are sensitive to different wavelengths of light defined by the specific rhodopsin proteins that they express. The fruit fly Drosophila melanogaster has several different ommatidium types that can be localised to specific retinal regions, such as the dorsal rim area (DRA), or distributed stochastically in a mosaic across the retina, like the 'pale' and 'yellow' types. Variation in these ommatidia patterns very likely has important implications for the vision of insects and could underlie behavioural and environmental adaptations. However, despite the detailed understanding of ommatidia specification in D. melanogaster, the extent to which the frequency and distribution of the different ommatidium types vary between sexes, strains and species of Drosophila is not known.

Insulin/insulin-like growth factor signalling (IIS) has been described as one of the major pathways involved in growth control and homeostasis in multicellular organisms. Whereas its core components are well established, less is known about the molecular functions of IIS regulators. The adaptor molecule Lnk/SH2B has been implicated in IIS but the mechanism by which it promotes IIS activity has remained enigmatic.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.