Integrated transcriptomic and proteomic analyses reveal potential mechanisms linking thermal stress and depressed disease resistance in the turbot Scophthalmus maximus

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.

Integrated cell metabolomics and network pharmacology approach deciphers the anti-testosterone deficiency mechanisms of Bushen Zhuanggu Tang

Testicular dysfunction is distinguished by a deficiency in testosterone levels, which can be attributed to the occurrence of oxidative stress injury in Leydig cells. The empirical prescription known as Bushen Zhuanggu Tang, developed by a highly experienced traditional Chinese medicine practitioner with six decades of clinical expertize, aligns with the traditional Chinese medicine principle of "kidney governing bone". Researchers have demonstrated that the administration of BSZGT can effectively enhance testosterone production. The objective of this study is to investigate the potential anti-testicular dysfunction effects of BSZGT and elucidate its underlying mechanism in an in vitro setting. Specifically, the impact of oxidative stress induced by H2O2 on the activity and testosterone levels of Leydig cells (TM3) was examined. Furthermore, the utilization of UPLC-QE-Qrbitrap-MS enabled the identification of the involvement of BSZGT in various metabolic pathways, including arginine biosynthesis, amino acyl-tRNA biosynthesis, Alanine, aspartate and glutamine metabolism, and Citrate Cycle, through the modulation of 25 distinct metabolites. Additionally, a network pharmacological analysis was conducted to investigate the pivotal protein targets associated with the therapeutic effects of BSZGT. The findings demonstrate the identification of six key proteins (CYP19A1, CYP1B1, ALOX5, ARG1, XDH, and MPO) that play a significant role in augmenting testicular function through their involvement in the ovarian steroid production pathway. In summary, our study presents a comprehensive research methodology that combines cell metabonomics and network pharmacology to enhance the discovery of new therapeutic agents for TD.

pparβ regulates lipid catabolism by mediating acox and cpt-1 genes in Scophthalmus maximus under heat stress

Peroxisome proliferator-activated receptor β (pparβ) is a key gene-regulating lipid metabolism pathway, but its function in turbot remains unclear. In this study, the CDS of pparβ was cloned from kidney for the first time. The CDS sequence length was 1533 bp encoding 510 amino acids. Structural analysis showed that the pparβ protein contained a C4 zinc finger and HOLI domain, suggesting that the pparβ gene of turbot has high homology with the PPAR gene of other species. The high expression patterns of pparβ, acox, and cpt-1 at high temperatures, as shown through qPCR, indicated that high temperatures activated the transcriptional activity of pparβ and increased the activity of the acox and cpt-1 genes. The expression of acox and cpt-1 was significantly inhibited when pparβ was downregulated using RNAi technology and inhibitor treatments, suggesting that pparβ positively regulated acox and cpt-1 expression at high temperatures and, thus, modulates lipid catabolism activity. These results demonstrate that pparβ is involved in the regulation of lipid metabolism at high temperatures and expand a new perspective for studying the regulation of lipid metabolism in stress environments of teleost.

Whole Transcriptome Analysis of Differentially Expressed Genes in Cultured Nile Tilapia (O. niloticus) Subjected to Chronic Stress Reveals Signaling Pathways Associated with Depressed Growth

Chronic stress is a serious threat to aquaculture as it lowers fish growth performance and compromises fish welfare. The exact mechanism by which growth is retarded is, however, not clearly understood. This study sought to elucidate the gene expression profiles associated with chronic stress in cultured Nile tilapia (Oreochromis niloticus) reared for 70 days at different ammonia concentrations and stocking densities. Fish in the treatment groups showed negative growth, while the controls showed positive allometric growth. The specific condition factor (Kn) ranged from 1.17 for the controls to 0.93 for the ammonia and 0.91 for the stocking density treatments. RNA was extracted from muscle tissue using TRIzol followed by library construction and Illumina sequencing. Comparative transcriptome analysis revealed 209 differentially expressed genes (DEGs) (156 up- and 53 down-regulated) in the ammonia and 252 DEGs (175 up- and 77 down-regulated) in the stocking density treatment. In both treatments, 24 and 17 common DEGs were up- and down-regulated, respectively. DEGs were significantly enriched in six pathways associated with muscle activity, energy mobilization and immunity. The heightened muscular activity consumes energy which would otherwise have been utilized for growth. These results bring to fore the molecular mechanisms underlying chronic stress’ suppression of growth in cultured Nile tilapia.

Administration time-dependent effects of poly (I:C) on antioxidant and immune responses along the diurnal time scale in zebrafish

The circadian clock of vertebrates regulates many biological processes, including the immune system. This paper investigated whether responsiveness to poly (I:C), a synthetic analog of double-stranded RNA used as an immunostimulant, exhibits day/night differences in zebrafish. Fish were intraperitoneally (IP) injected with either phosphate-buffered saline (PBS) or poly (I:C) at two different time points: "Zeitgeber Time" (ZT) 4 (day) and ZT16 (night). Then, 6 h later, fish were euthanized, and tissue samples (skin, liver and kidney) were collected. A control group (intact fish) was also sampled at the same time points. The effect of poly (I:C) on the expression of antioxidant and immune genes was time-of-day-dependent, and the response was stronger following poly (I:C) administration in the day than at night. Time-dependent differences were observed for some genes in the PBS and control groups. However, these differences were tissue-specific. In liver, almost all the genes were affected by time of day. In kidney, poly (I:C) affected the expression of all the gene markers regardless of administration time. These findings highlight the importance of considering the time to administer poly (I:C) when evaluating the fish immune response.

