Clustered Fox genes in lophotrochozoans and the evolution of the bilaterian Fox gene cluster

Evolution of the expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

The Fox Gene Repertoire in the Annelid Owenia fusiformis Reveals Multiple Expansions of the foxQ2 Class in Spiralia

Fox genes are a large and conserved family of transcription factors involved in many key biological processes, including embryogenesis and body patterning. Although the role of Fox genes has been studied in an array of model systems, comprehensive comparative studies in Spiralia-a large clade of invertebrate animals including molluscs and annelids-are scarce but much needed to better understand the evolutionary history of this gene family. Here, we reconstruct and functionally characterize the Fox gene complement in the annelid Owenia fusiformis, a slow evolving species and member of the sister group to all remaining annelids. The genome of O. fusiformis contains at least a single ortholog for 20 of the 22 Fox gene classes that are ancestral to Bilateria, including an ortholog of the recently discovered foxT class. Temporal and spatial expression dynamics reveal a conserved role of Fox genes in gut formation, mesoderm patterning, and apical organ and cilia formation in Annelida and Spiralia. Moreover, we uncover an ancestral expansion of foxQ2 genes in Spiralia, represented by 11 paralogs in O. fusiformis. Notably, although all foxQ2 copies have apical expression in O. fusiformis, they show variable spatial domains and staggered temporal activation, which suggest cooperation and sub-functionalization among foxQ2 genes for the development of apical fates in this annelid. Altogether, our study informs the evolution and developmental roles of Fox genes in Annelida and Spiralia generally, providing the basis to explore how regulatory changes in Fox gene expression might have contributed to developmental and morphological diversification in Spiralia.

Gene activation of metazoan Fox transcription factors at the onset of metamorphosis in the marine demosponge Amphimedon queenslandica

Transcription factors encoded by the Forkhead (Fox) gene family have diverse, sometimes conserved, regulatory roles in eumetazoan development, immunity, and physiology. Although this gene family includes members that predate the origin of the animal kingdom, the majority of metazoan Fox genes evolved after the divergence of animals and choanoflagellates. Here, we characterize the composition, structure, and expression of Fox genes in the marine demosponge Amphimedon queenslandica to better understand the origin and evolution of this family. The Fox gene repertoire in A. queenslandica appears to be similar to the ancestral metazoan Fox gene family. All 17 A. queenslandica Fox genes are differentially expressed during development and in adult cell types. Remarkably, eight of these, all of which appear to be metazoan-specific, are induced within just 1 h of larval settlement and commencement of metamorphosis. Gene co-expression analyses suggest that these eight Fox genes regulate developmental and physiological processes similar to their roles in other animals. These findings are consistent with Fox genes playing deeply ancestral roles in animal development and physiology, including in response to changes in the external environment.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

Phylogenetic analysis of forkhead transcription factors in the Panarthropoda

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.

Placozoan fiber cells: mediators of innate immunity and participants in wound healing

Placozoa is a phylum of non-bilaterian marine animals. These small, flat organisms adhere to the substrate via their densely ciliated ventral epithelium, which mediates mucociliary locomotion and nutrient uptake. They have only six morphological cell types, including one, fiber cells, for which functional data is lacking. Fiber cells are non-epithelial cells with multiple processes. We used electron and light microscopic approaches to unravel the roles of fiber cells in Trichoplax adhaerens, a representative member of the phylum. Three-dimensional reconstructions of serial sections of Trichoplax showed that each fiber cell is in contact with several other cells. Examination of fiber cells in thin sections and observations of live dissociated fiber cells demonstrated that they phagocytose cell debris and bacteria. In situ hybridization confirmed that fiber cells express genes involved in phagocytic activity. Fiber cells also are involved in wound healing as evidenced from microsurgery experiments. Based on these observations we conclude that fiber cells are multi-purpose macrophage-like cells. Macrophage-like cells have been described in Porifera, Ctenophora, and Cnidaria and are widespread among Bilateria, but our study is the first to show that Placozoa possesses this cell type. The phylogenetic distribution of macrophage-like cells suggests that they appeared early in metazoan evolution.

FGF signaling acts on different levels of mesoderm development within Spiralia

FGF signaling is involved in mesoderm induction in members of deuterostomes (e.g. tunicates, hemichordates), but not in flies and nematodes, in which it has a role in mesoderm patterning and migration. However, we need comparable studies in other protostome taxa in order to decipher whether this mesoderm-inducing function of FGF extends beyond the lineage of deuterostomes. Here, we investigated the role of FGF signaling in mesoderm development in three species of lophophorates, a clade within the protostome group Spiralia. Our gene expression analyses show that the mesodermal molecular patterning is conserved between brachiopods and phoronids, but the spatial and temporal recruitment of transcription factors differs significantly. Moreover, the use of the inhibitor SU5402 demonstrates that FGF signaling is involved in different steps of mesoderm development, as well as in morphogenetic movements of gastrulation and axial elongation. Our findings suggest that the mesoderm-inducing role of FGF extends beyond the group of deuterostomes.

