The effects of nutritional level and body condition score on cytokines in seminal plasma of beef bulls

Introduction High quality semen is essential for reproductive efficiency. Nutrition and environmental factors impact the ejaculate components, like cytokines, that are essential for pregnancy establishment. We hypothesized that differing nutritional periods and body condition scores would affect seminal plasma cytokine concentrations of bulls. Methods Mature Angus bulls (n=11) were individually housed and randomly assigned to either over-fed (n=5) or restricted (n=6) treatment pathways. Bulls were fed different volumes of a single ration creating 8 individual treatment periods. Body weight and body condition scores were taken every 14 d to manage intake volumes. Ejaculates were collected every 84 d to determine seminal plasma cytokine profiles. A complete randomized design was used to evaluate seminal plasma cytokines after each nutritional treatment. Initial cytokine concentrations and volume of the ejaculate were included as covariates. Results All cytokines returned to initial concentrations following maintenance treatments at an ideal body condition score of 6. Nutritional treatments affected (P ≤ 0.05) IFN-γ, IL-8, MIP-1α, MIP-1β, TNF-α, IL-1β, and VEGF-A. However, TNF-α, IFN-γ, and MIP-1α had the greatest impact on cytokine profile. Discussion Nutritional levels and adiposity altered seminal plasma cytokine concentrations which could potentially impact the inflammatory balance of the uterus and the immune responses necessary for pregnancy establishment.


Introduction
A cow-calf producer's primary income is highly dependent on bull fertility since profits are based on the number of calves born each year (Taylor and Field, 1995). Monitoring bull fertility is critical to maximizing reproductive efficiency in beef production (Butler et al., 2020). One infertile bull is more detrimental than one infertile cow due to the number of offspring that bulls produce in a given breeding season (Kastelic, 2013). Bull management for maximizing fertility includes proper nutrition to produce bulls that fit industry demands. However, bull fertility and reproductive success can become limiting when issues such as inadequate nutrition, injuries, or infections are present (Kastelic, 2013). Correct nutritional management of bulls is critical due to sires experiencing variable periods of nutrition throughout the year. Therefore, understanding the complete impacts of nutrition on seminal plasma (SP) and male fertility could help increase reproductive efficiency.
Seminal plasma is the non-cellular portion of the ejaculate composed of cytokines, amino acids, enzymes, hormones, ions, sugars, lipids, antioxidants, and proteins that act as a nutritional, protection, and transport media for spermatozoa as it travels through the female reproductive tract (Juyena and Stelletta, 2012). Moreover, SP moieties can be impacted by many factors including nutrition (Binder et al., 2015;Watkins et al., 2018;Schjenken et al., 2021). Restricted diets can decrease the amount of SP volume produced in each ejaculate (Singh et al., 2018), while overfeeding increases the production of insulin-like growth factor 1, potentially influencing sperm production and libido of rams (Fourie et al., 2004;Selvaraju et al., 2012). Proinflammatory cytokines are also elevated when high energy diets are fed (Eckel and Ametaj, 2016). Therefore, nutritional levels, both over and under appropriate nutrition, will influence SP components and potentially reproductive efficiency.
Cytokines are a diverse group of signaling proteins involved in a multitude of immunological functions (Chen et al., 2018). Cytokines are known to be pleiotropic and synergistic, acting in cascading pathways to create a strong biological effect in a given tissue (Kany et al., 2019). Recently, SP was found to target female tissues in humans, mice, and cattle to elicit an immunological response needed for pregnancy establishment (Robertson et al., 2009;Bromfield, 2014). Specifically, cytokine signaling activates structural modifications in female tissues through the recruitment of leukocytes, macrophages, dendritic cells (Robertson et al., 1997), and regulatory T cells (Tsuda et al., 2019). The cytokines within SP cause cascading signals in the endometrium to promote the inflammatory immune response that will in turn facilitate pregnancy establishment, embryo tolerance, and fetal development in many species including in humans, mice, and cattle (Tremellen et al., 1998;Robertson, 2005). Furthermore, complex cytokine networks have important roles in a wide range of reproductive processes such as maternalfetal interaction by interferon-t, uterine expansion by Interleukin (IL)-1b and IL-8 (Orsi and Tribe, 2008), and cervical remodeling by IL-10 to prepare for parturition (Van Engelen et al., 2009). At copulation, SP stimulates the release of IL-10 in the female reproductive tract, potentially for the immunological balance needed to avoid rejection of the spermatozoa and subsequent embryo (Denison et al., 1999). The balance of pro-and anti-inflammatory cytokines is crucial for the survival of the spermatozoa in the female tract and the tolerance of the fetus. Since nutrition may impact the overall inflammatory state of SP and immune responses within the female after intercourse, it is crucial to understand how different diets affect cytokines within the SP. Therefore, the hypothesis of the current study is that differing levels of nutrition and body condition score (BCS) will alter cytokine concentrations within the SP of beef bulls.

Materials and methods
All experimental procedures were approved by the University of Tennessee Institutional Animal Care and Use Committee.

Experimental design and sample collection
The experimental design and sample collection (Figure 1) was previously reported in detail (Harrison et al., 2022). Briefly, mature Angus bulls (n = 11; Body Weight [BW] = 738 kg; BCS = 6; Age = 4 years) were purchased from a single farm in South Dakota, individually housed (2.44-m by 12.19-m paddock), and provided ad-libitum water, mineral supplement (CO-OP Supreme Cattle Mineral; Tennessee Farmers Cooperative; Lavergne, TN), and a single ration at varying levels to achieve targeted BW and BCS goals. The diet consisted of 35% grass hay, 35% cracked corn, 20% dried distillers grains, and 10% soybean meal resulting in a diet that was 80.3% TDN, 1.92 Mcal/kg NEm, 1.28 Mcal/kg NEg, and 19% crude protein. Prior to treatment assignment, all bulls were required to pass a breeding soundness exam (BSE). Pre-treatment average scrotal circumference was 38.52 ± 1.45 cm (34 cm minimum), forward progressive motility averaged 45.8 ± 7.3% (30% minimum), morphological defects were 19.9 ± 5.9%, 5.9 ± 3.9%, 0.64 ± 1.02%, and 26.5 ± 6.1% for head, mid piece, tail, and total (70% minimum), respectively (Harrison et al., 2022). The experimental design consisted of two nutritional pathways (over-fed [n = 5] and restricted [n = 6]) and four nutritional planes (gain, loss, maintenance at a suboptimal BCS, and maintenance at the ideal BCS) to achieve eight individual treatments. The treatments included: gain to obesity (GO), maintenance at an obese BCS (BCS = 8; MO), loss after obesity (LO), maintenance at ideal adiposity (BCS = 6) after obesity (IMO), loss from nutrient restriction (LR), maintenance at a low BCS (BCS = 4; MR), gain after nutrient restriction (GR), and maintenance at ideal adiposity after nutrient restriction (IMR). Individual non-shrunk BW and BCS (1 = emaciated and 9 = obese, (Wagner et al., 1988)) was taken every 14 d to monitor changes within bulls. Semen was collected every 84 d for cytokine concentrations. The ejaculate was centrifuged in the lab for 5 min at 2500 × G to pellet spermatozoa at the bottom of the tube. The SP was then aliquoted into 2 mL microcentrifuge tubes and stored at -80°C until cytokine analyses could be conducted.

Statistical analyses
A complete randomized design was implemented in GLIMMIX procedure of SAS 9.4 (SAS Institute; Cary, NC) to determine the effects of nutritional treatments on cytokine concentrations. The experimental unit was individual bull. The fixed effect was treatment (i.e. nutritional period) and the random effect was bull within treatment. Initial cytokine concentrations prior to the onset of treatments as well as the volume of each ejaculate were included as covariates for cytokine analyses but were removed from the model if found to be nonsignificant (P > 0.05). Normality was assessed by the Shapiro-Wilk statistic > 0.8 and the Kolmogorov-Smirnov test > 0.2. Cytokine concentrations that were not normally distributed were log transformed to achieve normality and were reported as back transformed means and standard errors. Interleukin-4, IL-6, IP-10, and MCP-1 were not able to achieve normality with or without transformation and were removed from further concentration analysis. Means were considered different when P ≤ 0.05 with a tendency at P ≤ 0.10.
For cytokine profile analyses, log-transformed cytokine concentrations were analyzed using MetaboAnalyst 5.0 (Pang et al., 2021) to identify differences in cytokine profiles amongst nutritional treatments. The chemometrics analysis using partial least squares discriminant analysis (PLS-DA) was utilized to evaluate nutritional treatments. Variable importance in projections (VIP), known as a weighted sum of squares of the partial least square loadings, was also performed on cytokines within PLS-DA to determine which cytokines were most influential in the resulting profile.

