Transcriptome profiles of the clam Meretrix petechialis hepatopancreas in response to Vibrio infection

Integrated Vibrio load variation and transcriptome profiles provide new insights into the defensive response of Cyclina sinensis under Vibrio parahaemolyticus infection

Vibrio parahaemolyticus is the dominant pathogen in mariculture, leading to the bivalves' mass summer mortality. This study investigated the antibacterial defense mechanism of the Cyclina sinensis against V. parahaemolyticus. The immersion challenge revealed that the daily mortality of C. sinensis increased and then decreased gradually, and the mass mortality occurred about a week after V. parahaemolyticus infection. The Vibrio load in the hepatopancreas of C. sinensis was dramatically increased at 12-24 hpi (hours post-infection) and then declined significantly at 2-7 dpi (days post-infection). RNA-Seq generated 866 differentially expressed genes (DEGs), and KEGG analyses enriched multiple innate immune-related and metabolic-related pathways. The expression levels of nine immune-related and metabolic-related DEGs were significantly changed after Vibrio infection, and their temporal expression patterns were multiple. Our results indicated that the immunity and metabolic responses might be reprogrammed to protect the host against pathogens at the early infection phase. This study would expand our knowledge of the pathogenesis mechanisms of clams infected with Vibrio and provide a theoretical basis for healthy shellfish cultivation.

Growth Performance, Dietary Enzyme Profiling, and Antioxidant-Induced Immunity of Crassostrea belcheri Fed With Microalgal Diet

Diseases, overfishing, and habitat loss are constantly reducing oyster populations. Moreover, environmental pollution and natural disasters also hinder offshore shellfish cultivation. Therefore, filter-feeding bivalves could be cultivated on land-based system to maintain the culture condition and secure food safety. Four different diets including Chaetoceros gracilis (CG), Tetraselmis chuii (TC), mixture of CG and TC (CG/TC), and fresh seawater without feed (Con) were trialed for Crassostrea belcheri spat in this research. After 35 days of culture, the highest survival rate (SR), volume, and weight increment (WI) of oyster and improved water quality appeared in CG/TC. In the same manner, CG/TC diet exhibited greater lipase (LPS), pepsin (PES), catalase (CAT), alkaline phosphatase (AKP), and hydrogen peroxide (H2O2) activity. Conversely, amylase (AMS), acid phosphatase (ACP), superoxide dimutase (SOD), malondialdehyde (MDA), and lysozymes (LZMs) activities were substantially higher in CG diet compared to TC, CG/TC, and Con. No significant differences were observed among CG, TC, and CG/TC for antioxidant capacity (AOC). In this investigation, mixed algal diet had excellent results for the growth and development of oyster, whereas unialgal diet improved immunity and AOC to survive in unfavorable conditions. This observation help elucidates the knowledge on microalgal diet influenced immune modulation and health of marine bivalves in the scenario of land-based farming.

The combination of high temperature and Vibrio infection worsens summer mortality in the clam Meretrix petechialis by increasing apoptosis and oxidative stress

The interaction between environmental factors and Vibrio in bivalves is not well understood, despite the widely held belief that pathogen infection and seawater temperature significantly impact summer mortality. In the present study, we conducted simulated experiments to explore the effects of high temperature and Vibrio infection on the clam Meretrix petechialis. The survival curve analysis revealed that the combined challenge of high temperature and Vibrio infection (31°C-vibrio) led to significantly higher clam mortality compared to the groups exposed solely to Vibrio (27°C-vibrio), high temperature (31°C-control), and the control condition (27°C-control). Furthermore, PCoA analysis of 11 immune genes indicated that Vibrio infection predominated during the incubation period, with a gradual equilibrium between these factors emerging during the course of the infection. Additionally, our investigations into apoptosis and autophagy processes exhibited significant induction of mTOR and Bcl2 of the 31°C-vibrio group in the early challenge stage, followed by inhibition in the later stage. Oxidative stress analysis demonstrated a substantial additive effect on malondialdehyde (MDA) and glutathione (GSH) content in the combined challenge group compared to the control group. Comparative transcriptome analysis revealed a significant increase in differentially expressed genes related to immunity, such as complement C1q-like protein, C-type lectin, big defensin, and lysozyme, in the 31°C-vibrio group, suggesting that the synergistic effect of high temperature and Vibrio infection triggers more robust antibacterial immune responses. These findings provide critical insights for understanding the infection process and uncovering the causes of summer mortality.

