Journal of Diabetes Research
Volume 2015 (2015), Article ID 538732, 7 pages
http://dx.doi.org/10.1155/2015/538732

Research Article

Influence of Cardiorespiratory Fitness on PPARG mRNA Expression Using Monozygotic Twin Case Control

Marcos Roberto Queiroga,¹ Ricardo Augusto Barbieri,² Sandra Aires Ferreira,¹ André Ducati Luchessi,³ Rosario Dominguez Crespo Hirata,⁴ Mario Hiroyuki Hirata,⁴ and Eduardo Kokubun²

¹Departamento de Educação Física da Universidade Estadual do Centro-Oeste (UNICENTRO), Rua Simeão Camargo Varela de Sá, 03, 85040-080 Guarapuava, PR, Brazil
²Departamento de Educação Física da Universidade Estadual Paulista (UNESP), Avenida 24 A 1515, 13506-900 Rio Claro, SP, Brazil
³Departamento de Análises Clínicas e Toxicológicas da Universidade Federal do Rio Grande do Norte (UFRN), Rua General Gustavo Cordeiro de Farias, 384, 59012-570 Natal, RN, Brazil
⁴Departamento de Análises Clínicas e Toxicológicas da Universidade de São Paulo (USP), Avenida Professor Lineu Prestes, 580, 05508-000 São Paulo, SP, Brazil

Received 22 October 2014; Revised 22 February 2015; Accepted 1 March 2015

Academic Editor: Norman Cameron

Copyright © 2015 Marcos Roberto Queiroga et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The influence of cardiorespiratory fitness (VO₂max) on anthropometric variables and PPARG mRNA expression was investigated. Monozygotic twin pairs aged 11–18 years were grouped into discordant (D) and concordant (C) high and low VO₂max groups. VO₂max was determined by progressive maximal exercise test on treadmill with gas exchange analysis. Body mass (BM), BMI, waist circumference (WC), triceps (TR), and subscapular (SB) skinfold thicknesses were measured. Twins from the discordant group had differences in VO₂max values (D-high = versus D-low =  mL·kg⁻¹·min⁻¹, ), while no differences were found in the concordant group (C-high = versus C-low =  mL·kg⁻¹·min⁻¹, ). In discordant group, VO₂max was negatively correlated with TR + SB (, ) and positively correlated with PPARG expression in leukocytes (, ). Moreover, PPARG expression was directly correlated with BM (, ) and height (, ). In concordant twins, VO₂max was inversely correlated with BM (, ), BMI (, ), WC (, ), and TR + SB (, ). Twins D-high had 1.78-fold greater PPARG expression when compared with twins D-low (). In conclusion, the cardiorespiratory fitness may modulate PPARG expression in childhood and adolescence, independently of the genetic background.

1. Introduction

Peroxisome proliferator-activated receptors (PPARs) are involved in the regulatory response processes of lipid and glucose metabolisms, adipocyte differentiation, and inflammatory response [1]. The three peroxisome proliferator-activated receptors (PPARs) α, β (δ), and γ constitute a distinct subfamily in the superfamily of nuclear receptors. PPARs are ligand-regulated transcription factors that control gene expression by binding to specific response elements within promoters. Studies have shown that PPARs are involved in the physiopathology of complex metabolic disorders, such as those related to metabolic syndrome, resulting in atherosclerosis and cardiovascular diseases [2].

PPARG regulates the expression of numerous genes involved in controlling cellular energy homeostasis, lipid metabolism, insulin action, and adipocyte differentiation [3]. In a study using cell culture, Uruno et al. [4] showed the inhibitory effect of PPARG agonists, controlling blood pressure, on angiotensin-II induced aldosterone synthase expression and aldosterone secretion. An association of PPARG gene variants with insulin resistance, type 2 diabetes, obesity, and hypertension has been shown [5, 6]. Additionally, PPARG activation is associated with increased consumption of free fatty acids by skeletal muscle and adipose tissue, inhibition of the expression of proinflammatory genes (TNFα) that trigger atherosclerosis, and increased insulin sensitivity [7]. Atherosclerosis is usually preceded by endothelial dysfunction, whereas agonists have been reported to improve the function of these cells in type 2 diabetic patients [8] and nondiabetic patients with coronary artery disease [9].

Besides the role played by PPARG activation in regulating the metabolism of glucose and lipids, a possible regulatory role in the expression of proinflammatory adipokines has also been suggested [1, 6]. Additionally, the renal protective effect of PPARG agonists against nondiabetic renal disease [10] and anticancer activities of PPARG agonist was recently reported [11].

