Journal of Applied Physiology Ad Instruments
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 82: 1802-1810, 1997;
8750-7587/97 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Loos, R.
Right arrow Articles by Vlietinck, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Loos, R.
Right arrow Articles by Vlietinck, R.

Journal of Applied Physiology
Vol. 82, No. 6, pp. 1802-1810, June 1997
EXERCISE AND MUSCLE

Gender-specific regional changes in genetic structure of muscularity in early adolescence

R. Loos1, M. Thomis2, H. H. Maes3, G. Beunen2, A. L. Claessens2, C. Derom1, E. Legius1, R. Derom1, and R. Vlietinck1

1 Center for Human Genetics and 2 Center for Physical Development Research, Katholieke Universiteit Leuven, B-3000 Leuven, Belgium; and 3 Department of Human Genetics, Medical College of Virginia, Richmond, Virginia 23298

ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
ACKNOWLEDGEMENTS
FOOTNOTES
REFERENCES


ABSTRACT

Loos, R., M. Thomis, H. H. Maes, G. Beunen, A. L. Claessens, C. Derom, E. Legius, R. Derom, and R. Vlietinck. Gender-specific regional changes in genetic structure of muscularity in early adolescence. J. Appl. Physiol. 82(6): 1802-1810, 1997.---Genetic and environmental influences on muscle circumference measurements of the extremities were estimated in 105 pairs of twins between 10 and 14 yr of age. Four circumferences, extended upper arm (EAC), forearm (FC), thigh (TC), and calf (CC), were measured. Univariate model fitting revealed that the largest part (87-95%) of the variance for all circumferences at most ages was explained by additive genetic factors. Sex differences were observed for some age categories. Multivariate analyses showed a different pattern evolving according to age and gender. In boys from 10 to 12 yr of age, one general genetic factor influenced all four circumferences. With increasing age, an arm-leg model emerged, one genetic factor influencing the arm and another genetic factor the leg circumferences. In young girls one genetic factor loaded on the proximal (EAC,TC) and another on the distal (FC,CC) circumferences. With subjects at age 14 yr, an arm-leg model was observed. High genetic correlations indicated that genetic factors related to EAC, FC, TC, and CC did not act independently. The age- and gender-specific changes in the genetic structure suggest pubertal influences. This study shows that muscle circumferences are highly heritable characteristics and are therefore a promising starting point at which to locate their genes. Gene mapping could validate the gender-specific change of the genetic structure with age and region.

skeletal muscle circumferences; twins; genetic model fitting; univariate and multivariate genetic analysis; gender and age effects


INTRODUCTION

SKELETAL MUSCLE SIZE is an indicator for undernutrition, aging, and muscular and neuromuscular diseases and is important in a variety of sports disciplines (9, 10, 19). Circumferences of extended upper arm (EAC) and calf (CC) are the most commonly used sites to estimate muscularity, followed by the circumferences of the forearm (FC) and thigh (TC). Most often, they are investigated by means of anthropometric or radiographic measurement. Muscle is the major tissue composing the circumferences.

Genetic influence on extremity circumferences is suggested by family studies, most often by comparing the correlations between parents and offspring and between full and half-siblings, controlling for age, gender, and rearing conditions (1, 13-15, 20, 27). It needs to be emphasized that the comparison of results of these studies is complicated by differences in sample characteristics (type of relationship, age, gender, size) and differences in the methods used to estimate heritability.

Little et al. (15) found evidence for different genetic influence depending on body region. However most studies (13, 14, 20, 27) found only small differences between the correlations of upper and lower extremity circumferences. Gender differences were observed by Kaur and Singh (14), Mueller and Malina (20), and Little et al. (15).

Comparison of monozygotic (MZ) and dizygotic (DZ) twins also provides insight into the possible genetic control of the trait. Twin studies indicate strong evidence for a genetic component of circumference measurements, with a heritability generally ranging between 0.53 and 0.75 (2, 4, 5, 8). Heritability coefficients are higher for CC than for upper arm circumference in adults (4), whereas only little difference was found in schoolchildren (2). Gender differences were found by some authors (5) but not by others (8).

All these studies investigated the genetic influence of each circumference separately. Until now, apparently no study ever examined the covariation between circumferences in a multivariate way nor assessed whether the genetic structure was similar in either gender or changed over time.

The purpose of this study was to estimate the genetic and environmental contributions to the variation of skeletal muscle characteristics, as assessed by four anthropometric circumferences of the extremities. The specific questions to be answered were 1) how large are the genetic and environmental influences on each of the four circumferences?; 2) are the same genes influencing the circumference in different regions?; 3) how stable is this pattern over time?; and 4) is the heritability the same in men and women and are the same genes expressed in men and women?

These questions were answered by univariate and multivariate model fitting.


METHODS

Subjects. One hundred and five twin pairs volunteered to participate in the Leuven Longitudinal Twin Study (16). All twins were reared together and were selected from the East Flanders Prospective Twin Study (6). The parents gave informed consent, and the project was approved by the Local Committee of Medical Ethics and the Committee of the National Scientific Fund. From 1985 to 1991, 10-yr-old twin pairs entered the study, and anthropometric data were collected longitudinally for 6 yr. At present, data are nearly complete until the age of 14 yr. Table 1 gives the sample size for the five yearly data points by age, gender, and zygosity. One pair skipped age 11 yr, one dropped out at age 11 yr, another dropped out at age 12 yr, and five pairs had not yet reached age 14 yr.

Table  1.   Number of twin pairs by zygostiy, gender, and age
Age, yr MZ
DZ
Total
Male Female Male Female Opposite sex

10 21 22 21 20 21 105
11 21 22 20 20 20 103
12 20 22 21 20 20 103
13 20 22 21 20 20 103
14 19 21 20 18 20 98

MZ, monozygotic; DZ, dizygotic.

