Origin of Pax and Six gene families in sponges: Single PaxB and Six1/2 orthologs in Chalinula loosanoffi

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

It takes Two: Discovery of Spider Pax2 Duplicates Indicates Prominent Role in Chelicerate Central Nervous System, Eye, as Well as External Sense Organ Precursor Formation and Diversification After Neo- and Subfunctionalization

Paired box genes are conserved across animals and encode transcription factors playing key roles in development, especially neurogenesis. Pax6 is a chief example for functional conservation required for eye development in most bilaterian lineages except chelicerates. Pax6 is ancestrally linked and was shown to have interchangeable functions with Pax2. Drosophila melanogaster Pax2 plays an important role in the development of sensory hairs across the whole body. In addition, it is required for the differentiation of compound eyes, making it a prime candidate to study the genetic basis of arthropod sense organ development and diversification, as well as the role of Pax genes in eye development. Interestingly, in previous studies identification of chelicerate Pax2 was either neglected or failed. Here we report the expression of two Pax2 orthologs in the common house spider Parasteatoda tepidariorum, a model organism for chelicerate development. The two Pax2 orthologs most likely arose as a consequence of a whole genome duplication in the last common ancestor of spiders and scorpions. Pax2.1 is expressed in the peripheral nervous system, including developing lateral eyes and external sensilla, as well as the ventral neuroectoderm of P. tepidariorum embryos. This not only hints at a conserved dual role of Pax2/5/8 orthologs in arthropod sense organ development but suggests that in chelicerates, Pax2 could have acquired the role usually played by Pax6. For the other paralog, Pt-Pax2.2, expression was detected in the brain, but not in the lateral eyes and the expression pattern associated with sensory hairs differs in timing, pattern, and strength. To achieve a broader phylogenetic sampling, we also studied the expression of both Pax2 genes in the haplogyne cellar spider Pholcus phalangioides. We found that the expression difference between paralogs is even more extreme in this species, since Pp-Pax2.2 shows an interesting expression pattern in the ventral neuroectoderm while the expression in the prosomal appendages is strictly mesodermal. This expression divergence indicates both sub- and neofunctionalization after Pax2 duplication in spiders and thus presents an opportunity to study the evolution of functional divergence after gene duplication and its impact on sense organ diversification.

Secreted frizzled related protein is a target of PaxB and plays a role in aquiferous system development in the freshwater sponge, Ephydatia muelleri

Canonical and non-canonical Wnt signaling, as well as the Pax/Six gene network, are involved in patterning the freshwater sponge aquiferous system. Using computational approaches to identify transcription factor binding motifs in a freshwater sponge genome, we located putative PaxB binding sites near a Secreted Frizzled Related Protein (SFRP) gene in Ephydatia muelleri. EmSFRP is expressed throughout development, but with highest levels in juvenile sponges. In situ hybridization and antibody staining show EmSFRP expression throughout the pinacoderm and choanoderm in a subpopulation of amoeboid cells that may be differentiating archeocytes. Knockdown of EmSFRP leads to ectopic oscula formation during development, suggesting that EmSFRP acts as an antagonist of Wnt signaling in E. muelleri. Our findings support a hypothesis that regulation of the Wnt pathway by the Pax/Six network as well as the role of Wnt signaling in body plan morphogenesis was established before sponges diverged from the rest of the metazoans.

PaxA, but not PaxC, is required for cnidocyte development in the sea anemone Nematostella vectensis

Background

Pax genes are a family of conserved transcription factors that regulate many aspects of developmental morphogenesis, notably the development of ectodermal sensory structures including eyes. Nematostella vectensis, the starlet sea anemone, has numerous Pax orthologs, many of which are expressed early during embryogenesis. The function of Pax genes in this eyeless cnidarian is unknown.

Results

Here, we show that PaxA, but not PaxC, plays a critical role in the development of cnidocytes in N. vectensis. Knockdown of PaxA results in a loss of developing cnidocytes and downregulation of numerous cnidocyte-specific genes, including a variant of the transcription factor Mef2. We also demonstrate that the co-expression of Mef2 in a subset of the PaxA-expressing cells is associated with the development with a second lineage of cnidocytes and show that knockdown of the neural progenitor gene SoxB2 results in downregulation of expression of PaxA, Mef2, and several cnidocyte-specific genes. Because PaxA is not co-expressed with SoxB2 at any time during cnidocyte development, we propose a simple model for cnidogenesis whereby a SoxB2-expressing progenitor cell population undergoes division to give rise to PaxA-expressing cnidocytes, some of which also express Mef2.

Discussion