It has been also shown that PPARG variant (PPARG C1A) is related to effects on the effectiveness of aerobic exercise training to increase aerobic fitness and insulin sensitivity [12]. Moreover, results suggest that the mitochondrial function associated with aerobic fitness and insulin resistance is deeply affected by the expression of coactivators (PGC-1α and PGC-1β) of the PPARG [13]. Therefore, PPARG may be a promising candidate for regulation by maximal oxygen uptake [14].

The cardiorespiratory fitness as defined by maximal oxygen uptake (VO₂max) is modulated largely by physical activity level and genotype. A low VO₂max has been shown to be associated with glucose intolerance, insulin resistance, type-2 diabetes mellitus [15, 16], and downregulation of the PPARG mRNA expression from skeletal muscle and peripheral mononuclear cells [17, 18]. However, it has not been demonstrated whether association between VO₂max and PPARG mRNA expression is influenced by familial or genetic background.

Nevertheless, studies involving the analysis of the interaction between physical activity and genetic aspects have increased considerably in recent years, but knowledge about the impact or effect of physical exercise on genome or highlighting the limitations of genetic background in physical performance is preliminary. It has been argued that the human genome has evolved over a period when high levels of physical activity were essential for survival, which supports the view that there is a link between exercise and gene expression regulation [19].

The direct influence of VO₂max on PPARG expression has not been previously examined, regardless of genetic background, but it can be evaluated using the discordant monozygotic (MZ) twin model (cotwin case control study design). This experimental model allows estimating the effect of one discordant environmental factor [20, 21]. MZ twins have the same inherited genes, and if they differ in a particular trait, the difference can be considered to be due to environmental factors [22]. Using this model, the influence of VO₂max can be estimated, isolated from additional covariates that usually confound other studies. The aim of this study was to investigate the influence of VO₂max on PPARG mRNA expression using MZ twin model case control.

2. Materials and Methods

2.1. Subjects

Ninety-eight twin pairs who attended private and public schools of Rio Claro, São Paulo, were invited to participate in this study. Among the twin pairs, 53 were girls and 45 were boys aged from 11 to 18 years. In 31 twin pairs, one or both twins declined participation and 13 pairs were not found. Fifty-four twin pairs (35 girls and 19 boys) attended clinical and evaluation sessions. Parental testing by microsatellite analysis was performed as previously described [23] and revealed that the sample evaluated had 38 MZ twin pairs and 16 dizygotic (DZ) twin pairs. In the study, only MZ twin pairs were included (24 girls and 14 boys, years old;  kg;  cm). The cotwin control design is a powerful research methodology for studying the effects of environmental risk factors on the development of disease. In the case-control study, disease discordant pairs are examined for differences in antecedent exposures [24].

MZ twin pairs, their parents, and/or guardians were previously informed about the experimental procedures and provided written consent for participation. The intervention protocols were approved by the Research Ethics Committee (protocol no. 5093) of the São Paulo State University (UNESP), according to the norms of Resolution 196/96 of the National Health Council on research involving humans.

2.2. Anthropometric and Physical Fitness Measurements

Body mass and height were measured barefoot in light clothing to the nearest 100 g and 0.1 cm, respectively. Body mass index (BMI) was calculated by dividing body mass (kg) by height (m²). Waist circumference was measured in triplicate, midway between the lowest rib and the superior border of the iliac crest, with a flexible tape. The average of the 3 waist circumference measurements was used in all analyses. Skinfold thicknesses were measured at the triceps (TR) and subscapular (SB) on the right side of the body with a Harpenden skinfold calipers (Holtain Ltd., Bryberian, UK). Three sets of measurements were taken to the nearest 0.2 mm at each site and the mean of the three values was used. Anthropometric measurements were performed according to conventional procedures [25].

Work-conducted maximal exercise test with gas exchange analysis (spiroergometry) was performed using treadmill model ATL Super (Inbrasport, Porto Alegre, RS, Brazil) with 1% inclination, in the morning and afternoon periods from 9:00 a.m. to 11:30 a.m. and from 2:00 p.m. to 6:00 p.m. in a room temperature chamber maintained between 20 and 25°C. After a 5-minute adaptation period with treadmill at different speeds (4 to 7 km/h), the individuals remained five minutes on the treadmill at rest in standing position. At rest and during exercise test, respiratory minute volume (VE), oxygen intake (VO₂), and carbon dioxide production (VCO₂) were continuously recorded by assessing pulmonary gas exchanges (metabolic analyzer VO2000 MedGraphics Medical Graphics Corp., St. Paul, MN). The equipment was calibrated prior to the beginning of each individual test. Then, the test was performed using protocol that provided initial speed of 4 km/h with progressive increase in workload of 1 km/h per minute until exhaustion (19-20/20 on the Borg scale for perceived exertion or respiratory quotient (RQ = VCO₂/VO₂) of at least 1.1). Verbal encouragement was given in the attempt to obtain maximum physical exertion.

