Skip to main content

Role of cytoskeleton-related proteins in the acrosome reaction of Eriocheir sinensis spermatozoa

Abstract

Cytoskeleton-related proteins are essential for cell shape maintenance and cytoskeleton remodeling. The spermatozoa of Eriocheir sinensis (Chinese mitten crab) have a unique cellular structure, and the mechanism of spermatozoal metamorphosis during the acrosome reaction is not well understood. In this study, the E. sinensis spermatozoa were induced using calcium ionophore A23187 to undergo the acrosome reaction in vitro, and the acrosome-reacting and fresh (non-reacting) spermatozoa were collected separately. The differential expression of cytoskeleton-related protein genes in acrosome-reacting and fresh spermatozoa of E. sinensis was analyzed by whole transcriptome sequencing and bioinformatics analysis, and PPI network and miRNA-mRNA regulation network were constructed to analyze their possible function and regulation mechanism. The results showed that numerous differentially expressed cytoskeleton-related protein genes, miRNAs and lncRNAs were found in acrosome-reacting and fresh spermatozoa of E. sinensis; 27 cytoskeleton-related protein genes were down regulated and 687 miRNAs were up regulated in acrosome-reacting spermatozoa; 147 miRNAs target these 27 cytoskeleton-related protein genes. In the PPI networks, RAC1, BCAR1, RDX, NCKAP1, EPS8, CDC42BPA, LIMK1, ELMO2, GNAI1 and OCRL were identified as hub proteins. These proteins are mainly involved in the regulation of cytoskeleton organization, actin cytoskeleton organization, microtubule skeleton organization and small GTPase-mediated signal transduction and other biological processes, and play roles in pathways such as actin cytoskeletal regulation and axon guidance. miR-9, miR-31 and two novel miRNAs in the miRNA-mRNA regulatory network are the core miRNAs targeting cytoskeleton-related protein genes. miR-9 targets and regulates OBSCN, CDC42BPA, ELMO2, BCAS3, TPR and OCRL; while miR-31 targets and regulates CDC42BPA and TPR. This study provides a theoretical basis for revealing the mechanism of acrosome reaction under the special spermatozoa morphology of E. sinensis.

Peer Review reports

Acrosome reaction is a necessary event prior to sperm-egg fusion and is essential for fertilization to occur [1]. The spermatozoa of Eriocheir sinensis (Chinese mitten crab) are distinct from those of other mammals, with a unique cellular structure that is characterized by the absence of flagella, an extra-large acrosome capsule and a "cup-shaped" spermatozoa nucleus [2, 3]. Therefore, the mechanism of the acrosome reaction has attracted much attention from researchers. Cytoskeleton consists of actin filaments, microtubules and intermediate filaments, which play a key role in cell morphology maintenance and cell motility [4]. Cytoskeleton-related proteins, as a class of proteins that regulate the dynamic stability of the cytoskeleton, include tubulin, actin, myosin, and many others structural unit accessory proteins. Currently, research on cytoskeleton-related proteins has focused on humans, mice, guinea pigs, bovines and other mammals [5,6,7,8,9]. Cytoskeleton-related proteins are involved in regulating cytoskeleton remodeling, determining and controlling cell morphology, and also in the transport of intracytoplasmic components [10,11,12]. In mammals, cytoskeletal remodeling is thought to be necessary for mammalian sperm capacitation and acrosome reaction [13, 14]. Similar to mammalian spermatozoa, E. sinensis spermatozoa also undergo dramatic morphological changes and cytoskeleton remodeling during the acrosome reaction. The acrosome reaction of E. sinensis spermatozoa include radial arm contraction, acrosome vesicle outgrowth, acrosome tubule protrusion and lamellar structure shedding [15]. In this series of changes, cytoskeleton-related proteins may play important roles. However, since the role of cytoskeleton-related proteins in the acrosome reaction of E. sinensis spermatozoa is less reported, it is not clear which cytoskeleton-related proteins are involved in the acrosome reaction of E. sinensis spermatozoa, nor do we know the processes involved and the manner in which they are regulated. Therefore, the exploration of expression changes in cytoskeleton-related proteins and their regulatory mechanisms during the acrosome reaction of E. sinensis spermatozoa can help to reveal the mechanism of spermatozoal metamorphosis during the acrosome reaction in the special spermatozoa morphology of E. sinensis.

Accordingly, we induced acrosome reaction of E. sinensis spermatozoa in vitro and collected acrosome-reacting spermatozoa and fresh (non-reacting) spermatozoa separately. Then, Whole transcriptome sequencing and bioinformatics analysis were used to analyze the differential expression of cytoskeleton-related protein genes in acrosome-reacting spermatozoa and fresh spermatozoa. We also constructed protein–protein interaction (PPI) network and miRNA-mRNA regulatory network to analyze their possible functional roles and regulatory mechanisms. This study can provide valuable information for the reproduction research of E. sinensis and provide a certain theoretical basis for animal developmental biology research in the future.

