Abstract: Understanding the general trends in developmental changes during animal evolution, which are often associated with morphological diversification, has long been a central issue in evolutionary developmental biology. Recent comparative transcriptomic studies revealed that gene expression profiles of mid-embryonic period tend to be more evolutionarily conserved than those in earlier or later periods. While the hourglass-like divergence of developmental processes has been demonstrated in a variety of animal groups such as vertebrates, arthropods, and nematodes, the exact mechanism leading to this mid-embryonic conservation remains to be clarified. One possibility is that the mid-embryonic period (pharyngula period in vertebrates) is highly prone to embryonic lethality, and the resulting negative selections lead to evolutionary conservation of this phase. Here, we tested this “mid-embryonic lethality hypothesis” by measuring the rate of lethal phenotypes of three different species of vertebrate embryos subjected to two kinds of perturbations: transient perturbations and genetic mutations. By subjecting zebrafish (Danio rerio), African clawed frog (Xenopus laevis), and chicken (Gallus gallus) embryos to transient perturbations, namely heat shock and inhibitor treatments during three developmental periods [early (represented by blastula and gastrula), pharyngula, and late], we found that the early stages showed the highest rate of lethal phenotypes in all three species. This result was corroborated by perturbation with genetic mutations. By tracking the survival rate of wild-type embryos and embryos with genetic mutations induced by UV irradiation in zebrafish and African clawed frogs, we found that the highest decrease in survival rate was at the early stages particularly around gastrulation in both these species. In opposition to the “mid-embryonic lethality hypothesis,” our results consistently showed that the stage with the highest lethality was not around the conserved pharyngula period, but rather around the early period in all the vertebrate species tested. These results suggest that negative selection by embryonic lethality could not explain hourglass-like conservation of animal embryos. This highlights the potential contribution of alternative mechanisms such as the diversifying effect of positive selections against earlier and later stages, and developmental constraints which lead to conservation of mid-embryonic stages.

Abstract: Egg size represents an important form of maternal effect determined by a complex interplay of long-term adaptation and short-term plasticity balancing egg size with brood size. Haplochromine cichlids are maternal mouthbrooders showing differential parental investment in different species, manifested in great variation in egg size, brood size and duration of maternal care. Little is known about maternally determined molecular characters of eggs in fishes and their relation to egg size and trophic specialization. Here we investigate maternal mRNA inputs of selected growth- and stress-related genes in eggs of mouthbrooding cichlid fishes adapted to different trophic niches from Lake Tanganyika, Lake Malawi, Lake Victoria and compare them to their riverine allies. We first identified two reference genes, atf7ip and mid1ip1, to be suitable for cross-species quantification of mRNA abundance via qRT-PCR in the cichlid eggs. Using these reference genes, we found substantial variation in maternal mRNA input for a set of candidate genes related to growth and stress response across species and lakes. We observed negative correlation of mRNA abundance between two of growth hormone receptor paralogs (ghr1 and ghr2) across all haplochromine cichlid species which also differentiate the species in the two younger lakes, Malawi and Lake Victoria, from those in Lake Tanganyika and ancestral riverine species. Furthermore, we found correlations between egg size and maternal mRNA abundance of two growth-related genes igf2 and ghr2 across the haplochromine cichlids as well as distinct clustering of the species based on their trophic specialization using maternal mRNA abundance of five genes (ghr1, ghr2, igf2, gr and sgk1). These findings indicate that variations in egg size in closely related cichlid species can be linked to differences in maternal RNA deposition of key growth-related genes. In addition, the cichlid species with contrasting trophic specialization deposit different levels of maternal mRNAs in their eggs for particular growth-related genes; however, it is unclear whether such differences contribute to differential morphogenesis at later stages of development. Our results provide first insights into this aspect of gene activation, as a basis for future studies targeting their role during ecomorphological specialization and adaptive radiation.

