Transcriptome analysis of Mytilus coruscus hemocytes in response to Vibrio alginnolyficus infection

Exposure to nanopollutants (nZnO) enhances the negative effects of hypoxia and delays recovery of the mussels' immune system

Aquatic environments face escalating challenges from multiple stressors like hypoxia and nanoparticle exposure, with impact of these combined stressors on mussel immunity being poorly understood. We investigated the individual and combined effects of short-term and long-term hypoxia and exposure to zinc oxide nanoparticles (nZnO) on immune system of the mussels (Mytilus edulis). Hemocyte functional traits (mortality, adhesion capacity, phagocytosis, lysosomal abundance, and oxidative burst), and transcript levels of immune-related genes involved in pathogen recognition (the Toll-like receptors, the complement system components, and the adaptor proteins MyD88) were assessed. Short-term hypoxia minimally affected hemocyte parameters, while prolonged exposure led to immunosuppression, impacting hemocyte abundance, viability, phagocytosis, and defensin gene expression. Under normoxia, nZnO stimulated immune responses of mussel hemocytes. However, combined nZnO and hypoxia induced more pronounced and rapid immunosuppression than hypoxia alone, indicating a synergistic interaction. nZnO exposure hindered immune parameter recovery during post-hypoxic reoxygenation, suggesting persistent impact. Opposing trends were observed in pathogen-sensing and pathogen-elimination mechanisms, with a positive correlation between pathogen-recognition system activation and hemocyte mortality. These findings underscore a complex relationship and potential conflict between pathogen-recognition ability, immune function, and cell survival in mussel hemocytes under hypoxia and nanopollutant stress, and emphasize the importance of considering multiple stressors in assessing the vulnerability and adaptability of mussel immune system under complex environmental conditions of anthropogenically modified coastal ecosystems.

Transcriptome analysis reveals tissue-specific responses of Mytilus unguiculatus to Vibrio alginolyticus infection

Mytilus unguiculatus is an important economic bivalve species with wide consumption and aquaculture value. Disease is one of the primary limiting factors in mussel aquaculture, thus understanding the response of different tissues of M. unguiculatus to pathogens is crucial for disease prevention and control. In this study, we investigated the physiological and transcriptomic responses of the gills, adductor muscle, and mantle of M. unguiculatus infected with Vibrio alginolyticus. The results showed that V. alginolyticus infection caused inflammation and tissue structure changes in the gill, adductor muscle and mantle of M. unguiculatus. Meanwhile, the activities of superoxide dismutase and catalase in the three tissues increased, while the total antioxidant capacity decreased, suggesting that M. unguiculatus have an activated defense mechanism against infection-induced oxidative stress, despite a compromised total antioxidant capacity. Transcriptomic studies reveal that infected M. unguiculatus exhibits upregulation of endocytosis, lysosome activity, cellular apoptosis, and immune-related signaling pathways, indicating that M. unguiculatus responds to pathogen invasion by upregulating efferocytosis. Compared with the gill and adductor muscle, the mantle had a higher level of mytimycin, mytilin and myticin, and the three tissues also increased the expression of mytimycin to cope with the invasion of pathogens. In addition, the analysis of genes related to taste transduction pathways and muscle contraction and relaxation found that after infection with V. alginolyticus, M. unguiculatus may reduce appetite by inhibiting taste transduction in the gill, while improving muscle contraction of the adductor muscle and keeping the shell closed, to resist further invasion of pathogens and reduce the risk of pathogen transmission in the population.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Integrative mRNA-miRNA interaction analysis associated with the immune response of Strongylocentrotus intermedius to Vibrio harveyi infection

Strongylocentrotus intermedius is one of the most economically valuable sea urchin species in China and has experienced mass mortality owing to outbreaks of bacterial diseases such as black mouth disease. This has caused serious economic losses to the sea urchin farming industry. To investigate the immune response mechanism of S. intermedius with different tube feet colors in response to Vibrio harveyi infection, we examined the different tube feet-colored S. intermedius under V. harveyi challenge and compared their transcriptome and microRNA (miRNA) profiles using RNA-Seq. We obtained 1813 differentially expressed genes (DEGs), 28 DE miRNAs, and 303 DE miRNA-DEG pairs in different tube feet-colored S. intermedius under V. harveyi challenge. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the most significant DEGs were associated with the Notch signaling and phagosome pathways. The target genes of immune-related miRNAs (miR-71, miR-184, miR-193) and genes (CALM1, SPSB4, DMBT, CSRP1) in S. intermedius were predicted and validated. This study provides insight into the molecular mechanisms that regulate genes involved in the immune response of S. intermedius infected with V. harveyi.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Comparative transcriptomic analysis revealed changes in multiple signaling pathways involved in protein degradation in the digestive gland of Mytilus coruscus during high-temperatures

As a result of global warming, the Mytilus coruscus living attached in the intertidal zone experience extreme and fluctuating changes in temperature, and extreme temperature changes are causing mass mortality of intertidal species. This study explores the transcriptional response of M. coruscus at different temperatures (18 °C, 26 °C, and 33 °C) and different times (0, 12, and 24 h) of action by analyzing the potential temperature of the intertidal zone. In response to high temperatures, several signaling pathways in M. coruscus, ribosome, endocytosis, endoplasmic reticulum stress, protein degradation, and lysosomes, interact to counter the adverse effects of high temperatures on protein homeostasis. Increased expression of key genes, including heat shock proteins (Hsp70, Hsp20, and Hsp110), Lysosome-associated membrane glycoprotein (LAMP), endoplasmic reticulum chaperone (BiP), and baculoviral IAP repeat-containing protein 7 (BIRC7), may further mitigate the effects of heat stress and delay mortality in M. coruscus. These results reveal changes in multiple signaling pathways involved in protein degradation during high-temperature stress, which will contribute to our overall understanding of the molecular mechanisms underlying the response of M. coruscus to high-temperature stress.

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.

