Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Phylogenomic analyses shed light on the relationships of chiton superfamilies and shell-eye evolution

Mollusca is the second-largest animal phylum with over 100,000 extant species representing eight classes. Across 1000 extant species in the class Polyplacophora, chitons have a relatively constrained morphology but with some notable deviations. Several genera possess "shell eyes", i.e., true eyes with a lens and retina that are embedded within the dorsal shells. The phylogeny of the major chiton clades is mostly well established, in a set of superfamily-level and higher level taxa supported by various approaches, including morphological studies, multiple gene markers, mitogenome-phylogeny, and phylotranscriptomic approaches. However, one critical lineage has remained unclear, namely Schizochiton which was controversially suggested as being the potential independent origin of chiton shell eyes. Here, with the draft genome sequencing of Schizochiton incisus (superfamily Schizochitonoidea) plus assemblies of transcriptome data from other polyplacophorans, we present phylogenetic reconstructions using both mitochondrial genomes and phylogenomic approaches with multiple methods. We found that phylogenetic trees from mitogenomic data are inconsistent, reflecting larger scale confounding factors in molluscan mitogenomes. However, a consistent and robust topology was generated with protein-coding genes using different models and methods. Our results support Schizochitonoidea as the sister group to other Chitonoidea in Chitonina, in agreement with the established classification. Combined with evidence from fossils, our phylogenetic results suggest that the earliest origin of shell eyes is in Schizochitonoidea, and that these structures were also gained secondarily in other genera in Chitonoidea. Our results have generated a holistic review of the internal relationship within Polyplacophora, and a better understanding of the evolution of Polyplacophora.

Octopod Hox genes and cephalopod plesiomorphies

Few other invertebrates captivate our attention as cephalopods do. Octopods, cuttlefish, and squids amaze with their behavior and sophisticated body plans that belong to the most intriguing among mollusks. Little is, however, known about their body plan formation and the role of Hox genes. The latter homeobox genes pattern the anterior-posterior body axis and have only been studied in a single decapod species so far. Here, we study developmental Hox and ParaHox gene expression in Octopus vulgaris. Hox genes are expressed in a near-to-staggered fashion, among others in homologous organs of cephalopods such as the stellate ganglia, the arms, or funnel. As in other mollusks Hox1 is expressed in the nascent octopod shell rudiment. While ParaHox genes are expressed in an evolutionarily conserved fashion, Hox genes are also expressed in some body regions that are considered homologous among mollusks such as the cephalopod arms and funnel with the molluscan foot. We argue that cephalopod Hox genes are recruited to a lesser extent into the formation of non-related organ systems than previously thought and emphasize that despite all morphological innovations molecular data still reveal the ancestral molluscan heritage of cephalopods.

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

RNA-Seq comparative study reveals molecular effectors linked to the resistance of Pinna nobilis to Haplosporidium pinnae parasite

With the intensification of maritime traffic, recently emerged infectious diseases have become major drivers in the decline and extinction of species. Since 2016, mass mortality events have decimated the endemic Mediterranean Sea bivalve Pinna nobilis, affecting ca. 100% of individuals. These events have largely been driven by Haplosporidium pinnae's infection, an invasive species which was likely introduced by shipping. While monitoring wild populations of P. nobilis, we observed individuals that survived such a mass mortality event during the summer of 2018 (France). We considered these individuals resistant, as they did not show any symptoms of the disease, while the rest of the population in the area was devastated. Furthermore, the parasite was not detected when we conducted a PCR amplification of a species-specific fragment of the small subunit ribosomal DNA. In parallel, the transcriptomic analysis showed evidence of some parasite RNA indicating that the resistant individuals had been exposed to the parasite without proliferating. To understand the underlying mechanisms of resistance in these individuals, we compared their gene expression with that of susceptible individuals. We performed de novo transcriptome assembly and annotated the expressed genes. A comparison of the transcriptomes in resistant and susceptible individuals highlighted a gene expression signature of the resistant phenotype. We found significant differential expressions of genes involved in immunity and cell architecture. This data provides the first insights into how individuals escape the pathogenicity associated with infection.

Gotta Go Slow: Two Evolutionarily Distinct Annelids Retain a Common Hedgehog Pathway Composition, Outlining Its Pan-Bilaterian Core

