Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 14 August 2025

Sec. Parasite and Host

Volume 15 - 2025 | https://doi.org/10.3389/fcimb.2025.1568349

This article is part of the Research TopicInteractions at the Viral-Host Nexus in animals: From Omics Insights to Immune ModulationView all 6 articles

Proteomics and phosphoproteomics of freshwater mollusk carcasses reveal novel insights as potential food source

Gongzhen Liu*Gongzhen Liu*Shengkang WangShengkang WangTielin WangTielin WangChangan ZhuChangan ZhuChanglong LiChanglong LiXinyi ZhangXinyi Zhang
  • College of Agriculture and Forestry, Linyi University, Linyi, Shandong, China

Cipangopaludina chinensis, Radix auricularia, and Nodularia douglasiae are common freshwater mollusks widely distributed in lakes and rivers across China. In this study, (DIA) Data - independent acquisition method was used to analyze the proteomics and phosphoproteomics among three species. The results showed that a total of 1,382 proteins were identified from mollusk carcasses, with 492 proteins were quantitively analyzed. These proteins are mainly involved in amino acid nutrition and metabolism, including glutamate dehydrogenase, glyceraldehyde-3-phosphate dehydrogenase, and tyrosinase copper-binding domain-containing protein. Bioinformatics analysis revealed that the identified proteins are primarily involved in carbon metabolism, microbial metabolism, ubiquitin-mediated proteolysis, and RAS signaling pathways. Notably, this study offers valuable insights into the differential protein expression of proteins across species. Furthermore, proteomics analysis revealed several novel proteins, which helps to improve our understanding of amino acid metabolism, protein diversity, and signaling pathways in these freshwater mollusks.

1 Introduction

Cipangopaludina chinensis (C. chinensis), Radix auricularia (R.auricularia), and Nodularia douglasiae (N.douglasiae) are three freshwater mollusk species widely distributed across Asian countries, including China, Japan, and Korea, as well as other parts of the world (Zhou et al., 2022a, Zhou et al., 2022a; Zhao et al, 2022; Choi et al., 2020; Schell et al., 2017). Mollusks inhabit a variety of freshwater environments such as rivers, lakes, rice fields, swamps, and ditches (Marrone et al., 2019). Taxonomically, both C. chinensis and R.auricularia belong to the phylum Mollusca and class Gastropoda; however, C. chinensis is classified under the family Viviparidae, while R.auricularia belongs to the family Lymnaeidae (Juhász and Lawton, 2022). In contrast, N.douglasiae belongs to the phylum Mollusca, class Bivalvia and family Unionidae, which is very similar to another relative species, Nodularia breviconcha (N.breviconcha) (Kim et al., 2021; Choi et al., 2020).

Aquatic mollusks play a vital role in ecosystems, significantly influencing food webs and nutrient cycling (Darrigran et al., 2025; Vasta et al., 2015). Economically, they contribute positively to aquaculture by serving as a food source and providing protein-rich nutrition (such as Viviparidae) (Liu et al., 2024), but they can also cause substantial economic losses through macrofouling and pose health risks, particularly in cases involving non-native species (Li et al., 2013; Darrigran et al., 2025). Freshwater mollusks are further recognized as key water quality indicators, natural biofilters, and agents of reproduction and dispersal. First, mollusks are an important part of the food chain in aquatic ecosystems by feeding on aquatic plants, plankton, organic detritus, and dead organisms (Deidda et al., 2021). Second, their sensitivity to environmental changes makes them effective bioindicators for monitoring water quality, including nutrient levels, pollution, and oxygen content (Elder and Collins, 1991; Mienis and Ashkenazi, 2011; Fernandes et al., 2011). Third, mollusks have a strong filtration capacity, enabling them to remove plankton and organic debris from water while metabolizing harmful substances and improving water quality (Giansante and Pelini, 2007). Additionally, they can form intricate symbiotic relationships with other aquatic organisms, enhancing the ecological stability of aquatic ecosystems (Chakraborty and Joy, 2020; de la Ballina et al., 2022). Finally, mollusks contribute to biodiversity by acting as dispersers and reproductive agents, spreading genetic material across diverse habitats and promoting species diversity (Wayne, 2001; Beirão et al., 2019; Gaudin-Zatylny et al., 2022).

C. Chinensis carcasses are an important food source for both humans and animals, and several recent studies have explored this species in greater depth. Single-molecule real-time (SMRT) sequencing technology has been used to analyze its full-length transcriptome, revealing that metabolic molecules, signal transduction proteins, and immune-related proteins are involved in drug metabolism and immune response (Zhou et al., 2022a). Meanwhile, high levels of aryl sulfatase and β-glucuronidase were identified in C. chinensis, both of which play a key role in the cleavage process of bound natural estrogens (C-NEs) (Zhao et al., 2022). Notably, a comparative analysis of shell proteins between Pomacea canaliculata (P. canaliculata) and C. chinensis showed that carbonic anhydrase was absent in both species, suggesting that freshwater snails may have a unique mechanism for regulating shell calcification and mineralization (Liu et al., 2023).

R. auricularia poses potential health risks, as it serves as an intermediate host for zoonotic trematodes. However, Radix auricularia also plays an important ecological role in aquatic sphere. Radix auricularia is an integral part of freshwater ecosystems, contributing to the balance of the food web, water quality maintenance, and the creation and modification of habitats. Its ecological role highlights the importance of conserving these organisms and their habitats for the overall health and functioning of freshwater environments. These include Clinostomum complanatum (Chung et al., 1998), Trichobilharzia franki (Juhász et al., 2022), and Petasiger exaeretus (Bashe and Ali, 2019). The draft genome of the pulmonate freshwater snail R.auricularia has been sequenced, providing a valuable foundation for genomic and population genetics research (Schell et al., 2017). Furthermore, phylogenetic analyses exploring the evolution of Toll-like receptors (TLRs) in this species revealed that, unlike many other mollusks, R. auricularia also possesses class I TLRs (Juhász et al., 2022).

N.douglasiae is a freshwater bivalve mollusk with distinct morphological characteristics that differs significantly from C. chinensis and R.auricularia. Investigating the genetic structure and diversity of N. douglasiae from the Yangtze River basin in China has provided important insights for the conservation of genetic diversity and effective management of riverine ecosystems (Liu et al., 2017). Additionally, the complete mitochondrial genomes of N.douglasiae and N.nipponensis have been sequenced from the Lake Biwa system in Japan (Mabuchi et al., 2021). Mitochondrial CO1 and 16S rRNA gene sequences from N. douglasiae have also been used to analyze phylogenetic relationships and population genetic structure, revealing close genetic similarity to N.breviconcha from the Korean Peninsula (Choi et al., 2020; Kim et al., 2021). Recent studies have further demonstrated that the kidneys of N.douglasiae improve water excretion efficiency in freshwater mussels through counter-current systems (Nikishchenko et al., 2022).

Although several studies have examined these three freshwater species, little research has been conducted at the protein level. The carcasses of all three species contain abundant proteins and proteases with potential for future applications. In this study, we aimed to analyze proteomic profile differences in mollusk carcasses among three species. Our results reveal novel proteins and differentially expressed proteins (DEPs) across the three species, offering valuable insights into protein composition, nutritional potential, and involvement in key signaling pathways.

2 Materials and methods

2.1 Samples

The three species (Cipangopaludina chinensis, Radix auricularia, and Nodularia douglasiae) were initially collected from the Yi River Basin in Shandong Province, China (latitude 34°22′–36°13′ N; longitude 117°24′–119°11′ E), a region characterized by a temperate monsoon climate.). We choose the same season in clean waters areas for three species. We used C.chinensis (2.8cm×1.5cm), R.auricularia (1.5cm×0.5cm) and N.douglasia (6.0cm×1.5cm) carcasses (muscle) respectively as protein samples for proteomic and phosphoproteomics analysis. All protein samples were taken from the same part of the carcass, after postmortem examination immediately transported on dry ice. The experimental groups were designated as Group A (C. chinensis), Group B (R. auricularia), and Group C (N. douglasiae), with each individual treated as a separate biological replicate.

2.2 High -performance liquid chromatography fractionation

Samples were concentrated, digested, and extracted, and the protein concentrations were determined using the Bradford assay. Protein sample preparation followed methods described previously (Liu et al., 2022).

2.3 LC−MS/MS analysis

For construction of the transition library, shotgun proteomics was performed using an Ultimate 3000 Ultra High-Performance Liquid Chromatography (UHPLC) system (Dionex) coupled with a Q Exactive HF mass spectrometer (Thermo Fisher Scientific). The detailed procedures were based on previously established protocols (Liu et al., 2022).

2.4 Phosphorylation analysis

For phosphopeptide enrichment, an appropriate volume of PuriMag Si-TiO2 magnetic beads was added to a 1.5 mL tube and washed three times with sample buffer (1 M glycolic acid in 80% acetonitrile [ACN] and 5% trifluoroacetic acid [TFA]). The sample buffer was then added and thoroughly mixed, followed by incubation at room temperature for 20 minutes. After incubation, the supernatant was transferred to a new tube. The beads were washed once with a sample buffer, twice with washing buffer (80% ACN, 1% TFA), and twice with an additional buffer (10% ACN, 0.2% TFA). Phosphopeptides were eluted with 5% ammonium hydroxide (NH4OH) and acidified with formic acid.

A second enrichment was performed using immobilized metal affinity chromatography (IMAC). IMAC beads (200 uL) were transferred to a 1.5 mL tube, centrifuged at 2,100 ×g for 30 seconds, and washed with 6% acetic acid. The beads were incubated in 0.1 M FeCl3 in 6% acetic acid for 2 hours, then washed three times with 6% acetic acid. A mixed buffer (250 mM acetic acid, 30% ACN) was added to a 50% slurry. The dried peptide sample was resuspended in this slurry at a concentration of 1 μg/μL and incubated at room temperature with vigorous shaking for 1 hour. After incubation, the supernatant was transferred to a new tube. The beads were washed twice with buffer and once with water. Phosphopeptides were then eluted with 5% NH4OH and acidified with formic acid. All enriched samples were desalted and vacuum freeze-dried. Finally, phosphopeptides enriched using PuriMag Si-TiO2 and IMAC were pooled and desalted.

2.5 Bioinformatic analysis

The identified proteins were classified using Gene Ontology (GO) annotation. The GO annotation proteome was derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/). The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate protein pathways. KEGG pathway mapping was carried out using the KEGG Automatic Annotation Server (KAAS).

3 Results

3.1 Identified proteins

Each individual carcass of C. chinensis, R. auricularia, and N. douglasiae was analyzed using proteomic and phosphoproteomic methods. phosphoproteomics. Cellular component (CC), molecular function (MF), and biological process (BP)—as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were used to examine protein expression differences among the three pairwise comparisons: Avs-B, Avs-C and Bvs-C.

A total of 1,382 proteins were identified, and 492 were quantitatively analyzed. Among these, 100 (62), 66 (101), and 38 (98) significantly upregulated (downregulated) differentially expressed proteins (DEPs) were detected in pairwise comparisons (A-vs-B, A-vs-C, and B-vs-C) (Figures 1a–f).

Figure 1
The image consists of six panels. Panel (a) displays a volcano plot for A vs B, showing log2 ratios against negative log10 p-values, with red and green dots indicating upregulated and downregulated values. Panel (b) shows a similar plot for A vs C. Panel (c) presents B vs C. Panels (a), (b), and (c) share a similar structure. Panel (d) is a bar graph depicting identified peptide distribution by length, with height decreasing as length increases. Panel (e) plots mass error distribution, showing peptide score versus peptide mass error in PPM. Panel (f) is a heatmap illustrating data correlations, with a color gradient indicating values.

Figure 1. Total identified proteins. (a-c) Protein expression ratios of A-vs-B, A-vs-C, and B-vs-C; (d) Distribution of identified peptides; (e) Mass error distribution; (f) Repeatability analysis across groups A-vs-B, A-vs-C, and B-vs-C.

3.2 Functional classification of DEPs

DEPs were defined based on a protein quantification threshold of fold change ≥ 1.2 and p-value ≤ 0.05. Functional classification was conducted for all regulated proteins, including assessments of cellular components, molecular functions, biological processes, and subcellular localization. Across all three comparisons (A-vs-B, A-vs-C, and B-vs-C), the DEPs were primarily classified into the following functional categories: cytoskeleton (17.92%–20.98%), ribosomal structure and biogenesis (11.19%–12.57%), signal transduction mechanisms (9.09%–10.93%), energy production and conversion (8.09%–11.19%), and osttranslational modification, protein turnover, and chaperones (6.29%–9.83%) (Figures 2a, d, g). Functional classification revealed clear differences among the three comparison groups. DEPs were predominantly involved in cellular components such as cells, cell parts, organelles, membranes, macromolecular complexes, and cell junctions (shown in green in Figures 2b, e, h). In terms of molecular function, these proteins were associated with binding, catalytic activity, and structural molecular activity (blue columns in Figures 2c, f, i). Biological process analysis showed that DEPs were involved in cellular processes, single-organism processes, metabolic processes, localization, regulation of biological processes (both positive and negative), response to stimuli, and overall biological regulation (orange column in Figure 2). Subcellular localization results indicated that most DEPs were concentrated in the cytosol (38.97%–49.38%), nucleus (17.96%–20.99%), and mitochondria (9.26%–10.29%). Notably, cytosolic proteins exhibited a nearly 10% difference in abundance across the three comparisons.

Figure 2
Grouped bar and doughnut charts show gene ontology classifications. Charts (a), (d), and (g) are bar graphs displaying GO distributions. Charts (b), (e), and (h) are stacked bar graphs with categories: biological process, cellular component, and molecular function. Charts (c), (f), and (i) are doughnut charts representing percentages of each category.

Figure 2. GO classification involved in cellular components, molecular functions, biological processes and subcellular location in three species. (a, d, g) All regulated classification in A -vs- B, A -vs- C, B -vs- C, respectively; (b, e, h) The category DEPs in cellular components, molecular function, biological processes in A -vs- B, A -vs- C, B -vs- C, respectively; (c, f, i) The DEPs subcellular location in A -vs- B, A -vs- C, B -vs- C, respectively. GO Annotation: Gene Ontology (GO) is a major bioinformatics initiative to unify the representation of gene and gene product attributes across all species.

3.3 Functional enrichment of DEPs

In contrast to the functional classification, a functional enrichment analysis of all 1,260 proteins highlighted distinct differences in protein domains, biological processes, and KEGG pathways across the three comparison groups.

Domain enrichment analysis revealed 158 DEPs in the Avs. B group, including 97 upregulated and 61 downregulated proteins. The majority of enriched domains were NAD(P)-binding domains and EF-hand domain EF-hand domains (Figure 3a). In the Avs. C group, 163 DEPs were identified—64 upregulated and 99 downregulated (Figure 3d) Tubulin—with enrichment in domains related to tubulin, including the Tubulin (FtsZ, GTPase) domain Tubulin (FtsZ, GTPase domain) and Tubulin (conserved site)conserved tubulin sites. In the Bvs. C group, 132 DEPs were detected (35 upregulated and 97 downregulated), with actin family domains being the most enriched (Figure 3g). actin family. Notably, the upregulated and downregulated DEPs exhibited distinct domain patterns across the three comparisons.