Immune responses to ZnO nanoparticles are modulated by season and environmental temperature in the blue mussels Mytilus edulis

Increased production and release of ZnO nanoparticles (nZnO) can cause toxic effects on marine ecosystems and aquatic organisms. However, nZnO toxicity and its modulation by common environmental stressors such as temperature are not yet fully understood. We examined the responses of immune cells (hemocytes) of the blue mussels (Mytilus edulis) exposed to different concentrations (0, 10, 100 μg l^-1) of nZnO or dissolved zinc combined with two temperatures (ambient (10 °C in winter and 15 °C in summer) and warming (+5 °C above ambient temperature)) in winter and summer for 21 days. In winter mussels, exposure to nZnO induced a strong transcriptomic response in multiple immune and inflammation-related genes, stimulated phagocytosis and hemocyte mortality yet suppressed adhesion capacity of hemocytes. In summer mussels, the immune cell responses to nZnO were blunted. The transcriptional responses of hemocytes to dissolved Zn were qualitatively similar but weaker than the responses to nZnO. In the absence of the toxic stress, +5 °C warming lead to dysregulation of the transcription of key immune-related genes in the summer but not the winter mussels. Seasonal warm acclimatization and additional warming in summer suppressed the nZnO-induced transcriptional upregulation of antimicrobial peptides, Toll-like receptors and the complement system. These findings demonstrate that nZnO act as an immunogen in M. edulis and indicate that +5 °C warming might have detrimental effect on innate immunity of the temperate mussel populations in summer when exposure to pathogens is especially high. Capsule: ZnO nanoparticles act as an immunotoxicant inducing a strong immune response in the mussels which is dysregulated by warming in summer but not in winter.

microRNAs facilitate comprehensive responses of Bathymodiolinae mussel against symbiotic and nonsymbiotic bacteria stimulation

Bathymodiolinae mussels are dominant species in cold seeps and hydrothermal vents and could harbor endosymbionts in gill bacteriocytes. However, mechanisms underlying the symbiosis have remained largely undisclosed for years. In the present study, the global expression pattern of immune-related genes and miRNAs were surveyed in Gigantidas platifrons during bacterial challenges using enriched symbiotic methane oxidation bacteria MOBs or nonsymbiotic Vibrio. As a result, multiple pattern recognition receptors were found differentially expressed at 12 h and 24 h post bacteria challenges and distinctly clustered between stimulations. Dozens of immune effectors along with signal transducers were also modulated simultaneously during MOB or Vibrio challenge. A total of 459 miRNAs were identified in the gill while some were differentially expressed post MOB or nonsymbiotic bacteria challenge. A variety of immune-related genes were annotated as target genes of aforesaid differentially expressed miRNAs. As a result, biological processes including the immune recognition, lysosome activity and bacteria engulfment were suggested to be dynamically modulated by miRNAs in either symbiotic or nonsymbiotic bacteria challenge. It was suggested that G. platifrons mussels could maintain a robust immune response against invading pathogens while establishing symbiosis with chemosynthetic bacteria with the orchestra of immune-related genes and miRNAs.

Effects of Ocean Acidification and Microplastics on Microflora Community Composition in the Digestive Tract of the Thick Shell Mussel Mytilus coruscus Through 16S RNA Gene Sequencing

Ocean acidification and microplastic pollution is a global environmental threat, this research evaluated the effects of ocean acidification and microplastics on mussel digestive tract microbial community. The 16S rRNA gene was sequenced to characterize the flora. Species diversity in the samples was assessed by clustering valid tags on 97% similarity. Bacteroidetes, Firmicutes and Proteobacteria were the three most abundant genera in the four groups, with Bacteroidetes showing the highest diversity. However, no differences in flora structure were evident under various treatments. Phylogenetic relationship analysis revealed Bacteroidetes and Firmicutes had the highest OTU diversity. The weighted UniFrac distance, principal coordinate analysis (PCoA), unweighted pair group method with arithmetic mean (UPGMA) cluster tree and analysis of molecular variance (AMOVA) evaluation results for all samples also showed that changes in pH and microplastics concentration did not significantly affect the microbial community structure in the mussel digestive tract. The results presented the no significant effects of ocean acidification and microplastics intake on mussel intestinal diversity.

Integrated analysis of miRNAome and transcriptome reveals miRNA-mRNA network regulation in Vibrio alginolyticus infected thick shell mussel Mytilus coruscus

The thick shell mussel Mytilus coruscus has developed into a model species for studying the interaction between molluscs and environmental stimuli. Herein, integrated analysis of miRNAome and transcriptome was performed to reveal miRNA-mRNA network regulation in Vibrio alginolyticus infected M. coruscus. There have detected some histological abnormalities in digestive gland and gills of V. alginolyticus challenged mussels, ascertaining the effective irritation by the present bacterial strain. A total of 265 novel miRNAs were finally predicted, of which 26 were differentially expressed miRNAs (DEMs). Additionally, 667 differentially expressed genes (DEGs) were detected, which may be potentially associated with innate immune response to V. alginolyticus infection. A regulatory network linked to 22 important pathways and 16 DEMs and 34 OGs was constructed. Some traditional immune-related signaling pathways such as toll-like receptor signaling pathway (TLR) signaling pathway, transforming growth factor-beta (TGF-beta) signaling pathway, peroxisome, phagosome, lysosome, mammalian target of rapamyoin (mTOR) signaling pathway were linked to specific miRNAs and genes in this network. Further, interactional relationship between certain miRNAs and TLR pathway was dissected, which the results predicted that a number of TLRs and TLR-associated signaling genes including TLR1, TLR2, TLR4, TLR6, IRAK1, TRAF6, MAPK, and IL-17 were negatively regulated by novel_miR_11, novel_miR_145, novel_miR_196, novel_miR_5, novel_miR_163 and novel_miR_217 in the TLR pathway. Additionally, interactional relationship between novel_miR_145 and TLR2 was validated by laboratory experiment. The integrated analysis of mRNA and microRNA deep sequencing data exhibited a sophisticated miRNA-mRNA regulation network in M. coruscus in response to V. alginolyticus challenge, which shed a new light on the underlying mechanism of molluscan confronting bacterial infection.