The role of PaxA in cnidocyte development among hydrozoans has not been studied, but the conserved role of SoxB2 in regulating the fate of a progenitor cell that gives rise to neurons and cnidocytes in Nematostella and Hydractinia echinata suggests that this SoxB2/PaxA pathway may well be conserved across cnidarians.

Ancestral and novel roles of Pax family genes in mollusks

Background

Pax genes are transcription factors with significant roles in cell fate specification and tissue differentiation during animal ontogeny. Most information on their tempo-spatial mode of expression is available from well-studied model organisms where the Pax-subfamilies Pax2/5/8, Pax6, and Paxα/β are mainly involved in the development of the central nervous system (CNS), the eyes, and other sensory organs. In certain taxa, Pax2/5/8 seems to be additionally involved in the development of excretion organs. Data on expression patterns in lophotrochozoans, and in particular in mollusks, are very scarce for all the above-mentioned Pax-subfamilies, which hampers reconstruction of their putative ancestral roles in bilaterian animals. Thus, we studied the developmental expression of Pax2/5/8, Pax6, and the lophotrochozoan-specific Paxβ in the worm-shaped mollusk Wirenia argentea, a member of Aplacophora that together with Polyplacophora forms the Aculifera, the proposed sister taxon to all primarily single-shelled mollusks (Conchifera).

Results

All investigated Pax genes are expressed in the developing cerebral ganglia and in the ventral nerve cords, but not in the lateral nerve cords of the tetraneural nervous system. Additionally, Pax2/5/8 is expressed in epidermal spicule-secreting or associated cells of the larval trunk and in the region of the developing protonephridia. We found no indication for an involvement of the investigated Pax genes in the development of larval or adult sensory organs of Wirenia argentea.

Conclusions

Pax2/5/8 seems to have a conserved role in the development of the CNS, whereas expression in the spicule-secreting tissues of aplacophorans and polyplacophorans suggests co-option in aculiferan skeletogenesis. The Pax6 expression pattern in Aculifera largely resembles the common bilaterian expression during CNS development. All data available on Paxβ expression argue for a common role in lophotrochozoan neurogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0919-x) contains supplementary material, which is available to authorized users.

The Pax gene family: Highlights from cephalopods

Pax genes play important roles in Metazoan development. Their evolution has been extensively studied but Lophotrochozoa are usually omitted. We addressed the question of Pax paralog diversity in Lophotrochozoa by a thorough review of available databases. The existence of six Pax families (Pax1/9, Pax2/5/8, Pax3/7, Pax4/6, Paxβ, PoxNeuro) was confirmed and the lophotrochozoan Paxβ subfamily was further characterized. Contrary to the pattern reported in chordates, the Pax2/5/8 family is devoid of homeodomain in Lophotrochozoa. Expression patterns of the three main pax classes (pax2/5/8, pax3/7, pax4/6) during Sepia officinalis development showed that Pax roles taken as ancestral and common in metazoans are modified in S. officinalis, most likely due to either the morphological specificities of cephalopods or to their direct development. Some expected expression patterns were missing (e.g. pax6 in the developing retina), and some expressions in unexpected tissues have been found (e.g. pax2/5/8 in dermal tissue and in gills). This study underlines the diversity and functional plasticity of Pax genes and illustrates the difficulty of using probable gene homology as strict indicator of homology between biological structures.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.

Comparative analyses of developmental transcription factor repertoires in sponges reveal unexpected complexity of the earliest animals

Developmental transcription factors (DTFs) control development of animals by affecting expression of target genes, some of which are transcription factors themselves. In bilaterians and cnidarians, conserved DTFs are involved in homologous processes such as gastrulation or specification of neurons. The genome of Amphimedon queenslandica, the first sponge to be sequenced, revealed that only a fraction of these conserved DTF families are present in demosponges. This finding was in line with the view that morphological complexity in the animal lineage correlates with developmental toolkit complexity. However, as the phylum Porifera is very diverse, Amphimedon's genome may not be representative of all sponges. The recently sequenced genomes of calcareous sponges Sycon ciliatum and Leucosolenia complicata allowed investigations of DTFs in a sponge lineage evolutionarily distant from demosponges. Surprisingly, the phylogenetic analyses of identified DTFs revealed striking differences between the calcareous sponges and Amphimedon. As these differences appear to be a result of independent gene loss events in the two sponge lineages, the last common ancestor of sponges had to possess a much more diverse repertoire of DTFs than extant sponges. Developmental expression of sponge homologs of genes involved in specification of the Bilaterian endomesoderm and the neurosensory cells suggests that roles of many DTFs date back to the last common ancestor of all animals. Strikingly, even DTFs displaying apparent pan-metazoan conservation of sequence and function are not immune to being lost from individual species genomes. The quest for a comprehensive picture of the developmental toolkit in the last common metazoan ancestor is thus greatly benefitting from the increasing accessibility of sequencing, allowing comparisons of multiple genomes within each phylum.