Figure 3
Nine charts labeled a to i, showing statistics of regulated enrichment using dot plots. Each chart represents different categories: Domain (a, d, g), Biological Process (b, e, h), and KEGG Pathway (c, f, i). The x-axis shows negative logarithm of Fisher's exact test p-value, y-axis lists categories, dot size indicates number of differential proteins, and color represents log2 fold enrichment. Color scale ranges from blue to red.

Figure 3. The functional enrichment of DEPs is involved in domains, biological processes, and KEGG in three species. (a, d, g) Enrichment domains in A -vs- B, A -vs- C, and B -vs- C, respectively; (b, e, h) Enrichment of biological processes in A -vs- B, A -vs- C, and B -vs- C, respectively; (c, f, i) Enrichment of KEGG in A -vs- B, A -vs- C, B -vs- C, respectively.

Biological process enrichment analysis showed that oxidation-reduction processes and energy derivation from the oxidation of organic compounds in the A-vs.-B group (Figure 3b). Compared with the A-vs-B group, the majority of DEPs in the A-vs-C group were involved in the cellular amide metabolic process, translation, amide biosynthetic process, and peptide biosynthetic process (Figure 3e). In the B-vs-C group, DEPs were primarily associated with the generation of precursor metabolites and energy, including the metabolism of glycosyl compounds, nucleosides, purines, ribonucleosides, and purine ribonucleosides, as well as energy derivation through the oxidation of organic compounds and muscle structure development (Figure 3h).

KEGG pathway enrichment analysis identified 414 proteins, with 85 DEPs in the A-vs,-B group—52 upregulated and 33 downregulated (Figure 3c). These DEPs were mainly associated with carbon metabolism, microbial metabolism in diverse environments, and the biosynthesis of antibiotics. In the A-vs-C group, 102 DEPs were found (36 upregulated and 66 downregulated), primarily involved in ribosome function and adrenergic signaling in cardiomyocytes (Figure 3f). For the B-vs. C group, 78 DEPs were enriched (18 upregulated and 60 downregulated) (Figure 3i), with dominant pathways including the citrate cycle (TCA cycle), dopaminergic synapse, insulin signaling pathway, alcoholism, insulin resistance, and plant–pathogen interactions.

Molecular function enrichment analysis of 983 proteins identified 131 DEPs in the A-vs-B group, comprising 78 upregulated and 53 downregulated proteins (Figure 4a). The top two enriched molecular functions were structural molecule activity and oxidoreductase activity.

Figure 4
Bar graphs showing different protein activities and their mappings, differentiated by color codes for Diffmapping, ProteinMapping, DiffNum, and ProteinNum. Graphs labeled a, b, and c display activities such as biotin carboxylase and magnesium ion binding with varied lengths indicating different mapping numbers.

Figure 4. Molecular Function Enrichment in A-vs-B, A-vs-C, and B-vs-C. Diffmapping (blue), ProteinMapping (orange), DiffNum (gray), and ProteinNum (yellow). DiffMapping is the number of differentially expressed proteins in a GO term, ProteinMapping is the number of quantitative proteins in a GO term, DiffNum is the number of differentially expressed proteins in all GO terms, and ProteinNum is the number of quantitative proteins in all GO terms.

In the A-vs-C group, 146 DEPs were identified—56 upregulated and 90 downregulated (Figure 4b). These DEPs were primarily enriched in structural molecule activity, nucleotide binding, nucleoside phosphate binding, ribonucleoside binding and ribonucleoside binding. For the B-vs-C group, 115 DEPs were detected, including 26 upregulated and 89 downregulated proteins (Figure 4c). The most significantly enriched molecular functions were heterocyclic compound binding and organic cyclic compound binding.

Functional enrichment–based clustering of protein groups was performed across cellular components, molecular functions, biological processes, domains, and KEGG pathways in all three comparison groups.

First, cellular component analysis revealed significant differences in the ribosome and macromolecular complex in the A-vs.-C group (Fisher’s exact p -value, P < 0.01), and in the small ribosomal subunit in the Bvs. C group (P<0.01) (Figure 5a).