Abstract: Understanding the evolution and development of morphological traits of the last common bilaterian ancestor is a major goal of the evo-devo discipline. The reconstruction of this “urbilaterian” is mainly based on comparative studies of common molecular patterning mechanisms in recent model organisms. The NK homeobox genes are key players in many of these molecular pathways, including processes regulating mesoderm, heart and neural development. Shared features seen in the expression patterns of NK genes have been used to determine the ancestral bilaterian characters. However, the commonly used model organisms provide only a limited view on the evolution of these molecular pathways. To further investigate the ancestral roles of NK cluster genes, we analyzed their expression patterns in the onychophoran Euperipatoides rowelli. We identified nine transcripts of NK cluster genes in E. rowelli, including single copies of NK1, NK3, NK4, NK5, Msx, Lbx and Tlx, and two copies of NK6. All of these genes except for NK6.1 and NK6.2 are expressed in different mesodermal organs and tissues in embryos of E. rowelli, including the anlagen of somatic musculature and the heart. Furthermore, we found distinct expression patterns of NK3, NK5, NK6, Lbx and Msx in the developing nervous system. The same holds true for the NKL gene NK2.2, which does not belong to the NK cluster but is a related gene playing a role in neural patterning. Surprisingly, NK1, Msx and Lbx are additionally expressed in a segment polarity-like pattern early in development—a feature that has been otherwise reported only from annelids. Our results indicate that the NK cluster genes were involved in mesoderm and neural development in the last common ancestor of bilaterians or at least nephrozoans (i.e., bilaterians to the exclusion of xenacoelomorphs). By comparing our data from an onychophoran to those from other bilaterians, we critically review the hypothesis of a complex “urbilaterian” with a segmented body, a pulsatile organ or heart, and a condensed mediolaterally patterned nerve cord.

Abstract: For early larvae of amphioxus, Kaji et al. (Zool Lett 2:2, 2016) proposed that mesoderm cells are added to the rim of the forming mouth, giving it the quality of a coelomoduct without homology to the oral openings of other animals. They depended in part on non-serial transmission electron microscopic (TEM) sections and could not readily put fine structural details into a broader context. The present study of amphioxus larvae is based largely on serial blockface scanning electron microscopy (SBSEM), a technique revealing TEM-level details within an extensive anatomical volume that can be reconstructed in three dimensions. In amphioxus larvae shortly before mouth formation, a population of compact mesoderm cells is present at the posterior extremity of the first left somite. As development continues, the more dorsal of these cells give rise to the initial kidney (Hatschek’s nephridium), while the more ventral cells become interposed between the ectoderm and endoderm in a localized region where the mouth will soon penetrate. SBSEM reveals that, after the mouth has opened, a majority of these mesoderm cells can still be detected, sandwiched between the ectoderm and endoderm; they are probably myoblasts destined to develop into the perioral muscles. SBSEM has provided the most accurate and detailed description to date of the tissues at the anterior end of amphioxus larvae. The present study supports the finding of Kaji et al. (2016) that the more dorsal of the cells in the posterior region of the first left somite give rise to the initial kidney. In contrast, the fate of the more ventral cells (called here the oral mesoderm) is less well understood. Although Kaji et al. (2016) implied that all of the oral mesoderm cells joined the rim of the forming mouth, SBSEM reveals that many of them are still present after mouth penetration. Even so, some of those cells go missing during mouth penetration and their fate is unknown. It cannot be ruled out that they were incorporated into the rim of the nascent mouth as proposed by Kaji et al. (2016). On the other hand, they might have degenerated or been shed from the larva during the morphogenetic interaction between the ectoderm and endoderm to form the mouth. The present SBSEM study, like Kaji et al. (2016), is based on static morphological data, and dynamic cell tracer experiments would be needed to decide among these possibilities.

Abstract: The axonal projection from the retina to the optic tectum maps visual information isomorphically from one to the other and serves as a model for the development of sensory maps more generally in the vertebrate brain. How or why this connection evolved is not known, nor why the midbrain is so important to the processing of visual information. Amphioxus is potentially informative here because its eye homolog, the frontal eye, also has a neural connection to a region of the brain now known to be homologous with the caudal diencephalon and midbrain. The frontal eye has only a one-dimensional receptor array, but simple alterations to the pattern and plane of cell division would have been sufficient to generate a structure more like the vertebrate retina. Accounting for the retinotectal map poses more of a problem. The hypothesis developed here is that this is best explained as a consequence of a prior association between the roof of the anterior nerve cord and an array of rhabdomeric photoreceptors, homologous with the Joseph cells of amphioxus, that were used by the common ancestor of amphioxus and vertebrates for detecting moving shadows. Hence, a rudimentary tectal map could have been present before the evolution of image-forming eyes and been coopted by them secondarily. Assuming the orientation of this map was fixed from the start relative to the external world, its retinal counterpart would have had to adjust to this to accommodate the image reversal that accompanies the conversion of a flat receptor array to a camera-type eye. Exploring this hypothesis further will require more information than is currently available on the Joseph cells, especially as to where and how their neural output is processed.

