Clam focal and systemic immune responses to QPX infection revealed by RNA-seq technology

A Comparative Study on Effects of Three Butyric Acid-Producing Additives on the Growth Performance, Non-specific Immunity, and Intestinal Microbiota of the Sea Cucumber Apostichopus japonicus

The providers of butyric acid, Clostridium butyricum (CB), sodium butyrate (SB), and tributyrin (TB), have been extensively studied as aquafeed additives in recent years. However, no comparative study has been reported on the probiotic effects of CB, SB, and TB as feed additives on sea cucumber (Apostichopus japonicus). A 63-day feeding trial was performed to assess the effects of dietary live cells of C. butyricum (CB group, the basal diet supplemented with 1% CB), sodium butyrate (SB group, the basal diet supplemented with 1% SB), and tributyrin (TB group, the basal diet supplemented with 1% TB) on the growth, non-specific immunity, and intestinal microbiota of A. japonicus with a basal diet group as the control. Results indicated that all three additives considerably increased A. japonicus growth, with dietary CB having the optimal growth-promoting effect. Of the seven non-specific enzyme parameters measured in coelomocytes of A. japonicus (i.e., the activities of phagocytosis, respiratory burst, superoxide dismutase, alkaline phosphatase, acid phosphatase, catalase, and lysozyme), dietary CB, SB, and TB considerably increased the activities of six, five, and six of them, respectively. The immune genes (Aj-p105, Aj-p50, Aj-rel, and Aj-lys) expression in the mid-intestine tissues of A. japonicus was significantly increased by all three additives. The CB group had the highest expression of all four genes. Additionally, the relative expression of Aj-p105, Aj-p50, and Aj-lys genes was significantly up-regulated in the three additive groups after stimulation with inactivated Vibrio splendidus. Dietary CB enhanced the intestinal microbial diversity and richness in A. japonicus while dietary TB decreased them. Meanwhile, dietary CB, SB, and TB significantly enhanced the abundance of Firmicutes, unclassified_f_Rhodobacteraceae, and Proteobacteria, respectively, while dietary CB and SB reduced the abundance of Vibrio. Dietary CB and SB improved the stability of microbial ecosystem in the intestine of A. japonicus. In contrast, dietary TB appeared to have a negative effect on the stability of intestinal microbial ecosystem. All three additives improved the intestinal microbial functions associated with energy production and immunity regulation pathways, which may contribute directly to growth promotion and non-specific immunity enhancement in A. japonicus. Collectively, in terms of enhancing growth and non-specific immunity, as well as improving intestinal microbiota, dietary live cells of C. butyricum exhibited the most effective effects in A. japonicus.

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.

The evolution of the metazoan Toll receptor family and its expression during protostome development

BACKGROUND

Toll-like receptors (TLRs) play a crucial role in immunity and development. They contain leucine-rich repeat domains, one transmembrane domain, and one Toll/IL-1 receptor domain. TLRs have been classified into V-type/scc and P-type/mcc TLRs, based on differences in the leucine-rich repeat domain region. Although TLRs are widespread in animals, detailed phylogenetic studies of this gene family are lacking. Here we aim to uncover TLR evolution by conducting a survey and a phylogenetic analysis in species across Bilateria. To discriminate between their role in development and immunity we furthermore analyzed stage-specific transcriptomes of the ecdysozoans Priapulus caudatus and Hypsibius exemplaris, and the spiralians Crassostrea gigas and Terebratalia transversa.

RESULTS

We detected a low number of TLRs in ecdysozoan species, and multiple independent radiations within the Spiralia. V-type/scc and P-type/mcc type-receptors are present in cnidarians, protostomes and deuterostomes, and therefore they emerged early in TLR evolution, followed by a loss in xenacoelomorphs. Our phylogenetic analysis shows that TLRs cluster into three major clades: clade α is present in cnidarians, ecdysozoans, and spiralians; clade β in deuterostomes, ecdysozoans, and spiralians; and clade γ is only found in spiralians. Our stage-specific transcriptome and in situ hybridization analyses show that TLRs are expressed during development in all species analyzed, which indicates a broad role of TLRs during animal development.

CONCLUSIONS

Our findings suggest that a clade α TLR gene (TLR-Ca) and a clade β/γ TLR gene (TLR-Cβ/γ) were already present in the cnidarian-bilaterian common ancestor. However, although TLR-Ca was conserved in cnidarians, TLR-Cβ/γ was lost during the early evolution of these taxa. Moreover, TLR-Cβ/γ duplicated to generate TLR-Cβ and TLR-Cγ in the lineage to the last common protostome-deuterostome ancestor. TLR-Ca, TLR-Cβ and TLR-Cγ further expanded generating the three major TLR clades. While all three clades radiated in several spiralian lineages, specific TLRs clades have been presumably lost in other lineages. Furthermore, the expression of the majority of these genes during protostome ontogeny suggests a likely role in development.

Aplysia Neurons as a Model of Alzheimer's Disease: Shared Genes and Differential Expression

Although Alzheimer’s disease (AD) is the most common form of dementia in the United States, development of therapeutics has proven difficult. Invertebrate alternatives to current mammalian AD models have been successfully employed to study the etiology of the molecular hallmarks of AD. The marine snail Aplysia californica offers a unique and underutilized system in which to study the physiological, behavioral, and molecular impacts of AD. Mapping of the Aplysia proteome to humans and cross-referencing with two databases of genes of interest in AD research identified 898 potential orthologs of interest in Aplysia. Included among these orthologs were alpha, beta and gamma secretases, amyloid-beta, and tau. Comparison of age-associated differential expression in Aplysia sensory neurons with that of late-onset AD in the frontal lobe identified 59 ortholog with concordant differential expression across data sets. The 21 concordantly upregulated genes suggested increased cellular stress and protein dyshomeostasis. The 47 concordantly downregulated genes included important components of diverse neuronal processes, including energy metabolism, mitochondrial homeostasis, synaptic signaling, Ca⁺⁺ regulation, and cellular cargo transport. Compromised functions in these processes are known hallmarks of both human aging and AD, the ramifications of which are suggested to underpin cognitive declines in aging and neurodegenerative disease.

