Nicola Bernabò*, Ilaria Saponaro, Mauro Mattioli and Barbara Barboni
Department of Comparative Biomedical Sciences, University of Teramo, Teramo, Italy
Received Date: October 18, 2011; Accepted Date: November 15, 2011; Published Date: November 18, 2011
Citation: Bernabò N, Saponaro I, Mattioli M, Barboni B (2012) Signaling Strategy in Spermatozoa Activation of Sea Urchin, C. elegans and Human: Three Different Players for the Same Melody. J Bioeng Biomed Sci S5:006. doi:10.4172/2155-9538.S5-006
Copyright: © 2012 Bernabò N, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Bioengineering & Biomedical Science
Human; C. elegans; Sea urchin; Spermatozoa; Sperm activation; Bioinformatics; Biological networks; Computational model
Many cellular systems undergo a process of activation during their metabolic activity. For instance the contraction of striated muscle, the release of synaptosomes from the presynaptic endplate, the degranulation (exocytosis) of sensitized mast cells, the activation of mammals eggs, the activation of visual transmission in mammalian rods, the activation of male gametes to reach the full fertilizing ability require an activation process. This last event, in particular, has attracted the attention of Researches since its discovery in 50’s because of its important implications both for basic (developmental biology, endocrinology, biochemistry) and applied science (andrology, male infertility, contraception). Thus sperm activation has been studied in several organisms, of different Phyla, and the data about this important topic are growing day by day. The availability of bioinformatics tools has provided to the Scientists new resources to investigate this phenomenon in its complexity. Overcoming the reductionist approach until now used, the use of bioinformatics resources makes possible the aggregation of molecular data in explicative models. The use of computational modeling to explore the biochemical architecture of mammalian sperm activation, thus allowing the quantitative analysis of this complex signaling system, has been recently introduced  and validated  by our group. In this context it could be very interesting to modelize and compare the sperm activation process in organisms showing completely different ecology, reproductive strategy and sexual behavior, as it is the case of sea urchin, C. elegans and Human. Sea urchins are members of the Phylum Echinodermata. They are dioecious, having separate male and female sexes, and have five gonads (fivefold symmetry) with a single duct, rising from the upper pole to open at a gonopore lying in one of the genital plates surrounding the anus . The gonads are lined with muscles underneath the peritoneum, that allow the animal to squeeze its gametes through the duct and into the surrounding sea water, where fertilization takes place . Sea urchin spermatozoa acquire the motility only once they reach sea water . To increase the success rate of fertilization the male gametes are attracted by the molecules dispersed by the homologous oocytes. The outer investment of oocytes, the egg jelly, contains short sperm-activating peptides (SAPs) that bind to specific sperm receptors and dramatically alter the metabolic rate and motility of spermatozoa, driving them towards the oocytes . When sperm encounter the egg jelly, the exocytosis of the acrosomal vesicle, the acrosome reaction (AR), occurs and the pHi-dependent polymerization of actin leads to the extension of acrosomal tubule, which exposes a new bindincovered membrane which will fuse specifically through a receptor to the egg [6,7]. The nematode C. elegans is an unsegmented, 1 mm long, vermiform and bilaterally symmetrical worm, living in the ground, having hermaphrodites and males. In particular the 99.95% of the individuals are hermaphrodite and the 0.05% of the population is composed by males. Males have a single-lobed gonad, a vas deferens, and a tail specialized for mating . Hermaphrodites have two ovaries, an oviduct, a spermatheca (a kind of chamber where oocytes are fertilized by sperm), and a single uterus. C. elegans eggs are laid by the hermaphrodite and, after hatching, they pass through four juvenile stages (L1–L4). Hermaphrodites produce all their sperm in the L4 stage (150 sperm per gonadal arm), then they switch over to produce oocytes. The sperm are stored in the same area of the gonad as the oocytes until the first oocyte pushes the sperm into the spermatheca. The male can inseminate the hermaphrodite, which will preferentially use male sperm (both types of sperm are stored in the spermatheca). When self-inseminated C. elegans will lay approximately 300 eggs, on the contrary when inseminated by a male, the number of eggs can exceed 1.000 . C. elegans spermatozoa are characterized by the absence of the flagellum and the lack of the acrosome, while they content many membranous vesicles, the Membranous Organelles (MOs). The motility is acquired by a process (the spermiogenesis or sperm activation), in which the spherical spermatid extends a pseudopod, conferring to the cell an amoeboid motility. This event is made possible by an important reorganization of the cytoskeleton (the C. elegans spermatozoa content MSP instead of actin) and of membrane involving the remodeling of specific microdomains .
In Human the reproduction is obtained by internal fertilization and the spermatozoa were ejaculated in vagina. Male gametes travel along the female genital tract until they reach the oviduct, where they reside for a relatively long period (from hours to days) thus acquiring the fertilizing ability. This process, called “the capacitation”, involves virtually all the cellular structures: the intracellular calcium concentration rises [10,11], the protein phosphorylation pattern changes [12,13], the actin cytoskeleton reorganizes , the plasma membrane and the outer acrosome membrane become more instable and gradually acquire the ability to fuse each other [15,16], the spermatozoa motility is hyperactivated ).
Aim of the present work was to carry out a comparison among the different biochemical events occurring during sperm activation in these organisms, adopting the graph theory formalism. The molecules that interact each other during sperm activation in sea urchin, C. elegans and Human were represented by a mathematical object called “graph” , constituted by a variable number of nodes (the molecules) linked by edges (the interactions) and originating a network, thus, the statistical analysis of the main topological parameters of the networks was carried out.
Database realization, network construction
The data about human sperm activation were obtained from the already realized database  after updating. Since at present a database containing the data about sperm activation of sea urchin (the available data are in particular referred to S. purpuratus,) and C. elegans does not exist, a new database was realized using Microsoft Office Excel 2003. The available information was obtained from peer-reviewed papers from PubMed (www.ncbi.nlm.nih.gov/pubmed/). As a reference were used the data published during latest 10 years. In some cases the record did not represent a single molecule but complex events, such as “membrane fusion” or “protein tyrosine phosphorylation” because all the single molecular determinants of the phenomenon are still unknown.
The database fields were:
Source molecule: representing the molecule source of interaction.
Interaction: representing the nature of interaction (activation, inhibition).
Target molecule: representing the molecule target of the interaction.
Biological function: representing the functional meaning or the contest of interaction (glycolysis, lipid remodelling, oxidative phosphorylation).
Species: representing the species in which the interaction was described.
Reference: representing the bibliographic source of information.
Notes: representing all the notation such as the presence of synonyms or the intracellular location, if relevant, or the explanation of complex cellular events.
The data, extracted from the databases, were used to build the networks using the Cytoscape 2.6.3 software (www.cytoscape.org). The network representing the elements common to the three networks was realized, using the Cytoscape plugin Network Analyzer (function “Compute intersection” of “Compare two networks” menu). This new network was called CE (Common Elements) network.
The networks were spatially represented using the Cytoscape Degree Sorted Circle Layout: all nodes with the same number of links are located together around the circle (Cytoscape User Manual, https://www.cytoscape.org/manual/Cytoscape2_6Manual.html). The node size was proportional to the number of connection and the node color gradient was dependent on the closeness centrality. This parameter is computed as: Cc(n) = 1/avg(L(n, m)), where L (n, m) is the length of the shortest path between two nodes n and m. The closeness centrality of each node ranges from 0 to 1 and it is a measure of how fast information spreads from a given node to the other nodes.
The statistical and topological analyses of networks were carried out, in agreement with Bernabò [1,2] considering the networks as directed by the Cytoscape plugin Network Analyzer (www.med.bioinf.mpi inf.mpg.de/netanalyzer/help/2.6.1/index.html).
The networks representing spermatozoa activation of sea urchin, C.elegans and Human spermatozoa are showed in Figure 1, 2 and 3 and their main topological parameters are tabulated in Table 1. In all cases the distribution of node linkages follows a power law, represented by the generic equation:
y = a x-b
|Sea urchin||C. elegans||Human|
|Avg. n° neighbours||2.740||2.620||2.662|
|Char. path length||8.128||7.816||6.546|
The number of nodes represent the total number of molecules involved; the number of edges represents the total number of interactions; the clustering coefficient is calculated as CI=2nI/k(k–1), where nI is the number of links connecting the kI neighbours of node I to each other; the network diameter is the largest distance between two nodes; the Averaged n° neighbours represents the mean number of connections of each node; the Char. path length gives the expected distance between two connected nodes.
Table 1: Main topological parameters of Sea urchin, C. elegans and Human spermatozoa activation networks.
The r, R2 and b coefficients of each network are tabulated in Table 2.
|Sea urchin||C. elegans||Human|
Table 2: Results of power law fitting of IN and OUT sea urchin, C. elegans and Human spermatozoa activation networks.
The clustering coefficient distribution does not follow a power law, as demonstrated by the results of power law fitting of clustering coefficient distribution (r = 0.598 in sea urchin; r = 0.100 in C. elegans; r = 0.183 in Human).
The most connected nodes (the hubs) in the networks (as shown in Table 3) were [Ca2+]i , cGMP, [H+]i, aCMP, ATP for sea urchin; [Ca2+]i, [H+]i, Motility, ATP, Vesicle fusion, NADH, NAD+, Pseudopod extension for C. elegans; [Ca2+]i, ATP, Tyr phosphorylation, PKA for Human.
|Sea urchin||C. elegans||Human|
|Node||N° of links||Node||N° of links||Node||N° of links|
Table 3: Most connected nodes (the hubs) of sea urchin, C.elegans and Human spermatozoa activation networks.
In all the three networks about the 40-45% of nodes have two links.
A new network showing the common elements present in the three networks has been created by computing the intersection of the existing networks (CE network). CE network is represented in Figure 4 and its topological parameters are showed in Table 4. The nodes belonging to this network are involved in energy metabolism (glicolysis and mitochondrial phosphorylative oxidation) or in signal transduction ([H+]i, PKA, PKC, [Ca2+]i).
|Avg. n° neighbours||2.508|
|Char. path length||8.158|
The number of nodes represents the total number of molecules involved; the number of edges represents the total number of interactions; the clustering coefficient is calculated as CI = 2nI/k(k– 1), where nI is the number of links connecting the kI neighbours of node I to each other; the network diameter is the largest distance between two nodes; the Averaged n° neighbours represents the mean number of connections of each node; the Char. path length gives the expected distance between two connected nodes.
Table 4: Main topological parameters of CE network.
This work was carried out to modelize and compare the signaling machinery of sperm activation processes in organisms belonging to different Phyla such as sea urchin, C. elegans and Human. To this aim, the biological networks representing these events were realized and their main topological parameters were analyzed. This bioinformatics-based approach is a powerful tool to describe biological organization and function of cellular components as well as to understand the principles driving their evolution  and has been already used to represent the biochemical events involved in human  and boar  spermatozoa activation.
First of all, it was evident that all the networks have virtually the same topological features (see Table 1). This finding is per se really interesting and poses an important question: why so different organisms, living in completely different environments and displaying different biological characteristics (see Table 5) share a common topology for sperm activation?
|Sea urchin||C. elegans||Human|
|Time for sperm activation||Seconds||Days||Hours to days|
Table 5: Main biological characteristics of reproduction and spermatozoa in sea urchin, C. elegans and Human
To answer this question it is necessary to evaluate the biological significance of the networks topological parameters. The most elementary characteristics are the node degree and the clustering coefficient. The first, also known as “connectivity”, k, indicates how many links the node has to other nodes. It makes possible to define the node degree distribution, P(k), which represents the probability that a selected node has exactly k links. The clustering coefficient CI = 2nI/k(k–1),where nI is the number of links connecting the kI neighbours of node I to each other, measures the network tendency to develop clusters of nodes: more the clustering coefficient is high, more the presence of clusters increases. On the basis of these topological indexes, the networks were classified in: random networks, scale free networks and hierarchical networks. In random networks, described by the Erdös–Rényi (ER) model, the node degrees follows a Poisson distribution, as a consequence the most of nodes have approximately the same number of links (that defines the network’s scale) and the clustering coefficient is independent of the nodes degree. The scale-free networks (Barabási– Albert , BA, model) are characterized by a power-law degree distribution of the number of links per node: a relatively small number of nodes is highly connected (the hubs) and the most of nodes is scarcely linked, thus it does not exist a “typical” node (scale free topology). In this case the clustering coefficient (CC) is independent of the number of links per node. In hierarchical networks the scale-free topology and the local clustering coexist: the clustering coefficient is higher in the most linked nodes and, consequently, its distribution follows a power law.
The analysis of topological parameters carried out on of sperm activation processes of sea urchin, C. elegans and showed that these networks are preferable to the scale-free networks, as demonstrated by the power law that links the number of edges to the node frequency and the dispersion of clustering coefficient in agreement with the BA model.
From a biological point of view, it is possible to speculate that this particular behavior could offer an important evolutionary advantage: the robustness against random failure. In fact a random perturbation, in most of cases, will involve the most frequent typology of nodes, i.e. the less connected ones, with negligible consequences on network architecture [1,20,21], that is on whole cellular function. In our model the probability that one of them will be involved in random damages is <5%.
Another interesting finding is the very low values of clustering coefficient which suggest that the architecture of signaling system privileges the rapidity of signal transduction instead the redundancy. Moreover the values of characteristic path length that measures how many links it is necessary to pass through to travel between two nodes, ranging between 6.5 and 8.1 are in line with this hypothesis. These two parameter, kept together with the finding that the higher percentage of nodes display two edges (one input and one output), seem to strengthen the idea that sperm activation is a highly efficient signaling process. In fact:
- each node receives an input signal and transfers an output response;
- any molecule can interact with any other in a small number of passages, thus, the loss of information due to the signal decrease is minimized and the signal efficiency is maximized;
- any local perturbation in signaling system could reach the whole network in a short time increasing the system responsiveness to intracellular and extracellular stimuli.
In the same context it is important to point out as the terminal events (such as “protein phosphorylation” or “membrane fusion” or “motility”) are highly linked. Reasonably this is due to the redundancy of biochemical signaling, as a safety strategy to overlap partial failure of the system. The most linked node was, in all the cases, the same: [Ca2+]i. This finding is not surprising, in fact, as in many other cellular systems, in spermatozoa activation Ca2+ behaves as a second messenger.
More in detail in sea urchin the Ca2+ enters the sperm cell through voltage dependent channels (Cav1.2 or 1.3) or through cAMP or cGMP (SpHCN1 or 2) gated channels. In the intracellular signalling cascade this ion is involved in the control of cAMP, cGMP and in the PIP2/IP3 pathway . In C. elegans [Ca2+]i is controlled by membrane channels and by sequestering the ion in intracellular stores and it is involved in PIP2/IP3 signaling pathways regulating the MOs fusion and the onset of motility . In Human it is known that during capacitation the [Ca2+]i increases and the capacitation does not take place in Ca2+ absence [10,23]. Four major Ca2+ clearance mechanisms are described in mammalian spermatozoa, two acting on the plasma membrane and two acting on intracellular organelles. The plasma membrane Ca2+- ATPase exports a cytoplasmic Ca2+ ion and imports one or two extracellular protons at the expense of ATP. When [Ca2+]i is elevated, the plasma membrane Na+- Ca2+ exchanger operates in forward mode exporting an intracellular Ca2+ ion and importing approximately three Na+ ions at the expense of the Na+ gradient [24,25]. The best characterized organellar clearance mechanism involves the sarcoplasmic-endoplasmic reticulum Ca2+- ATPase pumps and the mitochondrial Ca2+ uniporter . During the sperm capacitation the Ca2+ acts converting extracelluar stimuli to chemical response in a myriad of molecular system, such as, protein kynase C (PKC), protein kynase A (PKA), actin, and many others.
The other most connected nodes belong to different classes:
• Molecules involved in energetic balance, such as ATP. It is the main energetic source for spermatozoa, where the metabolic energy production derives exclusively from the glycolysis and from mitochondrial oxidative phosphorylation .
• [H+]i: in sea urchin and C. elegans the pH is a key element in control of many biochemical events, such as the membrane polarization, the rearrangement of cytoskeleton proteins (MSP) and the motility [9,27]. It is important to note that in Human and mammals spermatozoa the pH modification is an early event in the acquisition of full fertilizing ability .
• Molecules playing a key role in signal transduction and amplification: cAMP, cGMP. These two cyclic nucleotides are ubiquitous second messengers and concur to the control of cellular signal transduction by modulating the activity of membrane channels, kinases, phosphatases and many other molecules.
• Complex events whose biochemical determinants are not completely known (“Tyr phosphorylation” in Human, “Motility”, “Pseudopod extension”, “Vesicle fusion” in C.elegans).
The analysis of topological parameters of the CE networks, showed that the elements common to the three networks are those representing the scaffold of spermatozoa energetic metabolism (in all the cases the energy is provided, as ATP, by glicolysis and mitochondrial phosphorylative oxidation) or the signaling system [Ca2+]i, [H+]i, PKA, PKC.
From these data, it is evident that spermatozoa activation in sea urchin, C. elegans and Human recognizes similar biochemical determinants. This finding strengths that some signaling molecules and metabolic processes (such as [Ca2+]i, glycolysis, mitochondrial oxidative phosphorylation) have an ubiquitous role in spermatozoa signaling and biochemistry; in this context it will be interesting to compare the studied networks with the data concerning different cellular systems undergoing an activation process (for instance neurons, leucocytes, myocardiocytes)
More interesting is the evidence that the networks topology of the three networks is virtually identical displaying important features:
- robustness against random failure;
- signaling fastness, rapidity and efficiency.
Since, as Dobzhansky  famously stated, “nothing in biology makes sense except in the light of evolution”, it is possible to hypothesize that the similarity in the architecture of this process in such different organisms could be due to the maintenance of similar ancestral mechanisms that offer important evolutionary advantages.
In our opinion these findings in one hand strength the usefulness of bioinformatics modeling in the study of the biochemical asset of different organisms from a comparative and evolutionary point of view, in the other hand could contribute to the knowledge of the spermatozoa activation, i.e. of a process of pivotal importance for the survival of living beings.
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals