Biomarkers Correlate With Body Composition and Performance Changes Throughout the Season in Women's Division I Collegiate Soccer Players

Introduction

The physiological demands of soccer require athletes to possess high levels of aerobic capacity, as well as muscular strength, muscular power, speed, and speed-endurance (Kraemer et al., 2004; Stolen et al., 2005). Collegiate soccer players face unique challenges to maintaining these metrics throughout the competitive season. Specifically, multiple games per week, frequent practices, and travel requirements can increase these athletes' overall stress and recovery needs (Walker et al., 2019b). Moreover, changes in diet, environment, and sleep combined with the burdens of academic requirements that often accompany collegiate life can increase athletes' risk for injury and further hinder the ability to enhance aspects of sport performance that are required to succeed (Silva, 1990; Meeusen et al., 2013; Mann et al., 2015). Therefore, it is important to employ athlete monitoring techniques that assess the cumulative impact of these stressors, not simply the physical stress due to the sport, in order to better facilitate adequate recovery and optimal performance. Monitoring techniques that track changes in an athlete's physiological state may prove beneficial in determining an athlete's readiness to perform (Huggins et al., 2019; Walker et al., 2019b).

The implementation of systematic performance assessments is one important method of athlete monitoring that plays a key role in the determination of athlete readiness. Systematic assessments, more specifically, periodic maximal performance tests, are essential in determining an individual athlete's capabilities, as well as monitoring the effectiveness of training (Krustrup et al., 2005; Stolen et al., 2005; Sporis et al., 2009). For power-endurance athletes, these tests may include maximal oxygen consumption (VO_2max) tests that measures the highest flow of oxygen that can be used by the tissues in order to determine aerobic capacity (Taylor et al., 1955), maximal vertical jump (VJ) height to assess muscular power, and repetition maximum (RM) tests for muscular strength (Walker et al., 2019a). While maximal performance tests are optimal to determine player development, the potential for added fatigue incurred by the athletes as well as the time burden required to evaluate an entire team, makes these tests difficult to implement during the season in a team sport environment (Bangsbo and Lindquist, 1992; Twist and Highton, 2013; Halson, 2014; Thorpe et al., 2015). In a condensed collegiate season marked by high training volumes, adding an additional load of performance assessments periodically throughout the season may not be feasible, as coaches and support staff would be required to concede already limited sport specific training/practice times to implement testing. Therefore, it is important to incorporate other assessment methods that do not incur additional physical loads on the athletes, can be implemented without burdening coaches, and allow for continuous monitoring throughout the season (Thorpe et al., 2015). The use of workload tracking systems, such as heart rate (HR) and Global Navigation Satellite System (GNSS) metrics (or global positioning satellite (GPS) systems specific to the United States), throughout the season plays an important role in assessing recovery needs, detecting non-functional overreaching, reducing risk of injury, and ultimately maximizing athletic performance (Halson, 2014). While these workload tracking systems can detect an athlete's response to on-field training stressors, the use of biomarker monitoring offers the ability to track an individual player's physiological response to overall stress as well as to identify the balance between training and recovery (Lee et al., 2017; Pedlar et al., 2019).

Periodic biomarker assessments during the competitive soccer season offer an accurate, objective method for evaluating training related stress and recovery needs as well as overall athlete health (Bessa et al., 2016; Pedlar et al., 2019). The use of biomarkers may better indicate exercise related stress independent of other factors, such as motivation, which can confound the results of performance tests or other subjective measures (Bessa et al., 2016). Although there is no single definitive biomarker that can be used to monitor training progress or recovery status in team sport athletes (Marques-Jimenez et al., 2016; Lee et al., 2017), measuring biomarkers that are indicative of an athlete's physiological state may prove to be a useful tool in determining athlete readiness. Kraemer et al. (2004), proposed that if the physical demands of training are too high or there is insufficient recovery to meet the training demands, a catabolic state will persist throughout the season, resulting in impaired performance. Markers associated with catabolism, such as total and free cortisol (TCORT and FCORT), creatine kinase (CK), c-reactive protein (CRP), and interleukin-6 (IL-6), are often released in response to excessive strain or muscle damage (Lee et al., 2017). Elevations in these catabolic markers have been shown to promote inflammation and muscle breakdown and chronically high values may indicate impaired recovery (Lee et al., 2017) thus affecting performance outcomes (Kraemer et al., 2004). However, if training demands are met with adequate recovery, an anabolic environment will result in the maintenance or improvement of performance (Kraemer et al., 2004). For example, the anabolic hormone testosterone, promotes an anticatabolic effect and may improve athletic performance through direct effects on muscle growth, enhanced muscle repair after exercise, accelerated neuromuscular coordination, and the regulation of body fat (Koundourakis and Margioris, 2019). Other anabolic hormones relevant to the female athlete such as, insulin-like growth factor-1 (IGF-1) and growth hormone (GH), have been shown to promote muscle growth and development of muscle tissue (Frystyk, 2010; Nindl and Pierce, 2010; Lee et al., 2017) and may provide additional insights in regards to muscular adaptations to training. Estradiol (E₂) is a form of the hormone estrogen, which is responsible for female sex characteristics and has been shown to have various beneficial effects on athletic performance through anti-inflammatory properties and improved muscle recovery and repair rates (Tiidus, 2003; Koundourakis and Margioris, 2019). Estrogen may exert a protective effect on muscle tissue during exercise by acting as an antioxidant, membrane stabilizer, and through receptor substrate characteristics (Cano Sokoloff et al., 2016; Koundourakis and Margioris, 2019). Additionally, alterations in E₂ have been shown to be reflective of an energy deficiency and suggest hypothalamic-pituitary-gonadal (HPG) axis dysfunction in the female athlete (Dueck et al., 1996; Vanheest et al., 2014). Overall, chronic elevations in anabolic hormones may lead to improvements in body composition and performance, whereas declines in these markers may reflect impaired potential for adaptations to training (Lee et al., 2017).