Hedgehog signaling is one of the key regulators of morphogenesis, cell differentiation, and regeneration. While the Hh pathway is present in all bilaterians, it has mainly been studied in model animals such as Drosophila and vertebrates. Despite the conservatism of its core components, mechanisms of signal transduction and additional components vary in Ecdysozoa and Deuterostomia. Vertebrates have multiple copies of the pathway members, which complicates signaling implementation, whereas model ecdysozoans appear to have lost some components due to fast evolution rates. To shed light on the ancestral state of Hh signaling, models from the third clade, Spiralia, are needed. In our research, we analyzed the transcriptomes of two spiralian animals, errantial annelid Platynereis dumerilii (Nereididae) and sedentarian annelid Pygospio elegans (Spionidae). We found that both annelids express almost all Hh pathway components present in Drosophila and mouse. We performed a phylogenetic analysis of the core pathway components and built multiple sequence alignments of the additional key members. Our results imply that the Hh pathway compositions of both annelids share more similarities with vertebrates than with the fruit fly. Possessing an almost complete set of single-copy Hh pathway members, lophotrochozoan signaling composition may reflect the ancestral features of all three bilaterian branches.

Epigenetic machinery is functionally conserved in cephalopods

Background

Epigenetic regulatory mechanisms are divergent across the animal kingdom, yet these mechanisms are not well studied in non-model organisms. Unique features of cephalopods make them attractive for investigating behavioral, sensory, developmental, and regenerative processes, and recent studies have elucidated novel features of genome organization and gene and transposon regulation in these animals. However, it is not known how epigenetics regulates these interesting cephalopod features. We combined bioinformatic and molecular analysis of Octopus bimaculoides to investigate the presence and pattern of DNA methylation and examined the presence of DNA methylation and 3 histone post-translational modifications across tissues of three cephalopod species.

Results

We report a dynamic expression profile of the genes encoding conserved epigenetic regulators, including DNA methylation maintenance factors in octopus tissues. Levels of 5-methyl-cytosine in multiple tissues of octopus, squid, and bobtail squid were lower compared to vertebrates. Whole genome bisulfite sequencing of two regions of the brain and reduced representation bisulfite sequencing from a hatchling of O. bimaculoides revealed that less than 10% of CpGs are methylated in all samples, with a distinct pattern of 5-methyl-cytosine genome distribution characterized by enrichment in the bodies of a subset of 14,000 genes and absence from transposons. Hypermethylated genes have distinct functions and, strikingly, many showed similar expression levels across tissues while hypomethylated genes were silenced or expressed at low levels. Histone marks H3K27me3, H3K9me3, and H3K4me3 were detected at different levels across tissues of all species.

Conclusions

Our results show that the DNA methylation and histone modification epigenetic machinery is conserved in cephalopods, and that, in octopus, 5-methyl-cytosine does not decorate transposable elements, but is enriched on the gene bodies of highly expressed genes and could cooperate with the histone code to regulate tissue-specific gene expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01404-1.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

Evolutionary Divergence and Radula Diversification in Two Ecomorphs from an Adaptive Radiation of Freshwater Snails

(1) Background: Adaptive diversification of complex traits plays a pivotal role in the evolution of organismal diversity. In the freshwater snail genus Tylomelania, adaptive radiations were likely promoted by trophic specialization via diversification of their key foraging organ, the radula. (2) Methods: To investigate the molecular basis of radula diversification and its contribution to lineage divergence, we used tissue-specific transcriptomes of two sympatric Tylomelania sarasinorum ecomorphs. (3) Results: We show that ecomorphs are genetically divergent lineages with habitat-correlated abundances. Sequence divergence and the proportion of highly differentially expressed genes are significantly higher between radula transcriptomes compared to the mantle and foot. However, the same is not true when all differentially expressed genes or only non-synonymous SNPs are considered. Finally, putative homologs of some candidate genes for radula diversification (hh, arx, gbb) were also found to contribute to trophic specialization in cichlids and Darwin’s finches. (4) Conclusions: Our results are in line with diversifying selection on the radula driving Tylomelania ecomorph divergence and indicate that some molecular pathways may be especially prone to adaptive diversification, even across phylogenetically distant animal groups.

Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across ten species of North American true toads (Bufonidae). Upon assembly and annotation of the ten transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g. steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the ten focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

HES and Mox genes are expressed during early mesoderm formation in a mollusk with putative ancestral features

The mesoderm is considered the youngest of the three germ layers. Although its morphogenesis has been studied in some metazoans, the molecular components underlying this process remain obscure for numerous phyla including the highly diverse Mollusca. Here, expression of Hairy and enhancer of split (HES), Mox, and myosin heavy chain (MHC) was investigated in Acanthochitona fascicularis, a representative of Polyplacophora with putative ancestral molluscan features. While AfaMHC is expressed throughout myogenesis, AfaMox1 is only expressed during early stages of mesodermal band formation and in the ventrolateral muscle, an autapomorphy of the polyplacophoran trochophore. Comparing our findings to previously published data across Metazoa reveals Mox expression in the mesoderm in numerous bilaterians including gastropods, polychaetes, and brachiopods. It is also involved in myogenesis in molluscs, annelids, tunicates, and craniates, suggesting a dual role of Mox in mesoderm and muscle formation in the last common bilaterian ancestor. AfaHESC2 is expressed in the ectoderm of the polyplacophoran gastrula and later in the mesodermal bands and in putative neural tissue, whereas AfaHESC7 is expressed in the trochoblasts of the gastrula and during foregut formation. This confirms the high developmental variability of HES gene expression and demonstrates that Mox and HES genes are pleiotropic.

Development and Interrogation of a Transcriptomic Resource for the Giant Triton Snail (Charonia tritonis)

Gastropod molluscs are among the most abundant species that inhabit coral reef ecosystems. Many are specialist predators, along with the giant triton snail Charonia tritonis (Linnaeus, 1758) whose diet consists of Acanthaster planci (crown-of-thorns starfish), a corallivore known to consume enormous quantities of reef-building coral. C. tritonis are considered vulnerable due to overexploitation, and a decline in their populations is believed to have contributed to recurring A. planci population outbreaks. Aquaculture is considered one approach that could help restore natural populations of C. tritonis and mitigate coral loss; however, numerous questions remain unanswered regarding their life cycle, including the molecular factors that regulate their reproduction and development. In this study, we have established a reference C. tritonis transcriptome derived from developmental stages (embryo and veliger) and adult tissues. This was used to identify genes associated with cell signalling, such as neuropeptides and G protein-coupled receptors (GPCRs), involved in endocrine and olfactory signalling. A comparison of developmental stages showed that several neuropeptide precursors are exclusively expressed in post-hatch veligers and functional analysis found that FFamide stimulated a significant (20.3%) increase in larval heart rate. GPCRs unique to veligers, and a diversity of rhodopsin-like GPCRs located within adult cephalic tentacles, all represent candidate olfactory receptors. In addition, the cytochrome P450 superfamily, which participates in the biosynthesis and degradation of steroid hormones and lipids, was also found to be expanded with at least 91 genes annotated, mostly in gill tissue. These findings further progress our understanding of C. tritonis with possible application in developing aquaculture methods.

The genome of the venomous snail Lautoconus ventricosus sheds light on the origin of conotoxin diversity

BACKGROUND

Venoms are deadly weapons to subdue prey or deter predators that have evolved independently in many animal lineages. The genomes of venomous animals are essential to understand the evolutionary mechanisms involved in the origin and diversification of venoms.

RESULTS

Here, we report the chromosome-level genome of the venomous Mediterranean cone snail, Lautoconus ventricosus (Caenogastropoda: Conidae). The total size of the assembly is 3.59 Gb; it has high contiguity (N50 = 93.53 Mb) and 86.6 Mb of the genome assembled into the 35 largest scaffolds or pseudochromosomes. On the basis of venom gland transcriptomes, we annotated 262 complete genes encoding conotoxin precursors, hormones, and other venom-related proteins. These genes were scattered in the different pseudochromosomes and located within repetitive regions. The genes encoding conotoxin precursors were normally structured into 3 exons, which did not necessarily coincide with the 3 structural domains of the corresponding proteins. Additionally, we found evidence in the L. ventricosus genome for a past whole-genome duplication event by means of conserved gene synteny with the Pomacea canaliculata genome, the only one available at the chromosome level within Caenogastropoda. The whole-genome duplication event was further confirmed by the presence of a duplicated hox gene cluster. Key genes for gastropod biology including those encoding proteins related to development, shell formation, and sex were located in the genome.

CONCLUSIONS

The new high-quality L. ventricosus genome should become a reference for assembling and analyzing new gastropod genomes and will contribute to future evolutionary genomic studies among venomous animals.

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

Constructing a de novo transcriptome and a reference proteome for the bivalve Scrobicularia plana: Comparative analysis of different assembly strategies and proteomic analysis

Scrobicularia plana is a coastal and estuarine bivalve widely used in ecotoxicological studies. However, the underlying molecular mechanisms for S. plana pollutant responses are hardly known due to the lack of molecular databases. Thus, in this study we present a holistic approach to assess a robust reference transcriptome and proteome of this clam. A mixture of control and metal-exposed individuals was used for mRNA isolation. Four sets of high quality filtered preprocessed reads were generated (two quality scores and two sequenced lengths) and assembled with Mira, Ray and Trinity algorithms. The sixty-four generated assemblies were refined, filtered and evaluated for their proteomic quality. Eight assemblies presented top Detonate scores but one was selected due to its compactness and biological representation, which was generated: (i) from the highest quality dataset (Q20L100), (ii) using Trinity algorithm with all k-mers (AtKa), (iii) removing redundancy by CD-HIT (RR80), and (iv) filtering out poor contigs (F), that was subsequently named Q20L100AtKaRR80F. S. plana proteomic analysis revealed 10,017 peptide groups that corresponded to 2066 proteins with a wide coverage of molecular functions and biological processes, confirming the strength of the database generated.

Non-collinear Hox gene expression in bivalves and the evolution of morphological novelties in mollusks

Hox genes are key developmental regulators that are involved in establishing morphological features during animal ontogeny. They are commonly expressed along the anterior-posterior axis in a staggered, or collinear, fashion. In mollusks, the repertoire of body plans is widely diverse and current data suggest their involvement during development of landmark morphological traits in Conchifera, one of the two major lineages that comprises those taxa that originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia). For most clades, and bivalves in particular, data on Hox gene expression throughout ontogeny are scarce. We thus investigated Hox expression during development of the quagga mussel, Dreissena rostriformis, to elucidate to which degree they might contribute to specific phenotypic traits as in other conchiferans. The Hox/ParaHox complement of Mollusca typically comprises 14 genes, 13 of which are present in bivalve genomes including Dreissena. We describe here expression of 9 Hox genes and the ParaHox gene Xlox during Dreissena development. Hox expression in Dreissena is first detected in the gastrula stage with widely overlapping expression domains of most genes. In the trochophore stage, Hox gene expression shifts towards more compact, largely mesodermal domains. Only few of these domains can be assigned to specific developing morphological structures such as Hox1 in the shell field and Xlox in the hindgut. We did not find traces of spatial or temporal staggered expression of Hox genes in Dreissena. Our data support the notion that Hox gene expression has been coopted independently, and to varying degrees, into lineage-specific structures in the respective conchiferan clades. The non-collinear mode of Hox expression in Dreissena might be a result of the low degree of body plan regionalization along the bivalve anterior-posterior axis as exemplified by the lack of key morphological traits such as a distinct head, cephalic tentacles, radula apparatus, and a simplified central nervous system.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Computationally predicted gene regulatory networks in molluscan biomineralization identify extracellular matrix production and ion transportation pathways

MOTIVATION

The molecular processes regulating molluscan shell production remain relatively uncharacterized, despite the clear evolutionary and societal importance of biomineralization.

RESULTS

Here we built the first computationally predicted gene regulatory network (GRN) for molluscan biomineralization using Antarctic clam (Laternula elliptica) mantle gene expression data produced over an age-categorized shell damage-repair time-course. We used previously published in vivo in situ hybridization expression data to ground truth gene interactions predicted by the GRN and show that candidate biomineralization genes from different shell layers, and hence microstructures, were connected in unique modules. We characterized two biomineralization modules of the GRN and hypothesize that one module is responsible for translating the extracellular proteins required for growing, repairing or remodelling the nacreous shell layer, whereas the second module orchestrates the transport of both ions and proteins to the shell secretion site, which are required during normal shell growth, and repair. Our findings demonstrate that unbiased computational methods are particularly valuable for studying fundamental biological processes and gene interactions in non-model species where rich sources of gene expression data exist, but annotation rates are poor and the ability to carry out true functional tests are still lacking.

AVAILABILITY AND IMPLEMENTATION