Results
Nutritional plane effects on BW, BCS, and gross fertility (i.e. motility and morphology) measurements have been previously reported (Harrison et al., 2022). Initial cytokine concentrations, from the sample prior to the onset of treatment, was included in the model as a covariate but only remained in the models, due to a significant impact (P ≤ 0.04), on IFN-g and TNF-a. The volume of the ejaculate, for each bull, was also included as a covariate; however, only found to be impactful (P = 0.04) for the analysis of IFN-g. All cytokine concentrations returned to values similar to initial concentrations during the IMO and IMR treatments with the exception of MIP-1a and VEGF-A. Nutritional treatments influenced (P < 0.05) IFN-g, IL-8, MIP-1a, MIP-1b, TNF-a, IL-1b and VEGF-A. Interferon-g concentrations were influenced by nutritional treatments (P < 0.01), with concentrations decreasing the greatest amount during the IMO treatment with concentrations similar to the IMR treatment. Both IMO and IMR concentrations were similar to concentrations prior to the onset of treatments. The greatest concentration occurred during the MO treatment when bulls on the over-fed pathway were maintained at a BCS of 8 ( Figure 2A). Nutritional treatments also affected (P < 0.0001) IL-8 ( Figure 2B) and MIP-1a (P < 0.0001; Figure 2C) similarly with the IMO and IMR treatments having the greatest concentrations and the lowest concentrations during the GR treatment. However, IL-8 concentrations returned to levels similar to initial cytokine concentrations during the IMO and IMR treatments. Conversely, MIP-1a did not return to concentrations observed at the onset of treatments. Concentrations of MIP-1b were greatest (P < 0.01) during the IMR treatment with similar concentrations for the IMO treatment, and both treatments not different from concentrations prior to the onset of treatments. The most reduced concentrations for MIP-1b occurred during the GR and LO treatments which is the dietary restriction period of both pathways ( Figure 2D). The pro-inflammatory cytokines, TNF-a (P < 0.0001; Figure 2E) and IL-1b (P = 0.05; Figure 2F), followed similar trends with the greatest concentrations occurring during the IMO and IMR treatments, of which both were not different from initial concentrations. Tumor necrosis factor-a decreased in concentration but was not different from the remaining treatments. However, IL-1b had the lowest concentrations during the LO, MR and GR treatments. The angiogenic cytokine, VEGF-A, was influenced by nutritional treatments (P ≤ 0.02) and had the greatest increase in concentration during the IMO and MO treatments when bulls on the over-fed pathway were on diets which supported maintenance of BW and BCS. Similar concentrations did occur during the IMR and LO treatments; however, VEGF-A did not return to levels observed prior to the onset of nutritional treatments. The greatest decrease in VEGF-A concentrations occurred during the LR treatment which was similar to concentrations during the GO treatment ( Figure 2G).
To establish cytokine profiles for each nutritional treatment, a PLS-DA was created within MetaboAnalyst 5.0. Nutritional treatments appeared to be impactful on SP cytokine profiles with overlap trends occurring between some groups ( Figure 3A). Specifically, the IMO, IMR and initial nutritional treatments which were similar via overlapping ellipses but completely separate ellipses from the GR, LO, LR, MO and MR treatments. The only profile that stretched over both of the groups was the SP profile for the GO treatment. Variable importance in the projection scores were generated to determine the cytokine levels of contribution to these profiles during each nutritional period. All of the cytokines were included within the PLS-DA but the cytokines with < 0.5 VIP score included: IL-6, MCP-1, VEGF-A, IL-4, IL-8, MIP-1b and IL-10. These cytokines were deemed relatively non-impactful to the overall profile of the SP ( Figure 3B). Cytokines identified as moderately impactful (VIP score = > 0.5 and < 1.0) included: IP-10, IL-17a, and IL-1a ( Figure 3B) and cytokines, (TNF-a, IFNg, MIP-1a, IL-1b and IL-36RA) with a VIP scores > 1, were highly influential on the cytokine profiles of each nutritional period ( Figure 3B). Tumor necrosis factor-a had the greatest impact (VIP score > 2) on the SP profile among all nutritional treatments and was greatest during the initial period prior to treatments and the lowest during the MO treatment ( Figure 3B). In contrast, IFN-g (VIP score > 1.5; Figure 3B) had the greatest impact during the MR treatment and the lowest during the IMO treatment. Macrophage inflammatory protein-1a (VIP score > 1.5; Figure 3B) had the most influence during the IMO treatment and the lowest impact during the GR treatment. The cytokines, IL-1b and IL-36RA, followed a similar prevalence pattern with greatest impacts (VIP score = 1.5) during the initial period and the least during the LO treatment for IL-1b and the MR treatment for IL-36RA; however, IL-1b had a greater overall impact ( Figure 3B).

Discussion
Sire management, through the control of nutrition and other environmental factors, is imperative to ensure optimal fertility and advance overall herd genetics. Nutritional effects have been reported to impact the ejaculate in boars, rams, and bulls (Brown, 1994). Nutrition can also play a role in influencing the components of SP, specifically cytokines (Eckel and Ametaj, 2016). Cytokines rarely act individually, more so as a network of highly influential, cascading protein molecules to control local and systemic inflammation (Dinarello, 1989). Many cytokines and chemokines have been reported to be altered by seminal plasma in the endometrium of cows (Ibrahim et al., 2019) including IL-8 and IL-1b. The current study found IL-8, IL-1b, IFN-g, MIP-1a, MIP-1b, TNF-a, and VEGF-A, to be impacted by nutritional levels and may be indicative of a role in reproductive success.
Interferon-g is known as a pleiotropic cytokine that can have both pro-and anti-inflammatory effects, and has been found in the uterus during early pregnancy (Murphy et al., 2009). Within the maternal endometrium, IFN-g is abundantly produced by uterine natural killer cells as well as trophoblasts to initiate endometrial vasculature remodeling, angiogenesis, and the maintenance of the placenta in mice (Murphy et al., 2009). Interestingly, IFN-g in SP has been linked to female infertility in humans outside of sperm and male parameters (Robertson et al., 2003). Concentrations of IFN-g in the current dataset were greatest during periods of nutritional perturbation. This may be indicative of the nutritional stress influencing IFN-g concentrations in bovine SP which could impact the ability to produce viable embryos and establish a pregnancy. In support of these results, the immune-regulatory concentrations of IFN-g in humans have generally been detected at low levels within SP; however, was substantially elevated when a disease or infection was present (Leutscher et al., 2005;Vanpouille et al., 2016). Indicating that fluctuations in cytokine concentrations occur in response to external and internal factors as the body attempts to return back to homeostasis throughout pregnancy. Tumor necrosis factor-a is an acute phase cytokine responsible for mediating acute inflammatory reactions to diseases and infection sites (Burger and Dayer, 2002). The pro-inflammatory properties of TNF-a are known to hinder sperm motility and functional capacity by increasing chemotactic activity and inducing the immuno-cascade effects of neutrophils (Hill et al., 1989). In contrast, TNF-a has been found to be essential for early pregnancy establishment (Toder et al., 2003). Yet the overproduction of TNF-a may cause early embryonic loss or implantation failure in humans (Saito et al., 2010;Alijotas-Reig et al., 2017). Bovine TNF-a also incites inflammation by promoting neutrophil recruitment when induced by lipopolysaccharides (Sohn et al., 2007). Interleukin-1b has similar effects on semen quality and establishment of pregnancy as TNF-a (Gruschwitz et al., 1996). Interluekin-1b is also responsible for activating the innate immune response by mediating recruitment and activation of macrophages and neutrophils (Ott et al., 2007). The highly correlated pro-inflammatory cytokines, TNF-a and IL-1b (Eggert-Kruse et al., 2007), had similar cytokine concentration trends in the current study, but were also found to be impactful on overall cytokine profiles during each nutritional period. The current study identified TNF-a to have the greatest impact (VIP score) on cytokine profiles within SP during different nutritional periods. These closely associated and impactful pro-inflammatory cytokines are potentially needed and, thus, play important roles in pregnancy establishment.
Chemokines are a family of small cytokines, mostly known for their role in chemotaxis (Hughes and Nibbs, 2018), inflammation, immune surveillance, and angiogenesis (Dimberg, 2010). Many chemokines are known to be pro-inflammatory, while others are thought to control cell migration for normal tissue growth and maintenance (Hughes and Nibbs, 2018). Interleukin-8 is a potent chemokine involved in leukocyte migration and cell activation for events associated with inflammation (Nederlof et al., 2017). Interleukin-8 decreased substantially during the GR treatment in comparison to the greatest concentration following the IMO and IMR treatments. In humans, IL-8 concentrations have been reported to be present in high concentrations within SP of healthy fertile men (Politch et al., 2007), and, thus, may indicate excessive BW gain in bulls is detrimental to fertility. In contrast, IL-8 has also been known to dramatically increase in response to bacterial and viral infections which may negatively affect spermatozoa resulting in poor sperm motility within the ejaculate (Eggert-A B

FIGURE 3
Analysis of cytokine profiles based on nutritional treaments: (A) Partial least squares discriminant analysis visualizing differences in inflammatory cytokines during differing nutritional treatments and (B) Variable importance in the projection for individual cytokines plot. A single ellipse represents a 95% confidence interval within that nutritional treatment. Kruse et al., 2001). Therefore, a certain level of this immunostimulatory chemokine from nutritional stress could result in potential pregnancy failures and diminished herd outcomes. The pro-inflammatory chemokines, MIP-1a and MIP-1b, are known as chemoattractants of monocytes and macrophages which contribute to the regulation of uterine macrophage populations of mice . However, MIP-1a and MIP-1b have not been extensively researched in cattle. The current study identified MIP-1a to have an influential impact (VIP score) on the ejaculate during different nutritional periods; however, MIP-1b was deemed relatively non-impactful on the cytokine profiles within SP. Previous research has reported increases in MIP-1a and MIP-1b concentrations are indicators of infection or disease (Garzino-Demo et al., 1999;O'Grady et al., 1999;Chaisavaneeyakorn et al., 2003). Macrophage inflammatory protein-1a demonstrated significant concentrations within decidual tissue increases within humans who underwent abortions (Freis et al., 2018) and overall pregnancy losses within the first 10 weeks of gestation (Park et al., 2010). On the male side, humans with HIV expressed lower concentrations of MIP-1a and MIP-1b released by cytotoxic T cells (Cocchi et al., 2000). Thus, concentrations of MIP-1a and MIP-1b may increase acutely but long-term stress from diseases, like HIV (Cocchi et al., 2000), or nutritional stress (current dataset) results in lower concentrations altering the impacts of these chemokines on physiological functions. The increase of the pro-inflammatory chemokines MIP-1a and MIP-1b within SP may be detrimental to pregnancy outcomes due to early embryo loss.
Immuno-suppressive cytokines, IL-10 and IL-36RA, play important roles in inhibiting the synthesis of pro-inflammatory cytokines (Zhang and An, 2007). More specifically to reproduction, anti-inflammatory cytokines create the overall immuno-suppressive state of the mucosal environment to aid in sperm survival within the oviduct (Torres-Poveda et al., 2014). In humans, concentrations of IL-10 were greater in healthy male individuals demonstrating the natural immuno-tolerant environment for sperm of the female reproductive tract (Camejo, 2003). The immuno-suppressive functions of IL-10 also prevent the rejection of the semi-allogenic fetus (Chatterjee et al., 2014) and promote conceptus attachment to the uterine endometrium in dairy and beef cattle (Odhiambo et al., 2009). Our study demonstrated IL-10 to be non-influential to the overall ejaculate and the concentrations were not influenced by differing periods of nutritional stress. This potentially indicates that the desired immuno-tolerant levels for sperm survival can be maintained even under nutritional stress. The antagonistic effects of IL-36RA, a member of the IL-1 super family, inhibits inflammation by inhibiting IL-36 from binding to IL-36R (Yi et al., 2016). The antagonistic effects from IL-36RA impede signals between the toll/ interleukin-1 receptor domains which inhibits NF-kB signaling cascade (Murrieta-Coxca et al., 2019). The roles of the IL-36 cytokine family remain mostly unknown in pregnancy; however, the immuno-suppressive effects of IL-36RA could promote a balanced uterine environment for successful pregnancy establishment. Our findings quantified IL-36RA to be influential on the cytokine profiles within SP under different nutritional periods; however, concentrations of IL-36RA were found to not be influenced by the nutritional treatments. This may demonstrate that IL-36RA has a more important role in creating the immuno-suppressive environment compared to IL-10 to increase sperm survival and to promote embryo establishment.
Angiogenic cytokines can induce endothelial cell activation and proliferation for the creation of new blood vessels from precursor cells (Ucuzian et al., 2010). Proper function of VEGF-A is critical during every step of placental growth and vascular formation to provide enough blood and nutrients to the growing fetus (Chen and Zheng, 2014). Thus, VEGF-A being present in the greatest concentrations within SP is not entirely surprising. Furthermore, research has demonstrated a decrease of VEGF-A in the SP of obese humans causing a negative effect on semen quality (Zhu et al., 2021). Our study agrees with these results since VEGF-A had one of the lowest VIP impacts of all the cytokines presented, more specifically during the LR and GO treatment. Angiogenic cytokine VEGF-A did not play a significant role in any of the cytokine profiles yet had the greatest overall concentrations of all cytokines evaluated. The drastically higher concentrations of VEGF-A in SP may indicate a role in angiogenesis in the developing fetus, causing the overall impact of this cytokine to be low. These concentrations decreased during the LR and GO treatments when bulls deviated from their initial basal BCS, demonstrating the potential decrease in semen quality which may be a cause of pregnancy failures.
The focus of the current study was to elucidate the impacts of nutrition. However, the authors admit that it is possible that some of the results could have been, in part, influenced by photoperiod. Most mammalian immune and endocrine function have been reported to be impacted by photoperiod to varying extents (Reiter, 1983;Dahl et al., 2002). The current study did take place over the course of an entire year and, thus, the bulls were exposed to different photoperiods. Scrotal circumference and spermatogenesis measurements were altered by nutrition but not impacted by time of year in which the sample was taken (Harrison et al., 2022). The cytokines quantified in the current dataset, MIP-1a, IL-1b, IL-8, and TNF-a, that were impacted by nutritional treatment responded, in varying extents, to nutritional perturbations similarly compared with ideal nutrition and BCS. The exceptions to this were the changes observed in IFN-g, MIP-1b, and VEGF-A concentrations but these changes also did not correlate with seasonal changes in photoperiod. Therefore, photoperiod does not appear to impact cytokines in seminal plasma but more direct research studies would need to be conducted to completely elucidate the photoperiodic effect.
In conclusion, cytokines and chemokines fluctuated between the nutritional periods, potentially in response to nutritional stress. The cytokines, TNF-a, IFN-g, and MIP-1a had the greatest impact on the overall profile of SP during nutritional periods. Moreover, this dataset demonstrated a similarity between cytokine profiles when the animal gains or loses BW or is maintained at an abnormal BCS. The cytokines: IFN-g, MIP-1b, IL-1b, IL-8, and TNF-a, during the IMO and IMR treatments, returned to the initial concentrations and profiles. Therefore, varying nutritional levels can influence the immunological substrates within SP of mature bulls which could potentially affect the sire's ability to successfully establish pregnancy. However, these changes are not permanent and will return under correct nutritional management. Further studies are still required to fully understand the long-term impacts of nutrition on SP, the uterine responses to differing cytokine profiles, and how these cytokines influence the establishment of pregnancy all of which could alter reproductive efficiency and success.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement
The animal study was reviewed and approved by University of Tennessee Animal Care and Use Committee 2713-0819.

Author contributions
TH conducted the study, collected all samples, analyzed the data, and was the primary author for the manuscript. EC, KB, and TA-S assisted in sample collection, data analyses, and provided editorial comments to the manuscript. RP assisted in cytokine analysis and provided editorial comments to the manuscript.
LGSc aided in study design, statistical analyses, and provided editorial comments.
LGSt, and FS provided study guidance and editorial comments. KM designed and supervised the study and data management as well as provided editorial comments to the manuscript. All authors reviewed and approved the final version of the manuscript.

Funding
This work was funded by the University of Tennessee Institute of Agriculture, University of Tennessee Department of Animal Science, and the United States Department of Agriculture, National Institute of Food and Agriculture, Capacity Grant No. 1019048.