ORIGINAL RESEARCH article

Front. Anim. Sci., 06 November 2025

Sec. Animal Physiology and Management

Volume 6 - 2025 | https://doi.org/10.3389/fanim.2025.1665783

A multi-omics approach identifies candidate biomarkers predictive of boar cryotolerance and conception rate

  • 1. US Department of Agriculture, Agricultural Research Service, Beltsville Agricultural Research Center (BARC), Animal Biosciences and Biotechnology Laboratory, Beltsville, MD, United States

  • 2. AcuFast, Breese, IL, United States

  • 3. Metabolite Profiling Facility, Bindley Bioscience Center, Purdue University, West Lafayette, IN, United States

Abstract

Introduction:

The swine industry underutilizes cryopreserved boar semen due to poor post-thaw viability and variable fertility outcomes. Current semen evaluation methods are retrospective and insufficient for selecting cryotolerant and fertile sires. The aim of this study was to evaluate multiple reaction monitoring (MRM) profiling as a tool for predicting semen cryotolerance and fertility outcomes.

Methods:

Lipidomic and metabolomic analyses using MRM profiling were applied to fresh and post-thaw ejaculates from 16 commercial Duroc boars with known conception rates (CR) from single-sire matings to identify candidate biomarkers predictive of field CR, post-thaw motility loss, and to determine whether CR markers identified in fresh semen persist post-thaw. Boars were classified by their cryotolerance delta score (CDS), which was calculated as the absolute change in motility between arrival at the cryopreservation laboratory and post-thaw, relative to the average loss in motility (low vs. high), and by field CR (low: 75–79%; mid: 80–89%; high: 90–95%).

Results:

Distinct lipid and metabolite profiles were associated with each phenotype, revealing 20 candidate markers with an AUC ≥ 0.800 (P < 0.05). Markers predictive of higher post-thaw motility loss included compounds producing MRMs tentatively attributed to saturated long-chain fatty acids and elevated metabolites such as kynurenine (AUC = 0.905). MRMs predictive of < 80% CR were attributed to elevated guanosine (AUC = 0.850) and olealdehyde (AUC = 0.815), whereas > 80% CR showed higher abundance of TG(45:4) (AUC = 0.967) and creatine (AUC = 0.800). Candidate markers for CR were distinct from those associated with motility loss and remained detectable in post-thaw samples.

Discussion:

These findings demonstrate that CR and post-thaw motility loss are governed by independent molecular traits and support the development of a multidimensional biomarker-based screening strategy to enhance fertility postthaw. This approach could enable AI centers to improve boar selection and cryopreservation outcomes, ultimately increasing the utility of frozen semen in swine breeding programs.

1 Introduction

With the growing global demand for pork, optimizing genetic progress in swine production has become increasingly critical. Artificial insemination (AI) with fresh, extended semen has been the foundation of genetic progress in swine. However, fresh doses rapidly decline in quality during storage and must be used within 3–5 days of collection, limiting the breeding window, creating bottlenecks, and increasing biosecurity risks (; ; ; ). Cryopreservation of boar semen offers clear advantages, including long-term storage, reduced transport frequency, and compatibility with technologies such as sex-sorted sperm. Despite decades of refinement, fewer than 1% of artificial inseminations use frozen boar semen, and these are primarily employed to reintroduce elite genetics, heritage breeds, or high-value individuals (). The primary barrier to adoption lies in the performance gap between frozen–thawed and fresh semen, where only 50–60% of sperm survive the freeze–thaw process, leading to impaired motility, morphology, and membrane function (; Yeste, 2015, Yeste, 2016). Cryopreservation also induces changes resembling premature capacitation and causes disorder of plasma-membrane lipids, further compromising fertilization capacity, cell signaling, and interaction with the female reproductive tract (; ). These impairments ultimately reduce the fertility of the breeding dose and shorten the lifespan of spermatozoa post-thaw compared with fresh semen (, ; ).

Another major hurdle for the integration of genetic-advancement methods such as frozen–thawed semen is the lack of readily available tools capable of predicting fertility outcomes before insemination. Evaluation of boar fertility currently requires ≥ 50 single-sire matings (), and routine semen-quality assessment for motility and morphology still allows subfertile boars to enter the breeding herd (; ; ). Industry data suggest that approximately 25% of boars achieve conception rates below 80%, highlighting a significant subpopulation that borders or falls below culling thresholds—emphasizing the need for more sensitive selection tools (). While effective, single-sire matings used for fertility testing can result in negative economic consequences for the sow farm and cannot feasibly be deployed for every boar. Therefore, tools that enable fertility-prediction screening for each sire—and thus predict frozen–thawed breeding-dose fertility—would be highly valuable to the swine industry.