Figure 5
Heatmap panels labeled (a) through (e) display gene ontology (GO) and KEGG pathway enrichment. Colors range from green to red, indicating -log10(Fisher's exact p-value) significance levels from less to more significant. Panels illustrate cellular components, molecular functions, biological processes, protein domains, and signaling pathways.

Figure 5. Functional enrichment-based clustering for protein groups in cellular components, molecular function, biological processes, domains, and KEGG clusters of A-vs-B, A-vs-C, and B-vs-C. From red to green (score 0-5), it means the KEGG (P<0.05 or P<0.01) -log10 (Fisher’s exact p -value).

Second, molecular function analysis showed significant enrichment in structural molecule activity, GTPase activity, and structural constituents of ribosomes in the Avs. C group (P < 0.01) (Figure 5b). Third, domain enrichment results indicated significant differences in tubulin-related domains, including the conserved tubulin site, Tubulin/FtsZ, and GTPase domains, also in the Avs. C group (P < 0.01) (Figure 5c). KEGG pathway analysis revealed significant enrichment in carbon metabolism in the A vs. B group (P < 0.05) (Figure 5d). Finally, biological process analysis showed significant enrichment in striated muscle cell differentiation and development, muscle cell development, and muscle structure development in the Bvs. C group (P < 0.01) (Figure 5e).

3.4 Phosphoproteomics and phosphorylated proteins

A total of 2,914 phosphorylation modification sites and 1,039 phosphorylated proteins were identified across the three pairwise comparisons (A-vs-B, A-vs-C, and B-vs-C) (Figures 6a–c). The identified phosphorylation motif sites for serine (Ser), threonine (Thr), and tyrosine (Tyr) residues are shown in Figure 6.

Figure 6
Nine-panel image displaying various charts. Panels a, d, and g are bar charts of GO classification based on graph type. Panels b, e, and h are multi-colored bar charts categorizing biological processes, cellular components, and molecular functions. Panels c, f, and i are doughnut charts showing proportions of cellular components, such as cytoplasm and mitochondrion. Each chart includes a legend indicating color categories.

Figure 6. Total identified phosphorylation proteins. (a-c) Phosphorylation Protein Ratio Distribution of A-vs-B, A-vs-C, and B-vs-C; (d) Identified Phosphorylated Peptides Distribution; (e) Mass Error Distribution; (f) Identified phosphorylation motif sites: Serine (S), Threonine (T), and Tyrosine (Y).

Notably, these phosphorylation sites are involved in the regulation of key signaling pathways, including ubiquitin-mediated proteolysis and Ras signaling. Phosphorylation predominantly occurred on the side chains of Ser, Thr, and Tyr residues in substrate proteins, with 1,558 Ser sites (25 motif types), 1,144 Thr sites (19 motif types), and 209 Tyr sites (8 motif types) identified (Figure 6f).

Differentially expressed phosphorylation proteins (DEPPs) were identified based on a protein quantification threshold of fold change ≥ 1.2 and p-value ≤ 0.05. Functional classification analysis included cellular components, molecular functions, biological processes, and subcellular localization.

DEPPs were primarily involved in signal transduction mechanisms (20.0%–25.0%) and in posttranslational modification, protein turnover, and chaperones (7.69%–15.0%) across the A-vs-B, 547 A-vs-C, and B-vs-C groups (Figures 7a, d, g). Functional classification revealed distinct differences among the three comparison groups. These DEPPs participated in cellular components such as cells, organelles, and membranes (green columns in Figures 7b, e, h) and were functionally associated with binding and catalytic activity (blue columns in Figures 7c, f, i).

Figure 7
Six graphs depicting different data distributions in a scientific analysis. Graphs a, b, and c show A-vs-B, A-vs-C, and B-vs-C peptide ratio distributions, respectively, with varying peptide counts. Graph d illustrates the distribution of protein modify site numbers. Graph e displays a mass error distribution scatter plot with peptide scores against peptide mass error in parts per million. Graph f shows a bar chart of identified phosphorylation peptide motif sites, categorized by type: S, T, Y, and Total, with varying counts.

Figure 7. COG classification of differentially expressed phosphorylation proteins in cellular components, molecular function, biological processes, and subcellular location in three species. (a, d, g) All regulated COG classification in A -vs- B, A -vs- C, B -vs- C, respectively; (b, e, h) The category in cellular components, molecular function, and biological processes in A -vs- B, A vs. C, and B -vs- C, respectively; (c, f, i) The category in subcellular location in A -vs- B, A vs. C, and B -vs- C, respectively.

Biological process classification indicated that most DEPPs were involved in cellular processes and single-organism processes (orange columns). Subcellular localization analysis showed that DEPPs were predominantly located in the nucleus (44.87%–54.90%), mitochondria (8.82%–12.82%), and plasma membrane (11.76%–19.23%). The majority of DEPPs localized to the nucleus and were associated with phosphorylation-mediated signal transduction.

3.5 Phosphoproteomics of DEPs

Functional enrichment–based clustering of phosphorylation-associated DEPs was analyzed across cellular components, 596 molecular functions, biological processes, protein domains, and 597 KEGG pathways in the three comparison groups. Cellular 598 component analysis showed significant enrichment in the cytosol 599 for the A-vs-B group (P < 0.05) and in the microtubule organizing 600 center for the B-vs-C group (P < 0.05) (Figure 8a). Molecular 601 function analysis revealed significant differences in nucleoside-triphosphatase activity, hydrolase activity acting on acid anhydrides, phosphorus-containing anhydrides, and pyrophosphatase activity in the A-vs-B group (P < 0.05) 605 (Figure 8b). Domain enrichment analysis showed significant involvement of P-loop-containing nucleotide-binding domains in the A-vs-B group, and kinesin motor domains and their conserved sites in the B-vs-C group (P < 0.05) (Figure 8c). Biological process analysis indicated significant enrichment in mitochondrial localization in the A-vs-B group, as well as cellular responses to metal ions and inorganic substances (P < 0.01) (Figure 8d). KEGG pathway enrichment revealed significant differences in the Ras signaling pathway in the Avs. B group and in ubiquitin-mediated proteolysis in the B-vs-C group (P < 0.05) (Figure 8e).

Figure 8
Five clustered heat maps representing different analyses with color coding from green to red indicating significance values. Each map depicts relationships among categories based on -log10(Fisher's exact p-value). Categories include cellular components, molecular functions, biological processes, domain structures, and KEGG pathways, with labels such as spindle, kinase binding, response to metal ion, Pleckstrin homology domain, and Ubiquitin mediated proteolysis. Comparisons show analyses between different groups labeled A, B, and C.

Figure 8. Phosphoproteomics cluster in cellular components, molecular function, biological processes, domains, and KEGG cluster of A-vs-B, A-vs-C, and B-vs-C. From red to green (score 0-5), it means the KEGG (P<0.05 or P<0.01) -log10(Fisher’s exact p value). KEGG Annotation: KEGG connects known information on molecular interaction networks, such as pathways and complexes (the “Pathway” database); information about genes and proteins generated by genome projects (including the gene database); and information about biochemical compounds and reactions (including compound and reaction databases).

4 Discussion

The diversity of freshwater mollusks in China has been well documented, including Cipangopaludina chinensis (C. chinensis), Cipangopaludina cathayensis (C. cathayensis), Cipangopaludina yunnaneusis (C. yunnaneusis), Cipangopaludina menglaensis (C. menglaensis), Nodularia douglasiae (N. douglasiae), Bellamya aeruginosa (B.aeruginosa), and Radix auricularia (R. auricularia) (Zhao et al., 2022; Zhou et al., 2022b). Previous studies have shown that specie differences can lead to proteomic variation under specific conditions (Russomanno et al., 2023; Paulus et al., 2022; Grabowsky et al., 2023). In the present study, we compared the proteomic and phosphoproteomic profiles of three freshwater mollusk species. Several novel proteins and phosphorylated proteins were identified and predicted to play key roles in molecular functions, biological processes, and signaling pathways. The aim of this study was to explore differential protein expression among C. chinensis, R.auricularia, and N.douglasiae. A total of 1,382 proteins, 2,914 modified sites, and 1,039 phosphorylated proteins were identified across the pairwise comparisons (C. chinensis vs. R. auricularia, C. chinensis-vs-N douglasiae, and R. auricularia vs. N. douglasiae). Our findings may provide valuable insights into the biological characteristics of these species and their potential susceptibility to microbial infections.

Whether mollusks found in nature can serve as safe and sustainable food sources remains a topic of ongoing discussion. Freshwater gastropods and bivalves have been consumed since at least the Mid–Late Neolithic period (ca. 5600–4500 BP) (Li et al., 2013). Today, species such as C.chinensis and Unio douglasiae (U.douglasiae) are still used as food or feed additives for both humans and animals, commonly found in rivers, reservoirs, rice fields, and experimental aquaculture systems (Li et al., 2013).

However, while freshwater mollusks can serve as a valuable protein source, they may also pose potential health risks. These include microbial co-infections with the host (Chen et al., 2021), latent parasitic infections (Chung et al., 1998; Juhász et al., 2022), and the bioaccumulation of heavy metals (Fang et al., 2001). However, there are still differing opinions on whether naturally sourced or farmed mollusks are more suitable for consumption. Some consumers prefer wild-caught mollusks for their distinctive flavor compared to cultivated varieties.

Notably, our proteomic results microbial metabolism in diverse environments revealed signatures of microbial metabolism in diverse environments across all three species, highlighting their complex interactions with environmental microbes and their potential as both a nutritional source and a model for ecological and food safety studies.

In the context of carbon neutrality, microbes and plants have received significant attention for their ability to store carbon (Pradeep Ram et al., 2016; Li et al., 2023; Song et al., 2024; Reis et al., 2022; Wang et al., 2023), whereas animals have largely been overlooked due to their relatively lower biomass and more limited distribution. However, a growing body of research indicates that the role of animals in regulating greenhouse gas emissions has been underestimated (Schmitz et al., 2018; Clare et al., 2019). According to a report by a research team at Yale University, animals can influence biogeochemical processes by 15% to 250% (Schmitz et al., 2018). In some regions, the scale of carbon uptake or release from specific animal species or groups may even rival fossil fuel emissions in the same region (Schmitz et al., 2018). Our findings revealed significant differences in carbon metabolism–related signaling pathways in the C. chinensis-vs.-R.auricularia comparison group (Gorman et al., 2023; Lind et al., 2017; Nissa et al., 2023). Notably, KEGG enrichment and functional clustering analyses showed that a substantial number of DEPs were involved in carbon metabolism, suggesting that these signaling pathways play a key role in the metabolic processes of freshwater mollusks.

Interestingly, microbial metabolism was also detected in our study, particularly in the C. chinensis-vs.-R. auricularia group. While we cannot fully exclude the possibility of microbial infection in natural hosts, it is important to consider that the growth and developmental environments of freshwater mollusks naturally contain a wide variety of microorganisms, including viruses, bacteria, and parasites (Lavelle and Sokol, 2020; Fiorucci et al., 2021; Gabriel and Ferguson, 2023; Wypych et al., 2021). Both C. chinensis and R. auricularia are univalve mollusks widely distributed in rivers and lakes. Notably, R.auricularia may pose a greater microbial risk in freshwater ecosystems due to its broader environmental exposure and potential to harbor pathogens (Schell et al., 2017; Juhász et al., 2022; Chung et al., 1998).

Phosphorylation plays a critical role in host cell signaling pathway transduction (Li et al., 2022; Chen et al., 2024; Bachtell et al., 2002; Heinrich et al., 1998). Our research showed that majority of DEPPs (20.0%–25.0%) were involved in signal transduction mechanisms. Moreover, these DEPPs were primarily localized to the nucleus (44.87%–54.90%), suggesting that the three freshwater species may preferentially mediate signaling pathways in the nucleus and cytoplasm. Generally, phosphorylation is a common form of post-translational modification (PTM) involving the addition of phosphate groups to proteins, usually catalyzed by kinases. Prior studies have shown that phosphorylation mainly occurs on three amino acid residues: serine (Takai et al., 2018; León-Vaz et al., 2021; Chen et al., 2013; Li et al., 2020), threonine (Hou et al., 2022), and tyrosine (Birkelund et al., 1994; Chen et al., 2017). In our study, several novel DEPPs and phosphorylation sites were identified and quantified in C. chinensis, R.auricularia, and N.douglasiae, which hold significance for understanding regulatory processes in freshwater mollusks.

Previous studies have reported the involvement of ubiquitin mediated proteolysis pathways in freshwater species. For example, transcriptome analysis revealed that the ubiquitin-mediated proteolysis pathway plays a key role in metabolism during testicular differentiation in C. carpio (Zhao et al., 2022). Skeletal muscle protein degradation, mediated through the ubiquitin proteasome system, is also crucial for energy supply, particularly under environmental stress such as crowding (Valenzuela et al., 2020). High-throughput RNA-seq analysis showed that upregulated transcripts in rainbow trout skeletal muscle were predominantly associated with ubiquitin -mediated proteolysis during Flavobacterium psychrophilum infection (Rivas-Aravena et al., 2019). Similarly, in Procambarus clarkii, oxidative stress from heavy metal exposure activated ubiquitination enzymes and proteasomes (Jiao et al., 2019). In our phosphoproteomics analysis, KEGG enrichment results revealed that the ubiquitin -mediated proteolysis pathway plays an important role in the R.auriculariavs. N.douglasiae comparison group, highlighting its relevance in stress response and protein turnover in freshwater mollusks.

In addition to ubiquitin-mediated proteolysis, Ras signaling pathways have also been implicated in freshwater species such as zebrafish, Carassius auratus, Macrobrachium rosenbergii, and Trachemys scripta. For example, angiotensin II within the Ras signaling pathway has been shown to regulate hemodynamic pressure overload –mediated cardiovascular pathogenesis in mammals using the zebrafish model (Joshi et al., 2021). In Macrobrachium rosenbergii, the Ras pathway plays a critical role in maintaining homeostasis in response to hypotonic stress at various developmental stages. Notably, in Trachemys scripta elegans, Ras signaling is involved in regulating spermatogenesis, inflammation, and apoptosis under salinity stress (Li et al., 2024). QTL analysis of alkaline tolerance in C. auratus revealed that Ras signaling may contribute to adaptation to extremely alkaline environments (Zhang et al., 2024). Our results similarly suggest that Ras signaling plays an important role in the C. chinensis vs. R.auricularia comparison group.

While our study presents some novel findings across different species, it is not without limitations. Initially, we intended to analyze the microbial communities present in freshwater mollusk carcasses. However, the volume and diversity of associated viruses, bacteria, and parasites presented a considerable workload, which could not be addressed within the current scope. We plan to explore microbial diversity in these mollusks in future studies.

Some readers may wonder why specific microbes were not tested. This is because the proteomic profiling in this study was conducted using databases relevant to mollusks (specifically snails), not databases of pathogenic microorganisms. Importantly, pathohistological observations revealed no visible tissue damage or presence of pathogenic microbes in the carcass tissues of any of the three mollusk species.

This study was primarily designed to explore the proteomic and phosphoproteomic profiles of C. chinensis, R. auricularia, and N. douglasiae. The results contribute new knowledge by identifying novel proteins and proteases in freshwater mollusks. Furthermore, our research highlights the importance of food safety awareness, particularly in the fishing and catering industries. Finally, bioinformatics analysis indicated that several DEPs and DEPPs may play critical roles in molecular functions, biological processes, and signaling pathways. The identification of these carcass proteins provides a new basis for understanding differential protein expression across natural environments and supports the development of novel strategies for food security in freshwater species.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: http://proteomecentral.proteomex-change.org, IPX0007512000; IPX0009869000; IPX0009870000.

Ethics statement

The manuscript presents research on animals that do not require ethical approval for their study.

Author contributions

GL: Writing – original draft, Writing – review & editing, Funding acquisition, Visualization. TW: Writing – original draft, Conceptualization. SW: Writing – review & editing, Conceptualization. CZ: Writing – review & editing. CL: Writing – review & editing. XZ: Writing – review & editing.

Funding

The author(s) declare financial support was received for there search and/or publication of this article. This study was supported by Linyi University High-level Talent Funding Support (No.Z6122016) and Development and Application of Mollusk Products (No. HX250133).

Acknowledgments

We would like to thank Monitor Helix MH BioTech Co.,Ltd provide for proteomic analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Generative AI was used in the creation of this manuscript. The AI was used for polished manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Bachtell, R. K., Tsivkovskaia, N. O., and Ryabinin, A. E. (2002). Alcohol-induced c-Fos expression in the Edinger-Westphal nucleus: pharmacological and signal transduction mechanisms. J. Pharmacol. Exp. Ther. 302, 516–524. doi: 10.1124/jpet.102.036046

PubMed Abstract | Crossref Full Text | Google Scholar

Bashe, S. K. and Ali, L. A. (2019). First molecular identification and phylogenetic tree of Petasiger exaeretus Dietz 1909 (Digenea: Echinostomatidae) from an intermediate host Radix auricularia (L. 1758) in Greater Zab river, Iraq. Annals Agric. Environ. medicine: AAEM 26, 504–507. doi: 10.26444/aaem/111831

PubMed Abstract | Crossref Full Text | Google Scholar

Beirão, J., Boulais, M., Gallego, V., O’Brien, J. K., Peixoto, S., Robeck, T. R., et al. (2019). Sperm handling in aquatic animals for artificial reproduction. Theriogenology. 133, 161–178. doi: 10.1016/j.theriogenology.2019.05.004

PubMed Abstract | Crossref Full Text | Google Scholar

Birkelund, S., Johnsen, H., and Christiansen, G. (1994). Chlamydia trachomatis serovar L2 induces protein tyrosine phosphorylation during uptake by HeLa cells. Infection Immun. 62, 4900–4908. doi: 10.1128/iai.62.11.4900-4908.1994

PubMed Abstract | Crossref Full Text | Google Scholar

Chakraborty, K. and Joy, M. (2020). High-value compounds from the mollusks of marine and estuarine ecosystems as prospective functional food ingredients: An overview. Food Res. Int. (Ottawa Ont.) 137, 109637. doi: 10.1016/j.foodres.2020.109637

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, C. L., Chien, S. C., Leu, T. H., Harn, H. I., Tang, M. J., and Hor, L. I. (2017). Vibrio vulnificus MARTX cytotoxin causes inactivation of phagocytosis-related signaling molecules in macrophages. J. Biomed. Sci. 24, 58. doi: 10.1186/s12929-017-0368-2

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, D., Li, X., Ni, L., Xu, D., Xu, Y., Ding, Y., et al. (2021). First experimental evidence for the presence of potentially toxic vibrio cholerae in snails, and virulence, cross-resistance and genetic diversity of the bacterium in 36 species of aquatic food animals. Antibiotics (Basel Switzerland) 10, 412. doi: 10.3390/antibiotics10040412

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Liang, H., Wu, Y., and Li, C. (2024). Phosphoproteomics changes due to allograft-induced stress responses of Pinctada fucata martensii. Comparative biochemistry and physiology. Part D Genomics Proteomics 49, 101153. doi: 10.1016/j.cbd.2023.101153

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, L., Wu, J., Pier, E., Zhao, Y., and Shen, Z. (2013). ). mTORC2-PKBα/Akt1 Serine 473 phosphorylation axis is essential for regulation of FOXP3 Stability by chemokine CCL3 in psoriasis. J. Invest. Dermatol. 133, 418–428. doi: 10.1038/jid.2012.333

PubMed Abstract | Crossref Full Text | Google Scholar

Choi, E. H., Kim, G., Cha, S. H., Lee, J. S., Ryu, S. H., Suk, H. Y., et al. (2020). Molecular phylogenetic, population genetic and demographic studies of Nodularia douglasiae and Nodularia breviconcha based on CO1 and 16S rRNA. Sci. Rep. 10, 16572. doi: 10.1038/s41598-020-72015-5

PubMed Abstract | Crossref Full Text | Google Scholar

Chung, D. I., Kong, H. H., and Joo, C. Y. (1998). Radix auricularia coreana: natural snail host of Clinostomum complanatum in Korea. Korean J. Parasitol. 36, 1–6. doi: 10.3347/kjp.1998.36.1.1

PubMed Abstract | Crossref Full Text | Google Scholar

Clare, C. E., Brassington, A. H., Kwong, W. Y., and Sinclair, K. D. (2019). One-carbon metabolism: linking nutritional biochemistry to epigenetic programming of long-term development. Annu. Rev. Anim. Biosci. 7, 263–287. doi: 10.1146/annurev-animal-020518-115206

PubMed Abstract | Crossref Full Text | Google Scholar

Darrigran, G., Belz, C., Carranza, A., Collado, G. A., Correoso, M., Fabres, A. A., et al. (2025). What do we know about non-native, invasive, and transplanted aquatic mollusks in south america? Biol. (Basel). 14, 151. doi: 10.3390/biology14020151

PubMed Abstract | Crossref Full Text | Google Scholar

Deidda, I., Russo, R., Bonaventura, R., Costa, C., Zito, F., and Lampiasi, N. (2021). Neurotoxicity in marine invertebrates: an update. Biology 10, 161. doi: 10.3390/biology10020161

PubMed Abstract | Crossref Full Text | Google Scholar

de la Ballina, N. R., Maresca, F., Cao, A., and Villalba, A. (2022). Bivalve haemocyte subpopulations: A review. Front. Immunol. 13. doi: 10.3389/fimmu.2022.826255

PubMed Abstract | Crossref Full Text | Google Scholar

Elder, J. F. and Collins, J. J. (1991). Freshwater mollusks as indicators of bioavailability and toxicity of metals in surface-water systems. Rev. Environ. contamination Toxicol. 122, 37–79. doi: 10.1007/978-1-4612-3198-1_2

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, Z. Q., Cheung, R. Y., and Wong, M. H. (2001). Heavy metal concentrations in edible bivalves and gastropods available in major markets of the Pearl River Delta. J. Environ. Sci. (China) 13. doi: 10.3389/fmicb.2016.00081

PubMed Abstract | Crossref Full Text | Google Scholar

Fernandes, D., Loi, B., and Porte, C. (2011). Biosynthesis and metabolism of steroids in mollusks. J. Steroid Biochem. Mol. Biol. 127, 189–195. doi: 10.1016/j.jsbmb.2010.12.009

PubMed Abstract | Crossref Full Text | Google Scholar

Fiorucci, S., Carino, A., Baldoni, M., Santucci, L., Costanzi, E., Graziosi, L., et al. (2021). Bile acid signaling in inflammatory bowel diseases. Digestive Dis. Sci. 66, 674–693. doi: 10.1007/s10620-020-06715-3

PubMed Abstract | Crossref Full Text | Google Scholar

Gabriel, C. L. and Ferguson, J. F. (2023). Gut microbiota and microbial metabolism in early risk of cardiometabolic disease. Circ. Res. 132, 1674–1691. doi: 10.1161/CIRCRESAHA.123.322055

PubMed Abstract | Crossref Full Text | Google Scholar

Gaudin-Zatylny, C., Corre, E., Zanuttini, B., Endress, M., Bernay, B., Pontin, J., et al. (2022). Identification of a new set of polypeptidic sex pheromones from cuttlefish (Sepia officinalis). Mar. Biotechnol. (NY). 24, 574–587. doi: 10.1007/s10126-022-10126-y

PubMed Abstract | Crossref Full Text | Google Scholar

Giansante, C. and Pelini, S. (2007). The use of geographic information systems in sea and freshwater ecosystems. Veterinaria italiana 43, 507–512.

PubMed Abstract | Google Scholar

Gorman, D., Beale, D. J., Crosswell, J., Stephenson, S. A., Shah, R. M., Hillyer, K. E., et al. (2023). Multiple-biomarkers show the importance of blue carbon to commercially important fishery species. Sci. total Environ. 881, 163162. doi: 10.1016/j.scitotenv.2023.163162

PubMed Abstract | Crossref Full Text | Google Scholar

Grabowsky, E. R., Saviola, A. J., Alvarado-Díaz, J., Mascareñas, A. Q., Hansen, K. C., Yates, J. R., 3rd, et al. (2023). Montane Rattlesnakes in México: Venoms of Crotalus tancitarensis and Related Species within the Crotalus intermedius Group. Toxins 15, 72. doi: 10.3390/toxins15010072

PubMed Abstract | Crossref Full Text | Google Scholar

Heinrich, P. C., Behrmann, I., Müller-Newen, G., Schaper, F., and Graeve, L. (1998). Interleukin-6-type cytokine signalling through the gp130/Jak/STAT pathway. Biochem. J. 334, 297–314. doi: 10.1042/bj3340297

PubMed Abstract | Crossref Full Text | Google Scholar

Hou, J., Dai, J., Chen, Z., Wang, Y., Cao, J., Hu, J., et al. (2022). Phosphorylation regulation of a histone-like HU protein from deinococcus radiodurans. Protein Pept. Lett. 29, 891–899. doi: 10.2174/0929866529666220819121911

PubMed Abstract | Crossref Full Text | Google Scholar

Jiao, T., Chu, X. H., Gao, Z. Q., Yang, T. T., Liu, Y., Yang, L., et al. (2019). New insight into the molecular basis of Fe (III) stress responses of Procambarus clarkii by transcriptome analysis. Ecotoxicology Environ. Saf. 182, 109388. doi: 10.1016/j.ecoenv.019.109388

PubMed Abstract | Crossref Full Text | Google Scholar

Joshi, B., Wagh, G., Kaur, H., and Patra, C. (2021). Zebrafish model to study angiotensin II-mediated pathophysiology. Biology 10, 1177. doi: 10.3390/biology10111177

PubMed Abstract | Crossref Full Text | Google Scholar

Juhász, A. and Lawton, S. P. (2022). Toll like receptors and their evolution in the lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts for zoonotic trematodes. Dev. Comp. Immunol. 127, 104297. doi: 10.1016/j.dci.2021.104297

PubMed Abstract | Crossref Full Text | Google Scholar

Juhász, A., Majoros, G., and Cech, G. (2022). Threat of cercarial dermatitis in Hungary: A first report of Trichobilharzia franki from the mallard (Anas platyrhynchos) and European ear snail (Radix auricularia) using molecular methods.International journal for parasitology. Parasites wildlife 18, 92–100. doi: 10.1016/j.ijppaw.2022.04.009

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, Y. R., Choi, H. K., Cheon, K. S., and Lee, H. J. (2021). Complete mitochondrial genome of a Korean endemic freshwater mussel Nodularia breviconcha (Bivalvia: Unionidae). Mitochondrial DNA. Part B Resour. 6, 79–81. doi: 10.1080/23802359.2020.1846473

PubMed Abstract | Crossref Full Text | Google Scholar

Lavelle, A. and Sokol, H. (2020). Gut microbiota-derived metabolites as key actors in inflammatory bowel disease. Nature reviews. Gastroenterol. Hepatol. 17, 223–237. doi: 10.1038/s41575-019-0258-z

PubMed Abstract | Crossref Full Text | Google Scholar

León-Vaz, A., Romero, L. C., Gotor, C., León, R., and Vigara, J. (2021). Effect of cadmium in the microalga Chlorella sorokiniana: A proteomic study. Ecotoxicology Environ. Saf. 207, 111301. doi: 10.1016/j.ecoenv.2020.111301

PubMed Abstract | Crossref Full Text | Google Scholar

Li, W., Luo, D., Adyel, T. M., Wu, J., Miao, L., and Hou, J. (2023). Dynamic responses of carbon metabolism of sediment microbial communities to Ag nanoparticles: Effects of the single and repeated exposure scenarios. Sci. total Environ. 870, 161891. doi: 10.1016/j.scitotenv.2023.161891

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Luo, Z., Zhang, C., Qu, X., Chen, M., Song, T., et al. (2020). Seasonal variation in the rhizosphere and non-rhizosphere microbial community structures and functions of camellia yuhsienensis hu. Microorganisms 8, 1385. doi: 10.3390/microorganisms8091385

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Tang, H., Zhao, W., Yang, Y., Fan, X., Zhan, G., et al. (2022). Study of dimorphism transition mechanism of tremella fuciformis based on comparative proteomics. J. fungi (Basel Switzerland) 8, 242. doi: 10.3390/jof8030242

PubMed Abstract | Crossref Full Text | Google Scholar

Li, F., Wu, N., Lu, H., Zhang, J., Wang, W., Ma, M., et al. (2013). Mid-Neolithic exploitation of mollusks in the Guanzhong Basin of Northwestern China: preliminary results. PloS One 8, e58999. doi: 10.1371/journal.pone.0058999

PubMed Abstract | Crossref Full Text | Google Scholar

Li, N., Zhang, Q., Dai, S., Rao, W., Shi, H., Ding, L., et al. (2024). Angiotensin-(1-7) plays an important role in regulating spermatogenesis in Trachemys scripta elegans under salinity stress. J. Exp. Biol. 227, jeb246742. doi. doi: 10.1242/jeb.246742

PubMed Abstract | Crossref Full Text | Google Scholar

Lind, M. V., Lauritzen, L., Pedersen, O., Vestergaard, H., Stark, K. D., Hansen, T., et al. (2017). Higher intake of fish and fat is associated with lower plasma s-adenosylhomocysteine: a cross-sectional study. Nutr. Res. (New York N.Y.) 46, 78–87. doi: 10.1016/j.nutres.2017.09.008

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Cao, Y., Xue, T., Wu, R., Zhou, Y., Zhou, C., et al. (2017). Genetic structure and diversity of Nodularia douglasiae (Bivalvia: Unionida) from the middle and lower Yangtze River drainage. PloS One 12, e0189737. doi: 10.1371/journal.pone.0189737

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Liu, C., and Huang, J. (2023). Characterization of the shell proteins in two freshwater snails Pomacea canaliculata and Cipangopaludina chinensis. Int. J. Biol. macromolecules 1), 124524. doi: 10.1016/j.ijbiomac.2023.124524

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, G., Miao, F., Wang, Y., Kou, J., Yang, K., Li, W., et al. (2022). Comparative proteomics analysis of Schistosoma japonicum developed in different Oncomelania snails as intermediate hosts. Front. Cell. infection Microbiol. 12. doi: 10.3389/fcimb.2022.959766

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, G., Zhong, K., Gong, S., Li, X., and Li, Y. (2024). Proteomics and phosphoproteomics reveal novel proteins involved in Cipangopaludina chinensis carcasses. Front. Chem. 12. doi: 10.3389/fchem.2024.1416942