The Nrf2 molecule trigger antioxidant defense against acute benzo(a)pyrene exposure in the thick shell mussel Mytilus coruscus

The NF-E2-related factor 2 (Nrf2), an ubiquitous, evolutionarily conserved transcription factor, acts as a major sensor of oxidative stress in cells. In the present study, a Nrf2 homolog was newly identified in the thick shell mussel Mytilus coruscus. Accordingly, its functional role in antioxidant defense in response to acute benzo(a)pyrene (Bap) exposure was assessed. The newly identified McNrf2 affiliated to traditional Nrf2 family through Blast, multiple alignment and phylogenetic analysis. After acute exposure to Bap, antioxidants including superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPx) and glutathine reductase (GR) were significantly induced in gills and digestive glands at both mRNA and enzymatic levels, and the expression of McNrf2 mRNA was also up-regulated. The analysis of correlating the expression of McNrf2 and the mRNA levels of these antioxidant genes showed positive ties, indicating that Nrf2 was needed for protracted induction of such genes. Further, the recombinant McNrf2 was produced through pET-32a prokaryotic system. After 50 μg/L Bap exposure, ROS generation and LPO level in gills of Nrf2 over-expressed mussels significantly decreased compared to Nrf2 wild-type mussels, as well as reduced ROS production in digestive glands. Collectively, these results show that Nrf2 pathway can provide protection from oxidative stress triggered by Bap in the thick shell mussel.

Transcriptome analysis and immune-related genes expression reveals the immune responses of Macrobrachium rosenbergii infected by Enterobacter cloacae

Macrobrachium rosenbergii is an important cultural species in China and other Southeast Asian countries. However, Enterobacter cloacae infection has caused a great economic loss in M. rosenbergii culture industry. The immune responses of M. rosenbergii to the E. cloacae infection is not fully characterized. To investigate the immune response of M. rosenbergii against E. cloacae, we performed transcriptome analysis of the M. rosenbergii hepatopancreas with and without E. cloacae infection using RNA-seq. After assembly and annotation, 29,731 high quality unigenes were obtained from RNA-seq data. Differential expression analysis revealed the existence of 2498 significantly differently expressed genes (DEGs) at 12 h post infection, with 1365 up-regulated and 1133 down-regulated genes. Among these DEGs, some well-known immune-related genes were up-regulated significantly, including C-type lectin 1, lectin 3, anti-lipopolysaccharide factor 2, Cu/Zn superoxide dismutase and heat shock protein 70. GO analysis demonstrated 24 biological process subcategories, 14 cellular component subcategories, and 12 molecular function subcategories that were enriched among these DEGs, and some DEGs were clustered into immune related subcategories such as immune system process, response to stimulus, biological adhesion, and antioxidant activity. These DEGs were enriched into 216 KEGG pathways including a core set of immune correlated pathways notably in phagosome and lysosome. In addition, 5 up-regulated and 5 down-regulated immune-related DEGs were selected for further validation by quantitative real-time PCR and the results showed consistence with the RNA-seq data. Additionally, the expression level of six selected immune-related genes (ALF2, CLEC1, LEC3, hemocyanin1, HSP70 and SOD) based on the transcriptomic data were monitored at different point of time in hepatopancreas, gill, hemolymph and intestine. Results revealed these immune-related genes were significantly up-regulated in different tissues from 6 to 24 h after E. cloacae infection. Overall, these results provided valuable information for further studying the immune response of M. rosenbergii against E. cloacae infection.

Transcriptome analysis of Catarina scallop (Argopecten ventricosus) juveniles treated with highly-diluted immunomodulatory compounds reveals activation of non-self-recognition system

Marine bivalve hatchery productivity is continuously challenged by apparition and propagation of new diseases, mainly those related to vibriosis. Disinfectants and antibiotics are frequently overused to prevent pathogen presence, generating a potential negative impact on the environment. Recently, the use of highly diluted compounds with immunostimulant properties in marine organisms has been trailed successfully to activate the self-protection mechanisms of marine bivalves. Despite their potential as immunostimulants, little is known about their way of action. To understand their effect, a comparative transcriptomic analysis was performed with Argopecten ventricosus juveniles. The experimental design consisted of four treatments formulated from pathogenic Vibrio lysates at two dilutions: [(T1) Vibrio parahaemolyticus and Vibrio alginolyticus 1D; (T2) V. parahaemolyticus and V. alginolyticus 7C]; minerals [(T3) PhA+SiT 7C], scorpion venom [(T4) ViT 31C]; and one control (C1) hydro-alcoholic solution (ethanol 1%). The RNA sequencing (RNAseq) analysis showed a higher modulation of differentially expressed genes (DEG) in mantle tissue compared to gill tissue. The scallops that showed a higher number of DEG related to immune response in mantle tissue corresponded to T1 (V. parahaemolyticus and V. alginolyticus lysate) and T3 (Silicea terra® - Phosphoric acid®). The transcriptome analysis allowed understanding some interactions between A. ventricosus juveniles and highly-diluted treatments.

Transcriptome analysis and histopathological observations of Geloina erosa gills upon Cr(VI) exposure

