Integrated application of transcriptomics and metabolomics provides insights into unsynchronized growth in pearl oyster Pinctada fucata martensii

Metabolic perturbation and oxidative damage induced by tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) and tris(2-ethylhexyl) phosphate (TEHP) on Escherichia coli through integrative analyses of metabolome

Organophosphate esters (OPEs) are one of the emerging environmental threats, causing the hazard to ecosystem safety and human health. Yet, the toxic effects and metabolic response mechanism after Escherichia coli (E.coli) exposed to TDCIPP and TEHP is inconclusive. Herein, the levels of SOD and CAT were elevated in a concentration-dependent manner, accompanied with the increase of MDA contents, signifying the activation of antioxidant response and occurrence of lipid peroxidation. Oxidative damage mediated by excessive accumulation of ROS decreased membrane potential and inhibited membrane protein synthesis, causing membrane protein dysfunction. Integrative analyses of GC-MS and LC-MS based metabolomics evinced that significant perturbation to the carbohydrate metabolism, nucleotide metabolism, lipids metabolism, amino acid metabolism, organic acids metabolism were induced following exposure to TDCIPP and TEHP in E.coli, resulting in metabolic reprogramming. Additionally, metabolites including PE(16:1(5Z)/15:0), PA(17:0/15:1(9Z)), PC(20:2(11Z,14Z)/12:0), LysoPC(18:3(6Z,9Z,12Z)/0:0) were significantly upregulated, manifesting that cell membrane protective molecule was afforded by these differential metabolites to improve permeability and fluidity. Overall, current findings generate new insights into the molecular toxicity mechanism by which E.coli respond to TDCIPP and TEHP stress and supply valuable information for potential ecological risks of OPEs on aquatic ecosystems.

Effects of ammonia nitrogen stress on the physiological, biochemical, and metabolic levels of the gill tissue of juvenile four-finger threadfin (Eleutheronema tetradactylum)

In this study, the impact of ammonia nitrogen stress on juvenile four-finger threadfin in pond culture was examined. The 96-hour median lethal concentration (LC50-96h) and safe concentration of ammonia nitrogen were assessed in juveniles with a body weight of 7.4 ± 0.6 g using ecotoxicological methods. The study design included a stress group exposed to LC50-96h levels of ammonia nitrogen and a control group without ammonia nitrogen exposure. To examine the physiological, biochemical, and metabolic effects of ammonia nitrogen on gill tissue, gill tissue samples were collected after 12, 24, 48, and 96 h of stress, with a resumption of treatment after 48 h. Compared to the control group, ammonia nitrogen adversely affected juvenile four-finger threadfin, with LC50-96h and safe concentration values of 20.70 mg/L and 2.07 mg/L, respectively. Exposure to ammonia nitrogen resulted in substantial gill damage, including fusion of lamellae, epithelial cell loss, and proliferation of chlorine-secreting cells. This tissue damage persisted even after a 48-h recovery period. Ammonia nitrogen stress triggered an increase in antioxidant enzyme activity (superoxide dismutase, catalase, and glutathione peroxidase) and malondialdehyde levels in gills, indicating oxidative stress from 12 h onwards. Although enzyme activity decreased over time, oxidative stress persisted even after recovery, suggesting an ongoing need for antioxidant defense. Metabolomics analysis showed significant alterations in 423 metabolites under ammonia nitrogen stress. Key metabolites such as L-arginine, taurine, 20-hydroxyarachidonic acid, 11,12-dihydroxy-5Z, 8Z, and 14Z eicosotrienic acid followed an increasing trend; uridine, adenosine, L-glutathione, and thymidine 5'-triphosphate followed a decreasing trend. These changes reflect metabolic adaptations to stress. In enriched metabolic pathways, the main differential pathways are membrane transport, lipid metabolism, and amino acid metabolism. After 48 h, significant differences were observed in 396 metabolites compared to the control group. Notably, L-arginine, choline, and L-histidine increased, while linoleic acid, adenosine, and glutathione decreased. Amino acid and lipid metabolism pathways were key affected pathways. Under ammonia nitrogen stress, juvenile four-finger threadfin increased the synthesis of unsaturated and saturated fatty acids to cope with low temperatures and bolster immune function by consuming spermidine. This adaptation helps to clear peroxides generated during fatty acid synthesis, thereby protecting cells from oxidative damage. This study provides insights for pond aquaculture and breeding of ammonia nitrogen-tolerant fish strains.

Comparative transcriptomics identifies genes underlying growth performance of the Pacific black-lipped pearl oyster Pinctada margaritifera

BACKGROUND

In bivalves, the rate at which organisms grow is a major functional trait underlying many aspects of their commercial production. Growth is a highly polygenic trait, which is typically regulated by many genes with small to moderate effects. Due to its complexity, growth variability in such shellfish remains poorly understood. In this study, we aimed to investigate differential gene expression among spat of the pearl oyster Pinctada margaritifera with distinct growth phenotypes.

RESULTS

We selected two groups of P. margaritifera spat belonging to the same F2 cohort based on their growth performance at 5.5 months old. Transcriptome profile analysis identified a total of 394 differentially expressed genes between these Fast-growing (F) and Slow-growing (S) phenotypes. According to functional enrichment analysis, S oysters overexpressed genes associated with stress-pathways and regulation of innate immune responses. In contrast, F oysters up-regulated genes associated with cytoskeleton activity, cell proliferation, and apoptosis. Analysis of genome polymorphism identified 16 single nucleotide polymorphisms (SNPs) significantly associated with the growth phenotypes. SNP effect categorization revealed one SNP identified for high effect and annotated for a stop codon gained mutation. Interestingly, this SNP is located within a gene annotated for scavenger receptor class F member 1 (SRF1), which is known to modulate apoptosis. Our analyses also revealed that all F oysters showed up-regulation for this gene and were homozygous for the stop-codon mutation. Conversely, S oysters had a heterozygous genotype and a reduced expression of this gene.

CONCLUSIONS

Altogether, our findings suggest that differences in growth among the same oyster cohort may be explained by contrasted metabolic allocation between regulatory pathways for growth and the immune system. This study provides a valuable contribution towards our understanding of the molecular components associated with growth performance in the pearl oyster P. margaritifera and bivalves in general.

Integrated transcriptomic and metabolomic analysis reveals the response of pearl oyster (Pinctada fucata martensii) to long-term hypoxia

The frequency at which organisms are exposed to hypoxic conditions in aquatic environments is increasing due to coastal eutrophication and global warming. To reveal the effects of long-term hypoxic stress on metabolic changes of pearl oyster, commonly known as Pinctada (Pinctada fucata martensii), the present study performed the integrated analysis of transcriptomics and metabolomics to investigate the global changes of genes and metabolites following 25 days hypoxia challenge. Transcriptome analysis detected 1108 differentially expressed genes (DEGs) between the control group and the hypoxia group. The gene ontology (GO) analysis of DEGs revealed that they are significantly enriched in functions such as "microtubule-based process", "histone (H3-K4, H3-K27, and H4-K20) trimethylation", "histone H4 acetylation", "kinesin complex", and "ATPase activity", and KEGG pathway functions, such as "DNA replication", "Apoptosis", and "MAPK signaling pathways". Metabolome analysis identified 68 significantly different metabolites from all identified metabolites, and associated with 25 metabolic pathways between the control and hypoxia groups. These pathways included aminoacyl-tRNA biosynthesis, arginine and proline metabolism, and phenylalanine metabolism. Our integrated analysis suggested that pearl oysters were subject to oxidative stress, apoptosis, immune inhibition, and neuronal excitability reduction under long-term hypoxic conditions. We also found a remarkable depression in a variety of biological functions under long-term hypoxia, including metabolic rates, biomineralization activities, and the repression of reorganization of the cytoskeleton and cell metabolism. These findings provide a basis for elucidating the mechanisms used by marine bivalves to cope with long-term hypoxic stress.

Comprehensive profiling of the metabolome in corn silage inoculated with or without Lactiplantibacillus plantarum using different untargeted metabolomics analyses

Silage fermentation is a complicated biochemical process involving interactions between microbes and metabolites. However, the overall metabolome feature of ensiled forage and its response to lactic acid bacteria inoculation is poorly understood. Hence, in this study metabolome profiles of whole-plant corn silage inoculated with or without Lactiplantibacillus plantarum were characterised via solid-phase microextraction/gas chromatography/mass spectrometry (SPME-GC-MS), gas chromatography/time-of-flight mass spectrometry (GC-TOF-MS), and Liquid chromatography/Q Exactive HFX mass spectrometry (LC-QE-MS/MS) analysis. There were 2087 identified metabolites including 1143 reliably identified metabolites in fresh and ensiled whole-plant corn. After ensiling, the increased metabolites in whole-plant corn were mainly composed of organic acids, volatile organic compounds (VOC), benzene and substituted derivatives, carboxylic acids and derivatives, fatty acyls, flavonoids, indoles and derivatives, organooxygen compounds (including amines and amides), phenols, pyridines and derivatives, and steroids and steroid derivatives, which includes neurotransmitters and metabolites with aromatic, antioxidant, anti-inflammatory, and antimicrobial activities. Phenylacetaldehyde was the most abundant aromatic metabolite after ensiling. L-isoleucine and oxoproline were the major free amino acids in silage. Ensiling markedly increased the relative abundances of 3-phenyllactic acid, chrysoeriol, 6-O-acetylaustroinulin, acetylcholine, γ-aminobutyric acid, pyridoxine, and alpha-linoleic acid. Inoculation with L. plantarum remarkably changed silage VOC composition, and essential amino acids, 3-phenyllactic acid, and cinnamaldehyde compared with untreated silage. The present study does not only provide a deeper insight into metabolites of the ensiled whole-plant corn but also reveals metabolites with specific biological functions that could be much helpful in screening novel lactic acid bacteria to well ensile forages. Inoculation with L. plantarum significantly affects the metabolome in ensiled whole-plant corn.

Metabolomic and transcriptomic profiling reveals the effect of dietary protein and lipid levels on growth performance in loach (Paramisgurnus dabryanus)

The subject of this study was to explore the optimum requirements of loach (Paramisgurnus dabryanus) regarding dietary proteins and lipids and discuss the underlying mechanism. We designed nine diets to determine the effects of different levels of dietary crude protein (CP: 30%, 35%, and 40%) and ether extract (EE: 6%, 10%, and 14%) on the growth performance and metabolism of P. dabryanus. In total, 2160 healthy P. dabryanus (5.19 ± 0.01 g) were divided into nine groups with four replications at 60 fish per barrel stocking density. The trial lasted for eight weeks. Serum and liver samples were gathered for metabolomic and transcriptomic analyses. The results showed that the specific growth rate of P. dabryanus in the CP40EE10 group was the fastest and notably higher than that in other groups (P< 0.05). Analysis of the metabolome results found that the mTOR signaling pathway, glycerophospholipid metabolism, D-arginine and D-ornithine metabolism were significantly enriched pathways in the CP40EE10 group compared with the other groups (P< 0.05). Moreover, the transcriptomic analysis of differentially expressed genes (DEGs) showed that the expression of ARG (arginase) involved in protein synthesis was significantly upregulated in the CP40EE10 group compared to the slowest growing group (P< 0.05). Additionally, the expression of SPLA2 (secretory phospholipase A2) involved in lipid metabolism and FBP (fructose-1,6-bisphosphatase) involved in glucose metabolism were all significantly downregulated in the CP30EE6 group compared with the CP40EE10 group (P< 0.05). Furthermore, the analysis of differentially expressed metabolites (DEMs) and DEGs co-enriched in the KEGG pathway revealed that the significantly enriched pathways were arginine and proline metabolism, glycerophospholipid metabolism, and glycolysis/gluconeogenesis in CP40EE10 compared with other groups (P< 0.05). We conclude that including 40% CP and 10% EE in the P. dabryanus diet could result in a better growth rate. We hypothesized from metabolomic and transcriptomic analyses that the CP40EE10 diet might promote the growth of P. dabryanus by promoting protein synthesis, lipid metabolism, and energy production.

Analysis of Volatile Compounds and Flavor Fingerprint Using Gas Chromatography-Ion Mobility Spectrometry (GC-IMS) on Crassostrea gigas with Different Ploidy and Gender

In this study, GC-IMS was used to analyze the volatile component and flavor profiles of Crassostrea gigas individuals of different ploidy and gender. Principal component analysis was used to explore overall differences in flavor profiles, and a total of 54 volatile compounds were identified. The total volatile flavor contents in the edible parts of tetraploid oysters were significantly higher than in diploid and triploid oysters. The concentrations of ethyl (E)-2-butenoate and 1-penten-3-ol were significantly higher in triploid oysters than in diploid and tetraploid oysters. In addition, the volatile compounds propanoic acid, ethyl propanoate, 1-butanol, butanal, and 2-ethyl furan were significantly higher in females than in males. The volatile compounds p-methyl anisole, 3-octanone, 3-octanone, and (E)-2-heptenal were present in higher levels in male than in female oysters. Overall, different ploidy and gender of oysters are connected with different sensory characteristics, providing new insights for understanding the flavor characteristics of oysters.

Transcriptome analysis of growth and shell color between two genetic variants of Corbicula fluminea with different shell colors

To clarify the molecular mechanism of the black and yellow shell coloration, we performed a transcriptome analysis of whole tissue of Corbicula fluminea in Hongze Lake (Jiangsu Province, China). After assembly, 335,247 unigenes were obtained, and 136,804 unigenes were functionally identified using public databases (NR, GO, KEGG, eggnog, and Swissprot). 1567 differentially expressed genes (DEGs) were detected through pairwise comparisons, of which 941 DEGs were up-regulated and 626 were down-regulated in the black-shelled clam. We compared the DEGs between two clams and identified some coloration-related genes. Notably, the black-shelled clam was larger than the yellow-shelled. We speculated that higher digestion and anabolic ability of black-shelled clam might lead to this phenomenon. In contrast, the yellow-shelled clam appeared to be more sensitive to environmental stress. The metabolic energy of the yellow-shelled clam was depleted to maintain or recover from stress, and provide less energy for growth. In summary, our finding provides a theoretical basis for the molecular mechanism of pigmentation and the difference of somatotype in bivalve, as well as promotes the future breeding of more elite varieties.

Integrated transcriptomic and metabolomic analysis sheds new light on adaptation of Pinctada fucata martensii to short-term hypoxic stress

Analyses of the transcriptome and metabolome were conducted to clarify alterations of key genes and metabolites in pearl oysters following exposure to short-term hypoxic treatment. We totally detected 209 DEGs between the control and hypoxia groups. Enrichment analysis indicated the enrichment of GO terms including "oxidation-reduction process", "ECM organization", "chaperone cofactor-dependent protein refolding", and "ECM-receptor interaction" KEGG pathway by the DEGs. In addition, between the two groups, a total of 28 SDMs were identified, which were implicated in 13 metabolic pathways, such as "phenylalanine metabolism", "D-amino acid metabolism", and "aminoacyl-tRNA biosynthesis". Results suggest that pearl oysters are exposed to oxidative stress and apoptosis under short-term hypoxia. Also, pearl oysters might adapt to short-term hypoxic treatment by increasing antioxidant activity, modulating immune and biomineralization activities, maintaining protein homeostasis, and reorganizing the cytoskeleton. The results of our study help unveil the mechanisms by which pearl oysters respond adaptively to short-term hypoxia.

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.

Responses of digestive metabolism to marine heatwaves in pearl oysters

Marine heatwaves (MHWs) have increased in intensity and frequency in global oceans, causing deleterious effects on many marine organisms and ecosystems they support. Bivalves are among the most vulnerable taxonomic groups to intensifying MHWs, yet little is known about the underlying mechanisms. Here, we investigated the impact of MHWs on the digestive metabolism of pearl oysters (Pinctada maxima). Two moderate and severe scenarios of MHWs were performed by increasing seawater temperature respectively from 24 °C to 28 °C and 32 °C for 3 days. When subjected to MHWs and with increasing intensity, pearl oysters significantly enhanced their digestive enzymatic activities, such as lipase and amylase. LC-MS-based metabolomics revealed negative responses in the lipid metabolism (e.g., steroid biosynthesis, glycerophospholipid metabolism, and sphingolipid metabolism), the amino acid metabolism (e.g., glutamate, histidine, arginine, and proline), and the B-vitamins metabolism. These findings indicate that the digestive metabolism of marine bivalves can likely succumb to intensifying MHWs events.

The Anti-Photoaging Activity of Peptides from Pinctada martensii Meat

Long-term exposure to ultraviolet-B (UVB) can cause photoaging. Peptides from Pinctada martensii meat have been shown to have anti-photoaging activities, but their mechanism of action is rarely studied. In this study, Pinctada martensii meat hydrolysates (PME) were prepared by digestive enzymes and then separated by ultrafiltration and Sephadex G-25 gel filtration chromatography to obtain a purified fraction (G2). The fraction G2 was identified by ultra-performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS), and peptide sequences were synthesized by solid-phase synthesis. The mechanism of anti-photoaging activities was investigated using a human immortalised epidermal (HaCaT) cell model. Results showed that peptides from Pinctada martensii meat increased UVB-induced cell viability and reduced the contents of interstitial collagenase (MMP-1) and matrix lysing enzyme (MMP-3) in HaCaT cells. Furthermore, the fraction of G2 significantly downregulated the expression of p38, EKR, JNK, MMP-1, and MMP-3 in HaCaT cells. The peptide sequences Phe-His (FH), Ala-Leu (AL), Met-Tyr (MY), Ala-Gly-Phe (AGF), and Ile-Tyr-Pro (IYP) were identified and synthesized. Besides, FH reduced the contents of MMP-1 and MMP-3 in HaCaT cells, combining them effectively in molecular docking analysis. Thus, peptides from Pinctada martensii meat showed anti-photoaging activities and might have the potential to be used as an anti-photoaging agent in functional foods.

