Molluscan Genomes Research Articles

Two draft genomes of enigmatic Solenogastres (Mollusca, Aplacophora) Epimenia babai and Neomenia megatrapezata

Abstract* Many molluscan genomes have been published to date, however only three are from representatives of the subphylum Aculifera (Polyplacophora, Caudofoveata, and Solenogastres), the sister taxon to all other molluscs. Currently, genomic resources are completely lacking for Solenogastres. This gap in knowledge hinders comparative and evolutionary studies. Here, we sequenced the genomes of the solenogaster aplacophorans Epimenia babai Salvini-Plawen, 1997 and Neomenia megatrapezata Salvini-Plawen & Paar-Gausch, 2004 using a hybrid approach combining Oxford Nanopore and Illumina reads. For E. babai, we produced a 628 Mbp haploid assembly (N50 = 413 Kbp, L50 = 370) that is rather complete with a BUSCO completeness score of 90.1% (82.0% single, 8.1% duplicated, 6.0% fragmented, and 3.9% missing). For N. megatrapezata, we produced a 412 Mbp haploid assembly (N50 = 132 Kbp, L50 = 881) that is also rather complete with a BUSCO completeness score of 85.1% (81.7% single, 3.4% duplicated, 8.1% fragmented, and 6.8% missing). Our annotation pipeline predicted 25,393 gene models for E. babai with a BUSCO score of 92.4% (80.5% single, 11.9% duplicated, 4.9% fragmented, and 2.7% missing) and 22,463 gene models for N. megatrapezata with a BUSCO score of 90.2% (81.0% single, 9.2% duplicated, 4.7% fragmented, and 5.1% missing). Phylogenomic analysis recovered Solenogastres as the sister taxon to Polyplacophora and Aculifera as the sister taxon to all other sampled molluscs with maximal support. These represent the first whole-genome resources for Solenogastres and will be valuable for future studies investigating this understudied group and molluscan evolution as a whole.

Molluscan Genomes Reveal Extensive Differences in Photopigment Evolution Across the Phylum.

In animals, opsins and cryptochromes are major protein families that transduce light signals when bound to light-absorbing chromophores. Opsins are involved in various light-dependent processes, like vision, and have been co-opted for light-independent sensory modalities. Cryptochromes are important photoreceptors in animals, generally regulating circadian rhythm, they belong to a larger protein family with photolyases, which repair UV-induced DNA damage. Mollusks are great animals to explore questions about light sensing as eyes have evolved multiple times across, and within, taxonomic classes. We used molluscan genome assemblies from 80 species to predict protein sequences and examine gene family evolution using phylogenetic approaches. We found extensive opsin family expansion and contraction, particularly in bivalve xenopsins and gastropod Go-opsins, while other opsins, like retinochrome, rarely duplicate. Bivalve and gastropod lineages exhibit fluctuations in opsin repertoire, with cephalopods having the fewest number of opsins and loss of at least 2 major opsin types. Interestingly, opsin expansions are not limited to eyed species, and the highest opsin content was seen in eyeless bivalves. The dynamic nature of opsin evolution is quite contrary to the general lack of diversification in mollusk cryptochromes, though some taxa, including cephalopods and terrestrial gastropods, have reduced repertoires of both protein families. We also found complete loss of opsins and cryptochromes in multiple, but not all, deep-sea species. These results help set the stage for connecting genomic changes, including opsin family expansion and contraction, with differences in environmental, and biological features across Mollusca.

De novo genome assembly and genome skims reveal LTRs dominate the genome of a limestone endemic Mountainsnail (Oreohelix idahoensis)

BackgroundCalcareous outcrops, rocky areas composed of calcium carbonate (CaCO3), often host a diverse, specialized, and threatened biomineralizing fauna. Despite the repeated evolution of physiological and morphological adaptations to colonize these mineral rich substrates, there is a lack of genomic resources for calcareous rock endemic species. This has hampered our ability to understand the genomic mechanisms underlying calcareous rock specialization and manage these threatened species.ResultsHere, we present a new draft genome assembly of the threatened limestone endemic land snail Oreohelix idahoensis and genome skim data for two other Oreohelix species. The O. idahoensis genome assembly (scaffold N50: 404.19 kb; 86.6% BUSCO genes) is the largest (~ 5.4 Gb) and most repetitive mollusc genome assembled to date (85.74% assembly size). The repetitive landscape was unusually dominated by an expansion of long terminal repeat (LTR) transposable elements (57.73% assembly size) which have shaped the evolution genome size, gene composition through retrotransposition of host genes, and ectopic recombination. Genome skims revealed repeat content is more than 2–3 fold higher in limestone endemic O. idahoensis compared to non-calcareous Oreohelix species. Gene family size analysis revealed stress and biomineralization genes have expanded significantly in the O. idahoensis genome.ConclusionsHundreds of threatened land snail species are endemic to calcareous rock regions but there are very few genomic resources available to guide their conservation or determine the genomic architecture underlying CaCO3 resource specialization. Our study provides one of the first high quality draft genomes of a calcareous rock endemic land snail which will serve as a foundation for the conservation genomics of this threatened species and for other groups. The high proportion and activity of LTRs in the O. idahoensis genome is unprecedented in molluscan genomics and sheds new light how transposable element content can vary across molluscs. The genomic resources reported here will enable further studies of the genomic mechanisms underlying calcareous rock specialization and the evolution of transposable element content across molluscs.

Evolutionary history of the NLR gene families across lophotrochozoans

NOD-like receptor (NLR) genes are critical innate immune receptors in animals and plants. Lophotrochozoans represent one of the most species-rich superphyla that includes molluscs, segmented worms, flatworms, bryozoans, and other invertebrates, which is crucial to our understanding of immune system evolution in bilaterians. However, NLRs have not been systematically described in lophotrochozoans. We annotated 185 NLRs in 29 lophotrochozoan genomes, and analyzed their domain organization, phylogenetic distribution, molecular evolution, and gene expression. We found that all the 24 molluscan genomes studied encoded no more than three NLRs. None of these molluscan NLRs represented an inducible expression pattern under the infection of eight pathogens; some molluscan NLRs showed developmental stage-specific expression patterns. Instead, 29 molluscan incomplete NLR (incNLR) genes, encoding for proteins absent in the NACHT domain were upregulated under pathogen infection. We also documented the species-specific expansion of NLRs in the clades Polychaeta and Pteriidae. Our study revealed that gene duplication, domain shuffling, gene loss, and novel expression pattern played important roles in the molecular evolution of NLRs.

Genome of the lepidopleurid chiton Hanleya hanleyi (Mollusca, Polyplacophora).

Mollusca is the second most species-rich phylum and includes animals as disparate as octopuses, clams, and chitons. Dozens of molluscan genomes are available, but only one representative of the subphylum Aculifera, the sister taxon to all other molluscs, has been sequenced to date, hindering comparative and evolutionary studies. To facilitate evolutionary studies across Mollusca, we sequenced the genome of a second aculiferan mollusc, the lepidopleurid chiton Hanleya hanleyi (Bean 1844), using a hybrid approach combining Oxford Nanopore and Illumina reads. After purging redundant haplotigs and removing contamination from this 1.3% heterozygous genome, we produced a 2.5 Gbp haploid assembly (>4X the size of the other chiton genome sequenced to date) with an N50 of 65.0 Kbp. Despite a fragmented assembly, the genome is rather complete (92.0% of BUSCOs detected; 79.4% complete plus 12.6% fragmented). Remarkably, the genome has the highest repeat content of any molluscan genome reported to date (>66%). Our gene annotation pipeline predicted 69,284 gene models (92.9% of BUSCOs detected; 81.8% complete plus 11.1% fragmented) of which 35,362 were supported by transcriptome and/or protein evidence. Phylogenomic analysis recovered Polyplacophora sister to all other sampled molluscs with maximal support. The Hanleya genome will be a valuable resource for studies of molluscan biology with diverse potential applications ranging from evolutionary and comparative genomics to molecular ecology.

Neuronal Transcriptome Analysis of a Widely Recognised Molluscan Model Organism Highlights the Absence of Key Proteins Involved in the De Novo Synthesis and Receptor-Mediation of Sex Steroids in Vertebrates

Over the last ten years, the interpretation of the presence of vertebrate sex steroids in molluscs has changed dramatically. Evidence has been accumulating that CYP11A and CYP19A genes (encoding cholesterol side-chain cleavage enzyme and aromatase), that are crucial for the biosynthesis of sex steroids in vertebrates, as well as key functional sex steroid receptors, are missing in molluscan genomes. To provide further evidence, we sequenced the whole transcriptome of the central nervous system of the great pond snail (Lymnaea stagnalis) and screened it for sequences homologous to those used in the generally accepted vertebrate sex steroidogenesis pathway as well as the known sex steroid receptor-related genes (such as CYP11A, CYP19A, 3β-HSD, nPR, and nAR). Our screening confirmed the absence of several key sequences that are essential to accomplish a full sex steroid biosynthesis pathway similar to that of vertebrates. There was also no evidence for nuclear sex steroid receptors. Our findings support the contention that molluscan endocrinology differs from the well-characterized vertebrate endocrine system.

The 'Shellome' of the Crocus Clam Tridacna crocea Emphasizes Essential Components of Mollusk Shell Biomineralization.

Molluscan shells are among the most fascinating research objects because of their diverse morphologies and textures. The formation of these delicate biomineralized structures is a matrix-mediated process. A question that arises is what are the essential components required to build these exoskeletons. In order to understand the molecular mechanisms of molluscan shell formation, it is crucial to identify organic macromolecules in different shells from diverse taxa. In the case of bivalves, however, taxon sampling in previous shell proteomics studies are focused predominantly on representatives of the class Pteriomorphia such as pearl oysters, edible oysters and mussels. In this study, we have characterized the shell organic matrix from the crocus clam, Tridacna crocea, (Heterodonta) using various biochemical techniques, including SDS-PAGE, FT-IR, monosaccharide analysis, and enzyme-linked lectin assay (ELLA). Furthermore, we have identified a number of shell matrix proteins (SMPs) using a comprehensive proteomics approach combined to RNA-seq. The biochemical studies confirmed the presence of proteins, polysaccharides, and sulfates in the T. crocea shell organic matrix. Proteomics analysis revealed that the majority of the T. crocea SMPs are novel and dissimilar to known SMPs identified from the other bivalve species. Meanwhile, the SMP repertoire of the crocus clam also includes proteins with conserved functional domains such as chitin-binding domain, VWA domain, and protease inhibitor domain. We also identified BMSP (Blue Mussel Shell Protein, originally reported from Mytilus), which is widely distributed among molluscan shell matrix proteins. Tridacna SMPs also include low-complexity regions (LCRs) that are absent in the other molluscan genomes, indicating that these genes may have evolved in specific lineage. These results highlight the diversity of the organic molecules – in particular proteins – that are essential for molluscan shell formation.

Benchmarking Oxford Nanopore read assemblers for high-quality molluscan genomes.

Choosing the optimum assembly approach is essential to achieving a high-quality genome assembly suitable for comparative and evolutionary genomic investigations. Significant recent progress in long-read sequencing technologies such as PacBio and Oxford Nanopore Technologies (ONT) has also brought about a large variety of assemblers. Although these have been extensively tested on model species such as Homo sapiens and Drosophila melanogaster, such benchmarking has not been done in Mollusca, which lacks widely adopted model species. Molluscan genomes are notoriously rich in repeats and are often highly heterozygous, making their assembly challenging. Here, we benchmarked 10 assemblers based on ONT raw reads from two published molluscan genomes of differing properties, the gastropod Chrysomallon squamiferum (356.6 Mb, 1.59% heterozygosity) and the bivalve Mytilus coruscus (1593 Mb, 1.94% heterozygosity). By optimizing the assembly pipeline, we greatly improved both genomes from previously published versions. Our results suggested that 40-50X of ONT reads are sufficient for high-quality genomes, with Flye being the recommended assembler for compact and less heterozygous genomes exemplified by C. squamiferum, while NextDenovo excelled for more repetitive and heterozygous molluscan genomes exemplified by M. coruscus. A phylogenomic analysis using the two updated genomes with 32 other published high-quality lophotrochozoan genomes resulted in maximum support across all nodes, and we show that improved genome quality also leads to more complete matrices for phylogenomic inferences. Our benchmarking will ensure efficiency in future assemblies for molluscs and perhaps also for other marine phyla with few genomes available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Potential of genomic technologies to improve disease resistance in molluscan aquaculture.