Gene Expression Data from the Moon Jelly, Aurelia, Provide Insights into the Evolution of the Combinatorial Code Controlling Animal Sense Organ Development

In Bilateria, Pax6, Six, Eya and Dach families of transcription factors underlie the development and evolution of morphologically and phyletically distinct eyes, including the compound eyes in Drosophila and the camera-type eyes in vertebrates, indicating that bilaterian eyes evolved under the strong influence of ancestral developmental gene regulation. However the conservation in eye developmental genetics deeper in the Eumetazoa, and the origin of the conserved gene regulatory apparatus controlling eye development remain unclear due to limited comparative developmental data from Cnidaria. Here we show in the eye-bearing scyphozoan cnidarian Aurelia that the ectodermal photosensory domain of the developing medusa sensory structure known as the rhopalium expresses sine oculis (so)/six1/2 and eyes absent/eya, but not optix/six3/6 or pax (A&B). In addition, the so and eya co-expression domain encompasses the region of active cell proliferation, neurogenesis, and mechanoreceptor development in rhopalia. Consistent with the role of so and eya in rhopalial development, developmental transcriptome data across Aurelia life cycle stages show upregulation of so and eya, but not optix or pax (A&B), during medusa formation. Moreover, pax6 and dach are absent in the Aurelia genome, and thus are not required for eye development in Aurelia. Our data are consistent with so and eya, but not optix, pax or dach, having conserved functions in sensory structure specification across Eumetazoa. The lability of developmental components including Pax genes relative to so-eya is consistent with a model of sense organ development and evolution that involved the lineage specific modification of a combinatorial code that specifies animal sense organs.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.

Phylogenetic analysis and expression patterns of Pax genes in the onychophoran Euperipatoides rowelli reveal a novel bilaterian Pax subfamily

Pax family genes encode a class of transcription factors that regulate various developmental processes. To shed light on the evolutionary history of these genes in Panarthropoda (Onychophora + Tardigrada + Arthropoda), we analyzed the Pax repertoire in the embryonic and adult transcriptomes of the onychophoran Euperipatoides rowelli. Our data revealed homologs of all five major bilaterian Pax subfamilies in this species, including Pax2/5/8, Pax4/6, Pox-neuro, Pax1/9/Pox-meso, and Pax3/7. In addition, we identified a new Pax member, pax-α, which does not fall into any other known Pax subfamily but instead clusters in the heterogenic Pax-α/β clade containing deuterostome, ecdysozoan, and lophotrochozoan gene sequences. These findings suggest that the last common bilaterian ancestor possessed six rather than five Pax genes, which have been retained in the panarthropod lineage. The expression data of Pax orthologs in the onychophoran embryo revealed distinctive patterns, some of which might be related to their ancestral roles in the last common panarthropod ancestor, whereas others might be specific to the onychophoran lineage. The derived roles include, for example, an involvement of pax2/5/8, pox-neuro, and pax3/7 in onychophoran nephridiogenesis, and an additional function of pax2/5/8 in the formation of the ventral and preventral organs. Furthermore, our transcriptomic analyses suggest that at least some Pax genes, including pax6 and pax-α, are expressed in the adult onychophoran head, although the corresponding functions remain to be clarified. The remarkable diversity of the Pax expression patterns highlights the functional and evolutionary plasticity of these genes in panarthropods.

Vertebrate cranial placodes as evolutionary innovations--the ancestor's tale

Evolutionary innovations often arise by tinkering with preexisting components building new regulatory networks by the rewiring of old parts. The cranial placodes of vertebrates, ectodermal thickenings that give rise to many of the cranial sense organs (ear, nose, lateral line) and ganglia, originated as such novel structures, when vertebrate ancestors elaborated their head in support of a more active and exploratory life style. This review addresses the question of how cranial placodes evolved by tinkering with ectodermal patterning mechanisms and sensory and neurosecretory cell types that have their own evolutionary history. With phylogenetic relationships among the major branches of metazoans now relatively well established, a comparative approach is used to infer, which structures evolved in which lineages and allows us to trace the origin of placodes and their components back from ancestor to ancestor. Some of the core networks of ectodermal patterning and sensory and neurosecretory differentiation were already established in the common ancestor of cnidarians and bilaterians and were greatly elaborated in the bilaterian ancestor (with BMP- and Wnt-dependent patterning of dorsoventral and anteroposterior ectoderm and multiple neurosecretory and sensory cell types). Rostral and caudal protoplacodal domains, giving rise to some neurosecretory and sensory cells, were then established in the ectoderm of the chordate and tunicate-vertebrate ancestor, respectively. However, proper cranial placodes as clusters of proliferating progenitors producing high-density arrays of neurosecretory and sensory cells only evolved and diversified in the ancestors of vertebrates.

