Skip to main content

Protein expression patterns and metal metabolites in a protogynous hermaphrodite fish, the ricefield eel (Monopterus albus)

Abstract

Background

The ricefield eel Monopterus albus undergoes a natural sex change from female to male during its life cycle, and previous studies have shown the potential mechanisms of this transition at the transcriptional and protein levels. However, the changes in protein levels have not been fully explored, especially in the intersexual stage.

Results

In the present study, the protein expression patterns in the gonadal tissues from five different periods, the ovary (OV), early intersexual stage gonad (IE), middle intersexual stage gonad (IM), late intersexual stage gonad (IL), and testis (TE), were determined by untargeted proteomics sequencing. A total of 5125 proteins and 394 differentially expressed proteins (DEPs) were detected in the gonadal tissues. Of the 394 DEPs, there were 136 between the OV and IE groups, 20 between the IM and IE groups, 179 between the IL and IM groups, and 59 between the TE and IL groups. Three candidate proteins, insulin-like growth factor 2 mRNA-binding protein 3 isoform X1 (Igf2bp3), triosephosphate isomerase (Tpi), and Cu-Zn superoxide dismutase isoform X1 [(Cu-Zn) Sod1], were validated by western blotting to verify the reliability of the data. Furthermore, metal metabolite-related proteins were enriched in the IL vs. IM groups and TE vs. IL groups, which had close relationships with sex change, including Cu2+-, Ca2+-, Zn2+- and Fe2+/Fe3+-related proteins. Analysis of the combined transcriptome data revealed consistent protein/mRNA expression trends for two metal metabolite-related proteins/genes [LOC109953912 and calcium Binding Protein 39 Like (cab39l)]. Notably, we detected significantly higher levels of Cu2+ during the sex change process, suggesting that Cu2+ is a male-related metal metabolite that may have an important function in male reproductive development.

Conclusions

In summary, we analyzed the protein profiles of ricefield eel gonadal tissues in five sexual stages (OV, IE, IM, IL, and TE) and verified the plausibility of the data. After preforming the functional enrichment of metal metabolite-related DEPs, we detected the contents of the metal metabolites Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+ at these five stages and screened for (Cu-Zn) Sod1 and Mmp-9 as possible key proteins in the sex reversal process.

Peer Review reports

Background

