Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica)

Exploration of Molecular Mechanisms of Immunity in the Pacific Oyster (Crassostrea gigas) in Response to Vibrio alginolyticus Invasion

Over the years, oysters have faced recurring mass mortality issues during the summer breeding season, with Vibrio infection emerging as a significant contributing factor. Tubules of gill filaments were confirmed to be in the hematopoietic position in Crassostrea gigas, which produce hemocytes with immune defense capabilities. Additionally, the epithelial cells of oyster gills produce immune effectors to defend against pathogens. In light of this, we performed a transcriptome analysis of gill tissues obtained from C. gigas infected with Vibrio alginolyticus for 12 h and 48 h. Through this analysis, we identified 1024 differentially expressed genes (DEGs) at 12 h post-injection and 1079 DEGs at 48 h post-injection. Enrichment analysis of these DEGs revealed a significant association with immune-related Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. To further investigate the immune response, we constructed a protein-protein interaction (PPI) network using the DEGs enriched in immune-associated KEGG pathways. This network provided insights into the interactions and relationships among these genes, shedding light on the underlying mechanisms of the innate immune defense mechanism in oyster gills. To ensure the accuracy of our findings, we validated 16 key genes using quantitative RT-PCR. Overall, this study represents the first exploration of the innate immune defense mechanism in oyster gills using a PPI network approach. The findings provide valuable insights for future research on oyster pathogen control and the development of oysters with enhanced antimicrobial resistance.

Transcriptional profiling of Bulinus globosus provides insights into immune gene families in snails supporting the transmission of urogenital schistosomiasis

Schistosomiasis, urogenital and intestinal, afflicts 251 million people worldwide with approximately two-thirds of the patients suffering from the urogenital form of the disease. Freshwater snails of the genus Bulinus (Gastropoda: Planorbidae) serve as obligate intermediate hosts for Schistosoma haematobium, the etiologic agent of human urogenital schistosomiasis. These snails also act as vectors for the transmission of schistosomiasis in livestock and wildlife. Despite their crucial role in human and veterinary medicine, our basic understanding at the molecular level of the entire Bulinus genus, which comprises 37 recognized species, is very limited. In this study, we employed Illumina-based RNA sequencing (RNAseq) to profile the genome-wide transcriptome of Bulinus globosus, one of the most important intermediate hosts for S. haematobium in Africa. A total of 179,221 transcripts (N50 = 1,235) were assembled and the benchmarking universal single-copy orthologs (BUSCO) was estimated to be 97.7%. The analysis revealed a substantial number of transcripts encoding evolutionarily conserved immune-related proteins, particularly C-type lectin (CLECT) domain-containing proteins (n = 316), Toll/Interleukin 1-receptor (TIR)-containing proteins (n = 75), and fibrinogen related domain-containing molecules (FReD) (n = 165). Notably, none of the FReDs are fibrinogen-related proteins (FREPs) (immunoglobulin superfamily (IgSF) + fibrinogen (FBG)). This RNAseq-based transcriptional profile provides new insights into immune capabilities of Bulinus snails, helps provide a framework to explain the complex patterns of compatibility between snails and schistosomes, and improves our overall understanding of comparative immunology.

Invertebrate C1q Domain-Containing Proteins: Molecular Structure, Functional Properties and Biomedical Potential

C1q domain-containing proteins (C1qDC proteins) unexpectedly turned out to be widespread molecules among a variety of invertebrates, despite their lack of an integral complement system. Despite the wide distribution in the genomes of various invertebrates, data on the structure and properties of the isolated and characterized C1qDC proteins, which belong to the C1q/TNF superfamily, are sporadic, although they hold great practical potential for the creation of new biotechnologies. This review not only summarizes the current data on the properties of already-isolated or bioengineered C1qDC proteins but also projects further strategies for their study and biomedical application. It has been shown that further broad study of the carbohydrate specificity of the proteins can provide great opportunities, since for many of them only interactions with pathogen-associated molecular patterns (PAMPs) was evaluated and their antimicrobial, antiviral, and fungicidal activities were studied. However, data on the properties of C1qDC proteins, which researchers originally discovered as lectins and therefore studied their fine carbohydrate specificity and antitumor activity, intriguingly show the great potential of this family of proteins for the creation of targeted drug delivery systems, vaccines, and clinical assays for the differential diagnosis of cancer. The ability of invertebrate C1qDC proteins to recognize patterns of aberrant glycosylation of human cell surfaces and interact with mammalian immunoglobulins indicates the great biomedical potential of these molecules.

Understanding Crassostrea virginica tolerance of Perkinsus marinus through global gene expression analysis

Disease tolerance, a host's ability to limit damage from a given parasite burden, is quantified by the relationship between pathogen load and host survival or reproduction. Dermo disease, caused by the protozoan parasite P. marinus, negatively impacts survival in both wild and cultured eastern oyster (C. virginica) populations. Resistance to P. marinus has been the focus of previous studies, but tolerance also has important consequences for disease management in cultured and wild populations. In this study we measured dermo tolerance and evaluated global expression patterns of two sensitive and two tolerant eastern oyster families experimentally challenged with distinct doses of P. marinus (0, 10⁶, 10⁷, and 10⁸ parasite spores per gram wet weight, n = 3-5 individuals per family per dose). Weighted Gene Correlation Network Analysis (WGCNA) identified several modules correlated with increasing parasite dose/infection intensity, as well as phenotype. Modules positively correlated with dose included transcripts and enriched GO terms related to hemocyte activation and cell cycle activity. Additionally, these modules included G-protein coupled receptor, toll-like receptor, and tumor necrosis factor pathways, which are important for immune effector molecule and apoptosis activation. Increased metabolic activity was also positively correlated with treatment. The module negatively correlated with infection intensity was enriched with GO terms associated with normal cellular activity and growth, indicating a trade-off with increased immune response. The module positively correlated with the tolerant phenotype was enriched for transcripts associated with "programmed cell death" and contained a large number of tripartite motif-containing proteins. Differential expression analysis was also performed on the 10⁸ dosed group using the most sensitive family as the comparison reference. Results were consistent with the network analysis, but signals for "programmed cell death" and serine protease inhibitors were stronger in one tolerant family than the other, suggesting that there are multiple avenues for disease tolerance. These results provide new insight for defining dermo response traits and have important implications for applying selective breeding for disease management.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

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.

Integrative Application of Transcriptomics and Metabolomics Provides Insights into Unsynchronized Growth in Sea Cucumber (Stichopus monotuberculatus)

Ever-increasing consumer demand for sea cucumbers mainly leads to huge damage to wild sea cucumber resources, including Stichopus monotuberculatus, which in turn exerts negative impacts on marine environments due to the lack of ecological functions performed by sea cucumbers. Aquaculture of sea cucumbers is an effective way to meet consumer demand and restore their resources. Unsynchronous growth is a prominent problem in the aquaculture of sea cucumbers which has concealed unelucidated molecular mechanisms until now. In this study, we carried out an integrative analysis of transcriptomics and metabolomics on fast-growing (SMF) and slow-growing (SMS) groups of S. monotuberculatus cultured in the same environmental conditions. The results revealed that a total of 2054 significantly differentially expressed genes (DEGs) were identified, which are mainly involved in fat digestion and absorption, histidine metabolism, arachidonic acid metabolism, and glutathione metabolism. 368 differential metabolites (DMs) were screened out between the SMF group and the SMS group; these metabolites are mainly involved in glycerophospholipid metabolism, purine metabolism, biosynthesis of unsaturated fatty acids, pyrimidine metabolism, arachidonic acid metabolism, and other metabolic pathways. The integrative analysis of transcriptomics and metabolomics of S. monotuberculatus suggested that the SMF group had a higher capacity for lipid metabolism and protein synthesis, and had a more frequent occurrence of apoptosis events, which are likely to be related to coping with environmental stresses. The results of this study provide potential values for the aquaculture of sea cucumbers which may promote their resource enhancement.