Comprehensive Multi-omics Approaches Provide Insights to Summer Mortality in the Clam Meretrix petechialis

Bivalve mass mortalities have been reported worldwide, which not only can be explained as a result of pathogen infection, but may reflect changes in environments. Although these episodes were often reported, there was limited information concerning the molecular responses to various stressors leading to summer mortality. In the present work, RNA sequencing (RNA-seq), tandem mass tagging (TMT)-based quantitative proteomics, and 16S rRNA sequencing were used to explore the natural outbreak of summer mortality in the clam Meretrix petechialis. We identified a total of 172 differentially expressed genes (DEGs) and 222 differentially expressed proteins (DEPs) in the diseased group compared to the normal group. The inconsistent expression profiles of immune DEGs/DEPs may be due to the immune dysregulation of the diseased clams. Notably, 11 solute carrier family genes were found among the top 20 down-regulated genes in the diseased group, indicating that weakened transmembrane transport ability might occur in the diseased clams. Integration analysis of transcriptomic and proteomic results showed that many metabolic processes such as "arginine and proline metabolism" and "tyrosine metabolism" were inhibited in the diseased group, suggesting metabolic inhibition. Moreover, 16S rRNA sequencing revealed that the microbial composition of clam hepatopancreas was disordered in the diseased group. The comparison of DEGs expression between the natural summer mortality event and an artificial challenge experiment involving both Vibrio infection and heat stress revealed 9/15 genes showing similar expression trends between the two conditions, suggesting that the summer mortality might be caused by a combination of high temperature and Vibrio infection. These results would deepen our understanding of summer mortality and provide candidate resistance markers for clam resistance breeding.

The physiological responses to titanium dioxide nanoparticles exposure in pearl oysters (Pinctada fucata martensii)

To evaluate the physiological responses to titanium dioxide nanoparticles exposure in pearl oysters (Pinctada fucata martensii), pearl oysters were exposed for 14 days to different levels (0.05, 0.5, and 5 mg/L) of nano-TiO2 suspensions, while a control group did not undergo any nano-TiO2 treatment. And then recovery experiments were performed for 7 days without nano-TiO2 exposure. At days 1, 3, 7, 14, 17, and 21, hepatopancreatic tissue samples were collected and used to examine the activities of protease, amylase, lipase, catalase (CAT), glutathione peroxidase (GPx), superoxide dismutase (SOD), lysozyme (LYS), alkaline phosphatase (AKP), and acid phosphatase (ACP). The microstructure of the nacreous layer in shell was also analyzed by scanning electron microscopy. Results showed that pearl oysters exposed to 5 mg/L of TiO2 nanoparticles had significantly lower protease, amylase, and lipase activities and significantly higher CAT, SOD, GPx, LYS, ACP, and AKP activities than control pearl oysters did even after 7-day recovery (P-values <0.05). Pearl oysters exposed to 0.5 mg/L or 0.05 mg/L of TiO2 nanoparticles had lower protease, amylase, and lipase activities and higher CAT, SOD, GPx, LYS, ACP, and AKP activities than control pearl oysters did during the exposure period. After 7-day recovery, no significant differences in protease, lipase, SOD, GPx, CAT, ACP, AKP, or LYS activities were observed between pearl oysters exposed to 0.05 mg/L of TiO2 nanoparticles and control pearl oysters (P-values >0.05). In the period from day 7 to day 14, indistinct and irregular nacreous layer crystal structure in shell was observed. This study demonstrates that TiO2 nanoparticles exposure influences the levels of digestion, immune function, oxidative stress, and biomineralization in pearl oysters, which can be partially and weakly alleviated by short-term recovery. These findings contribute to understanding the mechanisms of action of TiO2 nanoparticles in bivalves. However, studies should evaluate whether a longer recovery period can restore to their normal levels in the future.

