Expression and regulation of caudal in the lower cyclorrhaphan fly Megaselia

Spatial mapping of single cells in the Drosophila embryo from transcriptomic data based on topological consistency

The advancement in single-cell RNA sequencing technologies allow us to obtain transcriptome at single cell resolution. However, the original spatial context of cells, a crucial knowledge for understanding cellular and tissue-level functions, is often lost during sequencing. To address this issue, the DREAM Single Cell Transcriptomics Challenge launched a community-wide effort to seek computational solutions for spatial mapping of single cells in tissues using single-cell RNAseq (scRNA-seq) data and a reference atlas obtained from in situ hybridization data. As a top-performing team in this competition, we approach this problem in three steps. The first step involves identifying a set of most informative genes based on the consistency between gene expression similarity and cell proximity. For this step, we propose two different approaches, i.e., an unsupervised approach that does not utilize the gold standard location of the cells provided by the challenge organizers, and a supervised approach that relies on the gold standard locations. In the second step, a Particle Swarm Optimization algorithm is used to optimize the weights of different genes in order to maximize matches between the predicted locations and the gold standard locations. Finally, the information embedded in the cell topology is used to improve the predicted cell-location scores by weighted averaging of scores from neighboring locations. Evaluation results based on DREAM scores show that our method accurately predicts the location of single cells, and the predictions lead to successful recovery of the spatial expression patterns for most of landmark genes. In addition, investigating the selected genes demonstrates that most predictive genes are cluster specific, and stable across our supervised and unsupervised gene selection frameworks. Overall, the promising results obtained by our methods in DREAM challenge demonstrated that topological consistency is a useful concept in identifying marker genes and constructing predictive models for spatial mapping of single cells.

Spatial mapping of single cells in the Drosophila embryo from transcriptomic data based on topological consistency

The advancement in single-cell RNA sequencing technologies allow us to obtain transcriptome at single cell resolution. However, the original spatial context of cells, a crucial knowledge for understanding cellular and tissue-level functions, is often lost during sequencing. To address this issue, the DREAM Single Cell Transcriptomics Challenge launched a community-wide effort to seek computational solutions for spatial mapping of single cells in tissues using single-cell RNAseq (scRNA-seq) data and a reference atlas obtained from in situ hybridization data. As a top-performing team in this competition, we approach this problem in three steps. The first step involves identifying a set of most informative genes based on the consistency between gene expression similarity and cell proximity. For this step, we propose two different approaches, i.e., an unsupervised approach that does not utilize the gold standard location of the cells provided by the challenge organizers, and a supervised approach that relies on the gold standard locations. In the second step, a Particle Swarm Optimization algorithm is used to optimize the weights of different genes in order to maximize matches between the predicted locations and the gold standard locations. Finally, the information embedded in the cell topology is used to improve the predicted cell-location scores by weighted averaging of scores from neighboring locations. Evaluation results based on DREAM scores show that our method accurately predicts the location of single cells, and the predictions lead to successful recovery of the spatial expression patterns for most of landmark genes. In addition, investigating the selected genes demonstrates that most predictive genes are cluster specific, and stable across our supervised and unsupervised gene selection frameworks. Overall, the promising results obtained by our methods in DREAM challenge demonstrated that topological consistency is a useful concept in identifying marker genes and constructing predictive models for spatial mapping of single cells.

Linking gene regulation to cell behaviors in the posterior growth zone of sequentially segmenting arthropods

Virtually all arthropods all arthropods add their body segments sequentially, one by one in an anterior to posterior progression. That process requires not only segment specification but typically growth and elongation. Here we review the functions of some of the key genes that regulate segmentation: Wnt, caudal, Notch pathway, and pair-rule genes, and discuss what can be inferred about their evolution. We focus on how these regulatory factors are integrated with growth and elongation and discuss the importance and challenges of baseline measures of growth and elongation. We emphasize a perspective that integrates the genetic regulation of segment patterning with the cellular mechanisms of growth and elongation.

