The purplish bifurcate mussel Mytilisepta virgata gene expression atlas reveals a remarkable tissue functional specialization

A genome-based phylogeny for Mollusca is concordant with fossils and morphology

Extreme morphological disparity within Mollusca has long confounded efforts to reconstruct a stable backbone phylogeny for the phylum. Familiar molluscan groups-gastropods, bivalves, and cephalopods-each represent a diverse radiation with myriad morphological, ecological, and behavioral adaptations. The phylum further encompasses many more unfamiliar experiments in animal body-plan evolution. In this work, we reconstructed the phylogeny for living Mollusca on the basis of metazoan BUSCO (Benchmarking Universal Single-Copy Orthologs) genes extracted from 77 (13 new) genomes, including multiple members of all eight classes with two high-quality genome assemblies for monoplacophorans. Our analyses confirm a phylogeny proposed from morphology and show widespread genomic variation. The flexibility of the molluscan genome likely explains both historic challenges with their genomes and their evolutionary success.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves

BACKGROUND

Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution.

RESULTS

We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta.

CONCLUSIONS

Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

Genome-wide sequencing identifies a thermal-tolerance related synonymous mutation in the mussel, Mytilisepta virgata

The roles of synonymous mutations for adapting to stressful thermal environments are of fundamental biological and ecological interests but poorly understood. To study whether synonymous mutations influence thermal adaptation at specific microhabitats, a genome-wide genotype-phenotype association analysis is carried out in the black mussels Mytilisepta virgata. A synonymous mutation of Ubiquitin-specific Peptidase 15 (MvUSP15) is significantly associated with the physiological upper thermal limit. The individuals carrying GG genotype (the G-type) at the mutant locus possess significantly lower heat tolerance compared to the individuals carrying GA and AA genotypes (the A-type). When heated to sublethal temperature, the G-type exhibit higher inter-individual variations in MvUSP15 expression, especially for the mussels on the sun-exposed microhabitats. Taken together, a synonymous mutation in MvUSP15 can affect the gene expression profile and interact with microhabitat heterogeneity to influence thermal resistance. This integrative study sheds light on the ecological importance of adaptive synonymous mutations as an underappreciated genetic buffer against heat stress and emphasizes the importance of integrative studies at a microhabitat scale for evaluating and predicting the impacts of climate change.

First Insights into the Repertoire of Secretory Lectins in Rotifers

Due to their high biodiversity and adaptation to a mutable and challenging environment, aquatic lophotrochozoan animals are regarded as a virtually unlimited source of bioactive molecules. Among these, lectins, i.e., proteins with remarkable carbohydrate-recognition properties involved in immunity, reproduction, self/nonself recognition and several other biological processes, are particularly attractive targets for biotechnological research. To date, lectin research in the Lophotrochozoa has been restricted to the most widespread phyla, which are the usual targets of comparative immunology studies, such as Mollusca and Annelida. Here we provide the first overview of the repertoire of the secretory lectin-like molecules encoded by the genomes of six target rotifer species: Brachionus calyciflorus, Brachionus plicatilis, Proales similis (class Monogononta), Adineta ricciae, Didymodactylos carnosus and Rotaria sordida (class Bdelloidea). Overall, while rotifer secretory lectins display a high molecular diversity and belong to nine different structural classes, their total number is significantly lower than for other groups of lophotrochozoans, with no evidence of lineage-specific expansion events. Considering the high evolutionary divergence between rotifers and the other major sister phyla, their widespread distribution in aquatic environments and the ease of their collection and rearing in laboratory conditions, these organisms may represent interesting targets for glycobiological studies, which may allow the identification of novel carbohydrate-binding proteins with peculiar biological properties.

Diverse Localization Patterns of an R-Type Lectin in Marine Annelids