Transcriptome analysis reveals the diverse response of pearl oyster Pinctada fucata martensii after different PAMP stimulation

Bivalves have evolved effective strategies to combat different pathogens in the environment. They rely on innate immunity to deal with the invasion of various bacteria, viruses, and other microorganisms. However, the molecular mechanisms underlying the responses remain largely unknown. Herein, we constructed 21 transcriptomes of the hemocytes after lipopolysaccharide (LPS), peptidoglycan (PGN) and polyinosinic-polycytidylic acid (poly(I:C)) stimulation to investigate the molecular mechanisms underlying adaptations and plastic responses to different pathogen-related molecular patterns (PAMPs) in pearl oyster Pinctada fucata martensii. Transcriptome analysis revealed 1986-3427 responsive genes enriched in the major immune and cell cycle-related pathways at different times after PAMP stimulation, and the expression patterns of genes under these pathways are complex and diverse. Moreover, "lysosomes" were enriched 6 h after LPS and PGN stimulation, while "peroxisomes" were only enriched in poly(I:C) group. These results suggest different response strategies of pearl oyster to different PAMPs. Furthermore, we identified 261 pattern-recognition receptors (PRRs) including 4 retinoic acid-inducible gene I-like receptors, 38 NOD-like receptors, 83 Toll-like receptors, and 136 C-type lectins in the genome of P. f. martensii. The diverse expression patterns of these PRRs after different PAMP stimulation indicated that pearl oyster evolved complex and specific recognition systems due to tandem repeat and diverse domain combination, which may help pearl oyster cope with the different pathogens in the environment. The present study improved our understanding of the molecular response of pearl oyster to different PAMP stimulation.

Pooled resequencing of larvae and adults reveals genomic variations associated with Ostreid herpesvirus 1 resistance in the Pacific oyster Crassostrea gigas

The Ostreid herpesvirus 1 (OsHV-1) is a lethal pathogen of the Pacific oyster (Crassostrea gigas), an important aquaculture species. To understand the genetic architecture of the defense against the pathogen, we studied genomic variations associated with herpesvirus-caused mortalities by pooled whole-genome resequencing of before and after-mortality larval samples as well as dead and surviving adults from a viral challenge. Analysis of the resequencing data identified 5,271 SNPs and 1,883 genomic regions covering 3,111 genes in larvae, and 18,692 SNPs and 28,314 regions covering 4,863 genes in adults that were significantly associated with herpesvirus-caused mortalities. Only 1,653 of the implicated genes were shared by larvae and adults, suggesting that the antiviral response or resistance in larvae and adults involves different sets of genes or differentiated members of expanded gene families. Combined analyses with previous transcriptomic data from challenge experiments revealed that transcription of many mortality-associated genes was also significantly upregulated by herpesvirus infection confirming their importance in antiviral response. Key immune response genes especially those encoding antiviral receptors such as TLRs and RLRs displayed strong association between variation in regulatory region and herpesvirus-caused mortality, suggesting they may confer resistance through transcriptional modulation. These results point to previously undescribed genetic mechanisms for disease resistance at different developmental stages and provide candidate polymorphisms and genes that are valuable for understanding antiviral immune responses and breeding for herpesvirus resistance.

Evolutionary history of the NLR gene families across lophotrochozoans

NOD-like receptor (NLR) genes are critical innate immune receptors in animals and plants. Lophotrochozoans represent one of the most species-rich superphyla that includes molluscs, segmented worms, flatworms, bryozoans, and other invertebrates, which is crucial to our understanding of immune system evolution in bilaterians. However, NLRs have not been systematically described in lophotrochozoans. We annotated 185 NLRs in 29 lophotrochozoan genomes, and analyzed their domain organization, phylogenetic distribution, molecular evolution, and gene expression. We found that all the 24 molluscan genomes studied encoded no more than three NLRs. None of these molluscan NLRs represented an inducible expression pattern under the infection of eight pathogens; some molluscan NLRs showed developmental stage-specific expression patterns. Instead, 29 molluscan incomplete NLR (incNLR) genes, encoding for proteins absent in the NACHT domain were upregulated under pathogen infection. We also documented the species-specific expansion of NLRs in the clades Polychaeta and Pteriidae. Our study revealed that gene duplication, domain shuffling, gene loss, and novel expression pattern played important roles in the molecular evolution of NLRs.

Hydrocarbon exposure effect on energetic metabolism and immune response in Crassostrea virginica

Crassostrea virginica was exposed to different light crude oil levels to assess the effect on transcriptomic response and metabolic rate. The exposure time was 21 days, and levels of 100 and 200 μg/L were used, including a control. The most significant difference among treatments was the overexpression of several genes associated with energy production, reactive oxygen species (ROS) regulation, immune system response, and inflammatory response. Also, a hydrocarbon concentration-related pattern was identified in ROS regulation, with a gene expression ratio near 1.8:1 between 200 and 100 μg/L treatments. Statistical analysis showed no interaction effect for metabolic rate; however, significant differences were found for oil concentration and time factors, with a higher oxygen consumption at 200 μg/L. Our findings provide novel information about the metabolic response of C. virginica during hydrocarbons exposure. In addition, our results point out which biological processes should be investigated as targets for searching bioindicators.

Transcriptome analysis of the mantle tissue of Pinctada fucata with red and black shells under salinity stress

To understand the molecular responses of Pinctada fucata with different shell colors to salinity stress, we used transcriptome sequencing on the mantle of P. fucata with a black shell and red shell color under the salinity of 20, 35, and 50. The 414 and 2371 differentially expressed genes (DEGs) in P. fucata with a black shell under low- or high-salt stress, while there were 588 and 3009 DEGs in P. fucata with a red shell. KEGG pathway enrichment analysis showed that, under low salt stress, the DEGs of P. fucata with the black shell were significantly enriched in pathways MAPK signaling pathway, protein processing in endoplasmic reticulum, vitamin B6 metabolism, longevity regulating pathway-multiple species, estrogen signaling pathway and antigen processing and presentation, the DEGs of P. fucata with a red shell were significantly enriched in pathways vitamin B6 metabolism. Under high salt stress, the DEGs of P. fucata with a red shell were significantly enriched in pathways arginine biosynthesis. 11 DEGs were randomly selected for quantitative real-time PCR, and the results were consistent with the RNA-seq. In addition, under high salt stress, DEGs were enriched into some pathways related to osmotic regulation and immune defense of P. fucata with black shell and red shell, such as Glycolysis / Gluconeogenesis, AMPK signaling pathway, Beta-Alanine metabolism, Glycine, serine and threonine metabolism, MAPK signaling pathway and Phagosome. The study showed that high salt stress had a greater influence on P. fucata with two shell colors, and P. fucata with a black shell made a positive immune defense response. Our results will improve to further understand the salt tolerance mechanism of P. fucata with different shell colors.

Divergences of the RLR Gene Families across Lophotrochozoans: Domain Grafting, Exon-Intron Structure, Expression, and Positive Selection

Invertebrates do not possess adaptive immunity but have evolved a variety of unique repertoires of innate immune sensors. In this study, we explored the immune diversity and specificity of invertebrates based on the lophotrochozoan RLRs, a major component in antiviral immune recognition. By annotating RLRs in the genomes of 58 representative species across metazoan evolution, we explored the gene expansion of RLRs in Lophotrochozoa. Of note, the N-terminal domains of lophotrochozoan RLRs showed the most striking diversity which evolved independently by domain grafting. Exon–intron structures were revealed to be prevalent in the domain grafting of lophotrochozoan RLRs based on an analysis of sibling paralogs and orthologs. In more than half of the cases, the mechanism of ‘exonization/pseudoexonization’ led to the generation of non-canonical N-terminal domains. Transcriptomic studies revealed that many non-canonical RLRs display immune-related expression patterns. Two of these RLRs showed obvious evidence of positive selection, which may be the result of host defense selection pressure. Overall, our study suggests that the complex and unique domain arrangement of lophotrochozoan RLRs might result from domain grafting, exon–intron divergence, expression diversification, and positive selection, which may have led to functionally distinct lophotrochozoan RLRs.

