How to get ahead: the origin, evolution and function ofbicoid

Body plan evolvability: The role of variability in gene regulatory networks

Recent computational modeling of early fruit fly (Drosophila) development has characterized the degree to which gene regulation networks can be robust to natural variability. In the first few hours of development, broad spatial gradients of maternally derived transcription factors activate embryonic gap genes. These gap patterns determine the subsequent segmented insect body plan through pair-rule gene expression. Gap genes are expressed with greater spatial precision than the maternal patterns. Computational modeling of the gap-gap regulatory interactions provides a mechanistic understanding for this robustness to maternal variability in wild-type (WT) patterning. A long-standing question in evolutionary biology has been how a system which is robust, such as the developmental program creating any particular species' body plan, is also evolvable, i.e. how can a system evolve or speciate, if the WT form is strongly buffered and protected? In the present work, we use the WT model to explore the breakdown of such Waddington-type 'canalization'. What levels of variability will push the system out of the WT form; are there particular pathways in the gene regulatory mechanism which are more susceptible to losing the WT form; and when robustness is lost, what types of forms are most likely to occur (i.e. what forms lie near the WT)? Manipulating maternal effects in several different pathways, we find a common gap 'peak-to-step' pattern transition in the loss of WT. We discuss these results in terms of the evolvability of insect segmentation, and in terms of experimental perturbations and mutations which could test the model predictions. We conclude by discussing the prospects for using continuum models of pattern dynamics to investigate a wider range of evo-devo problems.

Role of maternal spiralian-specific homeobox gene SPILE-E in the specification of blastomeres along the animal-vegetal axis during the early cleavage stages of mollusks

Spiralians, one of the major clades of bilaterians, share a unique development known as spiralian development, characterized by the formation of tiers of cells called quartets, which exhibit different developmental potentials along the animal-vegetal axis. Recently, spiralian-specific TALE-type homeobox genes (SPILE) have been identified, some of which show zygotic and staggered expression patterns along the animal-vegetal axis and function in quartet specification in mollusks. However, it is unclear which maternal molecular components control the zygotic expression of these transcription factors. In this study, we focused on SPILE-E, a maternal transcription factor, and investigated its expression and function in mollusks. We found that the maternal and ubiquitous expression of SPILE-E in the cleavage stages is conserved in molluskan species, including limpets, mussels, and chitons. We knocked down SPILE-E in limpets and revealed that the expression of transcription factors specifically expressed in the first quartet (1q2 ; foxj1b) and second quartet (2q; SPILE-B) was abolished, whereas the macromere-quartet marker (SPILE-C) was ectopically expressed in 1q2 in SPILE-E morphants. Moreover, we showed that the expression of SPILE-A, which upregulates SPILE-B but represses SPILE-C expression, decreased in SPILE-E morphants. Consistent with changes in the expression pattern of the above transcription factors, SPILE-E-morphant larvae exhibited patchy or complete loss of expression of marker genes of ciliated cells and shell fields, possibly reflecting incomplete specification of 1q2 and 2q. Our results provide a molecular framework for quartet specification and highlight the importance of maternal lineage-specific transcription factors in the development and evolution of spiralians.

Evolutionary Change in Gut Specification in Caenorhabditis Centers on the GATA Factor ELT-3 in an Example of Developmental System Drift

Cells in a developing animal embryo become specified by the activation of cell-type-specific gene regulatory networks. The network that specifies the gut in the nematode Caenorhabditis elegans has been the subject of study for more than two decades. In this network, the maternal factors SKN-1/Nrf and POP-1/TCF activate a zygotic GATA factor cascade consisting of the regulators MED-1,2 → END-1,3 → ELT-2,7, leading to the specification of the gut in early embryos. Paradoxically, the MED, END, and ELT-7 regulators are present only in species closely related to C. elegans, raising the question of how the gut can be specified without them. Recent work found that ELT-3, a GATA factor without an endodermal role in C. elegans, acts in a simpler ELT-3 → ELT-2 network to specify gut in more distant species. The simpler ELT-3 → ELT-2 network may thus represent an ancestral pathway. In this review, we describe the elucidation of the gut specification network in C. elegans and related species and propose a model by which the more complex network might have formed. Because the evolution of this network occurred without a change in phenotype, it is an example of the phenomenon of Developmental System Drift.

Transcriptional activators in the early Drosophila embryo perform different kinetic roles

Combinatorial regulation of gene expression by transcription factors (TFs) may in part arise from kinetic synergy-wherein TFs regulate different steps in the transcription cycle. Kinetic synergy requires that TFs play distinguishable kinetic roles. Here, we used live imaging to determine the kinetic roles of three TFs that activate transcription in the Drosophila embryo-Zelda, Bicoid, and Stat92E-by introducing their binding sites into the even-skipped stripe 2 enhancer. These TFs influence different sets of kinetic parameters, and their influence can change over time. All three TFs increased the fraction of transcriptionally active nuclei; Zelda also shortened the first-passage time into transcription and regulated the interval between transcription events. Stat92E also increased the lifetimes of active transcription. Different TFs can therefore play distinct kinetic roles in activating the transcription. This has consequences for understanding the composition and flexibility of regulatory DNA sequences and the biochemical function of TFs. A record of this paper's transparent peer review process is included in the supplemental information.

Insect PRXamides: Evolutionary Divergence, Novelty, and Loss in a Conserved Neuropeptide System

The PRXamide neuropeptides have been described in both protostome and deuterostome species, including all major groups of the Panarthropoda. Best studied are the insect PRXamides consisting of three genes: pk/pban, capa, and eth, each encoding multiple short peptides that are cleaved post-translationally. Comparisons of genome and transcriptome sequences reveal that while retaining its fundamental ancestral organization, the products of the pk/pban gene have undergone significant change in the insect Order Diptera. Basal dipteran pk/pban genes are much like those of other holometabolous insects, while more crown species have lost two peptide coding sequences including the otherwise ubiquitous pheromone biosynthesis activating neuropeptide (PBAN). In the genomic model species Drosophila melanogaster, one of the remaining peptides (hugin) plays a potentially novel role in feeding and locomotor regulation tied to circadian rhythms. Comparison of peptide coding sequences of pk/pban across the Diptera pinpoints the acquisition or loss of the hugin and PBAN peptide sequences respectively, and provides clues to associated changes in life history, physiology, and/or behavior. Interestingly, the neural circuitry underlying pk/pban function is highly conserved across the insects regardless of the composition of the pk/pban gene. The rapid evolution and diversification of the Diptera provide many instances of adaptive novelties from genes to behavior that can be placed in the context of emerging selective pressures at key points in their phylogeny; further study of changing functional roles of pk/pban may then be facilitated by the high-resolution genetic tools available in Drosophila melanogaster.