Zygosity determination. In each twin gestation, the fetal membranes were examined and placental morphometry was performed. Zygosity was determined through sequential analysis (30). Monochorionic twins were classified as MZ. Genetic markers were determined in all same-gendered dichorionic twins. Alkaline placental phosphates were assayed by electrophoresis, and umbilical cord blood was used to determine the ABO, Rhesus, MNSs, Duffy, and Kell blood groups by routine methods. Five unlinked DNA restriction fragment-length polymorphisms were also studied. Same-gendered dichorionic twins with at least one different genetic marker were also classified as DZ. The probability of monozygosity on the basis of the genetic markers was calculated for all same-gendered dichorionic twins with identical genetic markers. All dichorionic twins of the same gender and with the same markers, reaching a probability of monozygosity of 0.90 or more, were considered MZ.

Variables. The measures of interest were EAC, FC, CC, and TC. All the circumferences were measured with a steel tape accurate to 1 mm and taken at the left side of the body by one of the two trained anthropometrists, whose interobserver reliability ranged between 0.94 and 0.98. The intraobserver reliability ranged between 0.97 and 0.99. A detailed description of the measurement is given by Claessens et al. (3) and of the reliability by Thomis et al. (29).

Statistical analysis. The normality of the distributions was tested by the Kolmogorov-Smirnoff two-sample test. Means were compared by t-tests and their variances by analysis of variance. The analyses and transformations were done with the SAS 6.04 computer package (SAS Institute, 1988).

Genetic analysis. The classic twin study, in which MZ and DZ twins are reared together in the same family, is one of the most powerful designs for estimating genetic and environmental effects. This approach is based on the fact that MZ twins are genetically identical and thus the correlation between genotypes (additive and dominance variance) is one, whereas DZ twins only share one-half of their genes on average, giving a correlation of 0.50 for additive variance and 0.25 for dominance variance. By definition, the common environment is shared by both individuals of a pair. It is assumed to be equal for both DZ and MZ twins reared together. Random environment is unique for each individual and also includes measurement errors.

The contribution of genetic and environmental influences to the individual differences in the circumference measurements was estimated by the method of genetic model fitting (22). With this method, the variation in the observed phenotype is decomposed into genetic and environmental variance. The genetic variance may be due to additive (A) or dominant (D) genetic influences. The environmental variance can be decomposed into common environmental factors (C), which are shared by both twins reared in the same family (e.g., diet, family attitude toward sports) and to the specific environmental factors (E), unique for each individual (e.g., subject-specific sports participation). Their influence on the phenotype is given by parameters a, d, c, and e, which are equivalent to the standardized regression coefficients of the phenotype on A, D, C, and E, respectively. The amount of variance is expressed as a proportion of the total phenotypic variance, explained by the main sources: additive (a2) or dominant (d2) genes and common (c2) and unique (e2) environment.

Hypotheses, based on a priori knowledge, were translated into mathematical models. Parameter values, predicted by the models, were optimized by maximum likelihood method. These values were compared with the observed data. When differences between predicted and observed data are small, it does not necessarily indicate that the hypothesis is correct but that it is accepted until an alternative model better fits the data.

In univariate analyses, the genetic and environmental influences responsible for variation in each of the four circumferences (EAC, FC, TC, CC) taken separately were identified. In the multivariate approach, the underlying genetic structure of the covariation between these four circumferences was determined, i.e., it was tested whether there is only one general genetic and/or environmental factor that contributes to the variation and covariation in the different circumference measures or whether two common genetic and/or environmental factors are needed. In the latter case two hypotheses were distinguished: 1) one common factor loads on proximal limbs (EAC and TC) and the other common factor on distal limbs (FC and CC), and 2) one common factor loads on the upper extremity (EAC and FC) and the other common factor on the lower extremity (TC and CC). All hypotheses tested whether additional genetic and/or environmental factors with a localized action were needed to explain the residual variance.

Univariate analysis. For each circumference measurement, correlations were computed for MZ male and female (MZM, MZF) and dizygotic male, female, and opposite-gender twin pairs (DZM, DZF, DZO). The comparison of the correlations of these five groups gave a first idea of the importance of the genetic and environmental influences. They also allowed us to compare the results of the present study with those reported in the literature.

Univariate genetic models were fitted separately to each of the four circumferences (EAC, FC, CC, TC). The goodness of fit of different models, each representing an alternative hypothesis, was tested. These models included A and/or D genetic factors and C and/or E environmental factors. The AE model assumes that variance is explained by only two factors, i.e., an additive genetic and/or a unique environmental factor. If dominance or shared environment is included, the models were referred to as ADE or ACE, respectively. It should be noted that the power of the sample is sufficient to detect additive genetic influences but may be too low to detect dominance or shared environmental effects, unless their effect accounts for >50% of the total variance. The CE and E models assume no genetic influence at all.

Gender differences can result from differences in magnitude of the genetic and/or environmental influences. This hypothesis was tested by comparing a model that equaled the genetic and environmental influences for boys and girls with a model that allowed for different estimates in either gender.

Neale and Cardon (22) provide a detailed description of univariate model fitting.

Multivariate analysis. It is known that most of the somatic characteristics are phenotypically correlated to some degree. Multivariate analysis estimates the genetic and environmental contributions to the covariance between the circumference measurements.

Multivariate models, called general models, tested whether only one genetic and only one environmental factor were responsible for the covariation between the four circumference measurements (EAC, FC, TC, CC) or whether two genetic and environmental factors better explained the observed covariation. To test the latter hypothesis, two types of models were tested: the proximal-distal models assume that one genetic and/or environmental factor explains covariation of the circumferences taken at the proximal regions of the limbs (EAC, TC) and one explains the covariation of the circumferences taken at the distal region of the limbs (FC, CC); the arm-leg models assume one factor is responsible for covariation of the circumferences of the upper extremity (EAC, FC) and one for covariation of the circumferences of the lower extremity (TC, CC). In summary, seven multivariate models were fitted to the data: the first model assumed the presence of only one general genetic (G) and one E factor to explain the covariation; the other six models included, alternatively, one G factor, one E factor, two common genetic factors [1 factor loading on proximal and 1 on distal circumferences or 1 on arm and 1 on leg circumferences (GpGd or GaGl)], and/or two common E factors [1 factor loading on proximal and 1 on distal circumferences or 1 on arm and 1 on leg circumferences (EpEd or EaEl)]. Additional variable-specific genetic (Gs) and environmental (Es) latent factors were included in all models to account for residual site-specific contributions of variation. In models with at least two common factors, factors were allowed to correlate. Alternative models, in which this correlation was constrained to zero, were also tested. Models for multivariate analysis were applied separately for men and women.

Assessment of models. Two criteria were used to evaluate models derived from different hypotheses: the goodness of fit and the parsimony.

For the univariate models, the overall goodness of fit of a model was assessed by the chi 2 goodness-of-fit statistic: the lower the chi 2, the better the fit. The significance of each latent factor (A, D, C, and E) was tested by comparing the full model to the model leaving out this factor. A significantly worse fit of this submodel compared with the full model indicated the significance of this factor. The second criterion was the parsimony, as measured by Akaike's Information Criterion (AIC) (22) of the models, which combines the chi 2 with the degrees of freedom (df) (22)
AIC = &khgr;<SUP>2</SUP> − 2 df
Models that fit data exactly have chi 2 values equal to zero, but the most parsimonious model (with the most df) will have the lowest AIC. Alternative univariate and multivariate models were compared by the AIC (22).

Because chi 2 does not provide an adequate assessment of the goodness of fit for multivariate models, the Tucker-Lewis Index (TLI) (22) was calculated, which takes parsimony into account. For this purpose, the traditional "null" model was fitted; this assumes no covariance at all but estimates only the variances of each variable of each individual (TLInull). Because of the special character of twin data, an additional null model, the diagonal model (TLIdiag), was fitted. This model estimates the within-variable across-twin covariances, in addi- tion to the variances of the variables. Parsimony is accounted for by dividing the chi 2 values of the tested model and the null model by their df: the higher the index the better the fit. The goodness of fit was further evaluated by fitting a saturated model, i.e., a double Cholesky model (22) (triangular decomposition for A and for E factors). Because this model is fully saturated, it provides an indication for the lowest limit of chi 2 obtainable with the data. Because multivariate models can be complex, alternative models were also compared by the AIC.

The statistical modeling package Mx (21) was used to evaluate the fitting of univariate and multivariate models. Results with a probability level <0.05 were considered as statistically significant.


RESULTS

Means and variances. All variables showed positive skewness. After a logarithmic transformation, all had a Gaussian distribution. Means and variances were not significantly different between boys and girls (Table 2), between MZ and DZ, or by birth rank. Compared with a singleton control group we examined (23, 26), subjects were within the expected range for their respective ages.

Table  2.   Descriptive statistics of circumference measurements by age and gender
Boys
Girls
Age, yr
Age, yr
10 11 12 13 14 10 11 12 13 14

EAC, cm
  Mean 18.8 19.5 20.4 21.4 22.9 19.4 20.1 21.0 22.0 22.9
  SD 1.5 1.5 1.7 1.9 2.0 1.9 2.0 2.2 2.3 2.2
FC, cm
  Mean 19.3 20.0 20.7 21.8 23.2 19.2 19.9 20.6 21.5 22.1
  SD 1.0 1.0 1.2 1.5 1.6 1.2 1.3 1.4 1.4 1.3
TC, cm
  Mean 39.7 41.8 43.7 45.6 47.9 41.6 43.8 46.1 48.7 50.8
  SD 3.0 3.2 3.7 3.8 3.8 3.8 4.2 4.6 4.8 4.2
CC, cm
  Mean 26.9 28.1 29.4 30.9 32.4 27.0 28.4 29.7 31.3 32.2
  SD 1.7 1.8 2.1 2.3 2.3 2.0 2.4 2.6 2.7 2.4

EAC, extended upper arm circumference; FC, forearm circumference; TC, thigh circumference; CC, calf circumference.

Univariate genetic analysis of each circumference. TABLE A1 shows the twin correlations for each variable and age of the five groups (MZM, MZF, DZM, DZF, DZO). MZ correlations for circumference measurements ranged, respectively, between 0.81 and 0.90 for EAC, 0.84 and 0.89 for FC, 0.86 and 0.95 for TC, and 0.72 and 0.93 for CC. The DZ correlations varied between 0.30 and 0.71 for EAC, 0.44 and 0.65 for FC, 0.37 and 0.76 for TC, and 0.33 and 0.65 for CC, respectively. The MZ correlations were quite stable across age, whereas the DZ correlation tended to vary more. Because MZ correlations were very high, a strong genetic component was expected. As DZ correlation exceeded one-half of the MZ correlation, a common environmental influence could be present as well.

The goodness-of-fit statistics of the most parsimonious models were satisfactory: the probability levels (P) of the chi 2 and the standardized parameter estimates are summarized in Table 3. The AE-model that incorporated additive genes and environment, unique to the individual members of a twin pair, provided the best explanation of the variation. Additive genes explained 87-95% of the variation of all circumferences, except for CC at age 14 yr, where a C factor explained 26%, and only 64% was left to be explained by genes. E factors accounted for 4-14%. Incorporation of differences due to gender (gAE) improved the fit in several cases, especially at younger ages. Between 10 and 12 yr of age, the variance in girls was ~1.2 times greater than in boys (gender heterogeneity).

Table  3.   Goodness-of-fit statistics and standardized parameter estimates of best fitting univariate model for circumferences by age
EAC
FC
TC
CC
Age, yr
Age, yr
Age, yr
Age, yr
10 11 12 13 14 10 11 12 13 14 10 11 12 13 14 10 11 12 13 14

Model gAE gAE gAE AE AE AE gAE AE AE gAE AE gAE gAE AE AE AE gAE gAE AE ACE
 chi 2 13.29 6.10 14.99 15.64 16.59 10.74 7.18 10.98 13.54 9.578 12.49 8.20 11.51 14.08 11.58 10.55 10.97 7.80 9.87 4.883