Analysis of Fox genes in Schmidtea mediterranea reveals new families and a conserved role of Smed-foxO in controlling cell death

The forkhead box (Fox) genes encode transcription factors that control several key aspects of development. Present in the ancestor of all eukaryotes, Fox genes underwent several duplications followed by loss and diversification events that gave rise to the current 25 families. However, few Fox members have been identified from the Lophotrochozoa clade, and specifically from planarians, which are a unique model for understanding development, due to the striking plasticity of the adult. The aim of this study was to identify and perform evolutionary and functional studies of the Fox genes of lophotrochozoan species and, specifically, of the planarian Schmidtea mediterranea. Generating a pipeline for identifying Forkhead domains and using phylogenetics allowed us the phylogenetic reconstruction of Fox genes. We corrected the annotation for misannotated genes and uncovered a new family, the QD, present in all metazoans. According to the new phylogeny, the 27 Fox genes found in Schmidtea mediterranea were classified into 12 families. In Platyhelminthes, family losses were accompanied by extensive gene diversification and the appearance of specific families, the A(P) and N(P). Among the newly identified planarian Fox genes, we found a single copy of foxO, which shows an evolutionary conserved role in controlling cell death.

Identification and expression profiles of Fox transcription factors in the Yesso scallop (Patinopecten yessoensis)

The forkhead box (Fox) gene family is a family of transcription factors that play important roles in a variety of biological processes in vertebrates, including early development and cell proliferation and differentiation. However, at present, studies on the mollusk Fox family are relatively lacking. In the present study, the Fox gene family of the Yesso scallop (Patinopecten yessoensis) was systematically identified. In addition, the expression profiles of the Fox gene family in early development and adult tissues were analyzed. The results showed that there were 26 Fox genes in P. yessoensis. Of the 26 genes, 24 belonged to 20 subfamilies. The Fox genes belonging to the I, Q1, R and S subfamilies were absent in P. yessoensis. The other 2 genes formed 2 independent clades with the Fox genes of other mollusks and protostomes. They might be new members of the Fox family and were named FoxY and FoxZ. P. yessoensis contained a FoxC-FoxL1 gene cluster similar in structure to that of Branchiostoma floridae, suggesting that the cluster might already exist in the ancestors of bilaterally symmetrical animals. The gene expression analysis of Fox showed that most of the genes were continuously expressed in multiple stages of early development, suggesting that Fox genes might be widely involved in the regulation of embryo and larval development of P. yessoensis. Nine Fox genes were specifically expressed in certain tissues, such as the nerve ganglia, foot, ovary, testis, and gills. For the 9 genes that were differentially expressed between the testis and ovary, their expression levels were analyzed during the 4 developmental stages of gonads. The results showed that FoxL2, FoxE and FoxY were highly expressed in the ovary during all developmental stages, while FoxZ was highly expressed in the testis during all developmental stages. The results suggested that these genes might play an important role in sex maintenance or gametogenesis. The present study could provide a reference for evolutionary and functional studies of the Fox family in metazoans.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

Evolution of the bilaterian mouth and anus

It is widely held that the bilaterian tubular gut with mouth and anus evolved from a simple gut with one major gastric opening. However, there is no consensus on how this happened. Did the single gastric opening evolve into a mouth, with the anus forming elsewhere in the body (protostomy), or did it evolve into an anus, with the mouth forming elsewhere (deuterostomy), or did it evolve into both mouth and anus (amphistomy)? These questions are addressed by the comparison of developmental fates of the blastopore, the opening of the embryonic gut, in diverse animals that live today. Here we review comparative data on the identity and fate of blastoporal tissue, investigate how the formation of the through-gut relates to the major body axes, and discuss to what extent evolutionary scenarios are consistent with these data. Available evidence indicates that stem bilaterians had a slit-like gastric opening that was partially closed in subsequent evolution, leaving open the anus and most likely also the mouth, which would favour amphistomy. We discuss remaining difficulties, and outline directions for future research.

Proteome and Transcriptome Analysis of Ovary, Intersex Gonads, and Testis Reveals Potential Key Sex Reversal/Differentiation Genes and Mechanism in Scallop Chlamys nobilis