Emerging developmental genetic model systems in holometabolous insects

The number of insect species that are amenable to functional genetic studies is growing rapidly and provides many new research opportunities in developmental and evolutionary biology. The holometabolous insects represent a disproportionate percentage of animal diversity and are thus well positioned to provide model species for a wide variety of developmental processes. Here we discuss emerging holometabolous models, and review some recent breakthroughs. For example, flies and midges were found to use structurally unrelated long-range pattern organizers, butterflies and moths revealed extensive pattern formation during oogenesis, new imaging possibilities in the flour beetle Tribolium castaneum showed how embryos break free of their extraembryonic membranes, and the complex genetics governing interspecies difference in head shape were revealed in Nasonia wasps.

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

Embryo development. A cysteine-clamp gene drives embryo polarity in the midge Chironomus

In the fruit fly Drosophila, head formation is driven by a single gene, bicoid, which generates head-to-tail polarity of the main embryonic axis. Bicoid deficiency results in embryos with tail-to-tail polarity and no head. However, most insects lack bicoid, and the molecular mechanism for establishing head-to-tail polarity is poorly understood. We have identified a gene that establishes head-to-tail polarity of the mosquito-like midge, Chironomus riparius. This gene, named panish, encodes a cysteine-clamp DNA binding domain and operates through a different mechanism than bicoid. This finding, combined with the observation that the phylogenetic distributions of panish and bicoid are limited to specific families of flies, reveals frequent evolutionary changes of body axis determinants and a remarkable opportunity to study gene regulatory network evolution.

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.

Quantitative system drift compensates for altered maternal inputs to the gap gene network of the scuttle fly Megaselia abdita

The segmentation gene network in insects can produce equivalent phenotypic outputs despite differences in upstream regulatory inputs between species. We investigate the mechanistic basis of this phenomenon through a systems-level analysis of the gap gene network in the scuttle fly Megaselia abdita (Phoridae). It combines quantification of gene expression at high spatio-temporal resolution with systematic knock-downs by RNA interference (RNAi). Initiation and dynamics of gap gene expression differ markedly between M. abdita and Drosophila melanogaster, while the output of the system converges to equivalent patterns at the end of the blastoderm stage. Although the qualitative structure of the gap gene network is conserved, there are differences in the strength of regulatory interactions between species. We term such network rewiring 'quantitative system drift'. It provides a mechanistic explanation for the developmental hourglass model in the dipteran lineage. Quantitative system drift is likely to be a widespread mechanism for developmental evolution.

SuperFly: a comparative database for quantified spatio-temporal gene expression patterns in early dipteran embryos

