A phylogenetic backbone for Bivalvia: an RNA-seq approach

Signatures of Extreme Longevity: A Perspective from Bivalve Molecular Evolution

Among Metazoa, bivalves have the highest lifespan disparity, ranging from 1 to 500+ years, making them an exceptional testing ground to understand mechanisms underlying aging and the evolution of extended longevity. Nevertheless, comparative molecular evolution has been an overlooked approach in this instance. Here, we leveraged transcriptomic resources spanning 30 bivalve species to unravel the signatures of convergent molecular evolution in four long-lived species: Margaritifera margaritifera, Elliptio complanata, Lampsilis siliquoidea, and Arctica islandica (the latter represents the longest-lived noncolonial metazoan known so far). We applied a comprehensive approach-which included inference of convergent dN/dS, convergent positive selection, and convergent amino acid substitution-with a strong focus on the reduction of false positives. Genes with convergent evolution in long-lived bivalves show more physical and functional interactions to each other than expected, suggesting that they are biologically connected; this interaction network is enriched in genes for which a role in longevity has been experimentally supported in other species. This suggests that genes in the network are involved in extended longevity in bivalves and, consequently, that the mechanisms underlying extended longevity are-at least partially-shared across Metazoa. Although we believe that an integration of different genes and pathways is required for the extended longevity phenotype, we highlight the potential central roles of genes involved in cell proliferation control, translational machinery, and response to hypoxia, in lifespan extension.

Comparative Neuroanatomy of Pediveliger Larvae of Various Bivalves from the Sea of Japan

Here, we describe the nervous system structures from pediveligers of eight bivalve species (Callista brevisiphonata, Mactromeris polynyma, Crenomytilus grayanus, Kellia japonica, Mizuhopecten yessoensis, and Azumapecten farreri) with different modes of life in their adult stages, corresponding to the ecological niches that they occupy (burrowing, cemented, byssally attached, and mobile forms). We have identified neuromorphological features of the central and peripheral nervous systems in larval bivalves. We show that the unpaired sensory apical organ is still present in pediveligers along with the developing paired cerebral ganglia characteristic of an adult mollusk. Pediveligers have the pleural ganglia connected to the pedal ganglia via the pedal nerve cords and to the visceral ganglia via the lateral nerve cords. We have found a number of structures of the peripheral nervous system whose presence varies between pediveligers of different species. Mactromeris, Callista, and Pododesmus have 5-HT-immunopositive stomatogastric neurons, whereas the Yesso and Farrer's scallops have an FMRFamide-immunopositive enteric nervous system. The innervation of the anterior part of the velum is connected to a system of the apical organ and cerebral ganglia, and the innervation of the posterior part is connected to the visceral ganglia. Most differences in the structure of the peripheral elements of the nervous system are species-specific and weakly depend on the ecological niche that pediveligers occupy.

A transcriptome-based phylogeny for Polynoidae (Annelida: Aphroditiformia)

Polynoidae is the most diverse radiation of Aphroditiformia and one of the most successful groups of all Annelida in terms of diversity and habitats colonized. With such an unmatched diversity, phylogenetic investigations have struggled to understand their evolutionary relationships. Previous phylogenetic analyses have slowly increased taxon sampling and employed methodologies, but despite their diversity and biological importance, large genomic sampling is limited. To investigate the internal relationships within Polynoidae, we conducted the first phylogenomic analyses of the group based on 12 transcriptomes collected from species inhabiting a broad array of habitats, including shallow and deep waters, as well as hydrothermal vents, anchialine caves and the midwater. Our phylogenomic analyses of Polynoidae recovered congruent tree topologies representing the clades Polynoinae, Macellicephalinae and Lepidonotopodinae. Members of Polynoinae and Macellicephalinae clustered in well supported and independent clades. In contrast, Lepidonotopodinae taxa were always recovered nested within Macellicephalinae. Though our sampling only covers a small proportion of the species known for Polynoidae, our results provide a robust phylogenomic framework to build from, emphasizing previously hypothesized relationships between Macellicephalinae and Lepidonotopodinae taxa, while providing new insights on the origin of enigmatic cave and pelagic lineages.

Castro LF, Froufe E. The Crown Pearl V2: an improved genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Contiguous assemblies are fundamental to deciphering the composition of extant genomes. In molluscs, this is considerably challenging owing to the large size of their genomes, heterozygosity, and widespread repetitive content. Consequently, long-read sequencing technologies are fundamental for high contiguity and quality. The first genome assembly of Margaritifera margaritifera (Linnaeus, 1758) (Mollusca: Bivalvia: Unionida), a culturally relevant, widespread, and highly threatened species of freshwater mussels, was recently generated. However, the resulting genome is highly fragmented since the assembly relied on short-read approaches. Here, an improved reference genome assembly was generated using a combination of PacBio CLR long reads and Illumina paired-end short reads. This genome assembly is 2.4 Gb long, organized into 1,700 scaffolds with a contig N50 length of 3.4 Mbp. The ab initio gene prediction resulted in 48,314 protein-coding genes. Our new assembly is a substantial improvement and an essential resource for studying this species' unique biological and evolutionary features, helping promote its conservation.

Molecular characterization of accripin11, a soluble shell protein with an acidic C-terminus, identified in the prismatic layer of the Mediterranean fan mussel Pinna nobilis (Bivalvia, Pteriomorphia)