Sex determination exhibits significant diversity and variability among animals, where it is often thought to be driven by so-called master sex-determining genes or multiple environmental factors [1,2,3,4]. Some fish have a long history of undergoing sex change, the process of transitioning from one sex to another [5,6,7]. Thus, hermaphroditic fish provide an interesting model for the study of the mechanisms of sex determination and differentiation in vertebrates. There are three kinds of hermaphroditic fish: protogynous (female-to-male sex change), such as the ricefield eel Monopterus albus [8], Lutjanus campechanus [9], wrasse (blue-headed wrasse Thalassoma bifasciatum and New Zealand spotted wrasse Notolabrus celidotus [10]; protandry (male-to-female sex change), e.g., Amphiprion bicinctus [11] and the black sea bream (Acanthopagrus schlegelii) [12]; and bidirectional, such as two genera of small gobies (Gobiodon and Paragobiodon) [13] and the chalk bass Serranus tortugarum [14].

The ricefield eel is an economically important fish in China and Southeast Asia and is also a good model species for studying sex change, as they naturally change from female to male [15]. After reaching sexual maturity, the ricefield eel is a female first, and after laying eggs, the ovary gradually changes into a testis. The sex change process of the ricefield eel includes five periods: ovary (OV), early intersexual stage gonad (IE), middle intersexual stage gonad (IM), late intersexual stage gonad (IL), and testis (TE) [16,17,18]. Based on the transcriptome data, a few genes involved in the regulation of the sex change in the ricefield eel, such as transforming growth factor beta 3 (tgfb3), follicle stimulating hormone receptor (fshr), calcium/calmodulin dependent protein kinase 4 (camk4) and calmodulin may be the key genes, which helps balance germ cell apoptosis and proliferation during sex change [19]. Most differentially expressed genes (DEGs; mRNAs, circRNAs, and lncRNAs) were primarily enriched in the intersexual stages (IE, IL, and IM) [19,20,21]. The DEGs are involved in endocytosis, autophagy, and p53 apoptosis signaling [19]. In addition, two-dimensional electrophoresis was used to analyze differentially expressed proteins (DEPs) during the three phases of sex change in the ricefield eel, including the ovary, ovotestis and testis [22]. DEPs in the ovotestis phase have two functions related to ovarian apoptosis and testicular differentiation [22]. Sex change is a complex biological process [3, 6]. Although the process of sex change are now described for involving in the behavioral, gonadal, and morphological modifications, the genetic cascade orchestrating this transformation needs deeper exploration [23]. At the transcriptional level, many genes were differentially expressed during the sex change in ricefield eel, especially during the IE, IM, and IL stages, such as forkhead box L2 (foxl2) and wnt family member 4 (wnt4) expression was high in OV and IE gonad, and SRY-box transcription factor 9 (sox9), doublesex and mab-3 related transcription factor 1 (dmrt1), dmrt2 and dmrt3 in the TE stage were more higher than those in the OV and intersexual stages [19]. The determination of the proteome in the previous study could not completely reveal the protein patterns during the sex change process [22].

To date, the biological process of sex change is not very clear. Therefore, it is necessary to comprehensively uncover the protein expression profiles during the sex change period. In our study, the protein profiles in gonadal tissues of ricefield eel in five sexual stages (OV, IE, IM, IL, and TE) were determined and quantified. Then, the functional enrichment of DEPs was analyzed, and the metal metabolites-related proteins were identified. Furthermore, western blotting was utilized to detect the expression levels of the DEPs. The contents of metal metabolites (Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+) were detected in tissues in these five sexual stages. The results here lay a foundation for the study of functional proteins during the process of sex change in ricefield eel and further provide references for the mechanism of sex change and fish hermaphroditism.

Results

Protein profiles during sex change

According to the histological changes in the gonads, the OV, IE, IM, IL, and TE stages were identified. The protein patterns in the five sexual stages were detected by iTRAQ mass spectrometry. The peptide spectrum matches (PSMs) were more than 95% credible, and each credible protein contained at least one unique peptide. Only credible PSMs and proteins were retained, and false discovery rate (FDR) verification was performed to remove peptides and proteins with FDR > 1%. The main peptide length was 4–20 amino acids (aa), and a peak was observed at 12 aa (Fig. 1A). The sizes of the enriched proteins were mainly 10–70 kDa with a peak at 20 kDa (Fig. 1B). The pros and cons of repeatability were evaluated by determining the coefficient of variance (CV) value. The CV values in the IL and IE periods were the smallest, showing good repeatability. The TE values were the largest, and the reproducibility was relatively poor (Fig. 1C).

A total of 5125 nonredundant proteins were identified through the five sexual stages (Additional file 2–3), including 4422 proteins in OV, 4430 proteins in IE, 4679 proteins in IM, 4196 proteins in IL, and 4868 proteins in TE. Only 3499 common proteins were detected in the five sexual stages (Fig. 1D). Four proteins were selected to verify the accuracy of the sequencing data by western blot (Fig. 1E-G and Additional file 1: Fig. S1p < 0.05), including V-type proton ATPase subunit G 1 (Vatpg1), triosephosphate isomerase (Tpi), insulin-like growth factor 2 mRNA-binding protein 3 isoform X1 (Igf2bp3), and Cu-Zn superoxide dismutase isoform X1 [(Cu-Zn) Sod1]. Among them, three proteins [Tpi, Igf2bp3, and (Cu-Zn) Sod1] had similar expression patterns in the sequencing results.

Fig. 1
figure 1

Characterization and identification of total proteins and western blot verification of four candidate proteins. A, peptide length. B, protein mass. C, coefficient of variation (CV). D, overlapping proteins. E, western blot verification. F, proteome data. G, western blot verification (grayscale processing). OV, ovary; IE, early intersex stage gonad; IM, middle intersex stage gonad; IL, late intersex stage gonad; TE, testis. The results are presented as the means ± SEMs. Different lowercase superscripted letters indicate significant differences (p < 0.05)

Verification and analysis of the DEPs during sex change

The clustering of the DEPs showed that the five sexual stages, OV, IE, IM, IL, and TE, exhibited distinct expression patterns (Fig. 2A). A total of 394 DEPs were screened. The most DEPs were found in the IL vs. IM group comparison (101 downregulated and 78 upregulated in IL, Fig. 2D and Additional file 7), followed by the IE vs. OV groups (101 upregulated and 35 downregulated in IE, Fig. 2B and Additional file 4) and TE vs. IL groups (31 downregulated and 28 upregulated in TE, Fig. 2E and Additional file 8). Twenty DEPs were found in the IM vs. IE groups (10 downregulated and 10 upregulated in IM, Fig. 2C and Additional file 6). However, there were no common DEPs in the four comparison groups (IE vs. OV, IM vs. IE, and IL vs. IM, and TE vs. IL). Two DEPs, noelin isoform X1 and pentraxin fusion protein-like, were significantly differentially expressed in three groups (IE vs. OV, IM vs. IE, and IL vs. IM). 19 DEPs were screened in two groups, including three DEPs between the IE vs. OV and IM vs. IE groups, 10 DEPs between the IE vs. OV and IL vs. IM groups, and 6 DEPs between the IE vs. OV and TE vs. IL groups (Fig. 2F).

Fig. 2
figure 2

Differentially expressed proteins (DEPs) during sex change. A, clustering of DEGs. B-E, volcano maps of DEPs (B, IE vs. OV. C, IM vs. IE. D, IL vs. IM. E, TE vs. IL). F, overlapping DEPs. OV, ovary; IE, early intersex stage gonad; IM, middle intersex stage gonad; IL, late intersex stage gonad; TE, testis

Functional enrichment of total proteins

The functional enrichment of the total 5125 proteins was carried out. A total of 1540 GO terms were enriched, including 699 terms in biological process, 235 terms in cellular component, and 604 terms in molecular function. The top five terms were enriched in protein binding (972 proteins), ATP binding (379 proteins), nucleic acid binding (251 proteins), oxidation‒reduction process (200 proteins), and zinc ion binding (183 proteins) (Fig. 3A and Additional file 9).

COG analysis suggested that 5125 proteins were enriched in a total of 25 classes, and the top five classes with the most abundant proteins were general function prediction only (422 proteins), posttranslational modification, protein turnover, chaperones (338 proteins), translation, ribosomal structure and biogenesis (331 proteins), signal transduction mechanisms (236 proteins), and carbohydrate transport and metabolism (148 proteins) (Fig. 3B and Additional file 10).

Furthermore, the top 10 enriched KEGG pathways were infectious diseases (604 proteins), cancers (429 proteins), signal transduction (595 proteins), transport and catabolism (398 proteins), immune system (394 proteins), translation, folding sorting and degradation (333 proteins), endocrine system (313 proteins), cellular community (289 proteins), carbohydrate metabolism (264 proteins), and neurodegenerative diseases (234 proteins) (Fig. 3C and Additional file 11).

In general, the overlapping functions from the GO, COG, and KEGG analyses were mainly concentrated in energy metabolism, signal transduction, and immune infection.

Fig. 3
figure 3

GO, COG, and KEGG analyses of proteins during sex change in the ricefield eel Monopterus albus. A, GO analysis. B, COG analysis. C, KEGG analysis. GO, Gene Ontology; COG, Clusters of Orthologous Groups; KEGG, Kyoto Encyclopedia of Genes and Genomes

GO analysis of sex-related DEPs

To analyze the functions of the DEPs during sex change, GO enrichment analysis was performed. The IL vs. IM group had the greatest number of GO terms (220 terms), with the most relevant terms in energy metabolism and ion metabolism, for example, ion binding, catabolic process, carbohydrate metabolic process, oxoacid metabolic process, and carboxylic acid metabolic process, etc. (Fig. 4A and Additional file 14). The IE vs. OV groups had 159 enriched terms, and the function of most the DEPs were related to posttranscriptional protein modifications, such as cytoplasm, lipid transporter activity, and lipid transport (Fig. 4A and Additional file 12). The TE vs. IL groups had 94 enriched terms, such as RNA binding, RNA processing and cell macromolecular complex assembly. Similarly, between the IL group and IM groups, there were abundant enrichments in the biological processes of ion metabolism and energy metabolism (Fig. 4A and Additional file 15). The DEPs in the IM vs. IE groups had the fewest enriched entries, and the most significant enrichment term was RNA binding (Fig. 4A and Additional file 13). Only one overlapping term was enriched in all four groups, which was immune response (Fig. 4B).

Importantly, the most abundant DEPs in the IM vs. IL groups were enriched in ion metabolism, such as ion binding, anion binding, metal ion homeostasis, and cellular metal ion homeostasis. These results suggest that ion metabolism may have important roles in the sex change process of ricefield eel. GO-Slim molecular function analysis was used to determine the functional classification. The main types of metal metabolites-related terms were molecular functions and biological processes. To further explore the potential functions of metal metabolite-related DEPs in sex change, we screened the relevant terms of four metal metabolites according to the GO database: Cu2+, Ca2+, Zn2+ and Fe2+/Fe3+ (Table 1). The DEPs of these metal metabolites were mainly enriched in the IM vs. IL group, with the calcium and iron terms being the most abundant. Cu2+ was enriched in only the IL vs. IM group, and Zn2+ had only one term for zinc ion binding. These data also suggest that metal metabolites are mainly related with the IM vs. IL comparison group, with an important role in male development.

Table 1 Metal metabolites-related terms in the ricefield eel Monopterus albus
Fig. 4
figure 4

GO and KEGG analyses of differentially expressed proteins (DEPs) during sex change in the ricefield eel Monopterus albus. A, significantly enriched GO terms. B, overlapping terms. C, significantly enriched KEGG pathways. D, overlapping KEGG pathways. OV, ovary; IE, early intersex stage gonad; IM, middle intersex stage gonad; IL, late intersex stage gonad; TE, testis. MF, molecular function; CC, cellular component; BP, biological process

KEGG enrichment of sex-related DEPs

KEGG enrichment analysis was performed on the DEPs, and significant terms (p < 0.05) were screened (Fig. 4C). Among them, the pathways with significant enrichment were mainly enriched in the IM vs. IE groups and IL vs. IM groups, and the 10 pathways with the most significant enrichment were Staphylococcus aureus infection, methane metabolism, carbon metabolism, complement and coagulation cascades, glyoxylate and dicarboxylate metabolism, carbon fixation in photosynthetic organisms, leishmaniasis, phagosome, pathogenic Escherichia coli infection, and biosynthesis of amino acids (Additional file 16–19). In addition, four overlapping pathways were enriched in all four groups, including the PPAR signaling pathway, peroxisome, legionellosis, and adipocytokine signaling pathway (Fig. 4D).

Four metal metabolites-related pathways were detected, including the HIF-1 signaling pathway, PI3K-Akt signaling pathway, TNF signaling pathway, and MAPK signaling pathway. Copper- and zinc-related pathways included two-component systems, Alzheimer’s disease, and mineral absorption. Among them, the number of enriched calcium-related pathways was the largest. Calcium-related pathways contained the insulin signaling pathway, mTOR signaling pathway, HIF-1 signaling pathway, and p53 signaling pathway. The significantly enriched DEPs and related pathways were screened (Table 2).

Table 2 DEPs in metal metabolite-related pathways and those with the most significant enrichment in the ricefield eel Monopterus albus

Identification of sex-related DEPs during sex change

The most significant KEGG pathways and GO terms in each comparison group based on p < 0.05 were screened. KEGG pathways were mainly enriched in the pentose phosphate pathway, Staphylococcus aureus infection, methane metabolism, carbon metabolism, and the TCA cycle. The main GO terms included cytoplasm, ion binding, and RNA binding. These functional annotations have been highly correlated with energy metabolism and ion binding. Then, several sex-related DEPs in these significant pathways and among these terms were identified, including matrix metalloproteinase-9 (Mmp-9), alcohol dehydrogenase 1-like isoform X1 (Adh1-l), Vatpg1, Tpi, and Igf2bp3 (Table 3).

Table 3 The most significantly enriched terms and pathway sex-related DEPs in the ricefield eel Monopterus albus

Characterization of metal metabolite -related DEPs during sex change

As mentioned above, a few metal metabolites-related signaling pathways related with sex changes were determined in the ricefield eel. To explore the expression patterns of these factors during sex change, a total of 111 nonredundant DEPs were identified from 15 iron-related proteins, 14 copper-related proteins, 10 zinc-related proteins, and 131 calcium-related proteins. To analyze the expression patterns of these 111 nonredundant DEPs more comprehensively, the transcriptome and proteome data from the five stages were analyzed based on the log2-fold change at both the protein level and mRNA level. Only four of these genes [LOC109952670, LOC109953418, LOC109953912, and calcium Binding Protein 39 Like (cab39l)] had consistent expression using RT‒qPCR. Two genes (LOC109953912 and cab39l) matched the transcriptome data (Fig. 5 and Additional file 1: Fig. S2p < 0.05).

Fig. 5
figure 5

Gene expression levels of metal metabolite-related DEPs from the RNA-seq and RT‒qPCR data during sex change in the ricefield eel Monopterus albus. A, the expression level of cab39l. B, the expression level of LOC109953912. qRT-PCR, quantitative real-time PCR; RNA-seq, RNA sequencing; OV, ovary; IE, early intersex stage gonad; IM, middle intersex stage gonad; IL, late intersex stage gonad; TE, testis. The results are presented as the means ± SEMs. Means marked with different letters were significantly different (p < 0.05)

Metal ion content during sex change

To investigate whether the levels of metal metabolites varied during sex change in the ricefield eel, the levels of Fe2+/Fe3+, Cu2+, Ca2+ and Zn2+ were measured during the five periods. The results showed that there were no significant changes in the levels of Fe2+/Fe3+ or Ca2+, while there were significant changes in the Cu2+ and Zn2+ levels in gonads during sex change (Fig. 6A). Importantly, the expression trends of Fe2+/Fe3+ and Ca2+ were opposite. In addition, the expression trend of Cu2+ increased with the process of sex change and reached the highest level in TE, which suggests that Cu2+ may be a male-biased metal metabolite. In contrast, the expression trend of Zn2+ increased and then decreased, and the content was the highest during the IE and TE periods. In addition, the levels of serum metal metabolites during sex change altered little for Fe2+/Fe3+ and Ca2+, while Cu2+ and Zn2+ remained significantly changed (Fig. 6B). However, Cu2+ was no longer elevated with the course of sex change, and changes in its content occurred mainly during intersex.

Fig. 6
figure 6

Trends in gonads and serum levels of metal metabolites during sex change in the ricefield eel Monopterus albus. A, levels of Fe2+/Fe3+, Cu2+, Ca2+ and Zn2+ in gonads. B, levels of Fe2+/Fe3+, Cu2+, Ca2+ and Zn2+ in serum. OV, ovary; IE, early intersex stage gonad; IM, middle intersex stage gonad; IL, late intersex stage gonad; TE, testis. Different lowercase letters indicate significant between-group differences (p < 0.05), while the same lowercase letters indicate that the between-group differences were not significant (p < 0. 05)

Discussion

Unlike the transcriptome, the proteomics-based approach for detection the expressed proteins in tissues or cell types has complemented the genome initiatives [24, 25]. To date, there has been a gradual increase in proteomic studies in fish gonads, which have revealed changes in the protein composition of the gonads, contributing to shedding light on the gonadal status. In the gill of tropical marbled eel, salinity acclimation proteins were analyzed by iTRAQ proteomic method [26]. The significant correlation between gonadal index (GI) and protein content might play key roles in changing the condition factor of gonads in sea urchin Strongylocentrotus nudus [27]. The proteome of carp sperm between the immobilized and activated spermatozoa mainly involved in ubiquitin‒proteasome pathways, glycolysis, the TCA cycle, and remodeling related to sperm energy metabolism and motility [28]. Furthermore, the proteomic studies during zebrafish oocyte maturation, 190 proteins were significantly changed between immature and mature oocytes, and several proteins such as vitellogenin (Vtg3) and 14-3-3 protein involved in oocyte development [29]. Therefore, proteomics is an important tool to study fish gonad development. Two-dimensional electrophoresis was used to isolate the DEP spots in ricefield eel [22]. Over 80 DEPs were identified in ricefield eel, in which there were five spots highly expressed in the ovary, three spots in the ovotestis, 16 in the testis, 24 in both the ovary and ovotestis, and 17 in both the testis and ovotestis [22]. Chromosome group 3 (Cbx3) was highly expressed in the testis, and a member of the RAS oncogene family (Rab37) was highly expressed in the ovary of ricefield eel [22]. However, two-dimensional electrophoretic proteome analysis has limitations, in which the protein cannot be completely acquired and fully utilized in a database, and the intersex stage of research is not systematic enough for ricefield eel.

In the present study, proteins from five periods of sex change in the ricefield eel (OV, IE, IM, IL, and TE) were comprehensively and systematically investigated. There were the most DEPs from the IM period to the TE period. Especially in the IL vs. IM group, clustered heatmap male-biased DEPs were abundant. Compared to previously reported studies, our study obtained 5362 total proteins and 335 DEPs (two-dimensional electrophoresis obtained only approximately 80 proteins [22]) and performed functional enrichment analysis. In addition, we have more details on the intersex stages (IE, IM, and IL). In our previous transcriptome study, we found that apoptosis-related genes were strongly related with sex change, and we screened for tgfb3, fshr, camk4, and calmodulin, which might be key genes [30]. Although some results on the process of sex change at the transcriptional level have been obtained, the relevant studies at the protein level need to be further improved. Therefore, high-throughput proteome sequencing was utilized to characterize the expression levels of large proteins to understand the molecular mechanisms of sex change in ricefield eel.

With GO and KEGG analyses, some proteins with significant sex differences were screened, including Vatpg1, Tpi, and Igf2bp3. In mouse epididymis and vas deferens, vacuolar H+ ATPase was necessary for establishing an acidic lumen pH environment that kept sperm at rest during sperm maturation and storage in these organs [31, 32]. Tpi was involved in the gene expression of retinol during sperm metabolism in rats [33]. Igf2bp3 was related to embryonic development and placenta formation in medaka [34]. In our present work, Vatpg1, Tpi and Igf2bp3 were the DEPs that were markedly increased in the OV to IM period. While the gonad lamella was gradually occupied by testicular tissues from the OV to IM period, the number of ovarian tissues were gradually reduced [35]. But Vatpg1, Tpi and Igf2bp3 were significantly reduced from IM to TE stage, and then returned to level of expression similar to the OV period. The gonadal lamella was gradually filled with testicular tissues and gradual reduction of ovarian tissues until it disappeared from the IM to TE period [35]. Above results suggest that Vatpg1, Tpi and Igf2bp3 might play important roles in oocyte disappearance and spermatogenesis during sex change of ricefield eel.

Previous studies have shown that metal metabolites, such as cadmium (Cd) [36], iron (Fe) [37, 38], copper (Cu) [39], lead (Pb) [40], and zinc (Zn) [41], play key roles in gonadal development, germ cell growth and maturation, and reproduction in mammals and fish. In our study, functional enrichment analysis of the DEPs focused on metal metabolite-related pathways and terms, such as Cu2+-, Ca2+-, Zn2+- and Fe2+/Fe3+-related pathways and terms. Among them, the abundance of calcium-related pathways and terms was the highest. In different vertebrate lineages, the calcium and redox pathways regulate sex determination, acting as the crucial missing link between sex and the environment [42]. Moreover, two metal metabolite-related DEPs were identified (LOC109953912 and cab39l) that were consistent at both transcriptomic and proteomic expression levels. These two DEPs (LOC109953912 and cab39l) were enriched in the Parkinson’s disease, ubiquitin-mediated proteolysis, gap junction, and phagosome pathogenic Escherichia coli infection pathways. These results suggest that metal metabolites may be involved in several biological functions during sex change in ricefield eel.

Chronic exposure to these metal metabolites impaired mammalian follicle production and fish ovarian development by interfering with reproductive hormones and reactive oxygen species (ROS) production, which in turn damaged various molecules, including proteins, lipids, and DNA, and disrupted antioxidant defenses, ion homeostasis and endoplasmic reticulum homeostasis [40, 43, 44]. Four metal metabolites (Cu2+, Ca2+, Zn2+ and Fe2+/Fe3+) were detected in our study. Among them, Ca2+ and Fe2+ did not vary greatly, but there was some correlation and antagonism between their levels because Ca2+ affects the absorption of Fe2+/Fe3+ [45]. In contrast, the content of Cu2+ tended to increase gradually during sex change in our study, suggesting that Cu2+ might be an important metal metabolite in spermatheca development. A previous study found that Cu might cause a decrease in antioxidant capacity and testicular spermatogenesis [38]. In addition, Cu exposure upregulated the expression of oxidative phosphorylation pathway genes (Cytochrome c oxidase, Sod1, and Gst) in the testis to induce ROS production, which was related with Cu toxicity in the testis [38, 46, 47]. Long-term exposure to high levels of Cu had a toxic effect on yellow catfish (Pelteobagrus fulvidraco), whereas at low doses and after relatively short-term exposures, serum steroid hormones and steroidogenesis-related gene levels were elevated [48]. These results reflect the disruption of the balance between antioxidants and oxidants after Cu exposure, affecting ovary and testis development. In contrast, the content of Zn2+ showed a trend of increasing, then decreasing, and then increasing again, and was the highest in the IE period, indicating that zinc ions have an important role in the process from OV to IM. Zn is an essential nutrient involved in many physiological processes, such as follicular development, immune response, homeostasis, oxidative stress, cell cycle progression, DNA replication, DNA damage repair, apoptosis, and aging [41]. Zn deficiency could lead to blocking the oocyte meiotic process, cumulus expansion, and follicle ovulation [41]. In flounder (Pseudopleuronectes americanus), females stored Zn in summer for winter gonadal development [49]. Thus, metal metabolites play important functions in regulating gonadal development by influencing the balance between antioxidants and oxidants, as well as in steroid hormone production.

(Cu-Zn) Sod1 and Mmp-9 were identified as potential proteins related with crosstalk between metal metabolites and sex-related proteins in ricefield eel. (Cu-Zn) Sod1 is an abundant copper- and zinc-containing protein [50]. Its primary function is to act as an antioxidant enzyme, lowering the steady-state concentration of superoxide, but when mutated, it can also cause disease [50]. In our data, (Cu-Zn) Sod1 was mainly enriched in the peroxisome pathway. Peroxisomes are highly dynamic cell organelles that play key roles in cellular lipid and hydrogen peroxide (H2O2) metabolism [51]. In the present work, the abundances of (Cu-Zn) Sod1 in the intersexual stages (especially during the IM stage) were significantly higher than those in the OV and TE stages, suggesting that (Cu-Zn) Sod1 had been associated with the sex change process in the ricefield eel. The expression level of Mmp-9 was progressively lower with sex change, suggesting that it is likely to be a male-related protein. Mmp-9 has been immunolocalized during the development of the follicular membrane and interstitium in rodents [52]. The level of Mmp-9 in sheep follicles was elevated from normal follicles to atresia follicles [53]. High expression levels of Mmp-9 were also accompanied by an increase in the rate of follicular apoptosis in patients [54]. These data suggest that (Cu-Zn) Sod1 and Mmp-9 are likely to influence gonadal development during sex change by affecting the level of ROS. In the future, we may perform more detailed functional studies on (Cu-Zn) Sod1 and Mmp-9 in the sex change process of ricefield eel.

Conclusion

In summary, we analyzed the protein profiles of ricefield eel gonadal tissues in five stages (OV, IE, IM, IL, and TE) and verified the plausibility of the data using western blotting. Then, we analyzed the functional enrichment of metal metabolite-related DEPs. Sex-related DEPs (Vatpg1, Tpi and Igf2bp3) may play important roles in oocyte disappearance and spermatogenesis during sex change. In addition, we detected the contents of the metal ions Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+ at these five stages and screened for (Cu-Zn) Sod1 and Mmp-9 as possible key proteins in the sex reversal process. Our data therefore nicely complement the protein data from ricefield eel. As mentioned above, through proteome analysis of the sex change stage of the ricefield eel, detailed genetic information about the mechanism of sex change has been more deeply described, including ovarian maturation and degeneration, gonadal differentiation, and sex changes.

Methods

Experimental animals

Wild ricefield eels (n = 210, body length = 39.42 ± 5.13 cm, and body weight = 48.53 ± 24.67 g) were purchased from a local market (Chengdu, Sichuan, China). The fish were maintained in the laboratory at a water temperature of 21.7 ± 2.5 °C under a photoperiod of 16 h light:8 h dark. Fish were decapitated after anesthesia with 0.02% tricaine buffer (80 µg/L) (Sigma, Saint Louis, MO, USA). Then, the gonads were collected and immediately stored in liquid nitrogen at − 80 °C. The gonads were divided into two parts. The first was fixed with Bouin’s solution for 24 h and then stored in 75% ethanol for determination of the gonadal developmental stage, while the second was immediately stored in liquid nitrogen at − 80 °C until RNA and protein extraction. Histological classification of the gonads, including OV, IE, IM, IL, and TE, has been described previously [16]. Three biological replicates were randomly selected from gonadal tissues at the OV, IE, IM, IL, and TE stages for proteome sequencing (Fig. 7).

Fig. 7
figure 7

Identification of the gonad developmental stages in the ricefield eel Monopterus albus. A, ovary (OV). B, early intersex gonad (IE). C, middle intersexual gonad (IM). D, late intersex gonad (IL). E, testis (TE). CAO, cortical alveolar oocyte; PGO, primary growth stage oocyte; EVO, early vitellogenesis oocyte; GR, gonadal ridge; SC, spermatocyte; ST, spermatid

Total protein extraction, isobaric tags for relative and absolute quantification (iTRAQ) peptide labeling, and liquid chromatography–tandem mass spectrometry (LC‒MS/MS) analysis

Each sample was ground individually in liquid nitrogen and lysed with PASP lysis buffer (100 mM NH4HCO3, 8 M urea, pH 8), followed by 5 min of ultrasonication on ice. The lysate was centrifuged at 12,000 × g for 15 min at 4 °C, and the supernatant was reduced with 10 mM DL-dithiothreitol for 1 h at 56 °C and subsequently alkylated with sufficient iodoacetamide for 1 h at room temperature in the dark. Then, each sample was completely mixed with 4 volumes of precooled acetone by vortexing and incubated at -20 °C for at least 2 h. The samples were then centrifuged at 12,000 × g for 15 min at 4 °C, and the precipitate was collected. After washing with 1 mL of cold acetone, the pellet was dissolved in dissolution buffer.

The volume of each protein sample was made up to 100 µL with DB dissolution buffer (8 M urea, 100 mM triethylammonium formate, pH8.5). Trypsin and 100 mM triethylammonium formate buffer were added, and the sample was digested at 37 °C for 4 h. Then, trypsin and CaCl2 were added, and the sample was digested overnight. Formic acid was mixed with the digested sample, which was adjusted to pH 3 and centrifuged at 12,000 × g for 5 min at room temperature. The supernatant was slowly loaded onto a C18 desalting column, which was washed with washing buffer (0.1% formic acid, 3% acetonitrile) 3 times, and then the sample was eluted with elution buffer (0.1% formic acid, 70% acetonitrile). The eluents of each sample were collected and lyophilized. Twenty microliters of 1 M TEAB buffer was added for reconstitution, and enough iTRAQ labeling reagent (dissolved in isopropanol) was added. The sample was mixed with shaking for 2 h at room temperature. Then, the reaction was stopped by adding 100 µL of 50 mM Tris-HCl (pH = 8). Equal volumes of each labeled sample were combined, mixed, desalted, and lyophilized. For multiple labeled groups, a common reference was created by pooling an equal quantity of each sample.

Samples were injected into homemade C18 Nano-well columns. The peptides were isolated using the linear gradients listed for the homemade analytical columns. The isolated peptides were analyzed by Q Exactive™ series mass spectrometers (Thermo Fisher) with a Nanosury Flex™ ion source (ESI), a spray voltage of 2.1 kV, and an ion transport capillary temperature of 320 °C. The following parameters were used: full scan range from m/z 407 to 1500, a resolution of 60,000 (at m/z 200), an automatic gain control target of 3 × 106, and a maximum ion implantation time of 20 ms. The top 40 precursors with the highest abundance in the full scan were selected for fragmentation by high-energy collision dissociation and analyzed by MS/MS at a resolution of 15,000 (at m/z 200), an automatic gain control target value of 5 × 104, a maximum ion implantation time of 45 ms, a normalized collision energy of 32%, an intensity threshold of 2.2 × 104, and a dynamic rejection parameter of 20 s.

Analysis of the DEPs

The relative quantitative values of the unique peptides were obtained based on the peak area from the original spectrum and peptide spectrum matches (PSMs) by Proteome Discoverer 2.2 (PD 2.2, Thermo). The relative quantitative value of each unique peptide could be corrected based on the quantitative information from all unique peptides in each protein. The protein quantitation results were statistically analyzed by T test, and the corresponding p value was calculated as a significance indicator. When FC ≥ 1.5 and p ≤ 0.05, the protein was upregulated, and when FC ≤ 0.67 and p ≤ 0.05, the protein was downregulated.

Gene Ontology (GO) functional analysis was conducted by the program InterProScan against nonredundant protein databases (including Pfam, PRINTS, ProDom, SMART, ProSite, and PANTHER) [55]. Then, the Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were utilized to identify the protein families and pathways. Volcano maps, cluster heatmaps, and GO, COG and KEGG enrichment analyses of the DEPs are presented [56].

Joint transcriptome-proteome analysis

Correlation analysis of the transcriptome and proteome data at the fold change of the two omics was carried out. Some metal-related genes/proteins from the transcriptome and proteome data were screened. Then, the expression patterns of these coidentified genes/proteins based on the log2-fold change values were analyzed. A log2-fold change > 0 indicated upregulation, and a log2-fold change < 0 indicated downregulation.

Western blot verification

Total protein was extracted from the ricefield eel gonads by the Tissue or Cell Total Protein Extraction Kit (Sangon Biotech C510003). All components were loaded with equal amounts of protein extract on a 12% polyacrylamide gel (standard gel, Bio-Rad), and the isolated proteins were transferred to a polyvinylidene fluoride membrane. The primary antibodies included those against the following proteins: Igf2bp3 polyclonal antibody (1:8000, Proteintech, 14642-1-AP, 86% homology), Tpi rabbit polyclonal antibody (1:4000, Proteintech, 10713-1-AP, 89% homology), Atp6v1g1 rabbit polyclonal antibody (1:2000, Proteintech, 16143-1-AP, 84% homology), and Sod1 rabbit pAb (1:1000, ABclonal, A0274, 70% homology). According to the molecular weight of proteins, specific bands for Igf2bp3 (64 kDa), Tpi (27 kDa), Vatpg1 (14 kDa) and Sod1 (16 kDa) were first detected in the gonad tissue of ricefield eel (Additional file 1: Fig S3). Therefore, membranes are cropped on the position of these proteins at scale of 10–20 kDa above and below and then the antibody incubations were performed. The protein concentration was determined by BCA (Beyotime, P0012S), and the band values were quantified by ImageJ software.

Quantitative real-time PCR (qRT‒PCR) verification

Total RNA was isolated from gonads in different developmental stages using TRIzol reagent (Invitrogen, Chicago, IL, USA) and then reverse-transcribed into cDNA using a RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, Waltham, MA, USA) according to the manufacturer’s protocol. qRT‒PCR was performed using a CFX Connect system (Bio-Rad, Chicago, IL, USA) in a final reaction volume of 10 µL comprising 5 µL of 2× SYBR Green Master Mix (TaKaRa, Dalian, China), 0.4 µL of each primer (10 µmol/L), 3.2 µL of nuclease-free water, and 1 µL of cDNA template. The cycling parameters were 95 °C for 5 min followed by 40 amplification cycles of 95 °C for 10 s, 59 °C for 15 s, and 72 °C for 20 s. The specificity of PCR amplification was confirmed by melting curve analysis, agarose gel electrophoresis, and sequencing of the PCR products. To minimize variation due to differences in cDNA loading, the expression levels of the target genes were normalized to the geometric mean expression levels of elongation factor-1 alpha (ef1α) [57]. 2-ΔΔCT was performed to identify the mRNA gene expression levels according to the cycle threshold values [58]. Differential expression levels were compared by the IBM SPSS Statistics 20 system at a significance level of less than 0.05 [59].