Changes in hormonal concentrations in conjunction with changes in exercise performance during the course of a soccer training program have been previously reported in male athletes (Filaire et al., 2003; Kraemer et al., 2004; Silva et al., 2008; Huggins et al., 2019). Over the course of a competitive season, Kraemer et al. (2004) demonstrated decrements in sprint speed, vertical jump height, and knee extensor strength. These downturns in performance coincided with elevated concentrations of circulating cortisol and reduced testosterone. Together these findings indicate a predominance of catabolic processes occurred during the season (Kraemer et al., 2004). In another study in professional soccer players, testosterone was found to have a direct relationship between vertical jump height and power output suggesting its beneficial effects on athletic performance, whereas a negative correlation was seen between endurance capacity and both serum cortisol and testosterone levels (Bosco et al., 1996; Koundourakis and Margioris, 2019). Silva et al. (2008) demonstrated increases in average team CK levels coincided with diminished team performance (win percentage) as well as stagnation in aerobic endurance (Silva et al., 2008). Similarly, Filaire et al. (2003) examined changes in various biomarkers throughout a season and found that decreases in glutamine and increases in cortisol overlapped with a decreased win percentage; however, immunological and hematological markers, VO_2max and body composition remained unchanged (Filaire et al., 2003). In a comprehensive study by Huggins et al. (2019), 92 separate biomarkers were analyzed at five time points throughout a collegiate men's competitive soccer season. Although fluctuations in biomarkers were observed, these fluctuations were interpreted to be typical changes that occur over the course of a season within a team setting (Huggins et al., 2019). Furthermore, the team experienced little change in player load between training blocks which coincided with slight decreases in GPS and accelerometry metrics as well as no significant changes in body composition, VO_2max, or HR at anaerobic threshold. In one of the only studies evaluating biomarker changes in female soccer athletes, researchers demonstrated significant declines in VO_2max, vertical jump, body weight and body fat percentage from pre- to post-season which occurred along with increases in various catabolic biomarkers including cortisol, IL-6, and CK (Walker et al., 2019b). The decrements in performance metrics in conjunction with sustained increases in catabolic biomarkers from preseason may be a result of the accumulated stress of the season thereby negatively impacting performance variables.

Although it is apparent that changes in biomarkers may occur during a female soccer season along with performance decrements, it is unclear if changes in specific anabolic and catabolic markers correlate to changes in performance outcomes. In addition, identifying and utilizing athlete-monitoring methods, such as biomarkers, that are associated with performance outcomes may allow for prompt player evaluation and may also be useful in designing individualized recovery strategies (de Hoya et al., 2016; Pedlar et al., 2019; Walker et al., 2019b). Therefore, the purpose of this study was to determine the impact of a National Collegiate Athletic Association (NCAA) Division I (DI) women's soccer season on anabolic and catabolic biomarkers in conjunction with various measures of fitness and performance. Additionally, this study sought to identify the relationship between biomarkers, performance, and body composition changes in female collegiate soccer athletes over the course of a competitive season. It was hypothesized that as a team, the athletes would experience increased catabolic markers and declines in anabolic markers throughout the season, and decrements in performance metrics would occur from pre- to post-season. However, when assessing individual changes in concentrations of biomarkers from baseline, greater change in catabolic markers would negatively correlate with performance outcomes and fat free mass (FFM), while changes in anabolic markers would be positively associated with performance and FFM.

Materials and Methods

Participants

Women's DI college soccer players (N = 21; Age = 20 ± 2 y; Weight = 66.3 ± 6.2 kg, Height = 169 ± 7 cm) participated in monitoring and assessments throughout the course of a competitive soccer season. The team consisted of six defenders, six forwards, seven midfielders, and two goal keepers. The 16-week season, which began in early August and finished in mid-November, corresponded with the end of summer period and lasted until mid-fall in the Northeast region of the United States. Athletes performed training and performance testing as part of regular team activities associated with their sports science program. During this time, athletes participated in a periodized strength and conditioning (S&C) program one time per week set by the team's strength coach in addition to on-field agility training sessions incorporated into practice. The strength training program focused on power and strength development with sessions consisting of 1–2 power movements followed by 3–5 strength movements and associated core work. Power movements consisted of exercises such as the high pull and push press. Strength movements consisted of exercises including back squat, split squat, bench press, incline dumbbell press, Romanian deadlifts, rack pulls, single arm dumbbell rows, and associated variations of these movements. Weekly loads were programmed based on preseason 3 RM strength testing results and consisted of 70–80% of their estimated 1RM (Brzycki, 1993). All subjects received clearance by the University's Sports Medicine staff prior to testing and at the start of the season. This research was approved, and written consent was waived by the Rutgers University Institutional Review Board for the Protection of Human Subjects (IRB#16-050M). All procedures performed were in accordance with the 1964 Declaration of Helsinki and its later amendments or comparable ethical standard.

Experimental Design and Procedures

This study sought to evaluate the effects of a competitive collegiate soccer season on biomarkers and performance metrics in order to determine the correlation between biomarkers, body composition, and performance outcomes in DI female soccer players. Biomarkers were analyzed prior to the start of preseason and every 4-weeks following the initial analysis. Maximal performance testing occurred prior to the start of the season and within 1-week following the final match. Body composition testing occurred prior to the start of preseason, at weeks 6, 10, 14, and within 1-week following the final match. See Figure 1 for a timeline of assessments.

FIGURE 1

Figure 1. Timeline of assessments. BC, Body Composition Testing; BD, Blood Draws. Strip boxes represents 2-weeks of preseason, dark gray boxes represent weeks of the regular season, black boxes represent tournament time, light gray boxes represent time between blood draws, white boxes represent time between body composition testing.

Performance Testing

