Zulfi Nur Amrina Rosyada1,2,3*
Berlin Pandapotan Pardede1,2
Ekayanti Mulyawati Kaiin2Muhammad Gunawan2Tulus Maulana2Syahruddin Said2Ligaya I. T. A Tumbelaka1Dedy Duryadi Solihin4
Mokhamad Fakhrul Ulum1
Bambang Purwantara1*- 1Division of Reproduction and Obstetrics, School of Veterinary Medicine and Biomedical Sciences, IPB University, Bogor, Indonesia
- 2Research Center for Applied Zoology, National Research and Innovation Agency (BRIN), Bogor, Indonesia
- 3Division of Veterinary Anatomy, Faculty of Veterinary Medicine, Universitas Airlangga, Surabaya, Indonesia
- 4Department of Biology, Faculty of Science, IPB University, Bogor, Indonesia
Proteins assist sperm mature, transit the female reproductive tract, and recognise sperm oocytes. Indigenous Indonesian bulls, Madura bulls, have not been studied for reproductive proteomics. As local Indonesian beef livestock, Madura cattle assist in achieving food security; hence, their number must be improved. Thus, the identification of molecular proteomics-based bull fertility biomarkers is needed. This study aimed to characterise the sperm fertility function of the superior Madura bull (Bos indicus × Bos Javanicus) spermatozoa proteome. Frozen semen from eight Madura superior bulls (Bos indicus × Bos javanicus) aged 4–8 years was obtained from the artificial insemination centre (AIC) in Singosari and Lembang. Madura superior bulls are those that have passed the bull breeding soundness evaluation. Frozen sperm were thawed and centrifuged at 3000 × g for 30 min. Proteins in sperm were characterised through proteomic analysis using liquid chromatography–tandem mass spectrometry (LC–MS/MS). The resulting gene symbols for each protein were then subjected to bioinformatics tools, including UniProt, DAVID, and STRING databases. Regarding sperm fertility, the analysis revealed that 15 proteins were identified in the sperm of Madura bulls. Amongst the identified proteins, the superior Madura bull sperm contained several motilities, energy-related proteins, and chaperone proteins. A substantial portion of characterised proteins are linked to metabolic pathways and the tricarboxylic acid (TCA) cycle, contributing to sperm energy production. In conclusion, the first in-depth proteome identification of sperm related to sperm quality and bull fertility of a unique indigenous Madura breed of Indonesia was performed using the LC–MS/MS proteomic method. These findings may serve as a reference point for further studies related to the functions of bovine sperm and biomarkers of fertility and sperm quality.
1 Introduction
Recent developments in proteomics have significantly impacted our understanding of how sperm become fertile (1). One of the most differentiated cells is the sperm head, which has a highly compacted chromatin structure and an enormous midpiece that contains the machinery required to interact and transfer paternal genetic and epigenetic information to the oocyte (2). Due to their high level of differentiation, sperm is suitable for studying proteomic sections, such as the sperm membrane, which is the most crucial part because of its role in interacting with the environment and the oocyte (3). The fusion of a sperm and an oocyte requires complex membrane modifications of the sperm (4). Spermatozoa proteins that regulate normal/abnormal sperm function have been identified through proteomic investigation (5).
Thus, bull fertility is linked to many proteins involved in spermatogenesis. The significance of specific proteins in controlling sperm quality and fertilization is unknown, although their availability and quantity may alter sperm fertility (6). According to Peddinti et al. (1), many bovine fertility protein markers have been reported. Major semen proteins, such as binding sperm protein (BSP) in Frisian Holstein cattle (6), zona pellucida binding protein (ZPBP) in Bali Polled bulls (7), A-kinase anchoring protein 3 in buffalo (8), and osteopontin in Limousin cattle (9), have been widely reported. The sperm proteome profile of mammals has been studied, including pigs (10), equines (11), sheep (12), and bovines (1). Nevertheless, the sperm proteome remains unknown in Madura bulls.
Madura cattle (Bos indicus × Bos javanicus) are native to Madura, Indonesia. Small-scale producers raise Madura cattle as working cattle to maintain the regional culture, which still values Sonok and Karapan cattle (13). In addition, Madura cattle are currently being developed as beef cattle using artificial insemination and reproduction technology. It also ensures food security for the Indonesian population. In particular, these cattle are highly valued for their unique reproductive traits, with regular calving even under low-input regimes and in dry and arid regions (14). However, it remains unknown whether the proteome level of Madura bull sperm affects sperm function and fertility. Proteomics offers various methods. However, complex protein analysis using LC–MS has become the preferred analytical technique for quantitative proteomics (15, 16).
Ultimately, acquiring knowledge about the proteome of spermatozoa can provide a holistic understanding of reproductive processes, not only specific to the Madura bull breed but also applicable to bovine species. This study aims to provide a comprehensive profiling of the sperm proteins associated with fertility function in the superior Madura bull (Bos indicus × Bos javanicus) spermatozoa proteome. Additionally, we used the literature to explain how Madura bull sperm proteins affect bull fertility. Further, proteomic approaches have made it possible to identify proteins that have the potential to function as indicators of male fertility, as well as proteins that are involved in the functional characteristics of sperm.
2 Materials and methods
2.1 Frozen semen samples
A total of 8 superior Madura bulls (Bos indicus × Bos javanicus) aged 4 to 8 years with spermatozoa motility >70%, based on secondary data from each AIC in Lembang and Singosari, Indonesia, were included in this study. Frozen semen from each of the eight Madura superior bulls was obtained from the National AIC in Singosari and Lembang Bank of Semen. Madura superior bulls are those that have passed the bull breeding soundness evaluation. The Animal Care and Use Committee excluded this study from its ethical review because artificial vaginal semen collection did not alter the physiology of animals. The study was directed by veterinarians from both institutions and followed SNI ISO 9001: 2015 No. 824 100 16072 at Lembang AIC and SNI ISO 9001: 2015 No. G.01-ID0139-VIII-2019 at Singosari AIC. The ethics committees of Lembang AIC and Singosari AIC provided ethical guidance and sanctions for the responsible collection of bull sperm. Additionally, an experienced bull technician collected sperm using an artificial vagina.
Semen samples from Lembang AIC were processed according to the AIC’s standard operating procedure (SOP), which used skim milk as the extender. The formulation of the skimmed milk diluent in a volume of 1000cc consisted of 100g of skimmed milk, 960cc of distilled water, and an antibiotic solution containing 3000000 IU of penicillin, 3g of streptomycin, and 30cc of distilled water. The ratio of skimmed milk diluent to antibiotic was 100:1. In contrast, the semen samples obtained from Singosari were cryopreserved utilising tris egg yolk as a cryoprotectant. The tris egg yolk extender was formulated using the following components: 20% egg yolk, 1.6% tris aminomethane, 1.4% lactose, 2.5% raffinose, 0.9% citric acid, and an antibiotic mixture consisting of 1000000 IU/L of penicillin, 1g/L of streptomycin, and 1g/L of distilled water. The AIC in Lembang and Singosari employ extenders consisting of skim milk and tris egg yolk supplemented with antibiotics. The extender utilised in every artificial insemination facility serves as a cryoprotectant, safeguarding sperm cells during the freezing process.
Additionally, it enhances semen volume, sustains sperm viability, and regulates sperm pH. Consequently, the inclusion of an extender is imperative in the manufacture of frozen semen to ensure the preservation of sperm quality. The cryopreserved sperm samples were transferred with a specialised transport container that maintained a temperature of −196°C through liquid nitrogen. The pieces were placed within the container for subsequent examination.
2.2 Sperm protein isolation
The collected semen samples of about ive straws of frozen semen from each of the eight bulls were washed in 2 mL of phosphate buffer saline and centrifuged twice (3000 × g for 10 min, 4°C) to clean and separate the spermatozoa, seminal plasma, and extender of frozen semen. The resulting sperm cell pellet was then resuspended in cell lysis buffer (2% SDS in 62.5 mM Tris–HCl, pH 6.8, 1.0 mM phenylmethanesulfonyl fluoride, and 23 mM benzidine as a protease inhibitor) and stored at −20°C for protein extraction. The sperm pellet was vortexed for 10 min before protein extraction. Next, sonication was used thrice for 20 s each to dissolve the protein. After centrifuging the cell lysates (10,000 × g, 10 min), protein lysates were isolated. The Bradford method (17) was used to estimate the total protein yield of the sperm protein lysate. The Bradford protocol was performed using the Coomassie protein kit instructions (Merck, Darmstadt, Germany). The data were processed using ThermoScan RE software version 3.2 Multiskan Go (Thermo Fisher Scientific, Waltham, Massachusetts, United States). Pools of sperm protein containing equal amounts from three samples were frozen at −20°C for later use. A summary of the experimental design is shown in Figure 1.
Figure 1. Schematic workflow for the proteomic analysis of Madura bull spermatozoa. (A) Sperm protein extraction. (B) Sperm protein identification.
2.3 Fractionation by one-dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis (1D-SDS–PAGE) and in-gel digestion
Following the separation of 35 μg of sperm protein (an equivalent quantity of protein amalgamated from three cases) by 12.5% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS–PAGE), the gel was stained with colloidal Coomassie brilliant blue R-250 (CAS number: 6104-59-2) from HiMedia. Subsequently, the gel bands were excised and digested using a previously described protocol (18). The proteins were reduced using 0.5 M dithiothreitol (DTT) at 55°C for 25 min, followed by alkylation with 14 mM iodoacetamide (IAA) for 40 min at room temperature without light. Gel fragments were treated with trypsin (Promega, Madison, Fitchburg, WI, USA) at a ratio of 1:50 (trypsin: protein) at 4°C for 30 min. The gel fragments were incubated for 18 h at 37°C in a ThermoMixer shaker incubator block (Eppendorf). The reaction in each vial was terminated by adding 1% trifluoroacetic acid (TFA). An extraction buffer (2% TFA in 20% acetonitrile) was used to extract peptide digests from the gel pieces, followed by 70% acetonitrile. The obtained peptide digests were subjected to a drying process and subsequently purified using a C18 spin column for desalination (Thermo Scientific, Pierce Biotechnology, N Meridian Rd., Rockford, IL, USA) and then preserved at −20°C until LC–MS/MS analysis.
2.4 Liquid chromatography–tandem mass spectrometry (LC–MS/MS) analysis
The peptides obtained from each sample were subjected to LC–MS/MS analysis using an Ultimate 3000 Nano LC system coupled with a Q-Exactive Plus Orbitrap HRMS system (Thermo Fisher Scientific, Bremen, Germany). The peptides were introduced onto a preanalytical column with dimensions of 75 μm ID, 15 cm length, and 100 pore size, packed with Acclaim PepMap C18 2 μm particles. The solvent used for loading was solvent A, consisting of 0.1% formic acid, and the flow rate was set at 300 nL/min. A Q-Exactive Plus Orbitrap mass analyser was used for subsequent analysis. The peptides were subsequently resolved on an analytical column (50 cm 75 m ID, Pep Map RSLC C18 2 m) at a flow rate of 300 nL/min for 120 min using an increasing gradient of 5–35% solvent B (98% acetonitrile and 0.1% formic acid). Peptide signals were collected using an LTQ-Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). The MS spectra within the 200–2000 m/z range were obtained using an Orbitrap analyser with a resolution of 30,000 (at m/z 400). Subsequently, ten precursor MS scans were conducted using collision-induced dissociation fragmentation at 35% normalised collision energy (19).
2.5 Protein database searching and bioinformatics analysis
Protein analysis was performed on proteins that contained at least one unique peptide. The exclusion of specific peptides based on the “minimum two-peptide rule may result in the perpetual disregard of genuine peptides, and proteins must have a sequence score HT > 0. Thus, we manually confirmed single peptide spectral match (PSM)-identified protein MS/MS spectra. Protein data were acquired and analysed using the UniProt bovine protein database (http://www.uniprot.org) and Proteome Discover version 2.2 software (Thermo Fisher Scientific). Furthermore, the discovered proteins from spermatozoa were also put into functional groups based on Gene Ontology (GO) using the DAVID web resource (https://david.ncifcrf.gov/tools.jsp). The interactions between the proteins were obtained using STRING version 11.0 (https://string-db.org/).
3 Result
3.1 Protein characterisation of Madura bull spermatozoa
LC–MS/MS identified Madura bull sperm protein markers of fertility and environmental compatibility. All the databases contained in the database were only Bos taurus origin proteins since Bos indicus × Bos javanicus proteins are limited. This study analysis identified 15 proteins in Madura bull (Bos indicus × Bos Javanicus) sperm using UniProt and DAVID software (Table 1). Fifteen proteins were identified that were associated with sperm function categories (Table 1), including adenosine triphosphate (ATP) synthesis activity, cellular metabolic processes, cilia/flagella, sperm motility, capacitation and acrosome reaction, sperm-egg fusion, spermatogenesis and fertilization, and chaperone proteins. Functional annotation of proteins from Madura bull spermatozoa was further carried out using DAVID software and divided into several categories: “Biological Process” (BP), “Cellular Component” (CC), and “Molecular Function” (MF) (Figure 2). Although different categories were created for each division, the most numerous were those involved in the generation of precursor metabolites and energy (16%), urine ribonucleotide metabolic processes (16%), ATP metabolic processes (13%), and pyruvate metabolic processes (10%) in the case of BP (Figure 2A); cytoplasm (24%), supramolecular fibre (21%), and cell surface and microtubule (10%) in CC (Figure 2B); and ATP-binding proteins (16%), ATPase activity (12%) and unfolded protein binding (8%) in MF (Figure 2C).
Table 1. Proteins identified in Madura bull sperm associated with sperm function.
Figure 2. Bar chart representing the gene ontology annotations of proteins identified in Madura bull spermatozoa according to biological processes (A), cellular components (B), and molecular functions (C).
3.2 Protein interaction networks and pathway enrichment analysis
All proteins identified in Madura bull spermatozoa were then searched using STRING software (version 11.0) for protein–protein interaction network analysis. Light blue lines show protein nodes, whereas pink lines reflect experimentally established interactions. The expected interactions are green for the gene neighbourhood, red for gene fusion, and dark blue for gene co-occurrence. Other protein relationships are shown by light green text mining, black co-expression, and blue protein homologies. A coloured node shows the query protein, and a white node leads the second shell of the interaction. Proteins with unknown 3D conformations had empty nodes, whereas those with known conformations had full nodes. The interactions between the 15 identified proteins are shown in Figure 3. Fifteen proteins were directly or indirectly connected through one or more interacting proteins, indicating functional links. We categorised proteins using the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway terminology to study Madura bull sperm pathways. As expected, many of the discovered proteins were involved in metabolic pathways. The discovered proteins were associated with the TCA cycle (Figure 4). The sperm head, mitochondria (mid-piece), and flagellum/tail expressed the majority of these proteins (Figure 5).
Figure 3. STRING protein–protein interaction network showing the interactions of the spermatozoa proteins identified in Madura bull sperm.
Figure 4. Pathway analysis using the KEGG pathway database identified proteins involved in the TCA cycle. The red star denotes the detected proteins.
Figure 5. Distribution of sperm protein of Madura bull (Bos indicus × Bos javanicus).
4 Discussion
Fifteen proteins found in the sperm of Madura bulls were found to be related to fertility based on LC–MS/MS analysis (Table 1). The utilisation of known molecular weight standards for calibration is a fundamental aspect of size-exclusion chromatography that provides insights into the distribution of molecular weights. The current study has demonstrated that Madura bull sperm exhibits a molecular weight distribution of proteins within the 40 to 180 kDa range. This study utilised GO annotation to ascertain the cellular localisation, biological processes, and molecular functions of proteins specific to the sperm of Madura bulls.
4.1 Protein related to sperm fertility function
Effect of sperm motility on semen quality Male infertility is often attributed to reduced sperm motility. Therefore, the motility proteins of the Madura bull breed must be identified. We found many of these proteins in the spermatozoa. SPAM1 is a transmembrane protein present in ejaculated bull spermatozoa. It originates from the testis and occupies the anterior head region, where the orientation of the C-terminus with the zonal-binding domain facilitates interaction with zona pellucida after the acrosome reaction (20). This protein has also been detected in male sperm (21). The roteomee analysis carried out in this study from the sperm of Madura Bull also revealed the existence of the hyaluronidase protein, also known as SPAM1, isoform 62.4 kDa.
The flagellum of mammalian spermatozoa has intricate supplemental structures around the core axoneme, including the outer dense fibers (ODFs). ODFs protect the sperm tail against shear stress that may arise during epididymal transit and ejaculation. Four major proteins (ODF1, ODF2, ODF3, and ODF4) are amongst more than 14 polypeptides produced by mammalian ODFs. In mice, disrupting ODF2 expression lowers sperm motility, which is consistent with the characteristics of asthenozoospermia (22). Moreover, the cervical midpiece junction requires the presence of ODF2, which comprises part of the centrosome and is separated from the flagellum (23). Our analysis identified ODF2 in Madura bull spermatozoa.
In addition, we found that enolase protein, previously discovered in the plasma membrane of sperm from Bali bulls (24), was present in sperm samples from Madura bulls. During glycolysis and gluconeogenesis, enolase reversibly converts 2-phosphoglycerate into phosphoenolpyruvate. This enzyme produces sperm motility energy (25). Xi et al. (26) found a significant positive correlation between ENO3 and sperm motility parameters in sturgeons. Similarly, He et al. (25) discovered ENO3 in the midpiece and tail of ram sperm, which strongly correlated with sperm motility. The presence of ENO3 expression in Madura cattle sperm and its correlation with sperm function, particularly motility, can be readily explained.
The present study revealed the detectability of AKAP 4 in Madura bull sperm using LC–MS/MS. AKAP4, a marker involved in flagellar structure and motion, is associated with the equine sperm tail (27) and is substantially preserved in mice (28), bulls (29), and humans (30). AKAP4 attaches cAMP-dependent protein kinase A (PKA) to different subcellular locations. In the sperm fibrous sheath, AKAP4 binds PKA and interacts with other proteins to regulate motility. AKAP4, a critical fibrous sheath protein, bundles glycolytic enzymes and phosphorylation-signalling cascade components to supply a localised source of ATP and govern flagellar motion, sperm motility, and hypermotility. Many glycolytic enzymes closely linked to the fibrous sheath at the primary component of the flagellum are required for sperm motility (27).
Tubulins were also found to be another type of protein that contributes to sperm motility that was uncovered in this study. Tubulins are proteins found in the microtubules of sperm (30). Two tubulin proteins were identified as present in the sperm of the studied Madura bull samples. These proteins are known as TUBA1A and TUBB4B. These proteins are necessary for cilia and flagella formation; hence, their importance cannot be overstated.
4.2 Chaperone proteins: heat shock proteins
Heat stress affects cattle productivity and reproduction. Heat stress damages developing spermatozoa. However, heat tolerance can be inherited by offspring (3). According to the relevant literature, Madura cattle are a native breed of beef cow thought to have developed on the dry and barren island of Madura (3). Madura Island receives 1600 mm of rain annually and is hot and arid (31). HSPs are linked to heat tolerance and reproductive performance (3, 32). In this study, HSPA9 and HSPA2 were found to be ubiquitously expressed. HSPA9, found in sperm mitochondria, significantly reduces the effects of heat stress in tropical cattle (33, 34). It is also known to regulate sperm motility through metabolic processes, ATP binding, and folding of mitochondrial proteins (35).
HSPA2 (also known as HSP70-2) was discovered in this study, which is consistent with previous findings in Australian Brahman bulls (Bos indicus) (36) and Zebu (37). HSP70-2 (heat shock protein 70–2) was found to be expressed in the sperm nucleus as well as in various other sperm organelles, including the mitochondria and flagellum (Figure 5). HSP70-2 is a component of the chromatin structure and contributes to gene regulation, particularly in the folding and unfolding of proteins. This HSPA2 sperm protein may also be crucial in protecting the sperm of tropical acclimated bulls such as Madura from heat or environmental stress. HSP70 is also an excellent predictor of thermotolerance and thermoresistance (38). The potential use of HSP70 as a biomarker for animal fertility and thermotolerance has been suggested (39). Therefore, HSP70-2 is recognised as a dual-functioning gene (3).
4.3 Pathway enrichment
The enriched pathways were also analysed for the identified proteins in the Madura bull spermatozoa. The present study revealed that 16% of the proteins were involved in the biosynthesis of precursor metabolites and energy. Only metabolic pathways were included in the list of the significant ways enriched by the proteins. Energy metabolism is essential for sperm development. Sperm require ATP, which is most likely to maintain morphological changes during the spermatid stage, for the degradation and synthesis of active proteins (40). ATP also serves as the molecular motor that provides the necessary energy for flagellar movement in all kinetics-related biochemical activities. The idea that mitochondria in the sperm midpiece drive mammalian sperm is evolving, as evidence suggests that glycolysis is the preferred metabolic pathway to maintain sperm motility in many species (41). Thus, the current study’s results are in line with earlier research showing that ATP synthases are hub proteins found in the sperm of rams that are comparatively less fertile and may regulate how much ATP is used by sperm throughout the fertilization process (42). Existing data indicate that ATP production in sperm is facilitated by both glycolysis and mitochondrial respiration. These processes are interdependent and regulate sperm function based on the availability of energy substrates in the surrounding environment. An alternative glycolysis pathway for ATP production is also found in stallion spermatozoa (43). This study presents novel findings that suggest that Madura sperm are dependent on mitochondrial ATP synthesis, specifically regarding the roles of ATPF51A and ATP5B in sperm function and fertility.
Further, glucose utilisation is restricted during spermatogenesis, and lactate and pyruvate are favoured as substrates for energy synthesis (19). The generation of ATP in mammalian sperm is attributed to the energy metabolism pathway. This pathway exhibits subcellular compartmentalisation, with oxidative phosphorylation (OXPHOS) predominantly occurring in the sperm midpiece and glycolysis in the principal piece (40). Glycolysis is a significant mechanism for facilitating the transport of ATP along the flagellum. Westhoff and Kamp (44), as well as Welch et al. (45), have reported the presence of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) in significant quantities within the fibrous sheath of sperm from different mammalian species, including humans. GAPDH is a NAD-dependent glycolytic enzyme responsible for facilitating the conversion of glyceraldehyde 3-phosphate (GAP) to 1,3-biphosphoglycerate (1,3 BPG). Notably, one of the isoforms of GAPDH, known as GAPDHS, has been identified in this study.
Moreover, Zini et al. (46) found that an augmentation in the activity of pyruvate dehydrogenase complexes can increase pyruvate-lactate usage, thus producing a more significant amount of NADH. NADH oxidase can utilise the NADH produced to generate reactive oxygen species (ROS) necessary for the capacitation process. NADH oxidase activity and ROS formation have also been associated with the motility of human spermatozoa. In this way, GAPDHS contributes significantly to energy metabolism. Elkina et al. (47) found GAPDHS in the testes of humans and rats; the enzyme is primarily present in the cytoplasm of all spermatogenic cells, whereas it is localised in the sperm tail in the epididymis. This finding provides support for this theory. This previous study by Kumar et al. (48) also established the involvement of extramitochondrial localised pyruvate dehydrogenase complex (PDHA) in the process of sperm capacitation, as well as the significance of pyruvate in the overall energy metabolism of mammalian sperm. The metabolic pathway encompassed the GAPDHS and PDHA2 proteins, both detected in the sperm of Madura bulls in the present study.
In agreement, energy metabolism plays a critical role in facilitating sperm function. The TCA cycle is the primary energy source for spermatids, although glycolytic and pentose phosphate pathways also play a role in energy synthesis in spermatozoa. The TCA cycle generates adenine, which is converted to ATP by the electron transport chain of the OXPHOS pathway. These proteins may be involved in creating the acrosome and activities there, both of which require energy provided by oxidative phosphorylation. ATP is then delivered to the microtubules responsible for supporting various sperm functions (49). In this study, two candidate proteins related to the tricarboxylic acid cycle pathway were PDHA2 and ACO2.
It was intriguing to discover that specific proteins are involved in so many energy metabolism processes, suggesting that a single protein may not perform glycolysis or OXPHOS as the primary actor in maintaining sperm functionality. The current study shows that proteomic methods may provide a reliable source for sperm protein detection and improve gene ontology comprehension. Therefore, we used high-throughput LC–MS/MS to develop a proteomic profile of Madura bull sperm connected to reproductive characteristics. Preserving superior native germplasm and preventing gene pool depletion requires studying protein differences in different breeds. Thus, investigating Madura bull protein composition can help us comprehend native Indonesian bull sperm biology. Finally, knowing the proteome of ejaculated spermatozoa can help us understand the reproductive process in Madura bull and bovine in general.
5 Conclusion
We examined the sperm proteomes of Madura cows (Bos indicus × Bos javanicus). The use of several methods to prepare samples has made it easier to identify proteins in sperm. Most of the proteins identified in this study are essential for bull breeding. We believe that information about the sperm proteome of the superior Madura bull, a breed with good reproductive traits, will speed up future studies on bull fertility. Additionally, this would enable the creation of molecular tools for accurately selecting bulls to conserve indigenous Indonesian cattle breeds and bovines.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was not required for the studies on animals in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.
Author contributions
ZR: Conceptualization, Data curation, Writing – original draft, Writing – review & editing, Formal analysis, Investigation, Methodology. BP: Data curation, Funding acquisition, Project administration, Validation, Writing – original draft, Writing – review & editing, Formal analysis, Investigation, Methodology, Visualization. EK: Writing – review & editing. MG: Writing – review & editing. TM: Writing – review & editing. LT: Writing – review & editing, Conceptualization, Supervision. DS: Conceptualization, Supervision, Writing – review & editing. MU: Conceptualization, Supervision, Writing – review & editing. BP: Conceptualization, Data curation, Funding acquisition, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing. SS: Methodology, Supervision, Formal analysis, Validation, Investigation, Funding acquisition, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by the Indonesian Ministry of Education, Culture, Research, and Technology through Fundamental Research grant number 18789/IT3.D10/PT.01.02/M/T/2023.
Acknowledgments
The authors sincerely thank the Singosari and Lembang AIC for providing bull data. We also thank Wiwit for his invaluable help in the laboratory.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Peddinti, D, Nanduri, B, Kaya, A, Feugang, JM, Burgess, SC, and Memili, E. Comprehensive proteomic analysis of bovine spermatozoa of varying fertility rates and identification of biomarkers associated with fertility. BMC Syst Biol. (2008) 2:1–13. doi: 10.1186/1752-0509-2-19
2. Thompson, SK, Kutchy, NA, Kwok, S, Rosyada, ZNA, Imumorin, IG, Purwantara, B, et al. Review: sperm: comparative morphology and function related to altered reproductive strategies and mammalian fertil. Professional animal. Scientist. (2018) 34:558–65. doi: 10.15232/pas.2018-01748
3. Rosyada, ZNA, Ulum, MF, Tumbelaka, LITA, Solihin, DD, Purwantara, B, and Memili, E. Implications of sperm heat shock protein 70-2 in bull fertility. Vet World. (2022) 15:1456–66. doi: 10.14202/vetworld.2022.1456-1466
4. Nowicka-Bauer, K, Malcher, A, Włoczkowska, O, Kamieniczna, M, Olszewska, M, and Kurpisz, MK. Evaluation of seminal plasma HSPA2 protein as a biomarker of human spermatogenesis status. Reprod Biol. (2022) 22:100597. doi: 10.1016/j.repbio.2021.100597
5. Rahman, MS, Kwon, WS, and Pang, MG. Prediction of male fertility using capacitation-associated proteins in spermatozoa. Mol Reprod Dev. (2017) 84:749–59. doi: 10.1002/mrd.22810
6. Somashekar, L, Selvaraju, S, Parthipan, S, and Ravindra, JP. Profiling of sperm proteins and association of sperm PDC-109 with bull fertility. Syst Biol Reprod Med. (2015) 61:376–87. doi: 10.3109/19396368.2015.1094837
7. Manabi Diansyah, A, Yusuf, M, Latief Toleng, A, Ihsan Andi Dagong, M, Maulana, T, Hasrin, H, et al. The sperms post-thawing quality and proteomic seminal plasma on fertility performance of Bali-polled bull. Adv Anim Vet Sci. (2022) 11:517–525. doi: 10.17582/journal.aavs/2023/11.4.517.525
8. Karanwal, S, Pal, A, Chera, JS, Batra, V, Kumaresan, A, Datta, TK, et al. Identification of protein candidates in spermatozoa of water buffalo (Bubalus bubalis) bulls helps in predicting their fertility status. Front cell. Dev Biol. (2023) 11:11. doi: 10.3389/fcell.2023.1119220
9. Monaco, E, Gasparrini, B, Boccia, L, De Rosa, A, Attanasio, L, Zicarelli, L, et al. Effect of osteopontin (OPN) on in vitro embryo development in cattle. Theriogenology. (2009) 71:450–7. doi: 10.1016/j.theriogenology.2008.08.012
10. Luongo, C, González-Brusi, L, Cots-Rodríguez, P, Izquierdo-Rico, MJ, Avilés, M, and García-Vázquez, FA. Sperm proteome after interaction with reproductive fluids in porcine: from the ejaculation to the fertilization site. Int J Mol Sci. (2020) 21:1–27. doi: 10.3390/ijms21176060
11. Guasti, PN, Souza, FF, Scott, C, Papa, PM, Camargo, LS, Schmith, RA, et al. Equine seminal plasma and sperm membrane: functional proteomic assessment. Theriogenology. (2020) 156:70–81. doi: 10.1016/j.theriogenology.2020.06.014
12. Pini, T, Leahy, T, Soleilhavoup, C, Tsikis, G, Labas, V, Combes-Soia, L, et al. Proteomic investigation of ram spermatozoa and the proteins conferred by seminal plasma. J Proteome Res. (2016) 15:3700–11. doi: 10.1021/acs.jproteome.6b00530
13. Rosyada, ZNA, Tumbelaka, LI, Ulum, MF, Solihin, DD, Kaiin, EM, Gunawan, M, et al. Meta data analysis of conception rate in relation to sperm motility in Madura superior bulls. IOP Conf Ser Earth Environ Sci. (2021) 902:012048. doi: 10.1088/1755-1315/902/1/012048
14. Rosyada, ZNA, Yoelinda, VT, Kaiin, EM, Gunawan, M, Ulum, MF, Tumbelaka, LITA, et al. Sperm osteopontin mRNA expression levels and its correlation on semen quality and fertility in Madura bulls. Biodiversitas. (2023) 24:563–70. doi: 10.13057/biodiv/d240165
15. Ramesha, KP, Mol, P, Kannegundla, U, Thota, LN, Gopalakrishnan, L, Rana, E, et al. Deep proteome profiling of semen of Indian indigenous Malnad Gidda (Bos indicus) cattle. J Proteome Res. (2020) 19:3364–76. doi: 10.1021/acs.jproteome.0c00237
16. Druart, X, and De, GS. Seminal plasma proteomes and sperm fertility seminal plasma proteomes and sperm fertility. Anim Reprod Sci. (2018) 194:33–40. doi: 10.1016/j.anireprosci.2018.04.061
17. Bradford, MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. (1976) 72:248–54. doi: 10.1016/0003-2697(76)90527-3
18. Prasad, TSK, Harsha, HC, Keerthikumar, S, Sekhar, NR, Selvan, LDN, Kumar, P, et al. Proteogenomic analysis of Candida glabrata using high resolution mass spectrometry. J Proteome Res. (2012) 11:247–60. doi: 10.1021/pr200827k
19. Fu, Q, Pan, L, Huang, D, Wang, Z, Hou, Z, and Zhang, M. Proteomic profiles of buffalo spermatozoa and seminal plasma. Theriogenology. (2019) 134:74–82. doi: 10.1016/j.theriogenology.2019.05.013
20. Martin-Deleon, PA. Germ-cell hyaluronidases: their roles in sperm function. Int J Androl. (2011) 34:e306–18. doi: 10.1111/j.1365-2605.2010.01138.x
21. Gómez-Torres, MJ, Sáez-Espinosa, P, Manzano-Santiago, P, Robles-Gómez, L, Huerta-Retamal, N, and Aizpurua, J. Sperm adhesion molecule 1 (SPAM1) distribution in selected human sperm by hyaluronic acid test. Biomedicine. (2022) 10:10102553. doi: 10.3390/biomedicines10102553
22. Zhao, W, Li, Z, Ping, P, Wang, G, Yuan, X, and Sun, F. Outer dense fibers stabilize the axoneme to maintain sperm motility. J Cell Mol Med. (2018) 22:1755–68. doi: 10.1111/jcmm.13457
23. Ito, C, Akutsu, H, Yao, R, Yoshida, K, Yamatoya, K, Mutoh, T, et al. Odf2 haploinsufficiency causes a new type of decapitated and decaudated spermatozoa, Odf2-DDS, in mice. Sci Rep. (2019) 9:14249. doi: 10.1038/s41598-019-50516-2
24. Sumarsono, T, Supriatna, I, Setiadi, MA, Agil, M, and Purwantara, B. Detection of plasma membrane alpha-enolase (ENO1) and its relationship with sperm quality of Bali cattle. Trop Anim Sci J. (2023) 46:36–42. doi: 10.5398/tasj.2023.46.1.36
25. He, Y, Li, H, He, J, and Zhao, X. Heavy ion radiation can promote greater motility and enolase protein expression in ram sperm in in vitro liquid storage. Anim Reprod Sci. (2014) 148:260–6. doi: 10.1016/j.anireprosci.2014.06.019
26. Xi, MD, Li, P, Du, H, Qiao, XM, Liu, ZG, and Wei, QW. Disaccharide combinations and the expression of enolase3 and plasma membrane Ca2+ ATPase isoform in sturgeon sperm cryopreservation. Reprod Domest Anim. (2018) 53:472–83. doi: 10.1111/rda.13134
27. Blommaert, D, Sergeant, N, Delehedde, M, Jouy, N, Mitchell, V, Franck, T, et al. Expression, localization, and concentration of A-kinase anchor protein 4 (AKAP4) and its precursor (proAKAP4) in equine semen: promising marker correlated to the total and progressive motility in thawed spermatozoa. Theriogenology. (2019) 131:52–60. doi: 10.1016/j.theriogenology.2019.03.011
28. Carrera, A, Gerton, G, and Moss, SB. The major fibrous sheath polypeptide of mouse sperm: structural and functional similarities to the A-kinase anchoring proteins. Dev Biol. (1994) 165:272–84. doi: 10.1006/dbio.1994.1252
29. Moss, SB, Turner, RMO, Burkert, KL, Butt, HV, and Gerton, GL. Conservation and function of a bovine sperm A-kinase anchor protein homologous to mouse AKAP82 1. Biol Reprod. (1999) 61:335–42. doi: 10.1095/biolreprod61.2.335
30. Kierszenbaum, AL. Genomic imprinting and epigenetic reprogramming: unearthing the garden of forking paths. Mol Reprod Dev. (2002) 63:269–72. doi: 10.1002/mrd.90011
31. Kutsiyah, F. Population dynamic and productivity of Madura cattle in Corservation area of Sapudi Island. Sains Peternakan. (2017) 15:70–7. doi: 10.20961/sainspet.v15i2.13160
32. Purandhar, K, Jena, PK, Prajapati, B, Rajput, P, and Seshadri, S. Understanding the role of heat shock protein isoforms in male fertility, aging and apoptosis. World J Mens Health. (2014) 32:123. doi: 10.5534/wjmh.2014.32.3.123
33. Aarif, O, Mir, N, Shergojry, S, Ramesha, KP, and Mir, A. Genetic polymorphisms within exon 8, 9, and 10 of heat shock protein (HSP) 90AA1 in deoni cattle. Anim Sci Report. (2014) 8:26–30.
34. Rosyada, Z, Pardede, B, Kaiin, E, Tumbelaka, L, Solihin, D, Purwantara, B, et al. Identification of HSP70-2 and PRM-1 mRNA, proteins, and analyses of their association with fertility using frozen-thawed sperm in Madura bulls. Animal. Bioscience. (2023) 36:1796–805. doi: 10.5713/ab.23.0142
35. Liang, J, Zheng, Y, Zeng, W, Chen, L, Yang, S, Du, P, et al. Proteomic profile of sperm in infertile males reveals changes in metabolic pathways. Protein J. (2021) 40:929–39. Available from:. doi: 10.1007/s10930-021-10013-w
36. Rego, JPA, Crisp, JM, Moura, AA, Nouwens, AS, Li, Y, Venus, B, et al. Seminal plasma proteome of electroejaculated Bos indicus bulls. Anim Reprod Sci. (2014) 148:1–17. doi: 10.1016/j.anireprosci.2014.04.016
37. Zhang, XG, Hong, JY, Yan, GJ, Wang, YF, Li, QW, and Hu, JH. Association of heat shock protein 70 with motility of frozen-thawed sperm in bulls. Czech JAnim. (2015) 60:256–62. doi: 10.17221/8239-CJAS
38. Ul, HF, Nawaz, A, Rehman, MS, Ali, MA, SMR, D, and Yang, C. Prospects of HSP70 as a genetic marker for thermo-tolerance and immuno-modulation in animals under climate change scenario. Animal. Nutrition. (2019) 5:340–50. doi: 10.1016/j.aninu.2019.06.005
39. Wang, J, Que, SQ, Liu, X, Jin, M, Xin, TR, Zou, ZW, et al. Characteristic and expression of Hsp70 and Hsp90 genes from Tyrophagus putrescentiae and their response to thermal stress. Sci Rep [Internet]. (2021) 11:1–12. Available from:. doi: 10.1038/s41598-021-91206-2
40. Amaral, A. Energy metabolism in mammalian sperm motility. WIREs Mechanisms of Disease. (2022) 14:1569. doi: 10.1002/wsbm.1569
41. Castellini, C, D’andrea, S, Cordeschi, G, Totaro, M, Parisi, A, Di Emidio, G, et al. Pathophysiology of mitochondrial dysfunction in human spermatozoa: Focus on energetic metabolism, oxidative stress and apoptosis. Antioxidants MDPI. (2021) 10. doi: 10.3390/antiox10050695
42. Hitit, M, Özbek, M, Ayaz-Guner, S, Guner, H, Oztug, M, Bodu, M, et al. Proteomic fertility markers in ram sperm. Anim Reprod Sci. (2021) 235:106882. doi: 10.1016/j.anireprosci.2021.106882
43. Davila, MP, Muñoz, PM, Bolaños, JMG, Stout, TAE, Gadella, BM, Tapia, JA, et al. Mitochondrial ATP is required for the maintenance of membrane integrity in stallion spermatozoa, whereas motility requires both glycolysis and oxidative phosphorylation. Reproduction. (2016) 152:683–94. doi: 10.1530/REP-16-0409
44. Westhoff, D, and Kamp, G. Glyceraldehyde 3-phosphate dehydrogenase is bound to the fibrous sheath of mammalian spermatozoa. J Cell Sci. (1997) 110:1821–9. doi: 10.1242/jcs.110.15.1821
45. Welch, JE, Brown, PL, O'Brien, DA, Magyar, PL, Bunch, DO, Mori, C, et al. Human glyceraldehyde 3-phosphate dehydrogenase-2 gene is expressed specifically in spermatogenic cells. J Androl. (2000) 21:328–38. doi: 10.1002/j.1939-4640.2000.tb02111.x
46. Zini, A, O’Bryan, MK, Israel, L, and Schlegel, PN. Human sperm NADH and NADPH diaphorase cytochemistry: correlation with sperm motility. Urology. (1998) 51:464–8. doi: 10.1016/S0090-4295(97)00631-6
47. Elkina, YL, Atroshchenko, MM, Bragina, EE, Muronetz, VI, and Schmalhausen, EV. Oxidation of glyceraldehyde-3-phosphate dehydrogenase decreases sperm motility. Biochemistry (Mosc). (2011) 76:268–72. doi: 10.1134/S0006297911020143
48. Kumar, V, Rangaraj, N, and Shivaji, S. Activity of pyruvate dehydrogenase a (PDHA) in Hamster spermatozoa correlates positively with Hyperactivation and is associated with sperm capacitation. Biol Reprod. (2006) 75:767–77. doi: 10.1095/biolreprod.106.053587
Keywords: fertility, food security, LC–MS/MS, Madura bulls, sperm proteins fertility, proteomic, sperm proteins
Citation: Rosyada ZNA, Pardede BP, Kaiin EM, Gunawan M, Maulana T, Said S, Tumbelaka LITA, Solihin DD, Ulum MF and Purwantara B (2023) A proteomic approach to identifying spermatozoa proteins in Indonesian native Madura bulls. Front. Vet. Sci. 10:1287676. doi: 10.3389/fvets.2023.1287676
Edited by:
Lucas Lima Verardo, Universidade Federal dos Vales do Jequitinhonha e Mucuri (UFVJM), BrazilReviewed by:
Mahak Singh, The ICAR Research Complex for North Eastern Hill Region (ICAR RC NEH), IndiaMustafa Hitit, Kastamonu University, Türkiye
Copyright © 2023 Rosyada, Pardede, Kaiin, Gunawan, Maulana, Said, Tumbelaka, Solihin, Ulum and Purwantara. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bambang Purwantara, purwantara@apps.ipb.ac.id; Zulfi Nur Amrina Rosyada, nur.amrina@fkh.unair.ac.id