VO₂max was collected breath by breath and the value adopted for analysis of data was recorded as the average oxygen uptake during 30 s that preceded test interruption. Heart rate was obtained by the Polar S810 heart rate monitor (Polar Electro, Finland). To apply the model of studies with MZ twin pairs (case-control model), the absolute intrapair difference in VO₂max ≥ 10 mL·kg⁻¹·min⁻¹ was used to define twin pairs discordant for VO₂max.

2.3. Blood Samples

Blood samples were collected from an antecubital vein using vacutainer tubes (Vacutainer Becton Dickinson Company, Plymouth, UK) exactly one week after cardiorespiratory fitness assessment (VO₂max), from the twin pairs under their parents’ care and authorization, at 7.30 a.m. after overnight fasting (10 to 12 h) and 30 min resting.

Approximately 150 μL of blood was transferred to QIAcard FTA Spots (Qiagen, Valencia, USA) for DNA extraction and parental testing by microsatellite analysis. For mRNA expression, 4 mL of blood was transferred to Paxgene tubes (Qiagen, Valencia, USA) to further total RNA extraction.

2.4. Total RNA Isolation

Total RNA was isolated immediately after blood collection according to manufacturer’s instructions with minor modifications using PAXgene Blood RNA kit (Qiagen, Valencia, USA). The entire procedure was carried out at room temperature. RNA purity was determined by spectrophotometry (A260/A280 ratio) using the equipment NanoDrop (NanoDrop Technologies, Inc., Wilmington, USA). RNA integrity was determined on BioAnalyzer (Agilent Technologies Inc., Santa Clara, USA) and RNA was regarded as intact when showing two distinct bands for 28S and 18S RNA.

2.5. Gene Expression Analysis by Real-Time PCR

The PPARG mRNA expression was measured by quantitative real-time PCR (q-PCR). cDNA was synthesized using the High Capacity Reagent (Applied Biosystems, Foster City, USA), according to manufacturer’s protocol using the thermocycler Veriti 96 Well Thermal Cycler (Applied Biosystems, Foster City, USA), and it was stored at −20°C.

The cDNA from total RNA of peripheral blood leukocytes was amplified by qPCR using TaqMan amplifying system and ABI 7500 equipment (Applied Biosystems, Foster City, USA). Primers and fluorophore-labeled probes specific for PPARG and reference gene glyceraldehyde-3 phosphate dehydrogenase (GAPD) (Table 1) were designed using Primer Express (Applied Biosystems, Foster City, USA) and Primer Premier 5.0 (Premier Biosoft, Palo Alto, USA) software.

Table 1: Primers and probes used for mRNA expression by qPCR.

Gene
Primers
Amplicon size

PPARG   
5′ GCAAAGGCGAGGGCG 3′
82 bp

5′ CCCATCATTAAGGAATTCATGTCAT 3′

5′ FAM CAGACAAATCACCATTCG MGB 3′

GAPD
5′ GGAAGGTGAAGGTCGGAGTCA 3′
229 bp

5′ CTGGAAGATGGTGATGGGATTTC 3′

5′ VIC TCAGCCTTGAGGGTGC MGB 3′

Note: bp: base par.

qPCR was performed under the following conditions: 20 μL of final reaction volume containing TaqMan master mix 1 X (AmpliTaq DNA Polymerase, dNTPs with dUTP, AmpErase UNG) (Applied Biosystems, Foster City, USA), 300 nM of primers, and 200 nM of labeled probes. The PCR amplification program consisted of (1) one cycle of 2 min at 50°C (UNG activation), (2) one cycle of 5 minutes at 95°C (UNG inactivation); (3) 40 cycles of 15 sec at 95°C (denaturation); (4) 1 min at 60°C (hybridization and extension). The fluorescence signals emitted by fluorophores of TaqMan probes were detected by ABI Prism 7500 (Applied Biosystems, Foster City, USA). Data were analyzed using the System 7500 SDS software (Applied Biosystems, Foster City, USA) that generates half-log curves of amplification signals and Ct values. The Ct values were used to calculate the relative mRNA expression of PPARG compared to GAPD (reference gene) using the ΔCt. It is a calculated formula: .

2.6. Parental Testing by Microsatellite Analysis