Omics technologies, particularly metabolomics and lipidomics, provide the sensitivity needed to resolve subtle molecular signals likely driving the prevalence of subfertility. Moreover, lipidomic and metabolomic markers are well suited for fertility-biomarker discovery because they capture the molecular phenotype most proximal to physiological conditions. Like fertility traits themselves, the metabolome and lipidome are highly responsive to environmental influences, making them sensitive indicators of reproductive potential (; ; ). In boars, lipidomic and metabolomic profiling of sperm and seminal plasma has revealed associations between energy and lipid metabolism, reproductive performance, capacitation, and cryotolerance (; Zhang et al., 2021; ; ) Differentially expressed metabolites linked to sperm function, motility, and membrane integrity are now being mapped across species, breeds, and sample types, suggesting the potential for robust cross-species biomarkers predictive of both fertility and post-thaw semen quality (; ; ). Recent work by has significantly advanced this area by identifying metabolomic biomarkers of sperm freezability in Duroc, Landrace, and Large White boars. The presence of lipidomic markers of cryotolerance has also been reported in Yorkshire boars (Zhang et al., 2023).

However, traditional exploratory approaches for lipidomic and metabolomic marker discovery are time-consuming, creating a bottleneck for herd-level validation of omic cryotolerance or fertility markers. Multiple reaction monitoring (MRM) profiling, an alternative biomarker-discovery approach, enables high-throughput screening of informative lipid and metabolite features through ion monitoring without full structural elucidation, thereby reducing instrument time and improving diagnostic efficiency (). MRM profiling has already been applied in several fertility-related contexts, including biomarker identification in prepubertal gilts, broiler-breeder roosters, and in distinguishing fertile from subfertile men (; ; ). In boars, MRM profiling has been used to examine age-related changes in the lipidome and to identify correlations with on-farm fertility measurements such as motility and morphology (). Beyond fertility, MRM profiling has also been extended to metabolic and developmental contexts, including biomarker discovery for polycystic ovary syndrome (PCOS) in women, bovine oocyte and embryo-quality screening, and assessment of perinatal diet effects on gilt reproductive development (; ; ; ; ; ; ; ). Taken together, we hypothesized that MRM profiling could elucidate highly informative features capable of predicting boar fertility and cryotolerance traits.

The overarching goal of this work is to develop biologically informed, pre-freeze selection tools that can be deployed within AI centers to optimize selection decisions in the future. Moreover, the ability to predict post-thaw semen quality—including fertility outcomes—would enhance the commercial adoption of frozen–thawed semen and support the broader implementation of reproductive technologies such as sex-sorted semen, thereby improving the global competitiveness of swine-breeding programs. The objective of this study was to evaluate MRM profiling as a tool for predicting semen cryotolerance and fertility outcomes. Specifically, we aimed to (1) identify potential biomarkers predictive of post-thaw motility loss and conception rate in fresh semen, (2) characterize informative features of post-thaw motility loss in post-thaw samples, and (3) assess how cryopreservation alters the abundance of conception-rate markers.

2 Materials and methods

2.1 Animals, sample collection, and conception rate classification

All samples for the study were sourced from ejaculates of 16 sexually mature Duroc boars enrolled in a commercial breeding program during routine collection. Briefly, boars were moved into a collection crate where preputial fluid was removed manually using a double-gloved hand. The sheath area and extended penis were cleaned and dried with a paper towel. The outer glove was removed, and manual collection was completed using the clean glove. The gel and pre-sperm fractions were collected into an inner bag that was detached and discarded prior to collecting the sperm-rich fraction in the collection mug. When collection of the sperm-rich fraction was complete, the collection mug was removed from under the collection area, ensuring only the sperm-rich portion of the ejaculate was collected for cryopreservation. The filter was removed, and the collection bag was closed using a band and sent to the laboratory via pneumatic tube.

Upon arrival at the stud laboratory, a 2 mL aliquot of the sperm-rich fraction was placed into a polypropylene conical tube (Falcon, USA) containing 10 mL of methanol (1:5, v/v). The sample was inverted five times to thoroughly mix the ejaculate and methanol (). Samples were packaged and shipped on cooling packs to Beltsville Agricultural Research Center (BARC) overnight, and frozen upon arrival at −80 °C until lipidomic extraction. The remainder of the ejaculate was evaluated using standard boar stud protocols for semen quality assessment using computer-assisted sperm analysis (CASA). Boars met or exceeded the boar stud’s minimum quality thresholds of 75% motile sperm and 75% normal morphology during routine evaluation. Semen was extended 1:1 (v/v) in AndroStar® Plus (Minitube, Verona, WI, USA) and shipped overnight on cooling packs to the AcuFast cryopreservation laboratory.

Conception rates (CR) for each boar were supported by historical records from single-sire matings produced during routine fertility screening. Sires that fall below 75% CR are considered subfertile and are culled after sufficient data are available (≥50 matings), which is reflected in our cohort (75-95%; ≥ 50: n=13; <50: n=3). With 13–17% of the breeding herd estimated to be subfertile () and the prevalence of <80% CR boars estimated to be ~25% of the herd () identifying and monitoring the remaining 8–12% that border cull thresholds would be valuable. To address this, we stratified boars into three operational bands reflecting herd-level management decisions and the cohort’s conception rate distribution: Low (75–79%; n = 6; “Watch-list”), Mid (80–89%; n=6; “Acceptable”), and High (90–95%; n = 4; “Excellent”). The three individual boars with <50 single-sire matings were distributed equally across CR groups: 47 (79%), 26 (85%), and 42 (95%).