A new perspective on endocrine disrupting effects of triphenyltin on marine medaka: From brain transcriptome, gut content metabolome and behavior

Triphenyltin (TPT) is an endocrine contaminant that is often detected in the environment. However, the mechanism of the effects of TPT on biological systems is not fully understood. Here we exposed marine medaka (Oryzias melastigma) to TPT for 21 days. Brain transcriptome, intestinal content metabolism group, and behavior analysis were carried out. Through the comprehensive analysis of multiomics for the in-depth understanding of the ways related to health improvement, we determined that the glycine-serine-threonine metabolic axis was most perturbed by TPT. Through behavioral analysis, it was found that there was behavioral hyperactivity in the exposed group; behavioral hyperactivity may be caused by the interference of TPT with the neuroendocrine system. In order to gain a full understanding of the impacts of TPT on human health, transcriptomic screening of differential genes and an impartial attitude based on bioinformatics were used. Gene-disease interaction analysis using the Comparative Toxicogenomics Database (CTD) revealed the possible effects of TPT on human health. Finally, based on these findings, the relevant adverse outcome pathway (AOP), which is the "epigenetic modification of PPARG leading to adipogenesis," was identified from AOP Wiki. Further research is required to validate the potential AOP of TPT.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

Fluoride exposure induces lysosomal dysfunction unveiled by an integrated transcriptomic and metabolomic study in bone marrow mesenchymal stem cells

Fluoride has received much attention for its predominant bone toxicity in the human body. However, the toxic mechanism of bone injury caused by fluoride exposure remains largely unclear. Bone marrow mesenchymal stem cells (BMSCs) are widely used as model cells for evaluating bone toxicity after environmental toxicant exposure. In this study, BMSCs were exposed to fluoride at 1, 2, and 4 mM for 24 h, and fluoride significantly inhibited cell viability at 2 and 4 mM. A multiomics analysis combining transcriptomics with metabolomics was employed to detect alterations in genes and metabolites in BMSCs treated with 2 mM fluoride. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of transcriptomics profiles identified "lysosomes" as the top enriched pathway, which was severely damaged by fluoride exposure. Lysosomal damage was indicated by decreases in the expression of lysosomal associated membrane protein 2 (LAMP 2) and cathepsin B (CTSB) as well as an increase in pH. Upregulation of the lysosome-related genes Atp6v0b and Gla was observed, which may be attributed to a compensatory lysosomal biogenesis transcriptional response. Interestingly, inhibition of glutathione metabolism was observed in fluoride-treated BMSCs at the metabolomic level. Moreover, an integrative analysis between altered genes, metabolites and lysosome signaling pathways was conducted. Palmitic acid, prostaglandin C2, and prostaglandin B2 metabolites were positively associated with Atp6v0b, a lysosome-related gene. Overall, our results provide novel insights into the mechanism responsible for fluoride-induced bone toxicity.

Integrated application of transcriptomics and metabolomics provides insights into condition index difference mechanisms in the Pacific oyster (Crassostrea gigas)

The condition index (CI) is an economically important tool for assessing the quality of oysters, such as the Pacific oyster Crassostrea gigas. However, little is known about the mechanisms that underlie differences in CI between different C. gigas populations. In this study, we integrated transcriptomic and metabolomic profiling to investigate the mechanisms that underlie the differences between high- and low-CI groups in one- and two-year-old populations of C. gigas. The results indicate that differences in CI were associated with the regulation of growth-related genes, the FoxO signaling pathway, and the complex regulation of carbohydrate, lipid, amino acid, and energy metabolism. Moreover, the mechanisms underlying these differences differed between the populations. This study is the first to elucidate the molecular and chemical mechanisms associated with CI, and the results will be helpful for breeding higher quality oysters.

Integrated application of multi-omics approach and biochemical assays provides insights into physiological responses to saline-alkaline stress in the gills of crucian carp (Carassius auratus)

Given the decline of freshwater resources in recent years, the accessible space for freshwater aquaculture is rapidly shrinking, and aquaculture in saline-alkaline water has become a critical approach to meet the rising demand. However, the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance remains unclear. Here, adult crucian carp (Carassius auratus) were exposed to 60 mmol/L NaHCO₃ for 30 days. It was observed that long-term carbonate alkalinity (CA) exposure not only caused gill oxidative stress but also changed the levels of several physiological parameters associated with ammonia transport, including blood ammonia, urea nitrogen (BUN), glutamine (Gln), and glutamine synthetase (GS). According to the metabolomics study, differential metabolites (DMs) engaged in various metabolic pathways, such as glycerophospholipid metabolism, sphingolipid metabolism, and arachidonic acid metabolism. In addition, transcriptomics data showed that differentially expressed genes (DEGs) were closely related to ammonia transport, apoptosis, and immunological response. In general, comprehensive multi-omics and biochemical analysis revealed that crucian carp might adopt Rh glycoprotein as a carrier to mediate ammonia transport and increase glutamine and urea synthesis under long-term high saline-alkaline stress to mitigate the adverse effects of blocked ammonia excretion. Simultaneously, saline-alkaline stress caused the destruction of the antioxidant system and the disorder of lipid metabolism in the crucian carp gills, which induced apoptosis and immunological response. To our knowledge, this is the first study to investigate fish's molecular and metabolic mechanisms under saline-alkaline stress using integrated metabolomics, transcriptomics, and biochemical assays. Overall, the results of this study provided new insights into the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance.