We present SuperFly (http://superfly.crg.eu), a relational database for quantified spatio-temporal expression data of segmentation genes during early development in different species of dipteran insects (flies, midges and mosquitoes). SuperFly has a special focus on emerging non-drosophilid model systems. The database currently includes data of high spatio-temporal resolution for three species: the vinegar fly Drosophila melanogaster, the scuttle fly Megaselia abdita and the moth midge Clogmia albipunctata. At this point, SuperFly covers up to 9 genes and 16 time points per species, with a total of 1823 individual embryos. It provides an intuitive web interface, enabling the user to query and access original embryo images, quantified expression profiles, extracted positions of expression boundaries and integrated datasets, plus metadata and intermediate processing steps. SuperFly is a valuable new resource for the quantitative comparative study of gene expression patterns across dipteran species. Moreover, it provides an interesting test set for systems biologists interested in fitting mathematical gene network models to data. Both of these aspects are essential ingredients for progress toward a more quantitative and mechanistic understanding of developmental evolution.

Heterochronic shifts in germband movements contribute to the rapid embryonic development of the coffin fly Megaselia scalaris

The coffin fly, Megaselia scalaris, is a species of medical and forensic importance and is increasingly being used for the study of genetics. Postmortem interval can be estimated based on the life stage of M. scalaris recovered from corpses, therefore many studies have addressed the duration of each life stage. These studies demonstrate that embryogenesis completes significantly faster in M. scalaris than in the congener Megaselia abdita and faster even than the 24 h needed for Drosophila melanogaster embryogenesis. However, until now it has been unclear if this increased speed is achieved by reducing developmental time across all embryonic stages or by the acceleration of individual stages and processes. Here I use time-lapse imaging to create a staging scheme for M. scalaris embryogenesis. Comparison of stages between D. melanogaster and both Megaselia species reveals that heterochronic shifts, simultaneous morphogenetic movements and compression of individual stages all contribute to the rapid development of M. scalaris.

Noncanonical expression of caudal during early embryogenesis in the pea aphid Acyrthosiphon pisum: maternal cad-driven posterior development is not conserved

Previously we identified anterior localization of hunchback (Aphb) mRNA in oocytes and early embryos of the parthenogenetic and viviparous pea aphid Acyrthosiphon pisum, suggesting that the breaking of anterior asymmetry in the oocytes leads to the formation of the anterior axis in embryos. In order to study posterior development in the asexual pea aphid, we cloned and analysed the developmental expression of caudal (Apcad), a posterior gene highly conserved in many animal phyla. We found that transcripts of Apcad were not detected in germaria, oocytes and embryos prior to the formation of the blastoderm in the asexual (viviparous) pea aphid. This unusual expression pattern differs from that of the existing insect models, including long- and short-germ insects, where maternal cad mRNA is passed to the early embryos and forms a posterior-anterior gradient. The first detectable Apcad expression occurred in the newly formed primordial germ cells and their adjacent blastodermal cells during late blastulation. From gastrulation onward, and as in other insects, Apcad mRNA is restricted to the posteriormost region of the germ band. Similarly, in the sexual (oviparous) oocytes we were able to identify anterior localization of Aphb mRNA but posterior localization of Apcad was not detected. This suggests that cad-driven posterior development is not conserved during early embryogenesis in asexual and sexual pea aphids.

Caudal is a negative regulator of the Anopheles IMD pathway that controls resistance to Plasmodium falciparum infection

Malaria parasite transmission depends upon the successful development of Plasmodium in its Anopheles mosquito vector. The mosquito's innate immune system constitutes a major bottleneck for parasite population growth. We show here that in Anopheles gambiae, the midgut-specific transcription factor Caudal acts as a negative regulator in the Imd pathway-mediated immune response against the human malaria parasite Plasmodium falciparum. Caudal also modulates the mosquito midgut bacterial flora. RNAi-mediated silencing of Caudal enhanced the mosquito's resistance to bacterial infections and increased the transcriptional abundance of key immune effector genes. Interestingly, Caudal's silencing resulted in an increased lifespan of the mosquito, while it impaired reproductive fitness with respect to egg laying and hatching.

Comparative transcriptomics of early dipteran development

Background

Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo).

Results

We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships.

Conclusions

We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

Life's attractors : understanding developmental systems through reverse engineering and in silico evolution

We propose an approach to evolutionary systems biology which is based on reverse engineering of gene regulatory networks and in silico evolutionary simulations. We infer regulatory parameters for gene networks by fitting computational models to quantitative expression data. This allows us to characterize the regulatory structure and dynamical repertoire of evolving gene regulatory networks with a reasonable amount of experimental and computational effort. We use the resulting network models to identify those regulatory interactions that are conserved, and those that have diverged between different species. Moreover, we use the models obtained by data fitting as starting points for simulations of evolutionary transitions between species. These simulations enable us to investigate whether such transitions are random, or whether they show stereotypical series of regulatory changes which depend on the structure and dynamical repertoire of an evolving network. Finally, we present a case study-the gap gene network in dipterans (flies, midges, and mosquitoes)-to illustrate the practical application of the proposed methodology, and to highlight the kind of biological insights that can be gained by this approach.

The gap gene network

