Journal of Applied Physiology Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 82: 959-967, 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 Thomis, M. A.
Right arrow Articles by Beunen, G. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Thomis, M. A.
Right arrow Articles by Beunen, G. P.

Journal of Applied Physiology
Vol. 82, No. 3, pp. 959-967, March 1997
EXERCISE AND MUSCLE

Multivariate genetic analysis of maximal isometric muscle force at different elbow angles

Martine A. Thomis1, Marc Van Leemputte1, Hermine H. Maes4, Cameron J. R. Blimkie5, Albrecht L. Claessens1, Guy Marchal2, Eustachius Willems1, Robert F. Vlietinck3, and Gaston P. Beunen1

1 Center for Physical Development Research, Department of Kinesiology, Faculty of Physical Education and Physiotherapy, B-3001 Leuven, and 2 Radiology Unit and 3 Center for Human Genetics, Faculty of Medicine, Katholieke Universiteit Leuven, B-3000 Leuven, Belgium; 4 Department of Human Genetics, Medical College of Virginia, Virginia Commonwealth University, Richmond, Virginia 23298; and 5 Department of Kinesiology, McMaster University, Hamilton, Ontario, Canada L8S 4K1

ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
ACKNOWLEDGEMENTS
FOOTNOTES
REFERENCES


ABSTRACT

Thomis, Martine A., Marc Van Leemputte, Hermine H. Maes, Cameron J. R. Blimkie, Albrecht L. Claessens, Guy Marchal, Eustachius Willems, Robert F. Vlietinck, and Gaston P. Beunen. Multivariate genetic analysis of maximal isometric muscle force at different elbow angles. J. Appl. Physiol. 82(3): 959-967, 1997.---The maximal isometric moment at five different elbow joint angles was measured in 25 monozygotic and 16 dizygotic male adult twin pairs (22.4 ± 3.7 yr). Genetic model fitting was used to quantify the genetic and environmental contributions to individual differences in isometric strength. Additive genetic factors explained 66-78% of the variance in maximal torque at 170-140-110 and 80° flexion (extension = 180°). At 50° flexion, common and subject-specific environmental factors contributed equally to the variation. The contribution of unique environmental factors concurs with the level of variability in muscle activation and (dis)-comfort of torque production in the specific angle. The relative contribution of lever arm and force-length relationship in torque varies according to the angle. Because these factors might be genetic, this variability is reflected in the genetic contribution at the extreme angles of 170 and 50°. Multivariate analyses suggested a general set of genes that control muscle area and isometric strength, together with a more specific strength factor. Genetic correlations were high (0.82-0.99). Genes responsible for arm-segment lengths did not contribute to muscle area nor to isometric strength.

heredity; genetic model-fitting analysis; twins; maximal voluntary contraction; muscle cross-sectional area


INTRODUCTION

THE GENETIC and environmental determinants of muscle strength as measured by isometric arm pull or handgrip tests have been studied by several authors. Transmission coefficients for isometric strength (handgrip and arm pull), on the basis of parent-child and sib-sib correlations from family studies, vary between 0.20 and 0.60 (12, 13, 18, 20-22, 25). Heritability estimates in twin studies are higher and vary between 0.60 and 0.80 (10, 11, 23). Absence of any genetic variation in isometric quadriceps strength is reported in one study (7) on the basis of a nonsignificant intrapair variance ratio of monozygotic (MZ) vs. dizygotic (DZ) twins. In 10-yr-old twins, 73% of the variance in arm pull is determined by additive genes, whereas 27% of the variance is explained by the unique environment (11). There is no evidence for gender differences in relative genetic and environmental contributions in the 10-yr-old subjects in this study. However, including strength measurements of the parents reveals a significant evidence for different genes acting in adults and in 10-year-old subjects.

Jones and Klissouras (6) demonstrate that individual differences in maximal isometric elbow flexion force at 100° (180° = straight arm) measured on a self-constructed ergometer in nine MZ and eight DZ male twins (11-17 yr) are highly determined by genetic factors. The within-pair variance F-ratio of MZ and DZ twins is significant, and the heritability estimate, computed by using the Holzinger index, is 83%.

Handgrip and arm pull are standardized and reliable tests for measuring static or isometric strength in large population-based surveys. However, from a biomechanical point of view, they provide only a crude measurement of nondynamic force, produced by the combined action of several muscles at the wrist, elbow, and shoulder.

Muscle force production is determined by general factors such as the size, number, and type of muscle fibers, pennation of the fibers, and type of contraction and activation level, which are independent of the joint angle. From previous studies, it is known that variation in physiological muscle cross-sectional area (MCSA; measured as limb circumferences or muscle widths) is moderately to highly determined by genetic factors (2, 4, 9, 11). Differences in limb-segment lengths have their impact on the measured torque, and twin studies have shown that individual differences in arm-segment length are largely explained by genetic factors (2, 10, 11). Additionally, torque at specific joint angles (19) is influenced by the force-length relationship (8) and moment arm with respect to the joint axis (lever arm) (1, 29). Muscle force production is, therefore, possibly determined by the same genetic factors that influence MCSA, arm-segment lengths, the force-length relationship, and moment arm. The relative genetic and environmental contributions of these factors to the resulting torque may vary considerably across joint angles and remain to be determined.

The method of genetic model fitting used in this study has become a standard in the analysis of twin data and been applied, in recent years, to the study of individual differences in several somatic characteristics (10, 11). Major advantages of this method over the classical heritability estimations are that 1) assumptions are explicitly stated in the model; 2) more than two groups of twins can be studied as well as extended family data; 3) a goodness-of-fit test of the model is given; 4) quantitative genetic parameters are estimated and SE can be calculated; 5) maximum likelihood solutions are obtained; 6) the fit of several alternative models can be compared; and 7) univariate and multivariate questions can be addressed. Basically, model fitting involves solving a series of simultaneous structural equations to estimate genetic and environmental parameters in such a way that the model predictions closely match the observed twin correlations.