PubMed Abstract | Crossref Full Text | Google Scholar

Mabuchi, K., Nishida, K., and Nakajima, N. (2021). Complete female-transmitted mitochondrial genomes of two freshwater mussels from the Lake Biwa system in Japan: Nodularia douglasiae and N. nipponensis. Mitochondrial DNA. Part B Resour. 6, 1536–1538. doi: 10.1080/23802359.2021.1914218

PubMed Abstract | Crossref Full Text | Google Scholar

Marrone, F., Nardi, G., Cianfanelli, S., GovediČ, M., Barra, S. A., Arculeo, M., et al. (2019). Diversity and taxonomy of the genus Unio Philipsson in Italy, with the designation of a neotype for Unio elongatulus C. Pfeiffer (molluska, Bivalvia, Unionidae). Zootaxa 4545, 339–374. doi: 10.11646/zootaxa.4545.3.2

PubMed Abstract | Crossref Full Text | Google Scholar

Mienis, H. K. and Ashkenazi, S. (2011). Lentic Basommatophora mollusks and hygrophilous land snails as indicators of habitat and climate in the Early-Middle Pleistocene (0.78 Ma) at the site of Gesher Benot Ya’aqov (GBY), Israel. J. Hum. Evol. 60, 328–340. doi: 10.1016/j.jhevol.2010.03.009