2.2 Cryopreservation and cryotolerance phenotype classification

Extended ejaculates arrived approximately 24 hours post-collection (referred to as check-in) and were evaluated for semen quality estimates including concentration, motility, and morphology. The freezing process was initiated following check-in evaluation. Briefly, semen was extended further to 3:1 (v/v) with AndroStar® Plus (Minitube, Verona, WI, USA) and centrifuged at 10,000 × g at 20 °C for 10 minutes. Seminal plasma was discarded, and the remaining sperm cell pellet was resuspended in a proprietary egg yolk-based cooling extender. Samples were cooled for 1.5 hr to 5 °C and diluted in another proprietary, cryoprotectant-supplemented extender to a final concentration of 450–500 million cells/mL. Motility was evaluated again prior to freezing for quality control purposes, and extended semen was packaged in 0.5 mL straws (IMV Technologies, Osseo, MN, USA) and frozen in a Digitcool programmable freezer (IMV Technologies, Osseo, MN, USA). Straws were submerged into liquid nitrogen and stored for a minimum of 1 day prior to post-thaw semen quality evaluation. A subset of frozen straws was shipped to BARC in a dry shipping container and stored in a liquid nitrogen dewar until thawed for lipidomic and metabolomic sample processing.

The same thawing protocol was used to evaluate post-thaw motility and prepare samples for lipidomic and metabolomic (omic) screening. Briefly, straws were thawed at 38 °C for 30 s in a water bath. Straw contents were deposited into an empty, pre-warmed (38 °C) tube, then extended with 5 mL of AndroStar® Plus and incubated at 38 °C for 20 min. Samples were evaluated for 20-min motility per standard straw quality control protocols, which are used to determine the number of straws for a breeding dose (2 billion motile cells). Samples for omic screening were processed in methanol using the 1:5 (v/v) method described above and stored at −80°C until lipid and metabolite extraction.

Cryotolerance phenotype groups were defined by the relative loss of motility during cryopreservation. A cryotolerance delta score (CDS) was calculated as the absolute change (Δ) in total motility between initial check-in (~24 hours post-collection) and 20 min post-thaw. To account for differences in starting motility, CDS values were normalized to the cohort mean and stratified by standard deviation. Boars with stable or minimally reduced motility post-thaw (Δ ≥ cohort mean; n = 7) were classified as high-CDS, whereas those with larger-than-average reduction in motility (Δ < cohort mean; n = 9) were classified as low-CDS. This approach ensured classification reflected relative cryotolerance performance within the cohort rather than arbitrary cutoffs. Group-level summary statistics for CDS and conception rate are provided in Table 1.

Table 1

Phenotype1Fresh conception rate, %Average change (Δ) in absolute motility, %Check-in motility2, %Post-thaw motility3, %
CDS4
High (n=7)83.6 (76.0-89.0)27.3 (15.7-37.3)64.0 (36.8-84.5)36.7 (15.9-47.2)
Low (n=9)85.2 (75.0-95.0)44.8 (39.9-51.5)73.1 (55.8-88.3)28.3 (15.9-39.7)
Conception Rate5
High (n=4)92.3 (90.0-95.0)44.8 (39.9-51.5)73.7 (55.8-88.3)28.9 (15.9-37.6)
Mid (n=6)86.5 (82.0-89.0)29.9 (15.7-43.0)64.4 (36.8-82.7)34.5 (15.9-47.1)
Low (n=6)77.3 (75.0-79.0)39.3 (17.6-47.3)70.8 (59.7-84.5)31.5 (19.6-47.2)

Averages and ranges of semen quality parameters, loss in post-thaw motility, and conception rate by cryotolerance delta score (CDS) and conception rate phenotype group.

1Boars enrolled in the study ranged from 14.3-36.9 months of age, with an average of 20.9 months. Age distribution within conception rate or cryotolerance groups was similar.

2Check-in motility refers to total motility measured following transportation to the cryopreservation laboratory from the boar stud (24 hours post-collection).

3Post-thaw total motility was measured 20 minutes post-thaw.

4Cryotolerance delta score (CDS) group was quantified as the change (Δ) in absolute total motility between initial check-in (sample arrival at the cryopreservation laboratory) and 20 minutes post-thaw, relative to the standard deviation from the group’s average Δ in absolute motility.

5Conception rates (CR) for each boar were supported by single-sire matings based on historical records, ranging from 75-95%. Thirteen boars had ≥ 50 matings, and three boars had <50 matings. Boars with <50 matings were distributed equally across CR groups: 47 (79% CR), 26 (85%), & 42 (95% CR).