The primary purpose of the present study was to quantify the relative genetic and environmental contributions to the observed variation in maximal isometric elbow flexion torque and its key determinants, MCSA and limb-segment lengths. Furthermore, it addressed the question of whether the same genes and environmental factors contributed to the variability in isometric torque at the different joint angles and the MCSA. A third purpose was to determine the degree of shared genetic variation among segment lengths, muscle size, and isometric torque.


METHODS

Subjects

The sample for this study was selected from within the region of Flemish Brabant, Belgium. Male volunteer twins aged 17-30 yr were included if they had equal physical activity profiles and had neither started nor stopped performing strength training during the last year. Two of 43 twins were excluded because of mental retardation and physical disability. Two subjects of this sample had strength-training experience but had not trained in the last year before the study. There were no significant differences in frequency of intense and moderate physical activity between first- and second-born twins or zygosity groups. Intrapair correlations were slightly higher in MZ twins than in DZ twins for these physical activity scores. Twenty-five MZ and 16 DZ male twin pairs, 22.4 ± 3.7 (SD) yr of age, volunteered to participate. Determination of zygosity was assessed by examination of the following genetic markers: ABO; rhesus (D, C, Cw, c, E, e); MNSs; and Duffy (a,b). Differences in two genetic markers were used to establish dizygosity. Subjects were fully informed of the measurement protocol before giving their written consent. The project was approved by the local medical ethics committee (Katholieke Universiteit Leuven).

Measurement of Maximal Isometric Strength With Promett Evaluation System

Maximal isometric strength was assessed by the maximal measured moment on an active programmable dynamometer (Promett) (27, 28). With this system, isometric, concentric, and eccentric contractions can be performed at different speeds, and movement amplitudes can be imposed by the dynamometer. Subjects were seated in a comfortable, standardized position, strapped to the chair with the right arm resting on the measurement device. They were asked to build up their maximal isometric strength and to hold this maximum for 3 s. The highest registered moment during this contraction was selected as the maximal isometric strength measure expressed in Newton meters (N · m). Measurements were made, consecutively, in the same order at these five different angles, at 170, 140, 110, 80, and 50° (with 180° = the arm in straight position), covering nearly the whole movement range at the elbow. Test-retest correlations ranged from 0.93 in the extreme angles to 0.97 in the middle angle (110°). The observer was able to evaluate each subject's maximal effort by visualized moment, and electromyographic signals were registered at biceps brachii, brachioradialis, brachialis, and triceps brachii muscles.

Anthropometric and Midarm MCSA Measurement

Weight was recorded as a general body-size measurement, and more specific dimensions of the arm, including forearm length (FAL), upper-arm length (UAL), and midarm MCSA, were taken by computed tomography. All arm-segment lengths were measured by a trained observer, technical errors of measurement were small (0.28-0.31 cm), and reliability coefficients were high (>0.98). Computed tomography-imaging scans were done at three positions of the upper arm to measure the mean cross-sectional arm muscle area (26). Starting from the midhumerus position, three scans were done, with a 3-cm interval in the direction of the hand. A fourth scan was taken at the second position, with the relaxed arm in 150° flexion. Tissue within the limits of -50 and +200 Houndsfield units was defined as muscle tissue. Both elbow flexors and extensors were included in this muscle area. Technical error of measurement for muscle area was 0.16 cm2, with a coefficient of reliability of 0.99. The mean MCSA over the four scans was used in the further analyses.

Statistical Analysis

The distributions of the maximal isometric strength outcome at each angle were tested for Gaussian normality by using the Shapiro-Wilk test. Birth-order effects and differences in means or variances between MZ and DZ twins were tested with t-tests and F-tests, respectively. Pearson correlations between first- and second-born twins were computed for MZ and DZ twins. Results were considered statistically significant if P < 0.05.

Univariate genetic analysis. The biometric approach using path-analytic models (14, 16) was applied to determine the relative genetic and environmental contributions to the observed variation in maximal isometric strength, MCSA, and arm lengths. Applied on a sample of MZ and DZ twins, the general quantitative genetic model decomposes the total observed phenotypic variance into genetic additive (A) and dominance (D) variances, environmental variation due to environmental factors shared by twins reared in the same family (C), and to non-shared environmental factors (E). The influence of these sources A-E on the phenotypic variation is given by parameters a-e, which are equivalent to the standardized regression coefficients of the phenotype on A-E, respectively. For a model including A, C, and E factors (Fig. 1), the phenotypic score is expressed as PT1 = aA1 + cC1 + eE1. In the twin model, subscript 1 is for the first-born twin and 2 is for the second-born twin. Squaring the factor loadings yields the variance (V) explained by each component (VA = a2, VD = d2, VC = c2, VE = e2). The contribution of genes and environment to the total variance is reported in the standardized form by dividing the specific variance component by the total phenotypic variance. In this model, it is also assumed that genetic and environmental factors do not correlate or interact and that there is no significant parental correlation for these characteristics.
Fig. 1. Univariate genetic standard path model. PT1 (moment 140°) and PT2 (moment 140°), observed static moment in first-born (twin 1) and second-born (twin 2) twins. a-e, Standardized regression coefficients of phenotype on factors A-E. Additive genetic factors (A) are correlated [1 in monozygotic (MZ) twins, 0.5 in dizygotic (DZ) twins], and common environmental factors (C) correlate (1 in both MZ and DZ twins). These latent factors, together with non-shared environmental factors (E), cause variation in phenotypes of these twins.
[View Larger Version of this Image (15K GIF file)]