Mantle Transcriptome Provides Insights into Biomineralization and Growth Regulation in the Eastern Oyster (Crassostrea virginica)

Growth of the eastern oyster Crassostrea virginica, a major aquaculture species in the USA, is highly variable and not well understood at molecular levels. As growth of mollusks is confined in shells constructed by the mantle, mantle transcriptomes of large (fast-growing) and small (slow-growing) eastern oysters were sequenced and compared in this study. Transcription was observed for 31,186 genes, among which 104 genes were differentially expressed between the large and small oysters, including 48 upregulated and 56 downregulated in large oysters. Differentially expressed genes (DEGs) included genes from diverse pathways highlighting the complexity of shell formation and growth regulations. Seventeen of the 48 upregulated DEGs were related to shell matrix formation, most of which were upregulated in large oysters, indicating that large oysters are more active in biomineralization and shell formation. Genomic and transcriptomic analyses identified 22 genes encoding novel polyalanine containing proteins (Pacps) with characteristic motifs for matrix function that are tandemly duplicated on one chromosome, all specifically expressed in mantle and at higher levels in large oysters, suggesting that these expanded Pacps play important roles in shell formation and growth. Analysis of sequence variation identified 244,964 SNPs with 328 associated with growth. This study provides novel candidate genes and markers for shell formation and growth, and suggests that genes related to shell formation are important for the complex regulation of growth in the eastern oyster and possibly other bivalve mollusks. Results of this study show that both transcriptional modulation and functional polymorphism are important in determining growth.

Multiomic Analysis Reveals Core Regulatory Mechanisms underlying Steroidal Glycoalkaloid Metabolism in Potato Tubers

Steroidal glycoalkaloids (SGAs) present in germinated potato tubers are toxic; however, the mechanisms underlying SGA metabolism are poorly understood. Therefore, integrated transcriptome, metabolome, and hormone analyses were performed in this study to identify and characterize the key regulatory genes, metabolites, and phytohormones related to glycoalkaloid regulation. Based on transcriptome sequencing of bud eyes of germinated and dormant potato tubers, a total of 6260 differentially expressed genes were identified, which were mainly responsible for phytohormone signal transduction, carbohydrate metabolism, and secondary metabolite biosynthesis. Two TCP14 genes were identified as the core transcription factors that potentially regulate SGA synthesis. Metabolite analysis indicated that 149 significantly different metabolites were detected, and they were enriched in metabolic and biosynthetic pathways of secondary metabolites. In these pathways, the α-solanine content was increased and the expression of genes related to glycoalkaloid biosynthesis was upregulated. Levels of gibberellin and jasmonic acid were increased, whereas that of abscisic acid was decreased. This study lays a foundation for investigating the biosynthesis and regulation of SGAs and provides the reference for the production and consumption of potato tubers.

Pyruvate kinase isoform M2 impairs cognition in systemic lupus erythematosus by promoting microglial synaptic pruning via the β-catenin signaling pathway

BACKGROUND

Neuropsychiatric systemic lupus erythematosus (NPSLE) is a severe complication, which involves pathological damage to the brain and cognitive function. However, its exact mechanism of action still remains unclear. In this study, we explored the role of microglia in the cognitive dysfunction of NPSLE mice. We also analyzed and compared the metabolites in the hippocampal tissues of the lupus model and control mice.

METHODS

MRL/MpJ-Fas^lpr (MRL/lpr) female mice were used as the NPSLE mouse model. Metabolomics was used to assess hippocampal glycolysis levels. Glucose, lactic acid, IL-6, and IL-1β of the hippocampus were detected by ELISA. Based on the glycolysis pathway, we found that pyruvate kinase isoform M2 (PKM2) in the hippocampus was significantly increased. Thus, the expression of PKM2 was detected by qRT-PCR and Western blotting, and the localization of PKM2 in microglia (IBA-1⁺) or neurons (NeuN⁺) was assessed by immunofluorescence staining. Flow cytometry was used to detect the number and phenotype of microglia; the changes in microglial phagocytosis and the β-catenin signaling pathway were detected in BV2 cells overexpressing PKM2. For in vivo experiments, MRL/lpr mice were treated with AAV9-shPKM2. After 2 months, Morris water maze and conditional fear tests were applied to investigate the cognitive ability of mice; H&E and immunofluorescence staining were used to evaluate brain damage; flow cytometry was used to detect the phenotype and function of microglia; neuronal synapse damage was monitored by qRT-PCR, Western blotting, and immunofluorescence staining.

RESULTS

Glycolysis was elevated in the hippocampus of MRL/lpr lupus mice, accompanied by increased glucose consumption and lactate production. Furthermore, the activation of PKM2 in hippocampal microglia was observed in lupus mice. Cell experiments showed that PKM2 facilitated microglial activation and over-activated microglial phagocytosis via the β-catenin signaling pathway. In vivo, AAV9-shPKM2-treated mice showed decreased microglial activation and reduced neuronal synapses loss by blocking the β-catenin signaling pathway. Furthermore, the cognitive impairment and brain damage of MRL/lpr mice were significantly relieved after microglial PKM2 inhibition.

CONCLUSION

These data indicate that microglial PKM2 have potential to become a novel therapeutic target for treating lupus encephalopathy.

Integrated analysis of physiological, transcriptomics and metabolomics provides insights into detoxication disruption of PFOA exposure in Mytilus edulis