Bivalve mollusks exhibit hermaphroditism and sex reversal/differentiation. Studies generally focus on transcriptional profiling and specific genes related to sex determination and differentiation. Few studies on sex reversal/differentiation have been reported. A combination analysis of gonad proteomics and transcriptomics was conducted on Chlamys nobilis to provide a systematic understanding of sex reversal/differentiation in bivalves. We obtained 4258 unique peptides and 93,731 unigenes with good correlation between messenger RNA and protein levels. Candidate genes in sex reversal/differentiation were found: 15 genes differentially expressed between sexes were identified and 12 had obvious sexual functions. Three novel genes (foxl2, β-catenin, and sry) were expressed highly in intersex individuals and were likely involved in the control of gonadal sex in C. nobilis. High expression of foxl2 or β-catenin may inhibit sry and activate 5-HT receptor and vitellogenin to maintain female development. High expression of sry may inhibit foxl2 and β-catenin and activate dmrt2, fem-1, sfp2, sa6, Amy-1, APCP4, and PLK to maintain male function. High expression of sry, foxl2, and β-catenin in C. nobilis may be involved in promoting and maintaining sex reversal/differentiation. The downstream regulator may not be dimorphic expressed genes, but genes expressed in intersex individuals, males and females. Different expression patterns of sex-related genes and gonadal histological characteristics suggested that C. nobilis may change its sex from male to female. These findings suggest highly conserved sex reversal/differentiation with diverged regulatory pathways during C. nobilis evolution. This study provides valuable genetic resources for understanding sex reversal/differentiation (intersex) mechanisms and pathways underlying bivalve reproductive regulation.

Phylogenetics of Lophotrochozoan bHLH Genes and the Evolution of Lineage-Specific Gene Duplicates

The gain and loss of genes encoding transcription factors is of importance to understanding the evolution of gene regulatory complexity. The basic helix–loop–helix (bHLH) genes encode a large superfamily of transcription factors. We systematically classify the bHLH genes from five mollusc, two annelid and one brachiopod genomes, tracing the pattern of bHLH gene evolution across these poorly studied Phyla. In total, 56–88 bHLH genes were identified in each genome, with most identifiable as members of previously described bilaterian families, or of new families we define. Of such families only one, Mesp, appears lost by all these species. Additional duplications have also played a role in the evolution of the bHLH gene repertoire, with many new lophotrochozoan-, mollusc-, bivalve-, or gastropod-specific genes defined. Using a combination of transcriptome mining, RT-PCR, and in situ hybridization we compared the expression of several of these novel genes in tissues and embryos of the molluscs Crassostrea gigas and Patella vulgata, finding both conserved expression and evidence for neofunctionalization. We also map the positions of the genes across these genomes, identifying numerous gene linkages. Some reflect recent paralog divergence by tandem duplication, others are remnants of ancient tandem duplications dating to the lophotrochozoan or bilaterian common ancestors. These data are built into a model of the evolution of bHLH genes in molluscs, showing formidable evolutionary stasis at the family level but considerable within-family diversification by tandem gene duplication.

Cleavage modification did not alter blastomere fates during bryozoan evolution

BACKGROUND

Stereotypic cleavage patterns play a crucial role in cell fate determination by precisely positioning early embryonic blastomeres. Although misplaced cell divisions can alter blastomere fates and cause embryonic defects, cleavage patterns have been modified several times during animal evolution. However, it remains unclear how evolutionary changes in cleavage impact the specification of blastomere fates. Here, we analyze the transition from spiral cleavage - a stereotypic pattern remarkably conserved in many protostomes - to a biradial cleavage pattern, which occurred during the evolution of bryozoans.

RESULTS

Using 3D-live imaging time-lapse microscopy (4D-microscopy), we characterize the cell lineage, MAPK signaling, and the expression of 16 developmental genes in the bryozoan Membranipora membranacea. We found that the molecular identity and the fates of early bryozoan blastomeres are similar to the putative homologous blastomeres in spiral-cleaving embryos.

CONCLUSIONS

Our work suggests that bryozoans have retained traits of spiral development, such as the early embryonic fate map, despite the evolution of a novel cleavage geometry. These findings provide additional support that stereotypic cleavage patterns can be modified during evolution without major changes to the molecular identity and fate of embryonic blastomeres.

The evolutionary origin of bilaterian smooth and striated myocytes

