New insights into the Manila clam and PAMPs interaction based on RNA-seq analysis of clam through in vitro challenges with LPS, PGN, and poly(I:C)

Prediction of influenza virus infection based on deep learning and peripheral blood proteomics: A diagnostic study

INTRODUCTION

Influenza viruses cause seasonal epidemics almost every year, and it is difficult to diagnose quickly and accurately. Machine learning and peripheral blood protein omics have brought new ideas to the research of clinical markers.

OBJECTIVES

Prediction of key molecular marker of influenza virus infection by the established machine learning model and peripheral blood protein omics.

METHODS

This study used the testing data of 850 patients (including influenza, COVID-19 and mixed infections) and 265 healthy individuals, to establish and validate a diagnostic prediction model for influenza infection and verified the potential value of this model in the differential diagnosis of influenza, COVID-19 and healthy people.

RESULTS

The overall analysis showed that there were significant differences in 9 clinical features in the influenza group. Principal component analysis can effectively group samples based on these clinical features. Based on the random forest model and LASSO regression model found that the selected features are clinical indicators that can accurately distinguish influenza patients. We performed proteome sequencing combined with machine learning and found a total of 26 DEPs. Through PPI and WGCNA analysis, we identified several genes related to the proportion of monocytes. We then analyzed the correlation of these factors with immune cell proportions and found that SAA1 and SAA2 were highly correlated with various vital immunocyte. ROC curve analysis shows that SERPINA3 can distinguish influenza, COVID-19, mixed infection and healthy people; SAA1 can distinguish COVID-19, mixed infection and healthy people; SAA2 can distinguish influenza and healthy people. In influenza, high expression of SERPINA3, SAA1, and SAA2 is associated with higher risk. Finally, we used the ELISA method to confirm that SAA2 protein can be used as an auxiliary diagnostic indicator for influenza infection.

CONCLUSIONS

Preliminary results showed that SAA2 is an important molecular marker specific to influenza infection.

RNA-seq analysis of LPS-induced immune priming in silkworms (Bombyx mori) and the role of cytochrome P450 detoxification system in the process

While immune priming has been identified in many invertebrates, the intricate mechanisms that drive this process in insects continue to be a subject of mystery. In this study, we exposed silkworm larvae to varying doses of lipopolysaccharide (LPS) to induce immune priming and assessed their survival upon challenge with Bacillus thuringiensis (Bt). Transcriptome analysis was performed to identify differentially expressed genes (DEGs) associated with immune priming. The role of CYP450 genes in this process was further explored using RNA interference (RNAi) to knockdown CYP9E2 and CYP6K1, followed by measurements of detoxification enzyme activities and reactive oxygen species (ROS) levels. We found that LPS exposure significantly increased silkworm survival rates upon Bt challenge, indicating the induction of immune priming. Transcriptome analysis revealed 549 DEGs, including a large number involved in detoxification, immunity, and metabolism, suggesting a complex regulatory network that encompasses immune responses and metabolic pathways. Functional enrichment and gene set enrichment analysis (GSEA) highlighted the activation of immune signaling pathways and the involvement of detoxification processes. Knockdown of CYP9E2 and CYP6K1 resulted in increased ROS levels, decreased detoxification enzyme activities, and reduced survival rates post-Bt challenge, implicating the critical role of these genes in immune priming and detoxification. Our findings demonstrate that LPS-induced immune priming in silkworms involves the upregulation of CYP450 genes, which play a critical role in detoxification and immune response modulation. The study provides insights into the molecular mechanisms of immune priming in insects and highlights the potential of CYP9E2 and CYP6K1 as targets for enhancing disease resistance and pest management in insects.

Transcriptomic and functional analysis of the antiviral response of mussels after a poly I:C stimulation

The study of mussels (Mytilus galloprovincialis) has grown in importance in recent years due to their high economic value and resistance to pathogens. Because of the biological characteristics revealed by mussel genome sequencing, this species is a valuable research model. The high genomic variability and diversity, particularly in immune genes, may be responsible for their resistance to pathogens found in seawater and continuously filtered and internalized by them. These facts, combined with the lack of proven mussel susceptibility to viruses in comparison to other bivalves such as oysters, result in a lack of studies on mussel antiviral response. We used RNA-seq to examine the genomic response of mussel hemocytes after they were exposed to poly I:C, simulating immune cell contact with viral dsRNA. Apoptosis and the molecular axis IRFs/STING-IFI44/IRGC1 were identified as the two main pathways in charge of the response but we also found a modulation of lncRNAs. Finally, in order to obtain new information about the response of mussels to putative natural challenges, we used VHSV virus (Viral Hemorrhagic Septicemia Virus) to run some functional analysis and confirm poly I:C's activity as an immunomodulator in a VHSV waterborne stimulation. Both, poly I:C as well as an injury stimulus (filtered sea water injection) accelerated the viral clearance by hemocytes and altered the expression of several immune genes, including IL-17, IRF1 and viperin.

