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.

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.

Recycling resources: silica of diatom frustules as a source for spicule building in Antarctic siliceous demosponges

Dissolved silicon (DSi) is biologically processed to produce siliceous skeletons by a variety of organisms including radiolarians, silicoflagellates, choanoflagellates, plants, diatoms and some animals. In the photic ocean, diatoms are dominant consumers over competing other silicifiers. In Antarctica, where DSi is not particularly limiting, diatoms and sponges coexist in high abundances. Interestingly, diatom ingestion by sponges is a regular feeding strategy there. Although it was known that the diatom organic nutrients are readily metabolized by the sponges, what happened to the inorganic diatom silica skeleton remained unexplored. Here, we have conducted a multi-analytical approach to investigate the processing of diatom silica and whether it is reconverted into sponge silica. We have documented widespread diatom consumption by several demosponges, identifying storage vesicles for the diatom-derived silica by electron microscopy and microanalysis. Diatom-consuming sponges showed upregulation of silicatein and silicase genes, which in addition to the δ 30Si values of their silica, supports that the sponges are converting the ingested diatom silica into sponge silica without much further Si fractionation. Our multidisciplinary approach suggests that the reutilization of diatom silica by sponges is a common feature among Antarctic sponges, which should be further investigated in other latitudes and in other silicifiers.

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.

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.

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.

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.

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.

Integrative systematics of clathrinid sponges: morphological, reproductive and phylogenetic characterisation of a new species of Leucetta from Antarctica (Porifera, Calcarea, Calcinea) with notes on the occurrence of flagellated sperm

Our study reports on the occurrence of a new species of Leucetta (Calcinea, Calcarea) from the Southern Ocean, Leucetta giribeti, sp. nov., collected in the shallow waters (15m depth) of Deception Island, South Shetland Islands. This new taxon is described based on a combination of morphological and molecular data, including the description of oocytes, embryos, larvae and sperm found in the choanosome. While female reproductive elements showed great similarities with those of other calcineans, sperm is reported here for the first time in the whole Calcinea subclass. Sperm cells are flagellated and possess a typical spermatic mid-piece, which is usually observed in cnidarians. In our phylogenetic analyses, we recovered Leucetta giribeti, sp. nov. as sister species of a clade formed by species of the genera Leucetta, Pericharax and Leucettusa. Although the clade in which Leucetta giribeti, sp. nov. is placed is supported by molecular and morphological features, we cannot propose a new genus due to uncertainties regarding the type species of the genus, Leucetta primigenia Haeckel, 1872. Our study reinforces the relevance of integrative approaches in the description of new taxa and contributes to resolving the poorly known reproductive patterns of Antarctic sponge species.

Mitochondrial genome and polymorphic microsatellite markers from the abyssal sponge Plenaster craigi Lim & Wiklund, 2017: tools for understanding the impact of deep-sea mining

The abyssal demosponge Plenaster craigi is endemic to the Clarion - Clipperton Zone (CCZ) in the NE Pacific, a region with abundant seafloor polymetallic nodules and of potential interest for mining. Plenaster craigi encrusts on these nodules and is an abundant component of the ecosystem. To assess the impact of mining operations, it is crucial to understand the genetics of this species, because its genetic diversity and connectivity across the area may be representative of other nodule-encrusting invertebrate epifauna. Here we describe and characterize 14 polymorphic microsatellite markers from this keystone species using Illumina MiSeq, tested for 75 individuals from three different areas across the CCZ, including an Area of Particular Environmental Interest (APEI-6) and two areas within the adjacent UK1 mining exploration area. The number of alleles per locus ranged from 3 to 30 (13.33 average alleles for all loci across areas). Observed and expected heterozygosity ranged from 0.909-0.048 and from 0.954-0.255, respectively. Several loci displayed significant deviation from the Hardy-Weinberg equilibrium, which appears to be common in other sponge studies. The microsatellite loci described here will be used to assess the genetic structure and connectivity on populations of the sponge across the CCZ, which will be invaluable for monitoring the impact of mining operations on its habitat. Also, we provide the annotated mitochondrial genome of P. craigi, compare its arrangement with other closely related species, and discuss the phylogenetic framework for the sponge after Maximum Likelihood and Bayesian Inference analyses using nucleotide and amino acid sequences data sets separately.

A novel role for Ets4 in axis specification and cell migration in the spider Parasteatoda tepidariorum