The mRNA levels of LOC109952670, LOC109953418, LOC109953912, and cab39l were normalized to the geometric mean expression levels of ef1α. The sequences of all the qRT‒PCR primers used in this study are provided in Table 4.

Table 4 Primers for real-time quantitative PCR analysis

Metal and enzyme activity assays

The contents of Zn2+ (E011-1-1), Cu2+ (E010-1-1), Ca2+ (C004-2-1), and Fe2+/Fe3+ (A039-2-1) were determined by their respective detection kits from Nanjing Jiancheng Bioengineering Institute. All procedures were performed according to the instructions.

Statistical analysis

The results are shown as the mean ± standard error (SE) and were analyzed by GraphPad Prism 8.0 software. Technical and biological triplicates were performed for each experiment. The differences in experimental data among groups were analyzed by one-way analysis of variance (ANOVA). A value of p < 0.05 was considered to indicate statistical significance.

Data availability

The datasets generated and analysed during the current study are available in the iProX database with accession number IPX0007324000 (https://www.iprox.cn//page/SCV017.html?query=IPX0007324000). The datasets analysed during this study are included in this published article and its supplementary information files. Please contact Zhi He (zhihe@sicau.edu.cn) if someone wants to request the data from this study.

References

  1. Pennell MW, Mank JE, Peichel CL. Transitions in sex determination and sex chromosomes across vertebrate species. Mol Ecol. 2018;27(19):3950–63.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Korpelainen H. Sex ratios and conditions required for environmental sex determination in animals. Biol Rev Camb Philos Soc. 1990;65(2):147–84.

    Article  CAS  PubMed  Google Scholar 

  3. Li XY, Mei J, Ge CT, Liu XL, Gui JF. Sex determination mechanisms and sex control approaches in aquaculture animals. Sci China Life Sci. 2022;65(6):1091–122.

    Article  PubMed  Google Scholar 

  4. Li XY, Gui JF. Diverse and variable sex determination mechanisms in vertebrates. Sci China Life Sci. 2018;61(12):1503–14.

    Article  PubMed  Google Scholar 

  5. Cutting A, Chue J, Smith CA. Just how conserved is vertebrate sex determination? Dev Dyn. 2013;242(4):380–7.

    Article  CAS  PubMed  Google Scholar 

  6. Ortega-Recalde O, Goikoetxea A, Hore TA, Todd EV, Gemmell NJ. The genetics and epigenetics of sex change in fish. Annu Rev Anim Biosci. 2020;8:47–69.

    Article  PubMed  Google Scholar 

  7. Todd EV, Ortega-Recalde O, Liu H, Lamm MS, Rutherford KM, Cross H, Black MA, Kardailsky O, Marshall Graves JA, Hore TA, et al. Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change. Sci Adv. 2019;5(7):eaaw7006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhao X, Luo M, Li Z, Zhong P, Cheng Y, Lai F, Wang X, Min J, Bai M, Yang Y et al. Chromosome-scale assembly of the Monopterus genome. Gigascience 2018, 7(5).

  9. Gallaway BJ, Szedlmayer ST, Gazey WJ. A life history review for red snapper in the gulf of mexico with an evaluation of the importance of offshore petroleum platforms and other artificial reefs. Rev Fish Sci. 2009;17(1):48–67.

    Article  Google Scholar 

  10. Thomas JT, Todd EV, Muncaster S, Lokman PM, Damsteegt EL, Liu H, Soyano K, Gléonnec F, Lamm MS, Godwin JR, et al. Conservation and diversity in expression of candidate genes regulating socially-induced female-male sex change in wrasses. PeerJ. 2019;7:e7032.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Casas L, Saborido-Rey F, Ryu T, Michell C, Ravasi T, Irigoien X. Sex change in clownfish: molecular insights from transcriptome analysis. Sci Rep. 2016;6:35461.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wu GC, Li HW, Huang CH, Lin HJ, Lin CJ, Chang CF. The testis is a primary factor that contributes to epigenetic modifications in the ovaries of the protandrous black porgy, Acanthopagrus schlegelii. Biol Reprod. 2016;94(6):132.

    Article  PubMed  Google Scholar 

  13. Yvonne S, De Mitcheson Min L. Functional hermaphroditism in teleosts. Fish Fisheries. 2008;9(1):1–43.

  14. Mary K, Hart Andrew W., Kratter Philip H, Crowley. Partner fidelity and reciprocal investments in the mating system of a simultaneous hermaphrodite. Behav Ecol. 2016;27(5):1471–79.

  15. Cheng H, He Y, Zhou R. Swamp eel (Monopterus albus). Trends Genet. 2021;37(12):1137–8.

    Article  CAS  PubMed  Google Scholar 

  16. He Z, Deng F, Ma Z, Zhang Q, He J, Ye L, Chen H, Yang D, He L, Luo J, et al. Molecular characterization, expression, and apoptosis regulation of siva1 in protogynous hermaphrodite fish ricefield eel (Monopterus albus). Fish Physiol Biochem. 2021;47(5):1585–96.

    Article  CAS  PubMed  Google Scholar 

  17. He Z, Wu Y, Xie J, Wang T, Zhang L, Zhang W. Growth differentiation factor 9 (Gdf9) was localized in the female as well as male germ cells in a protogynous hermaphroditic teleost fish, ricefield eel Monopterus albus. Gen Comp Endocrinol 2012, 178(2):355–362.

  18. He Z, Li Y, Wu Y, Shi S, Sun C, Deng Q, Xie J, Wang T, Zhang W, Zhang L. Differentiation and morphogenesis of the ovary and expression of gonadal development-related genes in the protogynous hermaphroditic ricefield eel Monopterus albus. J Fish Biol. 2014;85(5):1381–94.

    Article  CAS  PubMed  Google Scholar 

  19. He Z, Deng F, Yang D, He Z, Hu J, Ma Z, Zhang Q, He J, Ye L, Chen H, et al. Crosstalk between sex-related genes and apoptosis signaling reveals molecular insights into sex change in a protogynous hermaphroditic teleost fish, ricefield eel Monopterus albus. Aquaculture. 2022;552:737918.

    Article  CAS  Google Scholar 

  20. He Z, Ma Z, Yang D, Chen Q, He Z, Hu J, Deng F, Zhang Q, He J, Ye L, et al. Circular RNA expression profiles and CircSnd1-miR-135b/c-foxl2 axis analysis in gonadal differentiation of protogynous hermaphroditic ricefield eel Monopterus albus. BMC Genomics. 2022;23(1):552.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. He Z, Ye L, Yang D, Ma Z, Deng F, He Z, Hu J, Chen H, Zheng L, Pu Y, et al. Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus). BMC Genomics. 2022;23(1):450.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Sheng Y, Zhao W, Song Y, Li Z, Luo M, Lei Q, Cheng H, Zhou R. Proteomic analysis of three gonad types of swamp eel reveals genes differentially expressed during sex reversal. Sci Rep. 2015;5:10176.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Liu H, Todd EV, Lokman PM, Lamm MS, Godwin JR, Gemmell NJ. Sexual plasticity: a fishy tale. Mol Reprod Dev. 2017;84(2):171–94.

    Article  CAS  PubMed  Google Scholar 

  24. Beynon RJ. The dynamics of the proteome: strategies for measuring protein turnover on a proteome-wide scale. Brief Funct Genomic Proteomic. 2005;3(4):382–90.

    Article  CAS  PubMed  Google Scholar 

  25. Banks RE, Dunn MJ, Hochstrasser DF, Sanchez JC, Blackstock W, Pappin DJ, Selby PJ. Proteomics: new perspectives, new biomedical opportunities. Lancet. 2000;356(9243):1749–56.

    Article  CAS  PubMed  Google Scholar 

  26. Jia Y, Yin S, Li L, Li P, Liang F, Wang X, Wang X, Wang L, Su X. iTRAQ proteomic analysis of salinity acclimation proteins in the gill of tropical marbled eel (Anguilla marmorata). Fish Physiol Biochem. 2016;42(3):935–46.

    Article  CAS  PubMed  Google Scholar 

  27. Shang WH, Han JR, Yan JN, Du YN, Xu YS, Xue CF, Zhang TT, Wu HT, Zhu BW. Quantitative Proteome reveals variation in the Condition factor of Sea Urchin Strongylocentrotus nudus during the fishing season using an iTRAQ-based Approach. Mar Drugs. 2019;17(7):397.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Dietrich MA, Dietrich GJ, Mostek A, Ciereszko A. Motility of carp spermatozoa is associated with profound changes in the sperm proteome. J Proteom. 2016;138:124–35.

    Article  CAS  Google Scholar 

  29. Ge C, Lu W, Chen A. Quantitative proteomic reveals the dynamic of protein profile during final oocyte maturation in zebrafish. Biochem Biophys Res Commun. 2017;490(3):657–63.

    Article  CAS  PubMed  Google Scholar 

  30. Chi W, Gao Y, Hu Q, Guo W, Li D. Genome-wide analysis of brain and gonad transcripts reveals changes of key sex reversal-related genes expression and signaling pathways in three stages of Monopterus albus. PLoS ONE. 2017;12(3):e0173974.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Au CL, Wong PY. Luminal acidification by the perfused rat cauda epididymidis. J Physiol. 1980;309:419–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Pietrement C, Sun-Wada GH, Silva ND, McKee M, Marshansky V, Brown D, Futai M, Breton S. Distinct expression patterns of different subunit isoforms of the V-ATPase in the rat epididymis. Biol Reprod. 2006;74(1):185–94.

    Article  CAS  PubMed  Google Scholar 

  33. Russell DL, Kim KH. Expression of triosephosphate isomerase transcripts in rat testis: evidence for retinol regulation and a novel germ cell transcript. Biol Reprod. 1996;55(1):11–8.

    Article  CAS  PubMed  Google Scholar 

  34. Wu X, Zhang Y, Xu S, Chang Y, Ye Y, Guo A, Kang Y, Guo H, Xu H, Chen L, et al. Loss of Gsdf leads to a dysregulation of Igf2bp3-mediated oocyte development in medaka. Gen Comp Endocrinol. 2019;277:122–9.

    Article  CAS  PubMed  Google Scholar 

  35. Zhang Y, Zhang W, Yang H, Zhou W, Hu C, Zhang L. Two cytochrome P450 aromatase genes in the hermaphrodite ricefield eel Monopterus albus: mRNA expression during ovarian development and sex change. J Endocrinol. 2008;199(2):317–31.

    Article  CAS  PubMed  Google Scholar 

  36. Jackson TW, Baars O, Belcher SM. Gestational cd exposure in the CD-1 mouse sex-specifically disrupts essential metal ion homeostasis. Toxicol Sci. 2022;187(2):254–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Shao Y, Zhao H, Wang Y, Liu J, Zong H, Xing M. Copper-mediated mitochondrial fission/fusion is associated with intrinsic apoptosis and autophagy in the testis tissues of chicken. Biol Trace Elem Res. 2019;188(2):468–77.

    Article  CAS  PubMed  Google Scholar 

  38. Li Y, Chen H, Liao J, Chen K, Javed MT, Qiao N, Zeng Q, Liu B, Yi J, Tang Z, et al. Long-term copper exposure promotes apoptosis and autophagy by inducing oxidative stress in pig testis. Environ Sci Pollut Res Int. 2021;28(39):55140–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Cao J, Wang G, Wang T, Chen J, Wenjing G, Wu P, He X, Xie L. Copper caused reproductive endocrine disruption in zebrafish (Danio rerio). Aquat Toxicol. 2019;211:124–36.

    Article  CAS  PubMed  Google Scholar 

  40. Qu J, Niu H, Wang J, Wang Q, Li Y. Potential mechanism of lead poisoning to the growth and development of ovarian follicle. Toxicology. 2021;457:152810.

    Article  CAS  PubMed  Google Scholar 

  41. Liao X, Wu L, Yin D, Tian D, Zhou C, Liu J, Li S, Zhou J, Nie Y, Liao H, et al. The role of zinc in follicular development. Mol Biol Rep. 2023;50(5):4527–34.

    Article  CAS  PubMed  Google Scholar 

  42. Castelli MA, Whiteley SL, Georges A, Holleley CE. Cellular calcium and redox regulation: the mediator of vertebrate environmental sex determination? Biol Rev Camb Philos Soc. 2020;95(3):680–95.

    Article  PubMed  Google Scholar 

  43. Abdel-Khalek AA. Chronic exposure to water of lake qaroun induced metal-related testicular damage and endocrine disruption in male fish. Biol Trace Elem Res. 2018;185(1):197–204.

    Article  CAS  PubMed  Google Scholar 

  44. Qu J, Wang Q, Sun X, Li Y. The environment and female reproduction: potential mechanism of cadmium poisoning to the growth and development of ovarian follicle. Ecotoxicol Environ Saf. 2022;244:114029.

    Article  CAS  PubMed  Google Scholar 

  45. Lönnerdal B. Calcium and iron absorption–mechanisms and public health relevance. Int J Vitam Nutr Res. 2010;80(4–5):293–9.

    Article  PubMed  Google Scholar 

  46. Zhao D, Zhang X, Li X, Ru S, Wang Y, Yin J, Liu D. Oxidative damage induced by copper in testis of the red swamp crayfish Procambarus clarkii and its underlying mechanisms. Aquat Toxicol. 2019;207:120–31.

    Article  CAS  PubMed  Google Scholar 

  47. Roychoudhury S, Nath S, Massanyi P, Stawarz R, Kacaniova M, Kolesarova A. Copper-induced changes in reproductive functions: in vivo and in vitro effects. Physiol Res. 2016;65(1):11–22.

    Article  CAS  PubMed  Google Scholar 

  48. Zhang LH, Luo Z, Song YF, Shi X, Pan YX, Fan YF, Xu YH. Effects and mechanisms of waterborne copper exposure influencing ovary development and related hormones secretion in yellow catfish Pelteobagrus fulvidraco. Aquat Toxicol. 2016;178:88–98.

    Article  CAS  PubMed  Google Scholar 

  49. Fletcher GL, King MJ. Seasonal dynamics of Cu2+, Zn2+, Ca2+, and Mg2+ in gonads and liver of winter flounder (Pseudopleuronectes americanus): evidence for summer storage of Zn2+ for winter gonad development in females. Can J Zool. 1978;56(2):284–90.

    Article  CAS  Google Scholar 

  50. Valentine JS, Doucette PA, Zittin Potter S. Copper-zinc superoxide dismutase and amyotrophic lateral sclerosis. Annu Rev Biochem. 2005;74:563–93.

    Article  CAS  PubMed  Google Scholar 

  51. Aquaculture ReportsFransen M, Lismont C. Redox signaling from and to peroxisomes: progress, challenges, and prospects. Antioxid Redox Signal. 2019;30(1):95–112.

    Article  Google Scholar 

  52. García R, Ballesteros LM, Hernández-Pérez O, Rosales AM, Espinosa R, Soto H, Díaz de León L, Rosado A. Metalloproteinase activity during growth, maturation and atresia in the ovarian follicles of the goat. Anim Reprod Sci. 1997;47(3):211–28.

    Article  PubMed  Google Scholar 

  53. Huet C, Monget P, Pisselet C, Monniaux D. Changes in extracellular matrix components and steroidogenic enzymes during growth and atresia of antral ovarian follicles in the sheep. Biol Reprod. 1997;56(4):1025–34.

    Article  CAS  PubMed  Google Scholar 

  54. Ben-Shlomo I, Goldman S, Shalev E. Regulation of matrix metalloproteinase-9 (MMP-9), tissue inhibitor of MMP, and progesterone secretion in luteinized granulosa cells from normally ovulating women with polycystic ovary disease. Fertil Steril. 2003;79(Suppl 1):694–701.

    Article  PubMed  Google Scholar 

  55. Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.

    Article  PubMed  Google Scholar 

  57. Hu Q, Guo W, Gao Y, Tang R, Li D. Reference gene selection for real-time RT-PCR normalization in rice field eel (Monopterus albus) during gonad development. Fish Physiol Biochem. 2014;40(6):1721–30.

    Article  CAS  PubMed  Google Scholar 

  58. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  59. Namdar P, YekeFallah L, Jalalian F, Barikani A, Razaghpoor A. Improving door-to-balloon time for patients with Acute ST-Elevation myocardial infarction: a controlled clinical trial. Curr Probl Cardiol. 2021;46(3):100674.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the National Natural Science Foundation of China for supporting this experiment.