Rebuilding ships while at sea-Character individuality, homology, and evolutionary innovation

How novel traits originate in evolution is still one of the most perplexing questions in Evolutionary Biology. Building on a previous account of evolutionary innovation, I here propose that evolutionary novelties are those individualized characters that are not homologous to any characters in the ancestor. To clarify this definition, I here provide a detailed analysis of the concepts of "character individuality" and "homology" first, before addressing their role for our understanding of evolutionary innovation. I will argue (1) that functional as well as structural considerations are important for character individualization; and (2) that compositional (structural) and positional homology need to be clearly distinguished to properly describe the evolutionary transformations of hierarchically structured characters. My account will therefore integrate functional and structural perspectives and put forward a new multi-level view of character identity and transformation.

The GATA factor ELT-3 specifies endoderm in Caenorhabditis angaria in an ancestral gene network

Endoderm specification in Caenorhabditis elegans occurs through a network in which maternally provided SKN-1/Nrf, with additional input from POP-1/TCF, activates the GATA factor cascade MED-1,2→END-1,3→ELT-2,7. Orthologues of the MED, END and ELT-7 factors are found only among nematodes closely related to C. elegans, raising the question of how gut is specified in their absence in more distant species in the genus. We find that the C. angaria, C. portoensis and C. monodelphis orthologues of the GATA factor gene elt-3 are expressed in the early E lineage, just before their elt-2 orthologues. In C. angaria, Can-pop-1(RNAi), Can-elt-3(RNAi) and a Can-elt-3 null mutation result in a penetrant ‘gutless’ phenotype. Can-pop-1 is necessary for Can-elt-3 activation, showing that it acts upstream. Forced early E lineage expression of Can-elt-3 in C. elegans can direct the expression of a Can-elt-2 transgene and rescue an elt-7 end-1 end-3; elt-2 quadruple mutant strain to viability. Our results demonstrate an ancestral mechanism for gut specification and differentiation in Caenorhabditis involving a simpler POP-1→ELT-3→ELT-2 gene network.

Germline specification and axis determination in viviparous and oviparous pea aphids: conserved and divergent features

Aphids are hemimetabolous insects that undergo incomplete metamorphosis without pupation. The annual life cycle of most aphids includes both an asexual (viviparous) and a sexual (oviparous) phase. Sexual reproduction only occurs once per year and is followed by many generations of asexual reproduction, during which aphids propagate exponentially with telescopic development. Here, we discuss the potential links between viviparous embryogenesis and derived developmental features in the pea aphid Acyrthosiphon pisum, particularly focusing on germline specification and axis determination, both of which are key events of early development in insects. We also discuss potential evolutionary paths through which both viviparous and oviparous females might have come to utilize maternal germ plasm to drive germline specification. This developmental strategy, as defined by germline markers, has not been reported in other hemimetabolous insects. In viviparous females, furthermore, we discuss whether molecules that in other insects characterize germ plasm, like Vasa, also participate in posterior determination and how the anterior localization of the hunchback orthologue Ap-hb establishes the anterior-posterior axis. We propose that the linked chain of developing oocytes and embryos within each ovariole and the special morphology of early embryos might have driven the formation of evolutionary novelties in germline specification and axis determination in the viviparous aphids. Moreover, based upon the finding that the endosymbiont Buchnera aphidicola is closely associated with germ cells throughout embryogenesis, we propose presumptive roles for B. aphidicola in aphid development, discussing how it might regulate germline migration in both reproductive modes of pea aphids. In summary, we expect that this review will shed light on viviparous as well as oviparous development in aphids.

Neural development: Bicoid not as a morphogen

Bicoid is famous for its role in early embryo patterning of Drosophila by activating Hunchback expression to establish the anterior-posterior axis. A new study found that in a type of post-mitotic neuron Hunchback conversely activates Bicoid expression to regulate synapse targeting and locomotor behavior.

Hunchback activates Bicoid in Pair1 neurons to regulate synapse number and locomotor circuit function

Neural circuit function underlies cognition, sensation, and behavior. Proper circuit assembly depends on the identity of the neurons in the circuit (gene expression, morphology, synapse targeting, and biophysical properties). Neuronal identity is established by spatial and temporal patterning mechanisms, but little is known about how these mechanisms drive circuit formation in postmitotic neurons. Temporal patterning involves the sequential expression of transcription factors (TFs) in neural progenitors to diversify neuronal identity, in part through the initial expression of homeodomain TF combinations. Here, we address the role of the Drosophila temporal TF Hunchback and the homeodomain TF Bicoid in the assembly of the Pair1 (SEZ_DN1) descending neuron locomotor circuit, which promotes larval pausing and head casting. We find that both Hunchback and Bicoid are expressed in larval Pair1 neurons, Hunchback activates Bicoid in Pair1 (opposite of their embryonic relationship), and the loss of Hunchback function or Bicoid function from Pair1 leads to ectopic presynapse numbers in Pair1 axons and an increase in Pair1-induced pausing behavior. These phenotypes are highly specific, as the loss of Bicoid or Hunchback has no effect on Pair1 neurotransmitter identity, dendrite morphology, or axonal morphology. Importantly, the loss of Hunchback or Bicoid in Pair1 leads to the addition of new circuit partners that may underlie the exaggerated locomotor pausing behavior. These data are the first to show a role for Bicoid outside of embryonic patterning and the first to demonstrate a cell-autonomous role for Hunchback and Bicoid in interneuron synapse targeting and locomotor behavior.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Unexpected mutual regulation underlies paralogue functional diversification and promotes epithelial tissue maturation in Tribolium