Identification of rearing temperature-dependent host defense signaling against viral hemorrhagic septicemia virus infection

Viral hemorrhagic septicemia virus (VHSV) infection is associated with fatal outcomes in the aquaculture production of olive flounder (Paralichthys olivaceus). Olive flounders at low and high temperatures are known to be highly susceptible and resistant to VHSV infection, respectively. To study temperature-dependent innate immune activity, 4-aminobenzoic hydrazide (4-AH), a myeloperoxidase (MPO) inhibitor, was used to treat VHSV-infected olive flounders reared at a high temperature of 20 °C (20VI). Mortality, the MPO transcription, and the proteomic expression pattern of the 20VI group were then compared with those of groups of VHSV-infected flounders reared at 15 °C (15V) and 20 °C (20V). The cumulative mortality rate of the 20VI group was increased by 35% compared with that of the untreated 20V group. The MPO transcription was decreased 5.8-fold in 20VI than in 20V group. Its expression decreased further at a lower temperature and after exposure to VHSV. Histopathological analysis revealed necrosis of splenic tissue in 20VI and 15V, but not in 20V group. Based on clustering analysis, proteins with increased expression in 15V and 20VI groups were associated with viral mRNA translation and reproduction compared with those of 20V group. Increased expression of DHX58, MX1, and UBB was detected in 15V and 20VI groups, suggesting a role in triggering innate immune response. Unfortunately, these genes failed to induce the translocation of GLUT4 to the surface membrane from the intracellular location due to decreased expression of 14-3-3 proteins (YWHAB and YWHAZ) and microtubules (TUBA1A and TUBB4B). Suppression of glucose supply led to inactivation of MPO and suppression of MHC-I and MHC-II-linked immune activity, resulting in high viral infection and spread. In conclusion, this study highlights that defective GLUT4 translocation-dependent glucose uptake increases the mortality of VHSV-infected olive flounders by inhibiting MPO activity.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Transcriptome analysis reveals that high temperatures alter modes of lipid metabolism in juvenile turbot (Scophthalmus maximus) liver

With the increase in farming density and the continuously high summer temperatures against the background of global warming, high temperature stress has become a major challenge in fish farming. In this study, we simulated the high temperature environments (20 °C, 24 °C, and 28 °C) that may occur during turbot culture. High-throughput sequencing was used to investigate the lipid metabolism response patterns in juvenile turbot liver under high temperature stress. A total of 2067 differentially expressed genes (DEGs) were identified. KEGG analysis revealed that the DEGs were mainly associated with glycerophospholipid metabolism, steroid biosynthesis, glycerolipid metabolism, fatty acid metabolic pathways, and the PPAR signaling pathway. A regulatory network was constructed to further elucidate the transcriptional regulation of lipid metabolism. We speculated that high temperature activates PPAR signaling pathway through interaction with ligands such as fatty acids. On the one hand, the HMGCS1 gene in this pathway can inhibit sterol synthesis by down-regulating the expression of key genes in steroid biosynthesis pathway (SQLE, EBP, and DHCR24). On the other hand, the expression of ACSL1 in this pathway is significantly increased under high temperature, which may play an important role in regulating fatty acid metabolism. Moreover, we collected blood and detected changes in serum lipid parameters; the variation patterns were also consistent with our results. These findings reveal that lipid metabolism has an important regulatory role in stress resistance when turbot is exposed to high temperatures.

Multi-Omics Analysis Provides Novel Insight into Immuno-Physiological Pathways and Development of Thermal Resistance in Rainbow Trout Exposed to Acute Thermal Stress

In recent years, poikilothermic animals such as fish have increasingly been exposed to stressful high-temperature environments due to global warming. However, systemic changes in fish under thermal stress are not fully understood yet at both the transcriptome and proteome level. Therefore, the objective of this study was to investigate the immuno-physiological responses of fish under extreme thermal stress through integrated multi-omics analysis. Trout were exposed to acute thermal stress by raising water temperature from 15 to 25 °C within 30 min. Head-kidney and plasma samples were collected and used for RNA sequencing and two-dimensional gel electrophoresis. Gene enrichment analysis was performed: differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were identified to interpret the multi-omics results and identify the relevant biological processes through pathway analysis. Thousands of DEGs and 49 DEPs were identified in fish exposed to thermal stress. Most of these genes and proteins were highly linked to DNA replication, protein processing in the endoplasmic reticulum, cell signaling and structure, glycolysis activation, complement-associated hemolysis, processing of released free hemoglobin, and thrombosis and hypertension/vasoconstriction. Notably, we found that immune disorders mediated by the complement system may trigger hemolysis in thermally stressed fish, which could have serious consequences such as ferroptosis and thrombosis. However, antagonistic activities that decrease cell-free hemoglobin, heme, and iron might be involved in alleviating the side effects of thermally induced immuno-physiological disorders. These factors may represent the major thermal resistance traits that allow fish to overcome extreme thermal stress. Our findings, based on integration of multi-omics data from transcriptomics and proteomics analyses, provide novel insight into the pathogenesis of acute thermal stress and temperature-linked epizootics.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.