Gap genes are involved in segment determination during the early development of the fruit fly Drosophila melanogaster as well as in other insects. This review attempts to synthesize the current knowledge of the gap gene network through a comprehensive survey of the experimental literature. I focus on genetic and molecular evidence, which provides us with an almost-complete picture of the regulatory interactions responsible for trunk gap gene expression. I discuss the regulatory mechanisms involved, and highlight the remaining ambiguities and gaps in the evidence. This is followed by a brief discussion of molecular regulatory mechanisms for transcriptional regulation, as well as precision and size-regulation provided by the system. Finally, I discuss evidence on the evolution of gap gene expression from species other than Drosophila. My survey concludes that studies of the gap gene system continue to reveal interesting and important new insights into the role of gene regulatory networks in development and evolution.

Maternal activation of gap genes in the hover fly Episyrphus

The metameric organization of the insect body plan is initiated with the activation of gap genes, a set of transcription-factor-encoding genes that are zygotically expressed in broad and partially overlapping domains along the anteroposterior (AP) axis of the early embryo. The spatial pattern of gap gene expression domains along the AP axis is generally conserved, but the maternal genes that regulate their expression are not. Building on the comprehensive knowledge of maternal gap gene activation in Drosophila, we used loss- and gain-of-function experiments in the hover fly Episyrphus balteatus (Syrphidae) to address the question of how the maternal regulation of gap genes evolved. We find that, in Episyrphus, a highly diverged bicoid ortholog is solely responsible for the AP polarity of the embryo. Episyrphus bicoid represses anterior zygotic expression of caudal and activates the anterior and central gap genes orthodenticle, hunchback and Krüppel. In bicoid-deficient Episyrphus embryos, nanos is insufficient to generate morphological asymmetry along the AP axis. Furthermore, we find that torso transiently regulates anterior repression of caudal and is required for the activation of orthodenticle, whereas all posterior gap gene domains of knirps, giant, hunchback, tailless and huckebein depend on caudal. We conclude that all maternal coordinate genes have altered their specific functions during the radiation of higher flies (Cyclorrhapha).

Hox3/zen and the evolution of extraembryonic epithelia in insects

Insects have undergone dramatic evolutionary changes in extraembryonic development, which correlate with changes in the expression of the class-3 Hox gene zen. Here, we review the evolution of this gene in insects and point out how changes in zen expression may have affected extraembryonic development at the morphological and the genetic level.

Evidence for a composite anterior determinant in the hover fly Episyrphus balteatus (Syrphidae), a cyclorrhaphan fly with an anterodorsal serosa anlage

Most insect embryos develop from a monolayer of cells around the yolk, but only part of this blastoderm forms the embryonic rudiment. Another part forms extra-embryonic serosa. Size and position of the serosa anlage vary between species, and previous work raises the issue of whether such differences co-evolve with the mechanisms that establish anteroposterior (AP) polarity of the embryo. AP polarity of the Drosophila embryo depends on bicoid, which is necessary and sufficient to determine the anterior body plan. Orthologs of bicoid have been identified in various cyclorrhaphan flies and their occurrence seems to correlate with a mid-dorsal serosa or amnioserosa anlage. Here, we introduce with Episyrphus balteatus (Syrphidae) a cyclorrhaphan model for embryonic AP axis specification that features an anterodorsal serosa anlage. Current phylogenies place Episyrphus within the clade that uses bicoid mRNA as anterior determinant, but no bicoid-like sequence could be identified in this species. Using RNA interference (RNAi) and ectopic mRNA injection, we obtained evidence that pattern formation along the entire AP axis of the Episyrphus embryo relies heavily on the precise regulation of caudal, and that anterior pattern formation in particular depends on two localized factors rather than one. Early zygotic activation of orthodenticle is separated from anterior repression of caudal, two distinct functions which in Drosophila are performed jointly by bicoid, whereas hunchback appears to be regulated by both factors. Furthermore, we found that overexpression of orthodenticle is sufficient to confine the serosa anlage of Episyrphus to dorsal blastoderm. We discuss our findings in a phylogenetic context and propose that Episyrphus employs a primitive cyclorrhaphan mechanism of AP axis specification.