df 12 12 12 13 13 13 12 13 13 12 13 12 12 13 13 13 12 12 13 12
P 0.349 0.911 0.242 0.269 0.219 0.633 0.846 0.612 0.407 0.653 0.488 0.769 0.486 0.368 0.562 0.648 0.531 0.801 0.704 0.962
AIC  -10.71  -17.90  -9.01  -10.36  -9.414  -15.27  -16.82  -15.02  -12.46  -14.42  -13.51  -15.80  -12.49  -11.92  -14.42  -15.45  -13.03  -16.20  -16.13  -19.1
  a2 0.87 0.88 0.87 0.89 0.87 0.88 0.88 0.89 0.89 0.91 0.90 0.92 0.92 0.93 0.90 0.91 0.92 0.92 0.91 0.64
  c2 0.26
  e2 0.13 0.12 0.13 0.11 0.14 0.12 0.12 0.11 0.11 0.09 0.10 0.08 0.08 0.07 0.10 0.09 0.08 0.08 0.09 0.10
  gh 1.17 1.24 1.16 1.18 0.84 1.22 1.15 1.25 1.16

AE model includes an additive genetic factor and a specific environmental factor; gAE model assumes gender heterogeneity for the AE model; ACE model includes besides an additive genetic factor and a specific environmental factor, a common environmental factor; df, degrees of freedom; AIC, Aikake's Information Criterion; a2, c2, and: e2 proportion of total phenotypic variance, explained by main sources: additive genes and common and unique environmental factors, respectively; gh, gender heterogeneity (>1, genetic and environmental effect is more pronounced in girls than in boys; <1, genetic and environmental effect is less pronounced in girls than in boys).

Multivariate genetic analysis of covariation between circumferences taken at different regions of limbs. Three models emerged: general, proximal-distal, and arm-leg (Fig. 1). For each age, one model was selected on the basis of its goodness of fit and its parsimony (Table 4). All the selected common factor models showed an acceptable multivariate fit, as indicated by their TLInull and the TLIdiag, ranging from 0.90 to 0.98 for TLInull and from 0.88 to 0.94 for TLIdiag. At each age level and for both genders, Gs and Es factors were present. In boys between 10 and 12 yr old, the general model fitted best; i.e., one common genetic factor loaded on the four variables. At these ages, however, the arm-leg model with two common genetic factors (GaGl), one loading on the arm circumferences and the other on the leg circumferences, did not differ significantly (P > 0.05) from the selected, more parsimonious, general model, according to the chi 2 difference test. In girls 10-13 yr of age, two correlating common genetic factors were required, one for the proximal (EAC and TC) and one for the distal (FC and CC) limb circumferences. The chi 2 of the proximal-distal model was significantly (P < 0.05) lower than the chi 2 of the general model with only one common genetic factor. The correlation between the two common genetic factors, however, was high (0.95-0.97). At an older age, the best fitting model was the same in both genders. The arm-leg model contained two genetic factors, one common to the arm circumferences (EAC and FC) and one to the leg circumferences (TC and CC). Although the correlation between those two common genetic factors was high (0.93 and 0.92 in boys, 0.94 in girls), the fit was significantly (P < 0.05) better than for the general model with only one common genetic factor. At all age levels and for both genders, there was only one environmental factor common to all four variables.
Fig. 1. Best fitting multivariate models for upper arm circumference (EAC), forearm circumference (FC), thigh circumference (TC), and calf circumference (CC) by age. A: general model; this model assumes 1 common genetic factor (G), 1 common environmental factor (E), and for each variable a specific genetic (Gs) and specific environmental factor (Es). B: proximal-distal model; this model differs from the general model in that 2 common genetic factors are present, 1 loading on proximal circumferences (Gp) and 1 loading on distal circumferences (Gd). C: arm-leg model; this model differs from the general model in that 2 common genetic factors are present, 1 loading on arm circumferences (Ga) and 1 loading on leg circumferences (Gl). Numbers indicate estimated parameters.
[View Larger Version of this Image (28K GIF file)]

Table  4.   Explanatory factors of covariation among EAC, FC, TC, and CC and goodness-of-fit statistics of corresponding multivariate models by age and gender
Boys
Girls
Age, yr
Age, yr
10 11 12 13 14 10 11 12 13 14

Genetic factors
  Number of specific 4 4 4 4 4 4 4 4 4 4
  Type of common G G G GaGl GaGl GpGd GpGd GpGd GpGd GaGl
  r Between common 0.93 0.92 0.97 0.95 0.95 0.95 0.94
Environmental factors
  Number of specific 4 4 4 4 4 4 4 4 4 4
  Type of common E E E E E E E E E E
Model General General General Arm-leg Arm-leg Prox-dist Prox-dist Prox-dist Prox-dist Arm-leg
 chi 2 74.80 86.27 65.11 98.69 87.39 102.49 85.77 102.17 91.99 76.76
P 0.047 0.006 0.189 0.000 0.004 0.000 0.005 0.000 0.001 0.028
df 56 56 56 55 55 55 55 55 55 55
AIC  -37.20  -25.73  -46.89  -11.31  -22.61  -7.52  -24.23  -7.83  -18.02  -33.24
TLInull 0.95 0.93 0.98 0.90 0.92 0.92 0.95 0.92 0.93 0.95
TLIdiag 0.94 0.91 0.94 0.88 0.90 0.90 0.94 0.91 0.92 0.94

G, 1 general genetic factor loading on all 4 circumferences; GaGl, 2 genetic factors, 1 loading on arm circumferences and 1 on leg circumferences; GpGd, 2 genetic factors, 1 loading on proximal circumferences and 1 on distal circumferences; general model, 1 genetic and 1 environmental factor; arm-leg model, 1 factor explains covariation at upper extremity and 1 covariation at lower extremity; prox-dist, 1 factor explains covariation at proximal limb region and 1 at distal limb region; TLInull and TLIdiag, traditional null model and diagonal model fitted to Tucker-Lewis Index, respectively.