Data on twin 1 and twin 2 were summarized in 2 × 2 variance-covariance matrixes. The maximum likelihood estimation of parameters a-e was done in Mx (15). The goodness-of-fit between the observed measurements and the expected values on the basis of the model was assessed by chi 2. Low chi 2 indicates consistency of the model with the data, whereas high chi 2 values indicate poor fit. Seven alternative hypotheses about the genetic and environmental determinations of isometric strength were formulated. The simplest model is one in which all variation is explained by unique environmental factors that influence each individual separately, e.g., specific training status and unique lifetime events (model E). The significance of additional sources A and D (or C) was tested by likelihood-ratio tests and Akaike's Information Criterion (AIC = chi 2 - 2 · degrees of freedom) (30) by comparing the submodels to the full genetic model. The AIC combines the goodness-of-fit of a model (the discrepancy of expected to observed covariance matrixes) with its simplicity (the degrees of freedom of the model), resulting in a measure of parsimony. The reciprocal sibling interaction or phenotypic (P) interaction model (APE) tested whether the performance of one twin influenced the performance of his cotwin. An additional latent factor is only significant if the fit of the submodel without the parameter is significantly worse than the full genetic model. The maximum likelihood estimates of the parameters of the most parsimonious model were computed. These parameter estimates were expressed in percentages of the total variance explained by genetic and environmental factors.

Multivariate genetic analyses. To test whether the same genes influenced static strength at different joint angles, the contribution of genetic and environmental factors to the covariance between elbow flexion torques was estimated by using multivariate models. Mean midarm MCSA was also included in this multivariate analysis because this is a key physiological determinant and one of the most important covariates of static strength (5). Multivariate analyses of the isometric strength were performed only on the three middle angles at 140, 110, and 80°. These phenotypes were also modeled in this order because this was the test order. Strength measurements at 170 and 50° were less reliable due to the unfamiliar position in which the subjects were asked to perform a flexing moment and were, therefore, not included in the analyses.

A first model included four distinct genetic effects, A1-A4, and four distinct environmental effects, E1-E4 (Fig. 2). In a Cholesky or triangular decomposition of latent factors, all the genetic variation in MCSA was associated with factor A1, whereas paths 2-4 indicated shared genetic effects with all three strength measurements. Similarly, the genetic effects of factor A2 were associated with static strength at 140° ( path 5), and the other strength measurements ( paths 6 and 7) but not with MCSA. In other words, these are residual genetic effects after accounting for the first general genetic factor. The genetic effects of A3 were shared only by strength measurements at both 110 and 80° ( paths 8 and 9), whereas A4 represented the genetic effects that were unique for the variation in strength at 80° ( path 10). The same triangular decomposition was modeled for environmental factors that were unique to each individual but shared between characteristics. The model provided an indication for the lower limit of chi 2 obtainable with the data. This saturated model, however, did not provide a simple or biological explanation of the data; therefore, submodels were tested. These submodels had more restrictions on the number of genetic and environmental parameters. The most parsimonious model was selected on the basis of the lowest AIC and plausible biological interpretation. The standardized genetic and environmental variance and covariance components were computed as well as the standardized genetic and environmental correlations (16).
Fig. 2. Illustration of Cholesky decomposition of genetic and unique environmental effects. Genetic effects of muscle area (A1) are shared with genetic variation in static (Stat) strength at 140, 110, and 80°. A second distinct genetic factor (A2) accounts for genetic variation in static moment at 140, 110, and 80°. Genetic factor A3 presents genetic factors shared by static moments in 110 and 80° but not in 140° nor in muscle area variation. A4 is genetic factor that contributes uniquely to remaining genetic variation in static strength at 80° after genetic variation shared by other traits (in A1-A3) is partialized out. E1-E4, same structure in unique environmental determination. CSA, cross-sectional area.
[View Larger Version of this Image (32K GIF file)]

To determine the degree of shared genetic and environmental influences, a second set of multivariate analyses, including static strength at 110°, FAL, UAL, and midarm MCSA, was studied. Analyses started with the Cholesky decomposition model (order of variables: FAL-UAL-MCSA-static torque at 110°), and, subsequently, submodels with more restrictions on genetic and environmental parameters were also tested. Again, the most parsimonious model was selected on the basis of the lowest AIC and biological interpretation.


RESULTS

Descriptive Statistics

Data showed a Gaussian distribution, and there were no significant birth-order effects. The group of DZ twins was, although not statistically significant, smaller and weaker overall than the MZ group of twins (twins taken as subjects) (Table 1). The mean MCSA was significantly greater in MZ than in DZ twins; however, no significant differences were found in the means or variances of maximal isometric moments between MZ and DZ twins. Highest maximal isometric moments were obtained at 80 and 110° for both MZ and DZ groups, with decreasing maximal moments at more extreme angles.

Table 1. Descriptive statistics for maximal isometric moment at different angles: anthropometric measurements and muscle cross-sectional area for different groups


Variables Total Group (n = 82) Twins
MZ (n = 50) DZ (n = 32)