Total DNA was extracted from whole blood samples in QIAcard FTA Spots with QIAamp 96 DNA Blood Kit (QIAGEN, Valencia, USA) according to manufacturer’s recommendations. Microsatellite analysis was carried out by PCR using the Identifiler commercial kit, according to the manufacturer’s instructions (Applied Biosystems, Foster City, EUA). This assay allows the genotyping of 15 multiple informative genetic markers (CSF1PO, D2S1338, D3S1358, D7S820, D8S1179, D13S317, D16S539, D18S51, D19S433, D21S11, D5S51, FGA, TH01, TPOX, VWF, and AMELXloci), which allows for the identification of monozygotic and dizygotic twins [23].

2.7. Statistical Analyses

Data were analyzed with the use of descriptive statistics as means and standard deviation. Normality of variables was assessed by Shapiro-Wilks test. For normally distributed variables, one-way analysis of variance for multiple comparisons followed by a Bonferroni’s post hoc test was used to compare the means (anthropometric traits) between twin pairs with high and low VO₂max and Wilcoxon’s signed rank test for intrapair comparison of nonnormally distributed data (gene expression). Spearman’s rank coefficient was used to investigate correlations between VO₂max-discordant and concordant MZ twins. All statistical analyses were conducted using SPSS software version 20 (IBM SPSS), and alpha was set at .

3. Results

Twin pairs were distributed into discordant (10 girls and 8 boys aged years old) and concordant (38 girls and 20 boys aged years old) for VO₂max, and within each group, cotwins were divided into high and low VO₂max. Comparison of anthropometric characteristics for VO₂max-discordant and concordant pairs is shown in Table 2. There were no differences (between and within group) in age and anthropometry characteristics (body mass, height, BMI, waist circumference, and sum of skinfold TR + SB). As expected, cotwins from the discordant group (D-high versus D-low) showed difference in VO₂max values ( versus  mL·kg⁻¹·min⁻¹, ), while cotwins from the concordant group (C-high versus C-low) did not show difference in VO₂max ( versus  mL·kg⁻¹·min⁻¹, ). The analysis also revealed that discordant low cotwins (D-low) showed marginal difference in VO₂max values compared for concordant high twins (C-high) ( versus  mL·kg⁻¹·min⁻¹, ) (Table 2).

Table 2: Anthropometric characteristics of discordant and concordant MZ twins pairs for VO₂max.

Variables
Discordant pairs ()†
Concordant pairs ()†

D-high
D-low
C-high
C-low

Age (years)
13.9 ± 2.2
14.2 ± 2.1
—
—

(♂/♀)
18 (8/10)
58 (20/38)
—
—

VO₂max
45.9 ± 10.0
32.4 ± 10.6a
42.4 ± 9.2b
38.8 ± 9.8
3.690
0.016

Body mass (kg)
46.4 ± 9.0
46.2 ± 8.7
52.1 ± 14.1
53.7 ± 13.6
1.291
0.284

Height (cm)
155.7 ± 11.5
156.4 ± 11.0
158.8 ± 10.6
159.1 ± 10.6
0.335
0.800

BMI (kg/m²)
18.9 ± 1.4
18.7 ± 1.5
20.4 ± 4.2
21.0 ± 4.4
1.290
0.284

Waist (cm)
66.5 ± 4.4
65.7 ± 5.0
70.6 ± 11.7
71.8 ± 10.9
1.218
0.310

TR + SB (mm)
18.9 ± 6.2
19.4 ± 6.1
28.2 ± 16.6
29.7 ± 17.0
1.924
0.124

VO₂max: maximal oxygen uptake (mL·kg⁻¹·min⁻¹); BMI: body mass index; TR + SB: sum of skinfold triceps and subscapular. D-high and D-low: discordant pairs; C-high and C-low: concordant pairs. †VO₂max intrapair difference ≥10 mL·kg⁻¹·min⁻¹; data are shown as mean ± SD and compared by ANOVA test, with Bonferroni test for multiple comparisons. aP < 0.025, D-low versus D-high; bP < 0.050, D-low versus C-high.

Correlation analysis for VO₂max and PPARG mRNA expression with anthropometric characteristics in MZ twin’s pairs discordant and concordant for VO₂max is shown in Table 3. The VO₂max was negatively correlated with the sum of skinfolds (−0.540, ) and it had a strong positive correlation with PPARG mRNA expression (, ) in discordant MZ twins. Additionally, PPARG expression was positively correlated with body mass (, ) and height (, ).

Table 3: Correlation of VO₂max and PPARG mRNA expression with anthropometric characteristics in MZ twins pairs.

Variables
Discordant pairs ()
Concordant pairs ()

VO₂max
PPARG
VO₂max
PPARG

Body mass (kg)
0.064
0.714**
−0.290*
−0.163

Height (cm)
0.086
0.762*
0.126
0.135

BMI (kg/m²)
−0.109
0.214
−0.472**
−0.172

Waist (cm)
0.046
0.323
−0.426**
−0.201