Insect Hox3/zen genes represent an evolutionary hotspot for changes in function and copy number. Single orthologues are required either for early specification or late morphogenesis of the extraembryonic tissues, which protect the embryo. The tandemly duplicated zen paralogues of the beetle Tribolium castaneum present a unique opportunity to investigate both functions in a single species. We dissect the paralogues' expression dynamics (transcript and protein) and transcriptional targets (RNA-seq after RNAi) throughout embryogenesis. We identify an unexpected role of Tc-Zen2 in repression of Tc-zen1, generating a negative feedback loop that promotes developmental progression. Tc-Zen2 regulation is dynamic, including within co-expressed multigene loci. We also show that extraembryonic development is the major event within the transcriptional landscape of late embryogenesis and provide a global molecular characterization of the extraembryonic serosal tissue. Altogether, we propose that paralogue mutual regulation arose through multiple instances of zen subfunctionalization, leading to their complementary extant roles.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Arthropod segmentation

There is now compelling evidence that many arthropods pattern their segments using a clock-and-wavefront mechanism, analogous to that operating during vertebrate somitogenesis. In this Review, we discuss how the arthropod segmentation clock generates a repeating sequence of pair-rule gene expression, and how this is converted into a segment-polarity pattern by ‘timing factor’ wavefronts associated with axial extension. We argue that the gene regulatory network that patterns segments may be relatively conserved, although the timing of segmentation varies widely, and double-segment periodicity appears to have evolved at least twice. Finally, we describe how the repeated evolution of a simultaneous (Drosophila-like) mode of segmentation within holometabolan insects can be explained by heterochronic shifts in timing factor expression plus extensive pre-patterning of the pair-rule genes.

The Origin and Evolution of Maternal Genes

Proteins and RNA molecules are deposited into the developing egg by the mother. These gene products will drive the first stages of development and are coded by maternal genes. Maternal genes are essential, yet, despite their importance, their evolutionary dynamics is largely unknown. Here I review the current knowledge of maternal gene evolution. The evolutionary origin of maternal genes tends to be more recent than that of zygotic genes. Some studies support the theoretical prediction that maternal genes evolve faster than zygotic genes. However, most studies were done on a limited set of species and genes. I also discuss the way forward to understand the evolution of maternal genes by combining high-throughput genomics and theoretical evolutionary approaches.

Evolution of maternal control of axial patterning in insects

Positional and cell fate cues provided maternally to eggs are important factors in the development of many animals. The insects are a model clade where maternal establishment of embryonic axes is widespread and has been a topic of intense classical and molecular embryological analysis. Recently, significant progress has been made in revealing the molecular basis of some classical embryological experiments. In addition, observations of novel forms of maternal positional cues have been made. Finally, it has become increasingly clear that no maternal source of positional information acts alone without input and feedback from zygotic target genes to ensure precise and repeatable pattern formation in the early embryo. These advances will be discussed in the context of historical experiments, our current understanding of how positional cues can be generated, stored, and transmitted in insect ovaries and eggs, and how the nature of the cues can change in evolution.

Evolution of maternal and zygotic mRNA complements in the early Drosophila embryo

The earliest stages of animal development are controlled by maternally deposited mRNA transcripts and proteins. Once the zygote is able to transcribe its own genome, maternal transcripts are degraded, in a tightly regulated process known as the maternal to zygotic transition (MZT). While this process has been well-studied within model species, we have little knowledge of how the pools of maternal and zygotic transcripts evolve. To characterize the evolutionary dynamics and functional constraints on early embryonic expression, we created a transcriptomic dataset for 14 Drosophila species spanning over 50 million years of evolution, at developmental stages before and after the MZT, and compared our results with a previously published Aedes aegypti developmental time course. We found deep conservation over 250 million years of a core set of genes transcribed only by the zygote. This select group is highly enriched in transcription factors that play critical roles in early development. However, we also identify a surprisingly high level of change in the transcripts represented at both stages over the phylogeny. While mRNA levels of genes with maternally deposited transcripts are more highly conserved than zygotic genes, those maternal transcripts that are completely degraded at the MZT vary dramatically between species. We also show that hundreds of genes have different isoform usage between the maternal and zygotic genomes. Our work suggests that maternal transcript deposition and early zygotic transcription are remarkably dynamic over evolutionary time, despite the widespread conservation of early developmental processes.

Cell cycle regulator, small silencing RNA, and segmentation patterning gene expression in relation to embryonic diapause in the band-legged ground cricket

Insects enter diapause to synchronize their life cycle with biotic and abiotic conditions favorable for their development, reproduction, and survival. Adult females of the band-legged ground cricket Dianemobius nigrofasciatus (Orthoptera, Glyllidae) respond to environmental factors in autumn and lay diapause-destined eggs. The eggs arrest their development and enter diapause at a very early embryonic stage, specifically the cellular blastoderm. To elucidate the physiological mechanisms underlying this very early stage programmed developmental arrest, we investigated the cell division cycle as well as the expression of cell cycle regulators, small silencing RNAs, and segment patterning genes. The diapause embryo arrests its cell cycle predominantly at the G₀/G₁ phase. The proportion of cells in the S phase of the cell cycle abruptly decreased at the time of developmental arrest, but further changes of the G₀/G₁ and G₂/M were later observed. Thus, cell cycle arrest in the diapause embryo is not an immediate event, but it takes longer to reach the steady state. We further elucidated molecular events possibly involved in diapause preparation and entry. Downregulation of Proliferating cellular antigen (PCNA; a cell cycle regulator), caudal and pumilio (cad and pum; early segmentation genes) as well as P-element induced wimpy testis (piwi) (a small silencing RNA) prior to the onset of developmental arrest was notable. The downregulation of PCNA, cad and pum continued even after entry into developmental arrest. In contrast to upregulation in non-diapause eggs, Cyclin D (another cell cycle regulator) and hunchback, Krüppel, and runt (gap and pair-rule genes) were downregulated in diapause eggs. These molecular events may contribute to embryonic diapause of D. nigrofasciatus.