We have identified a novel shell protein, accripin11, as a major soluble component of the calcitic prisms of the fan mussel Pinna nobilis. Initially retrieved from a cDNA library, its full sequence is confirmed here by transcriptomic and proteomic approaches. The sequence of the mature protein is 103 residues with a theoretical molecular weight of 11 kDa and is moderately acidic (pI 6.74) except for its C-terminus which is highly enriched in aspartic acid. The protein exhibits a peculiar cysteine pattern in its central domain. The full sequence shares similarity with six other uncharacterized molluscan shell proteins from the orders Ostreida, Pteriida and Mytilida, all of which are pteriomorphids and produce a phylogenetically restricted pattern of nacro-prismatic shell microstructures. This suggests that accripin11 is a member of a family of clade-specific shell proteins. A 3D model of accripin11 was predicted with AlphaFold2, indicating that it possesses three short alpha helices and a disordered C-terminus. Recombinant accripin11 was tested in vitro for its ability to influence the crystallization of CaCO₃ , while a polyclonal antibody was able to locate accripin11 to prismatic extracts, particularly in the acetic acid-soluble matrix. The putative functions of accripin11 are further discussed in relation to shell biomineralization.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

Dynamic evolutionary history of spiralian-specific TALE homeobox genes in mollusks

Homeobox genes play essential roles in the early development of many animals. Although the repertoire of most homeobox genes, including three amino acid loop extension (TALE)-type homeobox genes, is conserved in animals, spiralian-TALE (SPILE) genes are a notable exception. In this study, SPILE genes were extracted from the genomic data of 22 mollusk species and classified into four clades (-A/C, -B, -D, and -E) to determine which SPILE genes exhibit dynamic repertoire changes. While SPILE-D and -E duplications were rarely observed, SPILE-B duplication was observed in the bivalve lineage and SPILE-A/C duplication was observed in multiple clades. Conversely, most or all SPILE genes were lost in cephalopods and in some gastropod lineages. SPILE gene expression patterns were also analyzed in multiple mollusk species using publicly available RNA-seq data. The majority of SPILE genes examined, particularly those in the A/C- and B-clades, were specifically expressed during early development, suggesting that most SPILE genes exert specific roles in early development. This comprehensive cataloging and characterization revealed a dynamic evolutionary history, including SPILE-A/C and -B gene duplications and the loss of SPILE genes in several lineages. Furthermore, this study provides a useful resource for studying the molecular mechanism of spiralian early development and the evolution of young and lineage-specific transcription factors.

Reinvestigating the phylogeny of Myriapoda with more extensive taxon sampling and novel genetic perspective

Background

There have been extensive debates on the interrelationships among the four major classes of Myriapoda-Chilopoda, Symphyla, Diplopoda, and Pauropoda. The core controversy is the position of Pauropoda; that is, whether it should be grouped with Symphyla or Diplopoda as a sister group. Two recent phylogenomic studies separately investigated transcriptomic data from 14 and 29 Myriapoda species covering all four groups along with outgroups, and proposed two different topologies of phylogenetic relationships.

Methods

Building on these studies, we extended the taxon sampling by investigating 39 myriapods and integrating the previously available data with three new transcriptomic datasets generated in this study. Our analyses present the phylogenetic relationships among the four major classes of Myriapoda with a more abundant taxon sampling and provide a new perspective to investigate the above-mentioned question, where visual genes' identification were conducted. We compared the appearance pattern of genes, grouping them according to their classes and the visual pathways involved. Positive selection was detected for all identified visual genes between every pair of 39 myriapods, and 14 genes showed positive selection among 27 pairs.

Results

From the results of phylogenomic analyses, we propose that Symphyla is a sister group of Pauropoda. This stance has also received strong support from tree inference and topology tests.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Expanding the Search for Sperm Transmission Elements in the Mitochondrial Genomes of Bivalve Mollusks

Doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA) in bivalve mollusks is one of the most notable departures from the paradigm of strict maternal inheritance of mtDNA among metazoans. Recently, work on the Mediterranean mussel Mytilus galloprovincialis suggested that a nucleotide motif in the control region of this species, known as the sperm transmission element (STE), helps protect male-transmitted mitochondria from destruction during spermatogenesis. Subsequent studies found similar, yet divergent, STE motifs in other marine mussels. Here, we extend the in silico search for mtDNA signatures resembling known STEs. This search is carried out for the large unassigned regions of 157 complete mitochondrial genomes from within the Mytiloida, Veneroida, Unionoida, and Ostreoida bivalve orders. Based on a sliding window approach, we present evidence that there are additional putative STE signatures in the large unassigned regions of several marine clams and freshwater mussels with DUI. We discuss the implications of this finding for interpreting the origin of doubly uniparental inheritance in ancestral bivalve mollusks, as well as potential future in vitro and in silico studies that could further refine our understanding of the early evolution of this unusual system of mtDNA inheritance.

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

Novel gene rearrangement in the mitochondrial genome of Siliqua minima (Bivalvia, Adapedonta) and phylogenetic implications for Imparidentia