TABLE A2 shows the standardized estimates for the best-fitting multivariate model by age and gender. The contribution of the common genetic factor(s) (parameters 1-4) was systematically lower in boys (0.60-0.79) than in girls (0.72-0.93). Inversely, the contribution of the specific genetic factors (parameters 5-8) was higher in boys (0.09-0.29) than in girls (null to 0.07). Summing the common genetic with the specific genetic factors, the overall genetic determination varied between 0.82 and 0.91 in boys and between 0.76 and 0.91 in girls. Total environment (summing specific and common factors) was responsible for 9-18% of the variance in boys and for 4-24% in girls. Similar to the genetic contribution, the environmental common factor (9-12) accounted for the largest part of the total environmental contribution, ranging for boys between 0.01 and 0.17 and for girls between 0.02 and 0.21.


DISCUSSION

The main findings of this longitudinal study are the high genetic determination of the circumferences and the gender-specific change of the multivariate genetic structure with age.

The striking similarity of all circumferences in MZ twins, indicated by the correlations ranging between 0.72 and 0.95, approachs the test-retest reliability of 0.97-0.99 we obtained. MZ twins resemble each other nearly as much as one individual measured on different occasions. This suggested a very strong genetic determination. The correlations in DZ twins, ranging from 0.30 to 0.76, showed also a familial relatedness. Because some of these DZ correlations exceeded one-half of the corresponding MZ correlations, environmental influences, common to both twin members, could be suspected. The presence of a strong genetic influence (87-95%), suggested by comparison of the MZ and DZ correlations, was confirmed by the model fitting. This genetic influence remained stable over age in both genders. A suspected common environmental factor, however, could not be demonstrated, except for CC at the age of 14 yr. In structural equation modeling, common environmental factors can only be identified when the effect or the sample size is relatively large.

The comparison of heritability estimates from previous studies with the present ones is difficult because of differences in sample composition and in the method of heritability calculation.

The high genetic contribution confirms the results found in most earlier twin studies (2, 4, 5). Heritabilities obtained from the family studies are known to be often lower than those obtained from twin studies. Because DZ twins can be considered to be similar to other first-degree relatives, the DZ correlations were used to compare our results with the sibling-sibling and parent-offspring correlations of previous family studies.

The parent-offspring correlations (13, 14, 27) of circumference measurements ranged between 0.17 and 0.42, which are lower than our DZ correlations. This might be due to a difference in generational and in environmental conditions between parents and offspring.

The DZ correlations in this study, ranging between 0.30 and 0.76, are similar but somewhat higher than the sibling-sibling correlations of schoolchildren (15, 20) and young adults (13, 14), ranging between 0.35 and 0.45.

The narrow heritability for the EAC (0.30) in 10-yr-old siblings estimated by Bouchard et al. (1), on the basis of path analysis, is much lower than the genetic contribution (0.87) in this study and in other twin studies (2, 4). Susanne (27) found also a low heritability coefficient (estimated with the method of Fisher) of only 0.47 for arm circumferences and CC. The genetic component was of the same magnitude for the circumferences of the upper and lower extremities. This confirmed the findings of most investigators (2, 13, 14, 20, 27). Only Clark (4) and Little et al. (15) found a higher genetic component for the lower extremity circumferences compared with the upper extremity circumferences. They did not mention, however, whether this difference was statistically significant. Gender differences were mainly apparent at younger ages. Additive genes contributed equally to the variance in both boys and girls, but because girls had a larger total variance, the genetic variance was higher in girls than in boys. This may be due to the difference in individual and gender-related timing of puberty. While the girls were already in puberty, most boys were still in preadolescence. This might result in a larger variance in girls than in boys. The gender difference may also be caused by the marked increase of subcutaneous fat in the extremities in girls (28). This influences the utility of limb circumference as a measure of muscularity for girls. A higher genetic influence for women on CC and EAC is also reported in adolescents (5) and in schoolchildren (15, 20). Little et al. (15) ascribe the gender difference to the environmental impact, which might be larger in boys than in girls. Because no common environmental factor was identified in this study, this hypothesized environmental impact could not be confirmed. Hewitt (11) and Hoshi et al. (12) did not observe any gender difference in calf muscularity by radiographic analysis, either.

Because of the method used, none of the earlier studies could identify a shared environmental component, although many of them (1, 8, 13-15, 20, 25, 27) assume that environment contributes to the variance of these variables. In the present study, the presence of a shared environmental factor was only apparent at age 14 yr.

To our knowledge, no comparable multivariate path analysis of the covariation between limb circumferences has been reported. This multivariate analysis revealed the importance of genetic factors, common to different regions. With increasing age, there is evidence that the genetic architecture that underlies the variation and covariation in limb circumferences changes. This is not identical in both genders.

In boys between 10 and 12 yr of age, covariation and variation of the variables were explained by the general model, assuming one common genetic factor influencing the circumferences taken at the four regions. The genetic architecture changes with age. This change could already be suspected because the arm-leg model, assuming two correlated common genetic factors, one for the arm circumferences and one for the leg circumferences, fitted the data equally well but was less parsimonious. At ages 13 and 14 yr, this tendency was confirmed because the fit of the arm-leg model was statistically significantly better. This might suggest that a different set of genes influences arm vs. leg circumferences. The strong correlation between the two common genetic factors indicated, however, that the same genes had a different influence on arms than on legs. The fact that the presence of two common genetic factors was more pronounced at ages 13 and 14 yr than at earlier ages, possibly reflected an age-associated effect of genes at different maturational stages, as suggested by Mueller and Malina (20) and Little et al. (15). Phenotypic observation might indicate this genetic structure because boys gain relatively more muscle mass in the arms compared with the legs (28).