Conservation of transcription factor binding specificities across 600 million years of bilateria evolution

Divergent morphology of species has largely been ascribed to genetic differences in the tissue-specific expression of proteins, which could be achieved by divergence in cis-regulatory elements or by altering the binding specificity of transcription factors (TFs). The relative importance of the latter has been difficult to assess, as previous systematic analyses of TF binding specificity have been performed using different methods in different species. To address this, we determined the binding specificities of 242 Drosophila TFs, and compared them to human and mouse data. This analysis revealed that TF binding specificities are highly conserved between Drosophila and mammals, and that for orthologous TFs, the similarity extends even to the level of very subtle dinucleotide binding preferences. The few human TFs with divergent specificities function in cell types not found in fruit flies, suggesting that evolution of TF specificities contributes to emergence of novel types of differentiated cells.

Pox neuro control of cell lineages that give rise to larval poly-innervated external sensory organs in Drosophila

The Pox neuro (Poxn) gene of Drosophila plays a crucial role in the development of poly-innervated external sensory (p-es) organs. However, how Poxn exerts this role has remained elusive. In this study, we have analyzed the cell lineages of all larval p-es organs, namely of the kölbchen, papilla 6, and hair 3. Surprisingly, these lineages are distinct from any previously reported cell lineages of sensory organs. Unlike the well-established lineage of mono-innervated external sensory (m-es) organs and a previously proposed model of the p-es lineage, we demonstrate that all wild-type p-es lineages exhibit the following features: the secondary precursor, pIIa, gives rise to all three support cells-socket, shaft, and sheath, whereas the other secondary precursor, pIIb, is neuronal and gives rise to all neurons. We further show that in one of the p-es lineages, that of papilla 6, one cell undergoes apoptosis. By contrast in Poxn null mutants, all p-es lineages have a reduced number of cells and their pattern of cell divisions is changed to that of an m-es organ, with the exception of a lineage in a minority of mutant kölbchen that retains a second bipolar neuron. Indeed, the role of Poxn in p-es lineages is consistent with the specification of the developmental potential of secondary precursors and the regulation of cell division but not apoptosis.

Evolution of the Pax-Six-Eya-Dach network: the calcisponge case study

Background

The Pax-Six-Eya-Dach network (PSEDN) is involved in a variety of developmental processes, including well documented roles in determination of sensory organs and morphogenesis in bilaterian animals. Expression of PSEDN components in cnidarians is consistent with function in sensory organ development. Recent work in demosponges demonstrated the presence of single homologs of Pax and Six genes, and their possible involvement in morphogenesis, but the absence of the remaining network components. Calcisponges are evolutionarily distant from demosponges, and the developmental toolkits of these two lineages differ significantly. We used an emerging model system, Sycon ciliatum, to identify components of the PSEDN and study their expression during embryonic and postembryonic development.

Results

We identified two Pax, three Six and one Eya genes in calcisponges, a situation strikingly different than in the previously studied demosponges. One of the calcisponge Pax genes can be identified as PaxB, while the second Pax gene has no clear affiliation. The three calcisponge Six genes could not be confidently classified within any known family of Six genes. Expression analysis in adult S. ciliatum demonstrated that representatives of Pax, Six and Eya are expressed in patterns consistent with roles in morphogenesis of the choanocyte chambers. Distinct paralogues of Pax and Six genes were expressed early in the development of the putative larval sensory cells, the cruciform cells. While lack of known photo pigments in calcisponge genomes precludes formal assignment of function to the cruciform cells, we also show that they express additional eumetazoan genes involved in specification of sensory and neuronal cells: Elav and Msi.

Conclusions

Our results indicate that the role of a Pax-Six-Eya network in morphogenesis likely predates the animal divergence. In addition, Pax and Six, as well as Elav and Msi are expressed during differentiation of cruciform cells, which are good candidates for being sensory cells of the calcaronean sponge larvae.

Connecting ears to eye muscles: evolution of a 'simple' reflex arc