The raw RNA-Seq data is freely available for download from NCBI SRA (Accession: PRJNA398984), the assembled and annotated transcriptome can be viewed and downloaded from molluscDB (ensembl.molluscdb.org) and in addition, the assembled transcripts, reconstructed GRN, modules and detailed annotations are all available as Supplementary Files.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The evo-devo of molluscs: Insights from a genomic perspective

Molluscs represent one of ancient and evolutionarily most successful groups of marine invertebrates, with a tremendous diversity of morphology, behavior, and lifestyle. Molluscs are excellent subjects for evo-devo studies; however, understanding of the evo-devo of molluscs has been largely hampered by incomplete fossil records and limited molecular data. Recent advancement of genomics and other technologies has greatly fueled the molluscan "evo-devo" field, and decoding of several molluscan genomes provides unprecedented insights into molluscan biology and evolution. Here, we review the recent progress of molluscan genome sequencing as well as novel insights gained from their genomes, by emphasizing how molluscan genomics enhances our understanding of the evo-devo of molluscs.

Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae

The transient larva-bearing biphasic life cycle is the hallmark of many metazoan phyla, but how metazoan larvae originated remains a major enigma in animal evolution. There are two hypotheses for larval origin. The 'larva-first' hypothesis suggests that the first metazoans were similar to extant larvae, with later evolution of the adult-added biphasic life cycle; the 'adult-first' hypothesis suggests that the first metazoans were adult forms, with the biphasic life cycle arising later via larval intercalation. Here, we investigate the evolutionary origin of primary larvae by conducting ontogenetic transcriptome profiling for Mollusca-the largest marine phylum characterized by a trochophore larval stage and highly variable adult forms. We reveal that trochophore larvae exhibit rapid transcriptome evolution with extraordinary incorporation of novel genes (potentially contributing to adult shell evolution), and that cell signalling/communication genes (for example, caveolin and innexin) are probably crucial for larval evolution. Transcriptome age analysis of eight metazoan species reveals the wide presence of young larval transcriptomes in both trochozoans and other major metazoan lineages, therefore arguing against the prevailing larva-first hypothesis. Our findings support an adult-first evolutionary scenario with a single metazoan larval intercalation, and suggest that the first appearance of proto-larva probably occurred after the divergence of direct-developing Ctenophora from a metazoan ancestor.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Remnants of ancestral larval eyes in an eyeless mollusk? Molecular characterization of photoreceptors in the scaphopod Antalis entalis

Background

Eyes have evolved and been lost multiple times during animal evolution, however, the process of eye loss has only been reconstructed in a few cases. Mollusks exhibit eyes as varied as the octopod camera eye or the gastropod cup eye and are ideal systems for studying the evolution of eyes, photoreceptors, and opsins.

Results

Here, we identify genes related to photoreceptor formation and function in an eyeless conchiferan mollusk, the scaphopod Antalis entalis, and investigate their spatial and temporal expression patterns during development. Our study reveals that the scaphopod early mid-stage trochophore larva has putative photoreceptors in a similar location and with a similar gene expression profile as the trochophore of polyplacophoran mollusks. The apical and post-trochal putative photoreceptors appear to co-express go-opsin, six1/2, myoV, and eya, while expression domains in the posterior foot and pavilion (posterior mantle opening) show co-expression of several other candidate genes but not go-opsin. Sequence analysis reveals that the scaphopod Go-opsin amino acid sequence lacks the functionally important lysine (K296; Schiff base) in the retinal-binding domain, but has not accumulated nonsense mutations and still exhibits the canonical G-protein activation domain.

Conclusions

The scaphopod Go-opsin sequence reported here is the only known example of a bilaterian opsin that lacks lysine K296 in the retinal-binding domain. Although this may render the Go-opsin unable to detect light, the protein may still perform sensory functions. The location, innervation, development, and gene expression profiles of the scaphopod and polyplacophoran apical and post-trochal photoreceptors suggest that they are homologous, even though the scaphopod post-trochal photoreceptors have degenerated. This indicates that post-trochal eyes are not a polyplacophoran apomorphy but likely a molluscan synapomorphy lost in other mollusks. Scaphopod eye degeneration is probably a result of the transition to an infaunal life history and is reflected in the likely functional degeneration of Go-opsin, the loss of photoreceptor shielding pigments, and the scarce expression of genes involved in phototransduction and eye development. Our results emphasize the importance of studying a phylogenetically broad range of taxa to infer the mechanisms and direction of body plan evolution.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

