Cross-sectional associations between cardiorespiratory fitness and NMR-derived metabolic biomarkers in children – the PANIC study

Objective: Cardiorespiratory fitness has been inversely associated with cardiovascular risk across the lifespan. Some studies in adults suggest that higher cardiorespiratory fitness is associated with cardioprotective metabolite profile, but the evidence in children is lacking. Therefore, we investigated the cross-sectional association of cardiorespiratory fitness with serum nuclear magnetic resonance derived metabolic biomarkers in children.

Methods: A population sample of 450 children aged 6–8 years was examined. Cardiorespiratory fitness was assessed by a maximal exercise test on a cycle ergometer and quantified as maximal power output normalised for lean body mass assessed by dual-energy X-ray absorbtiometry. Serum metabolites were assessed using a high throughput nuclear magnetic resonance platform. The data were analysed using linear regression analyses adjusted for age and sex and subsequently for body fat percentage (BF%) assessed by DXA.

Results: Cardiorespiratory fitness was directly associated with high density lipoprotein (HDL) cholesterol (β=0.138, 95% CI=0.042 to 0.135, p=0.005), average HDL particle diameter (β=0.102, 95% CI=0.004 to 0.199, p=0.041), and the concentrations of extra-large HDL particles (β=0.103, 95% CI=0.006 to 0.201, p=0.038), large HDL particles (β=0.122, 95% CI=0.025 to 0.220, p=0.014), and medium HDL particles (β=0.143, 95% CI=0.047 to 0.239, p=0.004) after adjustment for age and sex. Higher cardiorespiratory fitness was also associated with higher concentrations of ApoA1 (β=0.145, 95% CI=0.047 to 0.242, p=0.003), glutamine (β=0.161, 95% CI=0.064 to 0.257, p=0.001), and phenylalanine (β=0.187, 95% CI=0.091 to 0.283, p<0.001). However, only the direct associations of cardiorespiratory fitness with the concentrations of HDL cholesterol (β=0.114, 95% CI=0.018 to 0.210, p=0.021), medium HDL particles (β=0.126, 95% CI=0.030 to 0.223, p=0.010), ApoA1 (β=0.126, 95% CI=0.030 to 0.223, p=0.011), glutamine (β=0.147, 95% CI=0.050 to 0.224, p=0.003), and phenylalanine (β=0.217, 95% CI=0.122 to 0.311, p<0.001) remained statistically significant after further adjustment for BF%.

Conclusions: Higher cardiorespiratory fitness was associated with a cardioprotective biomarker profile in children. Most associations were independent of BF% suggesting that the differences in serum metabolites between children are driven by cardiorespiratory fitness and not adiposity.

Introduction

Low cardiorespiratory fitness has been associated with increased risk for cardiovascular diseases and cardiovascular events in adults (1, 2). Nevertheless, pathophysiological processes for atherosclerotic cardiovascular diseases start already in childhood. Accordingly, the early signs of atherosclerotic cardiovascular diseases such as increased arterial stiffness (3, 4) and increased carotid intima media thickness and distensibility (4, 5), have been detected already in children and adolescents. In addition, some evidence suggests that increasing cardiorespiratory fitness since youth could stagnate the development of early signs of atherosclerotic cardiovascular diseases over the life course (6). Nevertheless, cardiorespiratory fitness has been found to have weak if any associations with traditional cardiovascular disease risk factors, such as insulin resistance, dyslipidaemia, and blood pressure in children (7–10), leaving a knowledge gap on how cardiorespiratory fitness could contribute to cardiovascular health in a general population of children. Thus, while cardiorespiratory fitness may protect against atherosclerotic cardiovascular diseases, the mechanisms are not well understood, especially in children and more studies are warranted. Serum metabolomics provides a novel approach to deepen our understanding of the mechanisms underlying the health benefits of cardiorespiratory fitness.

Impaired lipid metabolism, such as increased low-density lipoprotein (LDL) cholesterol concentration and apolipoprotein B (ApoB), is well-known risk factor for atherosclerotic cardiovascular diseases (11). Metabolomics can provide further characterisation of lipids improving the cardiovascular risk estimation (12). For example, increased very-low-density lipoprotein (VLDL) concentration and reduced high density lipoprotein (HDL) concentration and size have been associated with an increased risk of atherosclerotic cardiovascular diseases (13–16). Furthermore, traditional lipid biomarkers, such as LDL cholesterol, HDL cholesterol, and triglycerides were not associated with carotid atherosclerosis or carotid-femoral pulse wave velocity, whereas a total HDL particle concentration and average HDL size were inversely associated with them in obese youth aged 18 years (17). The studies utilising metabolomics have also found that increased concentrations of branched-chain amino acids are associated with increased risk of atherosclerotic cardiovascular diseases (15, 18) and increased carotid artery intima media thickness, and exercise induced myocardial ischaemia (19) in adults. However, it should be noted that there are only few studies on the associations between metabolites and the measures of cardiovascular health among youth. However, one study found that none of the investigated metabolites were associated with carotid intima-media thickness or pulse wave velocity after adjustment for body mass index and blood pressure in children (20).

Cardiorespiratory fitness has been inversely associated with lipid metabolites and branched-chain amino acids in adults (21), but there are paucity of data in children. Nevertheless, one small study in adolescents found that cardiorespiratory fitness was inversely associated with amino acids valerate, glutamate, and tyrosine (22). However, cardiorespiratory fitness was assessed by maximal oxygen uptake (V̇O_2max) normalised for whole body mass, a measure of cardiorespiratory fitness strongly confounded by adiposity (23). Similarly, another study in young adults showed that higher cardiorespiratory fitness quantified, as estimated V̇O_2max normalised for whole body mass, was inversely associated with the concentrations of extra-large to small size VLDL, large to small size LDL, ApoB, and ApoB to apolipoprotein A1 (ApoA1) ratio and directly associated with the concentrations of extra-large to medium size HDL and ApoA1 (24). Furthermore, cardiorespiratory fitness has been inversely associated with branched-chain amino acids and phenylalanine, tyrosine and positively associated with glutamine (24). However, when body fat percentage (BF%) was controlled for, cardiorespiratory fitness was associated only with the concentrations of large HDL, ApoA1, and medium size VLDL (24) suggesting that adiposity is a strong confounder in these associations.