2.3 Lipid and metabolite extraction

Lipids were extracted from ejaculates using the Bligh and Dyer method (). Samples were thawed at room temperature, then sonicated for three minutes at 40 kHz in an ultrasonic bath (Bransonic® M2800, Emerson Electric Co., St. Louis, MO, USA). Two hundred microliters of sample were transferred to a new 1.7 mL tube (Axygen®, Corning, New York), and then 550 µL of methanol prepared with 1 mM butylated hydroxytoluene (BHT) and 250 µL of chloroform were added and vortexed for 10 s so that the final concentration of BHT in the solution was 550 nM. Samples were then incubated at 4 °C for 15 min. The addition of 250 µL ultrapure water and 250 µL chloroform induced the formation of two immiscible liquid phases: a lower organic phase (lipids) and an upper polar phase (metabolites). Samples were centrifuged at 10,000 x g for 5 minutes at 4 °C for better separation of the two-phase solution. The lower (lipid) phase was removed and transferred to a new 1.7 mL microcentrifuge tube, followed by transferring the upper (metabolite) phase to a separate 1.7 mL microcentrifuge tube. Samples were dried in a vacuum concentrator for 8 h and stored at −80 °C until ready for MRM profiling.

2.4 MRM profiling of lipids and metabolites

Dried pellets were resuspended in 200 µL of solvent containing acetonitrile, methanol, and 300 mM ammonium acetate (3:6.65:0.35, v/v/v), to yield a final concentration of 10 mM NH4Ac. A further 150× dilution of the sample in solvent was used for direct injection. Eight microliters from each sample were injected into the micro autosampler (G1377A) of a QQQ6410 triple quadrupole mass spectrometer (Agilent Technologies, San Jose, CA) equipped with an ESI ion source. A solvent solution containing acetonitrile with 1% formic acid at 10 µL/min was pumped between injections (CapPump G1376A, Agilent Technologies, San Jose, CA). Pure methanol was injected between samples to remove any remaining lipids from the previous injection.

Multiple reaction monitoring (MRM) profiling was performed in two phases, beginning with a discovery phase that involved screening a large database of MRMs in composite samples (one sample for each experimental group). Samples were pooled in the discovery phase based on CDS phenotype (“high” and “low”) for fresh and post-thaw samples to identify lipid classes of importance prior to individual sample interrogation. MRMs detected in the discovery phase within fresh and post-thaw samples were then used to interrogate individual samples during the screening phase. Samples were screened for the chemical classes acylcarnitine (AC), ceramide (Cer), free fatty acids (FA), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylinositol (PI), phosphatidylglycerol (PG), phosphatidylserine (PS), diacylglycerol (DG), and triacylglycerol (TG) using MRM profiling methods previously reported (). Processing of the initial chemical class data was completed using MSConvert20, which converted each set of profiling method data into mzML file format. Signal intensity for ions present in the samples was obtained using an in-house script. Ions with ion count values >30% in at least one of the samples compared with a blank within each profiling method were selected for the screening phase. Due to the large number of MRMs to be screened, the samples were interrogated using four lists of MRMs, which we refer to as methods:

Metabolomic analysis was completed during the screening phase and samples were analyzed in both the positive (Supplementary Table 5) and negative (Supplementary Table 6) ion modes. Question marks within metabolite names in the Supplementary Tables indicate unknown structural orientation of the molecule.

2.5 Data processing and statistical analysis