Siliqua minima (Gmelin, 1791) is an important economic shellfish species belonging to the family Pharidae. To date, the complete mitochondrial genome of only one species in this family (Sinonovacula constricta) has been sequenced. Research on the Pharidae family is very limited; to improve the evolution of this bivalve family, we sequenced the complete mitochondrial genome of S. minima by next-generation sequencing. The genome is 17,064 bp in length, consisting of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), and two ribosomal RNA genes (rRNA). From the rearrangement analysis of bivalves, we found that the gene sequences of bivalves greatly variable among species, and with closer genetic relationship, the more consistent of the gene arrangement is higher among the species. Moreover, according to the gene arrangement of seven species from Adapedonta, we found that gene rearrangement among families is particularly obvious, while the gene order within families is relatively conservative. The phylogenetic analysis between species of the superorder Imparidentia using 12 conserved PCGs. The S. minima mitogenome was provided and will improve the phylogenetic resolution of Pharidae species.

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

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

Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates

Bivalves represent valuable taxonomic group for aging studies given their wide variation in longevity (from 1–2 to >500 years). It is well known that aging is associated to the maintenance of Reactive Oxygen Species homeostasis and that mitochondria phenotype and genotype dysfunctions accumulation is a hallmark of these processes. Previous studies have shown that mitochondrial DNA mutation rates are linked to lifespan in vertebrate species, but no study has explored this in invertebrates. To this end, we performed a Bayesian Phylogenetic Covariance model of evolution analysis using 12 mitochondrial protein-coding genes of 76 bivalve species. Three life history traits (maximum longevity, generation time and mean temperature tolerance) were tested against 1) synonymous substitution rates (dS), 2) conservative amino acid replacement rates (Kc) and 3) ratios of radical over conservative amino acid replacement rates (Kr/Kc). Our results confirm the already known correlation between longevity and generation time and show, for the first time in an invertebrate class, a significant negative correlation between dS and longevity. This correlation was not as strong when generation time and mean temperature tolerance variations were also considered in our model (marginal correlation), suggesting a confounding effect of these traits on the relationship between longevity and mtDNA substitution rate. By confirming the negative correlation between dS and longevity previously documented in birds and mammals, our results provide support for a general pattern in substitution rates.

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

A polyvalent and universal tool for genomic studies in gastropod molluscs (Heterobranchia)

Molluscs are the second most diverse animal phylum and heterobranch gastropods present ~ 44,000 species. These comprise fascinating creatures with huge morphological and ecological disparity. Such great diversity comes with even larger phylogenetic uncertainty and many taxa have been largely neglected in molecular assessments. Genomic tools have provided resolution to deep cladogenic events but generating large numbers of transcriptomes/genomes is expensive and usually requires fresh material. Here we leverage a target enrichment approach to design and synthesize a probe set based on available genomes and transcriptomes across Heterobranchia. Our probe set contains 57,606 70mer baits and targets a total of 2,259 ultra-conserved elements (UCEs). Post-sequencing capture efficiency was tested against 31 marine heterobranchs from major groups, including Acochlidia, Acteonoidea, Aplysiida, Cephalaspidea, Pleurobranchida, Pteropoda, Runcinida, Sacoglossa, and Umbraculida. The combined Trinity and Velvet assemblies recovered up to 2,211 UCEs in Tectipleura, up to 1,978 in Nudipleura, and up to 1,927 in Acteonoidea, the latter two being the most distantly related taxa to our core study group. Total alignment length was 525,599 bp and contained 52% informative sites and 21% missing data. Maximum-likelihood and Bayesian inference approaches recovered the monophyly of all orders tested as well as the larger clades Nudipleura, Panpulmonata, and Euopisthobranchia. The successful enrichment of diversely preserved material and DNA concentrations demonstrate the polyvalent nature of UCEs, and the universality of the probe set designed. We believe this probe set will enable multiple, interesting lines of research, that will benefit from an inexpensive and largely informative tool that will, additionally, benefit from the access to museum collections to gather genomic data.

Mitogenomic architecture of the multivalent endemic black clam (Villorita cyprinoides) and its phylogenetic implications

The Indian black clam Villorita cyprinoides (Family: Cyrenidae), an extractive commercially exploited species with aquaculture importance contributing more than 70% of clam fishery in India, is endemic to the Indian peninsula. Currently, there is very sparse information, especially on the molecular data of Villorita. The present study aims to provide a comprehensive knowledge of mitogenome architecture and assess the phylogenetic status of Cyrenidae. This has resulted in reporting the first complete mitogenome of V. cyprinoides using next-generation sequencing technology. The A+T circular mitogenome was 15,880 bp long, exhibiting 13 protein-coding genes (PCGs) including ATP8 (absent in several bivalves), 22 transfer RNA, and two ribosomal RNA genes residing in the heavy strand in a clockwise orientation and a gene order akin to Corbicula fluminea. The molecular phylogeny inferred from a concatenated multi-gene sequence [14 mitochondrial (12 PCGs, rrnS and rrnL) and two nuclear genes (Histone H3, 18S rRNA)] from 47 representative species of superorder Imparidentia, clustered V. cyprinoides and Cyrenid clams to a single clade supporting the monophyly of Cyrenidae. The subsequent mitochondrial gene order analysis substantiates the close relationship of V. cyprinoides and C. fluminea, analogous to phylogenetic output. The multilocus tree topology calibrated with verified fossil data deciphered the origin and diversification of Cyrenid clams during late Triassic-early Jurassic. The data derived from this study shall contribute remarkably for further insights on cryptic species identification, molecular characterization of bivalve mitogenomes and mitochondrial evolutionary history of genus Villorita. Moreover, complete mitogenome can aid in potential marker development for assessing the genetic health of black clam populations.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Shedding light: a phylotranscriptomic perspective illuminates the origin of photosymbiosis in marine bivalves