Abstract: The larval nervous system of the solitary tunicate Ciona is a simple model for the study of chordate neurodevelopment. The development and connectivity of the Ciona motor ganglion have been studied in fine detail, but how this important structure develops in other tunicates is not well known. By comparing gene expression patterns in the developing MG of the distantly related tunicate Molgula occidentalis, we found that its patterning is highly conserved compared to the Ciona MG. MG neuronal subtypes in Molgula were specified in the exact same positions as in Ciona, though the timing of subtype-specific gene expression onset was slightly shifted to begin earlier, relative to mitotic exit and differentiation. In transgenic Molgula embryos electroporated with Dmbx reporter plasmids, we were also able to characterize the morphology of the lone pair of descending decussating neurons (ddNs) in Molgula, revealing the same unique contralateral projection seen in Ciona ddNs and their putative vertebrate homologs the Mauthner cells. Although Dmbx expression labels the ddNs in both species, cross-species transgenic assays revealed significant changes to the regulatory logic underlying Dmbx transcription. We found that Dmbx cis-regulatory DNAs from Ciona can drive highly specific reporter gene expression in Molgula ddNs, but Molgula sequences are not active in Ciona ddNs. This acute divergence in the molecular mechanisms that underlie otherwise functionally conserved cis-regulatory DNAs supports the recently proposed idea that the extreme genetic plasticity observed in tunicates may be attributed to the extreme rigidity of the spatial organization of their embryonic cell lineages.

Abstract: Current studies in evolutionary developmental biology are focused on the reconstruction of gene regulatory networks in target animal species. From decades, the scientific interest on genetic mechanisms orchestrating embryos development has been increasing in consequence to the fact that common features shared by evolutionarily distant phyla are being clarified. In 2011, a study across eumetazoan species showed for the first time the existence of a highly conserved non-coding element controlling the SoxB2 gene, which is involved in the early specification of the nervous system. This discovery raised several questions about SoxB2 function and regulation in deuterostomes from an evolutionary point of view. Due to the relevant phylogenetic position within deuterostomes, the sea urchin Strongylocentrotus purpuratus represents an advantageous animal model in the field of evolutionary developmental biology. Herein, we show a comprehensive study of SoxB2 functions in sea urchins, in particular its expression pattern in a wide range of developmental stages, and its co-localization with other neurogenic markers, as SoxB1, SoxC and Elav. Moreover, this work provides a detailed description of the phenotype of sea urchin SoxB2 knocked-down embryos, confirming its key function in neurogenesis and revealing, for the first time, its additional roles in oral and aboral ectoderm cilia and skeletal rod morphology. We concluded that SoxB2 in sea urchins has a neurogenic function; however, this gene could have multiple roles in sea urchin embryogenesis, expanding its expression in non-neurogenic cells. We showed that SoxB2 is functionally conserved among deuterostomes and suggested that in S. purpuratus this gene acquired additional functions, being involved in ciliogenesis and skeletal patterning.

Abstract: Although retinoic acid (RA) signaling plays a crucial role in the body patterning of chordates, its function in non-chordate invertebrates, other than its mediation of environmental cues triggering metamorphosis in cnidarians, is largely unknown. We investigated the role of RA signaling in the metamorphosis of starfish (Echinodermata). We found that exogenous RA treatment induced metamorphosis in starfish larvae. In contrast, inhibitors of RA synthesis and RA receptors suppressed metamorphosis triggered by attachment to a substrate. Gene expressions of the RA signaling component were detected in competent larvae. This study provides insight into the ancestral function of RA signaling, which is conserved in the metamorphosis of cnidarians and starfish.

Abstract: How genome complexity affects organismal phenotypic complexity is a fundamental question in evolutionary developmental biology. Previous studies proposed various contributing factors of genome complexity and tried to find the connection between genomic complexity and organism complexity. However, a general model to answer this question is lacking. Here, we introduce a ‘two-level’ model for the realization of genome complexity at phenotypic level. Five representative species across Protostomia and Deuterostomia were involved in this study. The intrinsic gene properties contributing to genome complexity were classified into two generalized groups: the complexity and age degree of both protein-coding and noncoding genes. We found that young genes tend to be simpler; however, the mid-age genes, rather than the oldest genes, show the highest proportion of high complexity. Complex genes tend to be utilized preferentially in each stage of embryonic development, with maximum representation during the late stage of organogenesis. This trend is mainly attributed to mid-age complex genes. In contrast, young genes tend to be expressed in specific spatiotemporal states. An obvious correlation between the time point of the change in over- and under-representation and the order of gene age was observed, which supports the funnel-like model of the conservation pattern of development. In addition, we found some probable causes for the seemingly contradictory ‘funnel-like’ or ‘hourglass’ model. These results indicate that complex and young genes contribute to organismal complexity at two different levels: Complex genes contribute to the complexity of individual proteomes in certain states, whereas young genes contribute to the diversity of proteomes in different spatiotemporal states. This conclusion is valid across the five species investigated, indicating it is a conserved model across Protostomia and Deuterostomia. The results in this study also support ‘funnel-like model’ from a new viewpoint and explain why there are different evo–devo relation models.

Abstract: The concerted activity of Meis and Hoxa11 transcription factors is essential for the subdivision of tetrapod limbs into proximo-distal (PD) domains; however, little is know about the evolution of this patterning mechanism. Here, we aim to study the expression of meis and hoxa11 orthologues in the median and paired rayed fins of zebrafish and in the lobed fins of the Australian lungfish. First, a late phase of expression of meis1.1 and hoxa11b in zebrafish dorsal and anal fins relates with segmentation of endochondral elements in proximal and distal radials. Second, our zebrafish in situ hybridization results reveal spatial and temporal changes between pectoral and pelvic fins. Third, in situ analysis of meis1, meis3 and hoxa11 genes in Neoceratodus pectoral fins identifies decoupled domains of expression along the PD axis. Our data raise the possibility that the origin of stylopod and zeugopod lies much deeper in gnathostome evolution and that variation in meis and hoxa11 expression has played a substantial role in the transformation of appendage anatomy. Moreover, these observations provide evidence that the Meis/Hoxa11 profile considered a hallmark of stylopod/zeugopod patterning is present in Neoceratodus.

Abstract: Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology. We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2. A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.

Abstract: The basic helix-loop-helix (bHLH) family of transcription factors is one of the largest superfamilies of regulatory transcription factors and is widely used in eukaryotic organisms. They play an essential role in a range of metabolic, physiological, and developmental processes, including the development of the nervous system (NS). These transcription factors have been studied in many metazoans, especially in vertebrates but also in early branching metazoan clades such as the cnidarians and sponges. However, currently very little is known about their expression in the most basally branching bilaterian group, the xenacoelomorphs. Recently, our laboratory has characterized the full complement of bHLH in the genome of two members of the Xenacoelomorpha, the xenoturbellid Xenoturbella bocki and the acoel Symsagittifera roscoffensis. Understanding the patterns of bHLH gene expression in members of this phylum (in space and time) provides critical new insights into the conserved roles of the bHLH and their putative specificities in this group. Our focus is on deciphering the specific roles that these genes have in the process of neurogenesis. Here, we analyze the developmental expression of the whole complement of bHLH genes identified in the acoel S. roscoffensis. Based on their expression patterns, several members of bHLH class A appear to have specific conserved roles in neurogenesis, while other class A genes (as well as members of other classes) have likely taken on more generalized functions. All gene expression patterns are described in embryos and early juveniles. Our results suggest that the main roles of the bHLH genes of S. roscoffensis are evolutionarily conserved, with a specific subset dedicated to patterning the nervous system: SrAscA, SrAscB, SrHes/Hey, SrNscl, SrSrebp, SrE12/E47 and SrOlig.

Abstract: The notochord has organizer properties and is required for floor plate induction and dorsoventral patterning of the neural tube This activity has been attributed to sonic hedgehog (shh) signaling, which originates in the notochord, forms a gradient, and autoinduces shh expression in the floor plate However, reported data are inconsistent and the spatiotemporal development of the relevant shh expression domains has not been studied in detail We therefore studied the expression dynamics of shh in rabbit, chicken and Xenopus laevis embryos (as well as indian hedgehog and desert hedgehog as possible alternative functional candidates in the chicken) Our analysis reveals a markedly divergent pattern within these vertebrates: whereas in the rabbit shh is first expressed in the notochord and its floor plate domain is then induced during subsequent somitogenesis stages, in the chick embryo shh is expressed in the prospective neuroectoderm prior to the notochord formation and, interestingly, prior to mesoderm immigration Neither indian hedgehog nor desert hedgehog are expressed in these midline structures although mRNA of both genes was detected in other structures of the early chick embryo In X laevis, shh is expressed at the beginning of gastrulation in a distinct area dorsal to the dorsal blastopore lip and adjacent to the prospective neuroectoderm, whereas the floor plate expresses shh at the end of gastrulation While shh expression patterns in rabbit and X laevis embryos are roughly compatible with the classical view of “ventral to dorsal induction” of the floor plate, the early shh expression in the chick floor plate challenges this model Intriguingly, this alternative sequence of domain induction is related to the asymmetrical morphogenesis of the primitive node and other axial organs in the chick Our results indicate that the floor plate in X laevis and chick embryos may be initially induced by planar interaction within the ectoderm or epiblast Furthermore, we propose that the mode of the floor plate induction adapts to the variant topography of interacting tissues during gastrulation and notochord formation and thereby reveals evolutionary plasticity of early embryonic induction