Age, yr 22.4 ± 3.7  22.28 ± 3.1  22.54 ± 4.5 
Height, cm 177.4 ± 7.7  178.2 ± 7.8  176.5 ± 7.6 
Weight, kg 71.1 ± 9.7  72.9 ± 9.4  68.1 ± 9.4 
Mean fl. + ext. arm-muscle area, cm2 49.23 ± 7.4  50.8 ± 6.7* 46.8 ± 7.8 
Forearm length, cm 25.6 ± 1.2  25.6 ± 1.2  25.6 ± 1.2 
Upper-arm length, cm 36.4 ± 1.9  36.4 ± 1.8  36.4 ± 2.1 
Maximal moment
  at 170°, N · m 33.9 ± 10.2  34.3 ± 10.5  33.3 ± 9.8 
  at 140°, N · m 41.4 ± 10.8  41.9 ± 11.1  40.5 ± 10.5 
  at 110°, N · m 51.2 ± 11.4  52.1 ± 11.7  49.6 ± 11.1 
  at 80°, N · m 53.7 ± 11.5  54.4 ± 12.4  52.8 ± 10.1 
  at 50°, N · m 41.6 ± 10.1  42.0 ± 11.2  41.0 ± 8.2

Values are means ± SD. MZ, monozygotic; DZ, dizygotic; fl, flexor; ext, extensor. * Significantly different means, MZ vs. DZ group (P < 0.05).

Correlation Analysis

Except for strength at the 50° angle, the MZ twin phenotypic correlations (Table 2) were higher than the DZ twin correlations at all angles, which indicated the presence of additive genetic variation. The difference between the theoretical degree of genetic relatedness in MZ twins (r = 1.0) and the observed MZ correlation gives an indication of the importance of unique environmental factors: these values varied between 20 and 30% for all angles, except for the 50° angle. If the DZ correlation was lower than one-half of the MZ correlation, as was the case for the mean arm muscle area and all static moments except at the 50° angle, nonadditive genetic factors can be expected to contribute to the observed variance. At 50°, the DZ twin correlation was higher than the MZ twin correlation, which suggests that the similarity found at this angle was determined by common environmental factors and unique environmental factors. DZ correlations were also higher than one-half of the MZ twin correlations for height, FAL, and UAL, which also suggested shared environmental effects.

Table 2. MZ and DZ twin correlations for anthropometric and maximal isometric strength measurements at different angles


Variable Twins
MZ (n = 25) DZ (n = 16)

Height, cm 0.96 0.64
Weight, kg 0.93 0.44
Mean fl. + ext. arm-muscle area, cm2 0.81  -0.01
Forearm length, cm 0.83 0.52
Upper-arm length, cm 0.84 0.72
Maximal moment
  at 170°, N · m 0.72 0.05
  at 140°, N · m 0.77 0.06
  at 110°, N · m 0.79 0.14
  at 80°, N · m 0.70 0.14
  at 50°, N · m 0.51 0.58

Univariate Genetic Analysis

For all angles, except for the maximal isometric strength at the smallest angle (50°), the AE model was the most parsimonious model to explain the variance in maximal isometric strength. Individual differences in the somatic dimensions were also best explained by a model including additive genetic effects and unique environmental effects. The percentages of explained variance of the factors in the most parsimonious model are presented in Fig. 3.
Fig. 3. Contribution of additive genetic, unique, and common environmental factors (in %) in the best fit model in maximal isometric moment at all angles and anthropometric measurements. M, muscle; FAL, forearm length; UAL, upper-arm length).
[View Larger Version of this Image (39K GIF file)]

Additive genetic factors (a2) explained 66-78% of the variance in torque, whereas unique environmental factors (e2) accounted for 22-34% of the variance. From the MZ and DZ correlations, nonadditive genetic factors were expected to contribute to the total variance; they were, however, not significant. The lower chi 2 value of the ADE model for angles 170-80° was not significantly lower than the chi 2 of the AE model; due to the difference in one degree of freedom, this resulted in a lower AIC for the AE model. To detect a rather small contribution of nonadditive genetic factors, as could be the case in the present study, sample sizes would need to be increased considerably (1,800 twin pairs are needed to reject an AE model with a power of 80% when an ADE model is the true model, with respective contributions of a2 = 0.50, d2 = 0.30, e2 = 0.20). The upper and lower 95% confidence interval of the heritability estimates were computed by maximum likelihood procedures (17). The lower confidence-interval limit ranged from 0.40 (80° angle) to 0.56 (110° angle) and the upper limit from 0.81 (80° angle) to 0.88 (110° angle).

Genetic factors seem to be absent in the determination of maximal isometric strength at the smallest angle (50°). The ACE and APE models, however, had the same chi 2 value as the CE model and fitted the data equally well, but the simpler CE (more degrees of freedom) model was preferred because of its parsimony. Common and unique environmental factors were of equal importance at this angle.

Both height and weight are highly genetically determined (a2 = 0.95 and 0.92, respectively). FAL and UAL also had high heritabilities (a2 = 0.84 and 0.86, respectively). The variation in mean MCSA was best explained by the APE model, with a high genetic component (a2 = 0.92) and the indication of a significant negative phenotypic interaction factor (-0.24). This small negative phenotypic interaction factor indicates that a large arm muscle area in one twin goes together with a smaller arm muscle area in the other twin.

Multivariate Genetic Analysis

MCSA-static torque at 140, 110, and 80°. A full Cholesky decomposition of both additive genetic and unique environmental factors of the covariance structure fit the data fairly well (chi 252 = 56.72, P = 0.30, AIC -47.28). Further testing of more restricted models resulted in the model presented in Fig. 4 (chi 260 = 62.81, P = 0.34, AIC = -55.19). Genetic factor A1 represented a general genetic factor shared by all characteristics, whereas A2 was a genetic factor common to all torque measurements. There was only one environmental factor common to all torque measurements (E2), and additional characteristic-specific environmental factors were also included in the model. The loadings of these variable-specific environmental influences could be equated between the strength measurements.
Fig. 4. Graphical representation of most parsimonious model explaining covariation in MCSA and static strength measurements at 3 angles. See text for explanation of symbols and paths.
[View Larger Version of this Image (25K GIF file)]