The dichotomy between smooth and striated myocytes is fundamental for bilaterian musculature, but its evolutionary origin is unsolved. In particular, interrelationships of visceral smooth muscles remain unclear. Absent in fly and nematode, they have not yet been characterized molecularly outside vertebrates. Here, we characterize expression profile, ultrastructure, contractility and innervation of the musculature in the marine annelid Platynereis dumerilii and identify smooth muscles around the midgut, hindgut and heart that resemble their vertebrate counterparts in molecular fingerprint, contraction speed and nervous control. Our data suggest that both visceral smooth and somatic striated myocytes were present in the protostome-deuterostome ancestor and that smooth myocytes later co-opted the striated contractile module repeatedly – for example, in vertebrate heart evolution. During these smooth-to-striated myocyte conversions, the core regulatory complex of transcription factors conveying myocyte identity remained unchanged, reflecting a general principle in cell type evolution.

DOI:http://dx.doi.org/10.7554/eLife.19607.001

Logics and properties of a genetic regulatory program that drives embryonic muscle development in an echinoderm

Evolutionary origin of muscle is a central question when discussing mesoderm evolution. Developmental mechanisms underlying somatic muscle development have mostly been studied in vertebrates and fly where multiple signals and hierarchic genetic regulatory cascades selectively specify myoblasts from a pool of naive mesodermal progenitors. However, due to the increased organismic complexity and distant phylogenetic position of the two systems, a general mechanistic understanding of myogenesis is still lacking. In this study, we propose a gene regulatory network (GRN) model that promotes myogenesis in the sea urchin embryo, an early branching deuterostome. A fibroblast growth factor signaling and four Forkhead transcription factors consist the central part of our model and appear to orchestrate the myogenic process. The topological properties of the network reveal dense gene interwiring and a multilevel transcriptional regulation of conserved and novel myogenic genes. Finally, the comparison of the myogenic network architecture among different animal groups highlights the evolutionary plasticity of developmental GRNs.

DOI:http://dx.doi.org/10.7554/eLife.07343.001

Muscles, bones, and blood vessels all develop from a tissue called the mesoderm, which forms early on in the development of an embryo. Networks of genes control which parts of the mesoderm transform into different cell types. The gene networks that control the development of muscle cells from the mesoderm have so far been investigated in flies and several species of animals with backbones. However, these species are complex, which makes it difficult to work out the general principles that control muscle cell development.

Sea urchins are often studied in developmental biology as they have many of the same genes as more complex animals, but are much simpler and easier to study in the laboratory. Andrikou et al. therefore investigated the ‘gene regulatory network’ that controls muscle development in sea urchins. This revealed that proteins called Forkhead transcription factors and a process called FGF signaling are crucial for controlling muscle development in sea urchins. These are also important factors for developing muscles in other animals.

Andrikou et al. then produced models that show the interactions between the genes that control muscle formation at three different stages of embryonic development. These models reveal several important features of the muscle development gene regulatory network. For example, the network is robust: if one gene fails, the network is connected in a way that allows it to still make muscle. This also allows the network to adapt and evolve without losing the ability to perform any of its existing roles.

Comparing the gene regulatory network that controls muscle development in sea urchins with the networks found in other animals showed that many of the same genes are used across different species, but are connected into different network structures. Investigating the similarities and differences of the regulatory networks in different species could help us to understand how muscles have evolved and could ultimately lead to a better understanding of the causes of developmental diseases.

DOI:http://dx.doi.org/10.7554/eLife.07343.002

Comparative Transcriptome Analysis Reveals Sex-Biased Gene Expression in Juvenile Chinese Mitten Crab Eriocheir sinensis

Sex-biased genes are considered to account for most of phenotypic differences between males and females. In order to explore the sex-biased gene expression in crab, we performed the whole-body transcriptome analysis in male and female juveniles of the Chinese mitten crab Eriocheir sinensis using next-generation sequencing technology. Of the 23,349 annotated unigenes, 148 were identified as sex-related genes. A total of 29 candidate genes involved in primary sex determination pathways were detected, indicating the sex determination cascade of the mitten crab might be more complex than previously supposed. Differential expression analysis showed 448 differentially expressed genes (DEGs) between the two transcriptomes. Most of DEGs were involved in processes such as metabolism and immunity, and not associated with obvious sexual function. The pathway predominantly enriched for DEGs were related to lysosome, which might reflect the differences in metabolism between males and females. Of the immune DGEs, 18 up-regulated genes in females were humoral immune factors, and eight up-regulated genes in males were pattern recognition receptors, suggesting sex differences of immune defense might exist in the mitten crab. In addition, two reproduction-related genes, vitellogenin and insulin-like androgenic gland factor, were identified to express in both sexes but with significantly higher level in males. Our research provides the first whole-body RNA sequencing of sex-specific transcriptomes for juvenile E. sinensis and will facilitate further studies on molecular mechanisms of crab sexual dimorphism.