Athletes reported to the IFNH Center for Health and Human Performance (CHHP) prior to the start of preseason (W0) (late July) and again within 1-week following the final competitive match (W17) (mid-November) to complete a battery of fitness tests to assess body composition, VO_2max, and muscular strength during a 1-week period. Subjects were instructed to arrive euhydrated, 2–3 h fasted, having abstained from caffeine, and having performed no exercise 24 h prior to testing. Body composition was assessed by air displacement plethysmography (Dempster and Aitkens, 1995) (BOD POD, COSMED, Concord, CA, USA) to determine percent body fat (%BF), and fat free mass (FFM) using the Brozek formula (Brozek et al., 1963). In addition to W0 and W17, body composition was also assessed at weeks 6, 10, and 14. Next, following a 10–15 min generalized dynamic warm up, subjects performed maximal countermovement vertical jumps with both hands-on-hips (VJ_HOH) and with arm swing (VJ) that was assessed using the Just Jump Mat (Probotics, Huntsville, AL, USA). Subject were given three attempts for each jump variation, with the highest jump for each recorded. The error of vertical jump height for females has been reported with a standard error of the mean of 1.7 cm and a CV of 4.4% (Nuzzo et al., 2011). After completing power testing, a treadmill-based maximal graded exercise test (GXT) was used to measure maximal aerobic capacity (VO_2max) via direct gas exchange using an indirect calorimeter (Quark CPET, COSMED, Concord, CA, USA). A speed-based protocol was used with stages that were MET-equated to the Bruce protocol. This protocol included 2 min stages at a constant 2% incline. The speeds were as follows: 6.4, 7.9, 10.0, 11.7, 13.7, 15.6, 17.1, 18.2, 19.8, 21.1 (km/h) (McFadden et al., 2020). Subjects continued the test with encouragement from lab staff until volitional fatigue. At least three of the following criteria were met verifying attainment of VO_2max: a leveling off or plateauing of VO₂ with an increase in workload, attainment of age predicted heart rate max, a respiratory exchange ratio >1.10, and/or an RPE ≥18 (Thompson et al., 2000). Heart rate was continuously monitored using a Polar S610 heart rate monitor (Polar Electro Co., Woodbury, NY, USA) to obtain maximal heart rate (HR_max). One player (n = 1) did not participate in post season VO_2max testing due to clearance from sport medicine.

On a separate day, maximal strength testing was performed within a group setting under the instruction and supervision of the team's strength and conditioning staff. Following a generalized dynamic warmup, 3RM tests were completed for squat (SQ), bench press (BP), and deadlift (DL) in that respective order and in accordance with NSCA guidelines (Haff and Triplett, 2016). Four players did not participate in post-season strength testing due to professional team opportunities (n = 1), job requirements (n = 1), or clearance from sports medicine (n = 2).

Biomarker Collection and Analysis

Athletes were instructed to report to the CHHP in a euhydrated condition immediately prior to the start of pre-season for blood draws (BD1), and every 4-weeks following (BD2, BD3, and BD4). All blood draws were taken at the CHHP between 0700 and 0830 h following an overnight fast, and BD2, BD3, and BD4 draws occurred ~18–36 h following a game (Walker et al., 2019b). Blood samples were obtained (antecubital fossa, 21G, BD Vacutainer, Safety-Lok) from the antecubital vein in the arm by an experienced phlebotomist while subjects were seated. Whole blood was collected and left to clot at room temperature. Samples were drawn into collection tubes consisting of anticoagulant (EDTA), or clot activator (gel free and Serum Separator Tubes). Following collection, blood was centrifuged for 10 min at 4,750 rpm (Allegra x-15R, Beckman Coulter, Brea, CA, USA). Plasma and serum samples were shipped (in containers designed to maintain approximately 20, 4, or −20°C depending on the analyte) to a processing facility (Quest Diagnostics, Inc., San Juan Capistrano, CA, USA) for analysis. Samples were run in duplicate and the coefficient of variation (CV) was between 0.5 and 10.0% for all biomarkers. Biomarkers analyzed included CK, CRP, TCORT, FCORT, E₂, GH, IGF-1, IL-6, free testosterone (FTEST), and total testosterone (TTEST).

Training Load

Athlete workloads were monitored at all practices, games, and strength sessions throughout the entirety of the 16-week season which consisted of 21 games, 63 practice sessions, and 15 strength sessions. Workload was monitored through HR and GPS tracking technology using the Polar TeamPro system (Polar Electro Co., Woodbury, NY, USA). Physiological attributes of the player obtained from laboratory testing (age, height, body weight, sex, VO_2max, HR_max, and ventilatory threshold) were used to program each individual player's monitor. Training load (TL), exercise energy expenditure (EEE), and distance covered (DIS) were used to quantify the team's total weekly workload (Ceesay et al., 1989). TL was calculated via an algorithm developed by Polar^TM based on physiological attributes of the player, obtained from laboratory testing, and physical workload measures comprised of HR, GPS, and accelerometry information (Ceesay et al., 1989; Giersch et al., 2018). Weekly team workloads are presented in Figures 2–4.

FIGURE 2

Figure 2. Total weekly training loads.

FIGURE 3

Figure 3. Total weekly distances.

FIGURE 4

Figure 4. Total weekly exercise energy expenditures.

Statistical Analysis

Pre- and post-season performance and body composition data were analyzed using repeated measures (RM) multivariate analysis of variances (MANOVAS) with RM analysis of variance univariate follow-ups using a listwise deletion factor. RM MANOVAS were used to assess biomarker changes throughout the season. For each univariate analysis, the Huynh-Feldt epsilon was examined for the general model to evaluate sphericity. If the Huynh-Feldt epsilon exceeded 0.75, sphericity was considered to have been met, and the unadjusted statistic was used. If epsilon was <0.75, the adjusted Huynh-Feldt statistic was used to test significance. Planned simple contrasts were conducted using the baseline values as the comparison term. Pairwise contrasts were included in the case of significant univariate findings using the least significant difference method. Delta area under the curve (DAUC) was calculated for the biomarkers and body composition variables using the trapezoidal method to account for seasonal changes adjusted for baseline. Pearson-product moment correlations (r) were used to assess the relationships between biomarker changes and changes in performance and body composition. All analyses were conducted using SPSS (version 23, IBM, Armonk, NY, USA) with significance set at p < 0.05. Trends were considered to be p ≤ 0.10. Shapiro-Wilks test was used to test normality. In the case of non-normality of the data, variables were log transformed for analysis. An r value of 0.36–0.67 were considered to represent a moderate correlation. An r value of 0.68–1.0 were considered to represent a strong correlation (Weber and Lamb, 1970). Cohen's d was used to calculate effect sizes (ES) from baseline (P1, BC1, and BD1). Using Cohen's conventions, ES of 0.20, 0.50, and 0.80 were considered indicative of small, medium, and large effects, respectively. Biomarker and performance variables are presented and expressed as means ± standard deviation. Percent change (Δ) is represented by means and standard error. All individual raw biomarker and performance data is provided in the Supplemental Material.