The heavy metal contamination like Cr(VI) has been increased by human activities and that threats the ecosystem health of mangrove areas. Bioindicator is an emerging tool in the environmental contamination assessment. The objective of this study was to investigate the Geloina erosa response mechanisms and sensitivities of several biomarkers in the Cr(VI) exposure and identify the G. erosa capability of being used as heavy metals bioindicator. In this study, G. erosa was exposed to 100 μmol·L^-1 Cr(VI) for 48 h. After transcriptome sequencing, a total of 134,817 unigenes were obtained, including 12,555 up-regulated and 18,829 down-regulated differentially expressed genes and were validated through quantitative real-time PCR. In addition, a total of 12,185 SSRs and 1,428,214 candidate SNPs were identified from all the G. erosa transcriptome libraries. Histopathology of the gill indicated the Cr(VI) exposure induced damage of the organ leading to its immunization, detoxification or apoptosis reactions. Among eight genes of the selected biomarkers, Calm, HSP70, CYP450, ATG5, TLR2, MYD88 and CASP8 were up-regulated, while TLR4 was down-regulated in response to the Cr(VI) exposure.

Interactive effects of salinity variation and exposure to ZnO nanoparticles on the innate immune system of a sentinel marine bivalve, Mytilus edulis

ZnO nanoparticles (nZnO) are released into the coastal environment from multiple sources, yet their toxicity to marine organisms is not well understood. We investigated the interactive effects of salinity (normal 15, low 5, and fluctuating 5-15) and nZnO (100 μg l-1) on innate immunity of the blue mussels Mytilus edulis from a brackish area of the Baltic Sea. Exposure to ionic Zn (100 μg l-1) was used to test whether the toxic effects of nZnO can be attributed to the potential release of Zn2+. Functional parameters and the expression of key immune-related genes were investigated in the mussels exposed to nZnO or ionic Zn under different salinity regimes for 21 days. nZnO exposures elevated hemocyte mortality, suppressed adhesion, stimulated phagocytosis, and led to an apparent increase in lysosomal volume. At salinity 15, nZnO suppressed the mRNA expression of the Toll-like receptors TLRb and c, C-lectin, and the complement system component C3q indicating impaired ability for pathogen recognition. In contrast, the mRNA levels of an antimicrobial peptide defensin increased during nZnO exposure at salinity 15. At fluctuating salinity (5-15), nZnO exposure increased expression of multiple immune-related genes in hemocytes including the complement system components C1 and C3q, and the Toll-like receptors TLRa, b and c. Low salinity (5) had strong immunosuppressive effects on the functional and molecular immune traits of M. edulis that overshadowed the effects of nZnO. The salinity-dependent modulation of immune response to nZnO cannot be attributed to the differences in the aggregation or solubility of nZnO, and likely reflects the interaction of the toxic effects of nanoparticles and physiological effects of the osmotic stress. These findings have implications for the environmental risk assessment of nanomaterials and the development of the context-specific biomarker baselines for coastal pollution monitoring.

Genome-wide identification, characterization of RLR genes in Yesso scallop (Patinopecten yessoensis) and functional regulations in responses to ocean acidification

Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), are crucial sensors with a conserved structure in cytoplasm, inducing the production of cytokines, chemokines and host restriction factors which mediate a variety of intracellular activities to interfere with distinct PAMPs (pathogen-associated molecular patterns) for eliminating pathogens in innate immune system. Although RLR genes have been investigated in most vertebrates and some invertebrates, the systematic identification and characterization of RLR genes have not been reported in scallops. In this study, four RLR genes (PY-10413.4, PY-10413.5, PY-443.7 and PY-443.8, designated PyRLRs) were identified in Yesso scallop (Patinopecten yessoensis) through whole-genome scanning through in silico analysis, including two pairs of tandem duplicate genes located on the same scaffold (PY-10413.4 and PY-10413.5, PY-443.7 and PY-443.8, respectively). Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of these genes. The expression profiles of PyRLRs were determined in all developmental stages, in healthy adult tissues, and in mantles that simulated ocean acidification (OA) exposure (pH = 6.5 and 7.5) at different time points (3, 6, 12 and 24 h). Spatiotemporal expression patterns suggested the functional roles of PyRLRs in all stages of development and growth of the scallop. Regulation expressions revealed PY-10413.4 and PY-10413.5 with one or two CARD(s) (caspase activation and recruitment domain) were up-regulated expressed at most time points, whereas PY-443.8 and PY-10413.4 without CARD were significantly down-regulated at each time points, suggesting functional differentiations in the two pairs of PyRLRs based on the structural differences in response to OA. Collectively, this study demonstrated gene duplication of RLR family genes and provide primary analysis for versatile roles in the response of the bivalve innate immune system to OA challenge.

Transcriptomic response to three osmotic stresses in gills of hybrid tilapia (Oreochromis mossambicus female × O. urolepis hornorum male)

Background

Osmotic stress is a widespread phenomenon in aquatic animal. The ability to cope with salinity stress and alkaline stress is quite important for the survival of aquatic species under natural conditions. Tilapia is an important commercial euryhaline fish species. What’s more tilapia is a good experimental material for osmotic stress regulation research, but the molecular regulation mechanism underlying different osmotic pressure of tilapia is still unexplored.

Results

To elucidate the osmoregulation strategy behind its hyper salinity, alkalinity and salinity-alkalinity stress of tilapia, the transcriptomes of gills in hybrid tilapia (Oreochromis mossambicus ♀ × O. urolepis hornorum ♂) under salinity stress (S: 25‰), alkalinity stress(A: 4‰) and salinity-alkalinity stress (SA: S: 15‰, A: 4‰) were sequenced using deep-sequencing platform Illumina/HiSeq-2000 and differential expression genes (DEGs) were identified. A total of 1958, 1472 and 1315 upregulated and 1824, 1940 and 1735 downregulated genes (P-value < 0.05) were identified in the salt stress, alkali stress and saline-alkali stress groups, respectively, compared with those in the control group. Furthermore, Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted in the significant different expression genes. In all significant DEGs, some of the typical genes involved in osmoregulation, including carbonic anhydrase (CA), calcium/calmodulin-dependent protein kinase (CaM kinase) II (CAMK2), aquaporin-1(AQP1), sodium bicarbonate cotransporter (SLC4A4/NBC1), chloride channel 2(CLCN2), sodium/potassium/chloride transporter (SLC12A2 / NKCC1) and other osmoregulation genes were also identified. RNA-seq results were validated with quantitative real-time PCR (qPCR), the 17 random selected genes showed a consistent direction in both RNA-Seq and qPCR analysis, demonstrated that the results of RNA-seq were reliable.