Genome of the Rusty Millipede, Trigoniulus corallinus, Illuminates Diplopod, Myriapod, and Arthropod Evolution

The increasing availability of genomic information from the Arthropoda continues to revolutionize our understanding of the biology of this most diverse animal phylum. However, our sampling of arthropod diversity remains uneven, and key clade such as the Myriapoda are severely underrepresented. Here we present the genome of the cosmopolitanly distributed Rusty Millipede Trigoniulus corallinus, which represents the first diplopod genome to be published, and the second example from the Myriapoda as a whole. This genomic resource contains the majority of core eukaryotic genes (94.3%), and key transcription factor classes that were thought to be lost in the Ecdysozoa. Mitochondrial genome and gene family (transcription factor, Dscam, circadian clock-driving protein, odorant receptor cassette, bioactive compound, and cuticular protein) analyses were also carried out to shed light on their states in the Diplopoda and Myriapoda. The ready availability of T. corallinus recommends it as a new model for evolutionary developmental biology, and the data set described here will be of widespread utility in investigating myriapod and arthropod genomics and evolution.

Mesodermal gene expression during the embryonic and larval development of the articulate brachiopod Terebratalia transversa

Background

Brachiopods undergo radial cleavage, which is distinct from the stereotyped development of closely related spiralian taxa. The mesoderm has been inferred to derive from the archenteron walls following gastrulation, and the primary mesoderm derivative in the larva is a complex musculature. To investigate the specification and differentiation of the mesoderm in the articulate brachiopod Terebratalia transversa, we have identified orthologs of genes involved in mesoderm development in other taxa and investigated their spatial and temporal expression during the embryonic and larval development of T. transversa.

Results

Orthologs of 17 developmental regulatory genes with roles in the development of the mesoderm in other bilaterian animals were found to be expressed in the developing mesoderm of T. transversa. Five genes, Tt.twist, Tt.GATA456, Tt.dachshund, Tt.mPrx, and Tt.NK1, were found to have expression throughout the archenteron wall at the radial gastrula stage, shortly after the initiation of gastrulation. Three additional genes, Tt.Pax1/9, Tt.MyoD, and Tt.Six1/2, showed expression at this stage in only a portion of the archenteron wall. Tt.eya, Tt.FoxC, Tt.FoxF, Tt.Mox, Tt.paraxis, Tt.Limpet, and Tt.Mef2 all showed initial mesodermal expression during later gastrula or early larval stages. At the late larval stage, Tt.dachshund, Tt.Limpet, and Tt.Mef2 showed expression in nearly all mesoderm cells, while all other genes were localized to specific regions of the mesoderm. Tt.FoxD and Tt.noggin both showed expression in the ventral mesoderm at the larval stages, with gastrula expression patterns in the archenteron roof and blastopore lip, respectively.

Conclusions

Expression analyses support conserved roles for developmental regulators in the specification and differentiation of the mesoderm during the development of T. transversa. Expression of multiple mesodermal factors in the archenteron wall during gastrulation supports previous morphological observations that this region gives rise to larval mesoderm. Localized expression domains during gastrulation and larval development evidence early regionalization of the mesoderm and provide a basis for hypotheses regarding the molecular regulation underlying the complex system of musculature observed in the larva.

Electronic supplementary material

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

Did homeobox gene duplications contribute to the Cambrian explosion?

The Cambrian explosion describes an apparently rapid increase in the diversity of bilaterian animals around 540-515 million years ago. Bilaterian animals explore the world in three-dimensions deploying forward-facing sense organs, a brain, and an anterior mouth; they possess muscle blocks enabling efficient crawling and burrowing in sediments, and they typically have an efficient 'through-gut' with separate mouth and anus to process bulk food and eject waste, even when burrowing in sediment. A variety of ecological, environmental, genetic, and developmental factors have been proposed as possible triggers and correlates of the Cambrian explosion, and it is likely that a combination of factors were involved. Here, I focus on a set of developmental genetic changes and propose these are part of the mix of permissive factors. I describe how ANTP-class homeobox genes, which encode transcription factors involved in body patterning, increased in number in the bilaterian stem lineage and earlier. These gene duplications generated a large array of ANTP class genes, including three distinct gene clusters called NK, Hox, and ParaHox. Comparative data supports the idea that NK genes were deployed primarily to pattern the bilaterian mesoderm, Hox genes coded position along the central nervous system, and ParaHox genes most likely originally specified the mouth, midgut, and anus of the newly evolved through-gut. It is proposed that diversification of ANTP class genes played a role in the Cambrian explosion by contributing to the patterning systems used to build animal bodies capable of high-energy directed locomotion, including active burrowing.

