Identification and classification of cis-regulatory elements in the amphipod crustacean Parhyale hawaiensis.

Emerging research organisms enable the study of biology that cannot be addressed using classical 'model' organisms. New data resources can accelerate research in such animals. Here, we present new functional genomic resources for the amphipod crustacean Parhyale hawaiensis, facilitating the exploration of gene regulatory evolution using this emerging research organism. We use Omni-ATAC-seq to identify accessible chromatin genome-wide across a broad time course of Parhyale embryonic development. This time course encompasses many major morphological events, including segmentation, body regionalization, gut morphogenesis and limb development. In addition, we use short- and long-read RNA-seq to generate an improved Parhyale genome annotation, enabling deeper classification of identified regulatory elements. We discover differential accessibility, predict nucleosome positioning, infer transcription factor binding, cluster peaks based on accessibility dynamics, classify biological functions and correlate gene expression with accessibility. Using a Minos transposase reporter system, we demonstrate the potential to identify novel regulatory elements using this approach. This work provides a platform for the identification of novel developmental regulatory elements in Parhyale, and offers a framework for performing such experiments in other emerging research organisms.

Dual Functions of labial Resolve the Hox Logic of Chelicerate Head Segments.

Despite an abundance of gene expression surveys, comparatively little is known about Hox gene function in Chelicerata. Previous investigations of paralogs of labial (lab) and Deformed (Dfd) in a spider have shown that these play a role in tissue maintenance of the pedipalp segment (lab-1) and in patterning the first walking leg identity (Dfd-1), respectively. However, extrapolations of these data across chelicerates are hindered by the existence of duplicated Hox genes in arachnopulmonates (e.g., spiders and scorpions), which have resulted from an ancient whole genome duplication (WGD) event. Here, we investigated the function of the single-copy ortholog of lab in the harvestman Phalangium opilio, an exemplar of a lineage that was not subject to this WGD. Embryonic RNA interference against lab resulted in two classes of phenotypes: homeotic transformations of pedipalps to chelicerae, as well as reduction and fusion of the pedipalp and leg 1 segments. To test for combinatorial function, we performed a double knockdown of lab and Dfd, which resulted in a homeotic transformation of both pedipalps and the first walking legs into cheliceral identity, whereas the second walking leg is transformed into a pedipalpal identity. Taken together, these results elucidate a model for the Hox logic of head segments in Chelicerata. To substantiate the validity of this model, we performed expression surveys for lab and Dfd paralogs in scorpions and horseshoe crabs. We show that repetition of morphologically similar appendages is correlated with uniform expression levels of the Hox genes lab and Dfd, irrespective of the number of gene copies.

Expression of Abdominal-B in the brine shrimp, Artemia franciscana, expands our evolutionary understanding of the crustacean abdomen.

The brine shrimp, Artemia franciscana, has a body plan composed of 11 thoracic segments, followed by 2 genital segments, and then 6 additional abdominal segments. Previous studies of Artemia reported that expression of the posterior-most Hox gene, Abdominal-B (Abd-B), is restricted to the genital segments and is not observed posteriorly in the abdomen at any developmental stage. This report was remarkable because it suggested that the Artemia abdomen posterior to the genital segments was a novel body region of 6 segments that bore no homology to any region in other crustaceans and was unique amongst arthropods in being a Hox-free segmented domain outside of the head. In this study, we used RT-PCR, antibody staining, and in situ hybridization on various stages of Artemia nauplii to show that Abd-B mRNA and protein are in fact expressed throughout the abdominal segments during Artemia development, but this expression later retracts to the two genital segments (G1, G2) and the T11 appendages. This suggests that Abd-B does play a role in specifying abdominal segment identity in all crustaceans that have been examined and suggests a common evolutionary origin for the crustacean abdomen.

A meta-analysis of butterfly structural colors: their color range, distribution and biological production.