Abstract: Segmentation, the subdivision of the major body axis into repeated elements, is considered one of the major evolutionary innovations in bilaterian animals. In all three segmented animal clades, the predominant segmentation mechanism is sequential segmentation, where segments are generated one by one in anterior–posterior order from a posterior undifferentiated zone. In vertebrates and arthropods, sequential segmentation is thought to arise from a clock-and-wavefront-type mechanism, where oscillations in the posterior growth zone are transformed into a segmental prepattern in the anterior by a receding wavefront. Previous evo-devo simulation studies have demonstrated that this segmentation type repeatedly arises, supporting the idea of parallel evolutionary origins in these animal clades. Sequential segmentation has been studied most extensively in vertebrates, where travelling waves have been observed that reflect the slowing down of oscillations prior to their cessation and where these oscillations involve a highly complex regulatory network. It is currently unclear under which conditions this oscillator complexity and slowing should be expected to evolve, how they are related and to what extent similar properties should be expected for sequential segmentation in other animal species. To investigate these questions, we extend a previously developed computational model for the evolution of segmentation. We vary the slope of the posterior morphogen gradient and the strength of gene expression noise. We find that compared to a shallow gradient, a steep morphogen gradient allows for faster evolution and evolved oscillator networks are simpler. Furthermore, under steep gradients, damped oscillators often evolve, whereas shallow gradients appear to require persistent oscillators which are regularly accompanied by travelling waves, indicative of a frequency gradient. We show that gene expression noise increases the likelihood of evolving persistent oscillators under steep gradients and of evolving frequency gradients under shallow gradients. Surprisingly, we find that the evolutions of oscillator complexity and travelling waves are not correlated, suggesting that these properties may have evolved separately. Based on our findings, we suggest that travelling waves may have evolved in response to shallow morphogen gradients and gene expression noise. These two factors may thus also be responsible for the observed differences between different species within both the arthropod and chordate phyla.

Abstract: Sall (Spalt-like) proteins are zinc-finger transcription factors involved in a number of biological processes. They have only been studied in a few model organisms, such as Drosophila melanogaster, Caenorhabditis elegans, Schmidtea mediterranea and some vertebrates. Further taxon sampling is critical to understand the evolution and diversification of this protein and its functional roles in animals. Using genome and transcriptome mining, we confirmed the presence of sall genes in a range of additional animal taxa, for which their presence had not yet been described. We show that sall genes are broadly conserved across the Bilateria, and likely appeared in the bilaterian stem lineage. Our analysis of the protein domains shows that the characteristic arrangement of the multiple zinc-finger domains is conserved in bilaterians and may represent the ancient arrangement of this family of transcription factors. We also show the existence of a previously unknown zinc-finger domain. In situ hybridization was used to describe the gene expression patterns in embryonic and larval stages in two species of snails: Crepidula fornicata and Lottia gigantea. In L. gigantea, sall presents maternal expression, although later on the expression is restricted to the A and B quadrants during gastrulation and larval stage. In C. fornicata, sall has no maternal expression and it is expressed mainly in the A, C and D quadrants during blastula stages and in an asymmetric fashion during the larval stage. Our results suggest that the bilaterian common ancestor had a Sall protein with at least six zinc-finger domains. The evolution of Sall proteins in bilaterians might have occurred mostly as a result of the loss of protein domains and gene duplications leading to diversification. The new evidence complements previous studies in highlighting an important role of Sall proteins in bilaterian development. Our results show maternal expression of sall in the snail L. gigantea, but not C. fornicata. The asymmetric expression shown in the ectoderm of the trochophore larva of snails is probably related to shell/mantle development. The observed sall expression in cephalic tissue in snails and some other bilaterians suggests a possible ancestral role of sall in neural development in bilaterians.