Conclusions

The present results would be helpful to elucidate the osmoregulation mechanism of aquatic animals adapting to saline-alkali challenge. This study provides a global overview of gene expression patterns and pathways that related to osmoregulation in hybrid tilapia, and could contribute to a better understanding of the molecular regulation mechanism in different osmotic stresses.

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).

Dynamic immune and metabolism response of clam Meretrix petechialis to Vibrio challenge revealed by a time series of transcriptome analysis

Meretrix petechialis is an important commercial aquaculture species in China. During the clam culture period, mass mortality events often occurred due to the Vibrio infection. In this paper, M. petechialis were challenged with Vibrio parahaemolyticus immersion to simulate a natural infection, and the infection process were divided into four phases including latency, prodrome, onset and recovery phases based on the clam mortality data. Then, the dynamic response of clams to Vibrio infection at different infection phases were investigated by transcriptome analysis. A total of 38,067 differentially expressed genes (DEGs) were identified at different infection phases. DEG annotations showed that immune-related and metabolism-related signaling pathways were enriched, indicating that immune defense and metabolism process play key roles during bacterial infection. Three kinds of expression pattern were classified by cluster analysis, including U-shape, L-shape and inverted V-shape. Anabolism and cellular growth proliferation related signaling pathways were repressed (long-lasting or transient) during bacterial infection. However, the immune related signaling pathways with different immune functions showed induction expression or repression expression against bacterial infection, which indicated that immune system take different strategies against bacterial infection. Furthermore, some signaling pathways such as PI3K-Akt signaling pathway both involved in immune defense and cell metabolism. This study provides a sight that the dynamic immunity and metabolic responses may be integrated to improve the host survival and shift more energy for immune defense.

Liquid biopsies for omics-based analysis in sentinel mussels

Liquid biopsy of plasma is a simple and non-invasive technology that holds great promise in biomedical research. It is based on the analysis of nucleic acid-based biomarkers with predictive potential. In the present work, we have combined this concept with the FTA technology for sentinel mussels. We found that hemocytes collected from liquid biopsies can be readily fixed on FTA cards and used for long-term transcriptome analysis. We also showed that liquid biopsy is easily adaptable for metagenomic analysis of bacterial profiles of mussels. We finally provide evidence that liquid biopsies contained circulating cell-free DNA (ccfDNA) which can be used as an easily accessible genomic reservoir. Sampling of FTA-fixed circulating nucleic acids is stable at room temperature and does not necessitate a cold-chain protection. It showed comparable performance to frozen samples and is ideally adapted for sampling in remote areas, most notably in polar regions threatened by anthropogenic activities. From an ethical point of view, this minimally-invasive and non-lethal approach further reduces incidental mortality associated with conventional tissue sampling. This liquid biopsy-based approach should thus facilitate biobanking activities and development of omics-based biomarkers in mussels to assess the quality of aquatic ecosystems.

Two novel CYP3A isoforms in marine mussel Mytilus coruscus: Identification and response to cadmium and benzo[a]pyrene

CYP3A enzymes play a crucial role in metabolic clearance of a variety of xenobiotics. However, their genetic information and function remain unclear in molluscs. In the present study, two novel CYP3A genes i.e. McCYP3A-1 and McCYP3A-2 were identified and characterized from the thick shell mussel Mytilus coruscus, and their tissue distribution as well as the response to cadmium (Cd) and benzo[a]pyrene (B[α]P) exposure were addressed using real time quantitative RT-PCR (qRT-PCR) and erythromycin N-demethylase (ERND) assay. McCYP3A-1 and McCYP3A-2 possess typically domains of CYP family such as helix-C, helix-I, helix-K, PERF and the heme binding domain as well as the characteristic domains of CYP3s including six SRS motifs. McCYP3A-1 and McCYP3A-2 transcripts were constitutively expressed in all examined tissues with high expression level in digestive glands, hepatopancreas and gonads. Upon B[α]P exposure, McCYP3A-1 and McCYP3A-2 mRNA expression in digestive glands showed a pattern of up-regulation followed by down-regulation, while under Cd exposure, showed a time-dependent induction profile. In addition, ERND activity, generally used as an indicator of CYP3, increased in a time-dependent manner after exposure to Cd and B[α]P. These results collectively indicated that McCYP3A-1 and McCYP3A-2 are CYP3A family member and may play a potential role in metabolic clearance of xenobiotics. Meanwhile, the current results may provide some baseline data to support McCYP3A-1 and McCYP3A-2 as candidate biomarkers for monitoring of PAHs and heavy metal pollution.

A novel invertebrate toll-like receptor is involved in TLR mediated signal pathway of thick shell mussel Mytilus coruscus

Toll-like receptors (TLRs) are the most well studied pattern recognition receptors (PRRs) that play a crucial role in both innate and adaptive immunity in animals. In the present study, a novel toll-like receptor (McTLRj) was identified and characterised in thick shell mussel Mytilus coruscus. McTLRj possessed a signal peptide, a transmembrane domain, leucine-rich repeats and an intracellular Toll/interleukin-1 receptor domain that were conserved in typical TLRs. McTLRj transcripts were constitutively expressed in all of the examined tissues with high expression level in immune-related tissues, and significantly induced in haemocytes upon live Vibrio alginolyticus, lipopolysaccharide, polyinosinic-polycytidylic acid and peptidoglycans challenge. The overexpression of the McTLRj TLR fragment in Drosophila S2 cells could induce the expression of Drosophila attacin A, drosomycin, cecropin A, and metchnikowin expression. The expression of McTLRj was obviously repressed by dsRNA-mediated RNA interference, and downstream TLR pathway factors, such as MyD88a, IRAK4, and TRAF6 were significantly repressed in McTLRj-silenced mussels upon LPS challenge. These results collectively indicated that McTLRj is a TLR family member that may play a potential PRR role in TLR-mediated signalling pathway. This research contributed to the clarification of innate immune response in molluscs.

