Oyster RNA-seq data support the development of Malacoherpesviridae genomics

The family of double-stranded DNA (dsDNA) Malacoherpesviridae includes viruses able to infect marine mollusks and detrimental for worldwide aquaculture production. Due to fast-occurring mortality and a lack of permissive cell lines, the available data on the few known Malacoherpesviridae provide only partial support for the study of molecular virus features, life cycle, and evolutionary history. Following thorough data mining of bivalve and gastropod RNA-seq experiments, we used more than five million Malacoherpesviridae reads to improve the annotation of viral genomes and to characterize viral InDels, nucleotide stretches, and SNPs. Both genome and protein domain analyses confirmed the evolutionary diversification and gene uniqueness of known Malacoherpesviridae. However, the presence of Malacoherpesviridae-like sequences integrated within genomes of phylogenetically distant invertebrates indicates broad diffusion of these viruses and indicates the need for confirmatory investigations. The manifest co-occurrence of OsHV-1 genotype variants in single RNA-seq samples of Crassostrea gigas provide further support for the Malacoherpesviridae diversification. In addition to simple sequence motifs inter-punctuating viral ORFs, recombination-inducing sequences were found to be enriched in the OsHV-1 and AbHV1-AUS genomes. Finally, the highly correlated expression of most viral ORFs in multiple oyster samples is consistent with the burst of viral proteins during the lytic phase

IL-17 signaling components in bivalves: Comparative sequence analysis and involvement in the immune responses

The recent discovery of soluble immune-regulatory molecules in invertebrates takes advantage of the rapid growth of next generation sequencing datasets. Following protein domain searches in the transcriptomes of 31 bivalve spp. and in few available mollusk genomes, we retrieved 59 domains uniquely identifying interleukin 17 (IL-17) and 96 SEFIR domains typical of IL-17 receptors and CIKS/ACT1 proteins acting downstream in the IL-17 signaling pathway. Compared to the Chordata IL-17 family members, we confirm a separate clustering of the bivalve domain sequences and a consistent conservation pattern of amino acid residues. Analysis performed at transcript and genome level allowed us to propose an updated view of the components outlining the IL-17 signaling pathway in Mytilus galloprovincialis and Crassostrea gigas (in both species, homology modeling reduced the variety of IL-17 domains to only two 3D structures). Digital expression analysis indicated more heterogeneous expression levels for the mussel and oyster IL-17 ligands than for IL-17 receptors and CIKS/CIKSL proteins. Besides, new qPCR analyses confirmed such gene expression trends in hemocytes and gills of mussels challenged with heat-killed bacteria. These results uphold the involvement of an ancient IL-17 signaling pathway in the bivalve immune responses and, likewise in humans, suggest the possibility of distinctive modulatory roles of individual IL-17s/IL-17 receptors. Overall, the common evidence of pro-inflammatory cytokines and inter-related intracellular signaling pathways in bivalves definitely adds complexity to the invertebrate immunity

Expansion and loss events characterized the occurrence of MIF-like genes in bivalves

Macrophage migration inhibitory factor (MIF) dynamically connects innate and adaptive immune systems in vertebrate animals, allowing highly orchestrated systemic responses to various insults. The occurrence of MIF-like genes in non-vertebrate organisms suggests its origin from an ancestral metazoan gene, whose function is still a matter of debate. In the present work, by analyzing available genomic and transcriptomic data from bivalve mollusks, we identified 137 MIF-like sequences, which were classified into three types, based on phylogeny and conservation of key residues: MIF, D-DT, and the lineage-specific type MDL. Comparative genomics revealed syntenic conservation of homologous genes at the family level, the loss of D-DT in the Ostreidae family as well as the expansion of MIF-like genes in the Mytilidae family, possibly underpinning the neofunctionalization of duplicated gene copies. In M. galloprovincialis, MIF and one D-DT were mostly expressed in haemocytes and mantle rim of untreated animals, while D-DT paralogs often showed very limited expression, suggesting an accessory role or their persistence as relict genes

FicD genes in invertebrates: A tale of transposons, pathogenic and integrated viruses

: Many gene families are shared across the tree of life between distantly related species because of horizontal gene transfers (HGTs). However, the frequency of HGTs varies strongly between gene families and biotic realms suggesting differential selection pressures and functional bias. One gene family with a wide distribution are FIC-domain containing enzymes (FicDs). FicDs catalyze AMPylation, a post-translational protein modification consisting in the addition of adenosine monophosphate to accessible residues of target proteins. Beside the well-known conservation of FicDs in deuterostomes, we report the presence of a conserved FicD gene ortholog in a large number of protostomes and microbial eukaryotes. We also reported additional FicD gene copies in the genomes of some rotifers, parasitic worms and bivalves. A few dsDNA viruses of these invertebrates, including White spot syndrome virus, Cherax quadricarinatus iridovirus, Ostreid herpesvirus-1 and the beetle nudivirus, carry copies of FicDs, with phylogenetic analysis suggesting a common origin of these FicD copies and the duplicated FicDs of their invertebrate hosts. HGTs and gene duplications possibly mediated by endogenous viruses or genetic mobile elements seem to have contributed to the transfer of AMPylation ability from bacteria and eukaryotes to pathogenic viruses, where this pathway could have been hijacked to promote viral infection

