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.

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.

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.

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.

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.

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.

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.

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.

Right across the tree of life: the evolution of left-right asymmetry in the Bilateria.

Directional left/right (LR) asymmetries, in which there are consistent, heritable differences in morphology between the left and right sides of bilaterally symmetrical organisms, are found in animals across the Bilateria. For many years, we have lacked evidence for shared mechanisms underlying their development. This led to the supposition that the mechanisms driving establishment of LR asymmetries, and consequently the asymmetries themselves, had evolved separately in the three major Superphyla that constitute the Bilateria. The recent discovery that the transforming growth factor-beta (TGF-B) ligand Nodal plays a role in the regulation of LR asymmetry in both Deuterostomia and Lophotrochozoa has reignited debate in this field, as it suggests that at least this aspect of the development of the LR axis is conserved. In this review, we discuss evidence for shared mechanisms of LR asymmetry establishment across the bilaterian tree of life and consider how these mechanisms might have diverged across the Metazoa over the last 500 million years or so of evolution. As well as the likelihood that Nodal is an ancestral mechanism for regulating LR asymmetry, we reemphasize cytoskeletal architecture as a potential shared mechanism underlying symmetry breaking. However, convergent evolution remains a distinct possibility and study of a wider diversity of species will be needed to distinguish between conserved and lineage-specific mechanisms.

Genomic sequence and experimental tractability of a new decapod shrimp model, Neocaridina denticulata.

The speciose Crustacea is the largest subphylum of arthropods on the planet after the Insecta. To date, however, the only publically available sequenced crustacean genome is that of the water flea, Daphnia pulex, a member of the Branchiopoda. While Daphnia is a well-established ecotoxicological model, previous study showed that one-third of genes contained in its genome are lineage-specific and could not be identified in any other metazoan genomes. To better understand the genomic evolution of crustaceans and arthropods, we have sequenced the genome of a novel shrimp model, Neocaridina denticulata, and tested its experimental malleability. A library of 170-bp nominal fragment size was constructed from DNA of a starved single adult and sequenced using the Illumina HiSeq2000 platform. Core eukaryotic genes, the mitochondrial genome, developmental patterning genes (such as Hox) and microRNA processing pathway genes are all present in this animal, suggesting it has not undergone massive genomic loss. Comparison with the published genome of Daphnia pulex has allowed us to reveal 3750 genes that are indeed specific to the lineage containing malacostracans and branchiopods, rather than Daphnia-specific (E-value: 10⁻6). We also show the experimental tractability of N. denticulata, which, together with the genomic resources presented here, make it an ideal model for a wide range of further aquacultural, developmental, ecotoxicological, food safety, genetic, hormonal, physiological and reproductive research, allowing better understanding of the evolution of crustaceans and other arthropods.

Effects of the toxic dinoflagellate, Alexandrium pacificum, on the marine diatom, Chaetoceros muelleri, and mussel (Perna canaliculus) sperm and hemocytes.

Harmful algal blooms (HABs) of Alexandrium pacificum have affected the Marlborough Sounds in New Zealand since 2010, posing a threat to green-lipped mussel (GLM, Perna canaliculus) farming. Previous studies have shown A. pacificum has negative effects GLM embryos and larvae. To further investigate these toxic mechanisms, in vitro bioassays were conducted on GLM spermatozoa, hemocytes, and the diatom, Chaetoceros muelleri. The three cell types were exposed to several treatments of A. pacificum for 2 h and responses were measured using flow cytometry and pulse amplitude-modulated fluorometry. Significant spermatozoa mortality was recorded in treatments containing A. pacificum cells or fragments, while hemocyte and C. muelleri mortality was recorded in cell-free treatments of A. pacificum which contained paralytic shellfish toxins (PSTs). Variation in sensitivity between cell types as well as the sublethal effects observed, emphasise the diverse toxic mechanisms of A. pacificum on co-occurring species in the environment.

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.

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.

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.

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.

Additive multiple k-mer transcriptome of the keelworm Pomatoceros lamarckii (Annelida; Serpulidae) reveals annelid trochophore transcription factor cassette.

Recent advances in both next-generation sequencing and assembly programmes have made the low-cost construction of transcriptome datasets for non-model species feasible, capable of yielding a raft of information even from less well-transcribed genes. Here we present the results of assemblies performed on a 51-bp paired end Illumina dataset derived from a mixed larval sample of the annelid Pomatoceros lamarckii at 24, 48 and 72 h post-fertilization. We used Oases to assemble 36.5 million paired end reads with k-mer sizes from 21 to 29, followed by amalgamation of assemblies, redundancy removal with Vmatch and TGICL and removal of contigs less than 500 bp in length. This resulted in a final assembly of 50,151 contigs, with a mean length of 1,221 bp and covering 61.3 Mbp. A total of 34,846 (69.4 %) of these returned a BlastX hit above a cutoff of 1.0e (-3), and 17,967 (35.8 %) were assigned at least one GO annotation using Blast2GO. We used the assembly to identify genes belonging to the homeobox superclass and the Fox, Sox and Tbx classes, recovering 37, 16, four and three genes, respectively. This included orthologues of genes previously unidentified in lophotrochozoans and protostomes. Our study illustrates the utility of such transcriptomic assembly methods as a gene discovery tool and greatly expands our knowledge of transcription factor genes in annelids in general and in this species in particular.

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.

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.