In girls of all ages, two common genetic factors were observed, although the loading of these factors at ages 10-13 yr differed significantly from the loading at age 14 yr. At ages 10-13 yr, the proximal-distal model fitted data best; one factor loaded on the proximal limb circumferences and another on the distal limb circumferences, whereas at age 14 yr, as in boys at this age, there was evidence for the arm-leg model; one factor loaded on the two arm circumferences and one on the two leg circumferences. This switch in girls was rather unexpected because the proximal-distal model in the younger subjects did not suggest a strong covariation between arm and leg circumferences. The accumulation of fat in the extremities may influence the circumference as an indicator for muscularity (28), which complicates the interpretation. Like in boys, the two factors correlated highly at all age levels. This again suggested that nearly the same set of genes was responsible for variation and covariation of the four variables but that they acted differently depending on the body region.

These age-associated changes and gender differences in genetic control could be influenced by several factors. First, the validity of circumference measurements as indicators of muscle tissue may result in gender differences; fat accumulation may be a confounding factor in girls. However, the pattern of age- and gender-associated variation in anthropometric estimates of limb musculature is similar to that for radiographic measurements (17). Second, genders may differ because of gender-specific changes in muscle mass and fat accumulation during puberty. Boys experience a large increase in muscle mass, more so in the upper arm than in the calf, resulting in a maximum increase in muscle mass ~3-6 mo after peak height velocity (PHV). Girls accumulate relatively more fat on the extremities and boys more on the trunk, especially during puberal years, when the typical gender-specific subcutaneous fat is accentuated (17, 18, 28). Furthermore, the difference in biological maturation of boys and girls at the same ages may resort to age-associated and gender-specific differences. On the average, girls experience PHV at ~12 yr of age and boys at ~14 yr of age (18). Consequently, girls are biologically two years more advanced than boys at the same age.

Given the timing of the changes in genetic architecture, especially in boys, it is tempting to associate these changes with hormonal secretions and the concomitant changes in muscle tissue. To demonstrate the common genetic architecture, multivariate analyses need to be carried out in which both hormonal secretion and muscle dimensions are incorporated (18).

Because circumferences are mainly determined by muscle (17, 29), our findings could be explained by the biology of this tissue. Two questions arise: 1) can the regional difference in genetic structure be explained by genes known from the literature to influence muscle tissue in different regions?; and 2) what factors could explain the temporal changes of the genetic structure in either gender?

Heritable neuromuscular diseases in childhood (7) show that the action of genes controlling skeletal muscle development is neither stable in time nor do these genes influence different regions in the same way. Some of these genes affect predominantly the muscles in certain body regions. Autosomal dominant myotonic dystrophy and different types of hereditary motor and sensory neuropathies have a more severe effect on the distal musculature than on the proximal. Spinal muscular atrophy, the Duchenne-type of muscular dystrophies, and limb-girdle muscular dystrophies affect the proximal musculature and, only in a later stage, the distal musculature. The facioscapulohumeral muscular dystrophy affects predominantly the muscles of the face, neck, and shoulder girdle. To some extent, this is what the proximal-distal model assumes. Furthermore, several genes on the X-chromosome have an impact on neuromuscular function, i.e., genes for Duchenne muscular dystrophy, spinobulbar muscular atrophy, and myotubular myopathy. Abnormalities or polymorphisms in these genes will have a different effect in male vs. female children (7). There is thus ample evidence from genetic diseases that several genes influence muscles in different regions of the body.

In both boys and girls, additive genetic factor(s), common to several regions, explained the largest part of the phenotypic variation. The genetic influence, specific for each region, was low, especially in girls. This suggested that there was no clear evidence for "circumference"-specific genes and even strengthened the assumption that genes coding for muscle circumferences did not act completely independently. Total genetic contribution (common and specific) to the variance of each phenotype accounted for up to 94% of the total variance, which confirmed the univariate results and suggested a strong genetic basis for the phenotypic variation in different muscle areas. The contribution of unique environment was low.

In summary, this study confirmed that muscle circumferences are some of the most heritable characteristics in humans and are therefore a promising starting point to locate their genes. The environmental influence, shared by both twin members, was too small to be demonstrable. Multivariate analysis of the covariance between different muscle regions showed that additive genetic factors influenced the muscle circumferences differently according to the region in boys and girls at different ages. The high correlation between the genetic factors suggested that muscle circumferences might be mediated to some extent by the same genes that act not completely independently. Gene mapping should validate this gender-specific change of the genetic structure with age and region.


ACKNOWLEDGEMENTS

We are very greatful to E. Feys for the anthropometric measurements and to Dr. R. Lysens for medical examinations.


FOOTNOTES

   The Leuven Longitudinal Twin and Family Study was supported by The Research Fund of the Katholieke Universiteit Leuven OT/86/80, The National Bank of Belgium, and the Fund of Scientific Medical Research (FGWO 3.0098.91). M. Thomis is supported by grants of the Fund of Scientific Research-Flanders (Belgium) as Researcher and H. H. Maes by National Heart, Lung, and Blood Institute Grant HL-48148 and the Carman Trust.

Address for reprint requests: R. Loos, Center of Human Genetics, Faculty of Medicine, Katholieke Universiteit Leuven, Herestraat 49, O&N6, B-3000 Leuven, Belgium (E-mail: ruth.loos{at}med.kuleuven.ac.be).

Received 5 November 1996; accepted in final form 12 February 1997.


APPENDIX

Table  A1.   Correlation of each circumference measurement between both twin members for each zygosity, gender, and age
Age, yr
10
11
12
13
14
EAC FC TC CC EAC FC TC CC EAC FC TC CC EAC FC TC CC EAC FC TC CC

MZM 0.86 0.84 0.86 0.88 0.86 0.85 0.90 0.92 0.87 0.86 0.89 0.90 0.87 0.88 0.86 0.87 0.90 0.89 0.91 0.87
MZF 0.86 0.85 0.92 0.91 0.84 0.88 0.92 0.72 0.86 0.89 0.91 0.91 0.88 0.88 0.95 0.93 0.81 0.89 0.88 0.90
DZM 0.47 0.56 0.68 0.37 0.57 0.60 0.76 0.60 0.54 0.51 0.74 0.50 0.33 0.44 0.66 0.57 0.30 0.49 0.69 0.55
DZF 0.48 0.49 0.43 0.33 0.57 0.63 0.50 0.51 0.55 0.65 0.45 0.54 0.41 0.65 0.37 0.44 0.52 0.61 0.44 0.59
DZO 0.71 0.60 0.56 0.60 0.60 0.60 0.49 0.58 0.51 0.50 0.55 0.56 0.70 0.58 0.61 0.62 0.69 0.60 0.52 0.65