The expanded inhibitor of apoptosis gene family in oysters possesses novel domain architectures and may play diverse roles in apoptosis following immune challenge

Background

Apoptosis plays important roles in a variety of functions, including immunity and response to environmental stress. The Inhibitor of Apoptosis (IAP) gene family of apoptosis regulators is expanded in molluscs, including eastern, Crassostrea virginica, and Pacific, Crassostrea gigas, oysters. The functional importance of IAP expansion in apoptosis and immunity in oysters remains unknown.

Results

Phylogenetic analysis of IAP genes in 10 molluscs identified lineage specific gene expansion in bivalve species. Greater IAP gene family expansion was observed in C. virginica than C. gigas (69 vs. 40), resulting mainly from tandem duplications. Functional domain analysis of oyster IAP proteins revealed 3 novel Baculoviral IAP Repeat (BIR) domain types and 14 domain architecture types across gene clusters, 4 of which are not present in model organisms. Phylogenetic analysis of bivalve IAPs suggests a complex history of domain loss and gain. Most IAP genes in oysters (76% of C. virginica and 82% of C. gigas), representing all domain architecture types, were expressed in response to immune challenge (Ostreid Herpesvirus OsHV-1, bacterial probionts Phaeobacter inhibens and Bacillus pumilus, several Vibrio spp., pathogenic Aliiroseovarius crassostreae, and protozoan parasite Perkinsus marinus). Patterns of IAP and apoptosis-related differential gene expression differed between the two oyster species, where C. virginica, in general, differentially expressed a unique set of IAP genes in each challenge, while C. gigas differentially expressed an overlapping set of IAP genes across challenges. Apoptosis gene expression patterns clustered mainly by resistance/susceptibility of the oyster host to immune challenge. Weighted Gene Correlation Network Analysis (WGCNA) revealed unique combinations of transcripts for 1 to 12 IAP domain architecture types, including novel types, were significantly co-expressed in response to immune challenge with transcripts in apoptosis-related pathways.

Conclusions

Unprecedented diversity characterized by novel BIR domains and protein domain architectures was observed in oyster IAPs. Complex patterns of gene expression of novel and conserved IAPs in response to a variety of ecologically-relevant immune challenges, combined with evidence of direct co-expression of IAP genes with apoptosis-related transcripts, suggests IAP expansion facilitates complex and nuanced regulation of apoptosis and other immune responses in oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08233-6.

Mantle Transcriptome Provides Insights into Biomineralization and Growth Regulation in the Eastern Oyster (Crassostrea virginica)

Growth of the eastern oyster Crassostrea virginica, a major aquaculture species in the USA, is highly variable and not well understood at molecular levels. As growth of mollusks is confined in shells constructed by the mantle, mantle transcriptomes of large (fast-growing) and small (slow-growing) eastern oysters were sequenced and compared in this study. Transcription was observed for 31,186 genes, among which 104 genes were differentially expressed between the large and small oysters, including 48 upregulated and 56 downregulated in large oysters. Differentially expressed genes (DEGs) included genes from diverse pathways highlighting the complexity of shell formation and growth regulations. Seventeen of the 48 upregulated DEGs were related to shell matrix formation, most of which were upregulated in large oysters, indicating that large oysters are more active in biomineralization and shell formation. Genomic and transcriptomic analyses identified 22 genes encoding novel polyalanine containing proteins (Pacps) with characteristic motifs for matrix function that are tandemly duplicated on one chromosome, all specifically expressed in mantle and at higher levels in large oysters, suggesting that these expanded Pacps play important roles in shell formation and growth. Analysis of sequence variation identified 244,964 SNPs with 328 associated with growth. This study provides novel candidate genes and markers for shell formation and growth, and suggests that genes related to shell formation are important for the complex regulation of growth in the eastern oyster and possibly other bivalve mollusks. Results of this study show that both transcriptional modulation and functional polymorphism are important in determining growth.

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.

The immune response of fairy shrimp Streptocephalus sirindhornae against bacterial black disease by de novo transcriptome analysis

To enhance genomic resources and to understand the molecular immune mechanisms underlying the response of fairy shrimp (Streptocephalus sirindhornae) to pathogens, we first performed a comparative gene transcription analysis from Aeromonas hydrophila-immunized shrimp and from a control group through RNA sequencing. Meanwhile, the differentially expressed genes (DEGs) were investigated, and a total of 46,958,894 clean reads were obtained and then assembled into 73,297 unigenes with an average length of 993 bp and an N50 of 1,458 bp. Unigenes were annotated by comparison with the NR/NT/KO/SwissProt/PFAM/GO and KOG databases, and 28,198 unigenes (38.47%) were annotated in at least one database. After a bacterial challenge, 143 and 287 genes were identified as markedly up- or downregulated, respectively, and 345 were associated with 142 pathways, including the classic immune-related apoptosis, toll-like receptor and MAPK signaling pathways. Moreover, ten differently expressed immune-related genes were confirmed by using quantitative real-time PCR. This study characterized a gene expression pattern for normal and Aeromonas hydrophila-immunized S. sirindhornae for the first time and shed new light on its molecular mechanisms, thus enabling the future efforts of disease control programs for this valuable aquaculture species.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.

Shotgun analysis to identify differences in protein expression between granulocytes and hyalinocytes of the European flat oyster Ostrea edulis

Recovery of wild populations of the European flat oyster Ostrea edulis is important for ecosystem health and conservation of this species, because native oyster populations have dramatically declined or disappeared in most European waters. Diseases have contributed to oyster decline and are important constrains for oyster recovery. Understanding oyster immune system should contribute to design effective strategies to fight oyster diseases. Haemocytes play a pivotal role in mollusc immune responses protecting from infection. Two main types of haemocytes, granulocytes and hyalinocytes, are distinguished in O. edulis. A study aiming to explore differential functions between both haemocyte types and, thus, to enrich the knowledge of Ostrea edulis immune system, was performed by comparing the proteome of the two haemolymph cell types, using a shotgun approach through liquid chromatography (LC) coupled to mass spectrometry (MS). Cells from oyster haemolymph were differentially separated by Percoll density gradient centrifugation. Shotgun LC-MS/MS performance allowed the identification of 145 proteins in hyalinocytes and 138 in the proteome of granulocytes. After a comparative analysis, 55 proteins with main roles in defence were identified, from which 28 were representative of granulocytes and 27 of hyalinocytes, plus 11 proteins shared by both cell types. Different proteins involved in signal transduction, apoptosis, oxidative response, processes related with the cytoskeleton and structure, recognition and wound healing were identified as representatives of each haemocyte type. Important signalling pathways in the immune response such as MAPK, Ras and NF-κβ seemed to be more relevant for granulocytes, while the Wnt signalling pathway, particularly relevant for wound healing, more relevant in hyalinocytes. The differences in proteins involved in recognition and in cytoskeleton and structure suggest differential specialisation in processes of phagocytosis and internalisation of pathogens between haemocyte types. Apoptosis seemed more active in granulocytes. The differences in proteins involved in oxidative response also suggest different redox processes in each cell type.

Phylogenetic analysis of the caspase family in bivalves: implications for programmed cell death, immune response and development

BACKGROUND

Apoptosis is an important process for an organism's innate immune system to respond to pathogens, while also allowing for cell differentiation and other essential life functions. Caspases are one of the key protease enzymes involved in the apoptotic process, however there is currently a very limited understanding of bivalve caspase diversity and function.

RESULTS