Abstract: Glycogen synthase kinase 3 (GSK3) regulates many cell fate decisions in animal development. In multicellular structures of the group 4 dictyostelid Dictyostelium discoideum, GSK3 promotes spore over stalk-like differentiation. We investigated whether, similar to other sporulation-inducing genes such as cAMP-dependent protein kinase (PKA), this role of GSK3 is derived from an ancestral role in encystation of unicellular amoebas. We deleted GSK3 in Polysphondylium pallidum, a group 2 dictyostelid which has retained encystation as an alternative survival strategy. Loss of GSK3 inhibited cytokinesis of cells in suspension, as also occurs in D. discoideum, but did not affect spore or stalk differentiation in P. pallidum. However, gsk3− amoebas entered into encystation under conditions that in wild type favour aggregation and fruiting body formation. The gsk3− cells were hypersensitive to osmolytes, which are known to promote encystation, and to cyst-inducing factors that are secreted during starvation. GSK3 was not itself regulated by these factors, but inhibited their effects. Our data show that GSK3 has a deeply conserved role in controlling cytokinesis, but not spore differentiation in Dictyostelia. Instead, in P. pallidum, one of many Dictyostelia that like their solitary ancestors can still encyst to survive starvation, GSK3 promotes multicellular development into fruiting bodies over unicellular encystment.

Abstract: In the last few years, accumulated information has indicated that the evolution of an extra-embryonic membrane in dipterans was accompanied by changes in the gene regulatory network controlled by the BMP/Dpp pathway, which is responsible for dorsal patterning in these insects. However, only comparative analysis of gene expression levels between distant species with two extra-embryonic membranes, like A. gambiae or C. albipunctata, and D. melanogaster, has been conducted. Analysis of gene expression in ancestral species, which evolved closer to the amnioserosa origin, could provide new insights into the evolution of dorsoventral patterning in dipterans. Here we describe the spatial expression of several key and downstream elements of the Dpp pathway and show the compared patterns of expression between Musca and Drosophila embryos, both dipterans with amnioserosa. Most of the analyzed gene showed a high degree of expression conservation, however, we found several differences in the gene expression pattern of M. domestica orthologs for sog and tolloid. Bioinformatics analysis of the promoter of both genes indicated that the variations could be related to the gain of several binding sites for the transcriptional factor Dorsal in the Md.tld promoter and Snail in the Md.sog enhancer. These altered expressions could explain the unclear formation of the pMad gradient in the M. domestica embryo, compared to the formation of the gradient in D. melanogaster. Gene expression changes during the dorsal–ventral patterning in insects contribute to the differentiation of extra-embryonic tissues as a consequence of changes in the gene regulatory network controlled by BMP/Dpp. In this work, in early M. domestica embryos, we identified the expression pattern of several genes members involved in the dorsoventral specification of the embryo. We believe that these data can contribute to understanding the evolution of the BMP/Dpp pathway, the regulation of BMP ligands, and the formation of a Dpp gradient in higher cyclorraphan flies.

Abstract: Morphogen signalling represents a key mechanism of developmental processes during animal development. Previously, several evolutionary conserved morphogen signalling pathways have been identified, and their players such as the morphogen receptors, morphogen modulating factors (MMFs) and the morphogens themselves have been studied. MMFs are factors that regulate morphogen distribution and activity. The interactions of MMFs with different morphogen signalling pathways such as Wnt signalling, Hedgehog (Hh) signalling and Decapentaplegic (Dpp) signalling are complex because some of the MMFs have been shown to interact with more than one signalling pathway, and depending on genetic context, to have different, biphasic or even opposing function. This complicates the interpretation of expression data and functional data of MMFs and may be one reason why data on MMFs in other arthropods than Drosophila are scarce or totally lacking. As a first step to a better understanding of the potential roles of MMFs in arthropod development, we investigate here the embryonic expression patterns of division abnormally delayed (dally), dally-like protein (dlp), shifted (shf) and secreted frizzled-related protein 125 (sFRP125) and sFRP34 in the beetle Tribolium castaneum, the spider Parasteatoda tepidariorum, the millipede Glomeris marginata and the onychophoran Euperipatoides kanangrensis. This pioneer study represents the first comprehensive comparative data set of these genes in panarthropods. Expression profiles reveal a high degree of diversity, suggesting that MMFs may represent highly evolvable nodes in otherwise conserved gene regulatory networks. Conserved aspects of MMF expression, however, appear to concern function in segmentation and limb development, two of the key topics of evolutionary developmental research.