Developmental and evolutionary data from vertebrates are beginning to elucidate the origin of the sensorimotor pathway that links gravity and motion detection to image-stabilizing eye movements--the vestibulo-ocular reflex (VOR). Conserved transcription factors coordinate the development of the vertebrate ear into three functional sensory compartments (graviception/translational linear acceleration, angular acceleration and sound perception). These sensory components connect to specific populations of vestibular and auditory projection neurons in the dorsal hindbrain through undetermined molecular mechanisms. In contrast, a molecular basis for the patterning of the vestibular projection neurons is beginning to emerge. These are organized through the actions of rostrocaudally and dorsoventrally restricted transcription factors into a 'hodological mosaic' within which coherent and largely segregated subgroups are specified to project to different targets in the spinal cord and brain stem. A specific set of these regionally diverse vestibular projection neurons functions as the central element that transforms vestibular sensory signals generated by active and passive head and body movements into motor output through the extraocular muscles. The large dynamic range of motion-related sensory signals requires an organization of VOR pathways as parallel, frequency-tuned, hierarchical connections from the sensory periphery to the motor output. We suggest that eyes, ears and functional connections subserving the VOR are vertebrate novelties that evolved into a functionally coherent motor control system in an almost stereotypic organization across vertebrate taxa.

The evolution and function of the Pax/Six regulatory network in sponges

Examining the origins of highly conserved gene regulatory networks (GRNs) will inform our understanding of the evolution of animal body plans. Sponges are believed to be the most ancient extant metazoan lineage, and as such, hold clues about the evolution of genetic programs deployed in animal development. We used the emerging freshwater sponge model, Ephydatia muelleri, to study the evolutionary origins of the Pax/Six/Eya/Dac (PSED) GRN. Orthologs to Pax and Six family members are present in E. muelleri and are expressed in endothelial cells lining the canal system as well as cells in the choanoderm. Knockdown of EmPaxB and EmSix1/2 by RNAi resulted in defects to the canal systems. We further show that PaxB may be in a regulatory relationship with Six1/2 in E. muelleri, thus demonstrating that a component of the PSED network was present early in metazoan evolution.

A NGS approach to the encrusting Mediterranean sponge Crella elegans (Porifera, Demospongiae, Poecilosclerida): transcriptome sequencing, characterization and overview of the gene expression along three life cycle stages

Sponges can be dominant organisms in many marine and freshwater habitats where they play essential ecological roles. They also represent a key group to address important questions in early metazoan evolution. Recent approaches for improving knowledge on sponge biological and ecological functions as well as on animal evolution have focused on the genetic toolkits involved in ecological responses to environmental changes (biotic and abiotic), development and reproduction. These approaches are possible thanks to newly available, massive sequencing technologies-such as the Illumina platform, which facilitate genome and transcriptome sequencing in a cost-effective manner. Here we present the first NGS (next-generation sequencing) approach to understanding the life cycle of an encrusting marine sponge. For this we sequenced libraries of three different life cycle stages of the Mediterranean sponge Crella elegans and generated de novo transcriptome assemblies. Three assemblies were based on sponge tissue of a particular life cycle stage, including non-reproductive tissue, tissue with sperm cysts and tissue with larvae. The fourth assembly pooled the data from all three stages. By aggregating data from all the different life cycle stages we obtained a higher total number of contigs, contigs with blast hit and annotated contigs than from one stage-based assemblies. In that multi-stage assembly we obtained a larger number of the developmental regulatory genes known for metazoans than in any other assembly. We also advance the differential expression of selected genes in the three life cycle stages to explore the potential of RNA-seq for improving knowledge on functional processes along the sponge life cycle.

Evolution of new characters after whole genome duplications: insights from amphioxus

Additional copies of genes resulting from two whole genome duplications at the base of the vertebrates have been suggested as enabling the evolution of vertebrate-specific structures such as neural crest, a midbrain/hindbrain organizer and neurogenic placodes. These structures, however, did not evolve entirely de novo, but arose from tissues already present in an ancestral chordate. This review discusses the evolutionary history of co-option of old genes for new roles in vertebrate development as well as the relative contributions of changes in cis-regulation and in protein structure. Particular examples are the FoxD, FGF8/17/18 and Pax2/5/8 genes. Comparisons with invertebrate chordates (amphioxus and tunicates) paint a complex picture with co-option of genes into new structures occurring both after and before the whole genome duplications. In addition, while cis-regulatory changes are likely of primary importance in evolution of vertebrate-specific structures, changes in protein structure including alternative splicing are non-trivial.

Six1 and Six1 cofactor expression is altered during early skeletal muscle overload in mice