In this work, we investigated the presence of caspase homologues using a combination of bioinformatics and phylogenetic analyses. We blasted the Crassostrea gigas genome for caspase homologues and identified 35 potential homologues in the addition to the already cloned 23 bivalve caspases. As such, we present information about the phylogenetic relationship of all identified bivalve caspases in relation to their homology to well-established vertebrate and invertebrate caspases. Our results reveal unexpected novelty and complexity in the bivalve caspase family. Notably, we were unable to identify direct homologues to the initiator caspase-9, a key-caspase in the vertebrate apoptotic pathway, inflammatory caspases (caspase-1, - 4 or - 5) or executioner caspases-3, - 6, - 7. We also explored the fact that bivalves appear to possess several unique homologues to the initiator caspase groups - 2 and - 8. Large expansions of caspase-3 like homologues (caspase-3A-C), caspase-3/7 group and caspase-3/7-like homologues were also identified, suggesting unusual roles of caspases with direct implications for our understanding of immune response in relation to common bivalve diseases. Furthermore, we assessed the gene expression of two initiator (Cg2A, Cg8B) and four executioner caspases (Cg3A, Cg3B, Cg3C, Cg3/7) in C. gigas late-larval development and during metamorphosis, indicating that caspase expression varies across the different developmental stages.

CONCLUSION

Our analysis provides the first overview of caspases across different bivalve species with essential new insights into caspase diversity, knowledge that can be used for further investigations into immune response to pathogens or regulation of developmental processes.

Single-molecule long-read (SMRT) transcriptome sequencing of Mercenaria mercenaria reveals a powerful anti-apoptotic system critical for air exposure endurance

Mercenaria mercenaria is an economically important clam species and exhibits an outstanding resistance to multiple environmental stressors. However, our understanding of their stress adaptability is limited due to a lack of genomic information, such as transcriptome resources. In this study, single-molecule long-read (SMRT) mRNA sequencing was performed to obtain the full-length gill transcriptome reference sequences of M. mercenaria under air exposure stress. In all, 14.5 G subreads were obtained and assembled into 64,603 unigenes, among which 50,613 were successfully annotated. Additionally, 56,295 SSRs, 1457 transcription factors, and 5924 lncRNAs were identified in M. mercenaria transcriptome. Furthermore, numerous apoptosis-related transcripts were identified according to Swiss-Prot annotation and their numbers were counted. We also found that most apoptosis-related transcripts exhibited typical domains of a certain protein family through conserved domain prediction. Additionally, eight typical sequences related to apoptosis pathway were detected by RT-PCR, with the aim to show the sequential variation of gene expression levels under air exposure. These results implied that the complicated apoptosis system, especially the powerful anti-apoptotic system was critical for M. mercenaria to endure air exposure.

Novel globular C1q domain-containing protein (PmC1qDC-1) participates in shell formation and responses to pathogen-associated molecular patterns stimulation in Pinctada fucata martensii

The C1q protein, which contains the globular C1q (gC1q) domain, is involved in the innate immune response, and is found abundantly in the shell, and it participates in the shell formation. In this study, a novel gC1q domain-containing gene was identified from Pinctada fucata martensii (P. f. martensii) and designated as PmC1qDC-1. The full-length sequence of PmC1qDC-1 was 902 bp with a 534 bp open reading frame (ORF), encoding a polypeptide of 177 amino acids. Quantitative real-time PCR (qRT-PCR) result showed that PmC1qDC-1 was widely expressed in all tested tissues, including shell formation-associated tissue and immune-related tissue. PmC1qDC-1 expression was significantly high in the blastula and gastrula and especially among the juvenile stage, which is the most important stage of dissoconch shell formation. PmC1qDC-1 expression was located in the outer epithelial cells of mantle pallial and mantle edge and irregular crystal tablets were observed in the nacre upon knockdown of PmC1qDC-1 expression at mantle pallial. Moreover, the recombined protein PmC1qDC-1 increased the rate of calcium carbonate precipitation. Besides, PmC1qDC-1 expression was significantly up-regulated in the mantle pallial at 6 h and was significantly up-regulated in the mantle edge at 12 h and 24 h after shell notching. The expression level of PmC1qDC-1 in mantle edge was significantly up-regulated at 48 h after LPS stimulation and was significantly up-regulated at 12 h, 24 h and 48 h after poly I:C stimulation. Moreover, PmC1qDC-1 expression was significantly up-regulated in hemocytes at 6 h after lipopolysaccharide (LPS) and poly I:C challenge. These findings suggest that PmC1qDC-1 plays a crucial role both in the shell formation and the innate immune response in pearl oysters, providing new clues for understanding the shell formation and defense mechanism in mollusk.

RNA-seq: The early response of the snail Physella acuta to the digenetic trematode Echinostoma paraensei

To analyze the response of the snail Physella acuta to Echinostoma paraensei, a compatible digenetic trematode, Illumina RNA-seq data were collected from snails with early infection (5 snails at 2 days post-exposure [DPE]) and established infection (4 snails, 8 DPE), and 7 control (unexposed) snails. A reference transcriptome (325,563 transcripts, including 98% of eukaryotic universal single-copy orthologs; BUSCO) and a draft P. acuta genome (employing available genomic Illumina reads; 799,945 scaffolds, includes 88% BUSCO genes) were assembled to guide RNA-seq analyses. Parasite exposure of P. acuta led to 10,195 differentially expressed (DE) genes at 2 DPE and 8,876 DE genes at 8 DPE with only 18% of up-regulated and 22% of down-regulated sequences shared between these time points. Gene ontology (GO) analysis yielded functional annotation of only 1.2% of DE genes but did not indicate major changes in biological activities of P. acuta between 2 and 8 DPE. Increased insights were achieved by analysis of expression profiles of 460 immune-relevant DE transcripts, identified by BLAST and InterProScan. Physella acuta has expanded gene families that encode immune-relevant domains, including CD109/TEP, GTPase IMAP, Limulus agglutination factor (dermatopontin), FReD (≥82 sequences with fibrinogen-related domains), and transcripts that combine C-type lectin (C-LECT) and C1q domains, novel among metazoa. Notably, P. acuta expressed sequences from these immune gene families at all time points, but the assemblages of unique transcripts from particular immune gene families differed between 2 and 8 DPE. The shift in profiles of DE immune genes, from early exposure to parasite establishment, suggests that compatible P. acuta initially respond to infection but switch to express immune genes that likely are less effective against E. paraensei but counter other types of (opportunistic) pathogens and parasites. We propose that the latter expression profile is part of an extended phenotype of E. paraensei, imposed upon P. acuta through parasite manipulation of the host, following successful parasite establishment in the snail after 2 DPE.

Transcriptome Analysis of Juvenile Tilapia (Oreochromis niloticus) Blood, Fed With Different Concentrations of Resveratrol

Oreochromis niloticus (genetically improved farmed tilapia, GIFT) often bites the root of Polygonum cuspidatum when it is used as a floating bed, and resveratrol (RES) is mainly accumulated in the root of P. cuspidatum. Blood acts as a pipeline for the fish immune system. Generating blood transcriptomic resources is crucial for understanding molecular mechanisms underlying blood immune responses. In this study, we determined the effects of RES administration on blood transcriptomic response in GIFT. With increasing RES concentration, 133 (0.025 vs. 0.05 g/kg RES), 155 (0.025 vs. 0.1 g/kg RES), and 123 (0.05 vs. 0.1 g/kg RES) genes were detected as significant differentially expressed genes (DEGs). Three and ninety-five shared significant DEGs were found to be enriched among the three (except 0.1 g/kg RES) and four groups (0, 0.025, 0.05, and 0.1 g/kg RES), respectively. To determine the relationship between mitochondrial regulation and RES supplementation, the results of RNA-Seq were analyzed and nine mitochondria-related genes (ATP synthase or mitochondrial-function-related genes) were verified. The results revealed the same expression pattern: cytochrome c isoform X2 (cox2), katanin p60 ATPase-containing subunit A1 isoform X1 (katna1), plasma membrane calcium-transporting ATPase 1-like (atp2b1) and GTP-binding protein A-like (gtpbpal) showed the highest expression in the 0.1 g/kg RES group, while NADH dehydrogenase [ubiquinone] iron-sulfur protein 2 mitochondrial (nad7), ATP synthase subunit beta, mitochondrial (atpb), ATP synthase subunit alpha, mitochondrial-like (atpal), ATP synthase subunit alpha, mitochondrial (atpa) and ATP-dependent Clp protease proteolytic subunit, mitochondrial (clpp) revealed a dose-dependent expression following RES supplementation. Blood Ca²⁺-ATPase activity, and malondialdehyde, glutathione, and ATP content were significantly increased in the 0.05 (except Ca²⁺-ATPase activity), 0.1 g/kg RES group when compared with the controls. Eighty-nine shared DGEs were mainly enriched in antigen processing and presentation, cell adhesion molecules and phagosome pathways, based on the comparison between previous reported hepatic and the present blood transcriptome. Our study demonstrated that RES supplementation might improve the resistance to metabolism dysfunction via mitochondrial energy synthesis and/or the respiratory chain (e.g., ATPase).