TR + SB (mm)
−0.540*
−0.357
−0.739**
−0.243

VO₂max/PPARG
0.952**
0.309

VO₂max: maximal oxygen uptake (mL·kg⁻¹·min⁻¹); BMI: body mass index; TR + SB: sum of skinfold triceps and subscapular. Spearman’s rank correlation tests **P = 0.01 level (2-tailed); *P = 0.05 level (2-tailed).

The concordant MZ twins demonstrated negative correlation of the VO₂max with body mass (, ), BMI (, ), waist circumference (, ), and sum of skinfolds (, ). However, in MZ group, there was no correlation between PPARG mRNA expression with anthropometric characteristics and VO₂max (Table 3).

Regarding the PPARG mRNA expression, it was observed that discordant cotwins with the high VO₂max values (D-High) had 1.78-fold greater PPARG mRNA expression when compared to cotwins with the low VO₂max values (D-low) () (Figure 1). In the concordant group, cotwins with low VO₂max values (C-low) had 1.22-fold greater PPARG mRNA expression when compared to cotwins with high VO₂max (C-high). However, there was no significant difference () (Figure 1).

Figure 1: PPARG mRNA expression in blood total leukocytes from VO₂max discordant and concordant MZ twins. mRNA expression is shown as ΔCt values in relation to GAPD (reference gene). Results between low and high VO₂max within discordant (D-low versus D-high) and concordant (C-low with C-high) pairs were compared by paired Wilcoxon test.

4. Discussion

The major finding from this study indicates that VO₂max is capable of affecting PPARG mRNA expression in VO₂max-discordant MZ twins regardless of genetic background. The results showed 30.6% higher VO₂max values in cotwins discordant for VO₂max (D-High) compared with cotwins with low VO₂max values (D-Low), but this difference is not associated with differences in anthropometric variables between twin’s pars discordant and concordant.

The correlation analysis was performed for anthropometric characteristics with VO₂max and PPARG mRNA expression in MZ twin’s pairs with discordant and concordant for VO₂max. According to results for VO₂max-discordant MZ twins, the VO₂max values were negatively correlated with sum of skinfolds and strong positive correlation with PPARG mRNA expression. The PPARG mRNA expression was positively correlated with body mass and height. For VO₂max-concordant MZ twins, the VO₂max values were negatively correlated with body mass, BMI, waist circumference, and sum of skinfolds. Interestingly for this group, there was no correlation between PPARG mRNA expression with anthropometric characteristics and VO₂max.

In previous studies, there was correlation between obesity (BMI) and VO₂max (, ) in young subjects [26] and a strong positive correlation between the ratio of PPARG (splice variants y2/y1) mRNA and the BMI (, ) in obese patients [27]. The expression of PPARG mRNA (splice variants y2) is increased in adipose tissue in obese men and women, and their direct correlation with BMI is unknown but it may be related to the expansion of adipose tissue mass [27].

The PPARG mRNA expression in circulating monocytes from healthy individuals was increased after participating in a program of cycling exercise that caused a significant increase in VO₂max [18]. Significant increase in leukocyte mRNA expression for PPARG was observed in healthy but previously sedentary individuals that participated in an 8-week low-intensity exercise program consisting of walking 10,000 steps, three times a week [28]. We have demonstrated that children and adolescents with high VO₂max showed higher PPARG mRNA expression. This result was reinforced by the strong correlation between VO₂max and PPARG mRNA expression. Although the cardiorespiratory fitness of the MZ twin pairs was evaluated only once (cross-sectional study), the method used to measure VO₂max (indirect calorimety based in gas exchange) is more accurrate than the physical activity evaluation using questionnaires and surveys, which have been questioned [29]. The absolute value of VO₂max is one of the indices of an individual’s cardiorespiratory fitness to transport oxygen to working muscles. Additionally the evidence suggests that physically active children have significantly better cardiorespiratory fitness levels than inactive children [30, 31].

In an animal study, the effect of acute and endurance exercise training on mRNA expression pattern of the different PPARs and ppara coactivator-1 alpha (PGC-1α) in muscles that largely rely on either glycolytic (plantaris) or oxidative (soleus) metabolism was analyzed [17]. Moreover, pparγ mRNA levels were the lowest in skeletal muscle which is the most responsive to changes in physical activity levels.

In spite of the important contribution of this study, there are some limitations such as the cut-off values of the discordant cardiorespiratory fitness and timing of blood sampling for RNA extraction. The cardiorespiratory fitness (based on studies of familial aggregation and twins) shows significant genetic effects (50–67%) [32, 33]; our results should be interpreted with caution. Among the assumptions of the study with twins, one of them admits that the phenotype carries great genetic background when showing low variance between MZ twins. In classical studies with MZ twin pairs, discordance is defined using criteria based on literature in which phenotype is a condition usually acquired [34]. However, even in studies in which variance allows establishing discordance between identical twins (i.e., one twin is diabetic/obese/hypertensive/sedentary and the other is not), environmental and/or genetic factors that may have significantly influenced this discrepancy are still questioned [35, 36]. In addition, it was observed that the behavior of the cardiorespiratory fitness is unstable during life, being influenced by several factors such as age, gender, and physical activity level. It should be emphasized that, besides the improvement promoted by physical exercise on VO₂max, the period corresponding to childhood and adolescence is marked by growth and gradual development of the cardiorespiratory system components that determine VO₂max (lungs, heart, and muscles) [37]. In this sense, possible discrepancies in VO₂max in this age group in genetically identical individuals would be more likely to be due to environmental factors (physical activity and/or exercise). In other words, the determinants that distinguish aerobic fitness of a MZ twin pair compared to his cotwin could be due to environmental factors because, in this model, the changes observed in VO₂max as a function of growth and development of the various systems would be under genetic control. Considering the increase in VO₂max caused by growth, reducing aerobic capacity at this stage is less likely than increasing it. Thus, ignoring errors from evaluation, evaluator, and measurement equipment, it is possible that discordant cotwin pair with higher VO₂max in this study is actually more active than his correspondent cotwin.

We should also highlight that differences in the PPARG mRNA expression observed between discordant MZ twins (D-high versus D-low) may be influenced by the timing of blood sample collection. However it is noteworthy that concordant twins (C-high versus C-low) underwent the same protocol, but did not show variability in PPARG expression. There is evidence showing that the PPARG mRNA expression increases significantly (1.5- to 2.5-fold) up to 3 hours after exercise, but returning to basal levels within 24 h of exercise [18]. Reinforcing this information, it was observed that the mRNA of the PPAR family of transcription factors is sensitive to exercise training in skeletal muscle but is not particularly sensitive to acute bouts of endurance exercise. These findings would suggest that changes in the PPAR mRNA expression are most likely mediated by the accumulation of repeated bouts of exercise as opposed to single bouts of activity [17].

Regarding the cut-off value of 10 mL·kg⁻¹·min⁻¹ adopted to establish discordance between MZ twin pairs in this study. In a longitudinal study with MZ and DZ twin pairs discordant for cardiorespiratory fitness (VO₂ peak of versus  mL·kg⁻¹·min⁻¹), the researchers used lower cut-off values (≈6 mL·kg⁻¹·min⁻¹) and emphasized that the differences within twin pairs were sufficient to identify metabolic disorders [38]. In another study, researchers adopted 9% as cut-off value to determine discordances (mean of %) in the VO₂max values of MZ twins. The results revealed that the more and less active twins showed, respectively, VO₂max values of versus  mL·kg⁻¹·min⁻¹ [39]. It is noteworthy that the relative difference between twins with higher and lower VO₂max (discordant) in this study ranged from 16.9% to 42.1%, with mean of % (data not shown). It is also noteworthy that the absolute average of these differences was  mL·kg⁻¹·min⁻¹ for discordant pairs and approximately 7.5 mL·kg⁻¹·min⁻¹ for the study of Hannukainen et al. [39].

In conclusion, the results of this study suggest that the cardiorespiratory fitness evaluated by VO₂max parameter may modulate the PPARG mRNA expression in childhood and adolescence, independently of the genetic background. Considering that the PPARG expression regulates lipid and glucose metabolism, insulin action and adipocyte differentiation the effect of cardiorespiratory fitness on PPARG may contribute to the prevention of future metabolic diseases in adults such as type 2 diabetes.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

The authors contributed equally to this work.

Acknowledgments

The authors thank all the volunteers who participated in this study and E. Kokubum, M. H. Hirata, and R. D. C. Hirata are recipients of fellowships from CNPq, Brasilia, Brazil. This paper was supported by Foundation for Research Support of São Paulo/SP/Brazil (FAPESP).

References

J. Berger and D. E. Moller, “The mechanisms of action of PPARs,” Annual Review of Medicine, vol. 53, pp. 409–435, 2002.

B. Desvergne, L. Michalik, and W. Wahli, “Be fit or be sick: peroxisome proliferator-activated receptors are down the road,” Molecular Endocrinology, vol. 18, no. 6, pp. 1321–1332, 2004.

J. Auwerx, “PPARgamma, the ultimate thrifty gene,” Diabetologia, vol. 42, no. 9, pp. 1033–1049, 1999.