Silver nanoparticles' impact on the gene expression of the cytosolic adaptor MyD-88 and the interferon regulatory factor IRF in the gills and digestive gland of mytilus galloprovincialis

Silver nanoparticles (AgNPs) have been reported as stressors for the bivalves' immune system at different regulatory levels, impacting the detection step and receptors, and other mediators, as well as effector molecules. However, studies on how AgNPs impact the transmission of signals from receptors and whether they have an effect on mediators and transcription factors are still scarce. This study aims to investigate the effect of 12 hours of in vivo exposure to 100 µg/L of AgNPs on the gene expression of the cytosolic adaptor Myeloid, the differentiation protein 88 (MgMyD88-b), and the interferon regulatory factor (Me4-IRF) in the gills and digestive gland of Mytilus galloprovincialis, before and after blocking two major uptake pathways of nanoparticles (clathrin- and caveolae-mediated endocytosis). The results illustrate a tissue-specific gene expression of the MgMyD88-b and the Me4-IRF in the gills and digestive gland of M. galloprovincialis. In the gills, AgNPs did not significantly impact the expression of the two genes. However, blocking the caveolae-mediated endocytosis decreased the expression of Me4-IRF. However, inhibition of clathrin-mediated endocytosis in the digestive gland recorded a significant decrease in the expression of MgMyD88-b. Overall, the inhibition of the AgNPs' uptake routes have highlighted their potential interference with the immune response through the studied mediators' genes, which need to be studied further in future investigations.

Identification of variants associated with hard clam, Mercenaria mercenaria, resistance to Quahog Parasite Unknown disease

Severe losses in aquacultured and wild hard clam (Mercenaria mercenaria) stocks have been previously reported in the northeastern United States due to a protistan parasite called QPX (Quahog Parasite Unknown). Previous work demonstrated that clam resistance to QPX is under genetic control. This study identifies single nucleotide polymorphism (SNP) associated with clam survivorship from two geographically segregated populations, both deployed in an enzootic site. The analysis contrasted samples collected before and after undergoing QPX-related mortalities and relied on a robust draft clam genome assembly. ~200 genes displayed significant variant enrichment at each sampling point in both populations, including 18 genes shared between both populations. Markers from both populations were identified in genes related to apoptosis pathways, protein-protein interaction, receptors, and signaling. This research begins to identify genetic markers associated with clam resistance to QPX disease, leading the way for the development of resistant clam stocks through marker-assisted selection.

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.

Immunotoxicity and transcriptome analysis of zebrafish embryos in response to glufosinate-ammonium exposure

Glufosinate-ammonium (Gla) is a broad-spectrum and non-selective herbicide that widely used in many countries worldwide, but the biological safety including potentially negative effects on aquatic organisms remains largely unknown. In this study, we investigated the immunotoxic effects of Gla exposure on zebrafish embryos. Firstly, Gla markedly decreased the survival rate and caused a series of morphological malformations in a dose-dependent manner. Meanwhile, the number of macrophages and neutrophils was substantially reduced upon Gla exposure. In addition, the levels of oxidative stress were changed and the antioxidant enzyme activities such as CAT and SOD were elevated with the increase of Gla concentrations. Secondly, comparative transcriptome analysis identified 1, 366 differentially expressed genes (DEGs) including 789 up-regulated and 577 down-regulated in zebrafish embryos after Gla exposure. KEGG pathway analysis revealed that metabolic pathways such as drug metabolism-cytochrome P450 was markedly regulated and proteolysis, oxidation-reduction process, and peptidase activity were significantly enriched by the GO analysis. Besides, 55 immunity-related genes were identified in the DEGs, and we found that the genes in the metabolism, redox and immunity display an unique expression profilings by clustering analysis. Finally, 8 inflammatory cytokines and chemokines were further confirmed and they were differentially regulated after Gla exposure. In summary, a global survey of zebrafish defense against glufosinate was performed, and a large number of gene expression levels regarding metabolism, redox, and immunity-related genes were acquired from RNA-Seq. This study provides valuable informations for future elucidating the molecular mechanism of herbicide induced immunotoxicity in aquatic ecosystems.

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.

Long-term affected flat oyster (Ostrea edulis) haemocytes show differential gene expression profiles from naïve oysters in response to Bonamia ostreae

European flat oyster (Ostrea edulis) production has suffered a severe decline due to bonamiosis. The responsible parasite enters in oyster haemocytes, causing an acute inflammatory response frequently leading to death. We used an immune-enriched oligo-microarray to understand the haemocyte response to Bonamia ostreae by comparing expression profiles between naïve (NS) and long-term affected (AS) populations along a time series (1 d, 30 d, 90 d). AS showed a much higher response just after challenge, which might be indicative of selection for resistance. No regulated genes were detected at 30 d in both populations while a notable reactivation was observed at 90 d, suggesting parasite latency during infection. Genes related to extracellular matrix and protease inhibitors, up-regulated in AS, and those related to histones, down-regulated in NS, might play an important role along the infection. Twenty-four candidate genes related to resistance should be further validated for selection programs aimed to control bonamiosis.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.