Perfluorooctanoic acid (PFOA), a persistent environmental contaminant, resists environmental degradation and bioaccumulates in food chains. Lots of literatures have proved that PFOA exposure could disrupt detoxifying function in a variety of organisms, however, it still remained poorly known about this in mollusk. Here, we examined physiological, transcriptomic, and metabolomic responses to PFOA in Mytilus edulis, a model organism frequently used in studies of aquatic pollution. We aimed to characterize PFOA-induced stress responses and detoxification mechanisms. PFOA exposure significantly altered antioxidant enzyme activity levels and the abundances of lipid peroxidation products. In addition, transcriptomic analysis indicated that several genes associated with oxidative stress and detoxication were differentially expressed after PFOA exposure. In combination, transcriptomic and metabolomic analyses showed that PFOA exposure disturbed several metabolic processes in M. edulis, including the lipid metabolism, amino acid metabolism, and carbohydrate metabolism etc. Molecular examination and enzymes assay of PFOA-exposed M. edulis after a 7-day depuration period still did not recover to control levels. The Pathway enrichment analysis proved that several pathways related to detoxification, such as c-Jun N-terminal kinase (JNK) and p38-dependent mitogen-activated protein kinase (MAPK) pathway, Peroxisome proliferator-activated receptor γ (PPARγ) pathway etc, were obviously affected. The present work verifies firstly PFOA disruption to molluscan detoxification and identifies the key pathways to understand the molecular mechanisms thereof. This study provides new insights into the detoxication disruption invoked in response to PFOA exposure in M. edulis.

Integrated transcriptomics and metabolomics analyses reveal benzo[a]pyrene enhances the toxicity of mercury to the Manila clam, Ruditapes philippinarum

Mercury (Hg²⁺) and benzo[a]pyrene (BaP) are ubiquitous and persistent pollutants with multiple toxicities in bivalve molluscs. Here, the toxicological responses in the gills of Manila clams, Ruditapes philippinarum, to Hg²⁺ (10 μg L^-1), BaP (3 μg L^-1), and their mixture were analysed using transcriptomics and metabolomics approaches. Comparisons of the transcriptomes and metabolomes of Hg²⁺-and/or BaP-treated clams with control animals revealed the involvement of the detoxification metabolism, immune defence, energy-related pathways, and osmotic regulation in the stress response of R. philippinarum. Exposure to Hg²⁺ alone primarily enhanced the detoxification and energy metabolic pathways by significantly increasing the expression of genes associated with heat-shock proteins and oxidative phosphorylation. However, co-exposure to Hg²⁺ and BaP caused greater immunotoxicity and disrupted detoxification metabolism, the TCA cycle, glycolysis, and ATP generation. The expression levels of cytochrome P450 1A1 (CYP1A1), multidrug resistance-associated protein 1 (MRP1), and myosin (MYO), and the activity of electron transport system (ETS) in gills were detected, supporting the underlying toxic mechanisms of Hg²⁺ and BaP. We suggest that the presence of BaP enhances the toxicity of Hg²⁺ by 1) hampering the detoxification of Hg²⁺, 2) increasing the immunotoxicity of Hg²⁺, and 3) constraining energy availability for clams.

Identification and Allelic Variants Associated With Cold Tolerance of PmPIAS in Pinctada fucata martensii

The protein inhibitor of activated STAT (PIAS) functions in diverse aspects, including immune response, cell apoptosis, cell differentiation, and proliferation. In the present study, the PIAS in the pearl oyster Pinctada fucata martensii was characterized. The sequence features of PmPIAS were similar to that of other PIAS sequences with PIAS typical domains, including SAP, Pro-Ile-Asn-Ile-Thr (PINIT), RLD domain, AD, and S/T-rich region. Homologous analysis showed that PmPIAS protein sequence showed the conserved primary structure compared with other species. Ribbon representation of PIAS protein sequences also showed a conserved structure among species, and the PINIT domain and RLD domain showed the conserved structure compared with the sequence of Homo sapiens. The expression pattern of PmPIAS in different tissues showed significant high expression in the gonad. PmPIAS also exhibited a significantly higher expression in the 1 and 2 days after cold tolerance stress (17°C) and showed its potential in the cold tolerance. The SNP analysis of the exon region of PmPIAS obtained 18 SNPs, and among them, 11 SNPs showed significance among different genotypes and alleles between cold tolerance selection line and base stock, which showed their potential in the breeding for cold tolerance traits.

Integrated GC-MS- and LC-MS-Based Untargeted Metabolomics Studies of the Effect of Vitamin D3 on Pearl Production Traits in Pearl Oyster Pinctada fucata martensii

Pearl oyster Pinctada fucata martensii is widely recognized for biomineralization and has been cultured for high-quality marine pearl production. To ascertain how dietary vitamin D3 (VD3) levels affect the features of pearl production by P. f. martensii and discover the mechanisms regulating this occurrence, five experimental diets with variable levels of VD3 were used with inclusion levels of 0, 500, 1,000, 3,000, and 10,000 IU/kg. The distinct inclusion levels were distributed into five experimental groups (EG1, EG2, EG3, EG4, and EG5). All the experimental groups were reared indoors except the control group (CG) reared at the sea. Pearl oysters, one year and a half old, were used in the grafting operation to culture pearls. During the growing period that lasted 137 days, EG3 had the highest survival rate, retention rate, and high-quality pearl rate. A similar trend was found for EG3 and CG with significantly higher pearl thickness and nacre deposition rates than other groups, but no significant differences were observed between them. A metabolomics profiling using GC–MS and LC–MS of pearl oysters fed with low quantities of dietary VD3 and optimal levels of dietary VD3 revealed 135 statistically differential metabolites (SDMs) (VIP > 1 and p < 0.05). Pathway analysis indicated that SDMs were involved in 32 pathways, such as phenylalanine metabolism, histidine metabolism, glycerophospholipid metabolism, alanine aspartate and glutamate metabolism, arginine and proline metabolism, glycerolipid metabolism, amino sugar and nucleotide sugar metabolism, and tyrosine metabolism. These results provide a theoretical foundation for understanding the impacts of VD3 on pearl production traits in pearl oyster and reinforce forthcoming prospects and application of VD3 in pearl oyster in aquaculture rearing conditions.