Six1 is a transcription factor that, along with cofactors (Eya1, Eya3, and Dach2), regulates skeletal muscle fiber-type and development. SIX1 (human) gene expression decreases after overload, but the time course of Six1 expression, if protein is affected, and if the response differs between muscles with differing phenotypes, is not known. Our purpose was to examine Six1 gene and protein expression and co-factor gene expression during the initiation of muscle overload, and determine if the muscle phenotype altered this response. The plantaris and soleus were functionally overloaded by synergistic ablation of the gastrocnemius, and Six1 gene and protein, and Six1 cofactor gene expression was measured. Six1 gene expression decreased at 1 day of overload 48 ± 9 and 47 ± 20 % (p < 0.01) in the plantaris and soleus. After 3 days of overload, Six1 protein expression increased 73 ± 17 and 168 ± 57 % in the plantaris and soleus (p < 0.05). After 1 day of overload, Dach2 gene expression decreased 56 ± 9 and 35 ± 3 % in both muscles (p < 0.001), while Eya1 decreased 33 ± 5 % only in the soleus (p < 0.01). Eya3 gene expression increased 127 ± 26 % (p < 0.05) and 76 ± 16 % (p < 0.05) in the plantaris and soleus, while Dach2 gene expression decreased 71 ± 4 % (p < 0.05) in the soleus after 3 days of overload. Six1 and Six1 co-factor expression is responsive to muscle overload in both fast and slow muscles. This indicates that this molecular program may affect overload adaptation regardless of muscle phenotype.

The function and developmental expression of alternatively spliced isoforms of amphioxus and Xenopus laevis Pax2/5/8 genes: revealing divergence at the invertebrate to vertebrate transition

Pax genes encode highly conserved transcription factors vital for metazoan development. Pax transcripts, particularly those in Group II (Pax2/5/8), are extensively alternatively spliced. This study compares the transcriptional activation capacity and developmental stage-specific expression of major isoforms of Group II Pax proteins in amphioxus (Branchiostoma floridae) and in Xenopus laevis. The comparison reveals considerable divergence of splice forms between the lineages, with the X. laevis Group II Pax genes (Pax2, Pax5, and Pax8) possessing a greater repertoire of regulated and functionally distinct splice forms than the single amphioxus gene (Pax2/5/8). Surprisingly, some apparently conserved splice forms are expressed at quite different levels during development in the two organisms and present different capacities to activate transcription. However, despite this divergence, the combinatorial transcriptional activation capacity of the isoforms present in early X. laevis and amphioxus development are broadly similar. This suggests that the some of the conserved functional roles, implied by the expression of Group II Pax genes in homologous tissues of amphioxus and X. laevis embryos, may depend upon the combination of isoforms expressed in a particular tissue at a particular time in development. Thus, during early development, the evolutionary constraint on the net effect of several isoforms co-expressed in a given tissue may be more strict than that on specific isoforms. This flexibility may facilitate the appearance of new exons and splicing patterns in the vertebrate duplicates, leading to isoforms with subtly distinct functions critical to the subsequent development of vertebrate-specific cell types and structures.

Chasing the urmetazoon: striking a blow for quality data?

The ever-lingering question: "What did the urmetazoan look like?" has not lost its charm, appeal or elusiveness for one and a half centuries. A solid amount of organismal data give what some feel is a clear answer (e.g. Placozoa are at the base of the metazoan tree of life (ToL)), but a diversity of modern molecular data gives almost as many answers as there are exemplars, and even the largest molecular data sets could not solve the question and sometimes even suggest obvious zoological nonsense. Since the problems involved in this phylogenetic conundrum encompass a wide array of analytical freedom and uncertainty it seems questionable whether a further increase in molecular data (quantity) can solve this classical deep phylogeny problem. This review thus strikes a blow for evaluating quality data (including morphological, molecule morphologies, gene arrangement, and gene loss versus gene gain data) in an appropriate manner.

Regulation of Six1 expression by evolutionarily conserved enhancers in tetrapods