Results

As a team, no changes in VO_2max, VJ, VJ_HOH, or DL performance were seen from W0 to W17 (p > 0.05). However, significant improvements in BP (ΔBP = 1.6 ± 0.7 kg; p = 0.02) and SQ (ΔSQ = 3.55 ± 1.08 kg; p = 0.01) performance were seen. Performance changes pre- to post-season are expressed in Table 1. No change in body weight was seen from W0 to W17 (p > 0.05) (Table 2). A significant decline in %BF was seen from W0 to W17 (Δ%BF = −1.4 ± 0.6%; p = 0.03) with a trend for an increase in FFM (ΔFFM = 0.59 ± 0.3 kg; p = 0.09) (Table 2).

TABLE 1

Table 1. Preseason to postseason performance changes.

TABLE 2

Table 2. Preseason to postseason body composition changes.

RM MANOVAS showed no significant time main effects with anabolic biomarkers (GH, IGF-1, TTEST, FTEST, E₂) over the course of the season (p > 0.05). However, significant time main effects were found for the catabolic biomarkers (TCORT, FCORT, CK, IL-6, CRP) (p = 0.001). Univariate follow-ups revealed significant changes in FCORT throughout the season (p = 0.01) (Changes from BD1 are presented in Table 3).

TABLE 3

Table 3. Biomarker Changes.

Biomarkers that correlate with performance and body composition variables are presented in Table 4. TCORT was negatively correlated with FFM (r = −0.48, p = 0.03). TCORT and FCORT also showed a positive correlation with %BF that approached significance (r = 0.39; r = 0.40; p = 0.08, respectively). TCORT was positively correlated with VO_2max (r = 0.47; p = 0.04) while a trend for a positive association with FCORT and VO_2max was seen (r = 0.39; p = 0.09). IL-6 was negatively correlated with BP (r = −0.53; p = 0.03) while IGF-1 and GH positively correlated with DL (r = 0.57; p = 0.03 and r = 0.59; p = 0.02, respectively). A trend for a positive correlation between CK and DL was seen (r = 0.49; p = 0.06). A positive correlation between TTEST and DL approached significance (r = 0.50; p = 0.07). No significant correlations were seen for biomarkers and VJ/VJ_HOH (p > 0.05). No significant correlations were seen between CRP, FTEST, and measures of body composition or performance (p > 0.05).

TABLE 4

Table 4. Biomarkers and performance/body composition correlations.

Within the biomarker analysis, TCORT positively correlated with FCORT (r = 0.63; p = 0.002). A negative correlation between TCORT and E₂ trended toward significance (r = −0.36; p = 0.10). GH was positively correlated with IGF-1 (r = 0.48; p = 0.03). IGF-1 positively correlated with FTEST (r = 0.60; p < 0.001) and E₂ (r = 0.49; p = 0.02). IL-6 positively correlated with CRP (r = 0.76; p < 0.001). The positive correlation between E₂ and FTEST approached significance (r = 0.45, p = 0.06). Correlations between biomarkers are presented in Table 5.

TABLE 5

Table 5. Biomarker correlations.

Within performance and body composition measures, a positive correlation was seen between body weight and BP (r = 0.56; p = 0.02). Trends were seen for a negative correlation between VO_2max and FFM (r = −0.37; p = 0.10) and a positive correlation with %BF (r = 0.42; p = 0.065). VJ_HOH positively correlated with VJ (r = 0.56; p = 0.01). Performance and body composition correlations are presented in Table 6.

TABLE 6

Table 6. Performance and body composition correlations.

Discussion

As a team, over the course of a 16-week soccer season, which was punctuated by high TL, EEE, and DIS, female collegiate soccer athletes experienced alterations in catabolic markers and a maintenance of anabolic markers. Additionally, improvements in body composition, increases in strength, and a preservation of both aerobic capacity and power metrics from pre- to post-season occurred in these athletes. After adjusting for individual baseline values, positive associations were seen between anabolic biomarkers (GH and IGF-1) and DL strength. Additionally, greater overall increases in markers of stress and inflammation (IL-6 and cortisol), were negatively associated with BP strength and FFM measures. Interestingly, net increases in TCORT from baseline were also positively correlated with VO_2max. These findings support a relationship between cumulative changes in anabolic and catabolic markers and changes in performance and body composition in female collegiate athletes throughout a competitive season, and therefore can be useful in the evaluation of athlete adaptation to training and readiness to perform.

In the current study, fluctuations were seen in catabolic markers at various timepoints indicating varying levels of stress occurred in these athletes throughout the season. The highest CK values occurred at the first two timepoints and declined as the season progressed. Though within typical female athlete reference ranges (Mougios, 2007), elevated CK levels at BD2 coincided with the highest TLs, EEE, and DIS of the season. Further, elevations in both IL-6 and CRP were apparent at BD2 which persisted before returning to baseline values at the end of the season. A lack of adequate recovery from soccer-related training stressors may result in a pro-inflammatory state exhibited by high concentrations of cytokines, such as IL-6 and CRP. IL-6 in particular, has been shown to increase in response to decreased muscle glycogen and muscle contraction (Pedersen and Febbraio, 2005), and is often associated with muscle damage (Bessa et al., 2016), diminished muscular strength, and physical function (Ferrucci et al., 2002; Stenholm et al., 2010). This may be evident in the negative correlation between IL-6 and BP in the current study. In addition to cytokines and CK, cortisol has been shown to be reflective of a catabolic state (Lee et al., 2017). In the current study, fluctuations in FCORT were seen throughout the season. Although no changes in TCORT values were noted, values at each time point were considered above the normal clinical reference ranges (127–569 nmol/L when assessed 0800–1,000 h) which indicates a catabolic state was evident in these athletes prior to the start of the season and remained as the season progressed. Chronic elevations in cortisol may make it difficult for an athlete to build or maintain muscle or recover from training (Lee et al., 2017), a relationship that was reflected in the negative correlations with TCORT and FFM and the moderate positive associations between TCORT, FCORT and %BF.