A. Uruno, K. Matsuda, N. Noguchi et al., “Peroxisome proliferator-activated receptor-γ suppresses CYP11B2 expression and aldosterone production,” Journal of Molecular Endocrinology, vol. 46, no. 1, pp. 37–49, 2011.

A. Sugawara, A. Uruno, M. Kudo, K. Matsuda, C. W. Yang, and S. Ito, “PPARγ agonist beyond glucose lowering effect,” The Korean Journal of Internal Medicine, vol. 26, no. 1, pp. 19–24, 2011.

I. Barroso, M. Gurnell, V. E. F. Crowley et al., “Dominant negative mutations in human PPARgamma associated with severe insulin resistance, diabetes mellitus and hypertension,” Nature, vol. 402, no. 6764, pp. 880–883, 1999.

C. Jiang, A. T. Ting, and B. Seed, “PPAR-γ agonists inhibit production of monocyte inflammatory cytokines,” Nature, vol. 391, no. 6662, pp. 82–86, 1998.

F. M. A. C. Martens, F. L. J. Visseren, E. J. P. de Koning, and T. J. Rabelink, “Short-term pioglitazone treatment improves vascular function irrespective of metabolic changes in patients with type 2 diabetes,” Journal of Cardiovascular Pharmacology, vol. 46, no. 6, pp. 773–778, 2005.

C. Staniloae, V. Mandadi, D. Kurian et al., “Pioglitazone improves endothelial function in non-diabetic patients with coronary artery disease,” Cardiology, vol. 108, no. 3, pp. 164–169, 2007.

B. H. Chung, S. W. Lim, K. O. Ahn et al., “Protective effect of peroxisome proliferator activated receptor gamma agonists on diabetic and non-diabetic renal diseases,” Nephrology, vol. 10, supplement 2, pp. S40–S43, 2005.

C. Blanquicett, J. Roman, and C. M. Hart, “Thiazolidinediones as anticancer agents,” Cancer Therapy, vol. 6, pp. 25–34, 2008.

N. Stefan, C. Thamer, H. Staiger et al., “Genetic variations in PPARD and PPARGC1A determine mitochondrial function and change in aerobic physical fitness and insulin sensitivity during lifestyle intervention,” The Journal of Clinical Endocrinology & Metabolism, vol. 92, no. 5, pp. 1827–1833, 2007.

J. St-Pierre, J. Lin, S. Krauss et al., “Bioenergetic analysis of peroxisome proliferator-activated receptor γ coactivators 1α and 1β (PGC-1α and PGC-1β) in muscle cells,” The Journal of Biological Chemistry, vol. 278, no. 29, pp. 26597–26603, 2003.

A. Lucia, F. Gómez-Gallego, I. Barroso et al., “PPARGC1A genotype (Gly482Ser) predicts exceptional endurance capacity in European men,” Journal of Applied Physiology, vol. 99, no. 1, pp. 344–348, 2005.

B. Nyholm, A. Mengel, S. Nielsen et al., “Insulin resistance in relatives of NIDDM patients: the role of physical fitness and muscle metabolism,” Diabetologia, vol. 39, no. 7, pp. 813–822, 1996.

M. Wei, L. W. Gibbons, T. L. Mitchell, J. B. Kampert, C. D. Lee, and S. N. Blair, “The association between cardiorespiratory fitness and impaired fasting glucose and type 2 diabetes mellitus in men,” Annals of Internal Medicine, vol. 130, no. 2, pp. 89–96, 1999.

E. E. Spangenburg, D. A. Brown, M. S. Johnson, and R. L. Moore, “Alterations in peroxisome proliferator-activated receptor mRNA expression in skeletal muscle after acute and repeated bouts of exercise,” Molecular and Cellular Biochemistry, vol. 332, no. 1-2, pp. 225–231, 2009.

A. W. Thomas, N. A. Davies, H. Moir et al., “Exercise-associated generation of PPARγ ligands activates PPARγ signaling events and upregulates genes related to lipid metabolism,” Journal of Applied Physiology, vol. 112, no. 5, pp. 806–815, 2012.

F. W. Booth, M. V. Chakravarthy, and E. E. Spangenburg, “Exercise and gene expression: physiological regulation of the human genome through physical activity,” The Journal of Physiology, vol. 543, no. 2, pp. 399–411, 2002.

K. Samaras, P. J. Kelly, M. N. Chiano, T. D. Spector, and L. V. Campbell, “Genetic and environmental influences on total-body and central abdominal fat: the effect of physical activity in female twins,” Annals of Internal Medicine, vol. 130, no. 11, pp. 873–882, 1999.

L. Mustelin, K. H. Pietiläinen, A. Rissanen et al., “Acquired obesity and poor physical fitness impair expression of genes of mitochondrial oxidative phosphorylation in monozygotic twins discordant for obesity,” The American Journal of Physiology—Endocrinology and Metabolism, vol. 295, no. 1, pp. E148–E154, 2008.