Lectins facilitate cell–cell contact and are critical in many cellular processes. Studying lectins may help us understand the mechanisms underlying tissue regeneration. We investigated the localization of an R-type lectin in a marine annelid (Perinereis sp.) with remarkable tissue regeneration abilities. Perinereis nuntia lectin (PnL), a galactose-binding lectin with repeating Gln-X-Trp motifs, is derived from the ricin B-chain. An antiserum was raised against PnL to specifically detect a 32-kDa lectin in the crude extracts from homogenized lugworms. The antiserum detected PnL in the epidermis, setae, oblique muscle, acicula, nerve cord, and nephridium of the annelid. Some of these tissues and organs also produced Galactose (Gal) or N-acetylgalactosamine (GalNAc), which was detected by fluorescent-labeled plant lectin. These results indicated that the PnL was produced in the tissues originating from the endoderm, mesoderm, and ectoderm. Besides, the localizing pattern of PnL partially merged with the binding pattern of a fluorescent-labeled mushroom lectin that binds to Gal and GalNAc. It suggested that PnL co-localized with galactose-containing glycans in Annelid tissue; this might be the reason PnL needed to be extracted with haptenic sugar, such as d-galactose, in the buffer. Furthermore, we found that a fluorescein isothiocyanate-labeled Gal/GalNAc-binding mushroom lectin binding pattern in the annelid tissue overlapped with the localizing pattern of PnL. These findings suggest that lectin functions by interacting with Gal-containing glycoconjugates in the tissues.

[Cell Function Research of β-Trefoil Lectins from Mytilidae]

Two novel β-trefoil lectins, MytiLec-1 and SeviL were found from mussels in the coast of Yokohama and Nagasaki. MytiLec-1 was purified from gill and mantle of Mytilus galloprovincialis. It was consisted of 149 amino acid residues and there was no similarity with any other proteins when it was discovered. We advocate for this "Mytilectin" as a new protein family because of their novelty of its primary structure and homologues were also found in other mussels. Glycan array analysis revealed that MytiLec-1 specifically bound to the Gb3 and Gb4 glycan which contained the α-galactoside. MytiLec-1 caused the apoptosis against the Burkitt's lymphoma cells through the interaction of Gb3 express in their cell surface. On the other hand, SeviL obtained from gill and mantle of Mytilisepta virgata showed the specific binding against GM1b, asialo GM1 and SSEA-4 which are known as glycosphingolipid glycan including the β-galactoside. In addition, SeviL was identified as R type lectin by confirmation of QXW motif within its primary structure. Messenger RNA of SeviL like R type lectins was also found among the musssels including Mytilus galloprovincialis. SeviL also showed the apoptosis against asialo GM1 expressing cells. To apply the anticancer lectin as a novel molecular target drug, primary structure of MytiLec-1 was analyzed to enhance the stabilization of confirmation by computational design technique. It was succeeded to produce a monomeric artificial β-trefoil lectin, Mitsuba-1 without losing the Gb3 binding ability. Comparison of biological function between Mitsuba-1 and MytiLec-1 is also described in this study.

Host-Endosymbiont Genome Integration in a Deep-Sea Chemosymbiotic Clam

Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate coevolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 Gb, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene families, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 Mb but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host’s high dependence on the symbiont for nutrition. Overall, the host–symbiont genomes show not only tight metabolic complementarity but also distinct signatures of coevolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

Insights into Antimicrobial Peptides from Limacus flavus Mucus

Antimicrobial peptides have been a major research subject since the rise of antimicrobial resistance as a major public health problem. These molecules are considered a potential therapeutic source of antibiotics with broad-spectrum activity against microorganisms. Two antimicrobial peptides were isolated from the mucus of the Limacus flavus slug. The mucus was obtained by thermal shock, lyophilized and extracted with acetic acid. The supernatant was prefractionated in Sep-Pak and shortly thereafter fractionated by reverse-phase high-performance liquid chromatography. The manually obtained fractions were used in antimicrobial and cytotoxic assays and finally subjected to mass spectrometry (MS/MS). Characterization was performed by bioinformatics analysis with the tool Peaks^®X + and by comparison with the NCBI and UniProt-SwissProt databases. Additionally, the physicochemical parameters of the samples were evaluated with online programs. Two fractions comtained antimicrobial peptides with the ability to inhibit Micrococcus luteus A270; both samples, LFMP-001 and LFMP-002, were hydrophilic molecules consisting of fewer than 20 residues. Comparison of the SDS-PAGE and Peaks^®X + data showed that both had Mw < 3 kDa. In summary, this study presents data on the isolation and characterization of antimicrobial peptides from a slug and shows their potential against gram-positive bacteria.

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.