Ocean acidification induces tissue-specific interactions with copper toxicity on antioxidant defences in viscera and gills of Asiatic hard clam Meretrix petechialis (Lamarck, 1818)

Toxicity of contaminants in organisms under ocean acidification (OA) has attracted increasing attention in ecotoxicological studies. This study investigated how pCO₂-driven OA affected waterborne copper (Cu) toxicity in antioxidant defences in viscera and gills of Asiatic hard clam Meretrix petechialis (Lamarck, 1818). Clams were continuously exposed to Cu at ambient relevant (0/no metal exposure, 10 and 50 μg L^-1) and polluted-high (100 μg L^-1) concentrations in unacidified (pH 8.10) and acidified (pH 7.70/moderate OA and 7.30/extreme OA) seawater for 21 days. Following coexposure, metal bioaccumulation and responses of antioxidant defence-related biomarkers to OA and Cu coexposure were investigated. Results showed that metal bioaccumulation was positively correlated with waterborne metal concentrations but was not notably influenced by OA conditions. Both Cu and OA affected the antioxidant responses to environmental stress. Additionally, OA induced tissue-specific interactions with Cu on antioxidant defences, varying with exposure conditions. In unacidified seawater, antioxidant biomarkers were activated to defend against oxidative stress induced by Cu and prevented clams from lipid peroxidation (LPO or MDA), but failed to defend against DNA damage (8-OHdG). OA exacerbated Cu toxicity in antioxidant defences and increased LPO levels in tissues. Gills and viscera adopted adaptive antioxidant defence strategies to manage oxidative stress, with the former being more vulnerable to oxidative stress than the latter. MDA and 8-OHdG were sensitive to OA and Cu exposure, respectively, and were useful bioindicators for assessing oxidative stress. Integrated biomarker response (IBR) and PCA can reflect the integrative responses of antioxidant biomarkers to environmental stress and illuminate the contributions of specific biomarkers to antioxidant defence strategies. The findings provided insights for understanding antioxidant defences against metal toxicity in marine bivalves under OA scenarios, which is essential into managing wild populations.

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.

Integrated RNA-seq and RNAi Analysis of the Roles of the Hsp70 and SP Genes in Red-Shell Meretrix meretrix Tolerance to the Pathogen Vibrio parahaemolyticus

The "Wanlihong" Meretrix meretrix (WLH-M) clam is a new variety of this species that has a red shell and stronger Vibrio tolerance than ordinary M. meretrix (ORI-M). To investigate the molecular mechanisms responsible for the WLH-M strain's tolerance to Vibrio, we challenged clams with Vibrio parahaemolyticus and then assessed physiological indexes and conducted transcriptome analysis and RNA interference experiments. The mortality, tissue bacterial load, and hemocyte reactive oxygen species level of ORI-M were significantly higher than those of WLH-M, whereas the content and activity of lysozyme were significantly lower. Gene Ontology functional annotation analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that immune and metabolic pathways were enriched in Vibrio-challenged clams. The expressions of the heat shock protein 70 (Hsp70) and serine protease (SP) genes, which are involved in antibacterial immunity, were significantly upregulated in WLH-M but not in ORI-M, while the expression of the kynurenine 3-monooxygenase gene, a proinflammatory factor, was significantly downregulated in WLH-M. RNA interference experiments confirmed that Hsp70 and SP downregulation could result in increased mortality of WLH-M. Therefore, we speculate that Hsp70 and SP may be involved in the antibacterial immunity of WLH-M in vivo. Our data provided a valuable resource for further studies of the antibacterial mechanism of WLH-M and provided a foundation for the breeding of pathogen-resistant strains.

Transcriptome analysis of immune-related genes in Sesarmops sinensis hepatopancreas in reaction to peptidoglycan challenge