Abstract: Wnt/β-catenin (or canonical) signalling pathway activity is necessary and used independently several times for specification of vegetal fate and endoderm, gut differentiation, maintenance of epithelium in adult intestine and the development of gut-derived organs in various vertebrate and non-vertebrate organisms. However, its conservation in later stages of digestive tract development still remains questionable due to the lack of detailed data, mainly from Spiralia. Here we characterize the Pdu-Tcf gene, a Tcf/LEF orthologue and a component of Wnt/β-catenin pathway from Platynereis dumerilii, a spiralian, marine annelid worm. Pdu-Tcf undergoes extensive alternative splicing in the C-terminal region of the gene generating as many as eight mRNA isoforms some of which differ in the presence or absence of a C-clamp domain which suggests a distinct DNA binding activity of individual protein variants. Pdu-Tcf is broadly expressed throughout development which is indicative of many functions. One of the most prominent domains that exhibits rather strong Pdu-Tcf expression is in the putative precursors of endodermal gut cells which are detected after 72 h post-fertilization (hpf). At day 5 post-fertilization (dpf), Pdu-Tcf is expressed in the hindgut and pharynx (foregut), whereas at 7 dpf stage, it is strongly transcribed in the now-cellularized midgut for the first time. In order to gain insight into the role of Wnt/β-catenin signalling, we disrupted its activity using pharmacological inhibitors between day 5 and 7 of development. The inhibition of Wnt/β-catenin signalling led to the loss of midgut marker genes Subtilisin-1, Subtilisin-2, α-Amylase and Otx along with a drop in β-catenin protein levels, Axin expression in the gut and nearly the complete loss of proliferative activity throughout the body of larva. At the same time, a hindgut marker gene Legumain was expanded to the midgut compartment under the same conditions. Our findings suggest that high Wnt/β-catenin signalling in the midgut might be necessary for proper differentiation of the endoderm to an epithelium capable of secreting digestive enzymes. Together, our data provide evidence for the role of Wnt/β-catenin signalling in gut differentiation in Platynereis.

Abstract: Both euarthropods and vertebrates have tripartite brains. Several orthologous genes are expressed in similar regionalized patterns during brain development in both vertebrates and euarthropods. These similarities have been used to support direct homology of the tripartite brains of vertebrates and euarthropods. If the tripartite brains of vertebrates and euarthropods are homologous, then one would expect other taxa to share this structure. More generally, examination of other taxa can help in tracing the evolutionary history of brain structures. Tardigrades are an interesting lineage on which to test this hypothesis because they are closely related to euarthropods, and whether they have a tripartite brain or unipartite brain has recently been a focus of debate. We tested this hypothesis by analyzing the expression patterns of six3, orthodenticle, pax6, unplugged, and pax2/5/8 during brain development in the tardigrade Hypsibius exemplaris—formerly misidentified as Hypsibius dujardini. These genes were expressed in a staggered anteroposterior order in H. exemplaris, similar to what has been reported for mice and flies. However, only six3, orthodenticle, and pax6 were expressed in the developing brain. Unplugged was expressed broadly throughout the trunk and posterior head, before the appearance of the nervous system. Pax2/5/8 was expressed in the developing central and peripheral nervous system in the trunk. Our results buttress the conclusion of our previous study of Hox genes—that the brain of tardigrades is only homologous to the protocerebrum of euarthropods. They support a model based on fossil evidence that the last common ancestor of tardigrades and euarthropods possessed a unipartite brain. Our results are inconsistent with the hypothesis that the tripartite brain of euarthropods is directly homologous to the tripartite brain of vertebrates.

Abstract: ETCHbox genes are eutherian-specific homeobox genes expressed during preimplantation development at a time when the first cell lineage decisions are being made. The mouse has an unusual repertoire of ETCHbox genes with several gene families lost in evolution and the remaining two, Crxos and Obox, greatly divergent in sequence and number. Each has undergone duplication to give a double homeodomain Crxos locus and a large cluster of over 60 Obox loci. The gene content differences between species raise important questions about how evolution can tolerate loss of genes implicated in key developmental events. We find that Crxos internal duplication occurred in the mouse lineage, while Obox duplication was stepwise, generating subgroups with distinct sequence and expression. Ectopic expression of three Obox genes and a Crxos transcript in primary mouse embryonic cells followed by transcriptome sequencing allowed investigation into their functional roles. We find distinct transcriptomic influences for different Obox subgroups and Crxos, including modulation of genes related to zygotic genome activation and preparation for blastocyst formation. Comparison with similar experiments performed using human homeobox genes reveals striking overlap between genes downstream of mouse Crxos and genes downstream of human ARGFX. Mouse Crxos and human ARGFX homeobox genes are paralogous rather than orthologous, yet they have evolved to regulate a common set of genes. This suggests there was compensation of function alongside gene loss through co-option of a different locus. Functional compensation by non-orthologous genes with dissimilar sequences is unusual but may indicate underlying distributed robustness. Compensation may be driven by the strong evolutionary pressure for successful early embryo development.