A year in the salt marsh: Seasonal changes in gill protein expression in the temperate intertidal mussel Geukensia demissa

Organisms living in temperate and polar regions experience extensive seasonal changes in the physical and biotic environment, including temperature, insolation, and food availability, among other factors. Sessile intertidal organisms respond to such seasonal fluctuations largely through physiological and biochemical means, because their behavioral responses are severely limited. In this study, we used a proteomic approach to examine changes in seasonal protein expression of gill from the intertidal mussel Geukensia demissa, a keystone species of the western Atlantic salt marsh, over the course of one year. Gill tissue of mussels collected in summer had the greatest number of proteins significantly increased in abundance (37 of 592 spots detected on two-dimensional polyacrylamide gels), although autumn mussels revealed a comparable proportion of up-regulated proteins (31 spots). In contrast, the number of proteins changing in abundance in winter and spring mussels were substantially smaller (15 and 9, respectively). Identification of these proteins revealed both expected and unanticipated changes to the proteome. Maintenance of gill cilia dominates in the summer when filter-feeding is most active, as evidenced by cytoskeletal proteins such as tektin-4 and tubulin isoforms; a signal of protection from heat stress is also present in summer (e.g., heat shock cognate 70). In autumn oxidative stress protection (peroxiredoxin-5 and manganese-containing superoxide dismutase) and aerobic ATP synthetic capacity (ATP synthase subunits a and delta) appear to increase. In winter a signal of cold-induced oxidative stress is apparent (Mn-SOD and NADP-dependent isocitrate dehydrogenase), perhaps in association with heavy metal toxicity and exposure to pathogens. Gill tissue from spring shows relatively little environmental acclimatization, other than a possible increase in protein synthesis capacity.

Purification and Functional Characterization of the Effects on Cell Signaling of Mytilectin: A Novel β-Trefoil Lectin from Marine Mussels

In the 2010s, a novel lectin family with β-trefoil folding has been identified in marine mussels from the family Mytilidae (phylum Mollusca). "MytiLec-1," the lectin described in this chapter, was the first member of this family to be isolated and characterized from the Mediterranean mussel Mytilus galloprovincialis, a commercially and ecologically important species, spread in marine coastal areas worldwide. MytiLec-1 bound to the sugar moiety of globotriose (Gb3: Galα1-4Galβ1-4Glc), an α-galactoside, leading to apoptosis of Gb3-expressing Burkitt's lymphoma cells. Although the primary structure of MytiLec-1 was quite unusual, its three-dimensional structure was arranged as a β-trefoil fold, which is the typical architecture of "Ricin B chain (or R)-type" lectins, which are found in a broad range of organisms. To date, MytiLec-1-like lectins have been exclusively found in a few species of the mollusk family Mytilidae (M. galloprovincialis, M. trossulus, M. californianus, and Crenomytilus grayanus) and in the phylum Brachiopoda. Transcriptome data revealed the presence of different structural forms of mytilectin in mussels, which included prototype and chimera-type proteins. The primary sequence of these lectins did not match any previously described known protein family, leading to their assignment to the new "mytilectin family." We here report the method of purification of this lectin and describe its use in cell biology.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Phylogeny and anatomy of marine mussels (Bivalvia: Mytilidae) reveal convergent evolution of siphon traits

Convergent morphology is a strong indication of an adaptive trait. Marine mussels (Mytilidae) have long been studied for their ecology and economic importance. However, variation in lifestyle and phenotype also make them suitable models for studies focused on ecomorphological correlation and adaptation. The present study investigates mantle margin diversity and ecological transitions in the Mytilidae to identify macroevolutionary patterns and test for convergent evolution. A fossil-calibrated phylogenetic hypothesis of Mytilidae is inferred based on five genes for 33 species (19 genera). Morphological variation in the mantle margin is examined in 43 preserved species (25 genera) and four focal species are examined for detailed anatomy. Trait evolution is investigated by ancestral state estimation and correlation tests. Our phylogeny recovers two main clades derived from an epifaunal ancestor. Subsequently, different lineages convergently shifted to other lifestyles: semi-infaunal or boring into hard substrate. Such transitions are correlated with the development of long siphons in the posterior mantle region. Two independent origins are reconstructed for the posterior lobules on the inner fold, which are associated with intense mucociliary transport, suggesting an important cleansing role in epifaunal mussels. Our results reveal new examples of convergent morphological evolution associated with lifestyle transitions in marine mussels.

Molecular Diversity of Mytilin-Like Defense Peptides in Mytilidae (Mollusca, Bivalvia)

The CS-αβ architecture is a structural scaffold shared by a high number of small, cationic, cysteine-rich defense peptides, found in nearly all the major branches of the tree of life. Although several CS-αβ peptides involved in innate immune response have been described so far in bivalve mollusks, a clear-cut definition of their molecular diversity is still lacking, leaving the evolutionary relationship among defensins, mytilins, myticins and other structurally similar antimicrobial peptides still unclear. In this study, we performed a comprehensive bioinformatic screening of the genomes and transcriptomes available for marine mussels (Mytilida), redefining the distribution of mytilin-like CS-αβ peptides, which in spite of limited primary sequence similarity maintain in all cases a well-conserved backbone, stabilized by four disulfide bonds. Variations in the size of the alpha-helix and the two antiparallel beta strand region, as well as the positioning of the cysteine residues involved in the formation of the C1-C5 disulfide bond might allow a certain degree of structural flexibility, whose functional implications remain to be investigated. The identification of mytilins in Trichomya and Perna spp. revealed that many additional CS-αβ AMPs remain to be formally described and functionally characterized in Mytilidae, and suggest that a more robust scheme should be used for the future classification of such peptides with respect with their evolutionary origin.

A GM1b/asialo-GM1 oligosaccharide-binding R-type lectin from purplish bifurcate mussels Mytilisepta virgata and its effect on MAP kinases

A 15‐kDa lectin, termed SeviL, was isolated from Mytilisepta virgata (purplish bifurcate mussel). SeviL forms a noncovalent dimer that binds strongly to ganglio‐series GM1b oligosaccharide (Neu5Acɑ2‐3Galβ1‐3GalNAcβ1‐4Galβ1‐4Glc) and its precursor, asialo‐GM1 (Galβ1‐3GalNAcβ1‐4Galβ1‐4Glc). SeviL also interacts weakly with the glycan moiety of SSEA‐4 hexaose (Neu5Acα2‐3Galβ1‐3GalNAcβ1‐3Galα1‐4Galβ1‐4Glc). A partial protein sequence of the lectin was determined by mass spectrometry, and the complete sequence was identified from transcriptomic analysis. SeviL, consisting of 129 amino acids, was classified as an R(icin B)‐type lectin, based on the presence of the QxW motif characteristic of this fold. SeviL mRNA is highly expressed in gills and, in particular, mantle rim tissues. Orthologue sequences were identified in other species of the family Mytilidae, including Mytilus galloprovincialis, from which lectin MytiLec‐1 was isolated and characterized in our previous studies. Thus, mytilid species contain lectins belonging to at least two distinct families (R‐type lectins and mytilectins) that have a common β‐trefoil fold structure but differing glycan‐binding specificities. SeviL displayed notable cytotoxic (apoptotic) effects against various cultured cell lines (human breast, ovarian, and colonic cancer; dog kidney) that possess asialo‐GM1 oligosaccharide at the cell surface. This cytotoxic effect was inhibited by the presence of anti‐asialo‐GM1 oligosaccharide antibodies. With HeLa ovarian cancer cells, SeviL showed dose‐ and time‐dependent activation of kinase MKK3/6, p38 MAPK, and caspase‐3/9. The transduction pathways activated by SeviL via the glycosphingolipid oligosaccharide were triggered apoptosis.

Database

Nucleotide sequence data have been deposited in the GenBank database under accession numbers MK434191, MK434192, MK434193, MK434194, MK434195, MK434196, MK434197, MK434198, MK434199, MK434200, and MK434201.

Organ transcriptomes of the lucinid clam Loripes orbiculatus (Poli, 1791) provide insights into their specialised roles in the biology of a chemosymbiotic bivalve

BACKGROUND