Molluscan aquaculture is a major contributor to global seafood production, but is hampered by infectious disease outbreaks that can cause serious economic losses. Selective breeding has been widely used to improve disease resistance in major agricultural and aquaculture species, and has clear potential in molluscs, albeit its commercial application remains at a formative stage. Advances in genomic technologies, especially the development of cost-efficient genomic selection, have the potential to accelerate genetic improvement. However, tailored approaches are required owing to the distinctive reproductive and life cycle characteristics of molluscan species. Transgenesis and genome editing, in particular CRISPR/Cas systems, have been successfully trialled in molluscs and may further understanding and improvement of genetic resistance to disease through targeted changes to the host genome. Whole-organism genome editing is achievable on a much greater scale compared to other farmed species, making genome-wide CRISPR screening approaches plausible. This review discusses the current state and future potential of selective breeding, genomic tools and genome editing approaches to understand and improve host resistance to infectious disease in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'.

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Extensive genome-wide duplications in the eastern oyster (Crassostrea virginica).

Genomic structural variation is an important source of genetic and phenotypic diversity, playing a critical role in evolution. The recent availability of a high-quality reference genome for the eastern oyster, Crassostrea virginica, and whole-genome sequence data of samples from across the species range in the USA, provides an opportunity to explore structural variation across the genome of this species. Our analysis shows significantly greater individual-level duplications of regions across the genome than that of most model vertebrate species. Duplications are widespread across all ten chromosomes with variation in frequency per chromosome. The eastern oyster shows a large interindividual variation in duplications as well as particular chromosomal regions with a higher density of duplications. A high percentage of duplications seen in C. virginica lie completely within genes and exons, suggesting the potential for impacts on gene function. These results support the hypothesis that structural changes may play a significant role in standing genetic variation in C. virginica, and potentially have a role in their adaptive and evolutionary success. Altogether, these results suggest that copy number variation plays an important role in the genomic variation of C. virginica. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

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'.

Structure, function and parallel evolution of the bivalve byssus, with insights from proteomes and the zebra mussel genome

The byssus is a structure unique to bivalves. Byssal threads composed of many proteins extend like tendons from muscle cells, ending in adhesive pads that attach underwater. Crucial to settlement and metamorphosis, larvae of virtually all species are byssate. By contrast, in adults, the byssus is scattered throughout bivalves, where it has had profound effects on morphological evolution and been key to adaptive radiations of epifaunal species. I compare byssus structure and proteins in blue mussels ( Mytilus ), by far the best characterized, to zebra mussels ( Dreissena polymorpha ), in which several byssal proteins have been isolated and sequenced. By mapping the adult byssus onto a recent phylogenomic tree, I confirm its independent evolution in these and other lineages, likely parallelisms with common origins in development. While the byssus is superficially similar in Dreissena and Mytilus , in finer detail it is not, and byssal proteins are dramatically different. I used the chromosome-scale D. polymorpha genome we recently assembled to search for byssal genes and found 37 byssal loci on 10 of the 16 chromosomes. Most byssal genes are in small families, with several amino acid substitutions between paralogs. Byssal proteins of zebra mussels and related quagga mussels ( D. rostriformis ) are divergent, suggesting rapid evolution typical of proteins with repetitive low complexity domains. Opportunities abound for proteomic and genomic work to further our understanding of this textbook example of a marine natural material. A priority should be invasive bivalves, given the role of byssal attachment in the spread of, and ecological and economic damage caused by zebra mussels, quagga mussels and others. This article is part of the Theo Murphy meeting issue ‘Molluscan genomics: broad insights and future directions for a neglected phylum’.

Sticky problems: extraction of nucleic acids from molluscs

Traditional molecular methods and omics-techniques across molluscan taxonomy increasingly inform biology of Mollusca. Recovery of DNA and RNA for such studies is challenged by common biological properties of the highly diverse molluscs. Molluscan biomineralization, adhesive structures and mucus involve polyphenolic proteins and mucopolysaccharides that hinder DNA extraction or copurify to inhibit enzyme-catalysed molecular procedures. DNA extraction methods that employ the detergent hexadecyltrimethylammoniumbromide (CTAB) to remove these contaminants importantly facilitate molecular-level study of molluscs. Molluscan pigments may stain DNA samples and interfere with spectrophotometry, necessitating gel electrophoresis or fluorometry for accurate quantification. RNA can reliably be extracted but the ‘hidden break’ in 28S rRNA of molluscs (like most protostomes) causes 18S and 28S rRNA fragments to co-migrate electrophoretically. This challenges the standard quality control based on the ratio of 18S and 28S rRNA, developed for deuterostome animals. High-AT content in molluscan rRNA prevents the effective purification of polyadenylated mRNA. Awareness of these matters aids the continuous expansion of molecular malacology, enabling work also with museum specimens and next-generation sequencing, with the latter imposing unprecedented demands on DNA quality. Alternative methods to extract nucleic acids from molluscs are available from literature and, importantly, from communications with others who study the molecular biology of molluscs. This article is part of the Theo Murphy meeting issue ‘Molluscan genomics: broad insights and future directions for a neglected phylum’.