Although greater net elevations in TCORT showed a negative correlation with FFM, interestingly, net elevations in TCORT were also positively correlated with VO_2max. This may be indicative of a power/endurance tradeoff (Kraemer et al., 1995) typically noted in times of increased aerobic activity that might be seen during high volumes of soccer-specific training. Intensive aerobic training can improve endurance performance reflected by improvements in VO_2max (Bresciani et al., 2011). However, if intensive aerobic training is not met with adequate recovery, altered hormonal responses including increases in cortisol may persist. The positive correlation between VO_2max and TCORT may be a result of an overload of aerobic training, although more research is needed to determine the underlying mechanism behind this observed relationship. It is also interesting to note that a trend was seen for a moderate positive association between CK and DL. This relationship may be reflective of the work completed by the athletes, as increases in CK may be indicative of muscle breakdown and remodeling which could result in improvements in strength if coupled with a proper anabolic stimulus.

As a team, anabolic markers remained relatively stable despite the fluctuations in catabolic markers. The maintenance of anabolic markers over the course of the season, coincided with the preservation of both power and endurance as well increases in strength. Greater concentrations of anabolic hormones may help to promote protein synthesis and counteract muscle damage that may occur during soccer (Hackney and Machado, 2012). Whereas, low levels of key anabolic hormones, such as GH and IGF-1, that play a role in regulating muscle growth (Frystyk, 2010) have been associated with low muscle mass and strength (Araujo et al., 2008; Stenholm et al., 2010). This is potentially evidenced by greater individual changes in GH and IGF-1 measures positively correlated with DL in this study. In addition, a moderate positive correlation occurred between DL and TTEST.

For the female athlete, one area of particular interest is E₂ and its relationship to performance outcomes. Although no correlations between E₂ and performance or body composition were seen, overall elevations in TCORT showed a moderate negative correlation with E₂. Lowered E₂ in conjunction with increased TCORT has been associated with energy deficiencies caused by an inadequate caloric intake relative to exercise energy expenditures (Mountjoy et al., 2018). Low energy availability is shown to alter levels of various hormones leading to disruptions in optimal health and performance (Vanheest et al., 2014; Mountjoy et al., 2018). Further, E₂ was positively correlated with IGF-1, while correlations between E₂ and FTEST approached significance. IGF-1 has been shown to reflect energy status, with a combination of low energy intake and high volumes of training shown to reduce IGF-1 (Zanker and Swaine, 2000). Ultimately, the positive correlation between increased IGF-1 and E₂ may be the result of a more favorable anabolic environment (Nindl and Pierce, 2010; Koundourakis and Margioris, 2019), while declines in IGF-1 in conjunction with E₂ may be a response to an energy deficit (Nemet et al., 2004; Vanheest et al., 2014; Joro et al., 2017) associated with the high caloric expenditures characteristic of soccer athletes. More research is warranted to determine the impact of energy availability on biomarker changes as well as performance outcomes in order to identify individuals in need of a nutrition intervention. It is worth noting that although biomarker monitoring did occur every 4-weeks to account for a typical menstrual cycle, the individual menstrual status of the athletes was not “controlled.” Given the nature of this applied study, researchers sought to determine a typical hormonal environment that exists for a collegiate team during the season. Additionally, since menstrual cycle does not control playing time, we believe that this is an accurate representative profile of a female collegiate soccer team.

In addition to lack of reported information on athlete's menstrual status, the authors acknowledge several other limitations to this current study which may have impacted biomarker results. First, diet was not controlled throughout the season. Researchers sought to determine typical biomarker responses that occur in these athletes without nutritional controls. Although the findings would have benefited from the addition of dietary information, accuracy and adherence challenges associated with dietary logs limit the feasibility of this method in a team setting. In addition, sleep and mood states have also been shown to be associated with physiological responses and performance outcomes (Morgan et al., 1987; Filaire et al., 2001; Halson, 2014). Evaluating subjective measures of sleep and mood may provide useful context to interpreting physiological responses to training. However, it is important to consider the burden of implementing additional testing procedures to collegiate athletes who already experience heightened responsibilities inherent with collegiate student athlete life. Future research may consider ways to address these concerns without overburdening the athlete. Further, although consistent with an average soccer team's roster, the small sample size in the current study may limit the interpretation of these results. Future research may be warranted to determine biomarker and performance associations in female soccer athletes with similar training requirements across multiply teams.

Finally, it is important to consider that increases and decreases in a single biomarker often do not occur independently as illustrated by various convergent hormonal pathways (Lee et al., 2017). This is further evident in the correlations amongst biomarkers in this study, such as the associations of CRP and IL-6, IGF-1, and GH, along with testosterone and IGF-1. Therefore, rather than looking at one biomarker to determine the effects of training load on an athlete, it may be more advantageous to analyze multiple anabolic and catabolic biomarkers together when assessing adaptations to a training program.