The lucinid clam Loripes orbiculatus lives in a nutritional symbiosis with sulphur-oxidizing bacteria housed in its gills. Although our understanding of the lucinid endosymbiont physiology and metabolism has made significant progress, relatively little is known about how the host regulates the symbiosis at the genetic and molecular levels. We generated transcriptomes from four L. orbiculatus organs (gills, foot, visceral mass, and mantle) for differential expression analyses, to better understand this clam's physiological adaptations to a chemosymbiotic lifestyle, and how it regulates nutritional and immune interactions with its symbionts.

RESULTS

The transcriptome profile of the symbiont-housing gill suggests the regulation of apoptosis and innate immunity are important processes in this organ. We also identified many transcripts encoding ion transporters from the solute carrier family that possibly allow metabolite exchange between host and symbiont. Despite the clam holobiont's clear reliance on chemosynthesis, the clam's visceral mass, which contains the digestive tract, is characterised by enzymes involved in digestion, carbohydrate recognition and metabolism, suggesting that L. orbiculatus has a mixotrophic diet. The foot transcriptome is dominated by the biosynthesis of glycoproteins for the construction of mucus tubes, and receptors that mediate the detection of chemical cues in the environment.

CONCLUSIONS

The transcriptome profiles of gills, mantle, foot and visceral mass provide insights into the molecular basis underlying the functional specialisation of bivalve organs adapted to a chemosymbiotic lifestyle.

Gene expression study in Bathymodiolus azoricus populations from three North Atlantic hydrothermal vent sites

The deep-sea hydrothermal vents are known as harsh environments, abundant in animal diversity surrounded by fluids with specific physiological and chemical composition. Bathymodiolus azoricus mussels are endemic species dwelling at hydrothermal vent sites and at distinct depth ranges. Mussels from Menez Gwen (MG), Lucky Strike (LS), Rainbow (Rb) were collected at 800 m, 1730 m and 2310 m depths respectively, along the Mid-Atlantic Ridge. Five different tissues including gill, digestive gland, mantle, adductor muscle and foot from MG, LS and Rb mussels were selected for gene expression analyses by qPCR. 30 genes were tested to investigate the level of immune and apoptotic gene expression among B. azoricus populations. Statistical analyses confirmed tissue-specific gene expression differences among the five tissues. The digestive gland tissue showed a higher transcriptional activity characterized by an up-regulation of gene activities, contrary to what was assessed in the adductor muscle tissue. Five categories included recognition, signaling, transcription, effector and apoptotic genes were analyzed in this study. The majority of genes differed in levels of expression between MG/LS and LS/Rb in the digestive gland. Our findings suggest that gene expression profiles are inherent to the tissue analyzed, thus implying an immune tissue-specificity controlling defense responses across B. azoricus mussel body as a whole.

De novo assembly, characterization of tissue-specific transcriptomes and identification of immune related genes from the scallop Argopecten purpuratus

The scallop Argopecten purpuratus is one of the most economically important cultured mollusks on the coasts from Chile and Peru but its production has declined, in part, due to the emergence of mass mortality events of unknown origin. Driven by this scenario, increasing progress has been made in recent years in the comprehension of immune response mechanisms in this species. However, it is still not entirely understood how different mucosal interfaces participate and cooperate with the immune competent cells, the hemocytes, in the immune defense. Thus, in this work we aimed to characterize the transcriptome of three tissues with immune relevance from A. purpuratus by next-generation sequencing and de novo transcriptome assembly. For this, 18 cDNA libraries were constructed from digestive gland, gills and hemocytes tissues of scallops from different immune conditions and sequenced by the Illumina HiSeq4000 platform. A total of 967.964.884 raw reads were obtained and 967.432.652 clean reads were generated. The clean reads were de novo assembled into 46.601 high quality contigs and 32.299 (69.31%) contigs were subsequently annotated. In addition, three de novo specific assemblies were performed from clean reads obtained from each tissue cDNA libraries for their comparison. Gene ontology (GO) and KEGG analyses revealed that annotated sequences from digestive gland, gills and hemocytes could be classified into both general and specific subcategory terms and known biological pathways, respectively, according to the tissue nature. Finally, several immune related candidate genes were identified, and the differential expression of tissue-specific genes was established, suggesting they could display specific roles in the host defense. The data presented in this study provide the first insight into the tissue specific transcriptome profiles of A. purpuratus, which should be considered for further research on the interplay between the hemocytes and mucosal immune responses.

