An atlas of spider development at single-cell resolution provides new insights into arthropod embryogenesis.

Spiders are a diverse order of chelicerates that diverged from other arthropods over 500 million years ago. Research on spider embryogenesis, particularly studies using the common house spider Parasteatoda tepidariorum, has made important contributions to understanding the evolution of animal development, including axis formation, segmentation, and patterning. However, we lack knowledge about the cells that build spider embryos, their gene expression profiles and fate. Single-cell transcriptomic analyses have been revolutionary in describing these complex landscapes of cellular genetics in a range of animals. Therefore, we carried out single-cell RNA sequencing of P. tepidariorum embryos at stages 7, 8 and 9, which encompass the establishment and patterning of the body plan, and initial differentiation of many tissues and organs. We identified 20 cell clusters, from 18.5 k cells, which were marked by many developmental toolkit genes, as well as a plethora of genes not previously investigated. We found differences in the cell cycle transcriptional signatures, suggestive of different proliferation dynamics, which related to distinctions between endodermal and some mesodermal clusters, compared with ectodermal clusters. We identified many Hox genes as markers of cell clusters, and Hox gene ohnologs were often present in different clusters. This provided additional evidence of sub- and/or neo-functionalisation of these important developmental genes after the whole genome duplication in an arachnopulmonate ancestor (spiders, scorpions, and related orders). We also examined the spatial expression of marker genes for each cluster to generate a comprehensive cell atlas of these embryonic stages. This revealed new insights into the cellular basis and genetic regulation of head patterning, hematopoiesis, limb development, gut development, and posterior segmentation. This atlas will serve as a platform for future analysis of spider cell specification and fate, and studying the evolution of these processes among animals at cellular resolution.

Annelid adult cell type diversity and their pluripotent cellular origins.

Many annelids can regenerate missing body parts or reproduce asexually, generating all cell types in adult stages. However, the putative adult stem cell populations involved in these processes, and the diversity of cell types generated by them, are still unknown. To address this, we recover 75,218 single cell transcriptomes of the highly regenerative and asexually-reproducing annelid Pristina leidyi. Our results uncover a rich cell type diversity including annelid specific types as well as novel types. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. We find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors, in piwi+ cells. Finally, lineage reconstruction analyses reveal computational differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids by single cell transcriptomics and suggest that a piwi+ cell population with a pluripotent stem cell signature is associated with adult cell type differentiation.

Exposure to Alexandrium spp. impairs the development of Green-lipped mussel (Perna canaliculus) embryos and larvae.

The green-lipped mussel (GLM) Perna canaliculus is an economically, ecologically, and culturally important species in Aotearoa New Zealand. Since 2011, harmful algal blooms (HABs) of Alexandrium spp. have occurred annually in the Marlborough Sounds, the largest GLM aquaculture region in New Zealand. Across a similar timeframe, there has been a severe reduction in wild spat (juvenile mussel) catch. This research investigated the effects of Alexandrium pacificum (which produces paralytic shellfish toxins; PSTs) and A. minutum (a non-producer of PSTs) on the development of four GLM larval life stages (gametes, embryos, D-stage and settlement). Early life stages of GLM were exposed to environmentally relevant concentrations of Alexandrium spp. as whole cell, lysate and filtrate treatments. A 48-h exposure of embryos to whole A. pacificum cells at 500 cells mL-1 caused lysis of embryos, severe abnormalities, and reduced development through to veliger (D-stage) larvae by 85%. GLM growth was impaired at cell concentrations as low as 250 cells mL-1 during a 4-day exposure of D-stage larvae to both Alexandrium spp. Exposure of GLM to both whole and lysed treatments of Alexandrium spp. at 500 cells mL-1 resulted in halved larval growth rates (2.00 µm day-1 vs 4.48 µm day-1 in the control) and growth remained impeded during a 4-day recovery period. Both A. pacificum and A. minutum were found to negatively impact D-larvae. Both whole-cell and lysed-cell treatments of A. pacificum had similar negative effects, suggesting that Alexandrium spp. toxicity to D-larvae is independent of PSTs. Additionally, cell membrane-free treatments of A. pacificum had no negative effects on embryo development, indicating that cell surface-associated bioactive compounds may be responsible for the observed negative effects during this early life stage. Conversely, non-PST-producing A. minutum was toxic to D-stage larvae but not to embryos; larval growth was reduced following a brief 1 h exposure of sperm to cell membrane-free treatments of A. pacificum. No effects were recorded in GLM larvae exposed during settlement, highlighting the potential for differences in susceptibility of early life stages to Alexandrium spp. exposure and the influence of exposure durations. In the wild, blooms of Alexandrium spp. can persist for several months, reaching cell densities higher than those investigated in the present study, and as such may be detrimental to the vulnerable early life stages of GLM.