Pathogen stimulations and immune cells synergistically affect the gene expression profile characteristics of porcine peripheral blood mononuclear cells

BACKGROUND

Pigs serve as a crucial source of protein in the human diet and play a fundamental role in ensuring food security. However, infectious diseases caused by bacteria or viruses are a major threat to effective global pig farming, jeopardizing human health. Peripheral blood mononuclear cells (PBMCs) are a mixture of immune cells that play crucial roles in immunity and disease resistance in pigs. Previous studies on the gene expression regulation patterns of PBMCs have concentrated on a single immune stimulus or immune cell subpopulation, which has limited our comprehensive understanding of the mechanisms of the pig immune response.

RESULTS

Here, we integrated and re-analyzed RNA-seq data published online for porcine PBMC stimulated by lipopolysaccharide (LPS), polyinosinic acid (PolyI:C), and various unknown microorganisms (EM). The results revealed that gene expression and its functional characterization are highly specific to the pathogen, identifying 603, 254, and 882 pathogen-specific genes and 38 shared genes, respectively. Notably, LPS and PolyI:C stimulation directly triggered inflammatory and immune-response pathways, while exposure to mixed microbes (EM) enhanced metabolic processes. These pathogen-specific genes were enriched in immune trait-associated quantitative trait loci (QTL) and eGenes in porcine immune tissues and were implicated in specific cell types. Furthermore, we discussed the roles of eQTLs rs3473322705 and rs1109431654 in regulating pathogen- and cell-specific genes CD300A and CD93, using cellular experiments. Additionally, by integrating genome-wide association studies datasets from 33 complex traits and diseases in humans, we found that pathogen-specific genes were significantly enriched for immune traits and metabolic diseases.

CONCLUSIONS

We systematically analyzed the gene expression profiles of the three stimulations and demonstrated pathogen-specific and cell-specific gene regulation across different stimulations in porcine PBMCs. These findings enhance our understanding of shared and distinct regulatory mechanisms of genetic variants in pig immune traits.

Physiological and transcriptomic analysis provides new insights into osmoregulation mechanism of Ruditapes philippinarum under low and high salinity stress

The Manila clam (Ruditapes philippinarum) is a commercially important marine bivalve, which inhabits the estuarine and mudflat areas. The osmoregulation is of great significance for molluscs adaptation to salinity fluctuations. In this study, we investigated the effects of low salinity (10 psu) and high salinity (40 psu) stress on survival and osmoregulation of the R. philippinarum. The results of physiological parameters showed that the ion (Na+, K+, Cl-) concentrations and Na+/K+-ATPase (NKA) activity of R. philippinarum decreased significantly under low salinity stress, but increased significantly under high salinity stress, indicating that there are differences in physiological adaptation of osmoregulation of R. philippinarum. In addition, we conducted the transcriptome analysis in the gills of R. philippinarum exposed to low (10 psu) and high (40 psu) salinity challenge for 48 h using RNA-seq technology. A total of 153 and 640 differentially expressed genes (DEGs) were identified in the low salinity (LS) group and high salinity (HS) group, respectively. The immune (IAP, TLR6, C1QL4, Ank3), ion transport (Slc34a2, SLC39A14), energy metabolism (PCK1, LDLRA, ACOX1) and DNA damage repair-related genes (Gadd45g, HSP70B2, GATA4) as well as FoxO, protein processing in endoplasmic reticulum and endocytosis pathways were involved in osmoregulation under low salinity stress of R. philippinarum. Conversely, the ion transport (SLC6A7, SLC6A9, SLC6A14, TRPM2), amino acid metabolism (GS, TauD, ABAT, ALDH4A1) and immune-related genes (MAP2K6, BIRC7A, CTSK, GVIN1), and amino acid metabolism pathways (beta-Alanine, Alanine, aspartate and glutamate, Glutathione) were involved in the process of osmoregulation under high salinity stress. The results obtained here revealed the difference of osmoregulation mechanism of R. philippinarum under low and high salinity stress through physiological and molecular levels. This study contributes to the assessment of salinity adaptation of bivalves in the context of climate change and provides useful information for marine resource conservation and aquaculture.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum)

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