Response of a novel selenium-dependent glutathione peroxidase from thick shell mussel Mytilus coruscus exposed to lipopolysaccharide, copper and benzo[α]pyrene

Glutathione peroxidase (GPx) plays an important antioxidant role in cellular defense against environmental stress. In the present study, a novel selenium-dependent glutathione peroxidase termed McSeGPx firstly identified in thick shell mussel Mytilus coruscus. McSeGPx consists of 197 amino acid residues, characterized with one selenocysteine residue encoded by an opal stop codon TGA, one selenocysteine insertion sequence (SECIS) in the 3' untranslated region (UTR), two active site motifs and one signature sequence motif. McSeGPx transcripts were constitutively expressed in all examined tissues, and were significantly induced in gills and digestive glands with the stimulations of lipopolysaccharide (LPS), copper (Cu) and benzo[α]pyrene (B[α]P). Additionally, rough increases in McSeGPx activity were detected in both tissues under the challenge of LPS, Cu and B[α]P. Collectively, these results suggested that McSeGPx affiliate to selenocysteine dependent GPx (SeGPx) family and might play an important role in mediating the environmental stressors and antioxidant response in M. coruscus.

A novel invertebrate toll-like receptor with broad recognition spectrum from thick shell mussel Mytilus coruscus

Toll-like receptors (TLRs) are a category of most well recognized pattern recognition molecules that act on a vital role in both innate and adaptive immunity. In the present study, a novel toll-like receptor (McTLRw) was identified and characterized in thick shell mussel Mytilus coruscus. McTLRw possesses one intracellular Toll/interleukin-1 (IL-1) receptor (TIR) domain, one transmembrane region (TM), one leucine rich repeat N-terminal domain (LRR_NT) and a few of leucine-rich repeats (LRRs), which all are common in typical TLRs. McTLRw transcripts were constitutively expressed in all examined tissues with higher expression levels in immune related tissues, and were significantly induced in haemocytes with the challenges of live Vibrio alginolyticus, lipopolysaccharide (LPS), peptidoglycans (PGN) and β-glucan (GLU), but not induced by polyinosinic-polycytidylic acid (poly I:C). rMcTLRw exhibited affinity to LPS, PGN and GLU while no affinity to poly I:C. Further, the downstream of TLR signaling pathway myeloid differentiation factor 88a (MyD88a), interleukin-1 receptor-associated kinase-4 (IRAK4) and tumor necrosis factor receptor-associated factor 6 (TRAF6) were significantly repressed in McTLRw silenced mussels while challenged with LPS. These results collectively indicated that McTLRw is one member of TLR family and involved in immune response to against invaders by taking participate in TLR mediated signaling pathway.

High individual variability in the transcriptomic response of Mediterranean mussels to Vibrio reveals the involvement of myticins in tissue injury

Mediterranean mussels (Mytilus galloprovincialis) are sessile filter feeders that live in close contact with numerous marine microorganisms. As all invertebrates, they lack an adaptive immune response and how these animals are able to respond to a bacterial infection and discriminate it from their normal microbiome is difficult to understand. In this work, we conducted Illumina sequencing of the transcriptome of individual mussels before and after being infected with Vibrio splendidus. The control mussels were injected with filtered seawater. We demonstrate that a great variability exists among individual transcriptomes and that each animal showed an exclusive repertoire of genes not shared with other individuals. The regulated genes in both the control and infected mussels were also analyzed and, unexpectedly, the sampling before the injection was considered a stress stimulus strong enough to trigger and modulate the response in hemocytes, promoting cell migration and proliferation. We found a clear response against the injection of filtered seawater, suggesting a reaction against a tissue injury in which the myticins, the most expressed antimicrobial peptides in mussel, appeared significantly up regulated. Functional experiments with flow cytometry confirmed the transcriptomic results since a significant alteration of hemocyte structures and a decrease in the number of hemocytes positive for myticin C were found only after a Vibrio infection and not observed when mussels were bled before, generating a tissue injury. Therefore, we report the involvement of myticins in the response to a danger signal such as a simple injection in the adductor muscle.

A novel interleukin-1 receptor-associated kinase-4 from thick shell mussel Mytilus coruscus is involved in inflammatory response

Interleukin-1 receptor-associated kinase-4 (IRAK4) is considered as the most upstream kinase of IRAKs and plays a vital role in Toll-like receptor/Interleukin-1 receptor (TLR/IL-1R) signal transduction. In the present study, IRAK4 from thick shell mussel Mytilus coruscus (McIRAK4) was identified and characterized. McIRAK4 showed the most similarity to its counterparts in bivalves. The conserved death domain (DD) and catalytic domain of serine/threonine kinases (STKc) were predicted in all examined IRAK4s. McIRAK4 transcripts were constitutively expressed in all examined tissues with the higher expression level in immune related tissues, and were significantly induced in haemocytes upon lipopolysaccharide (LPS) and polyinosinic-polycytidylic acid (poly I:C) challenge. Further, the expression of McIRAK4 was obviously repressed by dsRNA mediated RNA interference (RNAi), meanwhile the proinflammatory cytokines TNF-alpha and IL17 were down-regulated while the antiinflammatory cytokine TGF-β was up-regulated. Additionally, McIRAK4 showed a global cytoplasmic localization in HEK293T cells through fluorescence microscopy. These results collectively indicated that McIRAK4 is one member of IRAK4 subfamily and might play the potential signal transducer role in inflammatory response. The present study provides supplement for TLR-mediated signaling pathway triggered by pathogenic invasions in thick shell mussel, and contributes to the clarification of the innate immune response in molluscs.