There is limited evidence on the associations of cardiorespiratory fitness with serum nuclear magnetic resonance (NMR) derived metabolic biomarkers in children, although pathophysiological process leading to atherosclerotic cardiovascular diseases often begin in childhood. Therefore, we first investigated the associations of cardiorespiratory fitness with serum NMR-derived biomarkers related to lipoproteins, triglycerides, apolipoproteins, and amino acids in a general population of children. Second, we investigated whether BF% modifies these associations. Finally, because HDL characteristics potentially have specific effects on cardiovascular health and the determinants of cardiorespiratory fitness (25–27), we also investigated the associations of cardiorespiratory fitness with various HDL characteristics.

Methods

Participants

The present cross-sectional data are from the baseline assessments of the Physical Activity and Nutrition in Children (PANIC) Study, which is an 8-year physical activity and dietary intervention study and a long-term follow-up study in a population sample of children from the city of Kuopio, Finland (28). The Research Ethics Committee of the Hospital District of Northern Savo approved the study protocol in 2006 (Statement 69/2006). The parents or caregivers of the children gave their written informed consent, and the children provided their assent to participation. The PANIC study has been carried out in accordance with the principles of the Declaration of Helsinki as revised in 2008.

Altogether 736 children 6–8 years of age from primary schools of Kuopio were invited to participate in the baseline examination in 2007–2009. A total of 512 children, who represented 70% of those invited, participated in the baseline examinations. Six children were excluded from the study at baseline because of physical disabilities that could hamper participation in the intervention or no time or motivation to attend in the study. The participants did not differ in sex distribution, age, or body mass index standard deviation score (BMI-SDS) from all children who started the first grade in 2007–2009 based on data from the standard school health examinations performed for all Finnish children before the first grade (data not shown). Complete data on variables used in the analyses on the associations of cardiorespiratory fitness with serum metabolites were available for 450 children (217 girls, 233 boys).

Assessment of body size and composition

Whole body mass was measured twice with the children having fasted for 12 hours, emptied the bladder, and standing in light underwear by a calibrated InBody^® 720 bioelectrical impedance device (Biospace, Seoul, South Korea) to an accuracy of 0.1 kg. The mean of these two values was used in the analyses. Stature was measured three times with the children standing in the Frankfurt plane without shoes using a wall-mounted stadiometer to an accuracy of 0.1 cm. The mean of the nearest two values was used in the analyses. BMI was calculated by dividing body mass (kg) by body height (m) squared. BMI-SDS was calculated based on Finnish reference data (29). The prevalence of overweight and obesity was defined using the cut-off values provided by Cole et al. (30). Total fat mass, BF%, and lean mass (LM) were measured by the Lunar^® dual-energy X-ray absorptiometry device (GE Medical Systems, Madison, WI, USA) using standardised protocols (31).

Assessment of cardiorespiratory fitness

We assessed cardiorespiratory fitness by a maximal exercise test using an electromagnetically braked Ergoselect 200 K^® cycle ergometer coupled with a paediatric saddle module (Ergoline, Bitz, Germany) (32). The exercise test protocol included a 2.5-minute anticipatory period with the child sitting on the ergometer; a 3-minute warm-up period with a workload of 5 watts; a 1-minute steady-state period with a workload of 20 watts; an exercise period with an increase in the workload of 1 watt per 6 seconds until exhaustion, and a 4-minute recovery period with a workload of 5 watts.

The children were asked to keep the cadence stable and within 70–80 revolutions per minute. Exhaustion was defined as the inability to maintain the cadence above 65 revolutions per minute regardless of vigorous verbal exhortation. The exercise test was considered maximal by an experienced physician (TT) supervising the test, if objective and subjective criteria (heart rate >85% of predicted, sweating, flushing, inability to continue exercise test regardless of strong verbal encouragement) indicated maximal effort and maximal cardiovascular capacity. Heart rate was measured continuously during the last five minutes of the supine rest prior to commencing the exercise test protocol right through to the 5-minute supine post-exercise rest period using a 12-lead electrocardiogram registered by the Cardiosoft^® V6.5 Diagnostic System (GE Healthcare Medical Systems, Freiburg, Germany) and the highest heart rate during the test was defined as peak heart rate. Maximal power output (W_max) measured at the end of the exercise test divided by LM^-1 were used as a measure of cardiorespiratory fitness. We used W_max as a measure of cardiorespiratory fitness because we did not perform respiratory gas analyses at baseline and it has been found to be a good surrogate measure of cardiorespiratory fitness in children (33). W_max x kg of LM^-1 was no associated with LM (β=0.072, 95% CI=-0.013 to 0.156) and the sex interaction was not statistically significant (p=0.235) suggesting the validity of the scaling procedure.

Assessment of metabolic biomarkers

A research nurse took blood samples in the morning, after children had fasted overnight for at least 12 hours. Blood was immediately centrifuged and stored at a temperature of -75°C until biochemical analyses. The Nightingale high-throughput NMR platform was used to quantify serum metabolic biomarkers (34). The Nightingale NMR platform quantifies different metabolic biomarkers in absolute concentration units. Based on the aims of the study and previous literature, we selected metabolites related to lipoprotein, triglyceride, apolipoprotein, and amino acid metabolism and specifically focused on 49 HDL characteristics.