Annelid adult cell type diversity and their pluripotent cellular origins.

Annelids are a broadly distributed, highly diverse, economically and environmentally important group of animals. Most species can regenerate missing body parts, and many are able to reproduce asexually. Therefore, many annelids can generate all adult cell types in adult stages. However, the putative adult stem cell populations involved in these processes, as well as the diversity of adult cell types generated by them, are still unknown. Here, we recover 75,218 single cell transcriptomes of Pristina leidyi, a highly regenerative and asexually-reproducing freshwater annelid. We characterise all major annelid adult cell types, and validate many of our observations by HCR in situ hybridisation. Our results uncover complex patterns of regionally expressed genes in the annelid gut, as well as neuronal, muscle and epidermal specific genes. We also characterise annelid-specific cell types such as the chaetal sacs and globin+ cells, and novel cell types of enigmatic affinity, including a vigilin+ cell type, a lumbrokinase+ cell type, and a diverse set of metabolic cells. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. In these piwi+ cells, we also find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors. Finally, lineage reconstruction analyses reveal the existence of differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids for the first time and serve as a resource for studying annelid cell types and their evolution. On the other hand, our characterisation of a piwi+ cell population with a pluripotent stem cell signature will serve as a platform for the study of annelid stem cells and their role in regeneration.

A Small Change With a Twist Ending: A Single Residue in EGF-CFC Drives Bilaterian Asymmetry.

Asymmetries are essential for proper organization and function of organ systems. Genetic studies in bilaterians have shown signaling through the Nodal/Smad2 pathway plays a key, conserved role in the establishment of body asymmetries. Although the main molecular players in the network for the establishment of left-right asymmetry (LRA) have been deeply described in deuterostomes, little is known about the regulation of Nodal signaling in spiralians. Here, we identified orthologs of the egf-cfc gene, a master regulator of the Nodal pathway in vertebrates, in several invertebrate species, which includes the first evidence of its presence in non-deuterostomes. Our functional experiments indicate that despite being present, egf-cfc does not play a role in the establishment of LRA in gastropods. However, experiments in zebrafish suggest that a single amino acid mutation in the egf-cfc gene in at least the common ancestor of chordates was the necessary step to induce a gain of function in LRA regulation. This study shows that the egf-cfc gene likely appeared in the ancestors of deuterostomes and "protostomes", before being adopted as a mechanism to regulate the Nodal pathway and the establishment of LRA in some lineages of deuterostomes.

Detecting Structural Variants and Associated Gene Presence-Absence Variation Phenomena in the Genomes of Marine Organisms.

As complete genomes become easier to attain, even from previously difficult-to-sequence species, and as genomic resequencing becomes more routine, it is becoming obvious that genomic structural variation is more widespread than originally thought and plays an important role in maintaining genetic variation in populations. Structural variants (SVs) and associated gene presence-absence variation (PAV) can be important players in local adaptation, allowing the maintenance of genetic variation and taking part in other evolutionarily relevant phenomena. While recent studies have highlighted the importance of structural variation in Mollusca, the prevalence of this phenomenon in the broader context of marine organisms remains to be fully investigated.Here, we describe a straightforward and broadly applicable method for the identification of SVs in fully assembled diploid genomes, leveraging the same reads used for assembly. We also explain a gene PAV analysis protocol, which could be broadly applied to any species with a fully sequenced reference genome available. Although the strength of these approaches have been tested and proven in marine invertebrates, which tend to have high levels of heterozygosity, possibly due to their lifestyle traits, they are also applicable to other species across the tree of life, providing a ready means to begin investigations into this potentially widespread phenomena.

Single individual structural variant detection uncovers widespread hemizygosity in molluscs.

The advent of complete genomic sequencing has opened a window into genomic phenomena obscured by fragmented assemblies. A good example of these is the existence of hemizygous regions of autosomal chromosomes, which can result in marked differences in gene content between individuals within species. While these hemizygous regions, and presence/absence variation of genes that can result, are well known in plants, firm evidence has only recently emerged for their existence in metazoans. Here, we use recently published, complete genomes from wild-caught molluscs to investigate the prevalence of hemizygosity across a well-known and ecologically important clade. We show that hemizygous regions are widespread in mollusc genomes, not clustered in individual chromosomes, and often contain genes linked to transposition, DNA repair and stress response. With targeted investigations of HSP70-12 and C1qDC, we also show how individual gene families are distributed within pan-genomes. This work suggests that extensive pan-genomes are widespread across the conchiferan Mollusca, and represent useful tools for genomic evolution, allowing the maintenance of additional genetic diversity within the population. As genomic sequencing and re-sequencing becomes more routine, the prevalence of hemizygosity, and its impact on selection and adaptation, are key targets for research across the tree of life. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

ACME dissociation: a versatile cell fixation-dissociation method for single-cell transcriptomics.

Single-cell sequencing technologies are revolutionizing biology, but they are limited by the need to dissociate live samples. Here, we present ACME (ACetic-MEthanol), a dissociation approach for single-cell transcriptomics that simultaneously fixes cells. ACME-dissociated cells have high RNA integrity, can be cryopreserved multiple times, and are sortable and permeable. As a proof of principle, we provide single-cell transcriptomic data of different species, using both droplet-based and combinatorial barcoding single-cell methods. ACME uses affordable reagents, can be done in most laboratories and even in the field, and thus will accelerate our knowledge of cell types across the tree of life.

Striking parallels between dorsoventral patterning in Drosophila and Gryllus reveal a complex evolutionary history behind a model gene regulatory network.

Dorsoventral pattering relies on Toll and BMP signalling in all insects studied so far, with variations in the relative contributions of both pathways. Drosophila and the beetle Tribolium share extensive dependence on Toll, while representatives of more distantly related lineages like the wasp Nasonia and bug Oncopeltus rely more strongly on BMP signalling. Here, we show that in the cricket Gryllus bimaculatus, an evolutionarily distant outgroup, Toll has, like in Drosophila, a direct patterning role for the ventral half of the embryo. In addition, Toll polarises BMP signalling, although this does not involve the conserved BMP inhibitor Sog/Chordin. Finally, Toll activation relies on ovarian patterning mechanisms with striking similarity to Drosophila. Our data suggest two surprising hypotheses: (1) that Toll's patterning function in Gryllus and Drosophila is the result of convergent evolution or (2) a Drosophila-like system arose early in insect evolution and was extensively altered in multiple independent lineages.

Mitochondrial evolution in the Demospongiae (Porifera): Phylogeny, divergence time, and genome biology.

The sponge class Demospongiae is the most speciose and morphologically diverse in the phylum Porifera, and the species within it are vital components of a range of ecosystems worldwide. Despite their ubiquity, a number of recalcitrant problems still remain to be solved regarding their phylogenetic inter-relationships, the timing of their appearance, and their mitochondrial biology, the latter of which is only beginning to be investigated. Here we generated 14 new demosponge mitochondrial genomes which, alongside previously published mitochondrial resources, were used to address these issues. In addition to phylogenomic analysis, we have used syntenic data and analysis of coding regions to forge a framework for understanding the inter-relationships between Demospongiae sub-classes and orders. We have also leveraged our new resources to study the mitochondrial biology of these clades in terms of codon usage, optimisation and gene expression, to understand how these vital cellular components may have contributed to the success of the Porifera. Our results strongly support a sister relationship between Keratosa and (Verongimorpha + Heteroscleromorpha), contradicting previous studies using nuclear markers. Our study includes one species of Clionaida, and show for the first time support for a grouping of Suberitida+(Clionaida+(Tethyida + Poecilosclerida). The findings of our phylogenetic analyses are supported by in-depth examination of structural and coding-level evidence from our mitochondrial data. A time-calibrated phylogeny estimated the origin of Demospongiae in the Cambrian (~529 Mya), and suggests that most demosponge order crown-groups emerged in the Mesozoic. This work therefore provides a robust basis for considering demosponge phylogenetic relationships, as well as essential mitochondrial data for understanding the biological basis for their success and diversity.

Sponge microbiome stability during environmental acquisition of highly specific photosymbionts.

In this study, we used in situ transplantations to provide the first evidence of horizontal acquisition of cyanobacterial symbionts by a marine sponge. The acquisition of the symbionts by the host sponge Petrosia ficiformis, which was observed in distinct visible patches, appeared several months after transplantation and at different times on different sponge specimens. We further used 16S rRNA gene amplicon sequencing of genomic DNA (gDNA) and complementary DNA (cDNA) and metatranscriptomics to investigate how the acquisition of the symbiotic cyanobacterium Candidatus Synechococcus feldmannii perturbed the diverse microbiota associated with the host P. ficiformis. To our surprise, the microbiota remained relatively stable during cyanobacterial symbiont acquisition at both structural (gDNA content) and activity (cDNA expression) levels. At the transcriptomic level, photosynthesis was the primary function gained following the acquisition of cyanobacteria. Genes involved in carotene production and oxidative stress tolerance were among those highly expressed by Ca. S. feldmannii, suggesting that this symbiont may protect itself and its host from damaging light radiation.

Tracing animal genomic evolution with the chromosomal-level assembly of the freshwater sponge Ephydatia muelleri.

The genomes of non-bilaterian metazoans are key to understanding the molecular basis of early animal evolution. However, a full comprehension of how animal-specific traits, such as nervous systems, arose is hindered by the scarcity and fragmented nature of genomes from key taxa, such as Porifera. Ephydatia muelleri is a freshwater sponge found across the northern hemisphere. Here, we present its 326 Mb genome, assembled to high contiguity (N50: 9.88 Mb) with 23 chromosomes on 24 scaffolds. Our analyses reveal a metazoan-typical genome architecture, with highly shared synteny across Metazoa, and suggest that adaptation to the extreme temperatures and conditions found in freshwater often involves gene duplication. The pancontinental distribution and ready laboratory culture of E. muelleri make this a highly practical model system which, with RNAseq, DNA methylation and bacterial amplicon data spanning its development and range, allows exploration of genomic changes both within sponges and in early animal evolution.

The gene-rich genome of the scallop Pecten maximus.

BACKGROUND: The king scallop, Pecten maximus, is distributed in shallow waters along the Atlantic coast of Europe. It forms the basis of a valuable commercial fishery and plays a key role in coastal ecosystems and food webs. Like other filter feeding bivalves it can accumulate potent phytotoxins, to which it has evolved some immunity. The molecular origins of this immunity are of interest to evolutionary biologists, pharmaceutical companies, and fisheries management. FINDINGS: Here we report the genome assembly of this species, conducted as part of the Wellcome Sanger 25 Genomes Project. This genome was assembled from PacBio reads and scaffolded with 10X Chromium and Hi-C data. Its 3,983 scaffolds have an N50 of 44.8 Mb (longest scaffold 60.1 Mb), with 92% of the assembly sequence contained in 19 scaffolds, corresponding to the 19 chromosomes found in this species. The total assembly spans 918.3 Mb and is the best-scaffolded marine bivalve genome published to date, exhibiting 95.5% recovery of the metazoan BUSCO set. Gene annotation resulted in 67,741 gene models. Analysis of gene content revealed large numbers of gene duplicates, as previously seen in bivalves, with little gene loss, in comparison with the sequenced genomes of other marine bivalve species. CONCLUSIONS: The genome assembly of P. maximus and its annotated gene set provide a high-quality platform for studies on such disparate topics as shell biomineralization, pigmentation, vision, and resistance to algal toxins. As a result of our findings we highlight the sodium channel gene Nav1, known to confer resistance to saxitoxin and tetrodotoxin, as a candidate for further studies investigating immunity to domoic acid.

Warm temperatures, cool sponges: the effect of increased temperatures on the Antarctic sponge Isodictya sp.

Although the cellular and molecular responses to exposure to relatively high temperatures (acute thermal stress or heat shock) have been studied previously, only sparse empirical evidence of how it affects cold-water species is available. As climate change becomes more pronounced in areas such as the Western Antarctic Peninsula, both long-term and occasional acute temperature rises will impact species found there, and it has become crucial to understand the capacity of these species to respond to such thermal stress. Here, we use the Antarctic sponge Isodictya sp. to investigate how sessile organisms (particularly Porifera) can adjust to acute short-term heat stress, by exposing this species to 3 and 5 °C for 4 h, corresponding to predicted temperatures under high-end 2080 IPCC-SRES scenarios. Assembling a de novo reference transcriptome (90,188 contigs, >93.7% metazoan BUSCO genes) we have begun to discern the molecular response employed by Isodictya to adjust to heat exposure. Our initial analyses suggest that TGF-ß, ubiquitin and hedgehog cascades are involved, alongside other genes. However, the degree and type of response changed little from 3 to 5 °C in the time frame examined, suggesting that even moderate rises in temperature could cause stress at the limits of this organism's capacity. Given the importance of sponges to Antarctic ecosystems, our findings are vital for discerning the consequences of short-term increases in Antarctic ocean temperature on these and other species.

The utility of micro-computed tomography for the non-destructive study of eye microstructure in snails.

Molluscan eyes exhibit an enormous range of morphological variation, ranging from tiny pigment-cup eyes in limpets, compound eyes in ark clams and pinhole eyes in Nautilus, through to concave mirror eyes in scallops and the large camera-type eyes of the more derived cephalopods. Here we assess the potential of non-destructive micro-computed tomography (µ-CT) for investigating the anatomy of molluscan eyes in three species of the family Solariellidae, a group of small, deep-sea gastropods. We compare our results directly with those from traditional histological methods applied to the same specimens, and show not only that eye microstructure can be visualised in sufficient detail for meaningful comparison even in very small animals, but also that µ-CT can provide additional insight into gross neuroanatomy without damaging rare and precious specimens. Data from µ-CT scans also show that neurological innervation of eyes is reduced in dark-adapted snails when compared with the innervation of cephalic tentacles, which are involved in mechanoreception and possibly chemoreception. Molecular tests also show that the use of µ-CT and phosphotungstic acid stain do not prevent successful downstream DNA extraction, PCR amplification or sequencing. The use of µ-CT methods is therefore highly recommended for the investigation of difficult-to-collect or unique specimens.

Symbiosis, Selection, and Novelty: Freshwater Adaptation in the Unique Sponges of Lake Baikal.

Freshwater sponges (Spongillida) are a unique lineage of demosponges that secondarily colonized lakes and rivers and are now found ubiquitously in these ecosystems. They developed specific adaptations to freshwater systems, including the ability to survive extreme thermal ranges, long-lasting dessication, anoxia, and resistance to a variety of pollutants. Although spongillids have colonized all freshwater systems, the family Lubomirskiidae is endemic to Lake Baikal and plays a range of key roles in this ecosystem. Our work compares the genomic content and microbiome of individuals of three species of the Lubomirskiidae, providing hypotheses for how molecular evolution has allowed them to adapt to their unique environments. We have sequenced deep (>92% of the metazoan "Benchmarking Universal Single-Copy Orthologs" [BUSCO] set) transcriptomes from three species of Lubomirskiidae and a draft genome resource for Lubomirskia baikalensis. We note Baikal sponges contain unicellular algal and bacterial symbionts, as well as the dinoflagellate Gyrodinium. We investigated molecular evolution, gene duplication, and novelty in freshwater sponges compared with marine lineages. Sixty one orthogroups have consilient evidence of positive selection. Transporters (e.g., zinc transporter-2), transcription factors (aristaless-related homeobox), and structural proteins (e.g. actin-3), alongside other genes, are under strong evolutionary pressure in freshwater, with duplication driving novelty across the Spongillida, but especially in the Lubomirskiidae. This addition to knowledge of freshwater sponge genetics provides a range of tools for understanding the molecular biology and, in the future, the ecology (e.g., colonization and migration patterns) of these key species.

Trimitomics: An efficient pipeline for mitochondrial assembly from transcriptomic reads in nonmodel species.

Mitochondrial resources are of known utility to many fields of phylogenetic, population and molecular biology. Their combination of faster and slower-evolving regions and high copy number enables them to be used in many situations where other loci are unsuitable, with degraded samples and after recent speciation events.The advent of next-generation sequencing technologies (and notably the Illumina platform) has led to an explosion in the number of samples that can be studied at transcriptomic level, at relatively low cost. Here we describe a robust pipeline for the recovery of mitochondrial genomes from these RNA-sequencing resources. This pipeline can be used on sequencing of a variety of depths, and reliably recovers the protein coding and ribosomal gene complements of mitochondria from almost any transcriptomic sequencing experiment. The complete sequence of the mitochondrial genome can also be recovered when sequencing is performed in sufficient depth. We show the efficacy of our pipeline using data from eight nonmodel invertebrates of six disparate phyla. Interestingly, among our poriferan data, where microbiological symbionts are known empirically to make mitochondrial assembly difficult, this pipeline proved especially useful. Our pipeline will allow the recovery of mitochondrial data from a variety of previously sequenced samples, and add an additional angle of enquiry to future RNA-sequencing efforts, simplifying the process of mitochondrial genome assembly for even the most recalcitrant clades and adding these data to the scientific record for a range of future uses.

Population substructure and signals of divergent adaptive selection despite admixture in the sponge Dendrilla antarctica from shallow waters surrounding the Antarctic Peninsula.

Antarctic shallow-water invertebrates are exceptional candidates to study population genetics and evolution, because of their peculiar evolutionary history and adaptation to extreme habitats that expand and retreat with the ice sheets. Among them, sponges are one of the major components, yet population connectivity of none of their many Antarctic species has been studied. To investigate gene flow, local adaptation and resilience to near-future changes caused by global warming, we sequenced 62 individuals of the sponge Dendrilla antarctica along the Western Antarctic Peninsula (WAP) and the South Shetlands (spanning ~900 km). We obtained information from 577 double digest restriction site-associated DNA sequencing (ddRADseq)-derived single nucleotide polymorphism (SNP), using RADseq techniques for the first time with shallow-water sponges. In contrast to other studies in sponges, our 389 neutral SNPs data set showed high levels of gene flow, with a subtle substructure driven by the circulation system of the studied area. However, the 140 outlier SNPs under positive selection showed signals of population differentiation, separating the central-southern WAP from the Bransfield Strait area, indicating a divergent selection process in the study area despite panmixia. Fourteen of these outliers were annotated, being mostly involved in immune and stress responses. We suggest that the main selective pressure on D. antarctica might be the difference in the planktonic communities present in the central-southern WAP compared to the Bransfield Strait area, ultimately depending on sea-ice control of phytoplankton blooms. Our study unveils an unexpectedly long-distance larval dispersal exceptional in Porifera, broadening the use of genome-wide markers within nonmodel Antarctic organisms.

Probing recalcitrant problems in polyclad evolution and systematics with novel mitochondrial genome resources.

For their apparent morphological simplicity, the Platyhelminthes or "flatworms" are a diverse clade found in a broad range of habitats. Their body plans have however made them difficult to robustly classify. Molecular evidence is only beginning to uncover the true evolutionary history of this clade. Here we present nine novel mitochondrial genomes from the still undersampled orders Polycladida and Rhabdocoela, assembled from short Illumina reads. In particular we present for the first time in the literature the mitochondrial sequence of a Rhabdocoel, Bothromesostoma personatum (Typhloplanidae, Mesostominae). The novel mitochondrial genomes examined generally contained the 36 genes expected in the Platyhelminthes, with all possessing 12 of the 13 protein-coding genes normally found in metazoan mitochondrial genomes (ATP8 being absent from all Platyhelminth mtDNA sequenced to date), along with two ribosomal RNA genes. The majority presented possess 22 transfer RNA genes, and a single tRNA gene was absent from two of the nine assembled genomes. By comparison of mitochondrial gene order and phylogenetic analysis of the protein coding and ribosomal RNA genes contained within these sequences with those of previously sequenced species we are able to gain a firm molecular phylogeny for the inter-relationships within this clade. Our phylogenetic reconstructions, using both nucleotide and amino acid sequences under several models and both Bayesian and Maximum Likelihood methods, strongly support the monophyly of Polycladida, and the monophyly of Acotylea and Cotylea within that clade. They also allow us to speculate on the early emergence of Macrostomida, the monophyly of a "Turbellarian-like" clade, the placement of Rhabditophora, and that of Platyhelminthes relative to the Lophotrochozoa (=Spiralia). The data presented here therefore represent a significant advance in our understanding of platyhelminth phylogeny, and will form the basis of a range of future research in the still-disputed classifications within this taxon.

Delegating Sex: Differential Gene Expression in Stolonizing Syllids Uncovers the Hormonal Control of Reproduction.

Stolonization in syllid annelids is a unique mode of reproduction among animals. During the breeding season, a structure resembling the adult but containing only gametes, called stolon, is formed generally at the posterior end of the animal. When stolons mature, they detach from the adult and gametes are released into the water column. The process is synchronized within each species, and it has been reported to be under environmental and endogenous control, probably via endocrine regulation. To further understand reproduction in syllids and to elucidate the molecular toolkit underlying stolonization, we generated Illumina RNA-seq data from different tissues of reproductive and nonreproductive individuals of Syllis magdalena and characterized gene expression during the stolonization process. Several genes involved in gametogenesis (ovochymase, vitellogenin, testis-specific serine/threonine-kinase), immune response (complement receptor 2), neuronal development (tyrosine-protein kinase Src42A), cell proliferation (alpha-1D adrenergic receptor), and steroid metabolism (hydroxysteroid dehydrogenase 2) were found differentially expressed in the different tissues and conditions analyzed. In addition, our findings suggest that several neurohormones, such as methyl farnesoate, dopamine, and serotonin, might trigger stolon formation, the correct maturation of gametes and the detachment of stolons when gametogenesis ends. The process seems to be under circadian control, as indicated by the expression patterns of r-opsins. Overall, our results shed light into the genes that orchestrate the onset of gamete formation and improve our understanding of how some hormones, previously reported to be involved in reproduction and metamorphosis processes in other invertebrates, seem to also regulate reproduction via stolonization.