Molecular conservation of metazoan gut formation: evidence from expression of endomesoderm genes in Capitella teleta (Annelida)

BACKGROUND

Metazoan digestive systems develop from derivatives of ectoderm, endoderm and mesoderm, and vary in the relative contribution of each germ layer across taxa and between gut regions. In a small number of well-studied model systems, gene regulatory networks specify endoderm and mesoderm of the gut within a bipotential germ layer precursor, the endomesoderm. Few studies have examined expression of endomesoderm genes outside of those models, and thus, it is unknown whether molecular specification of gut formation is broadly conserved. In this study, we utilize a sequenced genome and comprehensive fate map to correlate the expression patterns of six transcription factors with embryonic germ layers and gut subregions during early development in Capitella teleta.

RESULTS

The genome of C. teleta contains the five core genes of the sea urchin endomesoderm specification network. Here, we extend a previous study and characterize expression patterns of three network orthologs and three additional genes by in situ hybridization during cleavage and gastrulation stages and during formation of distinct gut subregions. In cleavage stage embryos, Ct-otx, Ct-blimp1, Ct-bra and Ct-nkx2.1a are expressed in all four macromeres, the endoderm precursors. Ct-otx, Ct-blimp1, and Ct-nkx2.1a are also expressed in presumptive endoderm of gastrulae and later during midgut development. Additional gut-specific expression patterns include Ct-otx, Ct-bra, Ct-foxAB and Ct-gsc in oral ectoderm; Ct-otx, Ct-blimp1, Ct-bra and Ct-nkx2.1a in the foregut; and both Ct-bra and Ct-nkx2.1a in the hindgut.

CONCLUSIONS

Identification of core sea urchin endomesoderm genes in C. teleta indicates they are present in all three bilaterian superclades. Expression of Ct-otx, Ct-blimp1 and Ct-bra, combined with previously published Ct-foxA and Ct-gataB1 patterns, provide the most comprehensive comparison of these five orthologs from a single species within Spiralia. Each ortholog is likely involved in endoderm specification and midgut development, and several may be essential for establishment of the oral ectoderm, foregut and hindgut, including specification of ectodermal and mesodermal contributions. When the five core genes are compared across the Metazoa, their conserved expression patterns suggest that 'gut gene' networks evolved to specify distinct digestive system subregions, regardless of species-specific differences in gut architecture or germ layer contributions within each subregion.

Deregulated FOX genes in Hodgkin lymphoma

FOX genes encode transcription factors which regulate basic developmental processes during embryogenesis and in the adult. Several FOX genes show deregulated expression in particular malignancies, representing oncogenes or tumor suppressors. Here, we screened six Hodgkin lymphoma (HL) cell lines for FOX gene activity by comparative microarray profiling, revealing overexpression of FOXC1 and FOXD1, and reduced transcription of FOXN3, FOXO1, and FOXP1. In silico expression analyses of these FOX gene candidates in HL patient samples supported the cell line data. Chromosomal analyses demonstrated an amplification of the FOXC1 locus at 6p25 and a gain of the FOXR2 locus at Xp11, indicting genomic aberrations for their upregulation. Comparative expression profiling and ensuing stimulation experiments revealed implementation of the TGFβ- and WNT-signaling pathways in deregulation of FOXD1 and FOXN3. Functional analysis of FOXP1 implicated miR9 and miR34a as upstream regulators and PAX5, TCF3, and RAG2 as downstream targets. A similar exercise for FOXC1 revealed repression of MSX1 and activation of IPO7, both mediating inhibition of the B-cell specific homeobox gene ZHX2. Taken together, our data show that aberrantly expressed FOX genes and their downstream targets are involved in the pathogenesis of HL via deregulation of B-cell differentiation and may represent useful diagnostic markers and/or therapeutic targets.

Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes

BACKGROUND

Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis.

RESULTS

The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera.

CONCLUSIONS

The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

The Fox/Forkhead transcription factor family of the hemichordate Saccoglossus kowalevskii

Background

The Fox gene family is a large family of transcription factors that arose early in organismal evolution dating back to at least the common ancestor of metazoans and fungi. They are key components of many gene regulatory networks essential for embryonic development. Although much is known about the role of Fox genes during vertebrate development, comprehensive comparative studies outside vertebrates are sparse. We have characterized the Fox transcription factor gene family from the genome of the enteropneust hemichordate Saccoglossus kowalevskii, including phylogenetic analysis, genomic organization, and expression analysis during early development. Hemichordates are a sister group to echinoderms, closely related to chordates and are a key group for tracing the evolution of gene regulatory mechanisms likely to have been important in the diversification of the deuterostome phyla.