Butterfly scales are among the richest natural sources of optical nanostructures, which produce structural color and iridescence. Several recurring nanostructure types have been described, such as ridge multilayers, gyroids and lower lamina thin films. While the optical mechanisms of these nanostructure classes are known, their phylogenetic distributions and functional ranges have not been described in detail. In this Review, we examine a century of research on the biological production of structural colors, including their evolution, development and genetic regulation. We have also created a database of more than 300 optical nanostructures in butterflies and conducted a meta-analysis of the color range, abundance and phylogenetic distribution of each nanostructure class. Butterfly structural colors are ubiquitous in short wavelengths but extremely rare in long wavelengths, especially red. In particular, blue wavelengths (around 450 nm) occur in more clades and are produced by more kinds of nanostructures than other hues. Nanostructure categories differ in prevalence, phylogenetic distribution, color range and brightness. For example, lamina thin films are the least bright; perforated lumen multilayers occur most often but are almost entirely restricted to the family Lycaenidae; and 3D photonic crystals, including gyroids, have the narrowest wavelength range (from about 450 to 550 nm). We discuss the implications of these patterns in terms of nanostructure evolution, physical constraint and relationships to pigmentary color. Finally, we highlight opportunities for future research, such as analyses of subadult and Hesperid structural colors and the identification of genes that directly build the nanostructures, with relevance for biomimetic engineering.

Developmental, cellular and biochemical basis of transparency in clearwing butterflies.

The wings of butterflies and moths (Lepidoptera) are typically covered with thousands of flat, overlapping scales that endow the wings with colorful patterns. Yet, numerous species of Lepidoptera have evolved highly transparent wings, which often possess scales of altered morphology and reduced size, and the presence of membrane surface nanostructures that dramatically reduce reflection. Optical properties and anti-reflective nanostructures have been characterized for several 'clearwing' Lepidoptera, but the developmental processes underlying wing transparency are unknown. Here, we applied confocal and electron microscopy to create a developmental time series in the glasswing butterfly, Greta oto, comparing transparent and non-transparent wing regions. We found that during early wing development, scale precursor cell density was reduced in transparent regions, and cytoskeletal organization during scale growth differed between thin, bristle-like scale morphologies within transparent regions and flat, round scale morphologies within opaque regions. We also show that nanostructures on the wing membrane surface are composed of two layers: a lower layer of regularly arranged nipple-like nanostructures, and an upper layer of irregularly arranged wax-based nanopillars composed predominantly of long-chain n-alkanes. By chemically removing wax-based nanopillars, along with optical spectroscopy and analytical simulations, we demonstrate their role in generating anti-reflective properties. These findings provide insight into morphogenesis and composition of naturally organized microstructures and nanostructures, and may provide bioinspiration for new anti-reflective materials.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity.

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.

Unraveling the Tangled Skein: The Evolution of Transcriptional Regulatory Networks in Development.

The molecular and genetic basis for the evolution of anatomical diversity is a major question that has inspired evolutionary and developmental biologists for decades. Because morphology takes form during development, a true comprehension of how anatomical structures evolve requires an understanding of the evolutionary events that alter developmental genetic programs. Vast gene regulatory networks (GRNs) that connect transcription factors to their target regulatory sequences control gene expression in time and space and therefore determine the tissue-specific genetic programs that shape morphological structures. In recent years, many new examples have greatly advanced our understanding of the genetic alterations that modify GRNs to generate newly evolved morphologies. Here, we review several aspects of GRN evolution, including their deep preservation, their mechanisms of alteration, and how they originate to generate novel developmental programs.

Comprehensive analysis of Hox gene expression in the amphipod crustacean Parhyale hawaiensis.

Hox genes play crucial roles in establishing regional identity along the anterior-posterior axis in bilaterian animals, and have been implicated in generating morphological diversity throughout evolution. Here we report the identification, expression, and initial genomic characterization of the complete set of Hox genes from the amphipod crustacean Parhyale hawaiensis. Parhyale is an emerging model system that is amenable to experimental manipulations and evolutionary comparisons among the arthropods. Our analyses indicate that the Parhyale genome contains a single copy of each canonical Hox gene with the exception of fushi tarazu, and preliminary mapping suggests that at least some of these genes are clustered together in the genome. With few exceptions, Parhyale Hox genes exhibit both temporal and spatial colinearity, and expression boundaries correlate with morphological differences between segments and their associated appendages. This work represents the most comprehensive analysis of Hox gene expression in a crustacean to date, and provides a foundation for functional studies aimed at elucidating the role of Hox genes in arthropod development and evolution.