Identification and functional characterization of three myeloid differentiation factor 88 (MyD88) isoforms from thick shell mussel Mytilus coruscus

Myeloid differentiation factor 88 (MyD88) is a pivotal adapter protein that involved in interleukin-1 receptor/toll-like receptor (IL-1R/TLR) signal transduction, which could spur downstream cascades and eventually drawn into innate immune response. MyD88 has been extensively studied in vertebrates, however, the information ascribe to MyD88 in invertebrates is still very scarce especially its function annotation remains extremely obscure. At here, three novel MyD88 isoforms termed McMyD88a, McMyD88b and McMyD88c were firstly cloned from thick shell mussel Mytilus coruscus. McMyD88a, McMyD88b and McMyD88c shared domain topology containing the Death domain (DD) and TIR domain (TIR) with its counterparts in mammals. All three McMyD88s were ubiquitously expressed in examined tissues in thick shell mussel, with the higher expression levels in immune-related tissues such as haemocytes, gills and digestive glands. Upon Vibrio alginolyticus, polyinosine-polycytidylic acid (poly I:C) and lipopolysaccharide (LPS) challenge, McMyD88a, McMyD88b and McMyD88c transcripts were significantly induced in haemocytes despite of differential expression levels and responsive time points. Overexpression of McMyD88a, McMyD88b and McMyD88c showed a dose-dependent induction to NF-κB or ISRE in mammalian cell lines. Taken together, these results suggested that McMyD88a, McMyD88b and McMyD88c are members of MyD88 family and play potential roles in innate immune response to pathogenic invasions in thick shell mussel. Moreover, these results suggested indirectly the existence of a MyD88-dependent signaling pathway in thick shell mussel, and provide insight into the immunoregulatory effect in molluscs.

RNA-Seq Analysis for Assessing the Early Response to DSP Toxins in Mytilus galloprovincialis Digestive Gland and Gill

The harmful effects of diarrhetic shellfish poisoning (DSP) toxins on mammalian cell lines have been widely assessed. Studies in bivalves suggest that mussels display a resistance to the cytogenotoxic effects of DSP toxins. Further, it seems that the bigger the exposure, the more resistant mussels become. To elucidate the early genetic response of mussels against these toxins, the digestive gland and the gill transcriptomes of Mytilus galloprovincialis after Prorocentrum lima exposure (100,000 cells/L, 48 h) were de novo assembled based on the sequencing of 8 cDNA libraries obtained using an Illumina HiSeq 2000 platform. The assembly provided 95,702 contigs. A total of 2286 and 4523 differentially expressed transcripts were obtained in the digestive gland and the gill, respectively, indicating tissue-specific transcriptome responses. These transcripts were annotated and functionally enriched, showing 44 and 60 significant Pfam families in the digestive gland and the gill, respectively. Quantitative PCR (qPCR) was performed to validate the differential expression patterns of several genes related to lipid and carbohydrate metabolism, energy production, genome integrity and defense, suggesting their participation in the protective mechanism. This work provides knowledge of the early response against DSP toxins in the mussel M. galloprovincialis and useful information for further research on the molecular mechanisms of the bivalve resistance to these toxins.

Transcriptome analysis of Artemia sinica in response to Micrococcus lysodeikticus infection

To enhance genomic resources and understand the molecular immune mechanisms underlying the response topathogens, we first performed a comparative gene transcription analysis from Micrococcus lysodeikticus-immunized Artemia sinica and from a control group through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. In total, 80, 113, 984 clean reads were obtained and then assembled into 71,536 unigenes with an average length of 1115 bp and an N50 of 1783 bp. Unigenes were annotated by comparing against nr, Swiss-Prot\KEGG\ COG\ KOG\ GO and Pfam databases, and 27,689 unigenes (38.7%) were annotated in at least one database. After bacterial challenge, 183 and 298 genes were identified as remarkably up-regulated or down-regulated, respectively, amongst 481 were associated with 168 pathways, including classical immune-related pathways, such as 'Toll-like receptor signaling', 'the complement cascades', 'MAPK signaling pathway' and 'Apoptosis'. Besides, eight genes which were differently expressed immune-related were confirmed by using quantitative real-time PCR. This study characterized a gene expression pattern for normal and M. lysodeikticus -immunized A. sinica for the first time and sheds new light on the molecular mechanisms thus enabling future efforts on disease control programs in this valuable aquaculture species.

Transcriptomic analysis reveal an efficient osmoregulatory system in Siberian sturgeon Acipenser baeri in response to salinity stress

Sturgeons are euryhaline fish species that have developed specific mechanisms of osmotic and ion regulation to adapt to waters of varying salinity. For the aim to elucidate the osmoregulation strategy behind its high salinity tolerance of sturgeons, the transcriptomes of gills in Siberian sturgeon Acipenser baeri under salinity stress (30 ppt) were sequenced using deep-sequencing platform Illumina/HiSeq-2500 and differential expression genes (DEGs) were identified. A total of 167, 501, 278 clean reads were obtained and 280, 238 unigenes were composed of those clean reads with the mean length of 520nt, and the N50 of 630 bp. Unigenes Sequence alignment was implemented via KEGG, KOG, NT, NR, PFAM, Swiss-Prot, and GO databases. 62, 242 unigenes (22.21%) were annoated in at least one database. 11380 significantly differentially expressed unigenes were found, 6969 of which were up-regulated and 4411 were down-regulated by salinity stress. Amongst the top 20 KEGG pathways with the most amount of annotation sequences, some pathways such as glycerophospholipid metabolism, fatty-acid biosynthesis, glycolysis/gluconeogenesis, oxidative phosphorylation have been comprehensively proved to be relevant to osmoregulation. Despite of these, three possible osmoregulation-related signaling pathways as lipid metabolism related pathways, tight junction pathway and thyroid hormone signaling pathway have been widely analyzed in the current study. In all DEGs, some of the typical genes involved in osmoregulation, including calcium-transporting ATPase 4 (ATP2B4), Na⁺/K⁺-ATPase alpha subunit (α-NKA), potassium-transporting ATPase alpha chain 1 (ATP4A) and Ras GTPase-activating protein (RasGAP) etc were also identified. RNA-seq results were validated with quantitative real-time PCR (qPCR), the 12 selected genes showed a consistent direction in both DGE library and qPCR analysis, proving that the RNA-seq results are reliable. The present results would be helpful to elucidate the osmoregulation mechanism of aquatic animals adapting to salinity challenge.