Evolutionary coincidence of adaptive changes in exuperantia and the emergence of bicoid in Cyclorrhapha (Diptera)

The great radiation in the infraorder Cyclorrhapha involved several morphological and molecular changes, including important changes in anterior egg development. During Drosophila oogenesis, exuperantia (exu) is critical for localizing bicoid (bcd) messenger RNA (mRNA) to the anterior region of the oocyte. Because it is phylogenetically older than bcd, which is exclusive to Cyclorrhapha, we hypothesize that exu has undergone adaptive changes to enable this new function. Although exu has been well studied in Drosophila, there is no functional or transcriptional information about it in any other Diptera. Here, we investigate exu in the South American fruit fly Anastrepha fraterculus, a Cyclorrhapha of great agricultural importance that have lost bcd, aiming to understand the evolution of exu in this infraorder. We assessed its pattern of gene expression in A. fraterculus by analyzing transcriptomes from cephalic and reproductive tissues. A combination of next-generation data with classical sequencing procedures enabled identification of the structure of exu and its alternative transcripts in this species. In addition to the sex-specific isoforms described for Drosophila, we found that not only exu is expressed in heads, but this is mediated by two transcripts with a specific 5′UTR exon—likely a result from usage of a third promoter. Furthermore, we tested the hypothesis that exu is evolving under positive selection in Cyclorrhapha after divergence from lower Diptera. We found evidence of positive selection at two important exu domains, EXO-like and SAM-like, both involved with mRNA binding during bcd mRNA localization in Drosophila, which could reflect its cooptation for the new function of bcd mRNA localization in Cyclorrhapha.

A molecular view of onychophoran segmentation

This paper summarizes our current knowledge on the expression and assumed function of Drosophila and (other) arthropod segmentation gene orthologs in Onychophora, a closely related outgroup to Arthropoda. This includes orthologs of the so-called Drosophila segmentation gene cascade including the Hox genes, as well as other genetic factors and pathways involved in non-drosophilid arthropods. Open questions about and around the topic are addressed, such as the definition of segments in onychophorans, the unclear regulation of conserved expression patterns downstream of non-conserved factors, and the potential role of mesodermal patterning in onychophoran segmentation.

Hox genes, evo-devo, and the case of the ftz gene

The discovery of the broad conservation of embryonic regulatory genes across animal phyla, launched by the cloning of homeotic genes in the 1980s, was a founding event in the field of evolutionary developmental biology (evo-devo). While it had long been known that fundamental cellular processes, commonly referred to as housekeeping functions, are shared by animals and plants across the planet-processes such as the storage of information in genomic DNA, transcription, translation and the machinery for these processes, universal codon usage, and metabolic enzymes-Hox genes were different: mutations in these genes caused "bizarre" homeotic transformations of insect body parts that were certainly interesting but were expected to be idiosyncratic. The isolation of the genes responsible for these bizarre phenotypes turned out to be highly conserved Hox genes that play roles in embryonic patterning throughout Metazoa. How Hox genes have changed to promote the development of diverse body plans remains a central issue of the field of evo-devo today. For this Memorial article series, I review events around the discovery of the broad evolutionary conservation of Hox genes and the impact of this discovery on the field of developmental biology. I highlight studies carried out in Walter Gehring's lab and by former lab members that have continued to push the field forward, raising new questions and forging new approaches to understand the evolution of developmental mechanisms.

Evolution and development of the adelphophagic, intracapsular Schmidt's larva of the nemertean Lineus ruber

Background

The life cycle of many animals includes a larval stage, which has diversified into an astonishing variety of ecological strategies. The Nemertea is a group of spiralians that exhibits a broad diversity of larval forms, including the iconic pilidium. A pelagic planktotrophic pilidium is the ancestral form in the Pilidiophora, but several lineages exhibit deviations of this condition, mostly as a transition to pelagic lecithotrophy. The most extreme case occurs, however, in the Pilidiophoran Lineus ruber, which exhibits an adelphophagic intracapsular pilidium, the so-called Schmidt’s larva.

Results

We combined confocal laser scanning microscopy and gene expression studies to characterize the development and metamorphosis of the Schmidt’s larva of L. ruber. The larva forms after gastrulation, and comprises a thin epidermis, a proboscis rudiment and two pairs of imaginal discs from which the juvenile will develop. The cells internalized during gastrulation form a blind gut and the blastopore gives rise to the mouth of the larva and juvenile. The Schmidt’s larva eats other siblings that occupy the same egg capsule, accumulating nutrients for the juvenile. A gradual metamorphosis involves the differentiation of the juvenile cell types from the imaginal discs and the shedding of the larval epidermis. The expression of evolutionarily conserved anterior (foxQ2, six3/6, gsc, otx), endomesodermal (foxA, GATA456-a, twi-a) and posterior (evx, cdx) markers demonstrate that the juvenile retains the molecular patterning of the Schmidt’s larva. After metamorphosis, the juveniles stay over 20 days within the egg masses, until they are fully mature and hatch.

Conclusions

The evolution of the intracapsular Schmidt’s larva involved the loss of the typical feeding structures of the planktotrophic pilidium and a precocious formation of the imaginal discs, as also observed in other pelagic lecithotrophic forms. However, no special adaptations are observed related to adelphophagy. As in planktotrophic pilidium, the molecular mechanism patterning the juvenile is only active in the imaginal discs and not during the early development of the larva, suggesting two separate molecular programs during nemertean embryogenesis. Our results illuminate the diversification of larval forms in the Pilidiophora and Nemertea, and thus on the developmental mechanisms underlying metazoan larval evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0023-5) contains supplementary material, which is available to authorized users.

The evolutionary history of vertebrate cranial placodes II. Evolution of ectodermal patterning