Materials and methods

Experimental animals

Chinese mitten crabs were purchased from Ice King Aquatic Products Sales Co., Ltd. in Suqian City, Jiangsu Province, China ([SYXK (Su) 2018–0021]). A total of 30 male crabs with a body weight of 100-150 g, healthy limbs and active movement were selected and randomly divided into 2 groups. There were 18 in the experimental group (group JA) and 12 in the control group (group JO). Then they were fully aerated under tap water for 1 day. The water temperature was 24℃, pH was 6.8, and dissolved oxygen > 5.0 mg/L. The isolation and extraction of vas deferens tissues were performed at the Experimental Animal Center of Youjiang Medical University for Nationalities ([SYXK(GUI)2017–0004]).

Main reagents and instruments

Total RNA extraction kit (Beijing Solarbio Science & Technology Co., Ltd.), Ultramicro ultraviolet spectrophotometer (Shanghai Puyuan Instrument Co., Ltd.), single-double ultra-clean workbench (Jiangsu Sujing Group Co., Ltd.), biological microscope (Shanghai Leica Instruments Co., Ltd.), desktop high-speed refrigerated centrifuge (Eppendorf China Co., Ltd.), 15 ml glass tissue grinder (Guangxi Yanpu Biotechnology Co., Ltd.).

Acrosome-reacting spermatozoa and fresh spermatozoa preparation

A small amount (about 100 mg) of vas deferens tissue was taken from each Chinese mitten crab and quickly washed several times with calcium-free artificial seawater buffer of pH 8.1. The cleaned vas deferens tissues were mixed into one sample of every 6 animals in the experimental group and one sample of every 4 animals in the control group. Then the spermatozoa pods were gently squeezed out with forceps to prepare a spermatozoa pod suspension, and the supernatant was discarded after standing for 10 min; the spermatozoa pods precipitate was added to a 15 ml clean centrifuge tube and quickly resuspended with calcium-free artificial seawater buffer; then transferred to a glass tissue homogenizer to gently grind the spermatozoa pods about 40 ~ 50 times in order to induce the release of spermatozoa in the pods. Then, spermatozoa pod tissue fragments were removed after centrifugation at 500 r/min for 10 min, and spermatozoa were collected after centrifugation at 4000 r/min for 10 min. In the spermatozoa precipitate of the control group, 10 times the volume of RNA preservation solution was added and mixed, and permeated through 2 h of room temperature to prepare the fresh spermatozoa. In the spermatozoa precipitate of the experimental group, an appropriate amount of calcium ionophore A23187 was added, and then frozen overnight in a -80 °C refrigerator, and the acrosome reaction was observed the next day. After the acrosome reaction was successful, centrifuged at 4000 r/min for 10 min to obtain the precipitate of the acrosome-reacting spermatozoa. Then, 10 times the volume of RNA preservation solution was added to the precipitate of the acrosome-reacting spermatozoa, mixed well and infiltrated at room temperature for 2 h to serve as the acrosome-reacting spermatozoa for subsequent experiments.

RNA extraction, sequencing and data analysis

Take the prepared suspension of acrosome-reacting spermatozoa and fresh spermatozoa, three samples in each group, centrifuge at 4000 r/min for 10 min to discard the RNA preservation solution. Total RNA was extracted using a total RNA extraction kit and sent to Shenzhen Huada Gene Technology Co., Ltd. for whole transcriptome sequencing. Sequencing quality assurance and depth was performed as usual. The preliminary analysis of the raw sequencing data was done by Shenzhen Huada Gene Technology Co., Ltd.

Analysis of differentially expressed genes

Differentially expressed genes (DEGs), miRNAs and lncRNAs in acrosome-reacting spermatozoa and fresh spermatozoa from the sequencing data were analyzed. Differential folds were selected as Log2FC(JO/JA) ≥ 1 or ≤ -1, and P < 0.05 as the screening criteria for differential genes. The official gene names (Official Symbol) were annotated by comparison with the Chinese mitten crab reference genome, and the DEGs that were down regulated in acrosome-reacting spermatozoa were screened.

Screening of cytoskeleton-related proteins