The Six1 homeobox gene plays critical roles in vertebrate organogenesis. Mice deficient for Six1 show severe defects in organs such as skeletal muscle, kidney, thymus, sensory organs and ganglia derived from cranial placodes, and mutations in human SIX1 cause branchio-oto-renal syndrome, an autosomal dominant developmental disorder characterized by hearing loss and branchial defects. The present study was designed to identify enhancers responsible for the dynamic expression pattern of Six1 during mouse embryogenesis. The results showed distinct enhancer activities of seven conserved non-coding sequences (CNSs) retained in tetrapod Six1 loci. The activities were detected in all cranial placodes (excluding the lens placode), dorsal root ganglia, somites, nephrogenic cord, notochord and cranial mesoderm. The major Six1-expression domains during development were covered by the sum of activities of these enhancers, together with the previously identified enhancer for the pre-placodal region and foregut endoderm. Thus, the eight CNSs identified in a series of our study represent major evolutionarily conserved enhancers responsible for the expression of Six1 in tetrapods. The results also confirmed that chick electroporation is a robust means to decipher regulatory information stored in vertebrate genomes. Mutational analysis of the most conserved placode-specific enhancer, Six1-21, indicated that the enhancer integrates a variety of inputs from Sox, Pax, Fox, Six, Wnt/Lef1 and basic helix-loop-helix proteins. Positive autoregulation of Six1 is achieved through the regulation of Six protein-binding sites. The identified Six1 enhancers provide valuable tools to understand the mechanism of Six1 regulation and to manipulate gene expression in the developing embryo, particularly in the sensory organs.

Cellular and molecular processes leading to embryo formation in sponges: evidences for high conservation of processes throughout animal evolution

The emergence of multicellularity is regarded as one of the major evolutionary events of life. This transition unicellularity/pluricellularity was acquired independently several times (King 2004). The acquisition of multicellularity implies the emergence of cellular cohesion and means of communication, as well as molecular mechanisms enabling the control of morphogenesis and body plan patterning. Some of these molecular tools seem to have predated the acquisition of multicellularity while others are regarded as the acquisition of specific lineages. Morphogenesis consists in the spatial migration of cells or cell layers during embryonic development, metamorphosis, asexual reproduction, growth, and regeneration, resulting in the formation and patterning of a body. In this paper, our aim is to review what is currently known concerning basal metazoans--sponges' morphogenesis from the tissular, cellular, and molecular points of view--and what remains to elucidate. Our review attempts to show that morphogenetic processes found in sponges are as diverse and complex as those found in other animals. In true epithelial sponges (Homoscleromorpha), as well as in others, we find similar cell/layer movements, cellular shape changes involved in major morphogenetic processes such as embryogenesis or larval metamorphosis. Thus, sponges can provide information enabling us to better understand early animal evolution at the molecular level but also at the cell/cell layer level. Indeed, comparison of molecular tools will only be of value if accompanied by functional data and expression studies during morphogenetic processes.

The physiology and molecular biology of sponge tissues

Sponges have become the focus of studies on molecular evolution and the evolution of animal body plans due to their ancient branching point in the metazoan lineage. Whereas our former understanding of sponge function was largely based on a morphological perspective, the recent availability of the first full genome of a sponge (Amphimedon queenslandica), and of the transcriptomes of other sponges, provides a new way of understanding sponges by their molecular components. This wealth of genetic information not only confirms some long-held ideas about sponge form and function but also poses new puzzles. For example, the Amphimedon sponge genome tells us that sponges possess a repertoire of genes involved in control of cell proliferation and in regulation of development. In vitro expression studies with genes involved in stem cell maintenance confirm that archaeocytes are the main stem cell population and are able to differentiate into many cell types in the sponge including pinacocytes and choanocytes. Therefore, the diverse roles of archaeocytes imply differential gene expression within a single cell ontogenetically, and gene expression is likely also different in different species; but what triggers cells to enter one pathway and not another and how each archaeocyte cell type can be identified based on this gene knowledge are new challenges. Whereas molecular data provide a powerful new tool for interpreting sponge form and function, because sponges are suspension feeders, their body plan and physiology are very much dependent on their physical environment, and in particular on flow. Therefore, in order to integrate new knowledge of molecular data into a better understanding the sponge body plan, it is important to use an organismal approach. In this chapter, we give an account of sponge body organization as it relates to the physiology of the sponge in light of new molecular data. We focus, in particular, on the structure of sponge tissues and review descriptive as well as experimental work on choanocyte morphology and function. Special attention is given to pinacocyte epithelia, cell junctions, and the molecules present in sponge epithelia. Studies describing the role of the pinacoderm in sensing, coordination, and secretion are reviewed. A wealth of recent work describes gene presence and expression patterns in sponge tissues during development, and we review this in the context of the previous descriptions of sponge morphology and physiology. A final section addresses recent findings of genes involved in the immune response. This review is far from exhaustive but intends rather to revisit for non-specialists key aspects of sponge morphology and physiology in light of new molecular data as a means to better understand and interpret sponge form and function today.

Deep phylogeny and evolution of sponges (phylum Porifera)