Cranial placodes are evolutionary innovations of vertebrates. However, they most likely evolved by redeployment, rewiring and diversification of preexisting cell types and patterning mechanisms. In the second part of this review we compare vertebrates with other animal groups to elucidate the evolutionary history of ectodermal patterning. We show that several transcription factors have ancient bilaterian roles in dorsoventral and anteroposterior regionalisation of the ectoderm. Evidence from amphioxus suggests that ancestral chordates then concentrated neurosecretory cells in the anteriormost non-neural ectoderm. This anterior proto-placodal domain subsequently gave rise to the oral siphon primordia in tunicates (with neurosecretory cells being lost) and anterior (adenohypophyseal, olfactory, and lens) placodes of vertebrates. Likewise, tunicate atrial siphon primordia and posterior (otic, lateral line, and epibranchial) placodes of vertebrates probably evolved from a posterior proto-placodal region in the tunicate-vertebrate ancestor. Since both siphon primordia in tunicates give rise to sparse populations of sensory cells, both proto-placodal domains probably also gave rise to some sensory receptors in the tunicate-vertebrate ancestor. However, proper cranial placodes, which give rise to high density arrays of specialised sensory receptors and neurons, evolved from these domains only in the vertebrate lineage. We propose that this may have involved rewiring of the regulatory network upstream and downstream of Six1/2 and Six4/5 transcription factors and their Eya family cofactors. These proteins, which play ancient roles in neuronal differentiation were first recruited to the dorsal non-neural ectoderm in the tunicate-vertebrate ancestor but subsequently probably acquired new target genes in the vertebrate lineage, allowing them to adopt new functions in regulating proliferation and patterning of neuronal progenitors.

Hunchback is required for abdominal identity suppression and germband growth in the parthenogenetic embryogenesis of the pea aphid, Acyrthosiphon pisum

Aphid, a short germband insect, displays an embryogenesis different from that of long germband insect species. Furthermore, the development of its parthenogenetic and viviparous embryo is different from that of the embryo resulting from sexual reproduction. To better understand the genetic regulation of this type of embryogenesis, the functions of hunchback in asexual Acyrthosiphon pisum were investigated by parental RNAi. Microinjection of Aphb double-stranded RNA yielded several defective phenotypes. Quantitative real-time PCR analysis revealed that these defects resulted from reduction of Aphb mRNA level in injected aphids. All these results suggested that the hb gene in parthenogenetic and viviparous Acyrthosiphon pisum was involved in abdominal identity suppression and germband growth as its homologue does in sexual insects.

Using the developmental gene bicoid to identify species of forensically important blowflies (Diptera: calliphoridae)

Identifying species of insects used to estimate postmortem interval (PMI) is a major subject in forensic entomology. Because forensic insect specimens are morphologically uniform and are obtained at various developmental stages, DNA markers are greatly needed. To develop new autosomal DNA markers to identify species, partial genomic sequences of the bicoid (bcd) genes, containing the homeobox and its flanking sequences, from 12 blowfly species (Aldrichina grahami, Calliphora vicina, Calliphora lata, Triceratopyga calliphoroides, Chrysomya megacephala, Chrysomya pinguis, Phormia regina, Lucilia ampullacea, Lucilia caesar, Lucilia illustris, Hemipyrellia ligurriens and Lucilia sericata; Calliphoridae: Diptera) were determined and analyzed. This study first sequenced the ten blowfly species other than C. vicina and L. sericata. Based on the bcd sequences of these 12 blowfly species, a phylogenetic tree was constructed that discriminates the subfamilies of Calliphoridae (Luciliinae, Chrysomyinae, and Calliphorinae) and most blowfly species. Even partial genomic sequences of about 500 bp can distinguish most blowfly species. The short intron 2 and coding sequences downstream of the bcd homeobox in exon 3 could be utilized to develop DNA markers for forensic applications. These gene sequences are important in the evolution of insect developmental biology and are potentially useful for identifying insect species in forensic science.

Variation and constraint in Hox gene evolution

Despite enormous body plan variation, genes regulating embryonic development are highly conserved. Here, we probe the mechanisms that predispose ancient regulatory genes to reutilization and diversification rather than evolutionary loss. The Hox gene fushi tarazu (ftz) arose as a homeotic gene but functions as a pair-rule segmentation gene in Drosophila. ftz shows extensive variation in expression and protein coding regions but has managed to elude loss from arthropod genomes. We asked what properties prevent this loss by testing the importance of different protein motifs and partners in the developing CNS, where ftz expression is conserved. Drosophila Ftz proteins with mutated protein motifs were expressed under the control of a neurogenic-specific ftz cis-regulatory element (CRE) in a ftz mutant background rescued for segmentation defects. Ftz CNS function did not require the variable motifs that mediate differential cofactor interactions involved in homeosis or segmentation, which vary in arthropods. Rather, CNS function did require the shared DNA-binding homeodomain, which plays less of a role in Ftz segmentation activity. The Antennapedia homeodomain substituted for Ftz homeodomain function in the Drosophila CNS, but full-length Antennapedia did not rescue CNS defects. These results suggest that a core CNS function retains ftz in arthropod genomes. Acquisition of a neurogenic CRE led to ftz expression in unique CNS cells, differentiating its role from neighboring Hox genes, rendering it nonredundant. The inherent flexibility of modular CREs and protein domains allows for stepwise acquisition of new functions, explaining broad retention of regulatory genes during animal evolution.

Conservation and variation in Hox genes: how insect models pioneered the evo-devo field

Evolutionary developmental biology, or evo-devo, broadly investigates how body plan diversity and morphological novelties have arisen and persisted in nature. The discovery of Hox genes in Drosophila, and their subsequent identification in most other metazoans, led biologists to try to understand how embryonic genes crucial for proper development have changed to promote the vast morphological variation seen in nature. Insects are ideal model systems for studying this diversity and the mechanisms underlying it because phylogenetic relationships are well established, powerful genetic tools have been developed, and there are many examples of evolutionary specializations that have arisen in nature in different insect lineages, such as the jumping leg of orthopterans and the helmet structures of treehoppers. Here, we briefly introduce the field of evo-devo and Hox genes, discuss functional tools available to study early developmental genes in insects, and provide examples in which changes in Hox genes have contributed to changes in body plan or morphology.

Evolutionary crossroads in developmental biology: the spider Parasteatoda tepidariorum