PubMed Abstract | Crossref Full Text | Google Scholar

Nikishchenko, V. E., Sayenko, E. M., and Dyachuk, V. A. (2022). First Immunodetection of Sensory and Nervous Systems of Parasitic Larvae (Glochidia) of Freshwater Bivalve Nodularia douglasiae. Frontiers in physiology. 13, 879540. doi: 10.3389/fphys.2022.879540

PubMed Abstract | Crossref Full Text | Google Scholar

Nissa, M. U., Pinto, N., Ghosh, B., Singh, U., Goswami, M., and Srivastava, S. (2023). Proteomic analysis of liver tissue reveals Aeromonas hydrophila infection mediated modulation of host metabolic pathways in Labeo rohita. J. Proteomics 279, 104870. doi: 10.1016/j.jprot.2023.104870

PubMed Abstract | Crossref Full Text | Google Scholar

Paulus, E., Brix, S., Siebert, A., Martínez Arbizu, P., Rossel, S., Peters, J., et al. (2022). Recent speciation and hybridization in Icelandic deep-sea isopods: An integrative approach using genomics and proteomics. Mol. Ecol. ,31, 313–330. doi: 10.1111/mec.16234

PubMed Abstract | Crossref Full Text | Google Scholar

Pradeep Ram, A. S., Colombet, J., Perriere, F., Thouvenot, A., and Sime-Ngando, T. (2016). Viral regulation of prokaryotic carbon metabolism in a hypereutrophic freshwater reservoir ecosystem (Villerest, France). Front. Microbiol. 7. doi: 10.3389/fmicb.2016.00081