Sesarmops sinensis is a dominant omnivorous crab species, which plays an important ecological function in salt marsh ecosystems. To better understand its immune system and immune related genes under pathogen infection, the transcriptome was analyzed by comparing the data of S. sinensis hepatopancreas stimulated by PBS and PGN. A set of assembly and annotation identified 39,039 unigenes with an average length of 1105 bp, obtaining 1300 differentially expressed genes (DEGs) in all, which included 466 remarkably up-regulated unigenes and 834 remarkably down-regulated unigenes. In addition, based on mensurable real time-polymerase chain reaction and high-throughput sequencing, several immune responsive genes were found to be markedly up-regulated under PGN stimulation. In conclusion, in addition to enriching the existing transcriptome data of S. sinensis, this study also clarified the immune response of S. sinensis to PGN stimulation, which will help us to further understand the crustacean's immune system.

Transcriptomic analysis reveals the immune changes associated with reproduction in the clam Meretrix petechialis

Substantial mortality and economic losses in marine mollusk culture has drawn considerable attention in recent years. The changes in immune status and environmental stress are thought to be the main causes of shellfish summer mortality. The reproduction and immune defense are both physiologically demanding processes, therefore, the immune status of mollusk is likely to be affected by reproduction during breeding. In present study, we performed transcriptome and gene expression analyses in the clam Meretrix petechialis pre-/post-spawning. DEGs enrichment analysis revealed important immune signaling pathways and key genes changed after spawning. Further analysis showed females up-regulated genes involved in apoptosis, TLR signal pathway and heat shock, whereas males down-regulated complement-related genes after spawning. Additionally, both genders of clams up-regulated its immune response level to against Vibrio infection after spawning revealed by the changes of four immune-related DEGs. The up-regulation of two marker genes at the transcription and protein levels further confirmed that pathogen reinforced the expression differences of immune-related genes between the two groups. Our study provides a new insight into the understanding of molecular mechanisms underlying reproduction influenced immune differences in M. petechialis.

TAF5L functions as transcriptional coactivator of MITF involved in the immune response of the clam Meretrix petechialis

TAF5L is a component of the P300/CBP-associated factor (PCAF) histone acetylase complex, which serves as a coactivator and takes part in basal transcription such as promoter recognition, complex assembly and transcription initiation. In our study, the full-length sequence of MpTAF5L was identified and characterized in the clam M. petechialis. Sequence analysis showed that the predicted MpTAF5L protein had a N-terminal TAF5-NTD2 domain and a C-terminal WD40-repeats domain. The annotation and evolutionary analysis revealed MpTAF5L had close evolutionary relationship with other invertebrate species. Tissue distribution analysis of TAF5L claimed that it was highly expressed in the mantle, adductor muscle, foot and hepatopancreas. The mRNA expression of MpTAF5L was significantly up-regulated after Vibrio parahaemolyticus challenge, indicating its involvement in the immune response of clam. Yeast two-hybrid assays verified that MpTAF5L can interact with MpMITF (a critical immune-related transcription factor), and our further research clarified this interaction depended upon the N-terminal TAF5-NTD2 domain of MpTAF5L. Moreover, the mRNA expression of MpBcl-2 (a target gene of MITF) was significantly decreased but the mRNA expression of MpMITF was not significantly changed after knockdown of MpTAF5L, which indicated the reduction of MpMITF regulating activity at the same time. These results revealed that MpTAF5L interacted with MpMITF and enhanced the activation of MpMITF, which plays roles in the immune defense against V. parahaemolyticus.

Characterization and expression analysis of immune-related genes in the red swamp crayfish, Procambarus clarkii in response to lipopolysaccharide challenge

To learn more about red swamp crayfish related genes in response to bacterial infections, we investigated immune-related genes induced by lipopolysaccharide (LPS) in the hepatopancreas using high-throughput sequencing method. In present the study, a total of 55,107 unigenes were identified, with an average length of 678 bp. A total of 2215 differentially expressed genes (DEGs) were found, including 669 up-regulated genes and 1546 down-regulated genes. The result of Gene ontology (GO) analysis revealed that 3017 DEGs were enriched in 19 biological process subcategories, 17 cellular component subcategories and 15 molecular function subcategories. The top 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that "ribosome" was the most abundant group, which had 34 DEGs. KEGG enrichment analysis identified several immune response pathways. Real-time quantitative reverse transcription-PCR (qRT-PCR) results exhibited that several immune responsive genes were greatly up-regulated following LPS stimulation as observed in the results of high-throughput sequencing. Overall, this study provides new insight into the immune defense mechanisms of P. clarkii against LPS infection.