BACKGROUND

Photosymbiotic associations between metazoan hosts and photosynthetic dinoflagellates are crucial to the trophic and structural integrity of many marine ecosystems, including coral reefs. Although extensive efforts have been devoted to study the short-term ecological interactions between coral hosts and their symbionts, long-term evolutionary dynamics of photosymbiosis in many marine animals are not well understood. Within Bivalvia, the second largest class of mollusks, obligate photosymbiosis is found in two marine lineages: the giant clams (subfamily Tridacninae) and the heart cockles (subfamily Fraginae), both in the family Cardiidae. Morphologically, giant clams show relatively conservative shell forms whereas photosymbiotic fragines exhibit a diverse suite of anatomical adaptations including flattened shells, leafy mantle extensions, and lens-like microstructural structures. To date, the phylogenetic relationships between these two subfamilies remain poorly resolved, and it is unclear whether photosymbiosis in cardiids originated once or twice.

RESULTS

In this study, we establish a backbone phylogeny for Cardiidae utilizing RNASeq-based transcriptomic data from Tridacninae, Fraginae and other cardiids. A variety of phylogenomic approaches were used to infer the relationship between the two groups. Our analyses found conflicting gene signals and potential rapid divergence among the lineages. Overall, results support a sister group relationship between Tridacninae and Fraginae, which diverged during the Cretaceous. Although a sister group relationship is recovered, ancestral state reconstruction using maximum likelihood-based methods reveals two independent origins of photosymbiosis, one at the base of Tridacninae and the other within a symbiotic Fraginae clade.

CONCLUSIONS

The newly revealed common ancestry between Tridacninae and Fraginae brings a possibility that certain genetic, metabolic, and/or anatomical exaptations existed in their last common ancestor, which promoted both lineages to independently establish photosymbiosis, possibly in response to the modern expansion of reef habitats.

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.

Molecular characterization of a whirlin-like protein with biomineralization-related functions from the shell of Mytilus coruscus

Mollusc shells are produced from calcified skeletons and have excellent mechanical properties. Shell matrix proteins (SMPs) have important functions in shell formation. A 16.6 kDa whirlin-like protein (WLP) with a PDZ domain was identified in the shell of Mytilus coruscus as a novel SMP. In this study, the expression, function, and location of WLP were analysed. The WLP gene was highly expressed and specifically located in the adductor muscle and mantle. The expression of recombinant WLP (rWLP) was associated with morphological change, polymorphic change, binding ability, and crystallization rate inhibition of the calcium carbonate crystals in vitro. In addition, an anti-rWLP antibody was prepared, and the results from immunohistochemistry and immunofluorescence analyses revealed the specific location of the WLP in the mantle, adductor muscle, and myostracum layer of the shell, suggesting multiple functions for WLP in biomineralization, muscle-shell attachment, and muscle attraction. Furthermore, results from a pull-down analysis revealed 10 protein partners of WLP in the shell matrices and a possible network of interacting WLPs in the shell. In addition, in this study, one of the WLP partners, actin, was confirmed to have the ability to bind WLP. These results expand the understanding of the functions of PDZ-domain-containing proteins in biomineralization and provide clues for determining the mechanisms of myostracum formation and muscle-shell attachment.

The quagga mussel genome and the evolution of freshwater tolerance

Freshwater dreissenid mussels evolved from marine ancestors during the Miocene ∼30 million years ago and today include some of the most successful and destructive invasive species of freshwater environments. Here, we sequenced the genome of the quagga mussel Dreissena rostriformis to identify adaptations involved in embryonic osmoregulation. We provide evidence that a lophotrochozoan-specific aquaporin water channel, a vacuolar ATPase subunit and a sodium/hydrogen exchanger are involved in osmoregulation throughout early cleavage, during which time large intercellular fluid-filled 'cleavage cavities' repeatedly form, coalesce and collapse, expelling excess water to the exterior. Independent expansions of aquaporins coinciding with at least five freshwater colonization events confirm their role in freshwater adaptation. Repeated aquaporin expansions and the evolution of membrane-bound fluid-filled osmoregulatory structures in diverse freshwater taxa point to a fundamental principle guiding the evolution of freshwater tolerance and provide a framework for future species control efforts.

Resolving the relationships of clams and cockles: dense transcriptome sampling drastically improves the bivalve tree of life

Bivalvia has been the subject of extensive recent phylogenetic work to attempt resolving either the backbone of the bivalve tree using transcriptomic data, or the tips using morpho-anatomical data and up to five genetic markers. Yet the first approach lacked decisive taxon sampling and the second failed to resolve many interfamilial relationships, especially within the diverse clade Imparidentia. Here we combine dense taxon sampling with 108 deep-sequenced Illumina-based transcriptomes to provide resolution in nodes that required additional study. We designed specific data matrices to address the poorly resolved relationships within Imparidentia. Our results support the overall backbone of the bivalve tree, the monophyly of Bivalvia and all its main nodes, although the monophyly of Protobranchia remains less clear. Likewise, the inter-relationships of the six main bivalve clades were fully supported. Within Imparidentia, resolution increases when analysing Imparidentia-specific matrices. Lucinidae, Thyasiridae and Gastrochaenida represent three early branches. Gastrochaenida is sister group to all remaining imparidentians, which divide into six orders. Neoheterodontei is always fully supported, and consists of Sphaeriida, Myida and Venerida, with the latter now also containing Mactroidea, Ungulinoidea and Chamidae, a family particularly difficult to place in earlier work. Overall, our study, by using densely sampled transcriptomes, provides the best-resolved bivalve phylogeny to date.

New data from Monoplacophora and a carefully-curated dataset resolve molluscan relationships

Relationships among the major lineages of Mollusca have long been debated. Morphological studies have considered the rarely collected Monoplacophora (Tryblidia) to have several plesiomorphic molluscan traits. The phylogenetic position of this group is contentious as morphologists have generally placed this clade as the sister taxon of the rest of Conchifera whereas earlier molecular studies supported a clade of Monoplacophora + Polyplacophora (Serialia) and phylogenomic studies have generally recovered a clade of Monoplacophora + Cephalopoda. Phylogenomic studies have also strongly supported a clade including Gastropoda, Bivalvia, and Scaphopoda, but relationships among these taxa have been inconsistent. In order to resolve conchiferan relationships and improve understanding of early molluscan evolution, we carefully curated a high-quality data matrix and conducted phylogenomic analyses with broad taxon sampling including newly sequenced genomic data from the monoplacophoran Laevipilina antarctica. Whereas a partitioned maximum likelihood (ML) analysis using site-homogeneous models recovered Monoplacophora sister to Cephalopoda with moderate support, both ML and Bayesian inference (BI) analyses using mixture models recovered Monoplacophora sister to all other conchiferans with strong support. A supertree approach also recovered Monoplacophora as the sister taxon of a clade composed of the rest of Conchifera. Gastropoda was recovered as the sister taxon of Scaphopoda in most analyses, which was strongly supported when mixture models were used. A molecular clock based on our BI topology dates diversification of Mollusca to ~546 MYA (+/- 6 MYA) and Conchifera to ~540 MYA (+/- 9 MYA), generally consistent with previous work employing nuclear housekeeping genes. These results provide important resolution of conchiferan mollusc phylogeny and offer new insights into ancestral character states of major mollusc clades.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Environmental transcriptomes of invasive dreissena, a model species in ecotoxicology and invasion biology

Dreissenids are established model species for ecological and ecotoxicological studies, since they are sessile and filter feeder organisms and reflect in situ freshwater quality. Despite this strong interest for hydrosystem biomonitoring, omics data are still scarce. In the present study, we achieved full de novo assembly transcriptomes of digestive glands to gain insight into Dreissena polymorpha and D. rostriformis bugensis molecular knowledge. Transcriptomes were obtained by Illumina RNA sequencing of seventy-nine organisms issued from fifteen populations inhabiting sites that exhibits multiple freshwater contamination levels and different hydrosystem topographies (open or closed systems). Based on a recent de novo assembly algorithm, we carried out a complete, quality-checked and annotated transcriptomes. The power of the present study lies in the completeness of transcriptomes gathering multipopulational organisms sequencing and its full availability through an open access interface that gives a friendly and ready-to-use access to data. The use of such data for proteogenomic and targeted biological pathway investigations purpose is promising as they are first full transcriptomes for this two Dreissena species.

A transcriptome-based phylogenetic study of hard ticks (Ixodidae)

Hard ticks are widely distributed across temperate regions, show strong variation in host associations, and are potential vectors of a diversity of medically important zoonoses, such as Lyme disease. To address unresolved issues with respect to the evolutionary relationships among certain species or genera, we produced novel RNA-Seq data sets for nine different Ixodes species. We combined this new data with 18 data sets obtained from public databases, both for Ixodes and non-Ixodes hard tick species, using soft ticks as an outgroup. We assembled transcriptomes (for 27 species in total), predicted coding sequences and identified single copy orthologues (SCO). Using Maximum-likelihood and Bayesian frameworks, we reconstructed a hard tick phylogeny for the nuclear genome. We also obtained a mitochondrial DNA-based phylogeny using published genome sequences and mitochondrial sequences derived from the new transcriptomes. Our results confirm previous studies showing that the Ixodes genus is monophyletic and clarify the relationships among Ixodes sub-genera. This work provides a baseline for studying the evolutionary history of ticks: we indeed found an unexpected acceleration of substitutions for mitochondrial sequences of Prostriata, and for nuclear and mitochondrial genes of two species of Rhipicephalus, which we relate with patterns of genome architecture and changes of life-cycle, respectively.

Biodiversity Assessment, DNA Barcoding, and the Minority Majority

The majority of species on Earth are in "under-studied" groups, and indeed probably the majority of species remain undiscovered and undescribed. Species are natural units of evolution, and they are formed from branching phylogenetic processes that have a mathematical structure. So it follows that we should be able to develop a set of general principles that describe global patterns of species groups, like genera. Understanding such patterns would lend considerable power to the approach of "taxonomic surrogacy." In environmental assessments, ecology, and paleontology, it is common to substitute genus-level or family-level identification where definitive species identification is impractical. Clarity and confidence in fundamental patterns, based on a robust null model for species and genus level diversity, can accelerate species discovery: there are more species in the tropics, species-poor genera are very common, large genera are rare. Much hope has been placed in DNA barcoding as an effective tool to increase the pace of species discovery, but it is abundantly clear that certain mitochondrial DNA (mtDNA) markers are more or less variable in different clades and universal threshold values are impractical to delimit species. This study further examines the patterns of divergence in one common mtDNA barcode fragment, cytochrome c oxidase subunit 1at the genus level. We compared pairwise divergence in this fragment between two animal clades that have similar species richness but different evolutionary histories: birds and bivalves. We analyzed quality controlled alignments of over 39,000 published sequences in 1223 genera. Median pairwise differences at the genus level are positively correlated with the species richness of a genus, and this is not dependent of the number of sequences sampled. Unsurprisingly, sequence divergence in vertebrates was far more constrained than in evolutionarily more ancient non-vertebrate clades. Differences among the groups examined highlight the need for DNA barcode approaches to be considered in the context of specific biological groups. Vertebrates are better studied, but not necessarily representative of the majority of biodiversity. A technique that provides powerful insights for vertebrate species may be ineffective for the majority of organisms.

Phylogenomic analyses and divergence time estimation of Elateroidea (Coleoptera) based on RNA-Seq data

Bioluminescence, the emission of visible light in a living organism, is an intriguing phenomenon observed in different species and environments. In terrestrial organisms, the bioluminescence is observed mainly in beetles of the Elateroidea superfamily (Coleoptera). Several phylogenetic studies have been used different strategies to propose a scenario for the origin and evolution of the bioluminescence within this group, however some of them showed incongruences, mainly about the relationship of the bioluminescent families. In order to increase the number of molecular markers available for Elateroidea species and to propose a more accurate phylogeny, with high supported topology, we employed the Next-Generation Sequencing (NGS) methodology to perform the RNA-Seq analysis of luminescent (Elateridae, Phengodidae, Rhagophthalmidae, and Lampyridae) and non-luminescent (Cantharidae) species of Neotropical beetles. We used the RNA-Seq data to construct a calibrated phylogeny of Elateroidea superfamily using a large number of nuclear molecular markers. The results indicate Lampyridae and Phengodidae/Rhagophthalmidae as sister-groups, suggesting that the bioluminescence evolved later in Elateridae than other families (Lampyridae, Phengodidae, and Rhagophthalmidae), and indicating the Upper Cretaceous as the period for the main diversification of Elateroidea bioluminescent species.

Contrasting Metatrochal Behavior of Mollusc and Annelid Larvae and the Regulation of Feeding While Swimming

Molluscan veliger larvae and some annelid larvae capture particulate food between a preoral prototrochal band of long cilia that create a current for both swimming and feeding and a postoral metatrochal band of shorter cilia that beat toward the prototroch. Larvae encountering satiating or noxious particles must somehow swim without capturing particles or else reject large numbers of captured particles. Because high rates of particle capture are inferred to depend on the beat of both ciliary bands, arrest of the metatroch could be one way to swim while reducing captures. Larvae in eight families of annelids arrest metatrochal cilia frequently during prototrochal beat, often over a large part of the metatrochal band and with the arrested cilia aligned near the beginning of the effective stroke. In contrast, metatrochs of veligers of gastropods and bivalves rarely arrested while the prototroch beat, and those arrests were more localized and variable in position. This difference in metatrochal arrest was unexpected under hypotheses of either a single origin of this feeding mechanism or multiple origins within each phylum. Although different in metatrochal arrests, larvae of both phyla can separate swimming from feeding while both prototroch and metatroch beat. One hypothesis explaining low rates of capture per encounter, without metatrochal arrest, is a change in adhesion of prototrochal cilia with algae. In a few observations, part of the velar edge was retained within the veliger's shell so that exposed prototrochal cilia contributed to swimming while the adjacent metatroch and food groove were sequestered.

Signal, Uncertainty, and Conflict in Phylogenomic Data for a Diverse Lineage of Microbial Eukaryotes (Diatoms, Bacillariophyta)

Diatoms (Bacillariophyta) are a species-rich group of eukaryotic microbes diverse in morphology, ecology, and metabolism. Previous reconstructions of the diatom phylogeny based on one or a few genes have resulted in inconsistent resolution or low support for critical nodes. We applied phylogenetic paralog pruning techniques to a data set of 94 diatom genomes and transcriptomes to infer perennially difficult species relationships, using concatenation and summary-coalescent methods to reconstruct species trees from data sets spanning a wide range of thresholds for taxon and column occupancy in gene alignments. Conflicts between gene and species trees decreased with both increasing taxon occupancy and bootstrap cutoffs applied to gene trees. Concordance between gene and species trees was lowest for short internodes and increased logarithmically with increasing edge length, suggesting that incomplete lineage sorting disproportionately affects species tree inference at short internodes, which are a common feature of the diatom phylogeny. Although species tree topologies were largely consistent across many data treatments, concatenation methods appeared to outperform summary-coalescent methods for sparse alignments. Our results underscore that approaches to species-tree inference based on few loci are likely to be misled by unrepresentative sampling of gene histories, particularly in lineages that may have diversified rapidly. In addition, phylogenomic studies of diatoms, and potentially other hyperdiverse groups, should maximize the number of gene trees with high taxon occupancy, though there is clearly a limit to how many of these genes will be available.

SNP hot-spots in the clam parasite QPX

Background

Quahog Parasite Unknown (QPX) is an opportunistic protistan pathogen of the clam Mercenaria mercenaria. Infections with QPX have caused significant economic losses in the Northeastern United States. Previous research demonstrated a geographic gradient for disease prevalence and intensity, but little information is available on the genetic diversity of the parasite throughout its distribution range. Also, QPX virulence factors are not well understood. This study addresses the occurrence of QPX genetic variants with a particular focus on functions involved in virulence and adaptation to environmental conditions.

Results

Analyses were performed using transcriptome-wide single-nucleotide polymorphism (SNP) of four QPX isolates cultured from infected clams collected from disparate locations along the Northeastern United States. For contig assembly and mapping, two different genome builds and four transcriptomes of the parasite were examined. Genomic variants appeared at a differential rate relative to sequenced transcripts at 20.18 and 22.55% occurrence under 1000 base pairs upstream and downstream protein domains respectively and at 57.26% rate in protein domain coding sequences. QPX strains shared 30.50% of the mutations and exhibited a preferential nucleotide substitution towards thymine. Sequence identity suggested relatedness between different QPX strains, with the parasite being possibly introduced to Virginia from the Massachusetts region during clam trading, while QPX could have been naturally present in New York. Diversity in virulence, temperature, and salinity domains suggested a common variability between strains, but with a preferential higher variation in local adaptation genes. This could explain differences in disease prevalence noted in different regions. Overall, the results supported views that this opportunistic parasite might be able to adapt to varying environmental conditions.

Conclusion

Relatedness and mutations between the four QPX strains suggested that variability in environmental-related functions favors parasite survival, potentially promoting resilience against stressful conditions. These findings are in agreement with the widespread presence of QPX in the environment. Although QPX levels are enzootic in most areas, an increase in disease outbreaks were often associated with seasonal changes in environmental conditions. A selection mediated by the parasitic life of QPX remains possible, but the effect of the environment on the biology of the parasite appears more obvious.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4866-8) contains supplementary material, which is available to authorized users.

In-depth proteomic analyses of Haliotis laevigata (greenlip abalone) nacre and prismatic organic shell matrix

Background

The shells of various Haliotis species have served as models of invertebrate biomineralization and physical shell properties for more than 20 years. A focus of this research has been the nacreous inner layer of the shell with its conspicuous arrangement of aragonite platelets, resembling in cross-section a brick-and-mortar wall. In comparison, the outer, less stable, calcitic prismatic layer has received much less attention. One of the first molluscan shell proteins to be characterized at the molecular level was Lustrin A, a component of the nacreous organic matrix of Haliotis rufescens. This was soon followed by the C-type lectin perlucin and the growth factor-binding perlustrin, both isolated from H. laevigata nacre, and the crystal growth-modulating AP7 and AP24, isolated from H. rufescens nacre. Mass spectrometry-based proteomics was subsequently applied to to Haliotis biomineralization research with the analysis of the H. asinina shell matrix and yielded 14 different shell-associated proteins. That study was the most comprehensive for a Haliotis species to date.

Methods

The shell proteomes of nacre and prismatic layer of the marine gastropod Haliotis laevigata were analyzed combining mass spectrometry-based proteomics and next generation sequencing.

Results

We identified 297 proteins from the nacreous shell layer and 350 proteins from the prismatic shell layer from the green lip abalone H. laevigata. Considering the overlap between the two sets we identified a total of 448 proteins. Fifty-one nacre proteins and 43 prismatic layer proteins were defined as major proteins based on their abundance at more than 0.2% of the total. The remaining proteins occurred at low abundance and may not play any significant role in shell fabrication. The overlap of major proteins between the two shell layers was 17, amounting to a total of 77 major proteins.

Conclusions

The H. laevigata shell proteome shares moderate sequence similarity at the protein level with other gastropod, bivalve and more distantly related invertebrate biomineralising proteomes. Features conserved in H. laevigata and other molluscan shell proteomes include short repetitive sequences of low complexity predicted to lack intrinsic three-dimensional structure, and domains such as tyrosinase, chitin-binding, and carbonic anhydrase. This catalogue of H. laevigata shell proteins represents the most comprehensive for a haliotid and should support future efforts to elucidate the molecular mechanisms of shell assembly.

Electronic supplementary material

The online version of this article (10.1186/s12953-018-0139-3) contains supplementary material, which is available to authorized users.

Nervous system development in the Pacific oyster, Crassostrea gigas (Mollusca: Bivalvia)

Background

Bivalves comprise a large, highly diverse taxon of invertebrate species. Developmental studies of neurogenesis among species of Bivalvia are limited. Due to a lack of neurogenesis information, it is difficult to infer a ground pattern for Bivalvia. To provide more comprehensive morphogenetic data on bivalve molluscs and relationships among molluscan clades, we investigated neurogenesis in the Pacific oyster, Crassostrea gigas, from the appearance of the first sensory cells to the formation of the larval ganglionic nervous system by co-immunocytochemistry of the neuronal markers FMRFamide or 5-HT and vesicular acetylcholine transporter (VAChT).

Results

Neurogenesis begins with the emergence of the apical serotonin-immunoreactive (5-HT-ir) sensory cells and paired sensory posttrochal dorsal and ventral FMRFamide-immunoreactive (FMRFamide-ir) cells at the early trochophore stage. Later, at the early veliger stage, the apical organ (AO) includes 5-HT-ir, FMRFamide-ir, and VAChT-ir cells. At the same stage, VAChT-ir cells appear in the posterior region of larvae and send axons towards the AO. Thus, FMRFamide-ir neurites and VAChT-ir processes form scaffolds for longitudinal neurite bundles develop into the paired ventral nerve cords (VNC). Later-appearing axons from the AO/CG neurons join the neurite bundles comprising the VNC. All larval ganglia appear along the VNC as paired or fused (epiathroid) clusters in late veliger and pediveliger larvae. We observed the transformation of the AO into the cerebral ganglia, which abundantly innervated the velum, and the transformation of ventral neurons into the pedal ganglia, innervating the foot, gills, and anterior adductor muscle. The visceral ganglia appear last in the pediveliger oyster and innervate the visceral mass and posterior adductor of premetamorphic larvae. In addition, a local FMRFamide-ir network was detected in the digestive system of pediveliger larvae. We identified VAChT-ir nervous elements in oyster larvae, which have not been observed previously in molluscs. Finally, we performed a morphology-based comparative analysis of neuronal structures among bivalve, conchiferan, and aculiferan species.