Dynamics of F-actin prefigure the structure of butterfly wing scales.

The wings of butterflies and moths consist of dorsal and ventral epidermal surfaces that give rise to overlapping layers of scales and hairs (Lepidoptera, "scale wing"). Wing scales (average length ~200 µm) are homologous to insect bristles (macrochaetes), and their colors create the patterns that characterize lepidopteran wings. The topology and surface sculpture of wing scales vary widely, and this architectural complexity arises from variations in the developmental program of the individual scale cells of the wing epithelium. One of the more striking features of lepidopteran wing scales are the longitudinal ridges that run the length of the mature (dead) cell, gathering the cuticularized scale cell surface into pleats on the sides of each scale. While also present around the periphery of other insect bristles and hairs, longitudinal ridges in lepidopteran wing scales gain new significance for their creation of iridescent color through microribs and lamellae. Here we show the dynamics of the highly organized F-actin filaments during scale cell development, and present experimental manipulations of actin polymerization that reveal the essential role of this cytoskeletal component in wing scale elongation and the positioning of longitudinal ribs.

Nodal signalling is involved in left-right asymmetry in snails.

Many animals display specific internal or external features with left-right asymmetry. In vertebrates, the molecular pathway that leads to this asymmetry uses the signalling molecule Nodal, a member of the transforming growth factor-beta superfamily, which is expressed in the left lateral plate mesoderm, and loss of nodal function produces a randomization of the left-right asymmetry of visceral organs. Orthologues of nodal have also been described in other deuterostomes, including ascidians and sea urchins, but no nodal orthologue has been reported in the other two main clades of Bilateria: Ecdysozoa (including flies and nematodes) and Lophotrochozoa (including snails and annelids). Here we report the first evidence for a nodal orthologue in a non-deuterostome group. We isolated nodal and Pitx (one of the targets of Nodal signalling) in two species of snails and found that the side of the embryo that expresses nodal and Pitx is related to body chirality: both genes are expressed on the right side of the embryo in the dextral (right-handed) species Lottia gigantea and on the left side in the sinistral (left-handed) species Biomphalaria glabrata. We pharmacologically inhibited the Nodal pathway and found that nodal acts upstream of Pitx, and that some treated animals developed with a loss of shell chirality. These results indicate that the involvement of the Nodal pathway in left-right asymmetry might have been an ancestral feature of the Bilateria.

Evolving specialization of the arthropod nervous system.

The diverse array of body plans possessed by arthropods is created by generating variations upon a design of repeated segments formed during development, using a relatively small "toolbox" of conserved patterning genes. These attributes make the arthropod body plan a valuable model for elucidating how changes in development create diversity of form. As increasingly specialized segments and appendages evolved in arthropods, the nervous systems of these animals also evolved to control the function of these structures. Although there is a remarkable degree of conservation in neural development both between individual segments in any given species and between the nervous systems of different arthropod groups, the differences that do exist are informative for inferring general principles about the holistic evolution of body plans. This review describes developmental processes controlling neural segmentation and regionalization, highlighting segmentation mechanisms that create both ectodermal and neural segments, as well as recent studies of the role of Hox genes in generating regional specification within the central nervous system. We argue that this system generates a modular design that allows the nervous system to evolve in concert with the body segments and their associated appendages. This information will be useful in future studies of macroevolutionary changes in arthropod body plans, especially in understanding how these transformations can be made in a way that retains the function of appendages during evolutionary transitions in morphology.

Using phylogenetically-informed annotation (PIA) to search for light-interacting genes in transcriptomes from non-model organisms.