Assessment of confounding factors

The research physician assessed pubertal status using the 5-stage scale described by Tanner (35, 36). The boys were defined as having entered clinical puberty if their testicular volume assessed by an orchidometer was ≥4 mL (stage ≥2). The girls were defined having entered clinical puberty if their breast development had started (stage ≥2). Maturity was assessed as the difference between the current age from the age at predicted peak height velocity and it was computed using a sex-specific formula described by Moore et al. (37). Physical activity energy expenditure (PAEE) and time accumulated in moderate to vigorous physical activity (MVPA) were assessed by a combined heart rate and movement monitor (38). We included all children who had device-assessed data on physical activity regardless of wearing time requirements used in our previous studies (38, 39). The mean (standard deviation, SD) wear time was 113 (40) hours for the 402 children who met the wear time requirements and 61 (33) hours for the 33 children who did not meet these requirements. PAEE, MVPA, BF%, cardiorespiratory fitness, maturity, or diet quality did not differ between these two groups (p>0.290). Diet quality was assessed by four day dietary records (28), and the overall diet quality was computed using the Finnish Children Healthy Eating Index (FCHEI) (41). Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) as a measure of insulin resistance was computed from the fasting plasma glucose and serum insulin as prescribed earlier (42). Missing data were replaced using the sample mean values.

Statistical methods

The analyses were performed using the Jamovi statistical software, version 2.2 (Jamovi project 2021). First, we investigated the differences in basic characteristics between girls and boys using the Welch’s -test for normally distributed continuous variables, the Mann-Whitney U -test for continuous variables with skewed distributions, and the χ² test for categorical variables. Second, we studied the associations of cardiorespiratory fitness with metabolites using the linear regression analyses adjusted for age and sex (Model 1) and additionally for BF% (Model 2). We also investigated whether sex or BF% modified the associations of cardiorespiratory fitness with serum metabolites by adding the interaction term of sex or BF% and cardiorespiratory fitness into the models. If the interaction was statistically significant, we analysed the associations of cardiorespiratory fitness with metabolites separately among children below and at or above the sex-specific median of BF% or among girls and boys. The data were corrected for multiple comparisons using the Benjamini-Hochberg false discovery rate (FDR) using the FDR value of 0.2 (FDR_0.2) and 0.1 (FDR_0.1). The data were further adjusted for maturity, PAEE, MVPA, FCHEI, or HOMA-IR. The data on the associations between cardiorespiratory fitness and amino acids were further adjusted for protein intake as a percentage of total energy. We did not adjust the data for clinical puberty because few children had signs of clinical puberty.

Results

Basic characteristics

Girls were shorter and more likely to have entered clinical puberty and had advanced maturity, higher BF%, and lower LM than boys (Table 1). Girls also had lower cardiorespiratory fitness and accumulated less PAEE and MVPA than boys.

TABLE 1

Table 1 Characteristics of participants.

Associations of cardiorespiratory fitness with metabolic biomarkers

Cardiorespiratory fitness was directly associated with HDL cholesterol concentration, average HDL diameter, and the concentrations of extra-large, large, and medium size HDL particles as well as ApoA1, glutamine, and phenylalanine after adjustment for age and sex (Table 2). These associations remained after FDR_0.2 correction, and even after FDR_0.1 correction, except the association between cardiorespiratory fitness and the concentration of extra-large HDL particles. The direct associations of cardiorespiratory fitness with HDL cholesterol concentration, the concentration of medium size HDL particles, ApoA1, glutamine, and phenylalanine remained statistically significant even after further adjustment for BF%. Further adjustment for maturity, PAEE, MVPA, FCHEI, HOMA-IR, or protein intake had no effect on the magnitude of these associations (data not shown).

TABLE 2

Table 2 Associations of cardiorespiratory fitness with NMR-derived metabolic biomarkers.

BF% modified the associations of cardiorespiratory fitness with LDL cholesterol concentration, average LDL diameter, and the concentrations of small HDL particles and ApoB (Table 2). Cardiorespiratory fitness was directly associated with average LDL diameter (β=0.198, 95% CI=0.062 to 0.334, p=0.004) and the concentration of small HDL particles (β=0.139, 95% CI=0.004 to 0.274, p=0.043) in children with higher BF%, but not in children with lower BF% (β=-0.074, 95% CI=-0.215 to 0.067, p=0.301, and β=-0.010, 95% CI=-0.130 to 0.151, p=0.887, respectively). Further adjustment for maturity, FCHEI, PAEE, MVPA, or HOMA-IR had no effect on the magnitude of these associations (data not shown). Other associations among children with higher or lower BF% were not statistically significant.

Associations of cardiorespiratory fitness with HDL characteristics

The associations of cardiorespiratory fitness with HDL characteristics are presented in Supplementary Table. Cardiorespiratory fitness was directly associated with phospholipids, cholesteryl esters, free cholesterol, and total lipids in HDL. Cardiorespiratory fitness was also directly associated with total concentration, concentration of total lipids, phospholipids, cholesterol, cholesteryl esters, and free cholesterol of medium HDL particles and with the concentration of total lipids and phospholipids in small HDL particles after adjustment for age, sex, and BF%. All but one of these associations remained after FDR_0.2 and FDR_0.1 correction. The only exception was the association between cardiorespiratory fitness and phospholipids in extra-large HDL particles that did not remain after FDR_0.1 correction. Further adjustments for maturity, PAEE, MVPA, and HOMA-IR had no effect on the magnitude of the associations. However, adjustment for FCHEI slightly attenuated the association of cardiorespiratory fitness with the concentration of total lipids (β=0.081, 95% CI=-0.016 to 0.179, p=0.103), phospholipids (β=0.092, 95% CI=-0.006 to 0.190, p=0.065), and free cholesterol (β=0.093, 95% CI=-0.005 to 0.192, p=0.063) in small HDL particles.