PubMed Abstract | Crossref Full Text | Google Scholar

Reis, P. C. J., Thottathil, S. D., and Prairie, Y. T. (2022). The role of methanotrophy in the microbial carbon metabolism of temperate lakes. Nat. Commun. 13, 43. doi: 10.1038/s41467-021-27718-2

PubMed Abstract | Crossref Full Text | Google Scholar

Rivas-Aravena, A., Fuentes-Valenzuela, M., Escobar-Aguirre, S., Gallardo-Escarate, C., Molina, A., and Valdés, J. A. (2019). Transcriptomic response of rainbow trout (Oncorhynchus mykiss) skeletal muscle to Flavobacterium psychrophilum. Comparative biochemistry and physiology. Part D Genomics Proteomics 31, 100596. doi: 10.1016/j.cbd.2019.100596

PubMed Abstract | Crossref Full Text | Google Scholar

Russomanno, G., Sison-Young, R., Livoti, L. A., Coghlan, H., Jenkins, R. E., Kunnen, S. J., et al. (2023). A systems approach reveals species differences in hepatic stress response capacity. Toxicological sciences: an Off. J. Soc. Toxicol. 196, 112–125. doi: 10.1093/toxsci/kfad085

PubMed Abstract | Crossref Full Text | Google Scholar