BACKGROUND: Tools for high throughput sequencing and de novo assembly make the analysis of transcriptomes (i.e. the suite of genes expressed in a tissue) feasible for almost any organism. Yet a challenge for biologists is that it can be difficult to assign identities to gene sequences, especially from non-model organisms. Phylogenetic analyses are one useful method for assigning identities to these sequences, but such methods tend to be time-consuming because of the need to re-calculate trees for every gene of interest and each time a new data set is analyzed. In response, we employed existing tools for phylogenetic analysis to produce a computationally efficient, tree-based approach for annotating transcriptomes or new genomes that we term Phylogenetically-Informed Annotation (PIA), which places uncharacterized genes into pre-calculated phylogenies of gene families. RESULTS: We generated maximum likelihood trees for 109 genes from a Light Interaction Toolkit (LIT), a collection of genes that underlie the function or development of light-interacting structures in metazoans. To do so, we searched protein sequences predicted from 29 fully-sequenced genomes and built trees using tools for phylogenetic analysis in the Osiris package of Galaxy (an open-source workflow management system). Next, to rapidly annotate transcriptomes from organisms that lack sequenced genomes, we repurposed a maximum likelihood-based Evolutionary Placement Algorithm (implemented in RAxML) to place sequences of potential LIT genes on to our pre-calculated gene trees. Finally, we implemented PIA in Galaxy and used it to search for LIT genes in 28 newly-sequenced transcriptomes from the light-interacting tissues of a range of cephalopod mollusks, arthropods, and cubozoan cnidarians. Our new trees for LIT genes are available on the Bitbucket public repository ( http://bitbucket.org/osiris_phylogenetics/pia/ ) and we demonstrate PIA on a publicly-accessible web server ( http://galaxy-dev.cnsi.ucsb.edu/pia/ ). CONCLUSIONS: Our new trees for LIT genes will be a valuable resource for researchers studying the evolution of eyes or other light-interacting structures. We also introduce PIA, a high throughput method for using phylogenetic relationships to identify LIT genes in transcriptomes from non-model organisms. With simple modifications, our methods may be used to search for different sets of genes or to annotate data sets from taxa outside of Metazoa.

Genetic basis of eye and pigment loss in the cave crustacean, Asellus aquaticus.

Understanding the process of evolution is one of the great challenges in biology. Cave animals are one group with immense potential to address the mechanisms of evolutionary change. Amazingly, similar morphological alterations, such as enhancement of sensory systems and the loss of eyes and pigmentation, have evolved multiple times in a diverse assemblage of cave animals. Our goal is to develop an invertebrate model to study cave evolution so that, in combination with a previously established vertebrate cave system, we can address genetic questions concerning evolutionary parallelism and convergence. We chose the isopod crustacean, Asellus aquaticus, and generated a genome-wide linkage map for this species. Our map, composed of 117 markers, of which the majority are associated with genes known to be involved in pigmentation, eye, and appendage development, was used to identify loci of large effect responsible for several pigmentation traits and eye loss. Our study provides support for the prediction that significant morphological change can be mediated through one or a few genes. Surprisingly, we found that within population variability in eye size occurs through multiple mechanisms; eye loss has a different genetic basis than reduced eye size. Similarly, again within a population, the phenotype of albinism can be achieved by two different genetic pathways--either by a recessive genotype at one locus or doubly recessive genotypes at two other loci. Our work shows the potential of Asellus for studying the extremes of parallel and convergent evolution-spanning comparisons within populations to comparisons between vertebrate and arthropod systems.

A novel expression domain of extradenticle underlies the evolutionary developmental origin of the chelicerate patella.

Neofunctionalization of duplicated gene copies is thought to be an important process underlying the origin of evolutionary novelty and provides an elegant mechanism for the origin of new phenotypic traits. One putative case where a new gene copy has been linked to a novel morphological trait is the origin of the arachnid patella, a taxonomically restricted leg segment. In spiders, the origin of this segment has been linked to the origin of the paralog dachshund-2 , suggesting that a new gene facilitated the expression of a new trait. However, various arachnid groups that possess patellae do not have a copy of dachshund-2 , disfavoring the direct link between gene origin and trait origin. We investigated the developmental genetic basis for patellar patterning in the harvestman Phalangium opilio , which lacks dachshund-2 . Here, we show that the harvestman patella is established by a novel expression domain of the transcription factor extradenticle . Leveraging this definition of patellar identity, we surveyed targeted groups across chelicerate phylogeny to assess when this trait evolved. We show that a patellar homolog is present in Pycnogonida (sea spiders) and various arachnid orders, suggesting a single origin of the patella in the ancestor of Chelicerata. A potential loss of the patella is observed in Ixodida. Our results suggest that the modification of an ancient gene, rather than the neofunctionalization of a new gene copy, underlies the origin of the patella. Broadly, this work underscores the value of comparative data and broad taxonomic sampling when testing hypotheses in evolutionary developmental biology.

Independent migration of cell populations in the early gastrulation of the amphipod crustacean Parhyale hawaiensis.

Cells are the principal component of tissues and can drive morphogenesis through dynamic changes in structure and interaction. During gastrulation, the primary morphogenetic event of early development, cells change shape, exchange neighbors, and migrate long distances to establish cell layers that will form the tissues of the adult animal. Outside of Drosophila, little is known about how changes in cell behavior might drive gastrulation among arthropods. Here, we focus on three cell populations that form two aggregations during early gastrulation in the crustacean Parhyale hawaiensis. Using cytoskeletal markers and lineage tracing we observe bottle cells in anterior and visceral mesoderm precursors as gastrulation commences, and find that both Cytochalasin D, an inhibitor of actin polymerization, and ROCKOUT, an inhibitor of Rho-kinase activity, prevent gastrulation. Furthermore, by ablating specific cells, we show that each of the three populations acts independently during gastrulation, confirming previous hypotheses that cell behavior during Parhyale gastrulation relies on intrinsic signals instead of an inductive mechanism.

The functional relationship between ectodermal and mesodermal segmentation in the crustacean, Parhyale hawaiensis.

In arthropods, annelids and chordates, segmentation of the body axis encompasses both ectodermal and mesodermal derivatives. In vertebrates, trunk mesoderm segments autonomously and induces segmental arrangement of the ectoderm-derived nervous system. In contrast, in the arthropod Drosophila melanogaster, the ectoderm segments autonomously and mesoderm segmentation is at least partially dependent on the ectoderm. While segmentation has been proposed to be a feature of the common ancestor of vertebrates and arthropods, considering vertebrates and Drosophila alone, it is impossible to conclude whether the ancestral primary segmented tissue was the ectoderm or the mesoderm. Furthermore, much of Drosophila segmentation occurs before gastrulation and thus may not accurately represent the mechanisms of segmentation in all arthropods. To better understand the relationship between segmented germ layers in arthropods, we asked whether segmentation is an intrinsic property of the ectoderm and/or the mesoderm in the crustacean Parhyale hawaiensis by ablating either the ectoderm or the mesoderm and then assaying for segmentation in the remaining tissue layer. We found that the ectoderm segments autonomously. However, mesoderm segmentation requires at least a permissive signal from the ectoderm. Although mesodermal stem cells undergo normal rounds of division in the absence of ectoderm, they do not migrate properly in respect to migration direction and distance. In addition, their progeny neither divide nor express the mesoderm segmentation markers Ph-twist and Ph-Even-skipped. As segmentation is ectoderm-dependent in both Parhyale and holometabola insects, we hypothesize that segmentation is primarily a property of the ectoderm in pancrustacea.

giant is a bona fide gap gene in the intermediate germband insect, Oncopeltus fasciatus.

Drosophila undergoes a form of development termed long germ segmentation, where all segments are specified nearly simultaneously so that by the blastoderm stage, the entire body plan has been determined. This mode of segmentation is evolutionarily derived. Most insects undergo short or intermediate germ segmentation, where only anterior segments are specified early, and posterior segments are sequentially specified during germband elongation. These embryological differences imply that anterior and posterior segments might rely upon different molecular mechanisms. In Drosophila, embryos mutant for giant show a gap in the anterior as well fusions of several abdominal segments. In Tribolium, a short germ beetle, giant is required for segmental identity, but not formation, in gnathal segments and also for segmentation of the entire abdomen. This raises the possibility that giant might not act as a gap gene in short and intermediate germ insects. Oncopeltus fasciatus is an intermediate germ insect that is an outgroup to the clade containing Drosophila and Tribolium. We cloned the Oncopeltus homolog of giant and determined its expression and function during segmentation. We find that Oncopeltus giant is a canonical gap gene in the maxillary and labial segments and also plays a gap-like role in the first four abdominal segments. Our results suggest that giant was a bona fide gap gene in the ancestor of these insects with this role being lost in the lineage leading towards Tribolium. This highlights the conservation of anterior patterning and evolutionary plasticity of the genetic regulation controlling posterior segmentation, even in short and intermediate germ insects.

A prominent requirement for single-minded and the ventral midline in patterning the dorsoventral axis of the crustacean Parhyale hawaiensis.

In bilaterians, establishing the correct spatial positioning of structures along the dorsoventral (DV) axis is essential for proper embryonic development. Insects such as Drosophila rely on the Dorsal activity gradient and Bone morphogenetic protein (BMP) signaling to establish cell fates along the DV axis, leading to the distinction between tissues such as mesoderm, neurogenic ectoderm and dorsal ectoderm in the developing embryo. Subsequently, the ventral midline plays a more restricted role in DV patterning by establishing differential cell fates in adjacent regions of the neurogenic ectoderm. In this study, we examine the function of the ventral midline and the midline-associated gene single-minded (Ph-sim) in the amphipod crustacean Parhyale hawaiensis. Remarkably, we found that Ph-sim and the ventral midline play a central role in establishing proper fates along the entire DV axis in this animal; laser ablation of midline cells causes a failure to form neurogenic ectoderm and Ph-sim RNAi results in severely dorsalized embryos lacking both neurogenic ectoderm and the appendage-bearing lateral ectoderm. Furthermore, we hypothesize that this role of midline cells was present in the last common ancestor of crustaceans and insects. We predict that the transition to a Dorsal-dependent DV patterning system in the phylogenetically derived insect lineage leading to Drosophila has led to a more restricted role of the ventral midline in patterning the DV axis of these insects.

Discovery of the bicycle gene family provides new insights into insect manipulation of plant development during gall induction.

Galls are complex structures that develop from plant tissue, providing protection and food for gall-forming organisms, such as insects or mites. However, the molecules used by insects or mites to manipulate plant development have proved elusive. A landmark study has tracked down a gene in a gall-forming aphid that controls whether galls on witch hazel are green or red. The 'green allele' is strongly expressed in aphid salivary glands and represses plant genes used for red color formation. Excitingly, the gene product is part of a large suite of proteins that aphids may use to interact with plant biology.

Analysis of snail genes in the crustacean Parhyale hawaiensis: insight into snail gene family evolution.

The transcriptional repressor snail was first discovered in Drosophila melanogaster, where it initially plays a role in gastrulation and mesoderm formation, and later plays a role in neurogenesis. Among arthropods, this role of snail appears to be conserved in the insects Tribolium and Anopheles gambiae, but not in the chelicerates Cupiennius salei and Achaearanea tepidariorum, the myriapod Glomeris marginata, or the Branchiopod crustacean Daphnia magna. These data imply that within arthropoda, snail acquired its role in gastrulation and mesoderm formation in the insect lineage. However, crustaceans are a diverse group with several major taxa, making analysis of more crustaceans necessary to potentially understand the ancestral role of snail in Pancrustacea (crustaceans + insects) and thus in the ancestor of insects as well. To address these questions, we examined the snail family in the Malacostracan crustacean Parhyale hawaiensis. We found three snail homologs, Ph-snail1, Ph-snail2 and Ph-snail3, and one scratch homolog, Ph-scratch. Parhyale snail genes are expressed after gastrulation, during germband formation and elongation. Ph-snail1, Ph-snail2, and Ph-snail3 are expressed in distinct patterns in the neuroectoderm. Ph-snail1 is the only Parhyale snail gene expressed in the mesoderm, where its expression cycles in the mesodermal stem cells, called mesoteloblasts. The mesoteloblasts go through a series of cycles, where each cycle is composed of a migration phase and a division phase. Ph-snail1 is expressed during the migration phase, but not during the division phase. We found that as each mesoteloblast division produces one segment's worth of mesoderm, Ph-snail1 expression is linked to both the cell cycle and the segmental production of mesoderm.