Figures
Abstract
Metabolites play essential roles in biological systems, but detailed identities and significance of the seminal plasma metabolome related to bull fertility are still unknown. The objectives of this study were to determine the comprehensive metabolome of seminal plasma from Holstein bulls and to ascertain the potential of metabolites as biomarkers of bull fertility. The seminal plasma metabolome from 16 Holstein bulls with two fertility rates were determined by gas chromatography-mass spectrometry (GC-MS). Multivariate and univariate analyses of the data were performed, and the pathways associated with the seminal plasma metabolome were identified using bioinformatics approaches. Sixty-three metabolites were identified in the seminal plasma of all bulls. Fructose was the most abundant metabolite in the seminal fluid, followed for citric acid, lactic acid, urea and phosphoric acid. Androstenedione, 4-ketoglucose, D-xylofuranose, 2-oxoglutaric acid and erythronic acid represented the least predominant metabolites. Partial-Least Squares Discriminant Analysis (PLSDA) revealed a distinct separation between high and low fertility bulls. The metabolites with the greatest Variable Importance in Projection score (VIP > 2) were 2-oxoglutaric acid and fructose. Heat-map analysis, based on VIP score, and univariate analysis indicated that 2-oxoglutaric acid was less (P= 0.02); whereas fructose was greater (P= 0.02) in high fertility than in low fertility bulls. The current study is the first to describe the metabolome of bull seminal plasma using GC-MS and presented metabolites such as 2-oxoglutaric acid and fructose as potential biomarkers of bull fertility.
Citation:Velho ALC, Menezes E, Dinh T, Kaya A, Topper E, Moura AA, et al. (2018) Metabolomic markers of fertility in bull seminal plasma. PLoS ONE 13(4): e0195279. https://doi.org/10.1371/journal.pone.0195279
Editor:Carlos E. Ambrósio, Faculty of Animal Sciences and Food Engineering, University of São Paulo, BRAZIL
Received:December 10, 2017;Accepted:March 19, 2018;Published:April 10, 2018
Copyright:© 2018 Velho et al. This is an open access article distributed under the terms of theCreative Commons Attribution License,which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability:All relevant data are within the paper and its Supporting Information files.
Funding:This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2017-67016-26507 from the USDA National Institute of Food and Agriculture (NIFA). Partial funding was also provided by Mississippi Agricultural Experiment Station, Alta Genetics Inc., Hatch project under accession number 1005778 from the USDA NIFA, and by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) of Brazil. A. Velho and E. Menezes were funded by competitive fellowships from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) of Brazil.
Competing interests:The authors have declared that no competing interests exist.
Introduction
Male fertility relates to the capacity of an animal to produce spermatozoa with the ability to fertilize the oocyte, resulting in a living offspring. Fertility is affected by several factors, including management, nutrition, disease, stress, age, and genetics [1]. A decline in bull fertility affects the conception rate of herds, resulting in decreased production and profit. Therefore, the ability to predict bull fertility in advance offers enormous benefits for the economic success of livestock enterprise by improving pregnancy rates [2]. The “omics”approaches, such as genomics, transcriptomics and proteomics, have been used to ascertain molecular determinants of bull fertility. As a result, studies show that molecular compounds found in sperm and seminal plasma are significantly associated with bull fertility [3–6].
Seminal plasma is a complex mixture of secretions from testis, epididymis and accessory sex glands. Given the molecular contributions of the seminal plasma in sperm physiology, metabolites may affect downstream and complementary changes in gene and protein expressions and may be the underlying the key regulators of bull fertility [7]. Metabolites present in seminal plasma play several roles related to sperm function, such as energy production, motility, protection, pH control and regulation of metabolic activity [8].
Metabolomics is an emerging technique and has shown promise in identifying potential male fertility and infertility biomarkers [8–11]. Metabolomics represents the downstream of systems biology and has drawn significant interest for studying and understand fundamental biological processes related to reproduction [10,12] since it allows the identification and quantification of small molecules, such as amino acids, peptides, fatty acid, and carbohydrates in secretions, cells, tissues, and organs [13–15]. This then reveals information about metabolic reactions and mechanisms that can help identify potential biomarkers of phenotypes of interest [14–16]. There have been some studies on the potential functions of metabolites in sperm physiology and their specific roles in metabolic pathways for reproductive success. For example, metabolic profiling analysis of mouse spermatozoa showed that the action of glycolytic substrates was associated with tyrosine phosphorylation and energy production, which is essential for flagella motility [17]. In the boar species, a study compared the pathways of glycolysis and gluconeogenesis to produce lactate/pyruvate and observed that the route of gluconeogenesis was limited use to produce energy for the spermatozoa [18]. Recently, using NMR and GC-MS techniques to determine metabolites in human spermatozoa, a total of 42 metabolites were revealed [19]. Also, supplementation of metabolites into extender medium used for liquid storage of goat sperm showed that glucose and pyruvate supplementation were better as compared to lactate for maintenance of sperm motility [20].
Metabolomics approaches have been also used to identify potential fertility biomarkers in the seminal plasma of bulls [11] and men [10,20–26]. Moreover, a recent study reported that gas chromatography-mass spectrometry (GC-MS) can be an innovative method for fast and noninvasive diagnostic of man infertility [27]. Analyses of human seminal plasma by Raman spectroscopy allowed the identification of metabolites associated with asthenozoospermic patients as well [21]. In another study, GC-MS of seminal fluid detected lower concentrations of palmitic acid and oleic acid in healthy men as compared to asthenozoospermic individuals [26]. These findings suggest that metabolomic analyses are valid tools for identification of many classes of molecules associated with metabolic pathways essential to reproductive success [10,12].
Although metabolomic tools have been used to study man’s fertility [20–27] there are only a limited number of such studies in farm animals. Vast gaps in the knowledge base exist, including metabolite identities, their concentrations in seminal plasma of bulls, as well as the molecular mechanism of their involvement in fertility. As such, the present study was conducted to perform a comprehensive analysis of the seminal plasma metabolome from adult Holstein bulls using GC-MS. In addition, we tested the hypothesis that differences exist in the seminal fluid metabolome between sires of contrastingin vivofertility scores.
Materials and methods
Experimental design
Comprehensive metabolomics analysis of seminal plasma from Holstein bulls (n = 16) with contrastingin vivofertility categories was performed using GC-MS. Following the analysis of metabolome data, computational biology tools were employed to detect potential biomarkers for bulls of high (n = 8) and low (n = 8) fertility.
Sample collection and determination of bull fertility
Seminal plasma samples from 16 Holstein bulls with contrasting fertility phenotypes were provided by Alta Genetics (Watertown, WI, USA). All animals were raised under the same management conditions and received the same nutrition. Semen was collected with artificial vagina and seminal plasma was separated from sperm by centrifugation (700 ×g,4°C, 10 min). Supernatant seminal plasma was then transferred to a 2-mL microcentrifuge tube and centrifuged again (10,000 ×g,60 min, 4°C), as described previously [28]. After the second centrifugation, seminal plasma was aliquoted (100 μl) into a 2-mL Cryotube®(Sigma-Aldrich, St Louis, USA), snap-frozen in liquid nitrogen, and transported to Mississippi State University (MSU). At MSU, seminal plasma was stored at -80°C until preparation for GC-MS.
In the present study, every ejaculate collected from each bull was routinely evaluated by standard semen analysis methods [29]. The method to determine bull fertility used in our study is also similar to previous investigations about fertility biomarkers in bulls, conducted by several authors in the last decades [4,29–33]. The calculation of fertility scores was based on the actual conception rates confirmed by either veterinary palpation or ultrasound of cows inseminated with hundreds/thousands frozen-thawed semen straws from each bull. Sires used for this study were selected based on their fertility scores (Table 1), as previously described by Peddintiet al.[32]. Factors that influenced fertility performance of sires, e.g., breeding event, environmental factors and herd management were adjusted to determine reliable fertility scores using threshold models [34,35]. Using Probit.F90 software [36], fertility prediction of each sire was calculated according to the average conception of more than 300 breeding outcomes along with their percent deviation of conception rates. Based on this calculation, percent deviation of conception rates was used to categorize the fertility of bulls.
For this marker discovery study, bulls were selected among the ones performing conception rate of two standard deviation above and below the average of the bull population in the database. Thus, bulls which had percent difference of their conception rate above average were defined as high fertility (HF) with average of 3,174 breeding outcomes. Those bulls that had percent difference of their conception rate below average were classified as low fertility (LF) having 1,662 breeding outcomes (Table 1).
Sample preparation for gas chromatography-mass spectrometry analysis
Preparation of seminal plasma samples for GC-MS analysis was performed as described by Shiet al.[37], with modifications. Briefly, a 100 μL aliquot of samples or reference standards (10 μM; Sigma-Aldrich, St Louis, USA) was thawed and pipetted into a 2-mL polypropylene microcentrifuge tube (Sigma-Aldrich, St Louis, USA). Then, 150 μL of heptadecanoic acid in methanol (1 mg/mL, internal standard; Sigma-Aldrich, St Louis, USA) and 350 μL of methanol were then added to the tube. This mixture was vortexed vigorously for 1 min and was then centrifuged at 12,000 ×gat 4°C for 10 min. A volume of 100 μL of supernatant was transferred to a 2-mL amber glass vial (Agilent Technologies, Santa Clara, CA) and the solvent was evaporated to dryness in a TurboVap®LV evaporator (Biotage, Charlotte, NC) with a gentle stream of nitrogen at 45°C.The dried extract was suspended in 50 μL of methoxyamine hydrochloride in pyridine (20 mg/mL; Sigma-Aldrich, St Louis, USA), vortexed vigorously for 1 min, and heated in a water bath at 70°C for 1 h. The sample was then derivatized by adding 100 μL of N,O-Bis(trimethylsilyl)trifluoroacetamide with 1% trimethylchlorosilane (BSTFA + 1% TMCS; Sigma-Aldrich, St Louis, USA) and heated again in a water bath at 70°C for 1 h. Derivatives of metabolites were transferred to an amber glass vial having a fixed insert (Agilent Technologies, Santa Clara, CA) for GC-MS analysis.
Gas chromatography-mass spectrometry analysis
Samples and reference standards were analyzed using an Agilent 7890A GC System coupled to an Agilent 5975C inert XL MSD with triple-axis mass detector, an Agilent 7693 Series Autosampler, and a DB-5MS capillary column (30 m × 0.25 mm i.d. × 0.25 μm film thickness; Agilent Technologies, Santa Clara, CA). A volume of 1 μL of derivatized mixture was injected into the inlet heated at 280°C with 1:10 split ratio. Standard septum purge was performed after sample injection at 3 mL/min and helium carrier gas was at 1 mL/min constant flow rate. Transfer line, ion source, and quadrupole were heated at 250°C, 230°C, and 150°C, respectively. Oven was programmed initially at 70°C for 4 min, ramped up to 300ºC at 8°C/min, and then held at 300ºC for 5 min. Ionization was performed in an electron impact mode at 70 eV. Masses were scanned for full spectra from m/z 35 to 800 at 10,000 amu/s and 10.3 scans/s (m/z 0.2 step size). The solvent delay time was at 6 min.
Calculation and statistical analysis
All identified metabolites were categorized according to their chemical classes. All compounds were identified by their retention times and one target and two quantitative ions in comparison with mass spectra of authentic standards and mass spectra in the NIST Mass Spectral Search Program (NIST/EPA/NIH Mass Spectral Library, Version 2.0). Abundances of the target ions of each metabolite were divided by that of target ion of the internal standard (heptadecanoic acid) and the unit less ratios were used for statistical analysis [37].
Multivariate analysis of abundance ratios was conducted by uploading data to MetaboAnalyst 3.0 (http://www.metaboanalyst.ca) [38]. The generated data matrix were normalized to a constant sum, auto-scaled, and analyzed by PLS-DA. VIP based on the PLS-DA model was calculated to identify the potential biomarkers and those variables with VIP score of more than 1.5 were considered important for group separation [11,38,39]. Univariate analyses were conducted using a simple t-test to evaluate the statistical significance between high fertility and low fertility bulls atP≤ 0.05. Correlations among all seminal plasma metabolites were determined using Pearson's method (P≤ 0.05) [38,39]. In addition, Pearson's correlation (P <0.05) was to determine the strength of the associations between seminal plasma metabolites (abundance ratios) and fertility scores (% deviation of conception rates).
In silicoanalysis of metabolic networks
Bioinformatics tool was used to visualize metabolic networks and pathways. Metabolic networks related to the most abundant metabolites in bull seminal plasma, as well as for two metabolites with the highest VIP scores, were analyzed using Metscape version 3.1.2, which is a plug-in for Cytoscape version 3.2.1 (http://www.cytoscape.org). The Kyoto Encyclopedia of Genes and Genomes (KEGG) was used for compound identification [40].
Results
Metabolome of bull seminal plasma
A total of 73 peaks were integrated after GC-MS analysis of bull seminal plasma, regardless of fertility scores. Of these peaks, 63 metabolites were identified and categorized according to their major chemical classes, including amino acids, peptides/analogues, carbohydrates/carbohydrate conjugates, fatty acids/conjugates, steroids/steroid derivatives, nucleosides/nucleotides/analogues, and other organic and inorganic compounds (Table 2). A characteristic GC-MS chromatogram of bull seminal plasma and peaks of important metabolites were depicted inFig 1.Of the 63 metabolites identified in the current study, 24 compounds were authenticated by external standard references and 39 compounds were identified by probable match parameters of the NIST Mass Spectral Search Program. In addition, retention time, target ion, and two quantitative ions, and chemical structure of the derivative product of each metabolite were also used for identification and quantification purposes (Table 2).
Peaks of lactic acid, glycerol, phosphoric acid, citric acid, fructose, heptadecanoic acid (internal standard) and cholesterol are indicated.
Metabolites chemical classes
Amino acids, peptides and analogues were the major compounds found in the bovine seminal plasma. Twenty-one amino acids were detected, including glutamic acid, alanine, isoleucine, leucine, and serine, among others. Twenty organic compounds, including citric acid, lactic acid, urea, uric acid, and myo-inositol comprised the second major group. Twelve carbohydrates and their conjugate metabolites were identified in the bull seminal plasma including fructose, mannitol, sorbitol, glycerol, and ribitol. In addition, we detected fatty acids and conjugates (malic acid, hexadecanoic acid, methylmaleic acid and stearic acid), steroids and steroid derivatives (androstenedione and cholesterol) and few nucleosides, nucleotides, and analogues (inosine and 5-methyluridine) and inorganic compounds (borate and phosphoric acid;Table 2andFig 2).
Metabolites identified were categorized according to their chemical classes, defined as amino acids, peptides and analogues; carbohydrates and carbohydrate conjugates; fatty acids and conjugates; steroids and steroid derivatives; nucleosides, nucleotides, and analogues; others organic compounds and inorganic compounds.
Most and least predominant metabolites in the bull seminal plasma
Fructose was the most predominant metabolite in seminal plasma of all bulls. Other predominant metabolites, based on their abundance ratios, were citric acid, lactic acid, urea and phosphoric acid (Fig 3A). Androstenedione, 4-ketoglucose, D-xylofuranose, 2-oxoglutaric acid and erythronic acid were among the five least predominant metabolites identified in bull seminal plasma (Fig 3B).
(A) The five most abundant metabolites were fructose, citric acid, lactic acid, urea and phosphoric acid and (B) the five least metabolites were identified as androstenedione, 4-ketoglucose, D-xylofuranose, 2-oxoglutaric acid and erythronic acid. Abundance ratio of the metabolites was calculated by dividing abundance of target ions of metabolites by that of target ion of the internal standard. Error bars represent standard error of the mean.
Associations between seminal plasma metabolites and bull fertility
A multivariate analysis (Partial-Least Squares Discriminant Analysis; PLS-DA) of the seminal plasma metabolome indicated a distinct separation between high (HF) and low fertility (LF) bulls, as shown by PLS-DA score plot (Fig 4). In addition, metabolites with Variable Importance in Projection (VIP) score greater than 1.5 were identified as 2-oxoglutaric acid, fructose, phosphoric acid, D-mannitol, 4-ketoglucose, dulcitol and erythronic acid (Fig 5). Among these compounds, 2-oxoglutaric acid had the highest VIP score (VIP = 2.17), followed by fructose (VIP = 2.1). VIP score and the corresponding heat-map indicate the high or low abundance ratio of each metabolite in HF and LF bulls. Abundance ratios of 2-oxoglutaric acid, phosphoric acid, D-mannitol and dulcitol were lower in HF than in LF bulls. However, fructose, 4-ketoglucose and erythronic acid were more abundant in high fertility in comparison with low fertility bulls.
The plots indicate that a separation could be observed between HF and LF bulls. Supervised PLS-DA was obtained with 2 components. The explained variances are shown in parentheses.
The selected metabolites were those with VIP score of greater than 1.5. Heat map with red or green boxes on the right indicates high and low abundance ratio, respectively, of the corresponding metabolite in HF and LF bulls. VIP score was based on the PLS-DA model.
Based on univariate analysis, abundance ratios of 2-oxoglutaric acid (P= 0.02), ornithine (P= 0.03), L-leucine (P= 0.04) and D-mannitol (P= 0.04) were lower in HF than in LF bulls, whereas abundance ratio of fructose was greater (P= 0.02) in HF as compared to LF bulls (Fig 6). Fructose was positively correlated with amino-butyrolactone (r= 0.55;P= 0.005), oxalic acid (r= 0.69;P< 0.0005) and fumaric acid (r= 0.57;P< 0.005), and inversely associated with p-tolyl-beta-D-glucuronide (r= -0.51;P= 0.01), lactic acid (r= -0.54;P= 0.006), urea (r= -0.52;P= 0.009), 2-pyrrolidone-5-carboxylic acid (r= -0.51;P= 0.01), glutamic acid (r= -0.53;P= 0.008), phosphoric acid (r= - 0.54;P= 0.007) and D-glycero-D-gulo-heptose (r= -0.56;P= 0.004) (Fig 7).
(A) 2-oxoglutaric acid, (B) fructose, (C) ornithine, (D) L-leucine and (E) D-mannitol were significantly different (P≤ 0.05) between high (HF) and low (LF) bulls.
Based on Pearson’s method, abundance ratio of ornithine (r = 0.59; P = 0.02) was significantly associated with fertility score of LF bulls. In addition, there was also a positive and significant correlation between abundance ratio of D-mannitol (r = 0.53; P = 0.04) and fertility score of LF bulls.
Metabolic networks of seminal plasma metabolites
Metabolic networks were determined for fructose, citric acid, lactic acid, and urea (the most abundant seminal plasma metabolites) as well as for 2-oxoglutaric acid (metabolite with the highest VIP score). Based on fructose network panel (Fig 8A), D-sorbitol undergoes a reversible reaction to produce fructose, which is converted to beta-D-fructose-6-phosphate and D-fructose-1-phosphate. Citric acid can be interconverted to cis-aconitate, isocitrate, acetate and oxaloacetate. Then, non-reversible reactions yield acetyl-CoA and oxaloacetate (Fig 8B). Lactic acid is synthesized from S-lactoylglutathione or pyruvate in a non-reversible and reversible reaction, respectively (Fig 8C). Urea can be produced from L-arginine and allantoate, which is converted to urea and ureidoglycolate. L-arginine is hydrolyzed to urea and ornithine (Fig 8D). Furthermore, 2-oxoglutaric acid is involved in several reversible and non-reversible reactions, as it can be synthetized from L-glutamate, succinyl-CoA, 2-hydroxyglutarate, 2-methyl-3-oxopropanoate and 3-amino-2-methylpropanoate, among others. It can be converted to 3-carboxy-1-hydroxypropyl-ThPP, succinate and peptide-3-hidroxy-L-aspartate (Fig 8E). Reactions and enzymes of each metabolite network are presented inS1 Table.
(A) fructose, (B) citric acid, (C) lactic acid, (D) urea, (E) 2-oxoglutaric acid. Networks were generated using the Metscape plug-in for Cytoscape. Metabolites are shown in red hexagons. Compounds are represented as pink hexagons, reversible reactions as gray squares with purple text and non-reversible reactions as gray squares with orange text.
Discussion
In the present study, we performed GC-MS analysis to determine the metabolite profile of bull seminal plasma and identify potential biomarkers of fertility. Additionally, we used bioinformatics tools to reveal the networks and reactions in which bull seminal plasma metabolites might be involved. A previous study has investigated amino acid and fatty acid composition of bovine seminal plasma using GC-MS [41]. To the best of our knowledge, however, our study is in fact the first to conduct a comprehensive evaluation of the bull seminal plasma metabolome, which includes not only amino acids and fatty acids but also carbohydrates, nucleosides, steroids, organic and inorganic compounds by GC-MS. Moreover, we described the association of specific seminal plasma metabolites with bull fertility scores measuredin vivo.
Metabolites, products of metabolic reactions, appear in numerous biochemical pathways [42] and have been reported as potential biomarkers for male fertility [10,20,25–27]. At ejaculation, spermatozoa are suspended in the seminal plasma and it marks clear qualitative and quantitative differences in its biochemical composition [43]. Although the effects of the seminal plasma are still not clearly understood, the exposure of spermatozoa to small molecules such as metabolites can improve or hinder sperm fertilizing capacity, even if semen is diluted in extender during preservation of spermatozoa. The removal of seminal plasma is not necessary for bull sperm preservation and the timeline of the preservation protocol (e.g., a long cooling curve and equilibrium step) leads to a prolonged exposure of spermatozoa to seminal plasma molecules [44]. In addition, the interaction of metabolites with other molecules in the uterine environment impacts fertilization, implantation, and optimal fetal and placental developments [5]. Furthermore, metabolites such as amino acids, peptides, carbohydrates, fatty acids, steroids and nucleosides, among others, participate in physiologically important events, influencing energy production, motility, pH control, membrane protection and metabolic activity of the sperm [22,45–47].
As presently evaluated by GC-MS, the major chemical classes of metabolites in the bovine seminal plasma were defined as amino acids, peptides, and their analogues, followed by carbohydrates and carbohydrate conjugates. We identified 21 metabolites classified as amino acids, peptides, and their analogues in bull seminal plasma. Similarly, researchers have detected 20 [48] to 23 amino acids [41] in seminal plasma of bulls, using GC-MS. A large number of amino acids were also found in human sperm [19] and in goat epididymal fluid [47]. Besides playing essential roles as basic building block of proteins [49], amino acids participate in crucial steps of sperm biology, including protection and regulation of metabolic activity [47]. Amino acids protect ram sperm cells during cryopreservation by decreasing lipid peroxidation and injury caused by free-radicals [50]. Meanwhile, the presence of carbohydrates in mammalian seminal plasma is essential for sperm because these molecules are part of crucial energy production pathways [46], and glycolysis is used by mammalian sperm to obtain energy. During this event, seminal plasma glycolysable carbohydrates, such as fructose, are required for ATP production, which leads to increased respiratory activity to support optimum sperm motility and survival [51,52].
The most predominant metabolites of the bull seminal plasma were fructose, citric acid, lactic acid, urea and phosphoric acid; whereas androstenedione, 4-ketoglucose, D-xylofuranose, 2-oxoglutaric acid and erythronic acid were among the least abundant. Multivariate statistical analysis showed that fructose abundance separated HF from LF groups and it had the second highest VIP score. Fructose abundance ratio was higher (P= 0.02) in HF than in LF bulls. Fructose is the primary energy source for bull spermatozoa and the major carbohydrate in seminal plasma of these animals [52–54]. Fructose is produced in the seminal vesicles, under androgen stimulation [46], and it has been found in seminal plasma of several species, including buffalo [55], goat [55,56], ram [57], boar [58], human [59], and rabbit [60]. As revealed byin siliconetwork analysis, fructose is involved in fundamental pathways for energy production for the spermatozoa. Fructose can be produced from sorbitol by the action of sorbitol dehydrogenase, which is located on the sperm membrane. Once fructose enters the cell, it is converted to beta-D-fructose-6-phosphate and fructose-1-phosphate [53,61]. Furthermore, fructose can be metabolized to lactic acid, depending on various factors such as pH and temperature [52,53]. Smaller fructose concentration has been found in the seminal plasma of azoopermic, oligozoospermic, and idiopathic infertile man compared to fertile man [59]. Yousef et al. [60] suggested that a decrease in fructose concentration observed in the seminal plasma of rabbits intoxicated with aluminum chloride could be one of the factors reducing sperm motility. Therefore, less fructose abundance in bull seminal plasma reduces energy supply to sperms, negatively affecting their metabolism and, subsequently, male fertility.
As in bull seminal plasma, citric acid is also found in semen of other species, such as boar [62], human [26,59], and rabbit [63]. Citric acid helps control pH in boar semen and acts as a chelator for zinc, magnesium, and calcium [62]. The concentration of zinc, magnesium and calcium in human seminal plasma and their chelation can influence sperm metabolism, affecting sperm transport, acrosome reaction, and fertilization [64]. A recent study demonstrated that citric acid in seminal plasma is associated with bull fertility by potentially affecting sperm capacitation and acrosome reaction [11]. Lactic acid, like fructose and citric acid, is used as another important energy source for boar sperm [65] and lesser concentration of lactic acid has been detected in bull sperm with low viability [66]. Both urea and phosphoric acid were detected in the seminal plasma samples of our study. Urea has also been found in human seminal fluid [20,26,67], however, the main role it plays in semen is still unknown. Given that urea is an end product of protein metabolism [68] and that seminal fluid contains protein, it is possible that seminal plasma urea is the consequence of protein degradation. We also observed that phosphoric acid had the third highest VIP score, meaning that this metabolite is associated with lower fertility scores of the Holstein bulls. In spermatozoa, phosphoric acid can be a product of a reaction catalyzed by inorganic pyrophosphatase (PPA1). PPA1 catalyzes the hydrolysis of one molecule of inorganic pyrophosphate (PPi) to two molecules of phosphoric acid, leading to the release of energy. The transport of PPi from spermatozoa to the seminal plasma may be regulated by a transmembrane protein, called progressive ankylosis protein (ANKH). Therefore, the energy produced from the conversion of PPi to phosphoric acid could be utilized for sperm motility and during fertilization [69]. Moreover, inorganic phosphate can be resulted from the hydrolysis of ATP used for metabolism and sperm motility. Presence of greater abundance of inorganic phosphate indicated an increase in level of either motility or metabolism during ejaculation. At this point, it is unclear why greater inorganic phosphate was associated with LF group. We hypothesized that increased activity during ejaculation led to decreased reservation of energy, thereby decreasing fertility.
In the current study, propylene glycol was identified in bull seminal plasma. This finding is in agreement with previous studies showing that propylene glycol is present in human sperm [19], serum [70] and urine [71]. In addition, propylene glycol, a synthetic molecule, is found in pharmaceutical products and in food that are routinely used for animal feeding. As a matter of fact, propylene glycol has been used as feed additive for dairy cattle [72]. Therefore, we believe that the presence of propylene glycol in bull seminal plasma is not the result of contamination since this molecule has also been detected in human cells and secretions. However, it is possible that propylene glycol originally comes from feeds and it may also be the reason it was found in human and animal samples. The role of this compound in sperm fertility is currently unknown.
The compound identified as 2-oxoglutaric acid was one of the least abundant metabolites in bull seminal plasma but had the greatest VIP score, being less abundant in high fertility than in low fertility bulls. In agreement with our results, 2-oxoglutaric acid was found in greater concentration in the seminal plasma of men with asthenozoospermia than in that of the healthy ones [73]. It is known that 2-oxoglutaric acid can be synthesized from glutamate by 2-oxoglutarate aminotransferase [74], as also observed by the in silico analysis conducted in the present study. In boar semen, 2-oxoglutarate aminotransferase was mainly detected in the spermatozoa [75] and 2-oxoglutarate aminotransferase was released at high levels from boar sperm with low viability after semen freezing [76]. Other experiments proposed that high activity of 2-oxoglutarate aminotransferase in seminal plasma is an indicator of sperm membrane injury and poor semen quality in goats and bulls [77–79]. Therefore, all these findings agree with our work reporting the occurrence of higher abundance ratio of 2-oxoglutaric acid content in seminal plasma of low fertility bulls.
In recent decades,omicsapproaches have been used to identify several seminal plasma macromolecules associated with male fertility [29,30,80,81]. Using 2-D gels and mass spectrometry, an interesting study on the investigation of capacitation-related proteins in boar spermatozoa identified, afterin vitrocapacitation, a significant increase in tyrosine-phosphorylated proteins. The results of this study revealed differentially expressed proteins involved in cellular processes and biological regulation. This finding is important not only as related to semen physiology, but also to determinein vivofertility in the boar species [82]. Another proteomic work was conducted to evaluate differences in protein profiles between high- and low-litter sizes in boar spermatozoa in a field trial [83]. The study uncovered proteins potentially involved in the regulation of male fertility by contributing in capacitation, acrosome reaction, sperm-egg interaction, and fertilization processes. These results provide insights into the role played by these proteins in sperm physiology and male fertility. Fertility-related proteins in capacitated boar spermatozoa were also reported by Kwonet al.[84]. The study showed proteins differentially expressed between high- and low-litter size spermatozoa (> 3-folds) involved in sperm physiology. In addition, ras-related protein Rab-2A and cytochrome b-c1 complex subunit 1 expression negatively associated with litter size spermatozoa whereas cytochrome b-c1 complex subunit 2 was positively correlated with litter size spermatozoa. These results suggest that those proteins found in sperm following capacitation could be helpful as indicators of male fertility [84]. Proteins present in seminal plasma along with spermatozoa may contribute to improvement of fertilizing capacity of the sperm. In boar seminal plasma, fibronectin-1, one of the most abundant proteins in boar seminal plasma, was identified [85]. The presence of fibronectin-1 and its association with sperm function demonstrates that seminal plasma has potentially relevant biomarkers for cryoinjury [86,87]. Using bull epididymal spermatozoa as a model, a comprehensive proteomic study employed 2D electrophoresis technique and identified proteins associated with cryostress and their associated signaling pathways, including the ephrinR-actin pathway, the ROS metabolism pathway, actin cytoskeleton assembly, actin cytoskeleton regulation, and the guanylate cyclase pathway [88]. Another study showed evidences that the addition of cryoprotectant alters the bull epididymal sperm proteome [89]. The authors also observed that NADH dehydrogenase flavoprotein 2, f-actin-capping protein subunit beta, superoxide dismutase 2, and outer dense fiber protein 2 were associated with several important signaling pathways. These findings are useful to understand the mechanisms of cryoprotectant in protecting bull spermatozoa and for the development of novel extenders. Metabolomics technology emerged as a means of understanding physiological and pathological conditions in a large-scale manner, through the identification of metabolic substrates and products of a given biochemical system. Considering the metabolome as the metabolic state of a given physiologic status of a given fluid, cell, tissue, or organism, metabolomics is not only a complementary tool for understanding proteomics data, but also to gain insight into biochemical reaction networks, to understand mechanistically how the metabolites may affect male fertility as well as potential biomarker discovery.
In conclusion, we have demonstrated that fructose, citric acid, lactic acid, urea and phosphoric acid are the predominant metabolites in bull seminal plasma. We discovered a clear separation of metabolite profiles between high and low fertility bulls, with fructose and 2-oxoglutaric acid being potential candidates for biomarkers of bull fertility. Findings of the present study will help advance our current understanding of the multifactorial and complex processes related to the physiology of male fertility.
Supporting information
S1 Table.Reactions and enzymes of metabolite networks identified in bull seminal plasma.
https://doi.org/10.1371/journal.pone.0195279.s001
(DOCX)
Acknowledgments
This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2017-67016-26507 from the USDA National Institute of Food and Agriculture (NIFA). Partial funding was provided by Mississippi Agricultural Experiment Station, Alta Genetics Inc., Hatch project under accession number 1005778 from the USDA NIFA, and by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) of Brazil. A. Velho and E. Menezes were funded by competitive fellowships from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) of Brazil. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
- 1. Petherick JC. A review of some factors affecting the expression of libido in beef cattle, and individual bull and herd fertility. Applied Animal Behaviour Science. 2005;90(3–4):185–205.
- 2. Moce E, Graham JK. In vitro evaluation of sperm quality. Anim Reprod Sci. 2008;105(1–2):104–18. pmid:18178345.
- 3. Feugang JM, Rodriguez-Osorio N, Kaya A, Wang H, Page G, Ostermeier GC, et al. Transcriptome analysis of bull spermatozoa: implications for male fertility. Reprod Biomed Online. 2010;21(3):312–24. pmid:20638337.
- 4. de Oliveira RV, Dogan S, Belser LE, Kaya A, Topper E, Moura A, et al. Molecular morphology and function of bull spermatozoa linked to histones and associated with fertility. Reproduction. 2013;146(3):263–72. pmid:23904564.
- 5. Bromfield JJ. A role for seminal plasma in modulating pregnancy outcomes in domestic species. Reproduction. 2016;152(6):R223–R32. pmid:27601714.
- 6. Westfalewicz B, Dietrich MA, Mostek A, Partyka A, Bielas W, Nizanski W, et al. Analysis of bull (Bos taurus) seminal vesicle fluid proteome in relation to seminal plasma proteome. J Dairy Sci. 2017;100(3):2282–98. pmid:28041731.
- 7. Fernandez-Gago R, Dominguez JC, Martinez-Pastor F. Seminal plasma applied post-thawing affects boar sperm physiology: a flow cytometry study. Theriogenology. 2013;80(4):400–10. pmid:23756043.
- 8. Bieniek JM, Drabovich AP, Lo KC. Seminal biomarkers for the evaluation of male infertility. Asian J Androl. 2016;18(3):426–33. pmid:26975492; PubMed Central PMCID: PMCPMC4854096.
- 9. Goodacre R, Vaidyanathan S, Dunn WB, Harrigan GG, Kell DB. Metabolomics by numbers: acquiring and understanding global metabolite data. Trends Biotechnol. 2004;22(5):245–52. pmid:15109811.
- 10. Kovac JR, Pastuszak AW, Lamb DJ. The use of genomics, proteomics, and metabolomics in identifying biomarkers of male infertility. Fertil Steril. 2013;99(4):998–1007. pmid:23415969; PubMed Central PMCID: PMCPMC3652543.
- 11. Kumar A, Kroetsch T, Blondin P, Anzar M. Fertility-associated metabolites in bull seminal plasma and blood serum: 1H nuclear magnetic resonance analysis. Mol Reprod Dev. 2015;82(2):123–31. pmid:25640164.
- 12. Muñoz M, Uyar A, Correia E, Díez C, Fernandez-Gonzalez A, Caamaño JN, et al. Non-invasive assessment of embryonic sex in cattle by metabolic fingerprinting of in vitro culture medium. Metabolomics. 2014;10(3):443–51.
- 13. Oliver SG, Winson MK, Kell DB, Baganz F. Systematic functional analysis of the yeast genome. Trends Biotechnol. 1998;16(9):373–8. pmid:9744112.
- 14. Fiehn O. Metabolomics—the link between genotypes and phenotypes. Plant Mol Biol. 2002;48(1–2):155–71. pmid:11860207.
- 15. Dunn WB, Bailey NJ, Johnson HE. Measuring the metabolome: current analytical technologies. Analyst. 2005;130(5):606–25. pmid:15852128.
- 16. Patti GJ, Yanes O, Siuzdak G. Innovation: Metabolomics: the apogee of the omics trilogy. Nat Rev Mol Cell Biol. 2012;13(4):263–9. pmid:22436749; PubMed Central PMCID: PMCPMC3682684.
- 17. Goodson SG, Qiu Y, Sutton KA, Xie G, Jia W, O'Brien DA. Metabolic substrates exhibit differential effects on functional parameters of mouse sperm capacitation. Biol Reprod. 2012;87(3):75. pmid:22837480; PubMed Central PMCID: PMCPMC3464911.
- 18. Marin S, Chiang K, Bassilian S, Lee WN, Boros LG, Fernandez-Novell JM, et al. Metabolic strategy of boar spermatozoa revealed by a metabolomic characterization. FEBS Lett. 2003;554(3):342–6. pmid:14623091.
- 19. Paiva C, Amaral A, Rodriguez M, Canyellas N, Correig X, Ballesca JL, et al. Identification of endogenous metabolites in human sperm cells using proton nuclear magnetic resonance ((1) H-NMR) spectroscopy and gas chromatography-mass spectrometry (GC-MS). Andrology. 2015;3(3):496–505. pmid:25854681.
- 20. Qiao S, Wu W, Chen M, Tang Q, Xia Y, Jia W, et al. Seminal plasma metabolomics approach for the diagnosis of unexplained male infertility. PLoS One. 2017;12(8):e0181115. pmid:28797078; PubMed Central PMCID: PMCPMC5552325.
- 21. Gilany K, Moazeni-Pourasil RS, Jafarzadeh N, Savadi-Shiraz E. Metabolomics fingerprinting of the human seminal plasma of asthenozoospermic patients. Mol Reprod Dev. 2014;81(1):84–6. pmid:24254375.
- 22. Deepinder F, Chowdary HT, Agarwal A. Role of metabolomic analysis of biomarkers in the management of male infertility. Expert Rev Mol Diagn. 2007;7(4):351–8. pmid:17620044.
- 23. Hamamah S, Seguin F, Bujan L, Barthelemy C, Mieusset R, Lansac J. Quantification by magnetic resonance spectroscopy of metabolites in seminal plasma able to differentiate different forms of azoospermia. Hum Reprod. 1998;13(1):132–5. pmid:9512244.
- 24. Hamamah S, Seguin F, Barthelemy C, Akoka S, Le Pape A, Lansac J, et al. 1H nuclear magnetic resonance studies of seminal plasma from fertile and infertile men. J Reprod Fertil. 1993;97(1):51–5. pmid:8464025.
- 25. Gilany K, Mani-Varnosfaderani A, Minai-Tehrani A, Mirzajani F, Ghassempour A, Sadeghi MR, et al. Untargeted metabolomic profiling of seminal plasma in nonobstructive azoospermia men: A noninvasive detection of spermatogenesis. Biomed Chromatogr. 2017;31(8). pmid:28058728.
- 26. Tang B, Shang X, Qi H, Li J, Ma B, An G, et al. Metabonomic analysis of fatty acids in seminal plasma between healthy and asthenozoospermic men based on gas chromatography mass spectrometry. Andrologia. 2017;49(9). pmid:28124472.
- 27. Zhou X, Wang Y, Yun Y, Xia Z, Lu H, Luo J, et al. A potential tool for diagnosis of male infertility: Plasma metabolomics based on GC-MS. Talanta. 2016;147:82–9. pmid:26592580.
- 28. Rego JP, 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–2):1–17. pmid:24889044.
- 29. Killian GJ, Chapman DA, Rogowski LA. Fertility-associated proteins in Holstein bull seminal plasma. Biol Reprod. 1993;49(6):1202–7. pmid:8286602.
- 30. Moura AA, Chapman DA, Koc H, Killian GJ. Proteins of the cauda epididymal fluid associated with fertility of mature dairy bulls. J Androl. 2006;27(4):534–41. pmid:16582409.
- 31. Moura AA, Chapman DA, Killian GJ. Proteins of the accessory sex glands associated with the oocyte-penetrating capacity of cauda epididymal sperm from holstein bulls of documented fertility. Mol Reprod Dev. 2007;74(2):214–22. pmid:16941673.
- 32. Peddinti D, Nanduri B, Kaya A, Feugang JM, Burgess SC, Memili E. Comprehensive proteomic analysis of bovine spermatozoa of varying fertility rates and identification of biomarkers associated with fertility. BMC Syst Biol. 2008;2:19. pmid:18294385; PubMed Central PMCID: PMCPMC2291030.
- 33. Han Y, Penagaricano F. Unravelling the genomic architecture of bull fertility in Holstein cattle. BMC Genet. 2016;17(1):143. pmid:27842509; PubMed Central PMCID: PMCPMC5109745.
- 34. Zwald NR, Weigel KA, Chang YM, Welper RD, Clay JS. Genetic selection for health traits using producer-recorded data. I. Incidence rates, heritability estimates, and sire breeding values. J Dairy Sci. 2004;87(12):4287–94. pmid:15545392.
- 35. Zwald NR, Weigel KA, Chang YM, Welper RD, Clay JS. Genetic selection for health traits using producer-recorded data. II. Genetic correlations, disease probabilities, and relationships with existing traits. J Dairy Sci. 2004;87(12):4295–302. pmid:15545393.
- 36. Chang YM, Gianola D, Heringstad B, Klemetsdal G. Effects of trait definition on genetic parameter estimates and sire evaluation for clinical mastitis with threshold models. Animal Science. 2004;2004(79):355–63.
- 37. Shi S, Yi L, Yun Y, Zhang X, Liang Y. A combination of GC-MS and chemometrics reveals metabolic differences between serum and plasma. Analytical Methods. 2015;7(5):1751–7.
- 38. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acids Res. 2015;43(W1):W251–7. pmid:25897128; PubMed Central PMCID: PMCPMC4489235.
- 39. Xia J, Wishart DS. Web-based inference of biological patterns, functions and pathways from metabolomic data using MetaboAnalyst. Nat Protoc. 2011;6(6):743–60. pmid:21637195.
- 40. Karnovsky A, Weymouth T, Hull T, Tarcea VG, Scardoni G, Laudanna C, et al. Metscape 2 bioinformatics tool for the analysis and visualization of metabolomics and gene expression data. Bioinformatics. 2012;28(3):373–80. pmid:22135418; PubMed Central PMCID: PMCPMC3268237.
- 41. Holden SA, Fernandez-Fuertes B, Murphy C, Whelan H, O'Gorman A, Brennan L, et al. Relationship between in vitro sperm functional assessments, seminal plasma composition, and field fertility after AI with either non-sorted or sex-sorted bull semen. Theriogenology. 2017;87:221–8. pmid:27678515.
- 42. Dunn WB, Broadhurst DI, Atherton HJ, Goodacre R, Griffin JL. Systems level studies of mammalian metabolomes: the roles of mass spectrometry and nuclear magnetic resonance spectroscopy. Chem Soc Rev. 2011;40(1):387–426. pmid:20717559.
- 43. Einarsson S. Studies on the composition of epididymal content and semen in the boar. Acta Vet Scand Suppl. 1971;36:1–80. pmid:5290046.
- 44. Forero-Gonzalez RA, Celeghini EC, Raphael CF, Andrade AF, Bressan FF, Arruda RP. Effects of bovine sperm cryopreservation using different freezing techniques and cryoprotective agents on plasma, acrosomal and mitochondrial membranes. Andrologia. 2012;44 Suppl 1:154–9. pmid:22506813.
- 45. Williams AC, Ford WC. The role of glucose in supporting motility and capacitation in human spermatozoa. J Androl. 2001;22(4):680–95. pmid:11451366.
- 46. Juyena NS, Stelletta C. Seminal plasma: an essential attribute to spermatozoa. J Androl. 2012;33(4):536–51. pmid:22016346.
- 47. Patel AB, Srivastava S, Phadke RS, Govil G. Identification of low-molecular-weight compounds in goat epididymis using multinuclear nuclear magnetic resonance. Anal Biochem. 1999;266(2):205–15. pmid:9888977.
- 48. al-Hakim MK, Graham EF, Schmehl ML. Free amino acids and amino compounds in bovine seminal plasma. J Dairy Sci. 1970;53(1):84–8. pmid:5413074.
- 49. Papp G, Grof J, Menyhart J. The role of basic amino acids of the seminal plasma in fertility. Int Urol Nephrol. 1983;15(2):195–203. pmid:6629696.
- 50. Sangeeta S, Arangasamy A, Kulkarni S, Selvaraju S. Role of amino acids as additives on sperm motility, plasma membrane integrity and lipid peroxidation levels at pre-freeze and post-thawed ram semen. Anim Reprod Sci. 2015;161:82–8. pmid:26362050.
- 51. Ford WC. Glycolysis and sperm motility: does a spoonful of sugar help the flagellum go round? Hum Reprod Update. 2006;12(3):269–74. pmid:16407453.
- 52. Mann T. Studies on the metabolism of semen: 3. Fructose as a normal constituent of seminal plasma. Site of formation and function of fructose in semen. Biochem J. 1946;40(4):481–91. pmid:16748040; PubMed Central PMCID: PMCPMC1258384.
- 53. TE K T. M. Sorbitol metabolism in spermatozoa. Proc R Soc Lond B Biol Sci. 1959;151:226–46.
- 54. Liberda J, Kraus M, Ryslava H, Vlasakova M, Jonakova V, Ticha M. D-fructose-binding proteins in bull seminal plasma: isolation and characterization. Folia Biol (Praha). 2001;47(4):113–9. pmid:11508854.
- 55. Anand SR. The carbohydrates of buffalo and goat semen. J Reprod Fertil. 1973;32(1):97–100. pmid:4700809.
- 56. Mendoza G, White IG, Chow P. Studies of chemical components of Angora goat seminal plasma. Theriogenology. 1989;32(3):455–66. pmid:16726692.
- 57. Matsuoka T, Imai H, Asakuma S, Kohno H, Fukui Y. Changes of fructose concentrations in seminal plasma and glucose and testosterone concentrations in blood plasma in rams over the course of a year. J Reprod Dev. 2006;52(6):805–10. pmid:16988491.
- 58. Baronos S. Seminal carbohydrate in boar and stallion. J Reprod Fertil. 1971;24(2):303–5. pmid:5551420.
- 59. Jayaraman V, Ghosh S, Sengupta A, Srivastava S, Sonawat HM, Narayan PK. Identification of biochemical differences between different forms of male infertility by nuclear magnetic resonance (NMR) spectroscopy. J Assist Reprod Genet. 2014;31(9):1195–204. pmid:24965760; PubMed Central PMCID: PMCPMC4156941.
- 60. Yousef MI, El-Morsy AM, Hassan MS. Aluminium-induced deterioration in reproductive performance and seminal plasma biochemistry of male rabbits: protective role of ascorbic acid. Toxicology. 2005;215(1–2):97–107. pmid:16098653.
- 61. Cao W, Aghajanian HK, Haig-Ladewig LA, Gerton GL. Sorbitol can fuel mouse sperm motility and protein tyrosine phosphorylation via sorbitol dehydrogenase. Biol Reprod. 2009;80(1):124–33. pmid:18799757; PubMed Central PMCID: PMCPMC2804811.
- 62. Kamp G, Lauterwein J. Multinuclear magnetic resonance studies of boar seminal plasma. Biochim Biophys Acta. 1995;1243(1):101–9. pmid:7827097.
- 63. Williams J, Gladen BC, Schrader SM, Turner TW, Phelps JL, Chapin RE. Semen analysis and fertility assessment in rabbits: statistical power and design considerations for toxicology studies. Fundam Appl Toxicol. 1990;15(4):651–65. pmid:2086311.
- 64. Sorensen MB, Bergdahl IA, Hjollund NH, Bonde JP, Stoltenberg M, Ernst E. Zinc, magnesium and calcium in human seminal fluid: relations to other semen parameters and fertility. Mol Hum Reprod. 1999;5(4):331–7. pmid:10321804.
- 65. Paventi G, Lessard C, Bailey JL, Passarella S. In boar sperm capacitation L-lactate and succinate, but not pyruvate and citrate, contribute to the mitochondrial membrane potential increase as monitored via safranine O fluorescence. Biochem Biophys Res Commun. 2015;462(3):257–62. pmid:25956060.
- 66. Dills WL Jr., Bell LS, Onuma EK. Inhibitory effects of substrate analogs on lactate production from fructose and glucose in bovine spermatozoa. Biol Reprod. 1981;25(3):458–65. pmid:7306636.
- 67. Srivastava A, Chopra SK, Dasgupta PR. Biochemical analysis of human seminal plasma. II. Protein, non-protein nitrogen, urea, uric acid and creatine. Andrologia. 1984;16(3):265–8. pmid:6465552.
- 68. Newairy AS, Salama AF, Hussien HM, Yousef MI. Propolis alleviates aluminium-induced lipid peroxidation and biochemical parameters in male rats. Food Chem Toxicol. 2009;47(6):1093–8. pmid:19425229.
- 69. Yi YJ, Sutovsky M, Kennedy C, Sutovsky P. Identification of the inorganic pyrophosphate metabolizing, ATP substituting pathway in mammalian spermatozoa. PLoS One. 2012;7(4):e34524. pmid:22485177; PubMed Central PMCID: PMCPMC3317647.
- 70. Psychogios N, Hau DD, Peng J, Guo AC, Mandal R, Bouatra S, et al. The human serum metabolome. PLoS One. 2011;6(2):e16957. pmid:21359215; PubMed Central PMCID: PMCPMC3040193.
- 71. Bouatra S, Aziat F, Mandal R, Guo AC, Wilson MR, Knox C, et al. The human urine metabolome. PLoS One. 2013;8(9):e73076. pmid:24023812; PubMed Central PMCID: PMCPMC3762851.
- 72. NI N, KL I. Propylene glycol for dairy cows: A review of the metabolism of propylene glycol and its effects on physiological parameters, feed intake, milk production and risk of ketosis. Animal Feed Science and Technology. 2004;115(3–4):191–213.
- 73. Zhang X, Diao R, Zhu X, Li Z, Cai Z. Metabolic characterization of asthenozoospermia using nontargeted seminal plasma metabolomics. Clin Chim Acta. 2015;450:254–61. pmid:26342261.
- 74. Li Y, Hou X, Qi Q, Wang L, Luo L, Yang S, et al. Scavenging of blood glutamate for enhancing brain-to-blood glutamate efflux. Mol Med Rep. 2014;9(1):305–10. pmid:24220720.
- 75. van der Horst CJ. The occurrence of phenylalanine- and alpha-aminocaprylic acid-alpha-ketoglutarate transaminase in boar spermatozoa. Experientia. 1970;26(3):249. pmid:5417477.
- 76. Moore HD, Hibbitt KG. Fertility of boar spermatozoa after freezing in the absence of seminal vesicular proteins. J Reprod Fertil. 1977;50(2):349–52. pmid:881673.
- 77. Tuli RK, Holtz W. Effect of glycerolization procedure and removal of seminal plasma on post-thaw survival and got-release from Boer goat spermatozoa. Theriogenology. 1994;42(3):547–55. pmid:16727561.
- 78. Graham EF, Pace MM. Some biochemical changes in spermatozoa due to freezing. Cryobiology. 1967;4(2):75–84. pmid:5584798.
- 79. van der Horst CJ, Elgersma A, de Boer MJ, Bestebroer AC. Some biochemical parameters for qualification of bull semen. Tijdschr Diergeneeskd. 1979;104(8):suppl 97–108. pmid:34902.
- 80. D'Amours O, Frenette G, Fortier M, Leclerc P, Sullivan R. Proteomic comparison of detergent-extracted sperm proteins from bulls with different fertility indexes. Reproduction. 2010;139(3):545–56. pmid:19952166.
- 81. Soggiu A, Piras C, Hussein HA, De Canio M, Gaviraghi A, Galli A, et al. Unravelling the bull fertility proteome. Mol Biosyst. 2013;9(6):1188–95. pmid:23392320.
- 82. Kwon WS, Rahman MS, Lee JS, Kim J, Yoon SJ, Park YJ, et al. A comprehensive proteomic approach to identifying capacitation related proteins in boar spermatozoa. BMC Genomics. 2014;15:897. pmid:25315394; PubMed Central PMCID: PMCPMC4287242.
- 83. Kwon WS, Rahman MS, Lee JS, Yoon SJ, Park YJ, Pang MG. Discovery of predictive biomarkers for litter size in boar spermatozoa. Mol Cell Proteomics. 2015;14(5):1230–40. pmid:25693803; PubMed Central PMCID: PMCPMC4424395.
- 84. Kwon WS, Rahman MS, Ryu DY, Park YJ, Pang MG. Increased male fertility using fertility-related biomarkers. Sci Rep. 2015;5:15654. pmid:26489431; PubMed Central PMCID: PMCPMC4614854.
- 85. Vilagran I, Yeste M, Sancho S, Castillo J, Oliva R, Bonet S. Comparative analysis of boar seminal plasma proteome from different freezability ejaculates and identification of Fibronectin 1 as sperm freezability marker. Andrology. 2015;3(2):345–56. pmid:25678437.
- 86. Wennemuth G, Meinhardt A, Mallidis C, Albrecht M, Krause W, Renneberg H, et al. Assessment of fibronectin as a potential new clinical tool in andrology. Andrologia. 2001;33(1):43–6. pmid:11167518.
- 87. Gonzalez-Cadavid V, Martins JA, Moreno FB, Andrade TS, Santos AC, Monteiro-Moreira AC, et al. Seminal plasma proteins of adult boars and correlations with sperm parameters. Theriogenology. 2014;82(5):697–707. pmid:25034516.
- 88. Yoon SJ, Rahman MS, Kwon WS, Ryu DY, Park YJ, Pang MG. Proteomic identification of cryostress in epididymal spermatozoa. J Anim Sci Biotechnol. 2016;7:67. pmid:27895910; PubMed Central PMCID: PMCPMC5117493.
- 89. Yoon SJ, Rahman MS, Kwon WS, Park YJ, Pang MG. Addition of Cryoprotectant Significantly Alters the Epididymal Sperm Proteome. PLoS One. 2016;11(3):e0152690. pmid:27031703; PubMed Central PMCID: PMCPMC4816509.