Sex modified the associations of cardiorespiratory fitness and some HDL particle characteristics (p ≤ 0.05 for cardiorespiratory fitness x sex interactions). Cardiorespiratory fitness was inversely associated with the proportion of phospholipids in large HDL particles (β=-0.149, 95% CI=-0.279 to -0.020, p=0.024) and proportion of cholesterol in large HDL particles (β=0.133, 95% CI=0.003 to 0.263, p=0.045) among boys, but not in girls (β=0.118, 95% CI=-0.021 to 0.257, p=0.095 and β=-0.125, 95% CI=-0.263 to 0.014, p=0.077, respectively). However, the inverse associations of cardiorespiratory fitness with HDL characteristics in boys were not statistically significant after further adjustment for BF% (β=0.100, 95% CI=-0.023 to 0.222, p=0.110 and β=0.095, 95% CI=-0.035 to 0.226, p=0.152, respectively).

BF% modified the associations between cardiorespiratory fitness with characteristics of small HDL particles. Cardiorespiratory fitness was directly associated with free cholesterol in small HDL particles in children with higher BF% (β=0.146, 95% CI=0.011 to 0.282, p=0.035), but not in children with lower BF% (β=0.034, 95% CI=0.107 to 0.175, p=0.633). Further adjustments had no effect on the magnitude of this association (data not shown).

Discussion

We found that cardiorespiratory fitness was directly associated with serum HDL characteristics, ApoA1, glutamine, and phenylalanine in a general population of children. While adiposity explained part of the associations between cardiorespiratory fitness and HDL characteristics, most associations of cardiorespiratory fitness and metabolites were independent of adiposity, maturity, physical activity, diet quality, and insulin resistance. Cardiorespiratory fitness had weak if any associations with other lipoproteins, amino acids, or triglycerides. Furthermore, we found that adiposity modified the associations of cardiorespiratory fitness with average LDL diameter and the concentration of small HDL particles.

Higher cardiorespiratory fitness has been associated with higher HDL cholesterol, lower total cholesterol, lower total cholesterol to HDL cholesterol ratio, and lower triglycerides in children and adolescents (7), but the evidence remains mixed and thus prevents firm conclusions (8). Furthermore, higher V̇O_2max normalised for body mass has been associated with higher total cholesterol to HDL cholesterol ratio, while V̇O_2max normalised for lean body mass exhibited no association with blood lipids (9). This suggests that adiposity confounds the associations between cardiorespiratory fitness and blood lipids. However, traditional lipid measurements may not be sensitive enough to provide a full picture of the associations of cardiorespiratory fitness with atherogenic and cardioprotective blood lipid profile and previous studies have not used NMR metabolomics to assess blood HDL characteristics in children. We found consistent associations of cardiorespiratory fitness with HDL characteristics, but controlling for adiposity attenuated several of these associations. Nevertheless, several associations, particularly those with medium size HDL particle characteristics, were also independent of adiposity and other potential confounders. One reason for the positive associations between cardiorespiratory fitness and HDL characteristics could be the positive associations of HDL particles on cardiac and artery structures and functions. For example, HDL particles participate in endothelial vasodilatory functions by stimulating the release and production of nitric oxide from the artery wall (43), and arterial vasolidatory function is important for high cardiorespiratory fitness (44). Accordingly, a recent study in youth suggested that healthy arterial structures could enhance cardiorespiratory fitness (4). HDL may also directly or indirectly improve cardiac structures and functions and thereby increase cardiorespiratory fitness (27). Moreover, HDL has been found to improve cellular and mitochondrial functions (25, 26). Therefore, it is possible that HDL improves cardiovascular and skeletal muscle structures and functions that are central components of cardiorespiratory fitness. Finally, although the cross-sectional associations of cardiorespiratory fitness with HDL characteristics observed in our study were largely independent of current physical activity levels, these associations may reflect the effects of previous longer-term physical activity behaviour (45).

Children with higher adiposity exhibited a direct association between cardiorespiratory fitness and the concentration of small HDL particles. Duparc et al. (12) reported an inverse association between the concentration of small HDL particles and cardiovascular mortality. Small HDL particles have been suggested to provide cardioprotection through cholesterol efflux, antioxidant, anti-inflammatory, cytoprotective, and anti-thrombotic mechanisms (43, 46, 47). They also protect LDL from oxidation (47). Therefore, it is possible that higher concentration of small HDL particles along with higher cardiorespiratory fitness in children is an adaptive mechanism protecting against cardiovascular damage caused by increased adiposity. Consistently, we have previously found that physical activity is inversely associated with biomarkers of inflammation particularly in children with higher levels of adiposity (48). However, Shah et al. (17) found a direct association between the concentration of small HDL particles and pulse wave velocity in obese youth indicating a worse arterial health. It is possible that these latter observations reflect a decreased antioxidative capacity of small HDL particles observed in adults those with type 2 diabetes (49), suggesting that not only HDL particle size but also composition is important for arterial health. Nevertheless, it is yet to be investigated whether higher small HDL concentration improves cardiovascular health among children and adolescents with higher levels of adiposity and cardiorespiratory fitness compared to their peers with lower cardiorespiratory fitness.

In line with the studies in adults (24), we observed a direct association between cardiorespiratory fitness and ApoA1 concentration. ApoA1 is a main protein component of HDL and ApoA1 containing particles mediate the reverse cholesterol transport (50). Furthermore, a higher ApoA1 has been associated with better artery structure and function (51) suggesting that the mechanisms explaining these associations are probably similar than those of HDL characteristics (50). Furthermore, we did not find associations of cardiorespiratory fitness with LDL or VLDL characteristics. Previous studies in adults have observed such associations, but they have used measures of cardiorespiratory fitness that are confounded by adiposity, perhaps explaining the discrepancy between the results. Furthermore, we found that cardiorespiratory fitness was directly associated with average LDL diameter. Small dense LDL is a potent atherogenic lipoprotein (40, 52), and LDL particle size has been inversely associated with insulin resistance (17, 53). Therefore, it is possible that increased average LDL diameter would be seen as a cardioprotective response in children with higher BF% and increased insulin resistance (42). However, it remains unknown whether increased average LDL diameter would serve as a cardioprotective mechanisms between cardiorespiratory fitness and cardiovascular health.

We found that cardiorespiratory fitness was directly associated with amino acid glutamine in conjunction with some studies in adults (24). We also found, in contrast to the results by Kujala et al. (24), that the association was independent of BF%. One explanation why cardiorespiratory fitness is associated with glutamine may be that glutamine serves a precursor for l-arginine, that is used to synthetise a vasodilator nitric oxide (54, 55), which in turn is expected to improve cardiorespiratory fitness. It is also possible that an increased glutamine concentration reflects reduced systemic inflammation and oxidative stress (54, 55). The use of glutamine supplements has been associated with decreased markers of fatigue (56), that in turn may improve cardiorespiratory fitness. Taken together, these results suggest that cardiorespiratory fitness could be associated with cardioprotective and performance enhancing metabolic alterations related to amino acids metabolism. However, in contrast to adult studies, we found a positive association between cardiorespiratory fitness and phenylalanine (24) that has been directly associated with the risk of cardiovascular events (57). The reason for the direct association between cardiorespiratory fitness and phenylalanine in our study is unclear. However, the results of some studies suggest that l-phenylalanine could be associated with improved endothelial functions (58, 59), that could partly explain our finding. Nevertheless, more research is warranted to investigate the association between cardiorespiratory fitness and phenylalanine and its role in cardiovascular health since childhood.

While branched-chain amino acids have been associated with increased risk of cardiovascular disease in adults (15, 18), we did not find an association between cardiorespiratory fitness and these amino acids in children. This could be because branched-chain amino acids are more related to adiposity, as suggested by previous studies on the associations of adiposity and cardiorespiratory fitness with amino acids in adults (24). Therefore, our results together with previous studies suggest that the cardioprotective role of higher cardiorespiratory fitness is not explained by its effects on branched-chain amino acids.

The strengths of the present study include the valid and reproducible measurements of cardiorespiratory fitness using an exercise test until exhaustion, serum metabolic biomarkers using a NMR platform, and body composition using whole-body DXA in a population sample of children. We were also able to control for several confounding factors, including physical activity, diet quality, insulin resistance, and maturation. Nevertheless, we used W_max as a measure of cardiorespiratory fitness instead of directly measured V̇O_2max. Additional studies investigating whether the associations of maximal oxygen uptake with serum metabolites are similar to those of peak performance with these metabolites are warranted. We also investigated only the associations between cardiorespiratory fitness and NMR-derived metabolic biomarkers and studies investigating whether HDL characteristics mediate the associations between cardiorespiratory fitness and arterial structures and functions. Finally, our study was cross-sectional which limits our ability to make causal inferences.

In conclusion, we found that higher cardiorespiratory fitness was associated with several cardioprotective NMR-derived metabolic biomarkers, especially those of related to HDL characteristics, independent of several confounding factors. Future studies are warranted to investigate whether the identified metabolites and especially those of HDL characteristics mediate the cardioprotective effects of cardiorespiratory fitness on arterial structures and functions since childhood.

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 studies involving human participants were reviewed and approved by The Research Ethics Committee of the Hospital District of Northern Savo. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author contributions

EH, MHL, ML, US, TL participated the conception of the study. TT, AV, NL, US, TL collected the data. EH conducted the analyses and produced the first draft of the manuscript. All authors participated in drafting and revising the manuscript and provided significant intellectual contribution to the manuscript, and all authors approved the final version of the manuscript. All authors agree to be accountable for the work and to ensure that any questions relating to the accuracy and integrity of the paper are investigated and properly resolved. All authors contributed to the article and approved the submitted version.

Funding

The PANIC Study has financially been supported by Ministry of Education and Culture of Finland, Ministry of Social Affairs and Health of Finland, Research Committee of the Kuopio University Hospital Catchment Area (State Research Funding), Finnish Innovation Fund Sitra, Social Insurance Institution of Finland, Finnish Cultural Foundation, Foundation for Paediatric Research, Diabetes Research Foundation in Finland, Finnish Foundation for Cardiovascular Research, Juho Vainio Foundation, Paavo Nurmi Foundation, Yrjö Jahnsson Foundation, and the city of Kuopio. Moreover, the PhD students and postdoctoral researchers of The PANIC Study have been supported by Program for Clinical Research and Program for Health Sciences of Doctoral School of University of Eastern Finland, Finnish Doctoral Programs in Public Health, Päivikki and Sakari Sohlberg Foundation, Paulo Foundation, Jalmari and Rauha Ahokas Foundation, Aarne and Aili Turunen Foundation, Finnish Medical Foundation, Jenny and Antti Wihuri Foundation, Kuopio Naturalists’ Society, Olvi Foundation, and the city of Kuopio. EH has been funded by the Juho Vaino foundation and the Finnish Foundation for Cardiovascular Research. The sponsors had no role in designing the study, the collection, analysis, or interpretation of the data, the writing of the report, or the decision to submit the manuscript for publication.

Acknowledgments

We are grateful to the members of the PANIC research team for their contribution in acquisition of data. We are also indebted to all children and their parents participating in the PANIC 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.

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/fendo.2022.954418/full#supplementary-material

References

1. Kodama S, Saito K, Tanaka S, Maki M, Yachi Y, Asumi M, et al. Cardiorespiratory fitness as women: A meta a quantitative predictor of all-cause mortality and cardiovascular events in healthy men and -analysis. JAMA (2009) 301:2024. doi: 10.1001/jama.2009.681

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Myers J, McAuley P, Lavie CJ, Despres JP, Arena R, Kokkinos P. Physical activity and cardiorespiratory fitness as major markers of cardiovascular risk: Their independent and interwoven importance to health status. Progr Cardiovasc Dis (2015) 57:306–14. doi: 10.1016/j.pcad.2014.09.011

CrossRef Full Text | Google Scholar

3. Veijalainen A, Tompuri T, Haapala EA, Viitasalo A, Lintu N, Väistö J, et al. Associations of cardiorespiratory fitness, physical activity, and adiposity with arterial stiffness in children. Scandinavian J Med Sci Sports (2016) 26:943–50. doi: 10.1111/sms.12523

CrossRef Full Text | Google Scholar

4. Agbaje AO, Barker AR, Tuomainen TP. Cardiorespiratory fitness, fat mass, and cardiometabolic health with endothelial function, arterial elasticity, and stiffness. Med Sci Sports Exerc (2022) 54:141–52. doi: 10.1249/MSS.0000000000002757

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Pahkala K, Laitinen TT, Heinonen OJ, Viikari JSA, Rönnemaa T, Niinikoski H, et al. Association of fitness with vascular intima-media thickness and elasticity in adolescence. Pediatrics (2013) 132:e77–84. doi: 10.1542/peds.2013-0041

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ferreira I, Twisk JWR, Stehouwer CDA, Van Mechelen W, Kemper HCG. Longitudinal changes in V̇O2max: Associations with carotid IMT and arterial stiffness. Med Sci Sports Exerc (2003) 35:1670–8. doi: 10.1249/01.MSS.0000089247.37563.4B

PubMed Abstract | CrossRef Full Text | Google Scholar

7. García-Hermoso A, Ramírez-Vélez R, García-Alonso Y, Alonso-Martínez AM, Izquierdo M. Association of cardiorespiratory fitness levels during youth with health risk later in life: A systematic review and meta-analysis. JAMA Ped (2020) 174:952–60. doi: 10.1001/jamapediatrics.2020.2400

CrossRef Full Text | Google Scholar

8. Mintjens S, Menting MD, Daams JG, van Poppel MNM, Roseboom TJ, Gemke RJBJ. Cardiorespiratory fitness in childhood and adolescence affects future cardiovascular risk factors: A systematic review of longitudinal studies. Sports Med (2018) 48:2577–605. doi: 10.1007/s40279-018-0974-5

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lätt E, Mäestu J, Rääsk T, Jürimäe T, Jürimäe J. Cardiovascular fitness, physical activity, and metabolic syndrome risk factors among adolescent estonian boys: A longitudinal study. Am J Hum Biol (2016) 28:782–8. doi: 10.1002/ajhb.22866

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lang JJ, Tremblay MS, Ortega FB, Ruiz JR, Tomkinson GR. Review of criterion-referenced standards for cardiorespiratory fitness: what percentage of 1 142 026 international children and youth are apparently healthy? Br J Sports Med (2019) 53:953–8. doi: 10.1136/bjsports-2016-096955

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Baigent C, Koskinas KC, Casula M, Badimon L, Chapman MJ, Backer GGD, et al. 2019 ESC/EAS guidelines for the management of dyslipidaemias: lipid modification to reduce cardiovascular risk. Eur Heart J (2020) 41:111–88. doi: 10.1093/eurheartj/ehz455

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Duparc T, Ruidavets JB, Genoux A, Ingueneau C, Najib S, Ferrières J, et al. Serum level of HDL particles are independently associated with long-term prognosis in patients with coronary artery disease: The GENES study. Sci Rep (2020) 10:8138. doi: 10.1038/s41598-020-65100-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Iliou A, Mikros E, Karaman I, Elliott F, Griffin JL, Tzoulaki I, et al. Metabolic phenotyping and cardiovascular disease: an overview of evidence from epidemiological settings. Heart (2021) 107:1123–9. doi: 10.1136/heartjnl-2019-315615

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Kontush A. HDL particle number and size as predictors of cardiovascular disease. Front Pharmacol (2015) 6:218. doi: 10.3389/fphar.2015.00218

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ruiz-Canela M, Toledo E, Clish CB, Hruby A, Liang L, Salas-Salvadó J, et al. Plasma branched-chain amino acids and incident cardiovascular disease in the PREDIMED trial. Clin Chem (2016) 62:582–92. doi: 10.1373/clinchem.2015.251710

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Newgard CB. Metabolomics and metabolic diseases: Where do we stand? Cell Metab (2017) 25:43–56. doi: 10.1016/j.cmet.2016.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Shah AS, Davidson WS, Gao Z, Dolan LM, Kimball TR, Urbina EM. Superiority of lipoprotein particle number to detect associations with arterial thickness and stiffness in obese youth with and without prediabetes. J Clin Lipidol (2016) 10:610–8. doi: 10.1016/j.jacl.2016.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