The results of this study support a relationship between biomarkers, body composition, and performance outcomes in female collegiate soccer athletes. Despite the accumulated stressors, as a team, performance metrics were maintained or modestly improved in the case of SQ and BP performance. In addition, individual elevations in anabolic hormones from baseline were related to improvements in measures of strength, specifically DL performance, whereas greater net elevations in catabolic or inflammatory markers were related to declines in FFM and BP. Further, TCORT was positively correlated with VO_2max and negatively correlated with FFM indicating that an overload of aerobic soccer-specific training may have predominated in certain players during the season. It appears that biomarker monitoring may be a useful tool in team sports to determine an individual athlete's recovery needs and player readiness status. Monitoring changes in anabolic and catabolic markers throughout a season may provide an opportunity to intervene before decrements in performance occur. Changes in an athlete's physiological state can provide context for coaches and support staff to manipulate training variables in order to promote recovery or provide a greater stimulus to enhance adaptations needed to optimize performance. Additionally, biomarker monitoring may be useful to detect individual's physiological response to a given training load and serve as a guide to strength and conditioning coaches regarding periodization of strength training programs to optimize performance. Improved knowledge of the athlete's physiological state may allow for more informed decisions regarding the application of training loads by coaches and trainers to promote athlete performance and adaptation to training. This may include modifying volume and intensity of training, programmed rest and recovery, improved nutritional strategies, and substitution strategies that promote better player management (Marques-Jimenez et al., 2016). Additionally, the use of multiple biomarkers may better detect training stressors related to recovery status and performance outcomes rather than relying on performance tests alone. Through properly structured strength and conditioning programs in conjunction with appropriate assessment methods and application of training loads, it may be possible to avoid the power/endurance tradeoff and loss of FFM often seen over the course of the season in these athletes.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by Rutgers University Institutional Review Board for the Protection of Human Subjects. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

BM, AW, and SA conceived and designed the experiments. BM, AW, DS, BB, HC, MB, MA, and SA performed the experiments. BM, AW, BB, and SA analyzed the data. SA contributed materials/tools. BM, BB, MA, HC, and SA wrote the paper.

Funding

The authors declare that this study received funding from Quest Diagnostics, Grant Number 822880. The funder had the following involvement with the study: blood sample collection and analysis. The funder was not involved in study design, data analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.

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.

Acknowledgments

Special thanks to the Rutgers Women's Soccer Team.

Supplementary Material

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

References

Araujo, A. B., Travison, T. G., Bhasin, S., Esche, G. R., Williams, R. E., Clark, R. V., et al. (2008). Association between testosterone and estradiol and age-related decline in physical function in a diverse sample of men. J. Am. Geriatr. Soc. 56, 2000–2008. doi: 10.1111/j.1532-5415.2008.01965.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bangsbo, J., and Lindquist, F. (1992). Comparison of various exercise tests with endurance performance during soccer in professional players. Int. J. Sports Med. 13, 125–132. doi: 10.1055/s-2007-1021243

PubMed Abstract | CrossRef Full Text | Google Scholar

Bessa, A. L., Oliveira, V. N., Agostini, G. G., Oliveira, R. J., Oliveira, A. C., White, G. E., et al. (2016). Exercise intensity and recovery: biomarkers of injury, inflammation, and oxidative stress. J. Strength Cond. Res. 30, 311–319. doi: 10.1519/JSC.0b013e31828f1ee9

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosco, C., Tihanyi, J., and Viru, A. (1996). Relationships between field fitness test and basal serum testosterone and cortisol levels in soccer players. Clin. Physiol. 16, 317–322. doi: 10.1111/j.1475-097X.1996.tb00577.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bresciani, G., Cuevas, M. J., Molinero, O., Almar, M., Suay, F., Salvador, A., et al. (2011). Signs of overload after an intensified training. Int. J. Sports Med. 32, 338–343. doi: 10.1055/s-0031-1271764

PubMed Abstract | CrossRef Full Text | Google Scholar

Brozek, J., Grande, F., Anderson, J., and Keys, A. (1963). Densitometric analysis of body composition: revision of some quantitative assumptions. Ann. N. Y. Acad. Sci. 110, 113–140. doi: 10.1111/j.1749-6632.1963.tb17079.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Brzycki, M. (1993). Strength testing: predicting a one-rep max from reps-to-fatigue. J. Health Phys. Educ. Recreat. Dance 64, 88–90. doi: 10.1080/07303084.1993.10606684

CrossRef Full Text | Google Scholar

Cano Sokoloff, N., Misra, M., and Ackerman, K. E. (2016). Exercise, training, and the hypothalamic-pituitary-gonadal axis in men and women. Front. Horm Res. 47, 27–43. doi: 10.1159/000445154

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceesay, S., Prentice, A., Day, K., Murgatroyd, P., Goldberg, G., and Scott, W. (1989). The use of heart rate monitoring in the estimation of energy expenditure:validation study using indirect whole-body calorimetry. Br. J. Nutr. 61, 175–186. doi: 10.1079/BJN19890107

PubMed Abstract | CrossRef Full Text | Google Scholar

de Hoya, M., Cohen, D. D., sanudo, B., Carrasco, L., Alvarez-Mesa, A., del Ojo, J. J., et al. (2016). Influence of football match time-motion parameters on recovery time courses of muscle damage and jump ability. J. Sports Sci. 34, 1363–1370. doi: 10.1080/02640414.2016.1150603

PubMed Abstract | CrossRef Full Text | Google Scholar

Dempster, P., and Aitkens, S. (1995). A new air displacement method for the determination of human body composition. Med. Sci. Sports Exerc. 27, 1692–1697. doi: 10.1249/00005768-199512000-00017

PubMed Abstract | CrossRef Full Text | Google Scholar

Dueck, C. A., Manore, M. M., and Matt, K. S. (1996). Role of energy balance in athletic menstrual dysfunction. Int. J. Sport Nutr. 6, 165–190. doi: 10.1123/ijsn.6.2.165

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrucci, L., Penninx, B. W., Volpato, S., Harris, T. B., Bandeen-Roche, K., Balfour, J., et al. (2002). Change in muscle strength explains accelerated decline of physical function in older women with high interleukin-6 serum levels. J. Am. Geriatr. Soc. 50, 1947–1954. doi: 10.1046/j.1532-5415.2002.50605.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Filaire, E., Bernain, X., Sagnol, M., and Lac, G. (2001). Preliminary results on mood state, salivary testosterone:cortisol ratio and team performance in a professional soccer team. Eur. J. Appl. Physiol. 86, 179–184. doi: 10.1007/s004210100512