In model 1, the general genetic factor (A1) accounts for 41-83% of the variation in all characteristics (Table 3). The contribution of unique environmental factors to the variation of each trait was mostly explained by environmental factors that are unique to the specific trait and reached ~17-24% (paths 11, 15, 18, 20; Table 3). Very high genetic correlations were found among the three torque measurements (0.95-0.99), whereas there was a lower genetic relationship between MCSA and the torque measurements (0.76-0.90) (Table 3). The covariation among all traits was largely caused by genetic factors (86-97, Table 3), whereas environmental factors were only responsible for 3-14% of the covariation (Table 3).

Table 3. Model-fitting results in MCSA-torque at 140-110-80° in Cholesky model and most parsimonious model


Model 1 
Most Parsimonious Model
A1 A2 A3 A4 A1 A2

Proportion of explained variance by genetic factors
MCSA 0.83(1) 0.82(1)
140° 0.61(2) 0.15(5) 0.66(2) 0.06(5)
110° 0.57(3) 0.13(6) 0.02(8) 0.61(3) 0.09(6)
80° 0.41(4) 0.27(7) 0.03(9) 0.0(10) 0.50(4) 0.24(7)

MCSA 140° 110° 80° MCSA 140° 110° 80°

Genetic correlation for proportion of covariation explained by genetic factors
MCSA 0.90 0.89 0.76 0.95 0.93 0.82
140° 97 0.99 0.95 100 0.99 0.95
110° 97 87 0.97 100 87 0.97
80° 89 90 86 100 87 86

E1 E2 E3 E4 E2 Es1 Es2 Es3 Es4

Proportion of explained variance by environmental factors
MCSA 0.17(11) 0.18(11)
140° 0.002(12) 0.24(15) 0.10(15) 0.17(21)
110° 0.003(13) 0.05(16) 0.23(18) 0.12(16) 0.18(21)
80° 0.03(14) 0.02(17) 0.02(19) 0.22(20) 0.10(17) 0.16(21)

MCSA 140° 110° 80° MCSA 140° 110° 80°

Environmental correlation for proportion of covariation explained by environmental factors
MCSA 0.09 0.11 0.33 0.0 0.0 0.0
140° 3 0.43 0.28 0 0.38 0.37
110° 3 13 0.38 0 13 0.39
80° 11 10 14 0 12 14

Data show proportion of variance in each phenotype explained by genetic factor (A1-A4); bivariate genetic correlations and proportion of covariation explained by genetic factors; proportion of explained variance in each phenotype by unique environmental factors (E1-E4); and bivariate environmental correlations and proportion of covariation explained by environmental factors, respectively. MCSA, muscle cross-sectional area; nos. in parentheses, path coefficients from Fig. 2 (model 1) and Fig. 4 (most parsimomious model); Es1-Es4; variable-specific environmental factors.

Because genetic factors A3 and A4 did not contribute significantly to the variation in torque at 110 and 80°, they were left out of the most parsimonious model. In this second model, 50-66% of the variation in maximal isometric torque was explained by genetic factors that also explain 82% of the variation in muscle area (Table 3). Other genetic factors that only influence the three torque measurements are responsible for a smaller additional amount of variation in torque (0.06-0.24%). The environmental factor (E2) that influences all three torque measurements, like test conditions and so on, explained ~10-12% of the variance, and environmental factors that are unique to each measurement accounted for ~16-18% of variation (Table 3). The covariation of muscle area with torque measurements was completely explained by genetic factors because there was no environmental covariation present in this model (Table 3).

FAL-UAL-MCSA torque at 110°. Model 1 fitted the data only moderately well (chi 252 = 68.71, P = 0.06, AIC = -35.29). After inspection of the path coefficients of this model, a more restricted model was formulated as presented in Fig. 5. Genetic factor A1 represented genetic factors shared by both arm-length measurements, whereas A3 was a genetic factor common to both muscle area and static torque at 110°. A correlational path between both genetic factors was included. Environmental contributions to the characteristics seemed to be specific to each trait and were not shared between traits (chi 259 = 82.361, P = 0.024, AIC = -35.639).
Fig. 5. Graphical representation of most parsimonious model explaining covariance among FAL, UAL, MCSA, and isometric torque at 110° flexion. See text for explanation of symbols and paths.
[View Larger Version of this Image (25K GIF file)]

Inspection of the standardized genetic paths in model 1 indicated that genetic factors shared by all traits (A1) as well as A2 did not contribute significantly to the variation in MCSA and torque at 110° flexion (paths 3, 4, 6, 7; Table 4). The genetic correlation between arm-segment lengths (FAL/UAL) was high (0.77), and the genetic predisposition was highly correlated between arm muscle area and static torque (0.71) (Table 4). Genetic correlations among arm-segment lengths and muscle area and torque were, however, low. The environmental contributions to arm lengths, MCSA, and torque at 110° were trait specific and explained ~11-17% of the variance in each trait (Table 4).

Table 4. Model-fitting results in FAL-UAL-MCSA-torque at 110° in Cholesky model and most parsimonious model


Model 1 
Most Parsimonious Model
A1 A2 A3 A4 A1* A3 As1 A s2 As3 As4