DNA Damage and Transcriptional Changes in the Gills of Mytilus galloprovincialis Exposed to Nanomolar Doses of Combined Metal Salts (Cd, Cu, Hg)

[ENG]Aiming at an integrated and mechanistic view of the early biological effects of selected metals in the marine sentinel organism Mytilus galloprovincialis, we exposed mussels for 48 hours to 50, 100 and 200 nM solutions of equimolar Cd, Cu and Hg salts and measured cytological and molecular biomarkers in parallel. Focusing on the mussel gills, first target of toxic water contaminants and actively proliferating tissue, we detected significant dose-related increases of cells with micronuclei and other nuclear abnormalities in the treated mussels, with differences in the bioconcentration of the three metals determined in the mussel flesh by atomic absorption spectrometry. Gene expression profiles, determined in the same individual gills in parallel, revealed some transcriptional changes at the 50 nM dose, and substantial increases of differentially expressed genes at the 100 and 200 nM doses, with roughly similar amounts of up- and down-regulated genes. The functional annotation of gill transcripts with consistent expression trends and significantly altered at least in one dose point disclosed the complexity of the induced cell response. The most evident transcriptional changes concerned protein synthesis and turnover, ion homeostasis, cell cycle regulation and apoptosis, and intracellular trafficking (transcript sequences denoting heat shock proteins, metal binding thioneins, sequestosome 1 and proteasome subunits, and GADD45 exemplify up-regulated genes while transcript sequences denoting actin, tubulins and the apoptosis inhibitor 1 exemplify down-regulated genes). Overall, nanomolar doses of co-occurring free metal ions have induced significant structural and functional changes in the mussel gills: the intensity of response to the stimulus measured in laboratory supports the additional validation of molecular markers of metal exposure to be used in Mussel Watch programsWork granted by MIUR and Co.Ri.La. to PV. LV is currently supported by the FP7-KBBE-2010-4-266157 Bivalife project. Work in the laboratory of MPC is funded by grants to consolidated research groups (ref GIC07/26-IT-393-07, Basque Government), and to the unit of formation and research (UFI ref 11/37, University of the Basque Country) and projects Nanoretox (ref CP-FP 214478-2, UE 7th FP) and Nanocancer (ref CTM2009-13477, Spanish Ministry of Science and Innovation

An Evolutionary Perspective of Dopachrome Tautomerase Enzymes in Metazoans

Melanin plays a pivotal role in the cellular processes of several metazoans. The final step of the enzymically-regulated melanin biogenesis is the conversion of dopachrome into dihydroxyindoles, a reaction catalyzed by a class of enzymes called dopachrome tautomerases. We traced dopachrome tautomerase (DCT) and dopachrome converting enzyme (DCE) genes throughout metazoans and we could show that only one class is present in most of the phyla. While DCTs are typically found in deuterostomes, DCEs are present in several protostome phyla, including arthropods and mollusks. The respective DCEs belong to the yellow gene family, previously reported to be taxonomically restricted to insects, bacteria and fungi. Mining genomic and transcriptomic data of metazoans, we updated the distribution of DCE/yellow genes, demonstrating their presence and active expression in most of the lophotrochozoan phyla as well as in copepods (Crustacea). We have traced one intronless DCE/yellow gene through most of the analyzed lophotrochozoan genomes and we could show that it was subjected to genomic diversification in some species, while it is conserved in other species. DCE/yellow was expressed in most phyla, although it showed tissue specific expression patterns. In the parasitic copepod Mytilicola intestinalis DCE/yellow even belonged to the 100 most expressed genes. Both tissue specificity and high expression suggests that diverse functions of this gene family also evolved in other phyla apart from insects

Tracing the invertebrate herpesviruses in the global sequence datasets

The family of Malacoherpesviridae is currently represented by only two viruses infecting molluscs, Ostreid herpesvirus 1 (OsHV-1) and Haliotid herpesvirus 1 (HaHV-1), both causing detrimental infections in aquaculture species. Malacoherpesvirus-like sequences were also detected through genome sequencing projects in amphioxus (Branchiostoma species) and annelid worm (Capitella teleta), suggesting the existence of a hidden diversity of malacoherpesviruses in aquatic animals. Here, to extend the knowledge on malacoherpesvirus diversity, we searched for the presence of malacoherpesvirus relatives in genomic, transcriptomic and metagenomic datasets, including from the Tara Oceans expedition, and report 4 novel malacoherpesvirus-like genomes (MalacoHV1-4). Genomic analysis suggested gastropods and bivalves as the most probable hosts for these new malacoherpesviruses. Phylogenetic analysis based on the family B DNA polymerase placed the novel MalacoHV1 and MalacoHV3 as sister lineages of OsHV-1 and HaHV-1, respectively, whereas MalacoHV2 and MalacoHV4 showed higher divergence. The viral genome found associated with amphioxus together with MalacoHV4 formed a sister clade to the mollusc and annelid malacoherpesviruses, suggesting an early divergence of the two virus assemblages. In conclusion, although relatively rare in the available sequence databases, the previously undescribed malacoherpesviruses, MalacoHV1-4, circulate in aquatic ecosystems and should be considered as possible emerging viruses under changing environmental conditions

Chromosomal-level assembly of the blood clam, Scapharca (Anadara) broughtonii, using long sequence reads and Hi-C

Abstract Background The blood clam, Scapharca (Anadara) broughtonii, is an economically and ecologically important marine bivalve of the family Arcidae. Efforts to study their population genetics, breeding, cultivation, and stock enrichment have been somewhat hindered by the lack of a reference genome. Herein, we report the complete genome sequence of S. broughtonii, a first reference genome of the family Arcidae. Findings A total of 75.79 Gb clean data were generated with the Pacific Biosciences and Oxford Nanopore platforms, which represented approximately 86× coverage of the S. broughtonii genome. De novo assembly of these long reads resulted in an 884.5-Mb genome, with a contig N50 of 1.80 Mb and scaffold N50 of 45.00 Mb. Genome Hi-C scaffolding resulted in 19 chromosomes containing 99.35% of bases in the assembled genome. Genome annotation revealed that nearly half of the genome (46.1%) is composed of repeated sequences, while 24,045 protein-coding genes were predicted and 84.7% of them were annotated. Conclusions We report here a chromosomal-level assembly of the S. broughtonii genome based on long-read sequencing and Hi-C scaffolding. The genomic data can serve as a reference for the family Arcidae and will provide a valuable resource for the scientific community and aquaculture sector

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 μVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses

Insights into the innate immunity of the Mediterranean mussel Mytilus galloprovincialis

<p>Abstract</p> <p>Background</p> <p>Sessile bivalves of the genus <it>Mytilus </it>are suspension feeders relatively tolerant to a wide range of environmental changes, used as sentinels in ecotoxicological investigations and marketed worldwide as seafood. Mortality events caused by infective agents and parasites apparently occur less in mussels than in other bivalves but the molecular basis of such evidence is unknown. The arrangement of Mytibase, interactive catalogue of 7,112 transcripts of <it>M. galloprovincialis</it>, offered us the opportunity to look for gene sequences relevant to the host defences, in particular the innate immunity related genes.</p> <p>Results</p> <p>We have explored and described the Mytibase sequence clusters and singletons having a putative role in recognition, intracellular signalling, and neutralization of potential pathogens in <it>M. galloprovincialis</it>. Automatically assisted searches of protein signatures and manually cured sequence analysis confirmed the molecular diversity of recognition/effector molecules such as the antimicrobial peptides and many carbohydrate binding proteins. Molecular motifs identifying complement C1q, C-type lectins and fibrinogen-like transcripts emerged as the most abundant in the Mytibase collection whereas, conversely, sequence motifs denoting the regulatory cytokine MIF and cytokine-related transcripts represent singular and unexpected findings. Using a cross-search strategy, 1,820 putatively immune-related sequences were selected to design oligonucleotide probes and define a species-specific Immunochip (DNA microarray). The Immunochip performance was tested with hemolymph RNAs from mussels injected with <it>Vibrio splendidus </it>at 3 and 48 hours post-treatment. A total of 143 and 262 differentially expressed genes exemplify the early and late hemocyte response of the <it>Vibrio</it>-challenged mussels, respectively, with AMP trends confirmed by qPCR and clear modulation of interrelated signalling pathways.</p> <p>Conclusions</p> <p>The Mytibase collection is rich in gene transcripts modulated in response to antigenic stimuli and represents an interesting window for looking at the mussel immunome (transcriptomes mediating the mussel response to non-self or abnormal antigens). On this basis, we have defined a new microarray platform, a mussel Immunochip, as a flexible tool for the experimental validation of immune-candidate sequences, and tested its performance on <it>Vibrio</it>-activated mussel hemocytes. The microarray platform and related expression data can be regarded as a step forward in the study of the adaptive response of the <it>Mytilus </it>species to an evolving microbial world.</p