In Silico Analysis of Pacific Oyster (Crassostrea gigas) Transcriptome over Developmental Stages Reveals Candidate Genes for Larval Settlement

Following their planktonic phase, the larvae of benthic marine organisms must locate a suitable habitat to settle and metamorphose. For oysters, larval adhesion occurs at the pediveliger stage with the secretion of a proteinaceous bioadhesive produced by the foot, a specialized and ephemeral organ. Oyster bioadhesive is highly resistant to proteomic extraction and is only produced in very low quantities, which explains why it has been very little examined in larvae to date. In silico analysis of nucleic acid databases could help to identify genes of interest implicated in settlement. In this work, the publicly available transcriptome of Pacific oyster Crassostrea gigas over its developmental stages was mined to select genes highly expressed at the pediveliger stage. Our analysis revealed 59 sequences potentially implicated in adhesion of C. gigas larvae. Some related proteins contain conserved domains already described in other bioadhesives. We propose a hypothetic composition of C. gigas bioadhesive in which the protein constituent is probably composed of collagen and the von Willebrand Factor domain could play a role in adhesive cohesion. Genes coding for enzymes implicated in DOPA chemistry were also detected, indicating that this modification is also potentially present in the adhesive of pediveliger larvae.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

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.

Myticalins: A Novel Multigenic Family of Linear, Cationic Antimicrobial Peptides from Marine Mussels (Mytilus spp.)

The application of high-throughput sequencing technologies to non-model organisms has brought new opportunities for the identification of bioactive peptides from genomes and transcriptomes. From this point of view, marine invertebrates represent a potentially rich, yet largely unexplored resource for de novo discovery due to their adaptation to diverse challenging habitats. Bioinformatics analyses of available genomic and transcriptomic data allowed us to identify myticalins, a novel family of antimicrobial peptides (AMPs) from the mussel Mytilus galloprovincialis, and a similar family of AMPs from Modiolus spp., named modiocalins. Their coding sequence encompasses two conserved N-terminal (signal peptide) and C-terminal (propeptide) regions and a hypervariable central cationic region corresponding to the mature peptide. Myticalins are taxonomically restricted to Mytiloida and they can be classified into four subfamilies. These AMPs are subject to considerable interindividual sequence variability and possibly to presence/absence variation. Functional assays performed on selected members of this family indicate a remarkable tissue-specific expression (in gills) and broad spectrum of activity against both Gram-positive and Gram-negative bacteria. Overall, we present the first linear AMPs ever described in marine mussels and confirm the great potential of bioinformatics tools for the de novo discovery of bioactive peptides in non-model organisms.

Myticalins: A Novel Multigenic Family of Linear, Cationic Antimicrobial Peptides from Marine Mussels (Mytilus spp.)

The application of high-throughput sequencing technologies to non-model organisms has brought new opportunities for the identification of bioactive peptides from genomes and transcriptomes. From this point of view, marine invertebrates represent a potentially rich, yet largely unexplored resource for de novo discovery due to their adaptation to diverse challenging habitats. Bioinformatics analyses of available genomic and transcriptomic data allowed us to identify myticalins, a novel family of antimicrobial peptides (AMPs) from the mussel Mytilus galloprovincialis, and a similar family of AMPs from Modiolus spp., named modiocalins. Their coding sequence encompasses two conserved N-terminal (signal peptide) and C-terminal (propeptide) regions and a hypervariable central cationic region corresponding to the mature peptide. Myticalins are taxonomically restricted to Mytiloida and they can be classified into four subfamilies. These AMPs are subject to considerable interindividual sequence variability and possibly to presence/absence variation. Functional assays performed on selected members of this family indicate a remarkable tissue-specific expression (in gills) and broad spectrum of activity against both Gram-positive and Gram-negative bacteria. Overall, we present the first linear AMPs ever described in marine mussels and confirm the great potential of bioinformatics tools for the de novo discovery of bioactive peptides in non-model organisms.