18. McGarrah RW, Crown SB, Zhang GF, Shah SH, Newgard CB. Cardiovascular metabolomics. Circ Res (2018) 122:1238–58. doi: 10.1161/CIRCRESAHA.117.311002

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Magnusson M, Lewis GD, Ericson U, Orho-Melander M, Hedblad B, Engström G, et al. A diabetes-predictive amino acid score and future cardiovascular disease. Eur Heart J (2013) 34:1982–9. doi: 10.1093/eurheartj/ehs424

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Juonala M, Ellul S, Lawlor DA, Santos Ferreira DL, Carlin JB, Cheung M, et al. A cross-cohort study examining the associations of metabolomic profile and subclinical atherosclerosis in children and their parents: The child health CheckPoint study and Avon longitudinal study of parents and children. JAHA (2019) 8:e011852. doi: 10.1161/JAHA.118.011852

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Carrard J, Guerini C, Appenzeller-Herzog C, Infanger D, Königstein K, Streese L, et al. The metabolic signature of cardiorespiratory fitness: A systematic review. Sports Med (2022) 52:527–46. doi: 10.1007/s40279-021-01590-y

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Duft RG, Castro A, Bonfante ILP, Lopes WA, da Silva LR, Chacon-Mikahil MPT, et al. Serum metabolites associated with increased insulin resistance and low cardiorespiratory fitness in overweight adolescents. Nutr Metab Cardiovasc Dis (2022) 32:269–78. doi: 10.1016/j.numecd.2021.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Armstrong N, Welsman J. Clarity and confusion in the development of youth aerobic fitness. Front Physiol (2019) 10:979. doi: 10.3389/fphys.2019.00979

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kujala UM, Vaara JP, Kainulainen H, Vasankari T, Vaara E, Kyröläinen H. Associations of aerobic fitness and maximal muscular strength with metabolites in young men. JAMA Netw Open (2019) 2:e198265. doi: 10.1001/jamanetworkopen.2019.8265

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Lehti M, Donelan E, Abplanalp W, Al-Massadi O, Habegger KM, Weber J, et al. High-density lipoprotein maintains skeletal muscle function by modulating cellular respiration in mice. Circulation (2013) 128:2364–71. doi: 10.1161/CIRCULATIONAHA.113.001551

PubMed Abstract | CrossRef Full Text | Google Scholar

26. White CR, Datta G, Giordano S. High-density lipoprotein regulation of mitochondrial function. Adv Exp Med Biol (2017) 982:407–29. doi: 10.1007/978-3-319-55330-6_22

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mangner N, Scheuermann K, Winzer E, Wagner I, Hoellriegel R, Sandri M, et al. Childhood obesity: Impact on cardiac geometry and function. JACC: Cardiovasc Imaging (2014) 7:1198–205.

PubMed Abstract | Google Scholar

28. Eloranta AM, Lindi V, Schwab U, Kiiskinen S, Kalinkin M, Lakka HM, et al. Dietary factors and their associations with socioeconomic background in Finnish girls and boys 6–8 years of age: the PANIC study. Eur J Clin Nutr (2011) 65:1211–8. doi: 10.1038/ejcn.2011.113

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Saari A, Sankilampi U, Hannila ML, Kiviniemi V, Kesseli K, Dunkel L. New Finnish growth references for children and adolescents aged 0 to 20 years: Length/height-for-age, weight-for-length/height, and body mass index-for-age. Ann Med (2011) 43:235–48. doi: 10.3109/07853890.2010.515603

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Cole TJ, Bellizzi MC, Flegal KM, Dietz WH. Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ (2000) 320:1240. doi: 10.1136/bmj.320.7244.1240

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Tompuri TT, Lakka TA, Hakulinen M, Lindi V, Laaksonen DE, Kilpeläinen TO, et al. Assessment of body composition by dual-energy X-ray absorptiometry, bioimpedance analysis and anthropometrics in children: the physical activity and nutrition in children study. Clin Physiol Funct Imaging (2015) 35:21–33. doi: 10.1111/cpf.12118

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Lintu N, Tompuri T, Viitasalo A, Soininen S, Laitinen T, Savonen K, et al. Cardiovascular fitness and haemodynamic responses to maximal cycle ergometer exercise test in children 6–8 years of age. J Sports Sci (2014) 32:652–9. doi: 10.1080/02640414.2013.845681

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Dencker M, Thorsson O, Karlsson MK, Lindén C, Wollmer P, Andersen LB. Maximal oxygen uptake versus maximal power output in children. J Sports Sci (2008) 26:1397–402. doi: 10.1080/02640410802199789

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Soininen P, Kangas AJ, Würtz P, Suna T, Ala-Korpela M. Quantitative serum nuclear magnetic resonance metabolomics in cardiovascular epidemiology and genetics. Circ Cardiovasc Gen (2015) 8:192–206. doi: 10.1161/CIRCGENETICS.114.000216

CrossRef Full Text | Google Scholar

35. Marshall WA, Tanner JM. Variations in the pattern of pubertal changes in boys. Arch Dis Child (1970) 45:13–23. doi: 10.1136/adc.45.239.13

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Marshall WA, Tanner JM. Variations in pattern of pubertal changes in girls. Arch Dis Child. (1969) 44:291–303. doi: 10.1136/adc.44.235.291

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Moore S, McKay H, Macdonald H, Nettlefold L, Baxter-Jones A, Cameron N, et al. Enhancing a somatic maturity prediction model. Med Sci Sports Exerc (2015) 47:1755–64. doi: 10.1249/MSS.0000000000000588

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Collings PJ, Westgate K, Väistö J, Wijndaele K, Atkin AJ, Haapala EA, et al. Cross-sectional associations of objectively-measured physical activity and sedentary time with body composition and cardiorespiratory fitness in mid-childhood: The PANIC study. Sports Med (2017) 47:769–80. doi: 10.1007/s40279-016-0606-x

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Väistö J, Haapala EA, Viitasalo A, Schnurr TM, Kilpeläinen TO, Karjalainen P, et al. Longitudinal associations of physical activity and sedentary time with cardiometabolic risk factors in children. Scand J Med Sci Sports (2019) 29:113–23. doi: 10.1111/sms.13315

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Koba S, Yokota Y, Hirano T, Ito Y, Ban Y, Tsunoda F, et al. Small LDL-cholesterol is superior to LDL-cholesterol for determining severe coronary atherosclerosis. J Atheroscler Thromb (2008) 15:250–60. doi: 10.5551/jat.E572

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Kyttälä P, Erkkola M, Lehtinen-Jacks S, Ovaskainen ML, Uusitalo L, Veijola R, et al. Finnish Children healthy eating index (FCHEI) and its associations with family and child characteristics in pre-school children. Publ Health Nutr (2014) 17:2519–27. doi: 10.1017/S1368980013002772