Proportion of explained variance by genetic factors
FAL 0.84(1) 0.60(1) 0.25
UAL 0.52(2) 0.34(5) 0.83(2) 0.04
MCSA 0.01(3) 0.006(6) 0.81(8) 0.52(8) 0.30
110° 0.05(4) 0.02(7) 0.35(9) 0.35(10) 0.79(9) 0.0

FAL UAL MCSA 110° FAL UAL MCSA 110°

Genetic correlation for proportion of covariation explained by genetic factors
FAL 0.77 0.11 0.26 0.82 0.19 0.24
UAL 91 0.14 0.29 100 0.22 0.28
MCSA 97 87 0.71 100 100 0.79
110° 92 87 88 100 100 100

E1 E2 E3 E4 Es1 Es2 Es3 Es4

Proportion of explained variance by environmental factors
FAL 0.16(11) 0.15(11)
UAL 0.03(12) 0.11(15) 0.13(15)
MCSA 0.0(13) 0.003(16) 0.17(18) 0.18(18)
110° 0.002(14) 0.006(17) 0.04(19) 0.17(20) 0.21(20)

FAL UAL MCSA 110° FAL UAL MCSA 110°

Environmental correlation for proportion of covariation explained by environmental factors
FAL 0.46 0.01 0.09 0.0 0.0 0.0
UAL 9  -0.11 0.19 0 0.0 0.0
MCSA 3 13 0.40 0 0 0.0
110° 8 13 12 0 0 0