Transcriptomic analysis of Manila clam Ruditapes philippinarum under lipopolysaccharide challenge provides molecular insights into immune response

The Manila clam, Ruditapes philippinarum, is an economically important shellfish in marine aquaculture. A better understanding of the immune system in R. philippinarum will provide the basis for the development of strategies to mitigate the impact of infectious diseases affecting this species but can also be of relevance for other bivalves of commercial interest. In this study, the transcriptional response of the Manila clam under lipopolysaccharide (LPS) challenge was characterized using RNA sequencing. The transcriptomes of LPS challenged group of clams (LH1, LH2 and LH3), and the PBS control group (CH1, CH2 and CH3), were sequenced with the Illumina HiSeq platform. Compared with the unigene expression profile of the control group, 223 unigenes were up-regulated and 389 unigenes were down-regulated in the LPS challenged group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that signal transduction, defense response, and immune-related pathways such as Chemokine signaling pathway, Complement and coagulation cascades, NOD-like receptor signaling pathway, and Inflammatory mediator regulation of TRP channels in sensory system were the most highly enriched pathways among the genes that were differentially expressed under LPS challenge. This study present understanding of the molecular basis underpinning response to LPS challenge and provides useful information for future work on the molecular mechanism of pathogen resistance and immunity in Manila clam.

The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

De novo transcriptome of the whole-body of the gastropod mollusk Philomycus bilineatus, a pest with medical potential in China

Philomycus bilineatus is a highly common gastropod mollusk pest in China and is also utilized to treat infectious diseases. However, no genomic resources are available for this non-model species. In the present study, the transcriptomic analysis of P. bilineatus was completed. After sequencing using the next generation sequencing technology, 9.11 Gb of clean reads were obtained, which led to the assembly and annotation of 145,523 transcripts and 125,690 unigenes. Unigenes were functionally classified using Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 27,554 unigenes were assigned into 55 GO terms, 13,989 unigenes were differentiated into 26 KOG categories, and 16,368 unigenes were assigned to 229 KEGG pathways. Furthermore, 16,614 simple sequence repeats (SSRs), 38 olfactory genes, and 40 antimicrobial peptide/protein genes were identified. The transcriptome profile of P. bilineatus will provide a valuable genomic resource for further study, will promote the development of new pest management strategies through interference of chemosensory communication, and will support potential medicinal uses of this species.

Sources of Resistance to Fusarium solani and Associated Genomic Regions in Common Bean Diversity Panels

Common bean (Phaseolus vulgaris L.) production worldwide is hampered by Fusarium root rot (FRR), which is caused by Fusarium solani. Screening for FRR resistance on a large scale is notoriously difficult and often yields inconsistent results due to variability within the environment and pathogen biology. A greenhouse screening assay was developed incorporating multiple isolates of F. solani to improve assay reproducibility. The Andean (ADP; n = 270) and Middle American (MDP; n = 280) Diversity Panels were screened in the greenhouse to identify genetic factors associated with FRR resistance. Forty-seven MDP and 34 ADP lines from multiple market classes were identified as resistant to FRR. Greenhouse phenotyping repeatability was confirmed via five control lines. Genome-wide association mapping using ∼200k SNPs was performed on standard phenotyping score 1–9, as well as binary and polynomial transformation of score data. Sixteen and seven significant genomic regions were identified for ADP and MDP, respectively, using all three classes of phenotypic data. Most candidate genes were associated with plant immune/defense mechanisms. For the ADP population, ortholog of glucan synthase-like enzyme, senescence-associated genes, and NAC domain protein, associated with peak genomic region Pv08:0.04–0.18 Mbp, were the most significant candidate genes. For the MDP population, the peak SNPs Pv07:15.29 Mbp and Pv01:51 Mbp mapped within gene models associated with ethylene response factor 1 and MAC/Perforin domain-containing gene respectively. The research provides a basis for bean improvement through the use of resistant genotypes and genomic regions for more durable root rot resistance.

Mechanistic molecular responses of the giant clam Tridacna crocea to Vibrio coralliilyticus challenge

Vibrio coralliilyticus is a pathogen of coral and mollusk, contributing to dramatic losses worldwide. In our study, we found that V. coralliilyticus challenge could directly affect adult Tridacna crocea survival; there were dead individuals appearing at 6 h post infection, and there were 45.56% and 56.78% mortality rates in challenged groups after 36 h of infection. The apoptosis rate of hemocytes was significantly increased by 1.8-fold at 6 h after V. coralliilyticus injection. To shed light on the mechanistic molecular responses of T. crocea to V. coralliilyticus infection, we used transcriptome sequencing analysis and other relevant techniques to analyze T. crocea hemocytes at 0 h, 6 h, 12 h and 24 h after V. coralliilyticus challenge. Our results revealed that the total numbers of unigenes and DEGs were 195651 and 3446, respectively. Additional details were found by KEGG pathway enrichment analysis, where DEGs were significantly enriched in immune-related signaling pathways, such as the TLR signaling pathway, and some were associated with signaling related to apoptosis. Quantitative validation results illustrated that with exposure to V. coralliilyticus, the expression of TLR pathway members, TLR, MyD88, IRAK4, TRAF6, and IкB-α, were significantly upregulated (by 22.9-, 9.6-, 4.0-, 3.6-, and 3.9-fold, respectively) at 6 h. The cytokine-related gene IL-17 exhibited an increase of 6.3-fold and 10.5-fold at 3 h and 6 h, respectively. The apoptosis-related gene IAP1 was dramatically increased by 2.99-fold at 6 h. These results indicate that adult T. crocea could initiate the TLR pathway to resist V. coralliilyticus, which promotes the release of inflammatory factors such as IL-17 and leads to the activation of a series of outcomes, such as apoptosis. The response mechanism is related to the T. crocea immunoreaction stimulated by V. coralliilyticus, providing a theoretical basis for understanding T. crocea immune response mechanisms.

A novel tumor necrosis factor in the Pacific oyster Crassostrea gigas mediates the antibacterial response by triggering the synthesis of lysozyme and nitric oxide

Tumor necrosis factors (TNFs) are a group of multifunctional inflammatory cytokines involved in various pathological and immune processes. Recently, a few primitive TNFs have been characterized from molluscs, which play important roles in modulating cell apoptosis, phagocytosis and production of immune-related enzymes. In the present study, a novel TNF (named as CgTNF-2) with the activity to mediate antibacterial response was identified from the Pacific oyster Crassostrea gigas. The open reading frame of CgTNF-2 was of 783 bp encoding a putative polypeptide of 261 amino acids with a typical TNF domain. The deduced amino acid sequence of CgTNF-2 shared high identity with that of TNFs previously identified from other molluscs, such as 96.1% identity with that in oyster C. hongkongensis, 33.7% identity with that in scallop Mizuhopecten yessoensis and 33.0% identity with CgTNF-1 in oyster C. gigas. There were two distinct TNF branches of vertebrate and invertebrate in the phylogenetic tree, and CgTNF-2 was firstly clustered with TNF-14 from C. hongkongensis, and then clustered with other molluscan TNFs. The mRNA transcripts of CgTNF-2 were widely expressed in various oyster tissues, with the highest expression level in hemocytes. The expression level of CgTNF-2 increased significantly at 6 h (2.45-fold and 6.20-fold, respectively, p < 0.05) after peptidoglycan and lipopolysaccharides treatments, and peaked at 12 h (31.86-fold and 7.90-fold, respectively, p < 0.05). The recombinant protein of CgTNF-2 (rCgTNF-2) inhibited the growth of human alveolar basal epithelial (A549) cells at a concentration of 800 ng/mL. After the oysters received an injection of rCgTNF-2, the serum from those oysters exhibited significantly higher antibacterial activity compared to that from control group, evidenced by inhibiting the growth of Vibrio splendidus. Moreover, the lysozyme activity as well as the contents of nitric oxide in the oyster serum also increased significantly. The above results collectively suggested that CgTNF-2 was a novel member of bivalve TNF-α family, which could prompt the antibacterial activity by inducing the lysozyme activity and the production of nitric oxide in the innate immune response of oyster.