Functional and evolutionary perspectives on gill structures of an obligate air-breathing, aquatic snail

Ampullariids are freshwater gastropods bearing a gill and a lung, thus showing different degrees of amphibiousness. In particular, Pomacea canaliculata (Caenogastropoda, Ampullariidae) is an obligate air-breather that relies mainly or solely on the lung for dwelling in poorly oxygenated water, for avoiding predators, while burying in the mud during aestivation, and for oviposition above water level. In this paper, we studied the morphological peculiarities of the gill in this species. We found (1) the gill and lung vasculature and innervation are intimately related, allowing alternation between water and air respiration; (2) the gill epithelium has features typical of a transporting rather than a respiratory epithelium; and (3) the gill has resident granulocytes within intraepithelial spaces that may serve a role for immune defence. Thus, the role in oxygen uptake may be less significant than the roles in ionic/osmotic regulation and immunity. Also, our results provide a morphological background to understand the dependence on aerial respiration of Pomacea canaliculata. Finally, we consider these findings from a functional perspective in the light of the evolution of amphibiousness in the Ampullariidae, and discuss that master regulators may explain the phenotypic convergence of gill structures amongst this molluscan species and those in other phyla.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

Extensive conservation of the proneuropeptide and peptide prohormone complement in mollusks

As one of the most diverse groups of invertebrate animals, mollusks represent powerful models for neurobiological and developmental studies. Neuropeptides and peptide hormones are a heterogeneous class of signalling molecules involved in chemical communication between neurons and in neuroendocrine regulation. Here we present a fine-grained view of the molluscan neuropeptide and peptide hormone toolkit. Our results expand the distribution of several peptide families (e.g., prokineticin, insulin-related peptides, prohormone-4, LFRFamide) within Lophotrochozoa and provide evidence for an early origin of others (e.g., GNXQN/prohormone-2, neuroparsin). We identified a new peptide family broadly distributed among conchiferan mollusks, the PXRX family. We found the Wnt antagonist dickkopf1/2/4 ortholog in lophotrochozoans and nematodes and reveal that the egg-laying hormone family is a DH44 homolog restricted to gastropods. Our data demonstrate that numerous peptides evolved much earlier than previously assumed and that key signalling elements are extensively conserved among extant mollusks.

The evolution of molluscs

Molluscs are extremely diverse invertebrate animals with a rich fossil record, highly divergent life cycles, and considerable economical and ecological importance. Key representatives include worm-like aplacophorans, armoured groups (e.g. polyplacophorans, gastropods, bivalves) and the highly complex cephalopods. Molluscan origins and evolution of their different phenotypes have largely remained unresolved, but significant progress has been made over recent years. Phylogenomic studies revealed a dichotomy of the phylum, resulting in Aculifera (shell-less aplacophorans and multi-shelled polyplacophorans) and Conchifera (all other, primarily uni-shelled groups). This challenged traditional hypotheses that proposed that molluscs gradually evolved complex phenotypes from simple, worm-like animals, a view that is corroborated by developmental studies that showed that aplacophorans are secondarily simplified. Gene expression data indicate that key regulators involved in anterior-posterior patterning (the homeobox-containing Hox genes) lost this function and were co-opted into the evolution of taxon-specific novelties in conchiferans. While the bone morphogenetic protein (BMP)/decapentaplegic (Dpp) signalling pathway, that mediates dorso-ventral axis formation, and molecular components that establish chirality appear to be more conserved between molluscs and other metazoans, variations from the common scheme occur within molluscan sublineages. The deviation of various molluscs from developmental pathways that otherwise appear widely conserved among metazoans provides novel hypotheses on molluscan evolution that can be tested with genome editing tools such as the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeats-associated protein9) system.

Staggered Hox expression is more widespread among molluscs than previously appreciated

Hox genes are expressed along the anterior-posterior body axis in a colinear fashion in the majority of bilaterians. Contrary to polyplacophorans, a group of aculiferan molluscs with conserved ancestral molluscan features, gastropods and cephalopods deviate from this pattern by expressing Hox genes in distinct morphological structures and not in a staggered fashion. Among conchiferans, scaphopods exhibit many similarities with gastropods, cephalopods and bivalves, however, the molecular developmental underpinnings of these similar traits remain unknown. We investigated Hox gene expression in developmental stages of the scaphopod Antalis entalis to elucidate whether these genes are involved in patterning morphological traits shared by their kin conchiferans. Scaphopod Hox genes are predominantly expressed in the foot and mantle but also in the central nervous system. Surprisingly, the scaphopod mid-stage trochophore exhibits a near-to staggered expression of all nine Hox genes identified. Temporal colinearity was not found and early-stage and late-stage trochophores, as well as postmetamorphic individuals, do not show any apparent traces of staggered expression. In these stages, Hox genes are expressed in distinct morphological structures such as the cerebral and pedal ganglia and in the shell field of early-stage trochophores. Interestingly, a re-evaluation of previously published data on early-stage cephalopod embryos and of the gastropod pre-torsional veliger shows that these developmental stages exhibit traces of staggered Hox expression. Considering our results and all gene expression and genomic data available for molluscs as well as other bilaterians, we suggest a last common molluscan ancestor with colinear Hox expression in predominantly ectodermal tissues along the anterior-posterior axis. Subsequently, certain Hox genes have been co-opted into the patterning process of distinct structures (apical organ or prototroch) in conchiferans.

AmpuBase: a transcriptome database for eight species of apple snails (Gastropoda: Ampullariidae)

BACKGROUND

Gastropoda, with approximately 80,000 living species, is the largest class of Mollusca. Among gastropods, apple snails (family Ampullariidae) are globally distributed in tropical and subtropical freshwater ecosystems and many species are ecologically and economically important. Ampullariids exhibit various morphological and physiological adaptations to their respective habitats, which make them ideal candidates for studying adaptation, population divergence, speciation, and larger-scale patterns of diversity, including the biogeography of native and invasive populations. The limited availability of genomic data, however, hinders in-depth ecological and evolutionary studies of these non-model organisms.

RESULTS

Using Illumina Hiseq platforms, we sequenced 1220 million reads for seven species of apple snails. Together with the previously published RNA-Seq data of two apple snails, we conducted de novo transcriptome assembly of eight species that belong to five genera of Ampullariidae, two of which represent Old World lineages and the other three New World lineages. There were 20,730 to 35,828 unigenes with predicted open reading frames for the eight species, with N50 (shortest sequence length at 50% of the unigenes) ranging from 1320 to 1803 bp. 69.7% to 80.2% of these unigenes were functionally annotated by searching against NCBI's non-redundant, Gene Ontology database and the Kyoto Encyclopaedia of Genes and Genomes. With these data we developed AmpuBase, a relational database that features online BLAST functionality for DNA/protein sequences, keyword searching for unigenes/functional terms, and download functions for sequences and whole transcriptomes.

CONCLUSIONS

In summary, we have generated comprehensive transcriptome data for multiple ampullariid genera and species, and created a publicly accessible database with a user-friendly interface to facilitate future basic and applied studies on ampullariids, and comparative molecular studies with other invertebrates.

De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species

Background

Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date.

Results

We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species.

Conclusion

New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.

The gonadal transcriptome of the unisexual Amazon molly Poecilia formosa in comparison to its sexual ancestors, Poecilia mexicana and Poecilia latipinna

Background

The unisexual Amazon molly (Poecilia formosa) originated from a hybridization between two sexual species, the sailfin molly (Poecilia latipinna) and the Atlantic molly (Poecilia mexicana). The Amazon molly reproduces clonally via sperm-dependent parthenogenesis (gynogenesis), in which the sperm of closely related species triggers embryogenesis of the apomictic oocytes, but typically does not contribute genetic material to the next generation. We compare for the first time the gonadal transcriptome of the Amazon molly to those of both ancestral species, P. mexicana and P. latipinna.

Results

We sequenced the gonadal transcriptomes of the P. formosa and its parental species P. mexicana and P. latipinna using Illumina RNA-sequencing techniques (paired-end, 100 bp). De novo assembly of about 50 million raw read pairs for each species was performed using Trinity, yielding 106,922 transcripts for P. formosa, 115,175 for P. latipinna, and 133,025 for P. mexicana after eliminating contaminations. On the basis of sequence similarity comparisons to other teleost species and the UniProt databases, functional annotation, and differential expression analysis, we demonstrate the similarity of the transcriptomes among the three species. More than 40% of the transcripts for each species were functionally annotated and about 70% were assigned to orthologous genes of a closely related species. Differential expression analysis between the sexual and unisexual species uncovered 2035 up-regulated and 564 down-regulated genes in P. formosa. This was exemplary validated for six genes by qRT-PCR.