FAL, forearm length; UAL, upper-arm length; nos. in parentheses, path coefficients from from Fig. 2 (model 1) and Fig. 5 (most parsimonious model. * Correlation between genetic factors A1 and A3 is 0.28.

In the most parsimonious model, the genetic factor accounting for the genetic covariation in arm-segment lengths (A1) explained 60% of the variation in FAL and 83% of the variation in UAL (Table 4). Specific genetic factors accounted for another 25% of the variation in FAL. The shared genetic (A3) factor for MCSA and torque at 110° explained 79% of the variation in torque and 52% of the variation in MCSA (Table 4). Both shared genetic factors were correlated to a rather small extent (0.28), which meant that a small part of the genes, influencing arm-segment lengths, arm muscle size, and torque, are the same or act in the same pathways. All environmental variation was specific to each trait and accounted for 13-21% of the variation (Table 4).


DISCUSSION

Previous studies on the heritability of isometric strength differ from this study in strength evaluation, size, age, and gender of sample and use of heritability indexes. Heritability estimates in previous studies were calculated on the basis of correlations, intrapair variance ratios, and the Holzinger index, which fail to include estimates of the contribution of common environment, dominance, and so on, and may, therefore, be biased heritability estimations.

Maximal voluntary isometric torque of the elbow flexors at different angles was highly genetically determined in this group of adult men. Heritability estimates ranged from 66 to 78%, whereas unique environmental factors accounted for the rest of the variance. The confidence intervals around the heritability estimates were, however, rather large, probably due to the relatively small sample size of this study. Genetic effects seemed to be absent in the determination of maximal isometric strength at 50°, whereas common environmental effects explained one-half of the variance.

The heritability estimates in this study corresponded to those found in other twin studies measuring maximal static or dynamic strength by arm pull, handgrip, combined strength scores, or pull-ups (11, 23). These studies report heritability estimates ranging from 0.65 to 0.85. Jones and Klissouras (6) studied maximal isometric arm strength by a similar dynamometer in nine MZ and eight DZ male twins (11-17 yr). Using the Holzinger index, they estimated that genetic variation accounts for 83% of the variance in maximal arm flexor strength at 100°. This estimate is somewhat higher than the 78% genetic determination found for maximal isometric moment at 110° in this study. However, taking into consideration the confidence intervals around the heritability estimates associated with the sample size of this study, both estimates are quite similar.

Unique environmental effects, accounting for a part of the variance, can be interpreted as factors that contribute to the differences between members of a twin pair. Test-specific measurement error and motivational factors are included in this factor, although the latter could also be partly influenced by genes. Studies that demonstrate that the motivation or activation by electromyography is less stable at the 50 and 170° angles than at the middle angles is confirmed in our findings (28). Subjects also experience greater performance discomfort at 170 and 50° elbow flexion than at the middle angles of 140-80° flexion. We found that the degree of variability in muscle activation and performance discomfort was reflected by the general shape of unique environmental variation at the different angles (see Fig. 3), with the highest unique environmental impact at the extreme angles.

Although the factors mentioned above can explain some of the differences in unique environmental factors at different angles, the question remained as to why the genetic contribution was rather small or absent at the smallest angle (50°), whereas genetic factors were important (70%) at the 170° angle. In this regard, we can consider the dual origin of torque as a product of two factors: the moment lever arm and the force-length relationship (8, 19). Torque at 170° is relatively more influenced by the force-length factor, where the muscle is closer to its physiological resting length, than by the small moment arm. At 50° of flexion, the large moment arm has the greatest impact on the torque (1). Because genetic factors might be important in both the physiological resting length of the muscles and in the bone structure that largely determines the moment arm at 170° (2, 9-11), higher expression of genetic factors could be expected at this angle. The moment arm at 50°, however, is less dependent on bone structure, and individual variation in muscle origins and insertions will have smaller effects on the torque, resulting in lower expression of genetic factors.

Variation in physiological MCSA, univariately studied in families and twins by indicators of muscle mass (limb circumferences or muscle widths), is moderately to highly determined by genetic factors (2, 4, 9, 11), which was confirmed by the large genetic component (92%) for MCSA in this study. There was no clear biological explanation for the presence of a negative phenotypic interaction factor in MCSA in this study, and we mean that this factor rather reflects the significant difference in the mean and difference in variance in the MZ and DZ groups for this character.

In the multivariate analyses, the contribution of genetic and environmental factors to the observed covariation between isometric strength at different joint angles and MCSA was investigated and quantified. One common genetic factor explained the largest part of the covariation between MCSA and isometric torque at three elbow angles (genetic correlation 0.82-0.95). This suggested that, to a large extent, the same genes account for the variation in physiological cross-sectional area and isometric torque. These genetic effects could be the genes programming the total number of muscle fibers to develop and size of muscle fiber area, which, in turn, determine the individual differences in muscle mass (24), and, therefore, largely the variation in strength. A second common genetic factor explained a small part of the variation in all three torque measurements (6-24%). These genetic factors that were not shared with muscle mass could, in part, be genes coding for specific myofilament types and their proportion in a muscle or genetic factors regulating muscle innervation and coordination and energy supply and metabolism in the muscle during maximal contraction.

Unique environmental factors explained 12-14% of the observed covariation, and environmental correlations between the strength measurements ranged from 0.37 to 0.39. Unique environmental factors that affect the three strength measurements may reflect individual strength-training experience and Promett-related measurement errors.

The only study that reports genetic correlations between muscle area and strength was performed on 10-yr-old twins. A genetic correlation of 0.46 was found between an individual's factor score on a muscularity factor, on the basis of bone-width and extremity-circumference measurements and static strength measured by arm pull (10).

Our results suggested that genes played an important role in the variation of the FAL and UAL because estimated heritabilities were 84 and 86%, respectively. These heritabilities were higher than reported by Maes et. al. (10, 11), who found a significant influence of shared environmental factors in the variation of FAL in 10-yr-old twins. Because arm-segment lengths might influence absolute torque measurements at the elbow, both segment lengths and MCSA were included in a multivariate model to study the covariation with isometric torque. We found that subjects genetically predisposed to have longer arm segments did not seem to be predisposed to have a larger MCSA nor larger torque. Arm-segment lengths and muscle dimensions thus appear to be coded by a different set of genes. Subject-specific environmental factors, such as training status, had no general effect on the set of different length, muscle mass, and strength measurements but seem to be specific to arm-segment length, muscle mass, or strength.

In summary, we found that the variability in maximal isometric strength of the elbow flexors between 170 and 110° is highly determined by genetic factors. The contribution of unique environmental factors was in agreement with the level of variability in muscle activation and strength performance discomfort at specific angles. Furthermore, we found that the underlying relative contribution of lever arm and force-length relationship in the product of torque at different angles was reflected in the difference in genetic contribution at the extreme angles. With multivariate genetic model fitting, evidence was found for a general set of additive genetic factors that accounted for variation in muscle area and the covariation in isometric strength, together with a more specific strength factor. Specific environmental factors and measurement errors were shared among the strength measurements but not with muscle area. Furthermore, we found that genes coding for arm-segment lengths did not contribute to muscle area nor to isometric strength.


ACKNOWLEDGEMENTS

We thank I. Vallaey for assistance in Promett testing and training supervision, B. Staf for taking the blood samples, and D. Kellens for serological preparations. We also thank Dr. R. Andries for guidance during Promett evaluations and all the twins for their participation in the study.


FOOTNOTES

   This research was supported by grant 0T/92/27 of the Research Fund of the Katholieke Universiteit Leuven. M. Thomis is supported by the Fund For Scientific Research-Flanders (Belgium) as Researcher. H. Maes is supported by the Carman Trust and National Institutes of Health Grants HL-48148 and MH-45268.

Address for reprint requests: M. Thomis, Faculty of Physical Education and Physiotherapy, Katholieke Universiteit Leuven, Tervuursevest 101, B-3001 Leuven, Belgium.

Received 3 September 1996; accepted in final form 22 October 1996.


REFERENCES

1. An, K. N., F. C. Hui, B. F. Morrey, R. L. Linsheid, and E. Y. Chao. Muscles across the elbow joint: a biomechanical analysis. J. Biomech. 14: 659-669, 1981. [Medline]
2. 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]
3. Engström, L. M., and S. Fishbein. Physical capacity in twins. Acta Genet. Med. Gemellol. 26: 159-165, 1977. [Medline]
4. Hoshi, H., K. Ashizawa, M. Kouchi, and C. Koyama. On the intrapair similarity of Japanese monozygotic twins in some somatological traits. Okajimas Folia Anat. Jpn. 58: 675-686, 1982. [Medline]
5. Ikai, M., and T. Fukunaga. Calculation of muscle strength per unit cross-sectional area of human muscle by means of ultrasonic measurement. Int. Z. Angew. Physiol. Einschl. Arbeitsphysiol. 26: 26-32, 1968.
6. Jones, B., and V. Klissouras. Genetic variation in the strength-velocity relation of human muscle. In: Sport and Human Genetics, edited by R. M. Malina, and C. Bouchard. ChampaignIL: Human Kinetics, 1986, p. 155-163.
7. Karlsson, J., P. V. Komi, and J. H. Viitasalo. Muscle strength and muscle characteristics in monozygous and dizygous twins. Acta Physiol. Scand. 106: 318-325, 1979.
8. Kaufman, K. R., K. N. An, and E. Y. S. Chao. Incorporation of muscle architecture into the muscle length-tension relationship. J. Biomech. 22: 943-947, 1989. [Medline]
9. Livshits, G., I. Otremski, and E. Kobyliansky. Genetics of human body size and shape: complex segregation analysis. Ann. Hum. Biol. 22: 13-27, 1995. [Medline]
10. Maes, H. H. Univariate and Multivariate Genetic Analysis of Physical Characteristics of Twins and Parents (PhD dissertation). Katholieke Universiteit Leuven, Belgium, 1992.
11. Maes, H. H., G. Beunen, R. Vlietinck, J. Lefevre, C. Van den Bossche, A. Claessens, R. Derom, R. Lysens, R. Renson, J. Simons, and B. Vanden Eynde. Heritability of health- and performance-related fitness. Data from the Leuven Longitudinal Twin Study. In: Kinanthropometry, edited by W. Duquet, and J. Day. London: E & FN Spon, 1993, vol. IV, p. 140-149.
12. Malina, R. M., and W. H. Mueller. Genetic and environmental influences on the strength and motor performance of Philadelphia school children. Hum. Biol. 53: 163-79, 1981. [Medline]
13. Montoye, H. J., H. L. Metzner, and J. B. Keller. Familial aggregation of strength and heart rate response to exercise. Hum. Biol. 47: 17-36, 1975. [Medline]
14. Neale, M. C., A. C. Heath, J. K. Hewitt, L. J. Eaves, and D. W. Fulker. Fitting genetic models with Lisrel: hypothesis testing. Behav. Genet. 19: 37-49, 1989. [Medline]
15. Neale, M. C. Mx: Statistical Modeling (2nd ed.). Richmond, VA: Medical College of Virgnia, Department of Human Genetics, 1994.
16. Neale, M. C., and L. R. Cardon. Methodology for genetic studies of twins and families. Dordrecht, The Netherlands: Kluwer Academic, 1992, p. 496.
17. Neale, M. C., and M. B. Miller. The use of likelihood-based confidence intervals in genetic models. Behav. Gen. In press.
18. Pawlak, K. Heritability of strength and speed: methods of testing and evolution. In: Genetics of Psychomotor Traits in Man. Proceedings of the International Society of Sports Genetics and Somatology, edited by N. Wolanski, and A. Siniarski. Warsaw: College of Physical Education, 1984, p. 189-216.
19. Pérès, G., and B. Maton. Validity of the muscle equivalent concept in human muscular isometric contractions at different elbow angles. Electromyogr. Clin. Neurophysiol. 27: 135-143, 1987. [Medline]
20. Pérusse, L., C. Leblanc, and C. Bouchard. Inter-generation transmission of physical fitness in the Canadian population. Can. J. Sport Sci. 13: 8-14, 1988. [Medline]
21. Pérusse, L., C. Leblanc, A. Tremblay, C. Allard, G. Théiault, F. Landry, J. Talbot, and C. Bouchard. Familial aggregation in physical fitness, coronary heart disease risk factors and pulmonary function measurements. Prev. Med. 16: 607-615, 1987. [Medline]
22. Pérusse, L., G. Lortie, C. Leblanc, A. Tremblay, G. Thériault, and C. Bouchard. Genetic and environmental sources of variation in physical fitness. Ann. Hum. Biol. 14: 425-434, 1987. [Medline]
23. Reed, T., R. Fabsitz, J. Selby, and D. Carmelli. Genetic influences and hand grip strength norms in the NHLBI twin study male aged 59-69. Ann. Hum. Biol. 18: 425-432, 1991. [Medline]
24. Sale, D. G., J. D. MacDougall, S. E. Alway, and J. R. Sutton. Voluntary strength and muscle characteristics in untrained men and women and male bodybuilders. J. Appl. Physiol. 62: 1786-1793, 1987. [Abstract/Free Full Text]
25. Szopa, J. Familial studies on genetic determination of some manifestations of muscular strength in man. Genet. Pol. 23: 65-79, 1982.
26. 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.
27. Vande Broek, G., M. Van Leemputte, R. Andries, and E. J. Willems. Mechanical muscle properties after two types of plyometric training. In: Biomechanics in Sports, edited by A. Barabas, and G. Y. Fabian. Budapest: ITC Plantin, 1995, vol. XII, p. 98-101.
28. Van Leemputte, M., and E. J. Willems. EMG quantification and its application to the analysis of human movements. Med. Sports Sci. 25: 177-194, 1987.
29. Van Zuylen, E. J., A. Van Velzen, and J. J. Denier van der Gon. A biomechanical model for flexion torques of human arm muscles as a function of elbow angle. J. Biomech. 21: 183-190, 1988. [Medline]
30. Williams, L. J., and P. J. Holahan. Parsimony-based fit indices for multiple-indicator models: do they work? Struct. Eq. Mod. 1: 161-189, 1994.

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



This article has been cited by other articles:


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
J. Gerontol. A Biol. Sci. Med. Sci.Home page
A. Ortega-Alonso, N. L. Pedersen, U. M. Kujala, S. Sipila, T. Tormakangas, J. Kaprio, M. Koskenvuo, and T. Rantanen
A twin study on the heritability of walking ability among older women.
J. Gerontol. A Biol. Sci. Med. Sci., October 1, 2006; 61(10): 1082 - 1085.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. M. Roth, E. J. Metter, M. R. Lee, B. F. Hurley, and R. E. Ferrell
C174T polymorphism in the CNTF receptor gene is associated with fat-free mass in men and women
J Appl Physiol, October 1, 2003; 95(4): 1425 - 1430.
[Abstract] [Full Text] [PDF]


Home page
J. Gerontol. A Biol. Sci. Med. Sci.Home page
P. T. Katzmarzyk, N. Gledhill, L. Perusse, and C. Bouchard
Familial Aggregation of 7-Year Changes in Musculoskeletal Fitness
J. Gerontol. A Biol. Sci. Med. Sci., December 1, 2001; 56(12): B497 - 502.
[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 Thomis, M. A.
Right arrow