Transcriptome analysis provides new insights into the immune response of Ruditapes philippinarum infected with Vibrio alginolyticus

Manila clam (Ruditapes philippinarum) is a bivalve species with commercial value, but it is easily infected by pathogenic microorganisms in aquaculture, which restricts the shellfish industry. Notably, the impact of Vibrio alginolyticus on clam culture is obvious. In this study, RNA-seq was performed to analyze clam hepatopancreas tissue in 48 h (challenge group, G48h) and 96 h (challenge group, G96h) after infection with V. alginolyticus and 0 h after injection of PBS (control group, C). The results showed that a total of 1670 differentially expressed genes were detected in the G48h vs C group, and 1427 differentially expressed genes were detected in the G96h vs C group. In addition, KEGG analysis showed that DEGs were significantly enriched in pathways such as Lysosome and Mitophagy. Moreover, 15 immune related DEGs were selected for qRT-PCR analysis to verify the accuracy of RNA-seq, and the results showed that the expression level of DEGs was consistent with that of RNA-seq. Therefore, the results obtained in this study provides a preliminary understanding of the immune defense of R. philippinarum and molecular insights for genetic breeding of V. alginolyticus resistance in Manila clam.

Investigation of Immune Responses in Giant African Snail, Achatina immaculata, against a Two-Round Lipopolysaccharide Challenge

As one of the 100 most-threatening invasive alien species, the giant African snail (Achatina immaculata) has successfully invaded and established itself in most areas of southern China. Protection against recurrent pathogen infections is vital to biological invasion. Enhanced immune protection has been previously found in other invertebrates, but not in the unique immune system of the giant African snail. In the present study, the survival rate of the giant African snail was recorded following a second infection with lethal doses of Escherichia coli after a previous first injection using lipopolysaccharide (LPS), and the mechanism of immune enhancement was investigated by examining the cellular and transcriptomic response of the giant African snail after two successive stimuli using LPS. Snails injected first with LPS, sterilized physiologic (0.9%) saline (SPS), phosphate-buffered saline (PBS) or untreated (Blank) were rechallenged at 7d with E. coli (Ec), and were named as LPS + Ec, SPS + Ec, PBS + Ec, Ec, and Blank. The log-rank test shows the survival rate of the LPS + Ec group as significantly higher than that of other control groups after the second injection (p < 0.05). By performing cell counting and BrdU labeling on newly generated circulating hemocytes, we found that the total hemocyte count (THC) and the ratio of BrdU-positive cells to total cells increased significantly after primary stimulation with LPS and that they further increased after the second challenge. Then, caspase-3 of apoptosis protease and two antioxidant enzyme activities (CAT and SOD) increased significantly after infection, and were significantly higher in the second response than they had been in the first round. Moreover, transcriptome analysis results showed that 84 differentially expressed genes (DEGs) were expressed at higher levels in both the resting and activating states after the second immune response compared to the levels observed after the first challenge. Among them, some DEGs, including Toll-like receptor 4 (TLR4) and its downstream signaling molecules, were verified using qRT-PCR and were consistent with the transcriptome assay results. Based on gene expression levels, we proposed that these genes related to the TLR signaling cascade participate in enhanced immune protection. All results provide evidence that enhanced immune protection exists in the giant African snail.

Transcriptome analysis of largemouth bass (Micropterus salmoides) challenged with LPS and polyI:C

Largemouth bass (Micropterus salmoides) is a worldwide commercially important aquatic species. In recent years, pathogenic diseases cause great economic losses and hinder the industry of largemouth bass. To further understand the immune response against pathogens in largemouth bass, splenic transcriptome libraries of largemouth bass were respectively constructed at 12 h post-challenged with phosphate-buffered saline (PBS), lipopolysaccharide (LPS) and polyinosinic-polycytidylic acid (polyI:C) by using RNA sequencing technology (RNA-seq). RNA libraries were constructed using 9 RNA splenic samples isolated from three biological replicates of the three groups and sequenced on the DNBSEQ platform. A total number of 86,306 unigenes were obtained. Through pairwise comparisons among the three groups, we identified 11,295 different expression genes (DEGs) exhibiting significant differences at the transcript level. There were 7, 7, and 13 signal pathways were significantly enriched in LPS-PBS comparison, polyI:C-PBS comparison, and LPS-polyI:C comparison, respectively, indicating that the immune response to different pathogens was distinct in largemouth bass. To the best of our knowledge, this is the first report on the immune response of largemouth bass against different pathogen-associated molecular patterns (PAMPs) stimuli using transcriptomic analysis. Our results provide a valuable resource and new insights to understanding the immune characteristics of largemouth bass against different pathogens.