p38 MAPK is involved in the immune response to pathogenic Vibrio in the clam Meretrix petechialis

p38 mitogen-activated protein kinases (MAPKs) are involved in the response to various extracellular stimuli via regulating gene expression. In the present study, a p38 MAPK gene (MpP38) was identified from the clam Meretrix petechialis. The full-length cDNA of MpP38 measures 1,720 bp, consisting of a 134-bp 5'-UTR, a 1,095-bp ORF and a 491-bp 3'-UTR. Both the mRNA and protein expression levels of MpP38 increased after Vibrio challenge, implying that MpP38 is involved in clam immunity. Based on our previous study, a transcription factor activated by p38 MAPK, i.e. microphthalmia-associated transcription factor (MITF), participated in clam immunity by regulating the expression of phenoloxidase (PO). Coupled with other related reports, the mechanism underlying the involvement of MpP38 in clam immunity is most likely that pathogen stimuli induce the phosphorylation of p38 MAPK and thus activate MITF to regulate the expression of the immune-related gene PO. The results obtained in this study support this mechanism. First, we found that the MpP38 phosphorylation level increased in response to Vibrio challenge. Second, as revealed by a yeast two-hybrid assay, there was a direct interaction between MpP38 and MITF. Meanwhile, inhibiting the phosphorylation of MpP38 decreased the phosphorylation level of MpMITF, implying that MpP38 phosphorylation is required for MpMITF activation. Additionally, our results showed that there was a regulatory relationship between MpP38 phosphorylation level and PO expression level. With increased MpP38 phosphorylation level, the PO expression level was also increased after Vibrio challenge; when MpP38 phosphorylation was inhibited, the PO expression level was significantly decreased. This study describes the immune function of p38 MAPK in the clam for the first time and analyses its potential underlying mechanism, which will help to elucidate the immune mechanism in the clam M. petechialis.

Effect of vitamin D3 on immunity and antioxidant capacity of pearl oyster Pinctada fucata martensii after transplantation: Insights from LC-MS-based metabolomics analysis

Postoperative care is a critical step of pearl culture that ultimately determines culture success. To determine the effect of dietary vitamin D3 (VD3) levels on immunity and antioxidant capacity of pearl oyster Pinctada fucata martensii during postoperative care and explore the mechanisms behind this phenomenon, five isonitrogenous and isolipidic experimental diets were formulated by adding different levels of dietary VD3 (0, 500, 1000, 3000, and 10000 IU/kg), and the diets were fed to five experimental groups (EG1, EG2, EG3, EG4, and EG5) in turn and cultured indoors. The control group (CG) was cultured in the natural sea. Pearl oysters that were 1.5 years old were subjected to nucleus insertion. After culturing for 30 days, EG3 exhibited significantly higher survival rates than those in CG and EG5 (P < 0.05). Moreover, EG3 exhibited the highest activities of alkaline phosphatase, acid phosphatase, catalase, superoxide dismutase, and lysozyme. However, EG5 achieved the highest activities of glutathione peroxidase. Metabolomics-based profiling of pearl oysters fed with high levels of dietary VD3 (EG5) and optimum levels of dietary VD3 (EG3) revealed 76 significantly differential metabolites (SDMs) (VIP > 1 and P < 0.05). Pathway analysis indicated that SDMs were involved in 21 pathways. Furthermore, integrated key metabolic pathway analysis suggested that pearl oysters in EG5 regulated the pentose phosphate pathway, glutathione metabolism, sphingolipid metabolism, and arachidonic acid metabolism in response to stress generated from excessive VD3. These findings had significant implications on strengthening the future development and application of VD3 in aquaculture of pearl oyster P. f. martensii.

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.

MITF Regulates Downstream Genes in Response to Vibrio parahaemolyticus Infection in the Clam Meretrix Petechialis