Mobilizing molluscan models and genomes in biology.

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

MolluscDB: a genome and transcriptome database for molluscs.

As sequencing becomes more accessible and affordable, the analysis of genomic and transcriptomic data has become a cornerstone of many research initiatives. Communities with a focus on particular taxa or ecosystems need solutions capable of aggregating genomic resources and serving them in a standardized and analysis-friendly manner. Taxon-focussed resources can be more flexible in addressing the needs of a research community than can universal or general databases. Here, we present MolluscDB, a genome and transcriptome database for molluscs. MolluscDB offers a rich ecosystem of tools, including an Ensembl browser, a BLAST server for homology searches and an HTTP server from which any dataset present in the database can be downloaded. To demonstrate the utility of the database and verify the quality of its data, we imported data from assembled genomes and transcriptomes of 22 species, estimated the phylogeny of Mollusca using single-copy orthologues, explored patterns of gene family size change and interrogated the data for biomineralization-associated enzymes and shell matrix proteins. MolluscDB provides an easy-to-use and openly accessible data resource for the research community.This article is part of the Theo Murphy meeting issue ‘Molluscan genomics: broad insights and future directions for a neglected phylum’.

Pearls of wisdom-a Theo Murphy issue on molluscan genomics.

Pearls of wisdom-a Theo Murphy issue on molluscan genomics.

Whole-genome re-sequencing data to infer historical demography and speciation processes in land snails: the study of two Candidula sister species.

Despite the global biodiversity of terrestrial gastropods and their ecological and economic importance, the genomic basis of ecological adaptation and speciation in land snail taxa is still largely unknown. Here, we combined whole-genome re-sequencing with population genomics to evaluate the historical demography and the speciation process of two closely related species of land snails from western Europe, Candidula unifasciata and C. rugosiuscula. Historical demographic analysis indicated fluctuations in the size of ancestral populations, probably driven by Pleistocene climatic fluctuations. Although the current population distributions of both species do not overlap, our approximate Bayesian computation model selection approach on several speciation scenarios suggested that gene flow has occurred throughout the divergence process until recently. Positively selected genes diverging early in the process were associated with intragenomic and cyto-nuclear incompatibilities, respectively, potentially fostering reproductive isolation as well as ecological divergence. Our results suggested that the speciation between species entails complex processes involving both gene flow and ecological speciation, and that further research based on whole-genome data can provide valuable understanding on species divergence.This article is part of the Theo Murphy meeting issue ‘Molluscan genomics: broad insights and future directions for a neglected phylum’.

Molluscan mitochondrial genomes break the rules.

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Molluscan phylogenomics requires strategically selected genomes.

The extraordinary diversity in molluscan body plans, and the genomic mechanisms that enable it, remains one of the great questions of evolution. The eight distinct living taxonomic classes of molluscs are each unambiguously monophyletic; however, significant controversy remains about the phylogenetic relationships among those eight branches. Molluscs are the second-largest animal phylum, with over 100 000 living species with broad biological, economic and medical interest. To date, only around 53 genome assemblies have been accessioned to NCBI GenBank covering only four of the eight living molluscan classes. Furthermore, the molluscan taxa where partial or whole-genome assemblies are available are often aberrantly fast evolving or recently derived lineages. Characteristic adaptations provide interesting targets for whole-genome projects, in animals like the scaly-foot snail or octopus, but without basal-branching lineages for comparison, the context of recently derived features cannot be assessed. The currently available genomes also create a non-optimal set of taxa for resolving deeper phylogenetic branches: they are a small sample representing a large group, and those that are available come primarily from a rarefied pool. Thoughtful selection of taxa for future projects should focus on the blank areas of the molluscan tree, which are ripe with opportunities to delve into peculiarities of genome evolution, and reveal the biology and evolutionary history of molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.