Variation in global transcriptomic response to Perkinsus marinus infection among eastern oyster families highlights potential mechanisms of disease resistance

Dermo disease, caused by the protozoan parasite Perkinsus marinus, negatively impacts wild and cultured Eastern oyster populations, yet our knowledge of the mechanistic bases for parasite pathogenicity and the Eastern oyster's response to it is limited. To better understand host responses to the parasite and identify molecular mechanisms underlying disease-resistance phenotypes, we experimentally challenged two families exhibiting divergent Dermo-resistance phenotypes with the parasite, generated global expression profiles using RNAseq and identified differentially expressed transcripts between control and challenged oysters from each family at multiple time points post-parasite injection. The susceptible and resistant families exhibited strikingly different transcriptomic responses to the parasite over a 28-day time period. The resistant family exhibited a strong, focused, early response to P. marinus infection, where many significantly upregulated transcripts were associated with the biological processes "regulation of proteolysis" and "oxidation-reduction process." P. marinus virulence factors are mainly comprised of proteases that facilitate parasite invasion and weaken host humoral defenses, thus host upregulation of transcripts associated with negative regulation of proteolysis is consistent with a Dermo-resistant phenotype. In contrast, the susceptible family mounted a very weak, disorganized, initial response to the parasite. Few transcripts were differentially expressed between control and injected oysters, and no functional enrichment was detected among them. At the final 28 d time point 2450 differentially expressed transcripts were identified and were associated with either "G-protein coupled receptor activity" (upregulated) or "microtubule-based process" (downregulated). A handful of protease inhibitors were differentially expressed between control and injected susceptible oysters, but this function was not enriched in the susceptible data set. The differential expression patterns observed in this study provide valuable insight into the functional basis of Dermo resistance and suggest that the timing of expression is just as important as the transcripts being expressed.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

Transcriptomic analysis of oligochaete immune responses to myxosporeans infection: Branchiura sowerbyi infected with Myxobolus cultus

The Myxozoa are endoparasites characterized by a two-host life cycle that typically involves invertebrates and vertebrates as definitive and intermediate hosts, respectively. However, little is known about invertebrate-myxosporean interactions, particularly about patterns of host immune defense. We used RNA-sequencing to identify genes that are possibly involved in the immune responses of the oligochaete Branchiura sowerbyi naturally infected with Myxobolus cultus. De novo assembly of the B. sowerbyi transcriptome yielded 119,031 unigenes, with an average length of 896 bp and an N50 length of 1754 bp. Comparative transcriptome analysis revealed 4059 differentially expressed genes (DEGs) between M. cultus-infected and uninfected B. sowerbyi groups, including 3802 upregulated genes and 257 downregulated genes. Among the B. sowerbyi immune factors implicated in the responses to M. cultus infection, DEGs related to lectins, ubiquitin-mediated proteolysis, phagocytosis, oxidative-antioxidative responses, proteases, and protease inhibitors were upregulated. The expression of some immune-related molecules such as calmodulin, heat shock proteins, antimicrobial peptides, lysenin, and serum amyoid A protein were also significantly upregulated. The expression patterns of 14 immune-related DEGs identified by RNA-seq were validated by quantitative real-time polymerase chain reaction. This study is the first attempt to characterize the B. sowerbyi transcriptome and identify immune-related molecules possibly associated with M. cultus infection. It is also the first report of invertebrate host-myxosporean interactions at the transcriptomic level. Our results will facilitate the elucidation of adaptive evolution mechanisms of myxosporean parasites in the definitive host and the genetic basis for differences in resistance of invertebrate hosts of different genotypes to a myxosporean species.

How leachates from wasted cigarette butts influence aquatic life? A case study on freshwater mussel Anodontites trapesiali

There are several reports on the damage smoking causes to human health available in the literature, but little is known about the environmental and biological consequences from inappropriate cigarette butt (CB) disposal in urban and natural environments. The immunotoxic, morphotoxic and mutagenic potential of leachates from cigarette butts (LCB) diluted at environmentally relevant rates (LCB1x: 1.375%; LCB10x: 13.75%) was evaluated in adult representatives of the bivalve species Anodontites trapesialis, which was adopted as model organism. Type II hyalinocytes and granulocytes (phagocytic cells) frequency increased in the hemolymph of subjects exposed to the pollutant for 14 days. Based on this outcome, LCB chemical constituents did not induce immunotoxic effects. The treatments also did not seem to have any impact on the subjects' hemocitary morphometry parameters: diameter, area, perimeter, circularity and nucleus - cytoplasm ratio. However, subjects in groups LCB1x and LCB10x recorded a larger number of hyalinocytes with some nuclear abnormality such as micronucleus, blebbed nucleus, asymmetric constriction nucleus, and nuclear multilobulation and binucleation. The association between these abnormalities and the treatments was confirmed by the Cr, Ni, Pb, Zn, Mn and Na bioaccumulation in tissue samples of the bivalve models exposed to LCB. To the best of our knowledge, this is the first report on LCB mutagenicity in representatives of a freshwater bivalve group. Given the chemical complexity of the addressed pollutants, it is imperative to develop further investigations about the topic.

Transcriptomic response and hydrocarbon accumulation in the eastern oyster (Crassostrea virginica) exposed to crude oil

The adverse effect of crude oil on marine invertebrates is well known. To have a better understanding of its effects on marine invertebrates, Crassostrea virginica was exposed to different concentrations (50, 100 and 200 μg/L) of a mixture of super-light and light crude oil for two weeks, evaluating the transcriptomic response of the digestive gland using RNA-Seq and their accumulation in soft tissues. A total of 33,469,374 reads were assembled, which resulted in 61,356 genome assemblies ('Genes'). Trinotate was used for transcript annotation. At the end of this process, 86,409 transcripts were maintained, comprising a broad set of enzymes from xenobiotics metabolism, oxidative stress, stress and immune responses, and energetic metabolism. The enrichment analysis revealed a change in biological processes and molecular functions, finding from 100 to 200 μg/L. Moreover, the differential gene expression analysis showed a dose-dependent transcriptional response, generally up to 100 μg/L and in some cases up to 200 μg/L, which suggested that oysters' response decreased after 100 μg/L; the analysis of crude oil presence in soft tissues indicated that C. virginica is a suitable candidate for ecotoxicology. Finally, these results should contribute to expanding current genomic resources for C. virginica. Furthermore, they will help to develop new studies in aquatic toxicology focused on knowledge in depth of metabolic pathways, jointly with other approaches (such as proteomics) to allow obtaining a complete idea about the eastern oyster response to crude oil.

Applications of flow cytometry in molluscan immunology: Current status and trends

Flow cytometry (FCM) is routinely used in fundamental and applied research, clinical practice, and clinical trials. In the last three decades, this technique has also become a routine tool used in immunological studies of molluscs to analyse physical and chemical characteristics of haemocytes. Here, we briefly review the current implementation of FCM in the field of molluscan immunology. These applications cover a diverse range of practices from straightforward total cell counts and cell viability to characterize cell subpopulations, and further extend to analyses of DNA content, phagocytosis, oxidative stress and apoptosis. The challenges and prospects of FCM applications in immunological studies of molluscs are also discussed.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Extensive Tandem Duplication Events Drive the Expansion of the C1q-Domain-Containing Gene Family in Bivalves