Abstract: The anterior neuroectoderm (ANE) in many deuterostome embryos (echinoderms, hemichordates, urochordates, cephalochordates, and vertebrates) is progressively restricted along the anterior–posterior axis to a domain around the anterior pole. In the sea urchin embryo, three integrated Wnt signaling branches (Wnt/β-catenin, Wnt/JNK, and Wnt/PKC) govern this progressive restriction process, which begins around the 32- to 60-cell stage and terminates by the early gastrula stage. We previously have established that several secreted Wnt modulators of the Dickkopf and secreted Frizzled-related protein families (Dkk1, Dkk3, and sFRP-1/5) are expressed within the ANE and play important roles in modulating the Wnt signaling network during this process. In this study, we use morpholino and dominant-negative interference approaches to characterize the function of a novel Frizzled-related protein, secreted Frizzled-related protein 1 (sFRP-1), during ANE restriction. sFRP-1 appears to be related to a secreted Wnt modulator, sFRP3/4, that is essential to block Wnt signaling and establish the ANE in vertebrates. Here, we show that the sea urchin sFRP3/4 orthologue is not expressed during ANE restriction in the sea urchin embryo. Instead, our results indicate that ubiquitously expressed maternal sFRP-1 and Fzl1/2/7 signaling act together as early as the 32- to 60-cell stage to antagonize the ANE restriction mechanism mediated by Wnt/β-catenin and Wnt/JNK signaling. Then, starting from the blastula stage, Fzl5/8 signaling activates zygotic sFRP-1 within the ANE territory, where it works with the secreted Wnt antagonist Dkk1 (also activated by Fzl5/8 signaling) to antagonize Wnt1/Wnt8–Fzl5/8–JNK signaling in a negative feedback mechanism that defines the outer ANE territory boundary. Together, these data indicate that maternal and zygotic sFRP-1 protects the ANE territory by antagonizing the Wnt1/Wnt8–Fzl5/8–JNK signaling pathway throughout ANE restriction, providing precise spatiotemporal control of the mechanism responsible for the establishment of the ANE territory around the anterior pole of the sea urchin embryo.

Abstract: Tapeworms are agents of neglected tropical diseases responsible for significant health problems and economic loss. They also exhibit adaptations to a parasitic lifestyle that confound comparisons of their development with other animals. Identifying the genetic factors regulating their complex ontogeny is essential to understanding unique aspects of their biology and for advancing novel therapeutics. Here we use RNA sequencing to identify up-regulated signalling components, transcription factors and post-transcriptional/translational regulators (genes of interest, GOI) in the transcriptomes of Larvae and different regions of segmented worms in the tapeworm Hymenolepis microstoma and combine this with spatial gene expression analyses of a selection of genes. RNA-seq reads collectively mapped to 90% of the > 12,000 gene models in the H. microstoma v.2 genome assembly, demonstrating that the transcriptome profiles captured a high percentage of predicted genes. Contrasts made between the transcriptomes of Larvae and whole, adult worms, and between the Scolex-Neck, mature strobila and gravid strobila, resulted in 4.5–30% of the genes determined to be differentially expressed. Among these, we identified 190 unique GOI up-regulated in one or more contrasts, including a large range of zinc finger, homeobox and other transcription factors, components of Wnt, Notch, Hedgehog and TGF-β/BMP signalling, and post-transcriptional regulators (e.g. Boule, Pumilio). Heatmap clusterings based on overall expression and on select groups of genes representing ‘signals’ and ‘switches’ showed that expression in the Scolex-Neck region is more similar to that of Larvae than to the mature or gravid regions of the adult worm, which was further reflected in large overlap of up-regulated GOI. Spatial expression analyses in Larvae and adult worms corroborated inferences made from quantitative RNA-seq data and in most cases indicated consistency with canonical roles of the genes in other animals, including free-living flatworms. Recapitulation of developmental factors up-regulated during larval metamorphosis suggests that strobilar growth involves many of the same underlying gene regulatory networks despite the significant disparity in developmental outcomes. The majority of genes identified were investigated in tapeworms for the first time, setting the stage for advancing our understanding of developmental genetics in an important group of flatworm parasites.