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

Environmental influence on single methylation variation sites (SMVs) in the large yellow croaker (Larimichthys crocea): identification and correlation analysis

BACKGROUND

Larimichthys crocea is an important aquaculture species along the southeastern coast of China, with diverse environment and farming practices since artificial breeding, these different aquatic habitats are subject to significant variations in environmental factors that may involve modulation of gene expression through epigenetic mechanisms to enable species to survive and reproduce.

METHODS AND RESULTS

This study aimed to identify methylation variation sites (SMVs) in different sequence contexts (CG, CHG, and CHH) within populations of L. crocea in different habitats. All SMV sites were subjected to linear regression with environmental factors to identify candidate genes involved environmental stress. The results indicate a significant correlation between SMV sites and various environmental factors. For the wild populations in Jinmen and Zhanjiang, the primary environmental pressures for adapting are temperature and salinity. In contrast, for the domesticated populations in Zhoushan and farmed population in Xiangshan, the main environmental pressures are nitrate and dissolved oxygen. Furthermore, genes related to temperature adaptation in different aquatic environments were identified, including nr3c2, igf1, hsp70, trpm3, and fgf1. The gene rasa3 was found to be associated with pH adaptation, while genes such as atp6ap1lb, slc15a4, and gpr39 were linked to salinity, ammonia nitrogen, and dissolved oxygen. Research on the association between single methylation variation sites (SMVs) and environmental factors in aquatic organisms is scarce.

CONCLUSIONS

These results suggest that selection pressures can influence a significant proportion of methylation sites in this species, indirectly implying that epigenetic variation is not solely attributed to patterns of genetic variation, but is also closely linked to environmental differences. These results highlight the complex interactions between epigenetic regulation and environmental influences. Hence, this study provides preliminary evidence for a new perspective on the role of methylation patterns in L. crocea in environmental adaptation.

Identification and function analysis of Toll-like receptor 4 (TLR4) from Manila clam (Ruditapes philippinarum)

Toll-like receptor 4 (TLR4) is a pattern recognition receptor that activates innate immunity in response to pathogen infection. However, the role of TLR4 in pathogen-induced apoptosis and host immunity in mollusks remains largely unknown. In this study, the TLR4 of the Manila clam Ruditapes philippinarum (RpTLR4) was cloned. The open reading frame of RpTLR4 encodes a protein of 734 amino acids, containing a conserved TIR domain. Phylogenetic analysis revealed that RpTLR4 clusters closely with TLR4s from mollusks. RpTLR4 mRNA was detected in all tested tissues, with notably high expression in hemocytes (428-fold) and gills (657-fold). Subcellular localization showed that RpTLR4 is expressed on the cell membrane. qRT-PCR and western blot analyses demonstrated that RpTLR4 expression was induced in Manila clams after treatment with Vibrio parahaemolyticus. Overexpression of RpTLR4 significantly increased apoptosis levels and the expression of apoptosis-related genes. Conversely, silencing RpTLR4 markedly reduced the apoptosis rate in hemocytes induced by V. parahaemolyticus, indicating that V. parahaemolyticus-induced hemocyte apoptosis depends on RpTLR4 expression. Overall, these findings confirm that RpTLR4 plays a pro-apoptotic role in the response of Manila clams to V. parahaemolyticus infection. This study provides a theoretical foundation for understanding the molecular mechanisms underlying mollusk responses to pathogen infection.

Transcriptome-based network analysis reveals hub immune genes and pathways of hepatopancreas against LPS in Amphioctopus fangsiao

The hepatopancreas is the biggest digestive organ in Amphioctopus fangsiao (A. fangsiao), but also undertakes critical functions like detoxification and immune defense. Generally, pathogenic bacteria or endotoxin from the gut microbiota would be arrested and detoxified in the hepatopancreas, which could be accompanied by the inevitable immune responses. In recent years, studies related to cephalopods immune have been increasing, but the molecular mechanisms associated with the hepatopancreatic immunity are still unclear. In this study, lipopolysaccharide (LPS), a major component of the cell wall of Gram-negative bacteria, was used for imitating bacteria infection to stimulate the hepatopancreas of A. fangsiao. To investigate the immune process happened in A. fangsiao hepatopancreas, we performed transcriptome analysis of hepatopancreas tissue after LPS injection, and identified 2615 and 1943 differentially expressed genes (DEGs) at 6 and 24 h post-injection, respectively. GO and KEGG enrichment analysis showed that these DEGs were mainly involved in immune-related biological processes and signaling pathways, including ECM-receptor interaction signaling pathway, Phagosome signaling pathway, Lysosome signaling pathway, and JAK-STAT signaling pathways. The function relationships between these DEGs were further analyzed through protein-protein interaction (PPI) networks. It was found that Mtor, Mapk14 and Atm were the three top interacting DEGs under LPS stimulation. Finally, 15 hub genes involving multiple KEGG signaling pathways and PPI relationships were selected for qRT-PCR validation. In this study, for the first time we explored the molecular mechanisms associated with hepatopancreatic immunity in A. fangsiao using a PPI networks approach, and provided new insights for understanding hepatopancreatic immunity in A. fangsiao.