Results

Of the 22 Fox gene families that were likely present in the last common ancestor of all deuterostomes, S. kowalevskii has a single ortholog of each group except FoxH, which we were unable to detect, and FoxQ2, which has three paralogs. A phylogenetic analysis of the FoxQ2 family identified an ancestral duplication in the FoxQ2 lineage at the base of the bilaterians. The expression analyses of all 23 Fox genes of S. kowalevskii provide insights into the evolution of components of the regulatory networks for the development of pharyngeal gill slits (foxC, foxL1, and foxI), mesoderm patterning (foxD, foxF, foxG), hindgut development (foxD, foxI), cilia formation (foxJ1), and patterning of the embryonic apical territory (foxQ2).

Conclusions

Comparisons of our results with data from echinoderms, chordates, and other bilaterians help to develop hypotheses about the developmental roles of Fox genes that likely characterized ancestral deuterostomes and bilaterians, and more recent clade-specific innovations.

Discovery of sex-related genes through high-throughput transcriptome sequencing from the salmon louse Caligus rogercresseyi

Understanding the molecular underpinnings involved in the reproduction of the salmon louse is critical for designing novel strategies of pest management for this ectoparasite. However, genomic information on sex-related genes is still limited. In the present work, sex-specific gene transcription was revealed in the salmon louse Caligus rogercresseyi using high-throughput Illumina sequencing. A total of 30,191,914 and 32,292,250 high quality reads were generated for females and males, and these were de novo assembled into 32,173 and 38,177 contigs, respectively. Gene ontology analysis showed a pattern of higher expression in the female as compared to the male transcriptome. Based on our sequence analysis and known sex-related proteins, several genes putatively involved in sex differentiation, including Dmrt3, FOXL2, VASA, and FEM1, and other potentially significant candidate genes in C. rogercresseyi, were identified for the first time. In addition, the occurrence of SNPs in several differentially expressed contigs annotating for sex-related genes was found. This transcriptome dataset provides a useful resource for future functional analyses, opening new opportunities for sea lice pest control.

Reproduction-related genes in the pearl oyster genome

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes

The gastropod Patella vulgata is abundant on rocky shores in Northern Europe and a significant grazer of intertidal algae. Here we report the application of Illumina sequencing to develop a transcriptome from the adult mantle tissue of P. vulgata. We obtained 47,237,104 paired-end reads of 51 bp, trialled de novo assembly methods and settled on the additive multiple K method followed by redundancy removal as resulting in the most comprehensive assembly. This yielded 29,489 contigs of at least 500 bp in length. We then used three methods to search for candidate genes relevant to biomineralisation: searches via BLAST and Hidden Markov Models for homologues of biomineralising genes from other molluscs, searches for predicted proteins containing tandem repeats and searches for secreted proteins that lacked a transmembrane domain. From the results of these searches we selected 15 contigs for verification by RT-PCR, of which 14 were successfully amplified and cloned. These included homologues of Pif-177/BSMP, Perlustrin, SPARC, AP24, Follistatin-like and Carbonic anhydrase, as well as three containing extensive G-X-Y repeats as found in nacrein. We selected two for further verification by in situ hybridisation, demonstrating expression in the larval shell field. We conclude that de novo assembly of Illumina data offers a cheap and rapid route to a predicted transcriptome that can be used as a resource for further biological study.

Mesodermal gene expression in the acoel Isodiametra pulchra indicates a low number of mesodermal cell types and the endomesodermal origin of the gonads

Acoelomorphs are bilaterally symmetric small marine worms that lack a coelom and possess a digestive system with a single opening. Two alternative phylogenetic positions of this group within the animal tree are currently debated. In one view, Acoelomorpha is the sister group to all remaining Bilateria and as such, is a morphologically simple stepping stone in bilaterian evolution. In the other, the group is a lineage within the Deuterostomia, and therefore, has derived a simple morphology from a more complex ancestor. Acoels and the closely related Nemertodermatida and Xenoturbellida, which together form the Acoelomorpha, possess a very limited number of cell types. To further investigate the diversity and origin of mesodermal cell types we describe the expression pattern of 12 orthologs of bilaterian mesodermal markers including Six1/2, Twist, FoxC, GATA4/5/6, in the acoel Isodiametra pulchra. All the genes are expressed in stem cells (neoblasts), gonads, and at least subsets of the acoel musculature. Most are expressed in endomesodermal compartments of I. pulchra developing embryos similar to what has been described in cnidarians. Our molecular evidence indicates a very limited number of mesodermal cell types and suggests an endomesodermal origin of the gonads and the stem cell system. We discuss our results in light of the two prevailing phylogenetic positions of Acoelomorpha.