Conclusions

We identified more than 130 genes related to meiosis and reproduction within the apomictically reproducing P. formosa. Overall expression of these genes seems to be down-regulated in the P. formosa transcriptome compared to both ancestral species (i.e., 106 genes down-regulated, 29 up-regulated). A further 35 meiosis and reproduction related genes were not found in the P. formosa transcriptome, but were only expressed in the sexual species. Our data support the hypothesis of general down-regulation of meiosis-related genes in the apomictic Amazon molly. Furthermore, the obtained dataset and identified gene catalog will serve as a resource for future research on the molecular mechanisms behind the reproductive mode of this unisexual species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4382-2) contains supplementary material, which is available to authorized users.

Approaches to Macroevolution: 1. General Concepts and Origin of Variation

Approaches to macroevolution require integration of its two fundamental components, i.e. the origin and the sorting of variation, in a hierarchical framework. Macroevolution occurs in multiple currencies that are only loosely correlated, notably taxonomic diversity, morphological disparity, and functional variety. The origin of variation within this conceptual framework is increasingly understood in developmental terms, with the semi-hierarchical structure of gene regulatory networks (GRNs, used here in a broad sense incorporating not just the genetic circuitry per se but the factors controlling the timing and location of gene expression and repression), the non-linear relation between magnitude of genetic change and the phenotypic results, the evolutionary potential of co-opting existing GRNs, and developmental responsiveness to nongenetic signals (i.e. epigenetics and plasticity), all requiring modification of standard microevolutionary models, and rendering difficult any simple definition of evolutionary novelty. The developmental factors underlying macroevolution create anisotropic probabilities-i.e., an uneven density distribution-of evolutionary change around any given phenotypic starting point, and the potential for coordinated changes among traits that can accommodate change via epigenetic mechanisms. From this standpoint, "punctuated equilibrium" and "phyletic gradualism" simply represent two cells in a matrix of evolutionary models of phenotypic change, and the origin of trends and evolutionary novelty are not simply functions of ecological opportunity. Over long timescales, contingency becomes especially important, and can be viewed in terms of macroevolutionary lags (the temporal separation between the origin of a trait or clade and subsequent diversification); such lags can arise by several mechanisms: as geological or phylogenetic artifacts, or when diversifications require synergistic interactions among traits, or between traits and external events. The temporal and spatial patterns of the origins of evolutionary novelties are a challenge to macroevolutionary theory; individual events can be described retrospectively, but a general model relating development, genetics, and ecology is needed. An accompanying paper (Jablonski in Evol Biol 2017) reviews diversity dynamics and the sorting of variation, with some general conclusions.

Hox and ParaHox gene expression in early body plan patterning of polyplacophoran mollusks

Molecular developmental studies of various bilaterians have shown that the identity of the anteroposterior body axis is controlled by Hox and ParaHox genes. Detailed Hox and ParaHox gene expression data are available for conchiferan mollusks, such as gastropods (snails and slugs) and cephalopods (squids and octopuses), whereas information on the putative conchiferan sister group, Aculifera, is still scarce (but see Fritsch et al., 2015 on Hox gene expression in the polyplacophoran Acanthochitona crinita). In contrast to gastropods and cephalopods, the Hox genes in polyplacophorans are expressed in an anteroposterior sequence similar to the condition in annelids and other bilaterians. Here, we present the expression patterns of the Hox genes Lox5, Lox4, and Lox2, together with the ParaHox gene caudal (Cdx) in the polyplacophoran A. crinita. To localize Hox and ParaHox gene transcription products, we also investigated the expression patterns of the genes FMRF and Elav, and the development of the nervous system. Similar to the other Hox genes, all three Acr‐Lox genes are expressed in an anteroposterior sequence. Transcripts of Acr‐Cdx are seemingly present in the forming hindgut at the posterior end. The expression patterns of both the central class Acr‐Lox genes and the Acr‐Cdx gene are strikingly similar to those in annelids and nemerteans. In Polyplacophora, the expression patterns of the Hox and ParaHox genes seem to be evolutionarily highly conserved, while in conchiferan mollusks these genes are co‐opted into novel functions that might have led to evolutionary novelties, at least in gastropods and cephalopods.