Integrated evidence of transcriptional, metabolic, and intestinal microbiota changes in Ruditapes philippinarum due to perfluorooctanoic acid-induced immunotoxicity

Perfluorooctanoic acid (PFOA) is a toxic pollutant that bioaccumulates and is a significant public health concern due to its ubiquitous and persistent occurrence in global environments. Few studies have evaluated the adverse effects of PFOA on immune system, and this is particularly true for mollusks. Here, the PFOA-associated effects on immune system were evaluated in Ruditapes philippinarum using integrated analysis of metabolomes, microbiomes, and transcriptomes, providing evidence for possible mechanisms related to immunotoxicity. PFOA exposure caused clear variation in several important metabolites related to immune regulatory function within the haemolyph from R. philippinarum, while also altering key metabolic pathways, including those of lipids, unsaturated fatty acids (UFAs), and bile acids (BAs). After exposure to PFOAs, intestinal bacterial communities also clearly changed, with the predominant microflora becoming Mycoplasma and Bacteroidetes that are related to intestinal inflammation. Molecular analyses provided consistent results, wherein the expression of immune-related genes was significantly altered. Integration of the multi-'omics' analyses suggested that the TLR/MyD88/NF-kB pathway, along with PI3K-Akt-mTOR pathway, PPAR-mediated lipid metabolism and the autophagy signaling pathway, likely play important roles in initiating immunotoxic effects in R. philippinarum after PFOA exposure. These results provide further evidence that PFOA exposure can lead to immunologic dysfunction and also provide new insights into the mechanisms of PFAS alteration of bivalve immune function.

MicroRNA transcriptome analysis reveals the immune regulatory mechanism of Crassostrea hongkongesis against Vibrio harveyi infection

MicroRNAs (miRNAs) are small non-coding RNA molecules that modulate target-genes expression and play crucial roles in post-transcriptional regulation and immune system regulation. The Hong Kong oyster (Crassostrea hongkongesis), as the main marine aquaculture shellfish in the South China Sea, not only has high economic and ecological value, but also is an ideal model for conducting research on pathogen host interaction. Vibrio harveyi, a Gram negative luminescent marine bacterium, is widely distributed in coastal water environments and can cause large-scale death of C. hongkongesis. However, little in formation is available on the immune regulatory mechanisms of C. hongkongesis infected with V. harveyi. Therefore, we performed microRNA transcriptome analysis for elucidating the immunoregulation mechanism of C. hongkongesis infected with V. harveyi. The results show that a total of 308468208 clean reads and 288371159 clean tags were obtained. 222 differentially expressed miRNAs were identified. A total of 388 target genes that were differentially expressed and negatively correlated with miRNA expression were predicted by 222 DEmiRs. GO enrichment analysis of 388 DETGs showed that they were mainly enriched in the immune-related term of membrane-bounded vesicle, endocytic vesicle lumen, antigen processing and presentation of exogenous peptide antigen via MHC class I, antigen processing and presentation of peptide antigen via MHC class I, and other immune-related term. KEGG enrichment analysis showed that DETGs were mainly enriched in the Complement and coagulation cascades, Herpes simplex virus 1 infection, Bacterial invasion of epithelial cells, Antigen processing and presentation and NOD-like receptor signaling pathway. The 16 key DEmiRs and their target genes form a regulatory network for seven immune-related pathways. These results suggest that V. harveyi infection induces a complex miRNA response with wide-ranging effects on immune gene expression in the C. hongkongesis. This study explored the immune response of C. hongkongesis to V. harveyi infection at the level of miRNAs, which provides new ideas for the healthy culture and selective breeding of C. hongkongesis.

Comparative Transcriptome Analysis of the Response to Vibrio parahaemolyticus and Low-Salinity Stress in the Swimming Crab Portunus trituberculatus