A comprehensive integrated transcriptome and metabolome analyses to reveal key genes and essential metabolic pathways involved in CMS in kenaf

Numbers of critical genes and pathways were found from the levels of transcriptome and metabolome, which were useful information for understanding of kenaf CMS mechanism. Cytoplasmic male sterility (CMS) is a maternally inherited trait in higher plants that leads to the inability to produce or release functional pollen. However, there is lack of comprehensive studies to reveal the molecular basis of CMS occurrence in kenaf. Herein, we performed transcriptome and UPLC-MS-based metabolome analyses in the anthers of a CMS (UG93A) and its maintainer (UG93B) to sort out essential genes and metabolites responding to CMS in kenaf. Transcriptome characterized 7769 differentially expressed genes (DEGs) between these two materials, and pathway enrichment analysis indicated that these DEGs were involved mainly in pentose and glucuronate interconversions, starch and sucrose metabolism, taurine and hypotaurine metabolism. In the metabolome assay, a total of 116 significantly different metabolites (SDMs) were identified between the CMS and its maintainer line, and these SDMs were involved in eight KEGG pathways, including flavone and flavonol biosynthesis, glycerophospholipid metabolism, flavonoid biosynthesis, glycosylphosphatidylinositol-anchor biosynthesi. Integrated analyses of transcriptome and metabolome showed that 50 genes had strong correlation coefficient values (R² > 0.9) with ten metabolites enriched in six pathways; notably, most genes and metabolites of flavonoid biosynthesis pathways and flavone and flavonol biosynthesis pathways involved in flavonoids biosynthetic pathways were downregulated in CMS compared to those in maintainer. Taken together, the decreased accumulation of flavonoids resulted from the compromised biosynthesis pathways coupled with energy deficiency in the anthers may contribute largely to CMS in UG93A of kenaf.

Hepatopancreatic metabolomics shedding light on the mechanism underlying unsynchronized growth in giant freshwater prawn, Macrobrachium rosenbergii

The giant freshwater prawn, Macrobrachium rosenbergii (M. rosenbergii) as an important freshwater aquaculture species with high commercial value, exhibited unsynchronized growth. However, the potentially metabolic mechanism remains unclear. In this study, we used liquid chromatography tandem mass spectrometry (LC-MS/MS) to investigate the hepatopancreatic metabolic profiles of twenty giant freshwater prawns between the fast-growing group and slow-growing group. In the metabolomics assay, we isolated 8,293 peaks in positive and negative iron mode. Subsequently, 44 significantly differential metabolites were identified. Functional pathway analysis revealed that these metabolites were significantly enriched in three key metabolic pathways. Further integrated analysis indicated that glycerophospholipid metabolism and aminoacyl-tRNA biosynthesis have significant impact on growth performance in M.rosenbergii. Our findings presented here demonstrated the critical metabolites and metabolic pathways involved in growth performance, moreover provided strong evidence for elucidating the potentially metabolic mechanism of the unsynchronized growth in M. rosenbergii.

Regulatory Role of Retinoic Acid in Male Pregnancy of the Seahorse

Seahorses epitomize the exuberance of evolution. They have the unique characteristic of male pregnancy, which includes the carrying of many embryos in a brood pouch that incubates and nourishes the embryos, similar to the mammalian placenta. However, the regulatory networks underlying brood pouch formation and pregnancy remain largely unknown. In this study, comparative transcriptomic and metabolomic profiling on the lined seahorse Hippocampus erectus, with unformed, newly formed, and pregnant brood pouches identified a total of 141 and 2,533 differentially expressed genes together with 73 and 121 significantly differential metabolites related to brood pouch formation and pregnancy, respectively. Specifically, integrative omics analysis revealed that retinoic acid (RA) synthesis and signaling pathway played essential roles in the formation of the brood pouch and pregnancy. RA might function upstream of testosterone and progesterone, thereby directly influencing brood pouch formation by regulating the expression of fshr and cyp7a1. Our results also revealed that RA regulates antioxidant defenses, particularly during male pregnancy. Alternatively, pregnancy caused a consistent decrease in RA, canthaxanthin, astaxanthin, and glutathione synthetase, and an increase in susceptibility to oxidative stress, which may balance brood pouch development and reproduction in seahorses and pave the way to successful gestation.

Molecular Cloning and Polymorphism Analysis of PmFGF18 from Pinctada fucata martensii

Fibroblast growth factor 18 (FGF18) plays an important functional role in skeletal growth and development. The FGF18 gene was characterized in pearl oyster Pinctada fucata martensii (PmFGF18) with the full-length sequence containing an open reading frame of 714 bp encoding 237 amino acids. The domain analysis of PmFGF18 showed a distinctive FGF domain, with a high similarity to FGF18 protein sequences from Crassostrea gigas (43.35%) and C. virginica (37.43%). PmFGF18 expression was revealed in all analyzed tissues with a significantly higher expression level in the fast-growing group than the slow-growing group. The analysis of PmFGF18 polymorphism demonstrated 33 SNPs (single nucleotide polymorphisms) in the CDS and promoter region of PmFGF18 sequence. Association analysis revealed 19 SNPs (2 SNPs from CDS and 17 SNPs from the promoter region) associating significantly with growth traits. Among the associated SNPs, one SNP g.50918198 A > C was verified in the other breeding line. Therefore, PmFGF18 can be utilized as a candidate gene for growth, and its related SNPs could be used in selective breeding of P. f. martensii for the improvement of growth traits.

Variation of Microbiological and Biochemical Profiles of Laowo Dry-Cured Ham, an Indigenous Fermented Food, during Ripening by GC-TOF-MS and UPLC-QTOF-MS

