High-resolution gene expression data from blastoderm embryos of the scuttle fly Megaselia abdita

Genomic Resources for the Scuttle Fly Megaselia abdita: A Model Organism for Comparative Developmental Studies in Flies

The order Diptera (true flies) holds promise as a model taxon in evolutionary developmental biology due to the inclusion of the model organism, Drosophila melanogaster, and the ability to cost-effectively rear many species in laboratories. One of them, the scuttle fly Megaselia abdita (Phoridae) has been used in evolutionary developmental biology for 30 years and is an excellent phylogenetic intermediate between fruit flies and mosquitoes but remains underdeveloped in genomic resources. Here, we present a de novo chromosome-level assembly and annotation of M. abdita and transcriptomes of 9 embryonic and 4 postembryonic stages. We also compare 9 stage-matched embryonic transcriptomes between M. abdita and D. melanogaster. Our analysis of these resources reveals extensive chromosomal synteny with D. melanogaster, 28 orphan genes with embryo-specific expression including a novel F-box LRR gene in M. abdita, and conserved and diverged features of gene expression dynamics between M. abdita and D. melanogaster. Collectively, our results provide a new reference for studying the diversification of developmental processes in flies.

Mapping Selection within Drosophila melanogaster Embryo's Anatomy

We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.

Gap Gene Regulatory Dynamics Evolve along a Genotype Network

Developmental gene networks implement the dynamic regulatory mechanisms that pattern and shape the organism. Over evolutionary time, the wiring of these networks changes, yet the patterning outcome is often preserved, a phenomenon known as “system drift.” System drift is illustrated by the gap gene network—involved in segmental patterning—in dipteran insects. In the classic model organism Drosophila melanogaster and the nonmodel scuttle fly Megaselia abdita, early activation and placement of gap gene expression domains show significant quantitative differences, yet the final patterning output of the system is essentially identical in both species. In this detailed modeling analysis of system drift, we use gene circuits which are fit to quantitative gap gene expression data in M. abdita and compare them with an equivalent set of models from D. melanogaster. The results of this comparative analysis show precisely how compensatory regulatory mechanisms achieve equivalent final patterns in both species. We discuss the larger implications of the work in terms of “genotype networks” and the ways in which the structure of regulatory networks can influence patterns of evolutionary change (evolvability).

Maternal co-ordinate gene regulation and axis polarity in the scuttle fly Megaselia abdita

Axis specification and segment determination in dipteran insects are an excellent model system for comparative analyses of gene network evolution. Antero-posterior polarity of the embryo is established through systems of maternal morphogen gradients. In Drosophila melanogaster, the anterior system acts through opposing gradients of Bicoid (Bcd) and Caudal (Cad), while the posterior system involves Nanos (Nos) and Hunchback (Hb) protein. These systems act redundantly. Both Bcd and Hb need to be eliminated to cause a complete loss of polarity resulting in mirror-duplicated abdomens, so-called bicaudal phenotypes. In contrast, knock-down of bcd alone is sufficient to induce double abdomens in non-drosophilid cyclorrhaphan dipterans such as the hoverfly Episyrphus balteatus or the scuttle fly Megaselia abdita. We investigate conserved and divergent aspects of axis specification in the cyclorrhaphan lineage through a detailed study of the establishment and regulatory effect of maternal gradients in M. abdita. Our results show that the function of the anterior maternal system is highly conserved in this species, despite the loss of maternal cad expression. In contrast, hb does not activate gap genes in this species. The absence of this activatory role provides a precise genetic explanation for the loss of polarity upon bcd knock-down in M. abdita, and suggests a general scenario in which the posterior maternal system is increasingly replaced by the anterior one during the evolution of the cyclorrhaphan dipteran lineage.

The basic head-to-tail polarity of an animal is established very early in development. In dipteran insects (flies, midges, and mosquitoes), polarity is established with the help of so-called morphogen gradients. Morphogens are regulatory proteins that are distributed as a concentration gradient, often involving diffusion from a localised source. This graded distribution then leads to the concentration-dependent activation of different target genes along the embryo’s axis. We examine this process, which differs to a surprising extent between dipteran species, in the scuttle fly Megaselia abdita, and compare our results to the model organism Drosophila melanogaster. In this way, we not only gain insights into how the mechanisms that establish polarity function differently in different species, but also how the system has evolved since these two flies shared a common ancestor. Specifically, we pin down the main difference between Drosophila and Megaselia in the altered function of the maternal Hunchback morphogen gradient, which activates target genes in the former, but not the latter species, where it has been completely replaced by the Bicoid morphogen during evolution.