Effects of environmental concentrations of the fragrance amyl salicylate on the mediterranean mussel Mytilus galloprovincialis

Amyl salicylate (AS) is a fragrance massively used as a personal care product and following the discharged in wastewaters may end up in the aquatic environment representing a potential threat for the ecosystem and living organisms. AS was recently detected in water of the Venice Lagoon, a vulnerable area continuously subjected to the income of anthropogenic chemicals. The lagoon is a relevant area for mollusc farming, including the Mediterranean mussels (Mytilus galloprovincialis) having an important economic and ecological role. Despite high levels of AS occurred in water of the Lagoon of Venice, no studies investigated the possible consequences of AS exposures on species inhabiting this ecosystem to date. For the first time, we applied a multidisciplinary approach to investigate the potential effects of the fragrance AS on Mediterranean mussels. To reach such a goal, bioaccumulation, cellular, biochemical, and molecular analyses (RNA-seq and microbiota characterization) were measured in mussels treated for 7 and 14 days with different AS Venice lagoon environmental levels (0.1 and 0.5 μg L^-1). Despite chemical investigations suggested low AS bioaccumulation capability, cellular and molecular analyses highlighted the disruption of several key cellular processes after the prolonged exposures to the high AS concentration. Among them, potential immunotoxicity and changes in transcriptional regulation of pathways involved in energy metabolism, stress response, apoptosis and cell death regulations have been observed. Conversely, exposure to the low AS concentration demonstrated weak transcriptional changes and transient increased representation of opportunistic pathogens, as Arcobacter genus and Vibrio aestuarianus. Summarizing, this study provides the first overview on the effects of AS on one of the most widely farmed mollusk species.

Apextrin from Ruditapes philippinarum functions as pattern recognition receptor and modulates NF-κB pathway

Apextrin belongs to ApeC-containing proteins (ACPs) and features a signal-peptide, an N-terminal membrane attack complex component/perforin (MACPF) domain, and a C-terminal ApeC domain. Recently, apextrin-like proteins were identified as pattern recognition receptor (PRR), which recognize the bacterial cell wall component and participate in innate immunity. Here, an apextrin (Rpape) was identified and characterized in Ruditapes philippinarum. Our results showed that Rpape mRNA was significantly induced under bacterial challenges. The Rpape recombinant protein exhibited a significant inhibitory effect on gram-positive bacteria (Bacillus subtilis and Staphylococcus aureus) and bound with Vibrio anguillarum, S. aureus and B. subtilis. We found Rpape protein positively activated the NF-κB signaling cascade and increased the activity of Nitric oxide (NO). This study revealed the immunity role of apextrin in R. philippinarum and provided a reference for further study on the role of apextrin in bivalves.

Gene Co-Expression Network Analysis Reveals the Correlation Patterns Among Genes in Different Temperature Stress Adaptation of Manila Clam

The Manila clam (Ruditapes philippinarum) is one of the most important aquaculture species and widely distributed along the coasts of China, Japan, and Korea. Due to its wide distribution, it can tolerate a wide range of temperature. Studying the gene expression profiles of clam gills had found differentially expressed genes (DEGs) and pathway involved in temperature stress tolerance. A systematic study of cellular response to temperature stress may provide insights into the mechanism of acquired tolerance. Here, weighted gene co-expression network analysis (WGCNA) was carried out using RNA-seq data from gill transcriptome in response to high and low temperature stress. There are a total 32 gene modules, of which 18 gene modules were identified as temperature-related modules. Blue module was one significantly correlated with temperature which was associated with cellular metabolism, apoptosis pathway, ER stress, and others.

Molecular Mechanisms Underlying Vibrio Tolerance in Ruditapes philippinarum Revealed by Comparative Transcriptome Profiling