The microphthalmia-associated transcription factor (MITF) is a basic helix-loop-helix-leucine zipper protein that plays a key role in cell proliferation, survival and immune defense through the direct transcriptional control of downstream genes. We have found that MITF participates in the immune response to Vibrio parahaemolyticus infection in the clam Meretrix petechialis. In this study, we focused on how MITF functions in immunity. First, PO, CTSK, and BCL-2 were identified as the target genes of MpMITF in the clam by RNAi. EMSAs showed direct binding between the MpMITF protein and the E-box of the MpPO, MpCTSK, and MpBCL-2 promoters. Yeast one-hybrid assays also suggested that MpMITF could activate the expression of these three downstream genes. These results demonstrated that the transcriptional expression of MpPO, MpCTSK, and MpBCL-2 is directly regulated by MpMITF. Second, we analyzed the roles of MpPO, MpCTSK, and MpBCL-2 in clam immunity. The mRNA expression of MpPO, MpCTSK, and MpBCL-2 increased significantly after V. parahaemolyticus challenge, which implied that these genes might take part in the immune defense against V. parahaemolyticus challenge in clams. The purified recombinant proteins, MpPO and MpCTSK, inhibited the growth of V. parahaemolyticus. Additionally, the apoptosis rate of clam haemocytes rose significantly when the activity of MpBCL-2 was suppressed. These results revealed that MpPO, MpCTSK, and MpBCL-2 are involved in the immune defense against V. parahaemolyticus. This study supports the idea that the MpMITF pathway plays a key role in immune defense through the direct regulation of the downstream genes MpPO, MpCTSK, and MpBCL-2 in the clam, M. petechialis.

Transcriptome analysis and discovery of genes involved in immune pathways in Solen strictus (Gould, 1861) under Vibrio anguillarum

The Solen strictus (Gould, 1861) has been recognized as an important marine economic bivalve in Eastern and Southeast Asia. To gain a better understanding of the S. strictus immune system and its related genes in response to bacterial infections, we performed a comparative gene transcription analysis from S. strictus with Vibrio anguillarum through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. After assembly, a total of 195,774 transcripts with an average length of 996 bp were obtained. Total 153,038 unigenes were annotated in the nr, Swiss-Prot, KEGG, COG, KOG, GO and Pfam databases, and 56,597 unigenes (36.98%) were annotated in at least one database. After bacterial challenge, there were 1588 significant differentially expressed genes (DEGs) between the challenged and control groups, including 999 up-regulated and 589 down-regulated genes. All the DEGs were classified into three gene ontology categories, and allocated to 225 KEGG pathways. Immune-related genes were detected from immune system pathways among the top 20 enriched pathways, such as Toll-like receptor signaling, RIG-I-like receptor signaling and NOD-like receptor signaling pathway. In addition, 56,079 potential simple sequence repeats (SSRs) and 1,031,521 candidate single nucleotide polymorphisms (SNPs) were detected and identified in the S. strictus transcriptome. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on S. strictus. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

Response to different dietary carbohydrate and protein levels of pearl oysters (Pinctada fucata martensii) as revealed by GC-TOF/MS-based metabolomics

Land-based culturing can avoid the effects of environmental pollution and natural disasters, thus ensuring food safety for shellfish. However, food availability, in this case, is limited. To achieve the optimum balance of dietary carbohydrates and proteins and explore the mechanisms behind the phenomenon, we formulated five isoenergetic and isolipidic diets (C30P40, C35P35, C40P30, C45P25, and C50P20) with different levels of carbohydrates (C) and proteins (P). There were five experimental groups (C30P40, C35P35, C40P30, C45P25, and C50P20) and two control groups (CG1 and CG2). CG1 was fed with mixed powders of yeast and Chlorella sp., and CG2 was cultured in natural sea. After 60-day feeding, the highest rates of survival and absolute growth appeared in C45P25. C45P25 exhibited significantly higher activities of amylase, protease, alkaline phosphatase, acid phosphatase, superoxide dismutase, catalase, glutathione peroxidase, and phenoloxidase and significantly lower malondialdehyde content than C30P40, C35P35, C40P30, C50P20, and CG1. No significant differences were observed between C45P25 and CG2. Furthermore, the total antioxidant capacity of the pearl oysters in C45P25 was significantly higher than that in C30P40, C35P35, C40P30, and C50P20. On the basis of these results, the optimal balance of proteins and carbohydrates for pearl oysters was the C45P25 diet. Metabolomics-based profiling of the pearl oysters fed with high-carbohydrate/low-protein diet (C45P25) and low-carbohydrate/high-protein diet (C30P40) revealed 80 significantly different metabolites (VIP > 1 and P < 0.05). Furthermore, integrated key metabolic pathway analysis showed that C45P25 regulated starch and sucrose metabolism, alanine, aspartate and glutamate metabolism and glycine, serine and threonine metabolism to meet the energy demand and increase the glucogenic amino acid, thereby promoting protein synthesis and reducing fatty acid β-oxidation in comparison with C30P40. This finding helps elucidate the underlying mechanisms leading to the high-carbohydrate/low-protein diet characteristic of the optimal dietary carbohydrate and protein levels of P. f. martensii.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

