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

Transcriptome analysis reveals the DNA replication genes response to Vibrio anguillarum and NNV infection in Jinhu grouper (Epinephelus fuscoguttatus♀ × Epinephelus tukulal♂)

Vibrio anguillarum acts as an infectious agent in the aquaculture industry that causes a fatal hemolytic septicaemic disease in fish and shellfish. Viral nervous necrosis (VNN) disease seriously impacts the healthy development of the aquaculture industry. While the detrimental effects of V. anguillarum and NNV have been widely researched on freshwater and marine fish, whether they affect DNA replication in fish is unclear. In this study, we used Jinhu grouper as a model to investigate the influence of V. anguillarum and NNV on their DNA replication. By RNA-seq analysis in conjunction with gene set enrichment analysis, we found a significant upregulation of genes related to DNA replication in the liver and spleen of the Jinhu grouper after V. anguillarum infection while a prominent downregulation of genes related to DNA replication in the brain of the Jinhu grouper after NNV infection. We identified 27, 29 and 18 key genes involved in DNA replication that may respond to V. anguillarum and NNV infection and selected and identified six DNA replication genes (pole, pole2, pold1, pola1, pcna and mcm2). In addition, our results indicated that these genes exhibit tissue-specific expression patterns, with an increasing pattern of expression when V. anguillarum infected and a decreasing pattern when NNV infected. These results may provide valuable understanding on the underlying mechanisms of V. anguillarum and NNV infection in fish.

Integrated microRNA study and pathological analysis provides new insights into the immune response of Ruditapes philippinarum under Vibrio anguillarum challenge

Manila clam (Ruditapes philippenarum) is an important shellfish aquaculture product. The large-scale breeding of clams is often affected by V. anguillarum and causes large-scale death. However, the pathogenesis, immune response and metabolic pathway of V. anguillarum are still unclear. In this study, we found that the bacterial load in the hepatopancreas of R. philippinarum peaked at 48 h after V. anguillarum infection, and then gradually decreased, while the activity of lysozyme reached the peak at 12 h. Tissue section observation reveals that the infected hepatopancreas cells lost normal structure or necrosis. Additionally, six small RNA libraries were constructed using hepatopancreas of clams. A total of 15 differentially expressed (DE) microRNA (miRNA) were identified at 48 h after V. anguillarum infection, including 8 upregulated and 7 downregulated miRNAs. GO and KEGG enrichment results indicated the prediction of 48 known miRNAs and 127 new miRNAs, with functional annotation suggests that endocytosis pathway and bacterial recognition proteins may play key roles in immune response. The sequencing results were basically consistent with the qRT-PCR validation, indicating the accuracy of the data. This study provides a new idea to explore the immune regulation mechanism of shellfish after V. anguillarum infection, which brings important reference significance for modern immunological research.

Unravelling a Latent Pathobiome Across Coral Reef Biotopes

Previous studies on disease in coral reef organisms have neglected the natural distribution of potential pathogens and the genetic factors that underlie disease incidence. This study explores the intricate associations between hosts, microbial communities, putative pathogens, antibiotic resistance genes (ARGs) and virulence factors (VFs) across diverse coral reef biotopes. We observed a substantial compositional overlap of putative bacterial pathogens, VFs and ARGs across biotopes, consistent with the 'everything is everywhere, but the environment selects' hypothesis. However, flatworms and soft corals deviated from this pattern, harbouring the least diverse microbial communities and the lowest diversity of putative pathogens and ARGs. Notably, our study revealed a significant congruence between the distribution of putative pathogens, ARGs and microbial assemblages across different biotopes, suggesting an association between pathogen and ARG occurrence. This study sheds light on the existence of this latent pathobiome, the disturbance of which may contribute to disease onset in coral reef organisms.

Understanding the probiotic potential of Lactobacillus plantarum: Antioxidant capacity, non-specific immunity and intestinal microbiota improvement effects on Manila clam Ruditapes philippinarum