C1q-domain-containing (C1qDC) proteins are rapidly emerging as key players in the innate immune response of bivalve mollusks. Growing experimental evidence suggests that these highly abundant secretory proteins are involved in the recognition of microbe-associated molecular patterns, serving as lectin-like molecules in the bivalve proto-complement system. While a large amount of functional data concerning the binding specificity of the globular head C1q domain and on the regulation of these molecules in response to infection are quickly accumulating, the genetic mechanisms that have led to the extraordinary lineage-specific expansion of the C1qDC gene family in bivalves are still largely unknown. The analysis of the chromosome-scale genome assembly of the Eastern oyster Crassostrea virginica revealed that the 476 oyster C1qDC genes, far from being uniformly distributed along the genome, are located in large clusters of tandemly duplicated paralogs, mostly found on chromosomes 7 and 8. Our observations point out that the evolutionary process behind the development of a large arsenal of C1qDC lectin-like molecules in marine bivalves is still ongoing and likely based on an unequal crossing over.

Transcriptome analysis of juvenile genetically improved farmed tilapia (Oreochromis niloticus) livers by dietary resveratrol supplementation

Here we used RNA-Seq to explore the transcriptomic response and specific involvement of hepatic mRNA of juvenile Oreochromis niloticus (GIFT) as a result of dietary resveratrol supplementation (0.05 g/kg RES). More than 24,513,018 clean reads were reference genome guided assembly into 23,417 unigenes. 12,596 unigenes (29.64%) were annotated to GO database. There were 5, 179 and 1526 genes significantly differentially expressed genes at 15, 30 and 45 d respectively, and 8 KEGG pathways were enriched associated with this immune response. Hyperemia and compressed hepatic sinusoid, fibrosis of liver cell and abnormal hepatic epidermal cell revealed by H&E and SEM analysis respectively. Genes related with cytokine production (il12rb2, scfr), immune system (ig8l, hlfl, cd226, prf1l), autophagy regulation (atg4b), foxo signaling (ccnb2), steroid hormone biosynthesis (cyp3a40), fatty acid metabolism (scd1), metabolism (cacna1b) have been significantly decreased, while genes associated with such pathways above (leap-2, prdx4, mb, homer1, mif, sat1, cytbc1_8) and the pathway of protein processing in endoplasmic reticulum (cne1, tram1) have been significantly increased. These findings suggested RES activated some immune and biological process-related genes to enhance GIFT's innate immunity. It also suggested high concentration addition or long-time administration may bring negative effect in tilapia liver.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

Varieties of immunity activities and gut contents in tilapia with seasonal changes

We performed 16S rDNA sequencing of tilapia fecal samples to analyze changes in tilapia gut contents after cultivation of the fish in the presence of sandwich-like floating beds of Chinese medicinal herbs (5 and 10% planting-areas; 5% Polygonum cuspidatum). The interactive effects between water quality and blood and hepatic pro- and anti-inflammatory concentrations were also assessed. Our results showed that the water quality (i.e., NO₃^--N, NO₂^--N, TP removal rates) improved, and the abundance of Chloroflexi and Cyanobacteria increased. The abundance of Bacteroidetes, Verrucomicrobia, Saccharibacteria, and Actinobacteria showed both significant seasonal decreases and increases in the presence of P. cuspidatum (increases in August and decreases in July). Fish blood and hepatic IL-10 and IFN-γ levels (together with fish sampled in September) significantly increased in the P. cuspidatum group sampled in August, while those of TNF-α (10% sandwich-like, P. cuspidatum), IL-1β (P. cuspidatum), IL-8 (5% sandwich-like in September, S905S) significantly decreased. Heat shock proteins 60 and 70 levels significantly increased in the P. cuspidatum group, and complement C3 and C4 concentrations significantly increased in S905S. This study demonstrated that enhanced immunity through the regulation of pro- and anti-inflammatory proteins was sustained throughout development until harvest, particularly in fish grown with P. cuspidatum.

Integrated transcriptome analysis of immunological responses in the pearl sac of the triangle sail mussel (Hyriopsis cumingii) after mantle implantation

For pearl culture of bivalve Hyriopsis cumingii, implantation of the sabio may cause nucleus discharge and increased host death rates. We performed a transcriptome analysis of the pearl sac of H. cumingii for 30 days after mantle implantation; 293863 unigenes were obtained, and 27176 unigenes were identified using nr, nt, KO, Swiss-Prot, Pfam, GO, and KOG databases. We detected 4878 differentially expressed genes (DEGs) through pairwise comparisons. We speculated that the physical condition of the recipient mussels returned to normal in about one month; the period was divided into six vital phases (0, 2 h-6 h, 12 h-24 h, 48 h to 7 days, 14 days and 30 days) on the basis of the overall similarities in DEGs. We compared the DEGs between time points and identified key immune-related genes. Our findings provide information on the immunological reactions induced by implantation in pearl mussels.

MATria: a unified centrality algorithm

Background

Computing centrality is a foundational concept in social networking that involves finding the most “central” or important nodes. In some biological networks defining importance is difficult, which then creates challenges in finding an appropriate centrality algorithm.

Results

We instead generalize the results of any k centrality algorithms through our iterative algorithm MATRIA, producing a single ranked and unified set of central nodes. Through tests on three biological networks, we demonstrate evident and balanced correlations with the results of these k algorithms. We also improve its speed through GPU parallelism.

Conclusions

Our results show iteration to be a powerful technique that can eliminate spatial bias among central nodes, increasing the level of agreement between algorithms with various importance definitions. GPU parallelism improves speed and makes iteration a tractable problem for larger networks.

The transcriptomic expression of pattern recognition receptors: Insight into molecular recognition of various invading pathogens in Oyster Crassostrea gigas

Pattern recognition receptors (PRRs) are essential in recognizing specific pathogen-associated molecular patterns (PAMPs) on microbes and triggering responses to eliminate the invading pathogens. Previous genomic studies have revealed a great number of PRR genes in the Pacific oyster Crassostrea gigas, a sessile and filter-feeder marine bivalve belonging to the phylum Mollusca. On the survey of PRRs in the assembly oyster reference genome version 9, a total of 1084 PRRs were identified, which were composed of at least 12 gene families. Some of the gene families were significantly expanded, including C-type lectins (CTLs), fibrinogen-related proteins (FREPs), scavenger receptor cysteine-rich repeat protein (SRCRs), leucine-rich repeat (LRR)-only proteins (LRRops), and especially C1q domain-containing proteins (C1qDCs). The transcriptomic profiles of these abundant PRRs in response to PAMP treatments were investigated by RNA-Seq using the SOLiD EZ BeadTM system. Compared to the control library, there were 6,655, 7,273, 7,593, 6,830, 6687 and 8250 differentially expressed genes in the haemocytes of oysters in response to lipopolysaccharide (LPS) stimulation for 6 h, 12 h and 24 h, and peptidoglycan (PGN), glucan (GLU) and poly I:C (IC) stimulation for 12 h, respectively. After stimulation for 12 h, there were 134, 97, 114 and 159 genes up-regulated in the LPS, PGN, GLU and IC library, respectively. Most of the gene families involved in immune response towards PAMPs were C1qDCs, CTLs and FREPs, while only a few members of LRR and immunoglobin-containing proteins (LRRIGs), retinoic acid-inducible gene I [RIG-I]-like receptors (RLRs) and Toll like receptors (TLRs) were up-regulated. After LPS stimulation, the expression level of 258 non-redundant PRR genes in oyster haemocytes increased significantly with different expression pattern, and most of them were C1qDCs, CTLs, LRRops and FREPs. The transcriptomic analyses indicated that there was a dynamic and orchestrated specific expression regulation of numerous PRR genes in response to pathogen invasion. The expanded PRR gene family members were differentiated with more specific functional responses to certain PAMPs rather than the versatile ones. Based on the different expression pattern during the LPS stimulation, the oyster PRRs could be assigned into three consecutive steps in the response against pathogen invading. All the results would provide useful information for future studies of oyster PRRs and deep insight into the researches on invertebrate innate immunity.