MRM ion pairs were acquired for all methods except for Method 1, where fatty acids did not fragment without derivatization during MRM profiling. As a result, MRMs reported in Method 1 only contain the precursor ion m/z. Tentative lipid class attributions based on associated functional group and biological information were assigned to unidentified MRM ion pairs and precursor ions to increase potential biomarker discovery using lipid maps (http://www.lipidmaps.org/). For data quality purposes, MRM intensities less than 1.3-fold (lipids) or 1.5-fold (metabolites) below the blank were removed. Next, relative ion intensity was calculated by dividing intensity by the sum of intensities of all lipids within a sample by screening method. Relative intensities of MRM ion pairs were uploaded into MetaboAnalyst 6.0 () and normalized using autoscaling for data visualization, statistical analysis, and biomarker analysis. Statistical analyses included Student’s t-tests to compare cryotolerance phenotypes and one-way ANOVA with Fisher’s post hoc tests to assess differences across conception rate groups. Because conception rate is inherently a continuous variable, Pearson’s correlation analysis was also performed to evaluate the relationship between conception rate and the abundance of informative biomarkers. Data visualizations included hierarchical clustering heat maps and box plots of significant features, with an α-level of 0.05 based on nominal P-values.

Biomarker analyses included receiver operating characteristic (ROC) analysis comparing lower conception rate boars (<80%, n = 6) to all others (n = 10), reflecting the estimated field prevalence where ~25% of boars fall below 80% conception rate (). This approach enabled us to evaluate whether omics markers could identify borderline boars approaching the cull threshold, thereby demonstrating the sensitivity of the screening approach without requiring phenotype extremes. Predictive performance was visualized in ROC curves and quantified by area under the curve (AUC). The following AUC scale was used to evaluate potential biomarkers: excellent = 0.900–1.000, good = 0.800–0.899, fair = 0.700–0.799, poor = 0.600–0.699, and fail = 0.500–0.599 (). Lipids and metabolites were considered candidate biomarkers if they met either of the following criteria: (1) Pearson’s correlation coefficient (r) ≥ 0.6 with P ≤ 0.05, or (2) area under the ROC curve (AUC) ≥ 0.800 with P ≤ 0.05. Features with AUC ≥ 0.800 and P-values approaching significance (0.05 < P ≤ 0.06) were also retained due to their high discriminatory potential. Tentative attributions that share an MRM are reported with each potential compound separated by a forward slash.

3 Results

3.1 Potential biomarkers of cryotolerance score in fresh semen

Fresh semen profiling distinguished high and low cryotolerance score (CDS) boars. High-CDS samples showed enrichment of select phosphatidylcholines and sphingomyelins, whereas low-CDS boars were characterized by greater abundance of long-chain fatty acids and additional phospholipid and sphingomyelin features. Metabolite analysis revealed enrichment of multiple small molecules in high-CDS semen, including features detected in both positive and negative ion modes. A total of 13 candidate markers for CDS score were identified (AUC ≥ 0.810). Three features with tentative attributions as 1,2-dihydroxyheptadec-16-Yn-4-Yl acetate, heneicosanoic acid, and nonadecanoic acid achieved AUC > 0.900 and had the strongest predictive value (Table 2). Predictive performance and distribution patterns for the top biomarkers identified within screening methods for lipids and metabolites are visualized in Figure 1.

Table 2

MethodTentative id attribution (MRM)1Common nameAUC2P-valFC3 (low-high)
Fresh markers
Method 14FOH 19:3;031,2-Dihydroxyheptadec-16-Yn-4-Yl Acetate0.9050.01-0.28
FA 21:0Heneicosanoic acid0.9050.01-0.28
FA 19:0Nonadecanoic acid0.9050.02-0.24
FA 20:0Arachidic acid0.8890.01-0.27
Method 2PC(36:1)/PC(O-37:1)PC(36:1)/PC(O-37:1)0.8730.020.58
PC(34:4)/PC(O-35:4)/PC(dO-36:4)PC(34:4)/PC(O-35:4)/PC(dO-36:4)0.8730.02-0.25
SM(38:3;O2)SM(38:3;O2)0.8570.05-0.20
SM(40:1;O2)SM(40:1;O2)0.8100.020.30
PE(37:5)/PE(O-38:5)/PE(O-38:5)PE(37:5)/PE(O-38:5)/PE(O-38:5)0.8100.05-0.27
Metabolite (+) Ion ModeKynurenine (209.1 -> 94.1)Kynurenine0.9050.01-0.61
Cytidine (242.1 -> 110)Cytidine0.889<0.01-0.89
D-Biotin (243.1 -> 200.1)D-Biotin0.8410.02-0.45
Metabolite (-) Ion ModeHexose sugar (179.1 -> 71)Fructose or Glucose0.8570.03-0.21
Post-thaw markers
Method 1FA 38:0Octatriacontanoic acid0.8410.020.30
Method 2LPC(16:0)/PC(O-16:0)/LPC(O-17:0)LPC(16:0)/PC(O-16:0)/LPC(O-17:0)0.8570.01-0.24
LPC(22:6)LPC(22:6)0.8410.02-0.23
PC(36:7)/PC(35:0)/PC(O-36:0)PC(36:7)/PC(35:0)/PC(O-36:0)0.8410.03-0.20
Cer(d18:1/16:0)Cer(d18:1/16:0)0.8100.02-0.21
PE(40:4)PE(40:4)0.8100.03-0.16
PC(40:4)PC(40:4)0.8100.05-0.12
Method 4TG(51:9)/TG(50:2) NL 16:0TG(51:9)/TG(50:2) NL 16:00.9050.010.10
TG(52:4) NL 16:0TG(52:4) NL 16:00.8570.020.10
TG(53:14;O2)/TG(54:5) NL 18:1TG(53:14;O2)/TG(54:5) NL 18:10.8570.02-0.11
TG(48:2) NL 18:2TG(48:2) NL 18:20.8410.010.17

Candidate biomarkers associated with cryotolerance delta score (CDS) in fresh and post-thaw semen.

1Tentative ID attributions that share MRM ions or ion pairs are reported together and separated by a forward slash.

2AUC, area under the curve from receiver operating characteristic (ROC) analysis.

3FC, Fold change of low-CDS minus high-CDS.

4Compounds identified in Method 1 did not fragment during MRM profiling. Tentative attributions are based on the precursor ion rather than the MRM ion pair.

Figure 1

3.2 Characteristic features of cryotolerance phenotype in post-thaw samples

Because all boars were frozen in the same extender, major differences were not expected; however, post-thaw profiling revealed distinct signatures between CDS phenotypes. High-CDS samples showed enrichment of select fatty acids, such as FA 38:0, and triacylglycerols, whereas low-CDS ejaculates accumulated lysophospholipids, phosphatidylcholines, and ceramides (Table 2). These post-thaw features met biomarker thresholds (AUC ≥ 0.81), while no informative associations were identified for diacylglycerols or metabolites (all AUC < 0.80, P > 0.05; full list provided in Supplementary Table 7).

3.3 Candidate markers predictive of conception rate

Despite not selecting individuals based on fertility extremes, Pearson’s correlation and ROC analysis identified several features exhibiting strong associations with conception success. Fertility-linked metabolites included both positively and negatively correlated features, while distinct lipid classes also showed strong associations with conception outcomes. Thirteen markers exceeded our candidate biomarker thresholds (|r| ≥ 0.6, AUC ≥ 0.800; Table 3; Figure 2), highlighting their potential to capture fertility variation within commercially relevant ranges. Notably, TG(45:4) NL 20:0 emerged as the most robust conception rate predictor (r = 0.8, AUC = 0.967; Figure 2). These markers effectively stratified boars by fertility outcome but not by post-thaw cryotolerance score class (Tables 2; Table 3).

Table 3

MethodTentative id attribution (MRM)2Common nameAUC3P-valueR-value4Pearson’s p-value
Method 2SM(42:3;O2)SM(42:3;O2)0.8000.050.60.01
PE(34:0)/PE(O-35:0)PE(34:0)/PE(O-35:0)0.8000.050.60.03
PC(40:4)PC(40:4)0.8330.020.60.01
Method 4TG(53:7)/TG(52:0) NL 16:0TG(53:7)/TG(52:0) NL 16:00.8500.060.60.02
TG(37:0) NL 20:0TG(37:0) NL 20:00.8670.020.7<0.01
TG(45:4) NL 20:0TG(45:4) NL 20:00.967<0.010.8<0.01
Metabolite (+) Ion ModeCreatine (114.1 -> 43)Creatine0.8000.060.60.01
Guanosine (284.1 -> 152.2)Guanosine0.8500.04-0.60.02
Metabolite (-) Ion ModeOlealdehyde (265.2 -> 97)Olealdehyde0.8150.01-0.60.02
Fumaric acid (115 ->71)Fumaric acid0.8330.03-0.60.01
Pantothenic acid (218.1 -> 88)Pantothenic acid0.8330.04-0.70.01
Hexose sugar (179.1 -> 89)Fructose or Glucose or Mannose0.8330.04-0.60.03
6-Hydroxynicotinic acid (138 -> 94)6-Hydroxynicotinic acid0.8890.01-0.7<0.01

Candidate biomarkers predictive of conception rate in fresh semen1.

1Conception rates for each boar were based on historical records produced from ≥50 single sire matings.

2Tentative ID attributions that share MRM ions or ion pairs are reported together and separated by a forward slash.

3AUC, area under the curve from receiver operating characteristic (ROC) analysis comparing low conception rate (<80%, n = 6) versus all others (n = 10), chosen to reflect the estimated field prevalence in which ~25% of boars fall below 80% conception rate ().

4Pearson’s correlation coefficient describing the linear relationship between candidate biomarker abundance and conception rate (continuous variable).

Figure 2

3.4 Persistence of fertility markers post-thaw and relationship with cryotolerance

Important fertility-associated features originally identified in fresh semen persisted post-thaw, but their relative abundance by fresh CR phenotype was not maintained within boar (Figure 3). Several low-CDS boars exhibited higher relative levels of CR-linked lipids such as TG(45:4) NL 20:0. By contrast, some high-CDS boars displayed lower abundance of these same markers post-thaw. These molecular patterns demonstrate the biological independence of post-thaw motility and conception potential, highlighting that post-thaw motility loss is not a complete indicator of frozen-thawed fertility. This further supports integrating multi-trait, pre-thaw screening to identify sires that produce breeding doses that remain both motile and fertile following cryopreservation, thereby optimizing cryotolerance selection. Boars who retained relative abundance of fertility markers pre- and post-thaw, represented as a percentage of their cryotolerance phenotype, are summarized in Figure 4.

Figure 3

Figure 4

4 Discussion and conclusions

In the present study, our findings support the existence of mechanistically distinct lipidomic and metabolomic signatures for sperm cryotolerance and fertility. However, MRM profiling is an exploratory ion screening tool, and tentative attributions for lipids and metabolites must be validated using approaches such as UHPLC or HRMS to determine final concentration and chemical structure of these small molecules. To help support our exploratory results and biological interpretation, tentative attributions were cross-referenced to previous work that characterized the lipidome and metabolome of boar spermatozoa and seminal plasma. Tentative attributions in the present study are consistent with small molecules previously reported in boar samples characterized with UHPLC-MS, GC-MS, or LC-MS/MS (; Zhang et al., 2023; ; ).

4.1 Informative lipid features suggest post-thaw motility loss is related to degree of lipid saturation

High motility loss ejaculates (low-CDS) were enriched in saturated fatty acids, including nonadecylic (FA 19:0), arachidic (FA 20:0), and heneicosanoic acid (FA 21:0). Saturated lipids tend to pack tightly in the lipid bilayer, making the plasma membrane more rigid and less flexible to temperature fluctuations (). Because cryotolerance was characterized based on motility in the present study, these findings are consistent with reports that low-motility boars have a higher proportion of saturated fatty acids in fresh ejaculates (). Lipid profiles of sperm cells from poor freezing Yorkshire boars were enriched in saturated lipids, and their abundance was correlated with weaker plasma membrane integrity, greater mitochondrial dysfunction, and poorer post-thaw motility and viability (Zhang et al., 2023). In contrast, monounsaturated fatty acids appear protective, promoting more flexible membranes; enrichment of these species was a trend in our low motility loss (high-CDS) samples and consistent with reports in Yorkshire boars (Zhang et al., 2023). However, samples enriched with monounsaturated phosphatidylcholines and sphingomyelins were more cryotolerant, indicating that there is a balance required to survive colder temperatures where moderate unsaturation is favored.

Polyunsaturated fatty acids (PUFAs) add further nuance. PUFAs are characteristic of boar spermatozoa and have been associated with fertility in fresh semen (; ). Increased PUFA composition increases plasma membrane fluidity and is typically characteristic of cells more resistant to ice crystal damage (). However, while PUFAs aid membrane fusion and function, their high PUFA concentrations also render boar spermatozoa prone to oxidative degradation which is amplified during the cryopreservation process (). Peroxidized lipids disrupt membrane integrity, mitochondrial function, and cell viability—hallmarks of sperm subjected to cold shock and cryopreservation (). Consistent with our findings in saturated fatty acids, Zhang et al. (2023) observed that poor freezing were enriched with PC and PE lipids containing arachidonic acid. Although PUFAs such as α-linolenic acid (18:3), arachidonic acid (ARA 20:4), and docosahexaenoic acid (DHA 22:6) are beneficial for fertility of fresh semen, they may have an opposite effect during cryopreservation (). Together, these results suggest that cryosurvival depends on a balanced lipid composition: overly rigid membranes resist remodeling, while excessive unsaturation leaves sperm vulnerable to cold-induced peroxidation. Future studies should quantify specific PUFA-containing lipid species in boar spermatozoa and evaluate their inherent risk for cryodamage.

Similar patterns persisted in post-thaw lipidomic signatures. Ejaculates with high motility loss were enriched in polyunsaturated LPCs and PEs, including LPC(22:6), which has been linked to oxidative membrane damage in cattle and roe deer sperm (). Cer(d18:1/16:0), a lipid previously linked to increased membrane permeability and induction of apoptosis via mitochondrial pathways (; ) was also associated with high motility loss in the present study. Seminal plasma extracellular vesicles containing elevated ceramides have been associated with low motility in fresh semen (). Although processing protocols typically minimize seminal plasma inclusion (), our findings suggest that vesicle-derived lipids may still influence motility post-thaw and merit further study.

Tentative attributions for other characteristic features in post-thaw samples—such as FA 38:0 and the triacylglycerols—have not been previously identified in boar spermatozoa or seminal plasma. These compounds are likely components of the cryopreservation extender rather than post-thaw modifications of the sperm cell itself. Interestingly, even though all ejaculates were frozen in the same extender, the relative abundance of presumed extender components differed by CDS phenotype. If extender composition varies enough from straw-to-straw to influence cryotolerance and post-thaw sperm metabolism, optimizing extender distribution could represent a new area for improvement.

4.2 Metabolic indicators of oxidative stress differentiate conception rate and motility loss outcomes

Metabolomic profiling revealed that candidate markers for post-thaw motility loss and lower conception rates were associated with elevated metabolites indicative of oxidative stress and inflammation. Kynurenine, a product of tryptophan degradation that reflects oxidative stress and immune activation detrimental to sperm function, was higher in high motility loss ejaculates (, ). 6-hydroxynicotinic acid, another metabolite linked to the kynurenine pathway, was inversely correlated with conception rate—but not informative for post-thaw motility loss (AUC<0.6, P>0.1)—suggesting broader dysregulation of NAD metabolism in fresh semen (). Other stress-linked metabolites included cytidine, biotin, olealdehyde, and excess sugar, each consistent with oxidative or metabolic imbalance (; ; ; ; ). Notably, guanosine was elevated in lower conception rate ejaculates and is known to be the most susceptible spermatozoal nucleoside to oxidative damage, with direct mechanistic links to male infertility (; ). In contrast, creatine was the only metabolite positively associated with conception rate. Its role as a substrate for phosphocreatine-driven ATP buffering is well established, and in sperm this pathway is critical for maintaining motility and energy supply during capacitation and fertilization (; ). When added to insemination media, exogenous creatine improves motility, velocity, and enhances capacitation during in vitro fertilization in mice and swine (, ). Thus, elevated creatine may reflect an energetically competent metabolic profile that supports fertilization success. While these pathways must be validated in larger cohorts, the recurring theme is that sperm experiencing—or more susceptible to—oxidative or metabolic stress pre-freeze are less resilient to cryopreservation and achieve lower conception rates, although different metabolites appear to underlie each trait.

4.3 Limitations, future directions, and implications

Several limitations must be acknowledged. The sample size (n = 16) and use of a single ejaculate per boar limit generalizability and preclude evaluation of intra-individual consistency. Future studies should include repeated measures, larger cohorts, and multiple breeds to validate biomarkers across a variety of genetic backgrounds. Moreover, field fertility in this study was based on fresh semen; direct fertility trials using cryopreserved semen remain necessary to confirm functional competence. Environmental and nutritional influences on biomarker expression, as well as optimization of extender formulations to preserve favorable lipid profiles, also warrant investigation. Analytical limitations of MRM profiling must also be considered, including reduced fragmentation of underivatized fatty acids and occasional ambiguous attributions. While this underscores the exploratory nature of MRM profiling, reproducible precursor ions and the strong predictive power of identified MRMs justify further targeted validation using high-resolution methods.

This approach addresses key limitations of current workflows that evaluate cryopreservation and conception rate outcomes, which rely heavily on motility and morphology as proxies for reproductive potential. Given the relatively lower abundance of conception rate biomarkers in post-thaw samples from some high-cryotolerance boars, current selection criteria may result in motile breeding doses with relatively less fertility compared to other sires. Moreover, the persistence of fertility markers post-thaw highlights the potential for using molecular screening as an additional quality control layer beyond motility, ensuring that ejaculates not only meet minimum viability thresholds but also retain key fertilization-associated metabolites.

Another important finding was that high-conception boars often overlapped with low cryotolerance phenotypes. For example, polyunsaturated PCs and SMs predicted higher conception rates but were also linked to greater post-thaw motility loss. This decoupling of fertility and cryotolerance observed in lipids and metabolites has practical implications for sire selection, suggesting that optimal candidates for cryopreservation may fall within mid-range conception rates (80–89%), where lipid balance favors both fertility and resilience to freezing. Our findings emphasize the importance of integrated selection systems and strongly support the implementation of a multi-axis evaluation framework for boar semen quality.

4.4 Conclusions

This study delineates distinct lipidomic and metabolomic signatures predictive of post-thaw motility loss and conception success in boars. The molecular dissociation of these two traits provides a framework for reimagining semen quality assessment—not as a unidimensional construct based solely on motility, but as a multidimensional phenotype amenable to precision molecular diagnostics. The deployment of MRM-based biomarker panels in AI centers offers a transformative opportunity to enhance the utility of cryopreserved semen, minimize genetic loss, and improve reproductive efficiency in swine genetic programs.

Statements

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 author.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because boars were not collected for the sole purpose of this study. Samples were sourced from ejaculates collected from boars enrolled in a commercial breeding program as part of routine collections per standard husbandry protocols.

Author contributions

KM: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing. AM: Investigation, Project administration, Resources, Writing – review & editing. JB: Conceptualization, Investigation, Project administration, Resources, Writing – review & editing. KW: Conceptualization, Investigation, Project administration, Resources, Writing – review & editing. CF: Data curation, Investigation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This research was funded by the in-house USDA-ARS CRIS project number 8042-31000-111-00D to Kayla Mills.

Acknowledgments

All lipidomic and metabolomic data collection and analysis were performed at Purdue University’s Metabolite Profiling Facility in the Bindley Bioscience Center. We are also very appreciative for the generous contributions made by AcuFast of all boar ejaculates used in this study.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fanim.2025.1665783/full#supplementary-material

References

Summary

Keywords

boar, cryopreservation, lipidome, biomarker, metabolome

Citation

Mills K, Minton A, Berndtson J, Willenburg K and Ferreira CR (2025) A multi-omics approach identifies candidate biomarkers predictive of boar cryotolerance and conception rate. Front. Anim. Sci. 6:1665783. doi: 10.3389/fanim.2025.1665783

Received

14 July 2025

Accepted

14 October 2025

Published

06 November 2025

Volume

6 - 2025

Edited by

Muhammet Rasit Ugur, IVF Michigan Fertility Centers, United States

Reviewed by

Naseer A. Kutchy, University of Alabama at Birmingham, United States

Martin Ntawubizi, University of Rwanda, Rwanda

Updates

Copyright

*Correspondence: Kayla Mills,

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics