- Open Access
Polymorphism of the PPARD Gene and Dynamic Balance Performance in Han Chinese Children
Hereditasvolume 156, Article number: 15 (2019)
Athletic performances are complex traits with heritability of ~66%. Dynamic balance is one of the most important athletic performances, and there has been little studies for it in sports genomics. The candidate PPARD gene was reported to be able to affect muscle development for balance predisposition and influence the athletic performance including skiing triumph in the Caucasian population. This study aims to investigate whether the PPARD gene is a susceptibility gene for dynamic balance performance in Han Chinese children.
A total 2244 children were recruited and their balance beam performances were measured. Five polymorphisms in the PPARD gene were genotyped through the MassARRAY Sequenom platform. Rs2016520 exerted significant association with dynamic balance performance (minor allele C, P = 0.015, Pcorrected < 0.05) and was affirmed in a meta-analysis by combining previously reported Caucasian cohorts (OR = 1.57, 95% CI = [1.30, 1.91], P < 10 -5) . Another polymorphism, rs2267668, was also significantly associated with dynamic balance performance (minor allele G, P = 0.015, Pcorrected < 0.05). In the dichotomous study, 321 cases (61% boys and 39% girls) and 370 controls (49% boys and 51% girls) in our samples were selected as representatives, and the thresholds were the mean velocity (0.737 m/s) ± standard deviation (0.264 m/s), in which rs2016520-C and rs2267668-G still remained significant (CI =1.41 [1.11~1.79], P = 0.004, Pcorrected < 0.016; CI =1.45 [1.14~1.86], P = 0.002, Pcorrected < 0.016). In different genders, consistent OR direction was observed for each variant.
Our results suggested that the PPARD gene is associated with dynamic balance performance of human being, and further studies to reveal its etiology is strongly suggested.
The heritability of athletic performance is approximately 66%, fluctuating in different sport types . Generally, sports phenotypes are currently divided into explosive power, endurance, strength, flexibility, balance and neuromuscular coordination. Up to now, overall 120 SNPs were found to be associated with athletes’ athletic performance. Among these SNPs, 77 were classified to endurance propensity and 43 to explosive/power propensity . Genetic basis of balance still remains unclear in sports genomics. Sports genomics usually focuses on athletes through case-control studies, and aims to explain the structure and function of genomic regions which might result to better athletic performances [3, 4]. Recently, attention to athletic performance has gradually expanded to lower age groups [5,6,7,8]. Notably, basic athletic skills have been developed in early childhood and mature at about 5 to 8 years old , therefore, the studies on young children's motor phenotype make sense. Till now, there are far less studies on children than on adults regarding sports genomics. The ACTN3 gene was found to be associated with children’s explosive power in Han Chinese . NOS3 and PPARGC1A were associated with cardiorespiratory endurance in children with cystic fibrosis . For children, the studied genes are also limited in explaining either explosive/power or endurance propensity, which is the same situation as that in the adults. It is essential to study the genetic basis of other athletic performance which is also of significant importance but remains unclear.
Dynamic balance is the process of maintaining or restoring a stable position while performing a task, and it is vitally important for most human locomotion . The improvement of dynamic balance can enhance a series of athletic performance including vertical jump , shuttle-run  and downhill skiing . Dynamic balance relies on rapid continuous feedback of vision, vestibular and somatosensory structures to perform coordinated neuromuscular movements . Dynamic balance is one of the physical fitness components that could be health-related or performance-related . In this article, we focused on performance-related components with all the children recruited in healthy status. Sports genomics is promising in explaining individual difference for children's dynamic balance performance.
The transcription factor peroxisome proliferator-activated receptor delta (PPARD) gene, which encodes a nuclear hormone receptor PPARδ, is among the hottest studied gene in sports genomics. PPARδ is considered to exert important functions in the process of transcriptional repression and nuclear receptor signaling . It is actively expressed in adipose tissue and skeletal muscles, mainly in type 1 (slow twitch) muscle fibers . Mice studies suggested that PPARδ plays an important role in myelination of the corpus callosum, cell differentiation, epidermal cell proliferation, and lipid metabolism [20,21,22]. Functional studies suggested that the PPARD gene was involved in muscle development and adaptive response to fitness training [23,24,25,26]. In mice skeletal muscle, targeted expression of activated PPARδ triggered adaptation from type 2 muscle fibers to type 1 muscle fibers . Type 1 muscle fiber has a beneficial effect on posture and endurance related performance, while type 2 muscle fiber favors strength and power related performance . Large cohort studies have shown that the PPARD gene is associated with overall athletic performance in the Polish populations  and athletes’ skiing triumph in the Russian populations , in which dynamic balance plays a crucial role.. Latest research on the involvement of PPARδ with circadian rhythm inhibitors showed the enhancement of exercise performance in mice via biallelic knock out of PPARδ repressors i.e. CRY1 -/- and CRY2 -/- . However, the biological mechanism of PPARD gene is not so clear.
In this study, we hypothesized that the PPARD gene might contribute to the individual variance of dynamic balance performance and tried to figure out how it could influence dynamic balance performance. In total, 2244 children were recruited, and their balance beam performances were measured. Five single nucleotide polymorphisms (SNPs) of PPARD gene including rs11571504, rs2016520, rs2267668, rs2299869 and rs3798343 were genotyped. Concisely, we conducted a quantitative trait loci (QTL) analysis, a case-control study, and a meta-analysis with one previously published data in the Caucasian population (152 cases and 610 controls)  to assess the association between PPARD SNPs and dynamic balance performance in Han Chinese children.
Physical characteristics of all recruited children were summarized in Additional file 1: Table S1. The mean age of the children was 5.36 years old and the average BMI was 16.1 kg/m2. The mean velocity walking through the 3-meter long balance beam was 0.737 ± 0.264 m/s. In addition, since a previous study observed correlations between measures of dynamic balance and lower-extremity maximal strength in healthy individuals across the lifespan, in which vertical jump height was one of the measurements to represent maximal strength , we also measured vertical jump height of each children for the further analysis. The mean vertical jump height was 17.46 ± 5.46 cm. The balance beam performance and vertical jump height of the children approximately obeyed normal distributions (Additional file 1: Figure S1).
All five SNPs were in Hardy-Weinberg equilibrium (HW_P > 0.05 in Table 1). Fig. 1 illustrated the linkage disequilibrium between rs2267668 and rs2016520. Table 1 shows the quantitative trait loci results of five SNPs in the balance beam performance. Two SNPs (rs2267668 and rs2016520) were significantly associated with dynamic balance performance after multiple testing correction. The SNP rs2016520 showed significance in association with dynamic balance performance (minor allele C, BETA = 0.021 m/s, P = 0.015, PFDR = 0.038), and the significance remained after taking gender, weight, BMI and age as covariates to eliminate their effects respectively (minor allele C, BETA = 0.021 m/s, Pgender = 0.02, Pgender_FDR = 0.05, Pweight= 0.013, Pweight_FDR = 0.032, PBMI = 0.014, P BMI-FDR = 0.035; BETA = 0.016 m/s, Page = 0.049). Another SNP, rs2267668, was also significantly associated with dynamic balance performance (minor allele G, BETA = 0.023 m/s, P = 0.015, PFDR = 0.038), and P-values maintained significant after correction by gender, weight, BMI and age (minor allele G, BETA = 0.023 m/s, Pgender = 0.015, Pgender_FDR = 0.05, Pweight= 0.008, Pweight_FDR = 0.032; BETA = 0.022 m/s, PBMI = 0.013, P BMI-FDR = 0.035; BETA = 0.019 m/s, Page = 0.024). The β value in the regression model pointed out the direction of the affecting way. The effect allele rs2016520-C could be a favorable allele favorable for balance beam performance. In a similar way, the β value of rs2267668-G also implied its status as a favorable allele.
To preclude interference from lower-extremity maximal strength on the results, independent influence of the PPARD polymorphism on dynamic balance and lower-extremity maximal strength was explored, as shown in Table 2 and 3. Significance remained in Table 2 between dynamic balance performance and the PPARD polymorphism after corrected by the covariate of vertical jump height (P < 0.05 after corrections). In Table 3, no significance was observed between vertical jump performance and the PPARD polymorphism either before or after it was corrected by the covariate of dynamic balance performance (P > 0.05). Both Table 2 and 3 supported that the PPARD polymorphism influenced dynamic balance performance rather than lower-extremity maximal strength.
We performed gender specific analysis too. Consistent OR directions were observed for each variant in boys and girls groups, which implied the similar predisposition to balance beam performance in different genders. In Table 4, boys’ results showed rs2016520-C (BETA = 0.027 m/s, P = 0.033, PBMI = 0.012, P BMI-FDR = 0.03; BETA = 0.025 m/s, Pweight = 0.04, Page = 0.03) and rs2267668-G (BETA = 0.024 m/s, PBMI = 0.013, Page = 0.049) as favorable alleles for dynamic balance performance. In Table 5, girls’ results also displayed the favor of rs2016520-C (BETA = 0.015 m/s, Pweight= 0.012; BETA = 0.013 m/s, PBMI = 0.012; BETA = 0.014 m/s, Page = 0.047) and rs2267668-G (BETA = 0.023 m/s, Pweight= 0.012; BETA = 0.021 m/s, PBMI = 0.012) on dynamic balance performance. In comparison with the gender-neutral dataset (Table 1), the significance of rs2016520-C and rs2267668-G slightly decreased in the gender specific dataset (Table 4 and Table 5), the major cause of which was the decline of the sample size resulting from gender split.
To further verify the association between the PPARD gene and balance beam performance, the mean velocity (0.737 m/s) and ± standard deviation (0.264 m/s) were used as the thresholds to define the case and the control group. Concisely, 321 children(61% boys and 39% girls)were picked up as the case group for their velocity above the mean value plus standard deviation (average velocity ± standard deviation: 1.16 ± 0.14 m/s) and 370 children (49% boys and 51% girls) were selected as the control group for their velocity below the mean value minus standard deviation (average velocity ± standard deviation: 0.37 ± 0.08 m/s). Table 6 and Table 7 demonstrates the case-control analysis results. In 321 cases and 370 controls, the two SNPs (rs2016520, rs2267668) remained significantly associated with the phenotype of balance beam performance after multiple corrections, i.e. the genetic predisposition of rs2016520-C (CI =1.41 [1.11~1.79], P = 0.004, PFDR = 0.012 in Table 6, P covariate corrections < 0.05 in Table 7) and rs2267668-G (CI =1.45 [1.14~1.86], P = 0.002, P FDR = 0.012 in Table 6, P covariate corrections < 0.05 in Table 7) was consistent with the effect in the previous QTL analysis.
Additionally, we performed meta-analysis by integrating our case- control data with similar data from one previously reported study. A genomic study on the association of the PPARD polymorphism and skiing triumph among the Russian population was selected in our meta-analysis. Skiing is considered to be a sort of sport largely depended on dynamic balance [32,33,34]. In real skiing situation, there are movements like turning, vacating and diverting in the air, which requires dynamic balance to restore a stable position while performing the skiing task, therefore, dynamic balance is one of the main factors that restricts skiing triumph . This Russian study recruited 152 athletes with competitive standards in skiing (including alpine skiing, ski jumping, cross-country skiing 5-10 km and 15–50 km) as cases and 610 non-athletes as controls, and then conducted genotyping to do association study . Similar research method as ours were seen in their study, therefore we selected this study to carry out the meta-analysis. With the inaccessibility of detailed individual information in the Russian study, it was unable to conduct covariate corrections. The two populations were not significantly heterogeneous (heterogeneity P-value = 0.1, I2 = 63% in Fig. 2), so that the fixed-effect model was used to estimate the overall effect of rs2016520-C on dynamic balance performance . Although no significant heterogeneity was seen in the two populations, the influence of differences in age, weight, training status and other factors on the analyses and results was not completely excluded. Fig. 2 shows that the favor of dynamic balance was linked to rs2016520-C allele based on the fixed-effect model(OR = 1.57, 95% CI = [1.30, 1.91], P < 10 -5). In different races, the P-values remained significant and the OR values were in the same direction. Chinese children and European athletes share consistent genetic predisposition on dynamic balance related performance. People owning rs2016520-C tends to have better dynamic balance related performance.
Our findings suggested that common SNPs in the PPARD gene could play an important role in Children’s dynamic balance performance. In general, the rs2016520-C allele could promote dynamic balance performance concluded from QTL analysis (BETA > 0, P < 0.05 after correction as shown in Table 1), case-control study (CI > 1, P < 0.05 after correction as shown in Table 6 and 7) and meta-analysis confirmation with the Caucasian population (OR = 1.57, 95% CI = [1.30, 1.91], P < 10 -5 in Fig. 2). Consistent BETA directions were observed for each variant in different genders, implying the similar predisposition to balance beam performance. In both genders, rs2016520-C favors better dynamic balance performance than rs2016520-T (BETA > 0 in Table 4 and 5).Although significant age difference were observed in other studies , significance remained in this study after taking age as a covariate (P < 0.05 in Table 1, 4, 5 and 7), which might be due to the elementary mature of basic motor skills in 5-8 years old children . The previously reported associations between measures of dynamic balance and lower-extremity maximal strength  did not interfere the effect of the PPARD polymorphism on the dynamic balance performance in our study (P < 0.05 after correction in Table 2; P > 0.05 in Table 3). Studies showed that the function of neurophysiological mechanisms during dynamic balance control and maximal strength production is task-specific, although the activation of corticospinal pathways are similar [36, 37]. Furthermore, studies of spinal and corticospinal excitability showed different activation patterns in task-specific performance. For example, during the execution of strength-related isometric ankle dorsiflexion and plantarflexion, short-latency facilitation induced by transcranial magnetic stimulation was observed , yet it was not seen during the performance of a balance-related perturbed stance task .
Rs2016520, located in the 5’ untranslated region of exon 4, is the best clinically studied polymorphism in the PPARD gene. The 39% higher transcriptional activity of the PPARD promotor was observed within the mutant rs2016520-C owners than the common T allele ones . Given that the rs2016520-T mutate into the rare C allele, the bilateral nucleotides become consensus sequence for Sp-1 binding, which might increase the affinity of Sp-1 transcriptional factor and explain the up-regulation of PPARD expression for rs2016520-C owners . Moreover, studies showed that activation of PPARδ triggered adaptation in the skeletal muscle from type 2 muscle fibers (strength and power related) to type 1 muscle fibers (postural and endurance related), enabling the mutant mice to run a longer distance in comparison with the wildtype [23, 27]. Large cohort studies in human also showed consistency in the PPARD gene regarding association with elite athlete’s performances in the European populations [28, 29]. In particular, the rs2016520-C allele was associated with the triumph of skiing in the Russian population , in which dynamic balance could play a rather crucial role. In this study, we replicated a significant association of rs2016520-C with dynamic balance performance in the Han Chinese children. After multiple corrections and further confirmation among the different populations, our results suggest the rs2016520-C as a favorable allele for dynamic balance performance.
rs2267668-G was significantly associated with dynamic balance performance concluded from QTL analysis (BETA > 0, P < 0.05 after correction in Table 1) and case-control analysis (CI > 1, P < 0.05 after correction in Table 6 and 7). In the boys’ or girls’ dataset, BETA directions were observed consistent with gender-neutral ones ( BETA > 0 for rs2267668-G, shown in Table 1, 4 and 5), although the significance decreased due to the split of the sample size. In both genders, rs2267668-G favors better dynamic balance performance. Previous studies showed that PPARD rs2267668 A/G SNP had impacts on the improvement of mitochondrial function and aerobic physical fitness [41, 42]. A blunted increase in individual anaerobic threshold (i.e. worse aerobic physical fitness) was observed in rs2267668-G, which suggested the less effectiveness of aerobic exercise training, while rs2267668-G of cultured human myotubes displayed low skeletal muscle mitochondrial function . In a physical and dietary lifestyle intervention, the minor G allele of rs2267668 was associated with less reduction in nonvisceral adipose tissue mass as well as blunted increase in relative muscle volume of the leg . The G allele of rs2267668 showed better dynamic balance performance in Chinese children in our study. The influence effect of the rs2267668 site still needs further investigation from experiments. The marginal positive result of the correlation between rs2267668 and dynamic balance performance might be caused by the linkage disequilibrium with the 5'UTR functional rs2016520.
In summary, we found the significant association of one 5’UTR variant (rs2016520-C) and one intron variant (rs2267668-G) in the PPARD gene with the balance beam performance. Our results suggest that the PPARD gene is associated with dynamic balance performance of children, and further studies to reveal its etiology is strongly suggested.
Recruited volunteers’ dataset
With the purpose of the study clearly explained to the parents, children, and all informed consent accessible to the children’s guardians, our study was scrutinized and approved by the local Ethical Committee of Human Genetics. The research was conducted in accordance with institutional requirements and the Declaration of Helsinki Principles. Overall 2244 unrelated healthy Chinese children (1155 boys and 1089 girls; average age ± standard deviation: 5.36 ± 1.03 years old) were recruited. All the children were of Han Chinese origin and lived in Shanghai. Oral epithelial cells were collected through oral swabs in identified tubes for follow-up experiments.
All the tests of balance beam performance and vertical jump height were conducted and recorded by the Shanghai University of Sports based on the standard manual for Chinese national physique determination (child part) . We made efforts to ensure the children to fully understand the process of measurements and reflect their true levels. The performance of walking through the balance beam was to evaluate children's dynamic balance performance. The balance beam was 3 meters long, 10 centimeters wide and 30 centimeters high. A 20 centimeters long, 20 centimeters wide, and 30 centimeters high board were added for step to each end of the balance beam as shown in Additional file 1: Figure S2. During the test, the children were guided to stood on the board and faced the balance beam with arms lifted horizontally. When ‘start’ was heard, the children proceeded as fast as possible. Walking through time was recorded and then the average velocity was calculated for further research. As for the vertical jump test, the Vertec apparatus was used. It was of steel frame construction with horizontal vanes which were rotated out of the way by the hand to indicate the height reached. We firstly recorded the standing height of each child with one arm fully extended upward, then instructed them to jump upwards vigorously and touch the highest possible vane. The vertical jump height was the difference between the maximal jump height and standing height. The score was recorded in centimeters. Each test allowed the children to perform three times to pick out the best score for the analysis. The same equipment under the same condition was used in each time to eliminate systematical errors.
SNPs selection and genotyping
The PPARD gene is located on chromosome 6p21.2-p21.1 and it spans about 92.62 kb of DNA and contains several SNPs. We downloaded the PPARD gene information of the Southern Han Chinese population from the 1000 Genomes website and selected 5 common tag SNPs (rs3798343, rs2299869, rs2267668, rs11571504, and rs2016520) as shown in Additional file 1: Figure S3 and Additional file 1: Table S2 through the Haploview 4.2 Software (Ensembl GRCh37 2018) . To be exact, the downloaded marker information and linkage pedigree format were firstly uploaded to the Haploview 4.2 Software. Then, major thresholds were set as listed below: HW-P value cutoff at 0.05, minimal genotype percentage at 95%, maximal Mendel errors at 1, minor allele frequency (MAF) at 0.05, r2 threshold at 0.8 . For the linkage format, pairwise comparisons of markers over 500 kb apart were set to be ignored. Afterwards, ‘Pairwise tagging only’ model was selected and the five tag SNPs were then given by the software. TIANamp Swab DNA Kit (TIANGEN) was used to extract genomic DNA from oral epithelial cells. The genomic DNA was obtained from participants’ epithelial mouth cells using high yield oral swab genomic DNA extraction kit. 2 μl DNA solution of each child was used for genotyping. All SNPs were genotyped on the Sequenom Massarray platform in Bio-x institute. SNPs were all genotyped by a call rate ≥ 95%, which suggested the failure rate of genotyping below 5%. The genotyping rate of each SNP in each studied group was listed in the relevant table. The process of Sequenom Massarray SNP genotyping involved polymerase chain reaction (PCR), Shrimp alkaline phosphatase (SAP) reaction, single base extension, resin cleanup and detection on mass spectrometry .
We used SHEsisPlus platform [47,48,49] to carry out statistical analysis, including Hardy–Weinberg equilibrium (HWE), linkage disequilibrium (LD), quantitative trait loci(QTL)and case-control analysis for 5 SNPs in the PPARD gene. QTL analyses were firstly carried out to find out the genotypes that might be related to the quantitative trait. The multiple-testing procedure was used to eliminate the duplicate effects of 5 SNPs to control the false discovery rate . The covariate of gender, weight, BMI, age and vertical jump height were further used in the regression model to eliminate their effect on the balance beam performance based on SHEsisPlus [47,48,49]. After taking these factors as covariates to eliminate their effects on dynamic balance performance, the sports genomics studies could be more meaningfully interpreted. The statistical significance was set at the threshold P < 0.05. The OR value in the case-control analysis indicates the influence direction and intensity of the PPARD loci on the balance trait. CI values of the Han Chinese children are used to estimate the range of the children's overall parameters. In addition, we searched related literature to conduct meta-analysis, which was carried out on RevMan5.3 software under Cochran guidelines . The Z-test was used to detect significance and the heterogeneity Q test was used to assess whether effects differ among different populations. In this test, the heterogeneity P > 0.05 (Fig. 2) suggests that the two populations are a homogeneous collection. Therefore, the fixed effect Mantel-Haenszel model was used.
Utah residents with Northern and Western Europe
Han Chinese in South
Benjamini and Hochberg procedure for controlling the false discovery rate of hypothesis tests
Hardy- Weinberg equilibrium
minor allele frequency
- PPARD/ PPARδ:
peroxisome proliferator-activated receptors delta
Single nucleotide polymorphism
De Moor MHM, Spector TD, Cherkas LF, et al. Genome-wide linkage scan for athlete status in 700 British female DZ twin pairs[J]. Twin Res Hum Genet. 2007;10(6):812–20.
Ahmetov II, Fedotovskaya ON. Current Progress in Sports Genomics, Makowski G S, editor. Adv Clin Chem. 2015;70:247–314.
Pitsiladis Y, Wang G, Wolfarth B. Genomics of elite sporting performance: what little we know and necessary advances (vol 47, pg 550, 2013)[J]. Br J Sports Med. 2013;47(10):656.
Wang G, Tanaka M, Eynon N, et al.: The Future of Genomic Research in Athletic Performance and Adaptation to Training, Posthumus M, Collins M, editor, Genetics and Sports, 2nd Revised and Extended Edition, 2016: 55-67.
이영란. Comparisons of Talent Children in Sports and General Students of Elementary School on Physique, Isokinetic Strength and Anaerobic Exercise Capacity[J]. The Korean Society of Sports Science, 2014, 23(6): 1225-1234.
Mueller J, Boehm B, Elmenhorst J, et al. Reduction of exercise capacity in children from summer to winter is associated with lower sporting activity: a serial study[J]. Pediatr Res. 2013;74(4):439–42.
Song J, 유성호, 정현철, et al. Comparison of Skeletal Maturation, Bone Mineral Density, Physical Fitness, Aerobic and Anaerobic Capacity in Children Who Participated in Different Sports[J]. The Korean Journal of Growth and Development, 2016, 24(1): 37-44.
Lubans DR, Morgan PJ, Cliff DP, et al. Fundamental Movement Skills in Children and Adolescents Review of Associated Health Benefits[J]. Sports Medicine. 2010;40(12):1019–35.
Malina RM, Katzmarzyk PT. Physical activity and fitness in an international growth standard for preadolescent and adolescent children[J]. Food Nutr Bull. 2006;27(4):S295–313.
Zhang Q, Cao Y, Chen J, et al. ACTN3 is associated with children's physical fitness in Han Chinese[J]. Mol Genet Genomics. 2019;294(1):47–56.
Yvert T, Santiago C, Santana-Sosa E, et al. Physical-Capacity-Related Genetic Polymorphisms in Children with Cystic Fibrosis[J]. Pediatr Exerc Sci. 2015;27(1):102–12.
Winter D A, Patla A E, Frank J S J M P T. Assessment of balance control in humans[J], 1990, 16(1-2): 31-51.
Kean CO, Behm DG, Young WB. Fixed foot balance training increases rectus femoris activation during landing and jump height in recreationally active women[J]. Journal of Sports Science and Medicine. 2006;5(1):138–48.
Yaggie J A, Campbell B M J J O S, Research C. Effects of balance training on selected skills[J], 2006, 20(2): 422-428.
Malliou P, Amoutzas K, Theodosiou A, et al. Proprioceptive training for learning downhill skiing[J]. Percept Mot Skills. 2004;99(1):149–54.
Lawrence EL, Cesar G, Bromfield MR, et al. Strength, Multijoint Coordination, and Sensorimotor Processing Are Independent Contributors to Overall Balance Ability[J]. Biomed Research International. 2015.
Caspersen C J, Powell K E, Christenson G M. Physical activity, exercise, and physical fitness: definitions and distinctions for health-related research[J]. Public health reports (Washington, D.C. : 1974), 1985, 100(2): 126-31.
Berger J, Moller DE. The mechanisms of action of PPARs[J]. Annu Rev Med. 2002;53:409–35.
Loviscach M, Rehman N, Carter L, et al. Distribution of peroxisome proliferator-activated receptors (PPARs) in human skeletal muscle and adipose tissue: relation to insulin action[J]. Diabetologia. 2000;43(3):304–11.
Lee C-H, Olson P, Hevener A, et al. PPARδ regulates glucose metabolism and insulin sensitivity[J], 2006, 103(9): 3444-3449.
Feige JN, Gelman L, Michalik L, et al. From molecular action to physiological outputs: Peroxisome proliferator-activated receptors are nuclear receptors at the crossroads of key cellular functions[J]. Prog Lipid Res. 2006;45(2):120–59.
Wang YX, Lee CH, Tiep S, et al. Peroxisome-proliferator-activated receptor delta activates fat metabolism to prevent obesity[J]. Cell. 2003;113(2):159–70.
Luquet S, Lopez-Soriano J, Holst D, et al. Peroxisome proliferator-activated receptor delta controls muscle development and oxydative capability[J]. Faseb Journal, 2003, 17(13): 2299-+.
Jones N, Kiely J, Suraci B, et al. A genetic-based algorithm for personalized resistance training[J]. Biology of Sport. 2016;33(2):117–26.
Mahoney D J, Parise G, Melov S, et al. Analysis of global mRNA expression in human skeletal muscle during recovery from endurance exercise[J]. Faseb Journal, 2005, 19(8): 1498-+.
Domanska-Senderowska D, Snochowska A, Szmigielska P, et al. ANALYSIS OF THE PPARD GENE EXPRESSION LEVEL CHANGES IN FOOTBALL PLAYERS IN RESPONSE TO THE TRAINING CYCLE[J]. Balkan Journal of Medical Genetics. 2018;21(3):19–25.
Wang YX, Zhang CL, Yu RT, et al. Regulation of muscle fiber type and running endurance by PPAR delta[J]. Plos Biology. 2004;2(10):1532–9.
Maciejewska-Karlowska A, Hanson ED, Sawczuk M, et al. Genomic haplotype within the Peroxisome Proliferator-Activated Receptor Delta ( PPARD) gene is associated with elite athletic status[J]. Scand J Med Sci Sports. 2014;24(3):E148–55.
Akhmetov II, Astranenkova IV, Rogozkin VA. Association of PPARD gene polymorphism with human physical performance[J]. Mol Biol (Mosk). 2007;41(5):852–7.
Jordan S D, Kriebs A, Vaughan M, et al. CRY1/2 Selectively Repress PPAR delta and Limit Exercise Capacity[J]. Cell Metabolism, 2017, 26(1): 243-+.
Muehlbauer T, Gollhofer A, Granacher U. Associations Between Measures of Balance and Lower-Extremity Muscle Strength/Power in Healthy Individuals Across the Lifespan: A Systematic Review and Meta-Analysis[J]. Sports Medicine. 2015;45(12):1671–92.
Hongguang YYMA, Yantao L, et al. Relationship between Lower Limb Muscle Strength and Static Balance Ability of National Elite Men's Freestyle Skiing Aerials Athletes[J]. Journal of Shenyang Sport University. 2012.
Moxnes JF, Sandbakk O, Hausken K. A simulation of cross-country skiing on varying terrain by using a mathematical power balance model[J]. Open access journal of sports medicine. 2013;4:127–39.
Zhang Y, Kong L. The function of body banlance in downhill skiing [J]. Snow and ice sports. 2001;3:28–9 (in Chinese).
Furlan AD, Pennick V, Bombardier C, et al. 2009 Updated Method Guidelines for Systematic Reviews in the Cochrane Back Review Group[J]. Spine. 2009;34(18):1929–41.
Schubert M, Beck S, Taube W, et al. Balance training and ballistic strength training are associated with task-specific corticospinal adaptations[J]. European Journal of Neuroscience. 2008;27(8):2007–18.
Beck S, Taube W, Gruber M, et al. Task-specific changes in motor evoked potentials of lower limb muscles after different training interventions[J]. Brain Res. 2007;1179:51–60.
Morita H, Olivier E, Baumgarten J, et al. Differential changes in corticospinal and Ia input to tibialis anterior and soleus motor neurones during voluntary contraction in man[J]. Acta Physiol Scand. 2000;170(1):65–76.
Taube W, Schubert M, Gruber M, et al. Direct corticospinal pathways contribute to neuromuscular control of perturbed stance[J]. J Appl Physiol. 2006;101(2):420–9.
Skogsberg J, Kannisto K, Cassel TN, et al. Evidence that peroxisome proliferator-activated receptor delta influences cholesterol metabolism in men[J]. Arteriosclerosis Thrombosis and Vascular Biology. 2003;23(4):637–43.
Stefan N, Thamer C, Staiger H, et al. Genetic variations in PPARD and PPARGC1A determine mitochondrial function and change in aerobic physical fitness and insulin sensitivity during lifestyle intervention[J]. J Clin Endocrinol Metab. 2007;92(5):1827–33.
Thamer C, Machann J, Stefan N, et al. Variations in PPARD determine the change in body composition during lifestyle intervention: A whole-body magnetic resonance study[J]. J Clin Endocrinol Metab. 2008;93(4):1497–500.
General Administration of Sport of China. Standard manual for Chinese national physique determination (child part) [M]. Peoples Sports Publishing House, 2003. (in Chinese)
Barrett JC, Fry B, Maller J, et al. Haploview: analysis and visualization of LD and haplotype maps[J]. Bioinformatics. 2005;21(2):263–5.
De Bakker PIW, Burtt NP, Graham RR, et al. Transferability of tag SNPs in genetic association studies in multiple populations[J]. Nat Genet. 2006;38(11):1298–303.
Gabriel S, Ziaugra L, Tabbaa D. SNP Genotyping Using the Sequenom MassARRAY iPLEX Platform[J]; 2009.
Shi YY, He L. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci[J]. Cell Res. 2005;15(2):97–8.
Li Z, Zhang Z, He Z, et al. A partition-ligation-combination-subdivision EM algorithm for haplotype inference with multiallelic markers: update of the SHEsis (http://analysis.bio-x.cn)[J]. Cell Res. 2009;19(4):519–23.
Shen J, Li Z, Chen J, et al. SHEsisPlus, a toolset for genetic studies on polyploid species[J]. Sci Rep. 2016;6.
Benjamini Y, Hochberg Y J J O T R S S. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing[J], 1995, 57(1): 289-300.
All the participants are appreciated for their contribution and cooperation .
This work was supported by the 973 Program (2015CB559100), the National Key R&D Program of China (2016YFC1306903, 2017YFC0908105), the Natural Science Foundation of China (81421061, 81701321, 31571012, 81501154), the Program of Shanghai Subject Chief Scientist (15XD1502200), the National Program for Support of Top-Notch Young Professionals, Shanghai Key Laboratory of Psychotic Disorders (13dz2260500), the National Program for Support of Top-Notch Young Professionals to Y.S., Shanghai Hospital Development Center (SHDC12016115), Shanghai Mental Health Center (2016-fx-02), Shanghai Science and Technology Committee (17JC1402900, 17490712200) and shanghai municipal health commission (ZK2015B01, 201540114)
Availability of data and materials
The original datasets used and analyzed during the current study are available from the coresponding author on reasonable request.
Ethics approval and consent to participate
With the purpose of the study clearly explained to the parents,children and all informed consent accessible to children’s guardians, our study was scrutinized and approved by the local Ethical Committee of Human Genetics.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Physical characteristics of all recruited children. Table S2. Genetic information of five markers. Figure S1. Distribution of the balance beam and vertical jump performance in the Han Chinese children. Figure S2. Diagram of balance beam performance. Figure S3. Location of the five single nucleotide polymorphisms (SNPs) in the PPARD gene. (DOCX 365 kb)