CrossRef Full Text | Google Scholar

42. Haapala EA, Wiklund P, Lintu N, Tompuri T, Väistö J, Finni T, et al. Cardiorespiratory fitness, physical activity, and insulin resistance in children. Med Sci Sports Exercise. (2020) 52(5):1144–52. doi: 10.1249/MSS.0000000000002216

CrossRef Full Text | Google Scholar

43. Camont L, Chapman MJ, Kontush A. Biological activities of HDL subpopulations and their relevance to cardiovascular disease. Trends Mol Med (2011) 17:594–603. doi: 10.1016/j.molmed.2011.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Poole DC, Behnke BJ, Musch TI. The role of vascular function on exercise capacity in health and disease. J Physiol (2021) 599:889–910. doi: 10.1113/JP278931

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Kujala UM, Mäkinen VP, Heinonen I, Soininen P, Kangas AJ, Leskinen TH, et al. Long-term leisure-time physical activity and serum metabolome. Circulation (2013) 127:340–8. doi: 10.1161/CIRCULATIONAHA.112.105551

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Van Lenten BJ, Navab M, Shih D, Fogelman AM, Lusis AJ. The role of high-density lipoproteins in oxidation and inflammation. Trends Cardiovasc Med (2001) 11:155–61. doi: 10.1016/S1050-1738(01)00095-0

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kontush A, Chantepie S, Chapman MJ. Small, dense HDL particles exert potent protection of atherogenic LDL against oxidative stress. Arterioscler Thromb Vasc Biol (2003) 23:1881–8. doi: 10.1161/01.ATV.0000091338.93223.E8

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Haapala EA, Väistö J, Ihalainen JK, Tomaselli González C, Leppänen MH, Veijalainen A, et al. Associations of physical activity, sedentary time, and diet quality with biomarkers of inflammation in children. Eur J Sport Sci (2022) 22:906–15. doi: 10.1080/17461391.2021.1892830

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Nobécourt E, Jacqueminet S, Hansel B, Chantepie S, Grimaldi A, Chapman MJ, et al. Defective antioxidative activity of small dense HDL3 particles in type 2 diabetes: relationship to elevated oxidative stress and hyperglycaemia. Diabetologia (2005) 48:529–38. doi: 10.1007/s00125-004-1655-5

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Chistiakov DA, Orekhov AN, Bobryshev YV. ApoA1 and ApoA1-specific self-antibodies in cardiovascular disease. Lab Invest (2016) 96:708–18. doi: 10.1038/labinvest.2016.56

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Juonala M, Viikari JSA, Kähönen M, Solakivi T, Helenius H, Jula A, et al. Childhood levels of serum apolipoproteins b and a-I predict carotid intima-media thickness and brachial endothelial function in adulthood. J Am Coll Cardiol (2008) 52:293–9. doi: 10.1016/j.jacc.2008.03.054

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Martinez E, Martorell J, Riambau V. Review of serum biomarkers in carotid atherosclerosis. J Vasc Surg (2020) 71:329–41. doi: 10.1016/j.jvs.2019.04.488

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Festa A, D’Agostino R Jr, Mykkänen L, Tracy RP, Hales CN, Howard BV, et al. LDL particle size in relation to insulin, proinsulin, and insulin sensitivity. the insulin resistance atherosclerosis study. Diabetes Care (1999) 22:1688–93. doi: 10.2337/diacare.22.10.1688

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Durante W. The emerging role of l-glutamine in cardiovascular health and disease. Nutrients (2019) 11:2092. doi: 10.3390/nu11092092

CrossRef Full Text | Google Scholar

55. Chen J, Zhang S, Wu J, Wu S, Xu G, Wei D. Essential role of nonessential amino acid glutamine in atherosclerotic cardiovascular disease. DNA Cell Biol (2020) 39:8–15. doi: 10.1089/dna.2019.5034

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Coqueiro AY, Rogero MM, Tirapegui J. Glutamine as an anti-fatigue amino acid in sports nutrition. Nutrients (2019) 11:863. doi: 10.3390/nu11040863

CrossRef Full Text | Google Scholar

57. Würtz P, Havulinna AS, Soininen P, Tynkkynen T, Prieto-Merino D, Tillin T, et al. Metabolite profiling and cardiovascular event risk: A prospective study of three population-based cohorts. Circulation (2015) 131:774–85. doi: 10.1161/CIRCULATIONAHA.114.013116

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Heikal L, Starr A, Hussein D, Prieto-Lloret J, Aaronson P, Dailey LA, et al. L-phenylalanine restores vascular function in spontaneously hypertensive rats through activation of the GCH1-GFRP complex. JACC Basic Transl Sci (2018) 30:3. doi: 10.1016/j.jacbts.2018.01.015

CrossRef Full Text | Google Scholar

59. Mitchell BM, Dorrance AM, Webb RC. Phenylalanine improves dilation and blood pressure in GTP cyclohydrolase inhibition-induced hypertensive rats. J Cardiovasc Pharmacol (2004) 43:758–63. doi: 10.1097/00005344-200406000-00004

PubMed Abstract | CrossRef Full Text | Google Scholar