PubMed Abstract | CrossRef Full Text | Google Scholar

Filaire, E., Lac, G., and Pequignot, J. M. (2003). Biological, hormonal, and psychological parameters in professional soccer players throughout a competitive season. Percept. Mot. Skills 97, 1061–1072. doi: 10.2466/pms.2003.97.3f.1061

PubMed Abstract | CrossRef Full Text | Google Scholar

Frystyk, J. (2010). Exercise and the growth hormone-insulin-like growth factor axis. Med. Sci. Sports Exerc. 42, 58–66. doi: 10.1249/MSS.0b013e3181b07d2d

PubMed Abstract | CrossRef Full Text | Google Scholar

Giersch, G., Huggins, R., Benjamin, C., Adams, W., Belval, L., Curtis, R., et al. (2018). Validity and reliability of a shirt-based integrated gps sensor: 2759 board #42 June 12. Med. Sci. Sports Exerc. 50:673. doi: 10.1249/01.mss.0000538219.84744.41

CrossRef Full Text | Google Scholar

Hackney, A. C., and Machado, M. (2012). Hormonal status, creatine kinase and soccer: a need for research. Arq. Bras. Endocrinol. Metabol. 56, 683–685. doi: 10.1590/S0004-27302012000900014

PubMed Abstract | CrossRef Full Text | Google Scholar

Haff, G., and Triplett, N. T. (2016). Essentials of Strength Training and Conditioning, 4th Edn. Champaign, IL: Human Kinetics, 452–454.

Google Scholar

Halson, S. L. (2014). Monitoring training load to understand fatigue in athletes. Sports Med. 44(Suppl. 2), S139–S147. doi: 10.1007/s40279-014-0253-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Huggins, R. A., Fortunati, A. R., Curtis, R. M., Looney, D. P., West, C. A., Lee, E. C., et al. (2019). Monitoring blood biomarkers and training load throughout a collegiate soccer season. J. Strength Cond. Res. 33, 3065–3077. doi: 10.1519/JSC.0000000000002622

PubMed Abstract | CrossRef Full Text | Google Scholar

Joro, R., Uusitalo, A., DeRuisseau, K. C., and Atalay, M. (2017). Changes in cytokines, leptin, and IGF-1 levels in overtrained athletes during a prolonged recovery phase: a case-control study. J. Sports Sci. 35, 2342–2349. doi: 10.1080/02640414.2016.1266379

PubMed Abstract | CrossRef Full Text | Google Scholar

Koundourakis, N. E., and Margioris, A. N. (2019). The complex and bidirectional interaction between sex hormones and exercise performance in team sports with emphasis on soccer. Hormones 18, 151–172. doi: 10.1007/s42000-019-00115-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kraemer, W. J., French, D. N., and Paxton, N. J. (2004). Changes in exercise performance and hormonal concentrations over a big ten soccer season in starters and non starters. J. Strength Cond. Res. 18, 121–128. doi: 10.1519/00124278-200402000-00018

PubMed Abstract | CrossRef Full Text | Google Scholar

Kraemer, W. J., Patton, J. F., Gordon, S. E., Harman, E. A., Deschenes, M. R., Reynolds, K., et al. (1995). Compatibility of high-intensity strength and endurance training on hormonal and skeletal muscle adaptations. J. Appl. Physiol. 78, 976–989. doi: 10.1152/jappl.1995.78.3.976

PubMed Abstract | CrossRef Full Text | Google Scholar

Krustrup, P., Mohr, M., Ellingsgaard, H., and Bangsbo, J. (2005). Physical demands during an elite female soccer game: importance of training status. Med. Sci. Sports Exerc. 37, 1242–1248. doi: 10.1249/01.mss.0000170062.73981.94

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, E. C., Fragala, M. S., Kavouras, S. A., Queen, R. M., Pryor, J. L., and Casa, D. J. (2017). Biomarkers in sports and exercise: tracking health, performance, and recovery in athletes. J. Strength Cond. Res. 31, 2920–2937. doi: 10.1519/JSC.0000000000002122

PubMed Abstract | CrossRef Full Text | Google Scholar

Mann, J., Bryant, K., Johnstone, B., Ivey, P., and Sayers, S. (2015). The effect of physical and academic stress on illness and injury in division 1 college football players. J. Strength Cond. Res. 30, 20–25. doi: 10.1519/JSC.0000000000001055

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques-Jimenez, D., Calleja-Gonzalez, A. I., and Terrados, N. (2016). Relevant biochemical markers of recovery process in soccer. Arch Med Deportes 33, 404–412. doi: 10.2174/1875399X01710010052

CrossRef Full Text | Google Scholar

McFadden, B. A., Walker, A. J., Bozzini, B. N., Sanders, D. J., and Arent, S. M. (2020). Comparison of internal and external training loads in male and female collegiate soccer players during practices vs. games. J. Strength Cond. Res. 34, 969–974. doi: 10.1519/JSC.0000000000003485

PubMed Abstract | CrossRef Full Text | Google Scholar

Meeusen, R., Duclos, M., Foster, C., Fry, A., Gleeson, M., Nieman, D., et al. (2013). Prevention, diagnosis, and treatment of the overtraining syndrome: joint consensus statement of the european college of sport science and the american college of sports medicine. Med. Sci. Sports Exerc. 45, 186–205. doi: 10.1249/MSS.0b013e318279a10a

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, W., Brown, D., and Raglin, J. (1987). Psychological monitoring of overtraining and staleness. Br. J. Sports Med. 21, 107–114. doi: 10.1136/bjsm.21.3.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Mougios, V. (2007). Reference intervals for serum creatine kinase in athletes. Br. J. Sports Med. 41, 674–678. doi: 10.1136/bjsm.2006.034041

PubMed Abstract | CrossRef Full Text | Google Scholar

