De novo RNA-Seq analysis of the venus clam, Cyclina sinensis, and the identification of immune-related genes

In search of the Aplysia immunome: an in silico study

The immune repertoires of mollusks beyond commercially important organisms such as the pacific oyster Crassostrea gigas or vectors for human pathogens like the bloodfluke planorb Biomphalaria glabrata are understudied. Despite being an important model for neural aging and the role of inflammation in neuropathic pain, the immune repertoire of Aplysia californica is poorly understood. Recent discovery of a neurotropic nidovirus in Aplysia has highlighted the need for a better understanding of the Aplysia immunome. To address this gap in the literature, the Aplysia reference genome was mined using InterProScan and OrthoFinder for putative immune genes. The Aplysia genome encodes orthologs of all critical components of the classical Toll-like receptor (TLR) signaling pathway. The presence of many more TLRs and TLR associated adapters than known from vertebrates suggest yet uncharacterized, novel TLR associated signaling pathways. Aplysia also retains many nucleotide receptors and antiviral effectors known to play a key role in viral defense in vertebrates. However, the absence of key antiviral signaling adapters MAVS and STING in the Aplysia genome suggests divergence from vertebrates and bivalves in these pathways. The resulting immune gene set of this in silico study provides a basis for interpretation of future immune studies in this important model organism.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Comparative transcriptomic analysis of surf clams (Paphia undulate) infected with two strains of Vibrio spp. reveals the identity of key immune genes involved in host defense

BACKGROUND

Vibrio spp. is the major infection-producing marine bacteria in commercially important bivalve Paphia undulata. The host resistance is the major determining factor for the development of pathogenesis. To explore defense mechanisms, researchers have focused primarily on the study of differential expression of individual or specific groups of host immune genes during pathogen-challenge.

RESULTS

We compared the expression profile in the surf clams infected with avirulent V. alginolyticus and virulent V. parahaemolyticus to mark the possible molecular mechanisms of pathogenesis. Comparison of the differentially expressed genes between the two groups of Vibrio-infected clams revealed that the number of down-regulate genes in V. parahaemolyticus injected clams (1433) were significantly higher than the other group (169). Based on Gene Ontology classification, a large proportion of these down-regulate genes were found to be associated with cellular and molecular mechanisms for pathogen recognition, and immunity development thereby explaining the low survival rate for the V. parahaemolyticus-treated clams and suggesting a higher virulence of this bacterium towards the surf clams. Quantitative real-time PCR of 24 candidate genes related to immunity involving the JAK-STAT signaling pathway, complementary cascade, cytokine signaling pathway, oxidative stress, phagocytosis and apoptosis down regulated under V. parahaemolyticus infection, indicating compromised host defense. Furthermore, we could demonstrate a central role of JAK-STAT pathway in bacterial clearance. dsRNA mediated depletion of a clam STAT homolog gene results in dramatic increase in the infection by V. alginolyticus, a mildly pathogenic strain under control conditions.

CONCLUSIONS

The difference in gene expression profiles in surf clams treated with two Vibrio species with a differential pathogenicity to P. undulate and downstream molecular analysis could enlighten on the probable molecular mechanisms of the Vibrio pathogenesis and the virulence of V. parahaemolyticus in surf clams, which also benefits to develop new strategies for disease control in surf calm aquaculture.

Two toll-like receptors identified in the mantle of Mytilus coruscus are abundant in haemocytes

Toll-like receptors (TLRs) are a large family of pattern recognition receptors (PRRs) that play a critical role in innate immunity. TLRs are activated when they recognize microbial associated molecular patterns (MAMPs) of bacteria, viruses, or fungus. In the present study, two TLRs were isolated from the mantle of the hard-shelled mussel (Mytilus coruscus) and designated McTLR2 and McTLR3 based on their sequence similarity and phylogenetic clustering with Crassostrea gigas, CgiTLR2 and CgiTLR3, respectively. Quantitative RT-PCR analysis demonstrated that McTLR2 and McTLR3 were constitutively expressed in many tissues but at low abundance.

Evolutionary Origins of Toll-like Receptor Signaling