D. Boomsma, A. Busjahn, and L. Peltonen, “Classical twin studies and beyond,” Nature Reviews Genetics, vol. 3, no. 11, pp. 872–882, 2002.

M. R. Queiroga, R. A. Barbieri, S. A. Ferreira et al., “Glucose metabolism in discordant monozygotic twins for cardiorespiratory fitness,” Revista Paulista de Pediatria, vol. 31, no. 1, pp. 77–82, 2013.

F. B. Hu, J. Goldberg, D. Hedeker, and W. G. Henderson, “Modelling ordinal responses from co-twin control studies,” Statistics in Medicine, vol. 17, no. 9, pp. 957–970, 1998.

T. Lohman, A. Roche, and R. Martorell, Eds., Anthropometric Standardization Reference Manual, Human Kinetics, Champaign, Ill, USA, 1988.

P. Setty, B. V. Padmanabha, and B. R. Doddamani, “Correlation between obesity and cardio respiratory fitness,” International Journal of Medical Science and Public Health, vol. 2, no. 3, pp. 300–304, 2013.

A. J. Vidal-Puig, R. V. Considine, M. Jimenez-Liñan et al., “Peroxisome proliferator-activated receptor gene expression in human tissues: effects of obesity, weight loss, and regulation by insulin and glucocorticoids,” The Journal of Clinical Investigation, vol. 99, no. 10, pp. 2416–2422, 1997.

L. R. Butcher, A. Thomas, K. Backx, A. Roberts, R. Webb, and K. Morris, “Low-intensity exercise exerts beneficial effects on plasma lipids via PPARgamma,” Medicine and Science in Sports and Exercise, vol. 40, no. 7, pp. 1263–1270, 2008.

L. Qi, F. B. Hu, and G. Hu, “Genes, environment, and interactions in prevention of type 2 diabetes: a focus on physical activity and lifestyle changes,” Current Molecular Medicine, vol. 8, no. 6, pp. 519–532, 2008.

L. Aires, M. Pratt, F. Lobelo, R. M. Santos, M. P. Santos, and J. Mota, “Associations of cardiorespiratory fitness in children and adolescents with physical activity, active commuting to school, and screen time,” Journal of Physical Activity and Health, vol. 8, pp. S198–205, 2011.

P.-L. Hsieh, M.-L. Chen, C.-M. Huang, W.-C. Chen, C.-H. Li, and L.-C. Chang, “Physical activity, body mass index, and cardiorespiratory fitness among school children in Taiwan: a cross-sectional study,” International Journal of Environmental Research and Public Health, vol. 11, no. 7, pp. 7275–7285, 2014.

H. H. M. Maes, G. P. Beunen, R. F. Vlietinck et al., “Inheritance of physical fitness in 10-yr-old twins and their parents,” Medicine & Science in Sports & Exercise, vol. 28, no. 12, pp. 1479–1491, 1996.

C. Bouchard, P. An, T. Rice et al., “Familial aggregation of VO(2max) response to exercise training: results from the HERITAGE family study,” Journal of Applied Physiology, vol. 87, no. 3, pp. 1003–1008, 1999.

A. Vaag, J. E. Henriksen, S. Madsbad, N. Holm, and H. Beck-Nielsen, “Insulin secretion, insulin action, and hepatic glucose production in identical twins discordant for non-insulin-dependent diabetes mellitus,” The Journal of Clinical Investigation, vol. 95, no. 2, pp. 690–698, 1995.

S.-W. Guo, “Does higher concordance in monozygotic twins than in dizygotic twins suggest a genetic component?” Human Heredity, vol. 51, no. 3, pp. 121–132, 2001.

C. H. Hawkes, “Twin studies in medicine—what do they tell us?” QJM: Monthly Journal of the Association of Physicians, vol. 90, no. 5, pp. 311–321, 1997.

T. W. Rowland, Pediatric Laboratory Exercise Testing: Clinical Guidelines, Human Kinetics, Champaign, Ill, USA, 1993.

T. Leskinen, K. Waller, S. Mutikainen et al., “Effects of 32-year leisure time physical activity discordance in twin pairs on health (TWINACTIVE study): aims, design and results for physical fitness,” Twin Research and Human Genetics, vol. 12, no. 1, pp. 108–117, 2009.

J. C. Hannukainen, P. Nuutila, B. Ronald et al., “Increased physical activity decreases hepatic free fatty acid uptake: a study in human monozygotic twins,” The Journal of Physiology, vol. 578, no. 1, pp. 347–358, 2007.