Insights into bilaterian evolution from three spiralian genomes

Current genomic perspectives on animal diversity neglect two prominent phyla, the molluscs and annelids, that together account for nearly one-third of known marine species and are important both ecologically and as experimental systems in classical embryology^1–3. Here we describe the draft genomes of the owl limpet (Lottia gigantea), a marine polychaete (Capitella teleta) and a freshwater leech (Helobdella robusta), and compare them with other animal genomes to investigate the origin and diversification of bilaterians from a genomic perspective. We find that the genome organization, gene structure and functional content of these species are more similar to those of some invertebrate deuterostome genomes (for example, amphioxus and sea urchin) than those of other protostomes that have been sequenced to date (flies, nematodes and flatworms). The conservation of these genomic features enables us to expand the inventory of genes present in the last common bilaterian ancestor, establish the tripartite diversification of bilaterians using multiple genomic characteristics and identify ancient conserved long- and short-range genetic linkages across metazoans. Superimposed on this broadly conserved pan-bilaterian background we find examples of lineage-specific genome evolution, including varying rates of rearrangement, intron gain and loss, expansions and contractions of gene families, and the evolution of clade-specific genes that produce the unique content of each genome.

Evolutionary crossroads in developmental biology: annelids

Annelids (the segmented worms) have a long history in studies of animal developmental biology, particularly with regards to their cleavage patterns during early development and their neurobiology. With the relatively recent reorganisation of the phylogeny of the animal kingdom, and the distinction of the super-phyla Ecdysozoa and Lophotrochozoa, an extra stimulus for studying this phylum has arisen. As one of the major phyla within Lophotrochozoa, Annelida are playing an important role in deducing the developmental biology of the last common ancestor of the protostomes and deuterostomes, an animal from which >98% of all described animal species evolved.

A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia

Background

Homeobox genes are a superclass of transcription factors with diverse developmental regulatory functions, which are found in plants, fungi and animals. In animals, several Antennapedia (ANTP)-class homeobox genes reside in extremely ancient gene clusters (for example, the Hox, ParaHox, and NKL clusters) and the evolution of these clusters has been implicated in the morphological diversification of animal bodyplans. By contrast, similarly ancient gene clusters have not been reported among the other classes of homeobox genes (that is, the LIM, POU, PRD and SIX classes).

Results

Using a combination of in silico queries and phylogenetic analyses, we found that a cluster of three PRD-class homeobox genes (Homeobrain (hbn), Rax (rx) and Orthopedia (otp)) is present in cnidarians, insects and mollusks (a partial cluster comprising hbn and rx is present in the placozoan Trichoplax adhaerens). We failed to identify this 'HRO' cluster in deuterostomes; in fact, the Homeobrain gene appears to be missing from the chordate genomes we examined, although it is present in hemichordates and echinoderms. To illuminate the ancestral organization and function of this ancient cluster, we mapped the constituent genes against the assembled genome of a model cnidarian, the sea anemone Nematostella vectensis, and characterized their spatiotemporal expression using in situ hybridization. In N. vectensis, these genes reside in a span of 33 kb with the same gene order as previously reported in insects. Comparisons of genomic sequences and expressed sequence tags revealed the presence of alternative transcripts of Nv-otp and two highly unusual protein-coding polymorphisms in the terminal helix of the Nv-rx homeodomain. A population genetic survey revealed the Rx polymorphisms to be widespread in natural populations. During larval development, all three genes are expressed in the ectoderm, in non-overlapping territories along the oral-aboral axis, with distinct temporal expression.

Conclusion

We report the first evidence for a PRD-class homeobox cluster that appears to have been conserved since the time of the cnidarian-bilaterian ancestor, and possibly even earlier, given the presence of a partial cluster in the placozoan Trichoplax. Very similar clusters comprising these three genes exist in Nematostella and diverse protostomes. Interestingly, in chordates, one member of the ancestral cluster (homeobrain) has apparently been lost, and there is no linkage between rx and orthopedia in any of the vertebrates. In Nematostella, the spatial expression of these three genes along the body column is not colinear with their physical order in the cluster but the temporal expression is, therefore, using the terminology that has been applied to the Hox cluster genes, the HRO cluster would appear to exhibit temporal but not spatial colinearity. It remains to be seen whether the mechanisms responsible for the evolutionary conservation of the HRO cluster are the same mechanisms responsible for cohesion of the Hox cluster and other ANTP-class homeobox clusters that have been widely conserved throughout animal evolution.