Conclusions

We described the development of the nervous system during the larval development in Crassostrea gigas. These data greatly advance the currently limited understanding of neurodevelopment in bivalves and mollusks, which has hampered the generation of a ground pattern reconstruction of the last common ancestor of Mollusca. Our morphological data support phylogenomic data indicating a closer Bivalvia-Gastropoda sister group relationship than the Bivalvia-Scaphopoda (Diasoma) group relationship.

Electronic supplementary material

The online version of this article (10.1186/s12983-018-0259-8) contains supplementary material, which is available to authorized users.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

Towards a ground pattern reconstruction of bivalve nervous systems: neurogenesis in the zebra mussel Dreissena polymorpha

Bivalvia is a taxon of aquatic mollusks that includes clams, oysters, mussels, and scallops. Within heterodont bivalves, Dreissena polymorpha is a small, mytiliform, freshwater mussel that develops indirectly via a planktotrophic veliger larva. Currently, only a few studies on bivalve neurogenesis are available, impeding the reconstruction of a ground pattern in Bivalvia. In order to inject novel data into this discussion, we describe herein the development of the serotonin-like and α-tubulin-like immunoreactive (lir) neuronal components of D. polymorpha from the early trochophore to the late veliger stage. Neurogenesis starts in the early trochophore stage at the apical pole with the appearance of one flask-shaped serotonin-lir cell. When larvae reach the veliger stage, four flask-shaped serotonin-lir cells are present in the apical organ. At the same time, the anlagen of the cerebral ganglia start to form at the base of the apical organ. From the apical organ, one pair of cerebro-visceral connectives projects posteriorly and connects to a posterior larval sensory organ that contains serotonin- and α-tubulin-like flask-shaped cells. Additional, paired serotonin-lir neurites originate from the apical organ and project into the velum. One unpaired stomatogastric serotonin-lir cell develops ventrally to the stomach at the veliger stage. The low number of serotonin-lir cells in the apical organ of bivalve veligers is shared with larvae of basally branching gastropods and scaphopods and is thus considered a feature of the last common ancestor of Conchifera, while the overall simplicity of the larval neural architecture appears to be a specific trait of Bivalvia.

Transcriptomic responses of the endangered freshwater mussel Margaritifera margaritifera to trace metal contamination in the Dronne River, France

The freshwater pearl mussel Margaritifera margaritifera is one of the most threatened freshwater bivalves worldwide. In this study, we aimed (i) to study the processes by which water quality might affect freshwater mussels in situ and (ii) to provide insights into the ecotoxicological significance of water pollution to natural populations in order to provide necessary information to enhance conservation strategies. M. margaritifera specimens were sampled in two close sites located upstream or downstream from an illegal dumping site. The renal transcriptome of these animals was assembled and gene transcription determined by RNA-seq. Correlations between transcription levels of each single transcript and the bioaccumulation of nine trace metals, age (estimated by sclerochronology), and condition index were determined in order to identify genes likely to respond to a specific factor. Amongst the studied metals, Cr, Zn, Cd, and Ni were the main factors correlated with transcription levels, with effects on translation, apoptosis, immune response, response to stimulus, and transport pathways. However, the main factor explaining changes in gene transcription appeared to be the age of individuals with a negative correlation with the transcription of retrotransposon-related genes. To investigate this effect further, mussels were classified into three age classes. In young, middle-aged and old animals, transcription levels were mainly explained by Cu, Zn and age, respectively. This suggests differences in the molecular responses of this species to metals during its lifetime that must be better assessed in future ecotoxicology studies.

Co-option of bacteriophage lysozyme genes by bivalve genomes

Eukaryotes have occasionally acquired genetic material through horizontal gene transfer (HGT). However, little is known about the evolutionary and functional significance of such acquisitions. Lysozymes are ubiquitous enzymes that degrade bacterial cell walls. Here, we provide evidence that two subclasses of bivalves (Heterodonta and Palaeoheterodonta) acquired a lysozyme gene via HGT, building on earlier findings. Phylogenetic analyses place the bivalve lysozyme genes within the clade of bacteriophage lysozyme genes, indicating that the bivalves acquired the phage-type lysozyme genes from bacteriophages, either directly or through intermediate hosts. These bivalve lysozyme genes underwent dramatic structural changes after their co-option, including intron gain and fusion with other genes. Moreover, evidence suggests that recurrent gene duplication occurred in the bivalve lysozyme genes. Finally, we show the co-opted lysozymes exhibit a capacity for antibacterial action, potentially augmenting the immune function of related bivalves. This represents an intriguing evolutionary strategy in the eukaryote–microbe arms race, in which the genetic materials of bacteriophages are co-opted by eukaryotes, and then used by eukaryotes to combat bacteria, using a shared weapon against a common enemy.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.