Fermented foods have unique microbiota and metabolomic profiles that can support dietary diversity, digestion, and gut health of consumers. Laowo ham (LWH) is an example of an indigenous fermented food from Southwestern China that has cultural, ecological, economic, and health significance to local communities. We carried out ethnobiological surveys coupled with metagenomic and metabolomic analyses using GC-TOF-MS and UPLC-QTOF-MS to elucidate the microbiota and metabolic profiles of LWH samples at different ripening stages. The results from high-throughput sequencing showed a total of 502 bacterial genera in LWH samples with 12 genera of bacteria and 6 genera of fungi identified as dominant groups. This is the first study to our knowledge to report the bacteria of Lentibacillus and Mesorhizobium along with fungi Eremascus and Xerochrysium on a fermented meat product. Findings further revealed that the metabolite profiles among LWH samples were significantly different. In total, 27 and 30 metabolites from GC-TOF-MS and UPLC-QTOF-MS analysis, respectively, were annotated as highly discriminative metabolites. Among the differential compounds, the relative contents of most amino acids showed the highest in the LWH sample ripened for two years, while some metabolites with potential therapeutic effects such as levetiracetam were the most abundant in the LWH sample ripened for three years. The correlation analysis indicated that the dominant microbes were closely related to differential metabolites, highlighting the importance of their functional characterization. Findings indicate that the consumption of LWH contributes to microbiological and chemical diversity of human diets as well as suggests efficacy of combining GC-MS and LC-MS to study the metabolites in dry-cured meat products.

PacBio single molecule long-read sequencing provides insight into the complexity and diversity of the Pinctada fucata martensii transcriptome

BACKGROUND

The pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms.

RESULTS

Herein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment.

CONCLUSIONS

Our results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.

Evaluation of the Metabolic Effects of Hydrogen Sulfide on the Development of Bombyx mori (Lepidoptera: Bombycidae), Using Liquid Chromatography-Mass Spectrometry-Based Metabolomics

Hydrogen sulfide (H2S) is a highly poisonous gas with an unpleasant smell of rotten eggs. Previous studies of H2S have primarily focused on its effects on mammalian nervous and respiratory systems. In this study, silkworm developmental parameters and changes in metabolites in response to H2S exposure were investigated using a hemolymph metabolomic approach, based on liquid chromatography-mass spectrometry (LC-MS). The developmental parameters, body weight, cocoon weight, cocoon shell weight, and cocoon shell ratio, were noticeably increased following H2S exposure, with the greatest effects observed at 7.5-μM H2S. Metabolites upregulated under H2S exposure (7.5 μM) were related to inflammation, and included (6Z, 9Z, 12Z)-octadecatrienoic acid, choline phosphate, and malic acid, while hexadecanoic acid was downregulated. Identified metabolites were involved in biological processes, including pyrimidine, purine, and fatty acid metabolism, which are likely to affect silk gland function. These results demonstrate that H2S is beneficial to silkworm development and alters metabolic pathways related to spinning function and inflammation. The present study provides new information regarding the potential functions of H2S in insects and metabolic pathways related to this phenomenon.

Metabolome and transcriptome profiling reveals quality variation and underlying regulation of three ecotypes for Cistanche deserticola

Cistanche deserticola is a plant used both as food and medicine. We are interested in understanding how C. deserticola responds to environmental conditions. Samples were collected from three ecotypes grown in saline-alkali land, grassland and sandy land. Transcriptome and metabolome analysis were performed by using RNA-seq and LC-ESI-MS/MS. Among 578 metabolites identified, 218, 209 and 215 compounds were found differentially produced among the three ecotypes. Particularly, 2'-acetylacteoside, belonging to phenylethanoid glycosides (PhGs) is the most significantly differentially produced with a VIP > 0.5 and fold change > 2, representing a potential chemical marker to distinguish the three ecotypes. RNA-Seq analysis revealed 52,043 unigenes, and 947, 632 and 97 of them were found differentially expressed among the three ecotypes. Analysis of the correlation between the metabolome profiles and transcriptome profiles among three ecotypes identified that the 12 key genes related to PhGs biosynthesis were differentially expressed. Particularly, the expression of PAL, ALDH and GOT genes were significantly up-regulated in saline-alkali land compared to the other two. In summary, we found PhGs content was higher in saline-alkali land compared with other ecotypes. This is likely due to the up-regulation of the PhGs biosynthetic genes in response to the saline-alkali conditions.

Combined Application of High-Throughput Sequencing and Metabolomics Reveals Metabolically Active Microorganisms During Panxian Ham Processing

Panxian ham, a traditional Chinese dry-cured ham, is protected by national geographical indication. Similar to other fermented foods, the microbial population of dry-cured ham is pivotal to taste and flavor formation. This study aimed to establish the relationship between microorganisms and metabolites during the spontaneous fermentation of Panxian ham. Multivariate analysis based on metabolomics data revealed that continuous metabolic changes occurred during the entire fermentation process, with the most significant changes occurring in the initial stage of ripening. Thirty-one significantly different metabolites (SDMs) were identified as discriminant factor, and pathway analysis suggested that these metabolites were involved in 30 pathways, including alanine, aspartate, and glutamate metabolism; glycine, serine, and threonine metabolism; and arginine and proline metabolism. Microbial community analysis using the Illumina MiSeq platform indicated that the bacterial community was more complex than the fungal community, and their succession regulation differed during processing. At the genus level, 11 bacteria and five fungi were identified as core microbes, of which Staphylococcus was the dominant bacteria and Debaryomyces and Aspergillus were the dominant fungi. Further, statistical redundancy analysis (RDA) indicated that Staphylococcus, Debaryomyces, and Chromohalobacter promoted the production of amino and fatty acids; Cobetia and Aspergillus were associated with sugar metabolism, and Kushneria, Penicillium, and Yamadazyma were closely related with organic acids. These findings provide fundamental knowledge regarding the metabolically active microorganisms in Panxian ham, helping industrial processors to develop effective strategies for standardizing quality parameters.

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

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