The clam Ruditapes philippinarum is an important species in the marine aquaculture industry in China. However, in recent years, the aquaculture of R. philippinarum has been negatively impacted by various bacterial pathogens. In this study, the transcriptome libraries of R. philippinarum showing different levels of resistance to challenge with Vibrio anguillarum were constructed and RNA-seq was performed using the Illumina sequencing platform. Host immune factors were identified that responded to V. anguillarum infection, including C-type lectin domain, glutathione S-transferase 9, lysozyme, methyltransferase FkbM domain, heat shock 70 kDa protein, Ras-like GTP-binding protein RHO, C1q, F-box and BTB/POZ domain protein zf-C2H2. Ten genes were selected and verified by RT-qPCR, and nine of the gene expression results were consistent with those of RNA-seq. The lectin gene in the phagosome pathway was expressed at a significantly higher level after V. anguillarum infection, which might indicate the role of lectin in the immune response to V. anguillarum. Comparing the results from R. philippinarum resistant and nonresistant to V. anguillarum increases our understanding of the resistant genes and key pathways related to Vibrio challenge in this species. The results obtained here provide a reference for future immunological research focusing on the response of R. philippinarum to V. anguillarum infection.

Transcriptomic analysis provides insights into candidate genes and molecular pathways involved in growth of Manila clam Ruditapes philippinarum

Growth is one of the most important traits of aquaculture breeding programs. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Ruditapes philippinarum is an economically important marine bivalve. In order to gain insights into the molecular mechanisms to growth variability in marine shellfish, we conducted the transcriptome sequencing and examined the expression differences in growth-related gene and molecular pathways involved in growth trait of R. philippinarum. In this study, we investigated the molecular and gene expression differences in fast-growing and slow-growing Manila clam and focused on the analysis of the differential expression patterns of specific genes associated with growth by RNA-seq and qPCR analysis. A total of 61 differentially expressed genes (DEGs) were captured significantly differentially expressed, and were categorized into Ras signaling pathway, hedgehog signaling pathway, AMPK signaling pathway, p53 signaling pathway, regulation of actin cytoskeleton, focal adhesion, mTOR signaling pathway, VEGF signaling pathway, and TGF-beta signaling pathway. A total of 34 growth-related genes were validated significantly and up/downregulated at fast growing and slow growing of R. philippinarum. Functional enrichment analysis revealed the insulin signaling pathway, PI3K-Akt signaling pathway, and mTOR signaling pathway play pivotal roles in molecular function and regulation of growth trait in R. philippinarum. The growth-related genes and pathways obtained here provide important insights into the molecular basis of physiological acclimation, metabolic activity, and growth variability in marine bivalves.

Genome-wide investigation and expression analysis of MACPF gene family reveals its immune role in response to bacterial challenge of Manila clam

In this study, 18 MACPF genes (RpMACPF) were identified and classed into three types (Macrophage-expressed gene 1, Apextrin, and MACPF domain contain protein) based on gene structure and phylogenetic relationship in R. philippinarum. The length of RpMACPF proteins varied from 287 to 785 amino acids. The molecular weights and Theoretical PI values ranged from 3.2 kDa to 8.7 kDa and 4.7 to 8.6, respectively. RNA-seq data analysis revealed that 14 of 18 RpMACPF genes were highly expressed at the pediveliger larvae stage indicate RpMACPF might contribute to the early development and metamorphosis processes of the R. philippinarum. Besides, we found RpMACPF genes were significantly regulated by pathogen-associated molecular patterns (PAMPs) and Vibrio parahemolyticus, which indicates RpMACPF genes may play significant roles in response to invading pathogens. The results obtained in this work will provide valuable insight into the immune function of MACPF gene in R. philippinarum.

Comparative Proteomics of Ostreid Herpesvirus 1 and Pacific Oyster Interactions With Two Families Exhibiting Contrasted Susceptibility to Viral Infection