Funding

This research was supported by the National Natural Science Foundation of China [grant numbers 31972777, and 31402286].

Author information

Authors and Affiliations

Authors

Contributions

ZH: Conceptualization, formal analysis, supervision, project administration. FX: Conceptualization, methodology, investigation, statistical analysis, writing draft. DYY: Project administration, validation and editing. FQD: Conceptualization, methodology, supervision, writing, original and draft. WXD: Investigation and statistical analysis. ZDH: Resources and investigation. SQW: Investigation. QQC: Methodology and investigation. MQC: Investigation. HCW: Investigation. KG: Resources. JXX: Resources. ZTT: Investigation and resources. MWZ: Conceptualization and investigation. TMY: Conceptualization and investigation.

Corresponding author

Correspondence to Taiming Yan.

Ethics declarations

Ethics approval and consent to participate

All experimental protocols involved in fishes in this study were conducted in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Sichuan Agricultural University and ARRIVE guidelines (https://arriveguidelines.org). All procedures and investigations were reviewed and approved by the Animal Research and Ethics Committees of Sichuan Agricultural University and performed in accordance with the guidelines of the committee (Approval No.20190031).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, Z., Xiao, F., Yang, D. et al. Protein expression patterns and metal metabolites in a protogynous hermaphrodite fish, the ricefield eel (Monopterus albus). BMC Genomics 25, 500 (2024). https://doi.org/10.1186/s12864-024-10397-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10397-w

Keywords