Vibrio parahaemolyticus is one of the main pathogenic bacteria of Portunus trituberculatus and causes mass mortality of P. trituberculatus in aquaculture. In addition, low-salinity stimulation makes P. trituberculatus more susceptible to V. parahaemolyticus infections. In order to elucidate the molecular mechanism of resistance to V. parahaemolyticus in P. trituberculatus, comparative transcriptomic analysis of blood cells stimulated by low salinity and V. parahaemolyticus was carried out in this study. Transcriptome sequencing of low-salinity stress and pathogen infection at different time points was completed using Illumina sequencing technology. A total of 5827, 6432, 5362 and 1784 differentially expressed genes (DEGs) involved in pathways related to ion transport and immunoregulation were found under low-salinity stress at 12, 24, 48 and 72 h compared with the control at 0 h. In contrast, 4854, 4814, 5535 and 6051 DEGs, which were significantly enriched in Toll and IMD signaling pathways, were found at 12, 24, 48 and 72 h compared with the control at 0 h under V. parahaemolyticus infection. Among them, 952 DEGs were shared in the two treatment groups, which were mainly involved in apoptosis and Hippo signaling pathway. Cluster analysis screened 103 genes that were differentially expressed in two factors that were negatively correlated, including immunoglobulin, leukocyte receptor cluster family, scavenger receptor, macroglobulin and other innate-immune-related genes. These results provide data support for the analysis of the mechanisms of immunity to V. parahaemolyticus under low-salinity stress in P. trituberculatus and help to elucidate the molecular mechanisms by which environmental factors affect immunity.

Transient receptor potential (TRP) channels in the Manila clam (Ruditapes philippinarum): Characterization and expression patterns of the TRP gene family under heat stress in Manila clams based on genome-wide identification

In this study, we identified a total of 40 transient receptor potential genes (RpTRP) in Manila clam by genome-wide identification and classified them into four categories (TRPV, TRPA, TRPM, TRPC) based on gene structure and subfamily relationships. The protein length of RpTRP genes ranges from 281 amino acids to 1601 amino acids. Molecular weight and theoretical PI values range from 182.82 kDa to 32.43 kDa, respectively, with PI values between 5.17 and 9.25. By comparing the expression profiles of TRP genes during heat stress in Manila clams at different latitudes, we found that most genes in the TRP gene family were up-regulated in expression during heat challenge. Therefore, we determined that TRP genes have an important role in the heat stress of Manila clams. This work provides a basis for further studies on the molecular mechanisms of TRP-mediated heat tolerance in Manila clam and for explaining differences in heat tolerance in Manila clam at different latitudes through key differential TRP genes at the molecular level.

Comparative transcriptome analysis reveals immunotoxicology induced by three organic UV filters in Manila clam (Ruditapes philippinarum)

Benzophenone-3 (BP-3), 4-methyl-benzylidene camphor (4-MBC) and 2-ethyl-hexyl-4-trimethoxycinnamate (EHMC) are commonly used organic ultraviolet (UV) filters and are frequently detected in water environments. In the present study, we studied the potential adverse impacts of UV filter exposures in Ruditapes philippinarum by investigating transcriptomic profiles and non-specific immune enzyme activities. Transcriptome analysis showed that more genes were differentially regulated in EHMC-treated group, and down-regulated genes (2009) were significantly more than up-regulated ones (410) at day 7. Function annotation revealed that pathways "immune system", "cell growth and death" and "infectious diseases" were significantly enriched. Generally, combined qPCR and biochemical analyses demonstrated that short-term exposure to low dose of UV filters could activate immune responses, whereas the immune system would be restrained after prolonged exposure. Taken together, the present study firstly demonstrated the immunotoxicology induced by BP-3, 4-MBC and EHMC on R. philippinarum, indicating their potential threats to the survival of marine bivalves.

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.

Immunotoxicity of Perfluorooctanoic Acid to the Marine Bivalve Species Ruditapes philippinarum

Polyfluorinated alkylated substances are recognized as an important class of pollutants in marine environments. Bivalves are good model organisms for evaluating the toxicity of pollutants and monitoring marine environments. In the present study, immunotoxicity of perfluorooctanoic acid (PFOA) was investigated by measuring biomarkers of the immune profile of Ruditapes philippinarum. In bivalves, hemocytes are an important component of the immune system. Thus, hemocyte proliferation, phagocytosis, cell viability, and immune enzyme activities, which have been applied as marine pollution bioindicators, were identified and observed for changes after exposure to PFOA in R. philippinarum. Based on the integrated biomarker responses method, we selected five biomarkers to evaluate PFOA risk at the multibiomarker level. In addition, the histopathological alterations of hemocytes in bivalves were used as indexes of the response to environmental stress. The subcellular structure of the hemocytes in R. philippinarum changed significantly with PFOA exposure, including hemocyte and nucleus morphological changes, organelle dissolution, cytomembrane and karyotheca swelling, and cytoplasm vacuolization. The present study verifies PFOA immunotoxicity to R. philippinarum at different levels and the integrated assessment of stress levels caused by PFOA in marine environment. Our results will provide new insights into evaluating adverse effects of PFOA and monitoring marine ecosystem. Environ Toxicol Chem 2022;41:426-436. © 2021 SETAC.

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.