MZM, monozygotic male; MZF, monozygotic female; DZM, dizygotic male; FZF, dizygotic female; DZO, dizygotic opposite sex.

Table  A2.   Standardized parameter estimates for best fitting multivariate model by age and gender, with each number representing a parameter referred to in Fig. 1
Age, yr Model Genetic
Environmental
Common
Specific
Common
Specific
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

Boys
10 General 0.65 0.60 0.70 0.65 0.17 0.29 0.15 0.20 0.17 0.06 0.01 0.14 0.01 0.06 0.04 0.01
11 General 0.73 0.65 0.69 0.81 0.12 0.26 0.15 0.10 0.12 0.03 0.14 0.04 0.03 0.06 0.01 0.05
12 General 0.74 0.65 0.71 0.79 0.09 0.24 0.15 0.11 0.14 0.06 0.13 0.06 0.03 0.06 0.01 0.04
13 Arm-leg 0.77 0.68 0.76 0.76 0.08 0.18 0.10 0.12 0.13 0.08 0.09 0.09 0.01 0.06 0.05 0.03
14 Arm-leg 0.70 0.71 0.76 0.77 0.15 0.15 0.12 0.13 0.13 0.07 0.10 0.07 0.02 0.07 0.03 0.03
Girls
10 Prox-dist 0.84 0.86 0.88 0.88 0.03 0.05 0.02 0.04 0.10 0.02 0.05 0.06 0.03 0.02 0.05 0.06
11 Prox-dist 0.77 0.91 0.82 0.93 0.07 0.02 0.04 0.00 0.14 0.03 0.12 0.05 0.02 0.05 0.03 0.02
12 Prox-dist 0.72 0.88 0.79 0.91 0.07 0.04 0.03 0.00 0.20 0.05 0.16 0.07 0.01 0.03 0.01 0.02
13 Prox-dist 0.77 0.93 0.78 0.92 0.07 0.00 0.06 0.02 0.14 0.02 0.13 0.05 0.02 0.05 0.02 0.02
14 Arm-leg 0.72 0.86 0.83 0.79 0.04 0.05 0.07 0.02 0.21 0.04 0.09 0.13 0.03 0.05 0.01 0.06

Bold values indicate paths due to same factors.


REFERENCES

1. Bouchard, C., A. Demirjian, and R. Malina. Path analysis of family resemblance in physique. Stud. Phys. Anthropol. 6: 61-70, 1980.
2. Byard, P. J., K. Sharma, J. M. Rusell, and D. C. Rao. A family study of anthropometric traits in a Punjabi Community: II an investigation of familial transmission. Am. J. Phys. Anthropol. 64: 97-104, 1984 [Medline] .
3. Claessens, A. L., B. Vanden Eynde, R. Renson, and D. Van Gerven. The description of tests and measurements. In: Growth and Fitness of Flemish Girls: The Leuven Growth Study, edited by J. Simons, G. Beunen, R. Renson, A. Claessens, B. Vanreusel, and J. Lefevre. Champaign, IL: Human Kinetics, 1990, p. 21-39.
4. Clark, P. J. The heritability of certain anthropometric characters as ascertained from measurements of twins. Am. J. Hum. Genet. 8: 49-54, 1956. [Medline]
5. Da Rocha, F. J., F. M. Salzano, H. F. Pena, and S. M. Callegari. New studies on the heritability of anthropometric characteristics as ascertained from twins. Acta Genet. Med. Gemellol. 21: 125-134, 1972.
6. Derom, C., R. Derom, H. Vanden Berghe, and R. Vlietinck. Increased monozygotic twinning rate after ovulation induction. Lancet 1: 1236-1238, 1987 [Medline] .
7. Dubowitz, V. Muscle Disorders in Childhood (2nd ed.). London, UK: Saunders, 1995.
8. Dupae, E., E. Defrise-Gussenhoven, and C. Susanne. Genetic and environmental influences on body measurements of Belgian twins. Acta Genet. Med. Gemellol. 31: 139-144, 1982 [Medline] .
9. Going, S., D. Williams, and T. Lohman. Aging and body composition: biological changes and methodological issues. Exercise Sports Sci. Rev. 23: 411-458, 1995 [Medline] .
10. Guo, S., S. Salisbury, A. Roche, W. Cameron, Chunlea, and R. Siervogel. Cardiovascular disease risk factors and body composition: a review. Nutr. Res. Rev. 14: 1721-1777, 1994.
11. Hewitt, D. Sib resemblance in bone, muscle and fat measurements of the human calf. Hum. Genet. 22: 213-221, 1957.
12. Hoshi, H., K. Ashizawa, M. Kouchi, and C. Koyoma. On the intrapair similarity of Japanese monozygotic twins in some somatological traits. Okajimas Folia Anat. Jpn. 58: 675-686, 1982 [Medline] .
13. Kapoor, S., A. K. Kapoor, R. Bhalla, and I. P. Singh. Parent-offspring correlation for body measurements and subcutaneous fat distribution. Hum. Biol. 57: 141-150, 1985 [Medline] .
14. Kaur, D. P., and R. Singh. Parent-adult offspring correlations and heritability of body measurements in a rural Indian population. Ann. Hum. Biol. 8: 333-339, 1981 [Medline] .
15. Little, B. B., R. M. Malina, P. H. Buschang, and J. H. DeMoss. Sibling correlations for growth status in schoolchildren from a rural community in Oaxaca, Mexico. Ann. Hum. Biol. 14: 11-21, 1987 [Medline] .
16. Maes, H. Univariate and Multivariate Genetic Analysis of Physical Characteristics of Twins and Parents (PhD thesis). Leuven, Belgium: Katholieke Universiteit Leuven, 1992.
17. Malina, R. Regional body composition: age, gender, and ethnic variation. In: Human Body Composition, edited by A. F. Roche, S. B. Heymsfield, and T. G. Lohman. Champaign, IL: Human Kinetics, 1996, p. 217-255.
18. Malina, R., and C. Bouchard. Growth, Maturation, and Physical Activity. Champaign, IL: Human Kinetics, 1991.
19. Micozzi, M. Applications of anthropometry to epidemiologic studies of nutrition and cancer. Am. J. Hum. Biol. 2: 727-739, 1990.
20. Mueller, W. H., and R. M. Malina. Genetic and environmental influences on growth of Philadelphia black and white schoolchildren. Ann. Hum. Biol. 7: 441-448, 1980 [Medline] .
21. Neale, M. C. Mx: Statistical Modeling (2nd ed.). Richmond, VA: Dept. of Psychiatry, 1994.
22. Neale, M. C., and L. Cardon. Methodology for genetic studies of twins and families. Dordrecht, The Netherlands: Kluwer Academic, 1992.
23. Ostyn, M., J. Simons, G. Beunen, R. Renson, and D. Van Gerven (Editors). Somatic and Motor Development of Belgian Secondary Schoolboys: Norms and Standards. Leuven, Belgium: Katholieke Universiteit Leuven Press, 1980.
24. Pollock, M., J. Graves, and L. Garazella. The measurement of body composition. In: Physiological Assessment of Human Fitness, edited by P. Maud, and C. Foster. Champaign, IL: Human Kinetics, 1995, p. 167-204.
25. Sharma, J. C., and K. Sharma. Estimates of genetic variance for some selected morphometric characters: a twin study. Acta Genet. Med. Gemellol. 33: 509-514, 1984 [Medline] .
26. Simons, J., G. Beunen, R. Renson, A. L. Claessens, B. Vanreusel, and J. Lefevre. Growth and Fitness of Flemish Girls: The Leuven Growth Study. Champaign, IL: Human Kinetics, 1990, vol. 3 (Sport Sci. Monogr. Ser.)
27. Susanne, C. Heritability of anthropological characters. Hum. Biol. 49: 573-580, 1977 [Medline] .
28. Tanner, J. M., P. C. R. Hughes, and R. H. Whitehouse. Radiographically determined widths of bone muscle and fat in the upper arm and calf from age 3-18 years. Ann. Hum. Biol. 8: 495-517, 1981 [Medline] .
29. Thomis, M., A. L. Claessens, R. Vlietinck, G. Marchal, and G. Beunen. Accuracy of anthropometric estimation of muscle cross-sectional area of the arm in males. Am. J. Hum. Biol. 9: 73-86, 1997.
30. Vlietinck, R. Determination of Zygosity of Twins (PhD thesis). Leuven, Belgium: Katholieke Universiteit Leuven, 1986. 