Lactic acid bacteria (LAB) have beneficial effects on aquatic animals, improving their immune system and intestinal microbiota. Nevertheless, the probiotic effects of LAB on the Manila clam Ruditapes philippinarum remain poorly understood. Herein, the effects of administering Lactobacillus plantarum at final doses of 1 × 105 CFU/L (T5 group), 1 × 107 CFU/L (T7 group), and 1 × 109 CFU/L (T9 group) in the rearing water for eight weeks were evaluated for the antioxidant capacity, non-specific immunity, resistance to Vibrio parahaemolyticus infection, and intestinal microbiota of R. philippinarum. The rearing water without the addition of L. plantarum served as a control. The results showed that the T7 and T9 groups demonstrated a significant elevation in the disease resistance of clams against V. parahaemolyticus, in the activities of alkaline phosphatase and lysozyme in the hepatopancreas, and in the expression of antioxidant- and immune-related genes, including SOD, GPx, and GST. Meanwhile, the T7 group showed a significant enhancement in superoxide dismutase and catalase activities and CAT expression, while the T9 group experienced a remarkable elevation in reduced glutathione content. Only catalase activity was markedly elevated in the T5 group. The expression of SOD, CAT, GPx, and GST was significantly elevated in three treatment groups following the V. parahaemolyticus challenge. The T7 group exhibited a significant increase in intestinal microbiota richness. Significant increases were noted in Firmicutes abundance across all three treatment groups and in Actinobacteriota in the T5 and T7 groups. Additionally, the opportunistic pathogen Escherichia-Shigella abundance significantly decreased in three treatment groups. Furthermore, administration of 1 × 107 CFU/L L. plantarum enhanced the stability of the intestinal microecosystem, whereas a dose of 1 × 109 CFU/L might have a negative effect. The application of three doses of L. plantarum significantly enhanced intestinal microbiota functions related to the immune response and oxidative stress regulation, while a higher dose (1 × 109 CFU/L) might inhibit several functions. In conclusion, the application of L. plantarum in the rearing water exerted beneficial effects on the antioxidant capacity, non-specific immunity, resistance to V. parahaemolyticus, and the intestinal microbiota stability and functions of R. philippinarum. The beneficial effects of L. plantarum on R. philippinarum were dose-dependent, and the final dose of 1 × 107 CFU/L exhibited the optimal effects.

Manipulated C5aR1 over/down-expression associates with IL-6 expression during bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis)

The complement component 5a/complement component 5 receptor 1 (C5a/C5aR1) pathway plays a crucial role in the onset and development of inflammation, but relevant studies in fish are lacking. In this study, we successfully characterized the relationship between half-smooth tongue sole (Cynoglossus semilaevis) C5aR1 (CsC5aR1) and bacterial inflammation. First, we showed that the overexpression of CsC5aR1 significantly increased bacterial pathological damage in the liver and intestine, whereas inhibition attenuated the damage. The in vitro experiments suggested that CsC5aR1 was able to positively regulate the phagocytic activity and respiratory burst of tongue sole macrophages. In terms of both transcriptional and translational levels, overexpression/inhibition of CsC5aR1 was followed by a highly consistent up-regulation/decrease of its downstream canonical inflammatory factor interleukin-6 (CsIL-6). Furthermore, we stimulated macrophages by lipopolysaccharide (LPS) and lipoteichoic acid (LTA) and found a broad-spectrum response to bacterial infections by the C5a/C5aR1 complement pathway together with the downstream inflammatory factor CsIL-6. Subsequently, we directly elucidated that CsIL-6 is an indicator of C5a/C5aR1-mediated inflammation at different infection concentrations, different infectious bacteria (Vibrio anguillarum and Mycobacterium marinum), and different detection levels. These results might provide a new inflammation bio-marker for early warning of bacteria-induced hyperinflammation leading to fish mortality and a promising target for the treatment of bacterial inflammation in teleost.

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.

A Full-Length Transcriptome and Analysis of the NHL-1 Gene Family in Neocaridina denticulata sinensis