Toll-like receptors (TLRs) are transmembrane pattern recognition receptors that are best known for their roles in innate immunity for the detection of and defense against microbial pathogens. However, TLRs also have roles in many nonimmune processes, most notably development. TLRs direct both immune and developmental programs by activation of downstream signaling pathways, often by activation of the NF-κB pathway. There are two primary TLR subtypes: 1) TLRs with multiple cysteine clusters in their ectodomain (mccTLRs) and 2) TLRs with a single cysteine cluster in their ectodomain (sccTLRs). For some time, it has been known that TLRs and the biological processes that they control are conserved in organisms from insects to mammals. However, genome and transcriptome sequencing has revealed that many basal metazoans also have TLRs and downstream NF-κB signaling components. In this review, we discuss what is known about the structure, biological function, and downstream signaling pathways of TLRs found in phyla from Porifera through Annelida. From these analyses, we hypothesize that mccTLRs emerged in the phylum Cnidaria, that sccTLRs evolved in the phylum Mollusca, and that TLRs have dual immune and developmental biological functions in organisms as ancient as cnidarians.

Evolutionary and functional analysis of Cyclina sinensis c-Jun AP-1 gene in response to LPS stimulation

The transcription factor activator protein-1 (AP-1) plays an essential and critical role in the regulation of numerous downstream genes involved in various physiological and chemical responses. In this study, we identified a full-length cDNA of the c-Jun AP-1 gene (termed Csc-Jun) from the transcriptome library in Cyclina sinensis. The cDNA contains an 825-bp open reading frame that encodes a 274-amino acid protein sequence, including a characteristic Jun transcription factor domain and a highly conserved basic leucine zipper (bZIP) signature that shares 90% identity to that of Ruditapes philippinarum. Furthermore, a phylogenetic analysis using MrBayes and PhyML software (with Bayesian and maximum likelihood approaches, respectively) revealed that the c-Jun AP-1 family genes might be involved in adapting to various environments in different invertebrates. We implemented the PAML software with the maximum likelihood method to further select and verify the positive selection sites (PSSs) in the Mollusca c-Jun AP-1 genes, and we detected four PSSs located in the Jun transcription factor domain. In addition, a spatial expression analysis showed that the Csc-Jun cDNA transcript was ubiquitously expressed in all of the tested tissues and was strongly expressed in the hepatopancreas and weakly expressed in the tissues of the hemocytes, gill filaments, mantle and adductor muscle. Quantitative real-time PCR showed that the expression profiles of Csc-Jun were significantly upregulated at different times in all of the tested tissues when challenged with lipopolysaccharide (LPS). Furthermore, knockdown of Csc-Jun by RNA interference resulted in a higher mortality of C. sinensis following LPS exposure. Finally, we explored the function of the TLR13-MyD88 signaling pathway in the innate immunity of C. sinensis by RNA interference and immune challenges. The results revealed that the mRNA expression levels of Csc-Jun were all decreased (P < 0.01) in normal and stimulated C. sinensis hemocytes. These data collectively indicated that the c-Jun AP-1 gene might play vital roles in innate immunity and provide new evidence for the evolutionary patterns of innate immune genes in Mollusca.

Comparative and evolutionary analysis of an adapter molecule MyD88 in invertebrate metazoans

The myeloid differentiation factor 88 (MyD88) is an essential adapter in Toll-like receptor (TLR) signalling pathways, with TLR the first pattern-recognition receptor (PRR) that was discovered in Drosophila. In the present study, a MyD88 gene was identified and characterized from a commercially important shellfish, Scapharca subcrenata, including a DEATH domain and TIR domain conserved within other molluscs. Furthermore, comparative genomic evidence revealed that MyD88 was of different lengths and contained quantitative exon and intron regions, which might be involved in specific mechanisms. To further explore the phylogenetic relationships of invertebrate metazoan MyD88, we applied MrBayes and PhyML software to construct phylogenetic trees using Bayesian and maximum likelihood approaches, respectively, which suggested that the MyD88 of Arthropoda was closely related to lower invertebrates, in contrast to morphological taxonomy. Finally, we investigated the evolutionary patterns and location of positive selection sites (PSSs) in the MyD88 gene from Arthropoda, Mollusca and Insecta using PAML software with the maximum likelihood method. The data showed that positive selection sites were detected in these groups, and partial sites were located in the TIR domain but were not found in the DEATH domain. To summarize, in this study, we report on the diversification of MyD88 in invertebrate metazoans, the specific evolutionary position of Arthropoda MyD88, and the positive selection pressures on MyD88 of Arthropoda, Mollusca and Insecta. These results are a valuable contribution to understand and clarify the evolutionary pattern of TLR/MyD88 signalling pathways in invertebrate and vertebrate taxa.

The TLR13-MyD88-NF-κB signalling pathway of Cyclina sinensis plays vital roles in innate immune responses

Toll-like receptors, the best known pattern recognition receptors, play important roles in recognizing non-self molecules and binding pathogen-associated molecular patterns in the innate immune system. In the present research, the cDNA and protein characterization of the TLR signalling pathway genes including IRAK4, TRAK6 and IKKα (named CsIRAK4, CsTRAF6 and CsIKKα, respectively) with the typical motifs from Cyclina sinensis showed significant similarity with their homologues from other shellfish. Furthermore, the mRNA transcripts of these three genes are ubiquitously expressed in all tissues tested and are dominantly expressed in C. sinensis haemocytes (P < 0.05). Moreover, IRAK4, TRAK6 and IKKα cDNA expression levels were all up-regulated after injection with Vibrio anguillarum, Micrococcus luteus and poly I:C (P < 0.01) as shown by quantitative real-time PCR, indicating that they were involved in responding to pathogenic stimulation. We explored the function of the TLR13-MyD88-NF-κB signalling pathway in the innate immune responses of C. sinensis by RNA interference and immune challenges. The results suggested the mRNA expression patterns of CsMyD88, CsIRAK4, CsTRAF6, CsIKKα, CsIκB, CsNF-κB, CsC-LYZ and CsAMP were all down-regulated (P < 0.01) in normal and stimulated C. sinensis haemocytes, revealing the involvement of the TLR13-MyD88-NF-κB signalling pathway in innate immunity by positively adjusting internal signalling factors and immune-related genes. In summary, a TLR13-MyD88-NF-κB signalling pathway exists and plays vital roles in innate immune responses in C. sinensis. These findings collectively lay the foundation for studying the functional characterization of internal signalling factors and establishing a regulatory network for the TLR signalling pathway in molluscs.

Transcriptome analysis of Ruditapes philippinarum hepatopancreas provides insights into immune signaling pathways under Vibrio anguillarum infection

The Manila clam, Ruditapes philippinarum, is one of the most economically important aquatic clams that are harvested on a large scale by the mariculture industry in China. However, increasing reports of bacterial pathogenic diseases have had a negative effect on the aquaculture industry of R. philippinarum. In the present study, the two transcriptome libraries of untreated (termed H) and challenged Vibrio anguillarum (termed HV) hepatopancreas were constructed and sequenced from Manila clam using an Illumina-based paired-end sequencing platform. In total, 75,302,886 and 66,578,976 high-quality clean reads were assembled from 101,080,746 and 99,673,538 raw data points from the two transcriptome libraries described above, respectively. Furthermore, 156,116 unigenes were generated from 210,685 transcripts, with an N50 length of 1125 bp, and from the annotated SwissProt, NR, NT, KO, GO, KOG and KEGG databases. Moreover, a total of 4071 differentially expressed unigenes (HV vs H) were detected, including 903 up-regulated and 3168 down-regulated genes. Among these differentially expressed unigenes, 226 unigenes were annotated using KEGG annotation in 16 immune-related signaling pathways, including Toll-like receptor, NF-kappa B, MAPK, NOD-like receptor, RIG-I-like receptor, and the TNF and chemokine signaling pathways. Finally, 20,341 simple sequence repeats (SSRs) and 214,430 potential single nucleotide polymorphisms (SNPs) were detected from the H and HV transcriptome libraries. In conclusion, these studies identified many candidate immune-related genes and signaling pathways and conducted a comparative analysis of the differentially expressed unigenes from Manila clam hepatopancreas in response to V. anguillarum stimulation. These data laid the foundation for studying the innate immune systems and defense mechanisms in R. philippinarum.

Identification and expression analysis of IκB and NF-κB genes from Cyclina sinensis