Dual transcriptomic analysis of Ostreid herpesvirus 1 infected Scapharca broughtonii with an emphasis on viral anti-apoptosis activities and host oxidative bursts

The ark shell, Scapharca (Anadara) broughtonii, is an economically important marine shellfish species in Northwestern Pacific. Mass mortalities of ark shell adults related to Ostreid herpesvirus-1 (OsHV-1) infection have occurred frequently since 2012. However, due to the lack of transcriptomic resource of ark shells, the molecular mechanisms underpinning the virus-host interaction remains largely undetermined. In the present study, we resolved the dual transcriptome changes of OsHV-1 infected ark shell with Illumina sequencing. A total of 44 M sequence reads were generated, of which 67,119 reads were mapped to the OsHV-1 genome. De novo assembly of host reads resulted in 276,997 unigenes. 74,529 (26.90%), 47,653 (17.20%) and 19, 611 (7.07%) unigenes were annotated into GO, KOG and KEGG database, respectively. According to RSEM expression values, we identified 2998 differentially expressed genes (DEGs) between control and challenged groups, which included 2065 up-regulated unigenes and 933 down-regulated unigenes. Further analysis of functional pathways indicated that OsHV-1 could inhibit host cell apoptosis mainly by the up-regulation of inhibitor of apoptosis protein (IAP), and thus facilitating its successful replication. While host hemoglobins could induce oxidative burst by suppressing its peroxidase activity, and thus defense against OsHV-1 infection. Although we reported a narrow expression of the OsHV-1 genome compared to Crassostrea gigas infection, we highlighted several common viral genes highly expressed in the two hosts, suggesting an important functional role. This study offers insights into the pathogenesis mechanisms of OsHV-1 infection in bivalve mollusks of the Arcidae family.

Differential lncRNA expression profiles reveal the potential roles of lncRNAs in antiviral immune response of Crassostrea gigas

Long noncoding RNAs (lncRNAs) may play widespread roles in various biological processes. However, systematic profiles of lncRNAs in the biological responses of Pacific Oyster (Crassostrea gigas) to pathogen infection have not yet been demonstrated. Here, we have conducted an exhaustive comparative transcriptome analysis using a bioinformatics approach to exam the functions of lncRNAs response to Ostreid herpesvirus 1μVar (OsHV-1μVar) challenge. In total, 101 differentially expressed lncRNAs (DE-lncRNA) during OsHV-1μVar infections were identified. Compared with differentially expressed mRNAs (DE-mRNA), DE-lncRNAs are shorter in terms of overall length but longer in terms of exon length. These lncRNAs shared similar characteristics with previously reported invertebrate lncRNAs, such as relatively low GC content, low exon number and low sequence conservation, but low expression level were not observed. 20 DE-lncRNAs are typically co-expressed with their neighboring genes annotated as GO terms (GO: 0044237), indicating that these lncRNAs are involved in binding and cellular process functions in cis mode. The weighted gene co-expression network (WGCNA) analysis resulted in 15 modules. The highlighted blue module was specifically demonstrated a co-expression relationship between 14 DE-lncRNAs and 17 immune-related DE-mRNAs (IR-DE-mRNA). Three hub lncRNAs within this module were co-expressed with one hub IR-DE-mRNA involved in fibrinogen-related protein. It was speculated that lncRNAs is extensively involved in oyster antiviral innate immune system. The present study will facilitate subsequently experimental studies to unravel the function of lncRNAs in marine invertebrate response to pathogen infection.

Molecular and cellular characterization of apoptosis in flat oyster a key mechanisms at the heart of host-parasite interactions

Bonamia ostreae has been associated with the decline of flat oyster Ostrea edulis populations in some European countries. This obligatory intracellular parasite persists and multiplies into hemocytes. Previous in vitro experiments showed that apoptosis is activated in hemocytes between 1 h and 4 h of contact with the parasite. The flat oyster uses the apoptosis pathway to defend against B. ostreae. However, the parasite might be also able to modulate this response in order to survive in its host. In order to investigate this hypothesis the apoptotic response of the host was evaluated using flow cytometry, transmission electron microscopy and by measuring the response of genes involved in the apoptotic pathway after 4 h. In parallel, the parasite response was investigated by measuring the expression of B. ostreae genes involved in different biological functions including cell cycle and cell death. Obtained results allow describing molecular apoptotic pathways in O. edulis and confirm that apoptosis is early activated in hemocytes after a contact with B. ostreae. Interestingly, at cellular and molecular levels this process appeared downregulated after 44 h of contact. Concurrently, parasite gene expression appeared reduced suggesting that the parasite could inhibit its own metabolism to escape the immune response.

Identification of SNPs potentially related to immune responses and growth performance in Litopenaeus vannamei by RNA-seq analyses

Litopenaeus vannamei is one of the most important shrimp species for worldwide aquaculture. Despite this, little genomic information is available for this penaeid and other closely related taxonomic crustaceans. Consequently, genes, proteins and their respective polymorphisms are poorly known for these species. In this work, we used the RNA sequencing technology (RNA-seq) in L. vannamei shrimp evaluated for growth performance, and exposed to the White Spot Syndrome Virus (WSSV), in order to investigate the presence of Single Nucleotide Polymorphisms (SNPs) within genes related to innate immunity and growth, both features of great interest for aquaculture activity. We analyzed individuals with higher and lower growth rates; and infected (unhealthy) and non-infected (healthy), after exposure to WSSV. Approximately 7,000 SNPs were detected in the samples evaluated for growth, being 3,186 and 3,978 exclusive for individuals with higher and lower growth rates, respectively. In the animals exposed to WSSV we found about 16,300 unique SNPs, in which 9,338 were specific to non-infected shrimp, and 7,008 were exclusive to individuals infected with WSSV and symptomatic. In total, we describe 4,312 unigenes containing SNPs. About 60% of these unigenes returned GO blastX hits for Biological Process, Molecular Function and Cellular Component ontologies. We identified 512 KEGG unique KOs distributed among 275 pathways, elucidating the majority of metabolism roles related to high protein metabolism, growth and immunity. These polymorphisms are all located in coding regions, and certainly can be applied in further studies involving phenotype expression of complex traits, such as growth and immunity. Overall, the set of variants raised herein enriches the genomic databases available for shrimp, given that SNPs originated from nextgen are still rare for this relevant crustacean group, despite their huge potential of use in genomic selection approaches.

Metabolomics Study of Immune Responses of New Zealand Greenshell™ Mussels (Perna canaliculus) Infected with Pathogenic Vibrio sp

Vibrio coralliilyticus is a bacterial pathogen which can affect a range of marine organisms, such as corals, fish and shellfish, with sometimes devastating consequences. However, little is known about the mechanisms involved in the host-pathogen interaction, especially within molluscan models. We applied gas chromatography-mass spectrometry (GC-MS)-based metabolomics to characterize the physiological responses in haemolymph of New Zealand Greenshell™ mussels (Perna canaliculus) injected with Vibrio sp. DO1 (V. coralliilyticus/neptunius-like isolate). Univariate data analyses of metabolite profiles in Vibrio-exposed mussels revealed significant changes in 22 metabolites at 6 h post-infection, compared to non-exposed mussels. Among them, 10 metabolites were up-regulated, while 12 metabolites were down-regulated in infected mussels. Multivariate analyses showed a clear distinction between infected and non-infected mussels. In addition, secondary pathway analyses indicated perturbations of the host innate immune system following infection, including oxidative stress, inflammation and disruption of the TCA cycle, change in amino acid metabolism and protein synthesis. These findings provide new insights into the pathogenic mechanisms of Vibrio infection of mussels and demonstrate our ability to detect detailed and rapid host responses from haemolymph samples using a metabolomics approach.