Neocaridina denticulata sinensis has emerged as a promising model organism for basic studies in Decapod. However, the current transcriptome information on this species is based on next-generation sequencing technologies, which are limited by a short read length. Therefore, the present study aimed to generate a full-length transcriptome assembly of N. denticulata sinensis utilizing the PacBio Sequel Ⅱ platform. The resulting transcriptome assembly comprised 5831 transcripts with an N50 value of 3697 bp. Remarkably, 90.5% of these transcripts represented novel isoforms of known genes. The transcripts were further searched against the NR, SwissProt, KEGG, KOG, GO, NT, and Pfam databases. A total of 24.8% of the transcripts can be annotated across all seven databases. Additionally, 1236 alternative splicing events, 344 transcription factors, and 124 long non-coding RNAs (LncRNAs) were predicted. Based on the alternative splicing annotation results, a RING finger protein NHL-1 gene from N. denticulata sinensis (NdNHL-1) was identified. There are 15 transcripts in NdNHL-1. The longest transcript is 4995 bp in length and encodes a putative protein of 1665 amino acids. A phylogenetic analysis showed its close relationship with NHL-1 from other crustacean species. This report represents the full-length transcriptome of N. denticulata sinensis and will facilitate research on functional genomics and environmental adaptation in this species.

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.

Genome-Wide Identification and Characterization of MITF Genes in Ruditapes philippinarum and Their Involvement in the Immune Response to Vibrio anguillarum Infection

Studies have shown that the shellfish have innate immune system, which is a very important immune form of shellfish, and they rely on the innate immune system to resist diseases. As a transcription factor, Microphthalmia-associated transcription factor (MITF) plays a regulatory role in immune response and the shell color is also an important index for the breeding of excellent varieties of R. philippinarum. The research on immune response mechanism of RPMITFs can provide important reference data for the breeding of excellent clam varieties. In the genome of R. philippinarum, the RPMITF genes family of shell color-related gene family was selected as the target gene of this experiment. There are 12 RpMITF genes named RpMITF1, RpMITF2, RpMITF3, RpMITF4, RpMITF5, RpMITF6, RpMITF7, RpMITF8, RpMITF9, RpMITF10, RpMITF11, and RpMITF12. The open reading frame length is 639, 1233, 996, 1239, 675, 624, 816, 1365, 612, 1614, 1122, and 486 bp, encoding 212, 410, 331, 412, 224, 207, 271, 454, 203, 537, 373, and 161 aa, respectively. The predicted molecular weight range of amino acids is 18.85-62.61 kda, and the isoelectric point range is 5.26-9.44. Real-time quantitative PCR was used to detect the gene expression of RpMITF gene family in hepatopancreas tissues of two populations of Manila clam at 6 time points (0, 3, 6, 12, 24, and 48 h) after Vibrio anguillarum stress. The results show that RpMITF gene family was significantly expressed in hepatopancreas of two clam populations after V. anguillarum stress (P < 0.05).

Transcriptome Profile Analyses of Head Kidney in Roach (Rutilus rutilus), Common Bream (Abramis brama) and Their Hybrids: Does Infection by Monogenean Parasites in Freshwater Fish Reveal Differences in Fish Vigour among Parental Species and Their Hybrids?

Hybrid generations usually face either a heterosis advantage or a breakdown, that can be expressed by the level of parasite infection in hybrid hosts. Hybrids are less infected by parasites than parental species (especially F1 generations) or more infected than parental species (especially post-F1 generations). We performed the experiment with blood-feeding gill parasite Paradiplozoon homoion (Monogenea) infecting leuciscid species, Abramis brama and Rutilus rutilus, their F1 generation and two backcross generations. Backcross generations tended to be more parasitized than parental lines and the F1 generation. The number of differentially expressed genes (DEGs) was lower in F1 hybrids and higher in backcross hybrids when compared to each of the parental lines. The main groups of DEGs were shared among lines; however, A. brama and R. rutilus differed in some of the top gene ontology (GO) terms. DEG analyses revealed the role of heme binding and erythrocyte differentiation after infection by blood-feeding P. homoion. Two backcross generations shared some of the top GO terms, representing mostly downregulated genes associated with P. homoion infection. KEGG analysis revealed the importance of disease-associated pathways; the majority of them were shared by two backcross generations. Our study revealed the most pronounced DEGs associated with blood-feeding monogeneans in backcross hybrids, potentially (but not exclusively) explainable by hybrid breakdown. The lower DEGs reported in F1 hybrids being less parasitized than backcross hybrids is in line with the hybrid advantage.