Spiders belong to the chelicerates, which is an arthropod group that branches basally from myriapods, crustaceans and insects. Spiders are thus useful models with which to investigate whether aspects of development are ancestral or derived with respect to the arthropod common ancestor. Moreover, they serve as an important reference point for comparison with the development of other metazoans. Therefore, studies of spider development have made a major contribution to advancing our understanding of the evolution of development. Much of this knowledge has come from studies of the common house spider, Parasteatoda tepidariorum. Here, we describe how the growing number of experimental tools and resources available to study Parasteatoda development have provided novel insights into the evolution of developmental regulation and have furthered our understanding of metazoan body plan evolution.

Expression of pair rule gene orthologs in the blastoderm of a myriapod: evidence for pair rule-like mechanisms?

Background

A hallmark of Drosophila segmentation is the stepwise subdivision of the body into smaller and smaller units, and finally into the segments. This is achieved by the function of the well-understood segmentation gene cascade. The first molecular sign of a segmented body appears with the action of the pair rule genes, which are expressed as transversal stripes in alternating segments. Drosophila development, however, is derived, and in most other arthropods only the anterior body is patterned (almost) simultaneously from a pre-existing field of cells; posterior segments are added sequentially from a posterior segment addition zone. A long-standing question is to what extent segmentation mechanisms known from Drosophila may be conserved in short-germ arthropods. Despite the derived developmental modes, it appears more likely that conserved mechanisms can be found in anterior patterning.

Results

Expression analysis of pair rule gene orthologs in the blastoderm of the pill millipede Glomeris marginata (Myriapoda: Diplopoda) suggests that these genes are generally involved in segmenting the anterior embryo. We find that the Glomeris pairberry-1 ( pby-1) gene is expressed in a pair rule pattern that is also found in insects and a chelicerate, the mite Tetraynchus urticae. Other Glomeris pair rule gene orthologs are expressed in double segment wide domains in the blastoderm, which at subsequent stages split into two stripes in adjacent segments.

Conclusions

The expression patterns of the millipede pair rule gene orthologs resemble pair rule patterning in Drosophila and other insects, and thus represent evidence for the presence of an ancestral pair rule-like mechanism in myriapods. We discuss the possibilities that blastoderm patterning may be conserved in long-germ and short-germ arthropods, and that a posterior double segmental mechanism may be present in short-germ arthropods.

Anterior-posterior patterning in early development: three strategies

The anterior-posterior (AP) axis is the most ancient of the embryonic axes and exists in most metazoans. Different animals use a wide variety of mechanisms to create this axis in the early embryo. In this study, we focus on three animals, including two insects (Drosophila and Tribolium) and a vertebrate (zebrafish) to examine different strategies used to form the AP axis. While Drosophila forms the entire axis within a syncytial blastoderm using transcription factors as morphogens, zebrafish uses signaling factors in a cellularized embryo, progressively forming the AP axis over the course of a day. Tribolium uses an intermediate strategy that has commonalities with both Drosophila and zebrafish. We discuss the specific molecular mechanisms used to create the AP axis and identify conserved features.

Novel function of Distal-less as a gap gene during spider segmentation

Despite many aspects of the regulation of segmentation being conserved among arthropods, the evolution of novel gene functions has played an important role in the evolution of developmental regulation and the emergence of new segmental structures. Moreover the study of such novel gene functions can be informative with respect to the patterns and direction of evolutionary changes in developmental programs. The homeobox gene Distal-less (Dll) is known for its conserved function in appendage development in metazoans. In arthropods, Dll is required for the specification of distal appendage structures. Here we describe a novel and unexpected role of Dll in the spider Achaearanea tepidariorum. We detect At-Dll transcripts not only in the appendages, but unexpectedly also in an anterior domain during early development, prior to the specification of the limb primordia. A similar early Dll domain is present in the distantly related spider Pholcus phalangioides. In A. tepidariorum this early At-Dll expression is required for head segmentation. RNA interference results in spiders that lack either the first or the first and the second walking leg segments. The early At-Dll expression is also required for the activation of the segment polarity genes engrailed and hedgehog in this region. Our work identifies the Distal-less gene as a novel factor in anterior spider segmentation with a gap gene-like function. This novel role of Dll is interesting because Dll expression is reduced in this region in crustaceans and the homologous insect segment, the mandible segment, does not express Dll and does not require this gene for patterning. We therefore discuss the possible implications of our results for understanding the evolution and diversification of the mandible segment.

The development and segmentation of the head of the fly Drosophila is one of the best-studied examples of how tissues become genetically specified during embryonic development. However, the mechanisms for head segmentation vary considerably among the arthropods. This is on the one hand surprising because the head consists of the same series of segments in all arthropods. On the other hand, changes in gene regulatory networks are the basis for the evolution of novel morphologies and biodiversity. In this paper, we describe our study of the Distal-less gene in the spider Achaearanea tepidariorum. Distal-less is well-known for its function in appendage development, but here we show that in the spider it also has a novel function during head segmentation that is not found in Drosophila or other arthropods. In Achaearanea the Distal-less gene is necessary for the formation of the first walking-leg segment, which is homologous to the mandible segment of the head of other arthropods. Intriguingly, the mandible segment in other arthropods exhibits reduced or missing Distal-less expression. Thus, our results suggest that this difference in the role of Distal-less in the first walking-leg/mandible segment of spiders and other arthropods may underlie the diversification of this segment.

Mathematical modeling of gene expression: a guide for the perplexed biologist

The detailed analysis of transcriptional networks holds a key for understanding central biological processes, and interest in this field has exploded due to new large-scale data acquisition techniques. Mathematical modeling can provide essential insights, but the diversity of modeling approaches can be a daunting prospect to investigators new to this area. For those interested in beginning a transcriptional mathematical modeling project, we provide here an overview of major types of models and their applications to transcriptional networks. In this discussion of recent literature on thermodynamic, Boolean, and differential equation models, we focus on considerations critical for choosing and validating a modeling approach that will be useful for quantitative understanding of biological systems.

Comparative genomic analysis of Drosophila melanogaster and vector mosquito developmental genes