0161-7567/97 $5.00 Copyright © 1997 the American Physiological Society



This article has been cited by other articles:


Home page
Physiol. GenomicsHome page
G. De Mars, A. Windelinckx, W. Huygens, M. W. Peeters, G. P. Beunen, J. Aerssens, R. Vlietinck, and M. A. I. Thomis
Genome-wide linkage scan for contraction velocity characteristics of knee musculature in the Leuven Genes for Muscular Strength Study
Physiol Genomics, September 17, 2008; 35(1): 36 - 44.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
A. D. G. Baxter-Jones, J. C. Eisenmann, R. L. Mirwald, R. A. Faulkner, and D. A. Bailey
The influence of physical activity on lean mass accrual during adolescence: a longitudinal analysis
J Appl Physiol, August 1, 2008; 105(2): 734 - 741.
[Abstract] [Full Text] [PDF]


Home page
J. Med. Genet.Home page
G De Mars, A Windelinckx, W Huygens, M W Peeters, G P Beunen, J Aerssens, R Vlietinck, and M A I Thomis
Genome-wide linkage scan for maximum and length-dependent knee muscle strength in young men: significant evidence for linkage at chromosome 14q24.3
J. Med. Genet., May 1, 2008; 45(5): 275 - 283.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. J. Prior, S. M. Roth, X. Wang, C. Kammerer, I. Miljkovic-Gacic, C. H. Bunker, V. W. Wheeler, A. L. Patrick, and J. M. Zmuda
Genetic and environmental influences on skeletal muscle phenotypes as a function of age and sex in large, multigenerational families of African heritage
J Appl Physiol, October 1, 2007; 103(4): 1121 - 1127.
[Abstract] [Full Text] [PDF]


Home page
Pers Soc Psychol BullHome page
D. A. Frederick and M. G. Haselton
Why Is Muscularity Sexy? Tests of the Fitness Indicator Hypothesis
Pers Soc Psychol Bull, August 1, 2007; 33(8): 1167 - 1183.
[Abstract] [PDF]


Home page
J. Appl. Physiol.Home page
G. De Mars, A. Windelinckx, G. Beunen, C. Delecluse, J. Lefevre, and M. A. I. Thomis
Polymorphisms in the CNTF and CNTF receptor genes are associated with muscle strength in men and women
J Appl Physiol, May 1, 2007; 102(5): 1824 - 1831.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
W. Huygens, M. A. I. Thomis, M. W. Peeters, J. Aerssens, R. Vlietinck, and G. P. Beunen
Quantitative trait loci for human muscle strength: linkage analysis of myostatin pathway genes
Physiol Genomics, August 11, 2005; 22(3): 390 - 397.
[Abstract] [Full Text] [PDF]


Home page
Physiol. GenomicsHome page
W. Huygens, M. A. Thomis, M. W. Peeters, J. Aerssens, R. Janssen, R. F. Vlietinck, and G. Beunen
Linkage of myostatin pathway genes with knee strength in humans
Physiol Genomics, May 19, 2004; 17(3): 264 - 270.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Loos, R.
Right arrow Articles by Vlietinck, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Loos, R.
Right arrow Articles by Vlietinck, R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online