Integrative analysis highlights molecular and immune responses of tick Amblyomma americanum to Escherichia coli challenge

Ticks are ectoparasites that can transmit various pathogens capable of causing life-threatening illnesses in people and animals, making them a severe public health threat. Understanding how ticks respond to bacterial infection is crucial for deciphering their immune defense mechanisms and identifying potential targets for controlling tick-borne diseases. In this study, an in-depth transcriptome analysis was used to investigate the molecular and immune responses of Amblyomma americanum to infection caused by the microinjection of Escherichia coli. With an abundance of differentially expressed genes discovered at different times, the analysis demonstrated significant changes in gene expression profiles in response to E. coli challenge. Notably, we found alterations in crucial immune markers, including the antimicrobial peptides defensin and microplusin, suggesting they may play an essential role in the innate immune response. Furthermore, KEGG analysis showed that following E. coli exposure, a number of key enzymes, including lysosomal alpha-glucosidase, fibroblast growth factor, legumain, apoptotic protease-activating factor, etc., were altered, impacting the activity of the lysosome, mitogen-activated protein kinase, antigen processing and presentation, bacterial invasion, apoptosis, and the Toll and immune deficiency pathways. In addition to the transcriptome analysis, we constructed protein interaction networks to elucidate the molecular interactions underlying the tick's response to E. coli challenge. Hub genes were identified, and their functional enrichment provided insights into the regulation of cytoskeleton rearrangement, apoptotic processes, and kinase activity that may occur in infected cells. Collectively, the findings shed light on the potential immune responses in A. americanum that control E. coli infection.

Metabolomic analysis provides new insights into the heat-hardening response of Manila clam (Ruditapes philippinarum) to high temperature stress

The temperature has always been a key environmental factor in Manila clam (Ruditapes philippinarum) culture. In this study, the Manila clam was treated to different temperature pre-heat (28 °C, 30 °C) and gained heat tolerance after recover of 12 h, and a survival rate (14.7 %-49.1 %) advantage after high temperature challenge (30 and 32 °C). To further investigate the physiological and metabolism changes in Manila clam that had experienced a heat stress, non-targeted metabolomics (LC-MS/MS) was used to analyze the metabolic responses of gills in three group Manila clams during the heat challenge. Metabolic profiles revealed that high temperature caused changes in fatty acid composition, energy metabolism, antioxidant metabolites, hydroxyl compounds, and amino acids in heat-hardened clams compared to non-hardened clams. We found a number of significantly enriched pathways, including cAMP signaling pathway, serotonergic synapse, and biosynthesis of unsaturated fatty acids in heat-hardened Manila clam compared with non-hardened and untreated Manila clam. After a brief high temperature treatment, the physiological maintenance ability of Manila clam was improved. Combined with metabolomics analysis, heat hardening treatment may improve the energy metabolism and antioxidant ability of Manila clam. These results provide new insights into the cellular and metabolic responses of Manila clams following high temperature stress.

Genome-wide identification and analysis of HECT E3 ubiquitin ligase gene family in Ruditapes philippinarum and their involvement in the response to heat stress and Vibrio anguillarum infection

E3 ubiquitin ligase (E3s) plays an important role in ubiquitin proteasome pathway, proteins containing homologous E6-AP carboxyl terminus (HECT) domains. However, the role of HECT E3 ubiquitin ligase in mollusk was rarely explored. In this study, we performed a genome-wide analysis of the HECT domain-containing gene in Ruditapes philippinarum to identify and predict the structural and functional characterization of HECT genes in response to abiotic and biotic stress. A total of sixteen members of HECT gene family were identified and analyzed for the gene structure, phylogenetic relation, three-dimensional structure, protein interaction network, and expression patterns. Experimental results demonstrated that Rph.HUWE1, Rph.HECTD1, Rph.Ubr5 were significantly up-regulated in response to heat stress and bacterial challenge. Taken together, our data provide insights into the potential function of HECT E3 ligase in heat stress and Vibrio anguillarum infection.

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.