Genome sequencing projects have presented the opportunity for analysis of developmental genes in three vector mosquito species: Aedes aegypti, Culex quinquefasciatus, and Anopheles gambiae. A comparative genomic analysis of developmental genes in Drosophila melanogaster and these three important vectors of human disease was performed in this investigation. While the study was comprehensive, special emphasis centered on genes that 1) are components of developmental signaling pathways, 2) regulate fundamental developmental processes, 3) are critical for the development of tissues of vector importance, 4) function in developmental processes known to have diverged within insects, and 5) encode microRNAs (miRNAs) that regulate developmental transcripts in Drosophila. While most fruit fly developmental genes are conserved in the three vector mosquito species, several genes known to be critical for Drosophila development were not identified in one or more mosquito genomes. In other cases, mosquito lineage-specific gene gains with respect to D. melanogaster were noted. Sequence analyses also revealed that numerous repetitive sequences are a common structural feature of Drosophila and mosquito developmental genes. Finally, analysis of predicted miRNA binding sites in fruit fly and mosquito developmental genes suggests that the repertoire of developmental genes targeted by miRNAs is species-specific. The results of this study provide insight into the evolution of developmental genes and processes in dipterans and other arthropods, serve as a resource for those pursuing analysis of mosquito development, and will promote the design and refinement of functional analysis experiments.

Evolution of complex gene regulatory circuits by addition of refinements

How do complex gene regulatory circuits evolve? These circuits involve many interacting components, which work together to specify patterns of gene expression. They typically include many subtle mechanistic features, but in most cases it is unclear whether these features are essential for the circuit to work at all, or if instead they make a functional circuit work better. In the latter case, such a feature is here termed 'dispensable', and it is plausible that the feature has been added at a late stage in the evolution of the circuit. This review describes experimental tests of this question, using the phage λ gene regulatory circuit. Several features of this circuit are found to be dispensable, in the sense that the circuitry works without these features, though not as well as the wild type. In some cases, second-site suppressor mutations are needed to confer near-normal behavior in the absence of such a feature. These findings are discussed here in the context of a two-stage model for evolution of gene regulatory circuits. In this model, a circuit evolves by assembly of a primitive or basic form, followed by adjustment of parameters and addition of qualitatively new features. Pathways are suggested for the addition of such features to a more basic form. Selected examples in other systems are described. Some of the dispensable features of phage λ may be evolutionary refinements. Finding that a feature is dispensable, however, does not prove that it is a late addition - it is possible that it was essential early in evolution, and became dispensable as the circuit evolved. Conversely, a late addition might have become essential. As ongoing work provides additional examples of dispensable features, it may become clearer how often they represent refinements.

Gene expression during blow fly development: improving the precision of age estimates in forensic entomology

Forensic entomologists use size and developmental stage to estimate blow fly age, and from those, a postmortem interval. Since such estimates are generally accurate but often lack precision, particularly in the older developmental stages, alternative aging methods would be advantageous. Presented here is a means of incorporating developmentally regulated gene expression levels into traditional stage and size data, with a goal of more precisely estimating developmental age of immature Lucilia sericata. Generalized additive models of development showed improved statistical support compared to models that did not include gene expression data, resulting in an increase in estimate precision, especially for postfeeding third instars and pupae. The models were then used to make blind estimates of development for 86 immature L. sericata raised on rat carcasses. Overall, inclusion of gene expression data resulted in increased precision in aging blow flies.

Rapid evolution of a novel signalling mechanism by concerted duplication and divergence of a BMP ligand and its extracellular modulators

Gene duplication and divergence is widely considered to be a fundamental mechanism for generating evolutionary novelties. The Bone Morphogenetic Proteins (BMPs) are a diverse family of signalling molecules found in all metazoan genomes that have evolved by duplication and divergence from a small number of ancestral types. In the fruit fly Drosophila, there are three BMPs: Decapentaplegic (Dpp) and Glass bottom boat (Gbb), which are the orthologues of vertebrate BMP2/4 and BMP5/6/7/8, respectively, and Screw (Scw), which, at the sequence level, is equally divergent from Dpp and Gbb. It has recently been shown that Scw has arisen from a duplication of Gbb in the lineage leading to higher Diptera. We show that since this duplication event, Gbb has maintained the ancestral BMP5/6/7/8 functionality while Scw has rapidly diverged. The evolution of Scw was accompanied by duplication and divergence of a suite of extracellular regulators that continue to diverge together in the higher Diptera. In addition, Scw has become restricted in its receptor specificity: Gbb proteins can signal through the Type I receptors Thick veins (Tkv) and Saxophone (Sax), while Scw signals through Sax. Thus, in a relatively short span of evolutionary time, the duplication event that gave rise to Scw produced not only a novel ligand but also a novel signalling mode that is functionally distinct from the ancestral Gbb mode. Our results demonstrate the plasticity of the BMP pathway not only in evolving new family members and new functions but also new signalling modes by redeploying key regulators in the pathway.

BSQA: integrated text mining using entity relation semantics extracted from biological literature of insects

Text mining is one promising way of extracting information automatically from the vast biological literature. To maximize its potential, the knowledge encoded in the text should be translated to some semantic representation such as entities and relations, which could be analyzed by machines. But large-scale practical systems for this purpose are rare. We present BeeSpace question/answering (BSQA) system that performs integrated text mining for insect biology, covering diverse aspects from molecular interactions of genes to insect behavior. BSQA recognizes a number of entities and relations in Medline documents about the model insect, Drosophila melanogaster. For any text query, BSQA exploits entity annotation of retrieved documents to identify important concepts in different categories. By utilizing the extracted relations, BSQA is also able to answer many biologically motivated questions, from simple ones such as, which anatomical part is a gene expressed in, to more complex ones involving multiple types of relations. BSQA is freely available at http://www.beespace.uiuc.edu/QuestionAnswer.

A compartmental model for the bicoid gradient