Schell, T., Feldmeyer, B., Schmidt, H., Greshake, B., Tills, O., Truebano, M., et al. (2017). An annotated draft genome for Radix auricularia (Gastropoda, molluska). Genome Biol. Evol. 9, 0. doi: 10.1093/gbe/evx032

PubMed Abstract | Crossref Full Text | Google Scholar

Schmitz, O. J., Wilmers, C. C., Leroux, S. J., Doughty, C. E., Atwood, T. B., Galetti, M., et al. (2018). Animals and the zoogeochemistry of the carbon cycle. Sci. (New York N.Y.) 362, eaar3213. doi: 10.1126/science.aar3213

PubMed Abstract | Crossref Full Text | Google Scholar

Song, Y., Cao, X., Li, S. A., Li, Z., Grossart, H. P., and Ma, H. (2024). Human activities-impacted lake dissolved organic matter (DOM) affects phycosphere microbial diversity and DOM diversification via carbon metabolism. J. Environ. Manage. 367, 122011. doi: 10.1016/j.jenvman.2024.122011

PubMed Abstract | Crossref Full Text | Google Scholar

Takai, A., Eto, M., Hirano, K., Takeya, K., Wakimoto, T., and Watanabe, M. (2018). Protein phosphatases 1 and 2A and their naturally occurring inhibitors: current topics in smooth muscle physiology and chemical biology. J. Physiol. sciences: JPS ,68, 1–17. doi: 10.1007/s12576-017-0556-6

PubMed Abstract | Crossref Full Text | Google Scholar

Valenzuela, C. A., Ponce, C., Zuloaga, R., González, P., Avendaño-Herrera, R., Valdés, J. A., et al. (2020). Effects of crowding on the three main proteolytic mechanisms of skeletal muscle in rainbow trout (Oncorhynchus mykiss). BMC veterinary Res. 16, 294. doi: 10.1186/s12917-020-02518-w

PubMed Abstract | Crossref Full Text | Google Scholar

Vasta, G. R., Feng, C., Bianchet, M. A., Bachvaroff, T. R., and Tasumi, S. (2015). Structural, functional, and evolutionary aspects of galectins in aquatic mollusks: From a sweet tooth to the Trojan horse. Fish Shellfish Immunol. 46, 94–106. doi: 10.1016/j.fsi.2015.05.012

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Bai, J., Wen, L., Wang, W., Zhang, L., Liu, Z., et al. (2023). Phytotoxicity of microplastics to the floating plant Spirodela polyrhiza (L.): Plant functional traits and metabolomics. Environ. pollut. (Barking Essex: 1987) 322, 121199. doi: 10.1016/j.envpol.2023.121199

PubMed Abstract | Crossref Full Text | Google Scholar

Wayne, N. L. (2001). Regulation of seasonal reproduction in mollusks. J. Biol. Rhythms. 16, 391–402. doi: 10.1177/074873001129002097

PubMed Abstract | Crossref Full Text | Google Scholar

Wypych, T. P., Pattaroni, C., Perdijk, O., Yap, C., Trompette, A., Anderson, D., et al. (2021). Microbial metabolism of L-tyrosine protects against allergic airway inflammation. Nat. Immunol. 22, 279–286. doi: 10.1038/s41590-020-00856-3

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, L., Su, B., Huang, J., Zhang, L., Chang, Y., and Hu, G. (2024). Fine mapping of QTLs for alkaline tolerance in crucian carp (Carassius auratus) using genome-wide SNP markers. Genes 15, 751. doi: 10.3390/genes15060751

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Chen, K., Liu, F., Jiang, M., Chen, Z., Chen, H., et al. (2022). Dynamic gene expression and alternative splicing events demonstrate co-regulation of testicular differentiation and maturation by the brain and gonad in common carp. Front. Endocrinol. 12. doi: 10.3389/fendo.2021.820463

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, K. M., Liu, Z. H., Zhang, J., Zhong, S. S., and Dang, Z. (2022). Property of arylsulfatase and β-glucuronidase extracted from digestive tracts of Cipangopaludina chinensis and their cleavage performance on conjugated natural estrogens. Environ. Sci. pollut. Res. Int. 29, 64244–64251. doi: 10.1007/s11356-022-22260-0

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, K., Chen, Z., Du, X., Huang, Y., Qin, J., Wen, L., et al. (2022a). SMRT sequencing reveals candidate genes and pathways with medicinal value in cipangopaludina chinensis. Front. Genet. 13. doi: 10.3389/fgene.2022.881952

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, K., Qin, J., Pang, H., Chen, Z., Huang, Y., Li, W., et al. (2022b). Comparison of the composition and function of gut microbes between adult and juvenile Cipangopaludina chinensis in the rice snail system. PeerJ 10, e13042. doi: 10.7717/peerj.13042

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: freshwater, mollusks carcasses, proteomics, phosphoproteomics, cellular components, molecular functions, biological processes, differentially expressed proteins (DEPs)

Citation: Liu G, Wang S, Wang T, Zhu C, Li C and Zhang X (2025) Proteomics and phosphoproteomics of freshwater mollusk carcasses reveal novel insights as potential food source. Front. Cell. Infect. Microbiol. 15:1568349. doi: 10.3389/fcimb.2025.1568349

Received: 29 January 2025; Accepted: 19 June 2025;
Published: 14 August 2025.

Edited by:

Shahid Karim, University of Southern Mississippi, United States

Reviewed by:

Duminda Senevirathna, Uva Wellassa University, Sri Lanka
Amna Bajwa, University of Veterinary and Animal Sciences, Pakistan

Copyright © 2025 Liu, Wang, Wang, Zhu, Li and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Gongzhen Liu, Z29uZ3poZW5saXVAMTI2LmNvbQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.