The evolution of ependymin-related proteins

Background

Ependymins were originally defined as fish-specific secreted glycoproteins involved in central nervous system plasticity and memory formation. Subsequent research revealed that these proteins represent a fish-specific lineage of a larger ependymin-related protein family (EPDRs). EPDRs have now been identified in a number of bilaterian animals and have been implicated in diverse non-neural functions. The recent discoveries of putative EPDRs in unicellular holozoans and an expanded EPDR family with potential roles in conspecific communication in crown-of-thorns starfish suggest that the distribution and diversity of EPDRs is significantly broader than currently understood.

Results

We undertook a systematic survey to determine the distribution and evolution of EPDRs in eukaryotes. In addition to Bilateria, EPDR genes were identified in Cnidaria, Placozoa, Porifera, Choanoflagellatea, Filasterea, Apusozoa, Amoebozoa, Charophyta and Percolozoa, and tentatively in Cercozoa and the orphan group Malawimonadidae. EPDRs appear to be absent from prokaryotes and many eukaryote groups including ecdysozoans, fungi, stramenopiles, alveolates, haptistans and cryptistans. The EPDR family can be divided into two major clades and has undergone lineage-specific expansions in a number of metazoan lineages, including in poriferans, molluscs and cephalochordates. Variation in a core set of conserved residues in EPDRs reveals the presence of three distinct protein types; however, 3D modelling predicts overall protein structures to be similar.

Conclusions

Our results reveal an early eukaryotic origin of the EPDR gene family and a dynamic pattern of gene duplication and gene loss in animals. This research provides a phylogenetic framework for the analysis of the functional evolution of this gene family.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1306-y) contains supplementary material, which is available to authorized users.

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.

Identification of a gene encoding microphthalmia-associated transcription factor and its association with shell color in the clam Meretrix petechialis

The microphthalmia-associated transcription factor (MITF) is a master regulator of melanocyte development through the direct transcriptional control of related genes, e.g., the phenoloxidase gene. In this study, an MITF gene, MpMITF2, was identified in the clam Meretrix petechialis. The full-length cDNA of MpMITF2 was 2026 bp, and the molecular mass of the predicted protein was 42.6 kDa. A basic helix-loop-helix leucine zipper domain was detected in the deduced protein sequence, which can bind the E-box motif within the promoter of the downstream genes. The mRNA of MpMITF2 was more highly expressed in the mantle compared to the other four tissues. Furthermore, there was a significant difference in the mRNA expression of MpMITF2 among three clam strains with different shell colors. The protein level of MpMITF2 was also different among these strains. These results implied that MpMITF2 was associated with shell color formation in the clam M. petechialis. When the mRNA expression of MpMITF2 was knocked down, the new shell showed discontinuous pigment distribution, suggesting that the reduced expression of MpMITF2 influenced pigment synthesis. A gene encoding phenoloxidase (MpPO) was identified as related to the shell color of the clam and was also a putative downstream gene of MITF. Both the mRNA and protein levels of MpPO decreased significantly at 12 h post-MpMITF-suppression, suggesting that MpMITF2 is required for the expression of MpPO. Our results indicate the close relationships among MpMITF2, MpPO and shell color. This study implicates the role of MITF in shell color formation in the clam M. petechialis.