The anterior region of the Drosophila embryo is patterned by the concentration gradient of the homeodomain transcription factor bicoid (Bcd). The Bcd gradient was the first identified morphogen gradient and continues to be a subject of intense research at multiple levels, from the mechanisms of RNA localization in the oocyte to the evolution of the Bcd-mediated patterning events in multiple Drosophila species. Critical assessment of the mechanisms of the Bcd gradient formation requires biophysical models of the syncytial embryo. Most of the proposed models rely on reaction-diffusion equations, but their formulation and applicability at high nuclear densities is a nontrivial task. We propose a straightforward alternative in which the syncytial blastoderm is approximated by a periodic arrangement of well-mixed compartments: a single nucleus and an associated cytoplasmic region. We formulate a compartmental model, constrain its parameters by experimental data, and demonstrate that it provides an adequate description of the Bcd gradient dynamics.

Anterior development in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum: hunchback and orthodenticle expression

In the dipteran Drosophila, the genes bicoid and hunchback work synergistically to pattern the anterior blastoderm during embryogenesis. bicoid, however, appears to be an innovation of the higher Diptera. Hence, in some non-dipteran insects, anterior specification instead relies on a synergistic interaction between maternally transcribed hunchback and orthodenticle. Here we describe how orthologues of hunchback and orthodenticle are expressed during oogenesis and embryogenesis in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum. A. pisum hunchback (Aphb) mRNA is localized to the anterior pole in developing oocytes and early embryos prior to blastoderm formation - a pattern strongly reminiscent of bicoid localization in Drosophila. A. pisum orthodenticle (Apotd), on the other hand, is not expressed prior to gastrulation, suggesting that it is the asymmetric localization of Aphb, rather than synergy between Aphb and Apotd, that regulates anterior specification in asexual pea aphids.

Self-organization of intracellular gradients during mitosis

Gradients are used in a number of biological systems to transmit spatial information over a range of distances. The best studied are morphogen gradients where information is transmitted over many cell lengths. Smaller mitotic gradients reflect the need to organize several distinct events along the length of the mitotic spindle. The intracellular gradients that characterize mitosis are emerging as important regulatory paradigms. Intracellular gradients utilize intrinsic auto-regulatory feedback loops and diffusion to establish stable regions of activity within the mitotic cytosol. We review three recently described intracellular mitotic gradients. The Ran GTP gradient with its elaborate cascade of nuclear transport receptors and cargoes is the best characterized, yet the dynamics underlying the robust gradient of Ran-GTP have received little attention. Gradients of phosphorylation have been observed on Aurora B kinase substrates both before and after anaphase onset. In both instances the phosphorylation gradient appears to result from a soluble gradient of Aurora B kinase activity. Regulatory properties that support gradient formation are highlighted. Intracellular activity gradients that regulate localized mitotic events bare several hallmarks of self-organizing biologic systems that designate spatial information during pattern formation. Intracellular pattern formation represents a new paradigm in mitotic regulation.

Interpretation of the UPD/JAK/STAT morphogen gradient in Drosophila follicle cells

We are using Drosophila follicle cells to study the mechanisms that promote cell motility. Using genetics we identified a gene regulatory network that controls the dynamic pattern of activation of JAK/STAT in anterior follicle cells. Under the influence of a graded signal, Unpaired (UPD), JAK/STAT becomes activated first in a graded fashion. STAT, in turn, locally activates its own repressor, Apontic (APT), a new feedback regulator of JAK/STAT signaling. High levels of JAK/STAT also activate Slow Border Cells (SLBO), which undermines APT-mediated repression. In this way, cells that achieve a high JAK/STAT level maintain SLBO expression and form border cells, which then migrate out of the cell layer. Cells with lower JAK/STAT activity express more APT than SLBO, ultimately lose STAT activity, and remain in the follicular epithelium. To better understand how the graded signal is converted to an all-or-none decision to move or stay, we developed a mathematical model. Simulations using the model reproduce the observed dynamics of JAK/STAT expression in the wild type and in several mutant situations. By combining biological experiments and mathematical modeling, we can achieve a more sophisticated understanding of how cells interpret molecular gradients.

Dynamic gene expression is required for anterior regionalization in a spider

Patterning of a multicellular embryo requires precise spatiotemporal control of gene expression during development. The gradient of the morphogen bicoid regulates anterior regionalization in the syncytial blastoderm of Drosophila. However many arthropod embryos develop from a cellular blastoderm that does not allow the formation of transcription factor gradients. Here we show that correct anterior development of the cellularized embryo of the spider Achaearanea tepidariorum requires an anterior-to-posterior wave of dynamic gene expression for positioning the stripes of hairy, hedgehog, and orthodenticle expression. Surprisingly, this dynamic repositioning of the expression of these segmentation genes is blocked in orthodenticle(pRNAi) embryos and no anterior structures are specified in those embryos. Our data suggest that dynamic gene expression across a field of cells is required for anterior regionalization in spiders and provides an explanation for the problem of how positional values for anterior segmentation genes are specified via a morphogen-independent mechanism across a field of cells.

The loci of evolution: how predictable is genetic evolution?

Is genetic evolution predictable? Evolutionary developmental biologists have argued that, at least for morphological traits, the answer is a resounding yes. Most mutations causing morphological variation are expected to reside in the cis-regulatory, rather than the coding, regions of developmental genes. This "cis-regulatory hypothesis" has recently come under attack. In this review, we first describe and critique the arguments that have been proposed in support of the cis-regulatory hypothesis. We then test the empirical support for the cis-regulatory hypothesis with a comprehensive survey of mutations responsible for phenotypic evolution in multicellular organisms. Cis-regulatory mutations currently represent approximately 22% of 331 identified genetic changes although the number of cis-regulatory changes published annually is rapidly increasing. Above the species level, cis-regulatory mutations altering morphology are more common than coding changes. Also, above the species level cis-regulatory mutations predominate for genes not involved in terminal differentiation. These patterns imply that the simple question "Do coding or cis-regulatory mutations cause more phenotypic evolution?" hides more interesting phenomena. Evolution in different kinds of populations and over different durations may result in selection of different kinds of mutations. Predicting the genetic basis of evolution requires a comprehensive synthesis of molecular developmental biology and population genetics.