Potential trade-offs between biomineralization and immunity revealed by shell properties and gene expression profiles of two closely related Crassostrea species

Species of the Ostreidae family are key ecosystem engineers and many of them - including Crassostrea gigas and Crassostreavirginica - are commercially important aquaculture species. Despite similarities in their morphology and ecology, these two species differ in their ability to defend against pathogens, potentially reflecting species-specific differential specialization of hemocytes on immune defense versus biomineralization. To test this hypothesis, we investigated the expression levels of immune- and biomineralization-related genes as well as mineralogical and mechanical properties of the shells and the calcium sequestration ability of the hemocytes of C. gigas and C. virginica The expression of biomineralization-related genes was higher in C. virginica than in C. gigas in multiple tissues including the mantle edge and hemocytes, while the expression of immune genes was higher in the hemocytes of C. gigas Hemocytes of C. virginica contained more calcium (stored intracellularly as calcium carbonate mineral) compared with those of C. gigas Analysis of the adult shells showed that the crystallinity of calcite was higher and the laths of the foliated layer of the shell were thicker in C. virginica than in C. gigas Mechanically, the shells of C. virginica were stiffer, harder and stronger than those of C. gigas Taken together, our results show that the species-specific differences in physiology (such as disease resistance and exoskeleton properties) are reflected at the cellular and molecular levels in the differential specialization of hemocytes on potentially competing functions (immunity and biomineralization) as well as different expression profiles of other tissues involved in biomineralization (such as the mantle edge).

Molecular cloning and functional analysis of tumor necrosis factor receptor-associated factor 6 (TRAF6) in thick shell mussel, Mytilus coruscus

Tumor necrosis factor receptor-associated factor 6 (TRAF6) is one of the key adapter molecules in Toll-like receptor signal transduction that triggers downstream cascades involved in innate immunity. Despite of the well study in vertebrates, there is few data ascribe to this TRAF member in invertebrates, especially in bivalves. In the present study, a novel TRAF6 homologue termed McTRAF6 was firstly characterized in Mytilus coruscus. Like its counterparts in mammals, McTRAF6 shared the domain topology containing one RING domain, two zinc finger domains, one coiled-coil region and a MATH domain. McTRAF6 transcripts predominantly expressed in gills, digestive glands and hemocytes in M. coruscus, and were significantly up-regulated in hemocytes after challenge with lipopolysaccharide (LPS) and polyinosine-polycytidylic acid (poly I:C). Further, the subcellular localization in cytoplasm and the activation of Nk-κB or ISRE luciferase reporter by overexpressed McTRAF6 were identified in HEK293T cells. These results collectively indicate that McTRAF6 is a member of TRAF6 subfamily and plays a potential role in immune defense system against pathogenic agents invasions in thick shell mussel. To our knowledge, this is the first report on component of TLR signaling pathway in thick shell mussel, providing further evidence for the existence of TLR pathway in M. coruscus and contribute to clarify the innate immune system of thick shell mussel.

The receptor for activated C kinase 1 (RACK1) mediating immune response in thick shell mussel Mytilus coruscus

The receptor for activated C kinase 1 (RACK1) is a intracellular receptor for the protein kinase C family which mediates various biological processes. Here, a novel RACK1 gene termed Mc-RACK1 was identified from thick shell mussel, Mytilus coruscus. Mc-RACK1 shared typical RACK1 domains containing WD repeats, PKC phosphorylation sites, N-myristoylation sites, PKC activation sites, TK phosphorylation site and WD motifs. Mc-RACK1 was constitutively expressed in all examined tissues, and its expression in gills, haemocytes and digestive glands were significantly up-regulated upon LPS challenge. Mc-RACK1 showed a significantly down-regulated expression in gills and haemocytes at the early phase upon copper exposure. Mc-RACK1 in haemocytes was silenced after receiving its dsRNA, meanwhile, the increases of SOD and CAT activity were investigated. Further, Mc-RACK1 could activate the NF-κB and ISRE reporter in HEK-293T cells. These suggested that Mc-RACK1 had a deeper involvement in mollusc immunity, and played an important role in antioxidant system.

Molecular characterization of complement component 3 (C3) in Mytilus coruscus improves our understanding of bivalve complement system

Complement component 3 (C3) plays a central role in the complement system whose activation is essential for all the important functions performed by this system. Here, a novel C3 gene, termed Mc-C3, was identified from thick shell mussel (Mytilus coruscus). The deduced Mc-C3 protein possessed the characteristic structure features present in its homologs and contained the A2M_N_2, ANATO, A2M, A2M_comp, A2M_recep, and C345C domains, as well as the C3 convertase cleavage site, thioester motif, and conserved Cys, His, and Glu residues. Mc-C3 gene constitutively expressed in all examined tissues and predominantly expressed in immune-related tissues such as gills, hemocytes and hepatopancreas. After stimulation with lipopolysaccharide or Cu²⁺, the expression of Mc-C3 was significantly induced in gills. Further luciferase reporter assays showed the ability for activation of NF-κB signaling transduction of Mc-C3a. Taken together, these results show that C3 may play an essential role in the immune defense of M. coruscus. The present data therefore provide a more detailed insight into the functional activities of the bivalve complement system.