Mountjoy, M., Sundgot-Borgen, J., Burke, L., Ackerman, K. E., Blauwet, C., Constantini, N., et al. (2018). International Olympic Committee (IOC) consensus statement on relative energy deficiency in sport (RED-S): 2018 Update. Int. J. Sport Nutr. Exerc. Metab. 28, 316–331. doi: 10.1123/ijsnem.2018-0136

PubMed Abstract | CrossRef Full Text | Google Scholar

Nemet, D., Connolly, P. H., Pontello-Pescatello, A. M., Rose-Gottron, C., Larson, J. K., Galassetti, P., et al. (2004). Negative energy balance plays a major role in the IGF-I response to exercise training. J. Appl. Physiol. 96, 276–282. doi: 10.1152/japplphysiol.00654.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Nindl, B. C., and Pierce, J. R. (2010). Insulin-like growth factor 1 as a biomarker of health, fitness, and training status. Med. Sci. Sports Exerc. 42, 39–49. doi: 10.1249/MSS.0b013e3181b07c4d

CrossRef Full Text | Google Scholar

Nindl, B. C., and Pierce, J. R. (2010). Insulin-like growth factor 1 as a biomarker of health, fitness, and training status. Med. Sci. Sports Exerc. 42, 39–49. doi: 10.1249/MSS.0b013e3181b07c4d

CrossRef Full Text | Google Scholar

Nuzzo, J., Anning, J., and Scharfenberg, J. (2011). The reliability of three devices used for measuring vertical jump height. Strength Cond. 25, 2580–2590. doi: 10.1519/JSC.0b013e3181fee650

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, B. K., and Febbraio, M. (2005). Muscle-derived interleukin-6-a possible link between skeletal muscle, adipose tissue, liver, and brain. Brain Behav. Immun. 19, 371–376. doi: 10.1016/j.bbi.2005.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedlar, C. R., Newell, J., and Lewis, N. A. (2019). Blood biomarker profiling and monitoring for high-performance physiology and nutrition: current perspectives, limitations and recommendations. Sports Med. 49, 185–198. doi: 10.1007/s40279-019-01158-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, A., Santhiago, V., Papoti, M., and Gobatto, C. (2008). Psychological, biochemical and physiological responses of Brazilian soccer players during a training program. Sci. Sports 23, 66–72. doi: 10.1016/j.scispo.2007.10.010

CrossRef Full Text | Google Scholar

Silva, J. M. (1990). An analysis of the training stress syndrome in competitive athletics. J. Appl. Sport Psychol. 2, 5–20. doi: 10.1080/10413209008406417

CrossRef Full Text | Google Scholar

Sporis, G., Jukic, I., Ostojic, S. M., and Milanovic, D. (2009). Fitness profiling in soccer: physical and physiologic characteristics of elite players. J. Strength Cond. Res. 23, 1947–1953. doi: 10.1519/JSC.0b013e3181b3e141

PubMed Abstract | CrossRef Full Text | Google Scholar

Stenholm, S., Maggio, M., Lauretani, F., Bandinelli, S., Ceda, G. P., Di Iorio, A., et al. (2010). Anabolic and catabolic biomarkers as predictors of muscle strength decline: the InCHIANTI study. Rejuven. Res. 13, 3–11. doi: 10.1089/rej.2009.0891

PubMed Abstract | CrossRef Full Text | Google Scholar

Stolen, T., Chamari, K., Castagna, C., and Wisloff, U. (2005). Physiology of soccer: an update. Sports Med. 35, 501–536. doi: 10.2165/00007256-200535060-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, H., Buskirk, E., and Henschel, A. (1955). Maximal oxygen intake as an objective measure of cardio-respiratory performance. J. Appl. Physiol. 8, 73–80. doi: 10.1152/jappl.1955.8.1.73

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, W., Gordon, N., and Pescatello, L. (2000). ACSM's Guidelines for Exercise Testing and Prescription, 8th Edn. Philadelphia: Wolters Kluwer/Lippincott Williams & Wilkins Health.

PubMed Abstract | Google Scholar

Thorpe, R. T., Strudwick, A. J., Buchheit, M., Atkinson, G., Drust, B., and Gregson, W. (2015). Monitoring fatigue during the in-season competitive phase in elite soccer players. Int. J. Sports Physiol. Perform. 10, 958–964. doi: 10.1123/ijspp.2015-0004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tiidus, P. M. (2003). Influence of estrogen on skeletal muscle damage, inflammation, and repair. Exerc. Sport Sci. Rev. 31, 40–44. doi: 10.1097/00003677-200301000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Twist, C., and Highton, J. (2013). Monitoring fatigue and recovery in rugby league players. Int. J. Sports Physiol. Perform. 8, 467–474. doi: 10.1123/ijspp.8.5.467

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanheest, J. L., Rodgers, C. D., Mahoney, C. E., and De Souza, M. J. (2014). Ovarian suppression impairs sport performance in junior elite female swimmers. Med. Sci. Sports Exerc. 46, 156–166. doi: 10.1249/MSS.0b013e3182a32b72

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, A., McFadden, B., Sanders, D., Rabideau, M., Hofacker, M., and Arent, M. (2019b). Biomarker response to a competitive season in division I female soccer players. J. Strength Cond. Res. 33, 2622–2628. doi: 10.1519/JSC.0000000000003264

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, A. J., Arent, M. A., McFadden, B.A., and Arent, S. M. (2019a) Physical merformance testing. in Elite Soccer Players: Maximizing Performance Safety, eds R. M. Curtis, C. L. Benjamin, R. A. Huggins, D. J. Cases (London: Routledge)

Weber, J., and Lamb, D. (1970). Statistical and Research in Physical Education. St Louis, MO: CV Mosby Co.

Zanker, C. L., and Swaine, I. L. (2000). Responses of bone turnover markers to repeated endurance running in humans under conditions of energy balance or energy restriction. Eur. J. Appl. Physiol. 83, 434–440. doi: 10.1007/s004210000293

PubMed Abstract | CrossRef Full Text | Google Scholar