Organizers play important roles during the embryonic development of many animals. The most famous example is the Spemann organizer that sets up embryonic axes in amphibian embryos. In spiders, a group of BMP secreting mesenchymal cells (the cumulus) functions as an organizer of the dorsoventral axis. Similar to experiments performed with the Spemann organizer, transplantation of the cumulus is able to induce a secondary axis in spiders. Despite the importance of this structure, it is unknown which factors are needed to activate cumulus specific gene expression. To address this question, we performed a transcriptomic analysis of early embryonic development in the spider Parasteatoda tepidariorum. Through this work, we found that the transcription factor Pt-Ets4 is needed for cumulus integrity, dorsoventral patterning and for the activation of Pt-hunchback and Pt-twist expression. Furthermore, ectopic expression of Pt-Ets4 is sufficient to induce cell delamination and migration by inducing a mesoderm-like cell fate.

Deep, Staged Transcriptomic Resources for the Novel Coleopteran Models Atrachya menetriesi and Callosobruchus maculatus

Despite recent efforts to sample broadly across metazoan and insect diversity, current sequence resources in the Coleoptera do not adequately describe the diversity of the clade. Here we present deep, staged transcriptomic data for two coleopteran species, Atrachya menetriesi (Faldermann 1835) and Callosobruchus maculatus (Fabricius 1775). Our sampling covered key stages in ovary and early embryonic development in each species. We utilized this data to build combined assemblies for each species which were then analysed in detail. The combined A. menetriesi assembly consists of 228,096 contigs with an N50 of 1,598 bp, while the combined C. maculatus assembly consists of 128,837 contigs with an N50 of 2,263 bp. For these assemblies, 34.6% and 32.4% of contigs were identified using Blast2GO, and 97% and 98.3% of the BUSCO set of metazoan orthologs were present, respectively. We also carried out manual annotation of developmental signalling pathways and found that nearly all expected genes were present in each transcriptome. Our analyses show that both transcriptomes are of high quality. Lastly, we performed read mapping utilising our timed, stage specific RNA samples to identify differentially expressed contigs. The resources presented here will provide a firm basis for a variety of experimentation, both in developmental biology and in comparative genomic studies.

Deep, multi-stage transcriptome of the schistosomiasis vector Biomphalaria glabrata provides platform for understanding molluscan disease-related pathways

BACKGROUND

The gastropod mollusc Biomphalaria glabrata is well known as a vector for the tropical disease schistosomiasis, which affects nearly 200 million people worldwide. Despite intensive study, our understanding of the genetic basis of B. glabrata development, growth and disease resistance is constrained by limited genetic resources, constraints for which next-generation sequencing methods provide a ready solution.

METHODS

Illumina sequencing and de novo assembly using the Trinity program was used to generate a high-quality transcriptomic dataset spanning the entirety of in ovo development in schistosomiasis-free B. glabrata. This was subjected to automated (KEGG, BLAST2GO) and manual annotation efforts, allowing insight into the gene complements of this species in a number of contexts.

RESULTS

Excellent dataset recovery was observed, with 133,084 contigs produced of mean size 2219.48 bp. 80,952 (60.8 %) returned a BLASTx hit with an E value of less than 10^-3, and 74,492 (55.97 %) were either mapped or assigned a GO identity using the BLAST2GO program. The CEGMA set of core eukaryotic genes was found to be 99.6 % present, indicating exceptional transcriptome completeness. We were able to identify a wealth of disease-pathway related genes within our dataset, including the Wnt, apoptosis and Notch pathways. This provides an invaluable reference point for further work into molluscan development and evolution, for studying the impact of schistosomiasis in this species, and perhaps providing targets for the treatment of this widespread disease.

CONCLUSIONS

Here we present a deep transcriptome of an embryonic sample of schistosomiasis-free B. glabrata, presenting a comprehensive dataset for comparison to disease-affected specimens and from which conclusions can be drawn about the genetics of this widespread medical model. Furthermore, the dataset provided by this sequencing provides a useful reference point for comparison to other mollusc species, which can be used to better understand the evolution of this commercially, ecologically and medically important phylum.

Ancestral whole-genome duplication in the marine chelicerate horseshoe crabs

Whole-genome duplication (WGD) results in new genomic resources that can be exploited by evolution for rewiring genetic regulatory networks in organisms. In metazoans, WGD occurred before the last common ancestor of vertebrates, and has been postulated as a major evolutionary force that contributed to their speciation and diversification of morphological structures. Here, we have sequenced genomes from three of the four extant species of horseshoe crabs-Carcinoscorpius rotundicauda, Limulus polyphemus and Tachypleus tridentatus. Phylogenetic and sequence analyses of their Hox and other homeobox genes, which encode crucial transcription factors and have been used as indicators of WGD in animals, strongly suggests that WGD happened before the last common ancestor of these marine chelicerates >135 million years ago. Signatures of subfunctionalisation of paralogues of Hox genes are revealed in the appendages of two species of horseshoe crabs. Further, residual homeobox pseudogenes are observed in the three lineages. The existence of WGD in the horseshoe crabs, noted for relative morphological stasis over geological time, suggests that genomic diversity need not always be reflected phenotypically, in contrast to the suggested situation in vertebrates. This study provides evidence of ancient WGD in the ecdysozoan lineage, and reveals new opportunities for studying genomic and regulatory evolution after WGD in the Metazoa.

The phylogenetic utility and functional constraint of microRNA flanking sequences

MicroRNAs (miRNAs) have recently risen to prominence as novel factors responsible for post-transcriptional regulation of gene expression. miRNA genes have been posited as highly conserved in the clades in which they exist. Consequently, miRNAs have been used as rare genome change characters to estimate phylogeny by tracking their gain and loss. However, their short length (21-23 bp) has limited their perceived utility in sequenced-based phylogenetic inference. Here, using reference taxa with established phylogenetic relationships, we demonstrate that miRNA sequences are of high utility in quantitative, rather than in qualitative, phylogenetic analysis. The clear orthology among miRNA genes from different species makes it straightforward to identify and align these sequences from even fragmentary datasets. We also identify significant sequence conservation in the regions directly flanking miRNA genes, and show that this too is of utility in phylogenetic analysis, as well as highlighting conserved regions that will be of interest to other fields. Employing miRNA sequences from 12 sequenced drosophilid genomes, together with a Tribolium castaneum outgroup, we demonstrate that this approach is robust using Bayesian and maximum-likelihood methods. The utility of these characters is further demonstrated in the rhabditid nematodes and primates. As next-generation sequencing makes it more cost-effective to sequence genomes and small RNA libraries, this methodology provides an alternative data source for phylogenetic analysis. The approach allows rapid resolution of relationships between both closely related and rapidly evolving species, and provides an additional tool for investigation of relationships within the tree of life.

Evolution, divergence and loss of the Nodal signalling pathway: new data and a synthesis across the Bilateria

Since the discovery that the TGF-β signalling molecule Nodal and its downstream effector Pitx have a parallel role in establishing asymmetry between molluscs and deuterostomes the debate over the degree to which this signalling pathway is conserved across the Bilateria as a whole has been ongoing. Further taxon sampling is critical to understand the evolution and divergence of this signalling pathway in animals. Using genome and transcriptome mining we confirmed the presence of nodal and Pitx in a range of additional animal taxa for which their presence has not yet been described. In situ hybridization was used to show the embryonic expression of these genes in brachiopods and planarians. We show that both nodal and Pitx genes are broadly conserved across the Spiralia, and nodal likely appeared in the Bilaterian stem lineage after the divergence of the Acoelomorpha. Furthermore, both nodal and Pitx mRNA appears to be expressed in an asymmetric fashion in the brachiopod Terebratalia transversa. No evidence for the presence of a Lefty ortholog could be found in the non-deuterostome genomic resources examined. Nodal expression is asymmetric in a number of spiralian lineages, indicating a possible ancestral role of the Nodal/Pitx cascade in the establishment of asymmetries across the Bilateria.

Draft genome assemblies and predicted microRNA complements of the intertidal lophotrochozoans Patella vulgata (Mollusca, Patellogastropoda) and Spirobranchus (Pomatoceros) lamarcki (Annelida, Serpulida)

MicroRNAs (miRNA) are small non-coding RNAs that act post-transcriptionally to regulate gene expression levels. Some studies have indicated that microRNAs may have low homoplasy, and as a consequence the phylogenetic distribution of microRNA families has been used to study animal evolutionary relationships. Limited levels of lineage sampling, however, may distort such analyses. Lophotrochozoa is an under-sampled taxon that includes molluscs, annelids and nemerteans, among other phyla. Here, we present two novel draft genomes, those of the limpet Patella vulgata and polychaete Spirobranchus (Pomatoceros) lamarcki. Surveying these genomes for known microRNAs identifies numerous potential orthologues, including a number that have been considered to be confined to other lineages. RT-PCR demonstrates that some of these (miR-1285, miR-1287, miR-1957, miR-1983 and miR-3533), previously thought to be found only in vertebrates, are expressed. This study provides genomic resources for two lophotrochozoans and reveals patterns of microRNA evolution that could be hidden by more restricted sampling.

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.

Identification of putative ecdysteroid and juvenile hormone pathway genes in the shrimp Neocaridina denticulata

Although the sesquiterpenoid juvenile hormone (JH) and the steroidal ecdysteroids are of vital importance to the development and reproduction of insects, our understanding of the evolution of these crucial hormonal regulators in other arthropods is limited. To better understand arthropod hormone evolution and regulation, here we describe the hormonal pathway genes (e.g. those involved in hormone biosynthesis, degradation, regulation and signal transduction) of a new decapod model, the shrimp Neocaridina denticulata. The majority of known insect sesquiterpenoid and ecdysteroid pathway genes and their regulators are contained in the N. denticulata genome. In the sesquiterpenoid pathway, these include biosynthetic pathway components: juvenile hormone acid methyltransferase (JHAMT); hormone binding protein: juvenile hormone binding protein (JHBP); and degradation pathway components: juvenile hormone esterase (JHE), juvenile hormone esterase binding protein (JHEBP) and juvenile hormone epoxide hydrolase (JHEH), with the JHBP, JHEBP and JHEH genes being discovered in a crustacean for the first time here. Ecdysteroid biosynthetic pathway genes identified include spook, phantom, disembodied, shadow and CYP18. Potential hormonal regulators and signal transducers such as allatostatins (ASTs), Methoprene-tolerant (Met), Retinoid X receptor (RXR), Ecdysone receptor (EcR), calponin-like protein Chd64, FK509-binding protein (FKBP39), Broad-complex (Br-c), and crustacean hyperglycemic hormone/molt-inhibiting hormone/gonad-inhibiting hormone (CHH/MIH/GIH) genes are all present in the shrimp N. denticulata. To our knowledge, this is the first report of these hormonal pathways and their regulatory genes together in a single decapod, providing a vital resource for further research into development, reproduction, endocrinology and evolution of crustaceans, and arthropods in general.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

Components of the dorsal-ventral pathway also contribute to anterior-posterior patterning in honeybee embryos (Apis mellifera)

Background

A key early step in embryogenesis is the establishment of the major body axes; the dorsal-ventral (DV) and anterior-posterior (AP) axes. Determination of these axes in some insects requires the function of different sets of signalling pathways for each axis. Patterning across the DV axis requires interaction between the Toll and Dpp/TGF-β pathways, whereas patterning across the AP axis requires gradients of bicoid/orthodenticle proteins and the actions of a hierarchy of gene transcription factors. We examined the expression and function of Toll and Dpp signalling during honeybee embryogenesis to assess to the role of these genes in DV patterning.

Results

Pathway components that are required for dorsal specification in Drosophila are expressed in an AP-restricted pattern in the honeybee embryo, including Dpp and its receptor Tkv. Components of the Toll pathway are expressed in a more conserved pattern along the ventral axis of the embryo. Late-stage embryos from RNA interference (RNAi) knockdown of Toll and Dpp pathways had both DV and AP patterning defects, confirmed by staining with Am-sna, Am-zen, Am-eve, and Am-twi at earlier stages. We also identified two orthologues of dorsal in the honeybee genome, with one being expressed during embryogenesis and having a minor role in axis patterning, as determined by RNAi and the other expressed during oogenesis.

Conclusions

We found that early acting pathways (Toll and Dpp) are involved not only in DV patterning but also AP patterning in honeybee embryogenesis. Changes to the expression patterns and function of these genes may reflect evolutionary changes in the placement of the extra-embryonic membranes during embryogenesis with respect to the AP and DV axes.

How are comparative genomics and the study of microRNAs changing our views on arthropod endocrinology and adaptations to the environment?

As the last few decades of work has shown, precise regulation of biosynthesis and release of arthropod hormones is essential to cope with environmental stresses and challenges. In crustaceans and insects, the sesquiterpenoids methyl farnesoate (MF), farnesoic acid (FA) and juvenile hormone (JH) regulate many developmental, physiological, and reproductive processes. In this review, we discuss how comparative genomics has and will impact our views on arthropod endocrinology. We will also highlight the current knowledge of regulation of genes involved in arthropod hormone biosynthesis by microRNAs, and describe the potential insights into arthropod endocrinology, evolution, and adaptation that are likely to come from the study of microRNAs.