Massive mortality outbreaks affecting Pacific oysters (Crassostrea gigas) spat/juveniles are often associated with the detection of a herpesvirus called ostreid herpesvirus type 1 (OsHV-1). In this work, experimental infection trials of C. gigas spat with OsHV-1 were conducted using two contrasted Pacific oyster families for their susceptibility to viral infection. Live oysters were sampled at 12, 26, and 144 h post infection (hpi) to analyze host-pathogen interactions using comparative proteomics. Shotgun proteomics allowed the detection of seven viral proteins in infected oysters, some of them with potential immunomodulatoy functions. Viral proteins were mainly detected in susceptible oysters sampled at 26 hpi, which correlates with the mortality and viral load observed in this oyster family. Concerning the Pacific oyster proteome, more than 3,000 proteins were identified and contrasted proteomic responses were observed between infected A- and P-oysters, sampled at different post-injection times. Gene ontology (GO) and KEGG pathway enrichment analysis performed on significantly modulated proteins uncover the main immune processes (such as RNA interference, interferon-like pathway, antioxidant defense) which contribute to the defense and resistance of Pacific oysters to viral infection. In the more susceptible Pacific oysters, results suggest that OsHV-1 manipulate the molecular machinery of host immune response, in particular the autophagy system. This immunomodulation may lead to weakening and consecutively triggering death of Pacific oysters. The identification of several highly modulated and defense-related Pacific oyster proteins from the most resistant oysters supports the crucial role played by the innate immune system against OsHV-1 and the viral infection. Our results confirm the implication of proteins involved in an interferon-like pathway for efficient antiviral defenses and suggest that proteins involved in RNA interference process prevent viral replication in C. gigas. Overall, this study shows the interest of multi-omic approaches applied on groups of animals with differing sensitivities and provides novel insight into the interaction between Pacific oyster and OsHV-1 with key proteins involved in viral infection resistance.

Stimulation of Mytilus galloprovincialis Hemocytes With Different Immune Challenges Induces Differential Transcriptomic, miRNomic, and Functional Responses

Mediterranean mussels (Mytilus galloprovincialis) are marine bivalve molluscs with high resilience to biotic and abiotic stress. This resilience is one of the reasons why this species is such an interesting model for studying processes such as the immune response. In this work, we stimulated mussel hemocytes with poly I:C, β-glucans, and LPS and then sequenced hemocyte mRNAs (transcriptome) and microRNAs (miRNome) to investigate the molecular basis of the innate immune responses against these pathogen-associated molecular patterns (PAMPs). An immune transcriptome comprising 219,765 transcripts and an overview of the mussel miRNome based on 5,175,567 non-redundant miRNA reads were obtained. The expression analyses showed opposite results in the transcriptome and miRNome; LPS was the stimulus that triggered the highest transcriptomic response, with 648 differentially expressed genes (DEGs), while poly I:C was the stimulus that triggered the highest miRNA response, with 240 DE miRNAs. Our results reveal a powerful immune response to LPS as well as activation of certain immunometabolism- and ageing/senescence-related processes in response to all the immune challenges. Poly I:C exhibited powerful stimulating properties in mussels, since it triggered the highest miRNomic response and modulated important genes related to energy demand; these effects could be related to the stronger activation of these hemocytes (increased phagocytosis, increased NO synthesis, and increased velocity and accumulated distance). The transcriptome results suggest that after LPS stimulation, pathogen recognition, homeostasis and cell survival processes were activated, and phagocytosis was induced by LPS. β-glucans elicited a response related to cholesterol metabolism, which is important during the immune response, and it was the only stimulus that induced the synthesis of ROS. These results suggest a specific and distinct response of hemocytes to each stimulus from a transcriptomic, miRNomic, and functional point of view.

Molecular cloning and expression analysis of mc5r like genes (mc5rl) in Ruditapes philippinarum (Manila clam) after aerial exposure and low-temperature stress

The melanocortin-5 receptor (mc5r) plays an important role in exocrine function, lipid metabolism, obesity, and stress response in the vertebrate. However, the functions of the mc5r in mollusks have been rarely investigated. We cloned the full length of Ruditapes philippinarum mc5r like gene (mc5rl) and the sequence structure and phylogenetic relationship of mc5rl were analyzed. Besides, we detected the tissue distribution and the expression pattern of R. philippinarum mc5r like (mc5rl) genes after aerial exposure and low-temperature stress. The full-length cDNA of the mc5rl-1 was 2143 bp, consisting of a 1224 bp open reading frame encoding (ORF) 408 amino acids. Sequence and phylogenetic analyses revealed that the nucleotide and amino acid sequences of Manila clam mc5rl were highly homologous with mc5r of Crassostrea virginica, Crassostrea gigas, Mizuhopecten yessoensis, and Pecten maximus (32%-36%) and low homologous with vertebrates. The results of the distribution of mc5rl genes showed that mc5rl genes were dominant in the mantle, gonad, and hepatopancreas in R. philippinarum. The expression of mc5rl genes was significantly increased after aerial exposure and low-temperature stress in R. philippinarum in hepatopancreas. Aerial exposure and low-temperature stress could induce mc5rl expressed. Mc5rl might serve as a sensor and promote stress response in R. philippinarum. The cloning and expression characteristics of mc5rl will facilitate the investigation of its function in stress response and other physiological processes in R. philippinarum.