The Gene Ontology and GO Annotations database (https://www.ebi.ac.uk/QuickGO/) was searched for "sperm cytoskeleton proteins". Cytoskeleton-related GO terms of mouse, nematode, drosophilid fly and other model organisms as well as human were selected, and the cytoskeleton-related protein genes in each GO term were downloaded and extracted. All genes were summarized after removing duplicate genes. Further compared with the differential genes that are down-regulated expressed in the acrosome-reacting spermatozoa from our sequencing data, we finally screened the cytoskeleton-related protein genes in acrosome-reacting spermatozoa and fresh spermatozoa of E. sinensis.

Functional enrichment of cytoskeleton-related protein genes and PPI network construction

The GO Annotations of candidate genes were obtained by mapping candidate genes to background gene sets using the R package org.HS.eg.db (version 3.1.0), and enrichment analysis was performed using the R package ClusterProfiler (version 3.14.3) to obtain enrichment analysis results, setting the minimum gene set number to 5, the maximum gene set number to 5000, P < 0.05, and FDR < 0.01. The protein–protein interaction (PPI) network was constructed using the STRING database (https://cn.string-db.org/), setting the minimum confidence limit equal to 0.15 and selecting the default for other parameters.

miRNA screening and miRNA-mRNA regulatory network construction for targeting cytoskeleton-related protein genes

TargetScan (http://www.targetscan.org/vert_71/), miRanda (http://www.microrna.org/microrna/home.do/) and RNAhybrid (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/) were used to predict and analyze the miRNAs regulating the above cytoskeleton-related protein genes. The miRNA-mRNA regulatory network was constructed and visualized using Cytoscape software, and the hub genes in the network were analyzed using the cytoHubba plugin, by assigning and ordering each gene with the topological network algorithm.

Statistical analysis

The data of this study were analyzed using SPSS 23.0 statistical software. Continuous normal distribution of measurement data was statistically described by mean ± standard deviation. Independent samples t-test was used for analysis, and differences were considered statistically significant at P < 0.05.

Results

Acrosome-reacting spermatozoa and fresh spermatozoa

The spermatozoa of E. sinensis were induced to undergo acrosome reaction in vitro, and the results are shown below (Fig. 1). When acrosome reaction occurs in spermatozoa, the radial arm contracts, the apical cap breaks to form an aperture, the acrosome vesicle reverses, the acrosome tubule extends forward, and then the lamellar structure shrinks and disappears, finally forms the structures as shown in the figure (Fig. 1B), such as a Y-shaped fork; but fresh (non-reacting) spermatozoa are round (Fig. 1A). The diagram for the process of acrosome reaction occurring in the spermatozoa of E. sinensis was shown in Fig. 1C.

Fig. 1
figure 1

Acrosome-reacting spermatozoa and fresh spermatozoa in Eriocheir sinensis. A fresh spermatozoa; B acrosome-reacting spermatozoa; C Diagram for the process of acrosome reaction in spermatozoa. The process of acrosome reaction can be divided into four stages: (1) the contraction of radial arm and the breaking to form an aperture on apical cap; (2) the reversion of acrosome vesicle; (3) the extension of acrosome tubule; (4) the shrinking and disappearance of lamellar structures. RA: Radial arm; AC: Apical cap; AT: Acrosome tubule; LS: Lamellar structures; ML: Middle layer; FL: Fribrous layer; N: Nucleus

Analysis of differentially expressed genes between acrosome-reacting spermatozoa and fresh spermatozoa

Through whole transcriptome sequencing of acrosome-reacting spermatozoa and fresh spermatozoa, a total of 6510 up-regulated genes, 8028 down-regulated genes and 9942 non-differentially expressed genes were identified in the two samples. The Volcano Plot of DEGs shows these genes (Fig. 2). Among them, 849 mRNAs were down regulated in acrosome-reacting spermatozoa and up regulated in fresh spermatozoa; on the contrary, 687 miRNAs were up regulated in acrosome-reacting spermatozoa and down regulated in fresh spermatozoa.

Fig. 2
figure 2

Volcano Plot of differentially expressed genes (DEGs) in the sample of acrosome-reacting spermatozoa and fresh spermatozoa of Eriocheir sinensis. Up-regulated genes are in red, down-regulated genes are in blue, non-differentially expressed genes are in gray

Screening of cytoskeleton-related protein genes

After the search of Gene Ontology and GO Annotation database, a total of 10 GO terms directly related to the spermatozoa cytoskeleton were retrieved, and their annotated gene numbers are shown in Table 1. 1976 genes were obtained after removing duplicate genes. Further comparison with 849 genes that were down regulated in acrosome-reacting spermatozoa from the sequencing data, finally yielded 27 cytoskeleton-related protein genes (Fig. 3). They were EPS8, OCRL, DCP2, SSH1, CHMP5, CUL3, LIMK1, TAOK1, CORO6, NCKAP1, CYLD, CALD1, BCR, GNAI1, TPR, FKBP4, BCAS3, RAC1, RAE1, KAT2A, ELMO2, ALDOA BCAR1, CDC42BPA, OBSCN, RDX and CALR.

Table 1 Number of genes annotated by the 10 GO terms related to cytoskeleton in Gene Ontology and GO Annotation database
Fig. 3
figure 3

Intersection of cytoskeleton-related protein genes in the database and differentially expressed genes in whole transcriptome sequencing data. Gene set A comes from Gene Ontology and GO Annotation database, with a total of 1976 genes; Gene set B comes from whole transcriptome sequencing of acrosome-reacting spermatozoa and fresh spermatozoa of Eriocheir sinensis, with a total of 849 genes

Functional enrichment and PPI network of cytoskeleton-related proteins

The PPI network was constructed for 27 cytoskeleton-related proteins by STRING database, and the results showed that they had strong interactions (Fig. 4). Among them, 10 proteins, including RAC1, BCAR1, RDX, NCKAP1, EPS8, CDC42BPA, LIMK1, ELMO2, GNAI1 and OCRL, were the hub proteins in the network. Functional enrichment analysis of the 27 cytoskeleton-related protein genes using the R package showed that they mainly have molecular functions such as Ras GTPase binding, small GTPase binding, GTPase binding, Rho GTPase binding, tubulin binding, actin binding, protein serine/threonine kinase activity and cell adhesion molecule binding, etc.; mainly in the cell leading edge, lamellipodium, actin cytoskeleton, focal adhesion, spindle, microtubules, site of polarized growth, A band and M band of spermatozoa, etc.; mainly involved in the regulation of cytoskeleton organization, microtubule cytoskeleton organization, regulation of actin cytoskeleton organization, small GTPase mediated signal transduction, Ras protein signal transduction, regulation of plasma membrane bounded cell projection assembly, regulation of cell projection assembly, cell protrusion organization, regulation of mitotic spindle organization and other biological processes; the pathways involved are mainly enriched in actin cytoskeleton regulation and axon guidance (Fig. 5).

Fig. 4
figure 4

The protein–protein interaction (PPI) network of 27 cytoskeleton-related proteins

Fig. 5
figure 5

GO and KEGG enrichment analysis of 27 cytoskeleton-related protein genes. A all KEGG pathways in which the 27 cytoskeleton-related protein genes are involved [16]; B Top20 biological processes in which the 27 cytoskeleton-related protein genes are involved; C Top20 cellular components in which the 27 cytoskeleton-related protein genes are located; D all molecular functions of the 27 cytoskeleton-related protein genes

miRNAs targeting cytoskeleton-related protein genes and miRNA-mRNA regulatory network

Based on the three databases of TargetScan, miRanda and RNAhybrid, the miRNAs targeting cytoskeleton-related protein genes were comprehensively predicted. It was finally found that 147 miRNAs and the 27 cytoskeleton-related protein genes may have targeted regulatory relationships. They formed a total of 276 targeted regulatory relationship pairs, which formed a strong miRNA-mRNA regulatory network (Fig. 6). Among them, 4 miRNAs (miR-9, miR-31, novel_mir10 and novel_mir3) and 6 genes (OBSCN, CDC42BPA, ELMO2, BCAS3, TPR and OCRL) constituted the core of the network, and they interacted with each other (Fig. 7). The core miRNAs in the network, miR-9, miR-31, novel_mir10 and novel_mir3, all have target-directed actions on multiple cytoskeleton-related protein genes. miR-9 has a possible target regulation relationship with OBSCN, CDC42BPA, ELMO2, BCAS3, TPR and OCRL; miR-31 has a possible target regulation relationship with CDC42BPA and TPR. In addition, the miRNAs targeting the top10 core genes in the PPI network in this miRNA-mRNA regulatory network are shown in Table 2.

Fig. 6
figure 6

miRNA-mRNA regulatory network consisting of 27 cytoskeleton-related protein genes and 147 miRNAs. mRNA is indicated by red diamonds, miRNA is indicated by blue ovals; the miRNAs beginning with “novel” have not yet been named in Eriocheir sinensis

Fig. 7
figure 7

Interactions between 4 core miRNAs and 6 core cytoskeleton-related protein genes. mRNA is indicated by diamonds, miRNA is indicated by ovals; the miRNAs beginning with “novel” are unnamed miRNAs in Eriocheir sinensis

Table 2 The miRNAs targeting top10 genes in PPI network

Discussion

The unique cytoskeleton can not only regulate cell shape, structure, stability, plasticity and dynamic response to environmental stimuli, but also maintain the biological function of certain organs [17]. The spermatozoa of E. sinensis have a cytoskeleton structure that is significantly different from that of normal mammalian spermatozoa. The spermatozoa are nearly spherical, tailless, with uneven surface, numerous radiating arms, and a huge spherical acrosome, which compresses the nucleus into a cup shape [18, 19]. This special cytoskeletal structure determines that the process of the acrosome reaction is also different: contraction of the radial arms, outgrowth of the acrosome capsule, forward extension of the acrosome canal, and finally shedding of the lamellar structure [15]. Each of these steps is inseparable from cytoskeleton remodeling, and its spermatozoa cytoskeleton-related proteins should play a significant role. Using the cytoskeleton-related proteins of model organisms such as mouse, nematode, drosophilid fly and other model organisms as well as human as the background, we screened out cytoskeleton-related protein genes from the acrosome-reacting spermatozoa and fresh spermatozoa of E. sinensis. And we further analyzed their differential expression and function in the acrosome-reacting spermatozoa and fresh spermatozoa of E. sinensis, and the miRNAs targeting them. The study will help to better understand their possible regulatory mechanisms in the acrosome reaction of E. sinensis spermatozoa.

Through whole transcriptome sequencing, we found a very large number of differentially expressed genes between acrosome-reacting spermatozoa and fresh spermatozoa of E. sinensis, and many cytoskeleton-related protein genes were up regulated in fresh spermatozoa and down regulated in the acrosome reaction of E. sinensis spermatozoa. The reasons for the down-regulated expression of cytoskeleton-related protein genes in E. sinensis spermatozoa are poorly understood. microRNAs (miRNAs) are an important class of endogenous factors with the ability to regulate the expression of protein-coding genes. In recent years, it has been found that its role in testicular development and spermatogenesis has become more and more extensive [20, 21]. They involved in various biological processes such as male germ cell differentiation, gonadal oxidative stress, and testicular immune homeostasis through targeted inhibition of gene expression [22, 23], and can be considered as crucial regulators of biological processes [24]. We found numerous miRNAs differentially expressed between acrosome-reacting spermatozoa and fresh spermatozoa, many of which were up regulated in acrosome-reacting spermatozoa of E. sinensis. Through the prediction of miRNAs and the construction of miRNA-mRNA regulatory network, we found that 147 miRNAs and 27 cytoskeleton-related protein genes have possible target regulatory relationships. They formed 276 pairs of targeted regulatory relationships. Among the 276 targeted regulatory relationship pairs formed, miR-9 and OBSCN, CDC42BPA, ELMO2, BCAS3, TPR, and OCRL, and miR-31 and CDC42BPA and TPR constituted the core of the targeted regulatory relationship network. miR-9 has been reported to promote the proliferation, invasion and metastasis of various cancer cells [25,26,27] and also has a role in regulating spermatogenesis of spermatogonia [25]. It is evident that miR-9 has a significant pro-differentiation effect on cells, but its involvement in acrosome reaction has not been reported. miR-31 has been reported to regulate apoptosis and invasion of hepatocellular carcinoma HepG2 cells through the ROCK1/F-actin pathway [28]. It can be seen that miR-31 can regulate cytoskeleton-related protein genes involved in the physiological functions of cells, but its study in acrosome reaction has also not been mentioned yet. Acrosome reaction of spermatozoa is a special physiological process, especially the acrosome reaction of E. sinensis spermatozoa has a different reaction process from that of general mammals. In the present study, the target gene of miR-9, Obscurin, cytoskeletal calmodulin and titin-interacting RhoGEF (OBSCN), a cytoskeletal calmodulin, can mediate direct interactions between the sarcoplasmic reticulum and sarcoplasmic fibers by participating in the organization of sarcoplasmic fibers during cytoskeletal assembly [29]; CDC42 binding protein kinase alpha (CDC42BPA), which can mediate Cdc42-MRCK signaling to induce peripheral actin formation and promote cytoskeleton reorganization [30]; BCAS3 microtubule associated cell migration factor (BCAS3), a microtubule associated cell migration factor, has been reported to promote cytoskeletal reorganization by coupling microtubules and intermediate filaments [31]. However, whether the specific roles of these protein genes in the unique acrosome reaction of E. sinensis spermatozoa are different? We don't know yet. Therefore, how miR-9 and miR-31 targeting these cytoskeleton-related protein genes are involved in the acrosome reaction of E. sinensis spermatozoa deserve subsequent in-depth study.

We know that the biological functions of miRNAs mainly act through their target genes. miRNAs target genes, and these target genes perform specific functions due to their own special domains, thereby participating in biological processes. In our present study, we found that 27 cytoskeleton-related protein genes were down regulated in acrosome reaction of E. sinensis spermatozoa, and the proteins encoded by RAC1, BCAR1, RDX, NCKAP1, EPS8, CDC42BPA, LIMK1, ELMO2, GNAI1 and OCRL are hub proteins in the PPI network of 27 cytoskeleton-related proteins. Their functional enrichment analysis revealed that they mainly have molecular functions such as Ras GTPase binding, small GTPase binding, GTPase binding, Rho GTPase binding, tubulin binding, actin binding, protein serine/threonine kinase activity and cell adhesion molecule binding, etc. Ras-GTPase, an activating protein, is considered an important putative of Ras inactivation and plays an important role in tissue homeostasis [32]. Local GTPase activation has been reported to lead to the exposure of effector structural domains, which upon binding to effector proteins can initiate downstream signaling [33]. G protein subunit alpha i1 (GNAI1) which has GTPase enzymatic activity can bind to GTPase activating protein, and play an important role in signal transduction [34]. Small GTPases control the signal transduction of transforming growth factor beta (TGFβ), which plays an important role in receptor endocytosis and actin cytoskeleton remodeling [35]. RAS-related C3 botulinum toxin substrate 1 (RAC1) has a small GTPase-binding protein structural domain that binds to small GTPases and is involved in actin cytoskeleton remodeling in the apical region of the acrosome during guinea pig spermatozoa capacitation. The formation of the actin cytoskeleton in the apical region of the acrosome is considered a necessary event for spermatozoa capacitation and acrosome reaction in guinea pigs and is driven by Rac1 [36]. In this study, we also found that RAC1 is differentially expressed in the acrosome-reacting spermatozoa of E. sinensis, suggesting that RAC1 may also regulate actin cytoskeleton remodeling in E. sinensis spermatozoa through its own specific structural domain and participate in the process of acrosome reaction. In addition, the extension and contraction of membrane protrusions and adhesion are controlled by the Rho family of small GTPases [37]. Rho GTPase, as the structural basis of signaling, has been reported to be directly related to miRNAs and regulated by miRNAs and covalent modifications [38]. OCRL inositol polyphosphate-5-phosphatase (OCRL) has a Rho GTPase activating protein domain and can interact with Rho GTPase in the trans-Golgi network, by activating phosphatidylinositol 4, 5-bisphosphate 5-phosphatase, thereby affecting actin polymerization [39]. When the acrosome reaction of E. sinensis spermatozoa occurs, the acrosome capsule will be everted and the outer membrane of the acrosome will be stretched and contracted. Therefore, proteins with such domains may play an important role in acrosome reaction.

In summary, we induced acrosome reaction of E. sinensis spermatozoa in vitro and found a large number of differentially expressed mRNAs and miRNAs in acrosome-reacting spermatozoa and fresh spermatozoa. By comparing with model organisms such as mouse, nematode, drosophilid fly and human cytoskeleton-related proteins, we finally found 27 down-regulated cytoskeleton-related protein genes in the acrosome-reacting spermatozoa of E. sinensis. They are mainly involved in the regulation of cytoskeleton organization, actin cytoskeleton organization, microtubule skeleton organization and small GTPase mediated signal transduction, acting through actin cytoskeleton regulation and axon guidance pathways. RAC1, BCAR1, RDX, NCKAP1, EPS8, CDC42BPA, LIMK1, ELMO2, GNAI1 and OCRL are the hub proteins in the PPI network of 27 cytoskeleton-related proteins. In the strong miRNA-mRNA regulatory network formed, 4 miRNAs with 6 cytoskeleton-related protein genes formed the core of the regulatory network. miR-9 had possible target regulatory relationships with OBSCN, CDC42BPA, ELMO2, BCAS3, TPR and OCRL, and miR-31 with CDC42BPA and TPR. The spermatozoa of E. sinensis have a unique cytoskeletal structure and the mechanism of acrosome reaction is obviously different. In this study, we investigated the expression and functional roles of cytoskeleton-related protein genes in acrosome-reacting spermatozoa and predicted the miRNAs that may regulate them, which may provide a theoretical basis for further understanding the mechanism of the acrosome reaction in E. sinensis spermatozoa. However, further validation of the miRNAs targeting the expression of related genes is needed to better elucidate the mechanism of acrosome response in E. sinensis and provide more valuable information for its reproduction and developmental biology of other animals. However, the targeting of miRNAs to cytoskeleton-related protein genes still needs to be further validated in order to provide more valuable information for their reproduction and developmental biology of other animals.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the Sequence Read Archive (SRA), https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA807466.

References

  1. Siu KK, Serrao V, Ziyyat A, et al. The cell biology of fertilization: gamete attachment and fusion[J]. J Cell Biol. 2021;220(10):e202102146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Chen T, Mu S, Sun Z, et al. Spermiogenic histone transitions and chromatin decondensation in Decapoda[J]. Theriogenology. 2020;156:242–52.

    Article  CAS  PubMed  Google Scholar 

  3. Wei YL, Yang T, Kovacs T, et al. C-terminal kinesin motor es-KIFC1 regulates nuclear formation during spermiogenesis in Chinese mitten crab Eriocheir sinensis[J]. Gene. 2019;719:144074.

    Article  CAS  PubMed  Google Scholar 

  4. Svitkina T. Imaging cytoskeleton components by electron Microscopy[J]. Methods Mol Biol. 2022;2364:25–52.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Brener E, Rubinstein S, Cohen G, et al. Remodeling of the actin cytoskeleton during mammalian sperm capacitation and acrosome reaction[J]. Biol Reprod. 2003;68(3):837–45.

    Article  CAS  PubMed  Google Scholar 

  6. Roa-Espitia AL, Hernandez-Rendon ER, Baltierrez-Hoyos R, et al. Focal adhesion kinase is required for actin polymerization and remodeling of the cytoskeleton during sperm capacitation[J]. Biol Open. 2016;5(9):1189–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Romarowski A, Luque GM, La Spina FA, et al. Role of actin cytoskeleton during mammalian sperm acrosomal exocytosis[J]. Adv Anat Embryol Cell Biol. 2016;220:129–44.

    Article  PubMed  Google Scholar 

  8. Cabello-Agueros JF, Hernandez-Gonzalez EO, Mujica A. The role of F-actin cytoskeleton-associated gelsolin in the guinea pig capacitation and acrosome reaction[J]. Cell Motil Cytoskeleton. 2003;56(2):94–108.

    Article  CAS  PubMed  Google Scholar 

  9. Howes EA, Hurst SM, Jones R. Actin and actin-binding proteins in bovine spermatozoa: potential role in membrane remodeling and intracellular signaling during epididymal maturation and the acrosome reaction[J]. J Androl. 2001;22(1):62–72.

    CAS  PubMed  Google Scholar 

  10. Benoit B, Baillet A, Pous C. Cytoskeleton and Associated Proteins: Pleiotropic JNK Substrates and Regulators[J]. Int J Mol Sci. 2021;22(16):8375.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kamal A, Goldstein LS. Connecting vesicle transport to the cytoskeleton[J]. Curr Opin Cell Biol. 2000;12(4):503–8.

    Article  CAS  PubMed  Google Scholar 

  12. Gautel M. The sarcomeric cytoskeleton: who picks up the strain?[J]. Curr Opin Cell Biol. 2011;23(1):39–46.

    Article  CAS  PubMed  Google Scholar 

  13. Delgado-Buenrostro NL, Mujica A, Chiquete-Felix N, et al. Role of Wasp and the small GTPases RhoA, RhoB, and Cdc42 during capacitation and acrosome reaction in spermatozoa of English guinea pigs[J]. Mol Reprod Dev. 2016;83(10):927–37.

    Article  CAS  PubMed  Google Scholar 

  14. Breitbart H, Finkelstein M. Actin cytoskeleton and sperm function[J]. Biochem Biophys Res Commun. 2018;506(2):372–7.

    Article  CAS  PubMed  Google Scholar 

  15. Sun WJ, Zhu M, Wang YL, et al. ERK is involved in the process of acrosome reaction in vitro of the Chinese mitten crab, Eriocheir sinensis[J]. Mar Biotechnol (NY). 2015;17(3):305–16.

    Article  CAS  PubMed  Google Scholar 

  16. Kanehisa M, Furumichi M, Sato Y, et al. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;27:gkac963.

    Google Scholar 

  17. Schell C, Huber TB. The evolving complexity of the podocyte cytoskeleton[J]. J Am Soc Nephrol. 2017;28(11):3166–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Zhu M, Sun WJ, Wang YL, et al. P38 participates in spermatogenesis and acrosome reaction prior to fertilization in Chinese mitten crab Eriocheir sinensis[J]. Gene. 2015;559(2):103–11.

    Article  CAS  PubMed  Google Scholar 

  19. Zeng J, Peng S, Zhong S, et al. The Spermatozoal Ultrastructure of the Chinese Mitten Crab (Eriocheir sinensis)[J]. J Microsc Ultrastruct. 2019;7(4):181–4.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Wang YY, Duan SH, Wang GL, et al. Integrated mRNA and miRNA expression profile analysis of female and male gonads in Hyriopsis cumingii[J]. Sci Rep. 2021;11(1):665.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Walker WH. Regulation of mammalian spermatogenesis by miRNAs[J]. Semin Cell Dev Biol. 2021;15:S1084-9521(21)00118-X.

    Google Scholar 

  22. Alves M, Arruda RP, Batissaco L, et al. Changes in miRNA levels of sperm and small extracellular vesicles of seminal plasma are associated with transient scrotal heat stress in bulls[J]. Theriogenology. 2021;161:26–40.

    Article  CAS  PubMed  Google Scholar 

  23. Li T, Wang X, Luo R, et al. Integrating miRNA and mRNA Profiling to Assess the Potential miRNA-mRNA Modules Linked With Testicular Immune Homeostasis in Sheep[J]. Front Vet Sci. 2021;8:647153.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Rahbar S, Pashaiasl M, Ezzati M, et al. MicroRNA-based regulatory circuit involved in sperm infertility[J]. Andrologia. 2020;52(1):e13453.

    Article  PubMed  Google Scholar 

  25. Gao Y, Jia D, Hu Q, et al. Foxl3, a Target of miR-9, Stimulates Spermatogenesis in Spermatogonia During Natural Sex Change in Monopterus albus[J]. Endocrinology. 2016;157(11):4388–99.

    Article  PubMed  Google Scholar 

  26. Shu Z, Gao F, Xia Q, et al. MiR-9-5p promotes cell proliferation and migration of hepatocellular carcinoma by targeting CPEB3[J]. Biomark Med. 2021;15(2):97–108.

    Article  CAS  PubMed  Google Scholar 

  27. Zheng J, Shao S, Dai C, et al. miR-9-5p promotes the invasion and migration of endometrial stromal cells in endometriosis patients through the SIRT1/NF-kappaB pathway[J]. Int J Clin Exp Pathol. 2020;13(7):1859–66.

    PubMed  PubMed Central  Google Scholar 

  28. Zhang X, Xu L, Yang T. miR-31 Modulates Liver Cancer HepG2 Cell Apoptosis and Invasion via ROCK1/F-Actin Pathways[J]. Onco Targets Ther. 2020;13:877–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kontrogianni-Konstantopoulos A, Jones EM, Van Rossum DB, et al. Obscurin is a ligand for small ankyrin 1 in skeletal muscle[J]. Mol Biol Cell. 2003;14(3):1138–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wilkinson S, Paterson HF, Marshall CJ. Cdc42-MRCK and Rho-ROCK signalling cooperate in myosin phosphorylation and cell invasion[J]. Nat Cell Biol. 2005;7(3):255–61.

    Article  CAS  PubMed  Google Scholar 

  31. Joshi D, Inamdar MS. Rudhira/BCAS3 couples microtubules and intermediate filaments to promote cell migration for angiogenic remodeling[J]. Mol Biol Cell. 2019;30(12):1437–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. King PD, Lubeck BA, Lapinski PE. Nonredundant functions for Ras GTPase-activating proteins in tissue homeostasis[J]. Sci Signal. 2013;6(264):e1.

    Article  Google Scholar 

  33. Ten KJ, Hordijk PL. Targeting and localized signalling by small GTPases[J]. Biol Cell. 2007;99(1):1–12.

    Article  Google Scholar 

  34. Li ZW, Sun B, Gong T, et al. GNAI1 and GNAI3 Reduce Colitis-Associated Tumorigenesis in Mice by Blocking IL6 Signaling and Down-regulating Expression of GNAI2[J]. Gastroenterology. 2019;156(8):2297–312.

    Article  CAS  PubMed  Google Scholar 

  35. Kardassis D, Murphy C, Fotsis T, et al. Control of transforming growth factor beta signal transduction by small GTPases[J]. FEBS J. 2009;276(11):2947–65.

    Article  CAS  PubMed  Google Scholar 

  36. Ramirez-Ramirez D, Salgado-Lucio ML, Roa-Espitia AL, et al. Rac1 is necessary for capacitation and acrosome reaction in guinea pig spermatozoa[J]. J Cell Biochem. 2020;121(4):2864–76.

    Article  CAS  PubMed  Google Scholar 

  37. Nguyen LK, Kholodenko BN, von Kriegsheim A. Rac1 and RhoA: Networks, loops and bistability[J]. Small GTPases. 2018;9(4):316–21.

    Article  CAS  PubMed  Google Scholar 

  38. Liu M, Bi F, Zhou X, et al. Rho GTPase regulation by miRNAs and covalent modifications[J]. Trends Cell Biol. 2012;22(7):365–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Suchy SF, Nussbaum RL. The deficiency of PIP2 5-phosphatase in Lowe syndrome affects actin polymerization[J]. Am J Hum Genet. 2002;71(6):1420–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Thanks to the support of the National Natural Science Foundation of China (Grant No. 31960728) and the Guangxi Graduate Education Innovation Program Teaching Reform Project (Grant No. JGY2020164).

Funding

This research is supported by the National Natural Science Foundation of China (Grant Nos.31960728) and the Guangxi Graduate Education Innovation Program Teaching Reform Project (Grant Nos. JGY2020164).

Author information

Authors and Affiliations

Authors

Contributions

Yulian Tang, Lishuang Sun, and Shu Li prepared the samples and analyzed the experimental data. Huiting Liu, Lvjing Luo and Zhengyu Chen prepared all the experimental reagents and materials. Yulian Tang completed the manuscript, Genliang Li conceived and designed the study and provided guidance of experiment and paper writting. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Genliang Li.

Ethics declarations

Ethics approval and consent to participate

The animal we studied was Eriocheir sinensis. Eriocheir sinensis is an important economic crab in China, which can be bought in the Chinese market. Therefore, this study has been granted an exemption from requiring ethics approval. All the animals we purchased were adapted under standard laboratory conditions (ventilation room, water temperature 25 ± 1 °C, pH 6.8, dissolved oxygen ≥ 5.0 mg/L), and were free to obtain standard water and food. All procedures were carried out in the Experimental Animal Center of Youjiang Medical University for Nationalities (SYXK(GUI) 2017–0004) in accordance with the "Guidelines for Animal Care and Use" (China).

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.

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

Tang, Y., Sun, L., Li, S. et al. Role of cytoskeleton-related proteins in the acrosome reaction of Eriocheir sinensis spermatozoa. BMC Genom Data 24, 4 (2023). https://doi.org/10.1186/s12863-023-01112-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-023-01112-x

Keywords