With the increasing economic importance of Cyclina sinensis aquaculture, interest in its defense mechanisms against pathogenic infection has grown in recent years. Inhibitor of nuclear factor-kappaB (IκB) and nuclear factor-kappaB (NF-κB) are proteins with central roles in many important physiological and pathological processes, such as innate immune responses. In this study, we identified CsIκB and CsNF-κB genes from a C. sinensis transcriptome library. In healthy adult clams, CsIκB and CsNF-κB genes were widely expressed in various tissues and highly expressed in hemocytes. Further, the expression levels of these genes were significantly increased in hemocytes challenged by Vibrio anguillarum, Micrococcus luteus and poly I:C. Inhibition of CsMyD88 expression by RNAi technology significantly altered the mRNA expression patterns of CsIκB and CsNF-κB as measured using quantitative real-time PCR. These results collectively indicated that the NF-κB signaling pathway, including CsIκB and CsNF-κB genes, might be involved in early innate immune responses and may be regulated by a MyD88-dependent signaling pathway in C. sinensis.

Long Non-Coding RNAs (lncRNAs) of Sea Cucumber: Large-Scale Prediction, Expression Profiling, Non-Coding Network Construction, and lncRNA-microRNA-Gene Interaction Analysis of lncRNAs in Apostichopus japonicus and Holothuria glaberrima During LPS Challenge and Radial Organ Complex Regeneration

Long non-coding RNA (lncRNA) structurally resembles mRNA but cannot be translated into protein. Although the systematic identification and characterization of lncRNAs have been increasingly reported in model species, information concerning non-model species is still lacking. Here, we report the first systematic identification and characterization of lncRNAs in two sea cucumber species: (1) Apostichopus japonicus during lipopolysaccharide (LPS) challenge and in heathy tissues and (2) Holothuria glaberrima during radial organ complex regeneration, using RNA-seq datasets and bioinformatics analysis. We identified A. japonicus and H. glaberrima lncRNAs that were differentially expressed during LPS challenge and radial organ complex regeneration, respectively. Notably, the predicted lncRNA-microRNA-gene trinities revealed that, in addition to targeting protein-coding transcripts, miRNAs might also target lncRNAs, thereby participating in a potential novel layer of regulatory interactions among non-coding RNA classes in echinoderms. Furthermore, the constructed coding-non-coding network implied the potential involvement of lncRNA-gene interactions during the regulation of several important genes (e.g., Toll-like receptor 1 [TLR1] and transglutaminase-1 [TGM1]) in response to LPS challenge and radial organ complex regeneration in sea cucumbers. Overall, this pioneer systematic identification, annotation, and characterization of lncRNAs in echinoderm pave the way for similar studies and future genetic, genomic, and evolutionary research in non-model species.

Multiple I-Type Lysozymes in the Hydrothermal Vent Mussel Bathymodiolus azoricus and Their Role in Symbiotic Plasticity

The aim of this study was first to identify lysozymes paralogs in the deep sea mussel Bathymodiolus azoricus then to measure their relative expression or activity in different tissue or conditions. B. azoricus is a bivalve that lives close to hydrothermal chimney in the Mid-Atlantic Ridge (MAR). They harbour in specialized gill cells two types of endosymbiont (gram-bacteria): sulphide oxidizing bacteria (SOX) and methanotrophic bacteria (MOX). This association is thought to be ruled by specific mechanism or actors of regulation to deal with the presence of symbiont but these mechanisms are still poorly understood. Here, we focused on the implication of lysozyme, a bactericidal enzyme, in this endosymbiosis. The relative expression of Ba-lysozymes paralogs and the global anti-microbial activity, were measured in natural population (Lucky Strike--1700 m, Mid-Atlantic Ridge), and in in situ experimental conditions. B. azoricus individuals were moved away from the hydrothermal fluid to induce a loss of symbiont. Then after 6 days some mussels were brought back to the mussel bed to induce a re-acquisition of symbiotic bacteria. Results show the presence of 6 paralogs in B. azoricus. In absence of symbionts, 3 paralogs are up-regulated while others are not differentially expressed. Moreover the global activity of lysozyme is increasing with the loss of symbiont. All together these results suggest that lysozyme may play a crucial role in symbiont regulation.