Sponges (phylum Porifera) are a diverse taxon of benthic aquatic animals of great ecological, commercial, and biopharmaceutical importance. They are arguably the earliest-branching metazoan taxon, and therefore, they have great significance in the reconstruction of early metazoan evolution. Yet, the phylogeny and systematics of sponges are to some extent still unresolved, and there is an on-going debate about the exact branching pattern of their main clades and their relationships to the other non-bilaterian animals. Here, we review the current state of the deep phylogeny of sponges. Several studies have suggested that sponges are paraphyletic. However, based on recent phylogenomic analyses, we suggest that the phylum Porifera could well be monophyletic, in accordance with cladistic analyses based on morphology. This finding has many implications for the evolutionary interpretation of early animal traits and sponge development. We further review the contribution that mitochondrial genes and genomes have made to sponge phylogenetics and explore the current state of the molecular phylogenies of the four main sponge lineages (Classes), that is, Demospongiae, Hexactinellida, Calcarea, and Homoscleromorpha, in detail. While classical systematic systems are largely congruent with molecular phylogenies in the class Hexactinellida and in certain parts of Demospongiae and Homoscleromorpha, the high degree of incongruence in the class Calcarea still represents a challenge. We highlight future areas of research to fill existing gaps in our knowledge. By reviewing sponge development in an evolutionary and phylogenetic context, we support previous suggestions that sponge larvae share traits and complexity with eumetazoans and that the simple sedentary adult lifestyle of sponges probably reflects some degree of secondary simplification. In summary, while deep sponge phylogenetics has made many advances in the past years, considerable efforts are still required to achieve a comprehensive understanding of the relationships among and within the main sponge lineages to fully appreciate the evolution of this extraordinary metazoan phylum.

Expansion, diversification, and expression of T-box family genes in Porifera

Sponges are among the earliest diverging lineage within the metazoan phyla. Although their adult morphology is distinctive, at several stages of development, they possess characteristics found in more complex animals. The T-box family of transcription factors is an evolutionarily ancient gene family known to be involved in the development of structures derived from all germ layers in the bilaterian animals. There is an incomplete understanding of the role that T-box transcription factors play in normal sponge development or whether developmental pathways using the T-box family share similarities between parazoan and eumetazoan animals. To address these questions, we present data that identify several important T-box genes in marine and freshwater sponges, place these genes in a phylogenetic context, and reveal patterns in how these genes are expressed in developing sponges. Phylogenetic analyses demonstrate that sponges have members of at least two of the five T-box subfamilies (Brachyury and Tbx2/3/4/5) and that the T-box genes expanded and diverged in the poriferan lineage. Our analysis of signature residues in the sponge T-box genes calls into question whether "true" Brachyury genes are found in the Porifera. Expression for a subset of the T-box genes was elucidated in larvae from the marine demosponge, Halichondria bowerbanki. Our results show that sponges regulate the timing and specificity of gene expression for T-box orthologs across larval developmental stages. In situ hybridization reveals distinct, yet sometimes overlapping expression of particular T-box genes in free-swimming larvae. Our results provide a comparative framework from which we can gain insights into the evolution of developmentally important pathways.

Flexibly deployed Pax genes in eye development at the early evolution of animals demonstrated by studies on a hydrozoan jellyfish

Pax transcription factors are involved in a variety of developmental processes in bilaterians, including eye development, a role typically assigned to Pax-6. Although no true Pax-6 gene has been found in nonbilateral animals, some jellyfish have eyes with complex structures. In the cubozoan jellyfish Tripedalia, Pax-B, an ortholog of vertebrate Pax-2/5/8, had been proposed as a regulator of eye development. Here we have isolated three Pax genes (Pax-A, Pax-B, and Pax-E) from Cladonema radiatum, a hydrozoan jellyfish with elaborate eyes. Cladonema Pax-A is strongly expressed in the retina, whereas Pax-B and Pax-E are highly expressed in the manubrium, the feeding and reproductive organ. Misexpression of Cladonema Pax-A induces ectopic eyes in Drosophila imaginal discs, whereas Pax-B and Pax-E do not. Furthermore, Cladonema Pax-A paired domain protein directly binds to the 5' upstream region of eye-specific Cladonema opsin genes, whereas Pax-B does not. Our data suggest that Pax-A, but not Pax-B or Pax-E, is involved in eye development and/or maintenance in Cladonema. Phylogenetic analysis indicates that Pax-6, Pax-B, and Pax-A belong to different Pax subfamilies, which diverged at the latest before the Cnidaria-Bilateria separation. We argue that our data, showing the involvement of Pax genes in hydrozoan eye development as in bilaterians, supports the monophyletic evolutionary origin of all animal eyes. We then propose that during the early evolution of animals, distinct classes of Pax genes, which may have played redundant roles at that time, were flexibly deployed for eye development in different animal lineages.