Articles |
From the Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Tex.
Correspondence to Dr Anthony Comuzzie, Department of Genetics, Southwest Foundation for Biomedical Research, PO Box 760549, San Antonio, TX 78245-0549. E-mail agcom{at}darwin.sfbr.org
| Abstract |
|---|
|
|
|---|
G=0.33 to 0.87) and all were significantly greater than zero (P<.05), indicative of pleiotropic effects. However, we found substantial unique genetic variance for each of these traits even after accounting for the effects shared in common with all the remaining measures. We conclude that the genetic variation in these HDL phenotypes is a result of the action of common as well as unique genes.
Key Words: pleiotropy HDL apolipoproteins quantitative genetics multivariate
| Introduction |
|---|
|
|
|---|
Given the physical and functional relationships among these HDL measures, we can hypothesize that many of the observed phenotypic correlations are in fact due to an underlying structure of common genetic effects (ie, pleiotropy), with this situation being manifest as a pattern of significant
G among all of these various HDL measures. To test this hypothesis, we conducted a multivariate quantitative genetic analysis in which we decomposed the phenotypic correlations among a group of seven related HDL measures into the portion due to genes shared in common and that due to shared environment. By use of this approach, it is possible to directly assess the extent of pleiotropy underlying the relationships among HDL measures. At the same time, however, it is also possible to use this information to remove the common genetic influences on each of these measures so as to allow for the identification of the genetic elements unique to each of these measures.
| Methods |
|---|
|
|
|---|
At the clinic visit, a blood sample was obtained after an over-night fast. Plasma was divided into aliquots2 and stored at -80°C until used. In addition, information about age, sex, health status, and medications was also collected. An Institutional Review Board (University of Texas Health Science Center at San Antonio) approved the procedures, and all subjects gave informed consent.
Phenotypes
HDL Concentration Measures
ß-Lipoproteins were precipitated from plasma samples by use of dextran sulfate,3 and HDL cholesterol concentrations were assayed enzymatically in a Gilford SBA-300 Clinical Chemistry Analyzer with commercial reagents supplied by Boehringer-Mannheim Diagnostics. HDL-UC concentrations were measured by use of a similar kit (Wako) except that the enzyme cholesterol esterase was not added; the concentrations of HDL-EC were calculated as the difference between HDL-C and HDL-UC. The interassay coefficients of variation for control products in these assays were 6.2% for HDL-C and 9.0% for HDL-UC. Apolipoprotein concentrations were determined by a commercial laboratory (Medical Research Laboratories). Apo A-I concentrations were determined by nephelometry,4 5 and apo A-II concentrations were determined with competitive immunoassays.6 The interassay coefficients of variation for control products in these assays were 3.5% for apo A-I and 4.4% for apo A-II.
HDL Particle Size Measures
Size distributions of HDL particles were determined by electrophoresis of plasma in nondenaturing 3% to 31% polyacrylamide gradient gels, which were made in the laboratory as previously described.7 After electrophoretic separation, the proteins were transferred to nitrocellulose paper and detected by binding with sheep antiapo A-I and sheep antiapo A-II antibodies (both from Boehringer-Mannheim) as previously described.8 9 These antibodies were in turn bound by donkey anti-sheep IgG (Chemicon International, Inc), which was radioiodinated by the chloramine T method,10 and distributions were measured by densitometry of autoradiograms with an LKB Ultroscan Laser Densitometer with GSXL software (Pharmacia). The distribution of cholesteryl esters among HDLs was detected by staining with Sudan black B and densitometry as described.11 12 HDL absorbance profiles were decomposed by a curve-fitting procedure as suggested previously13 ; we developed in house a program to automatically fit curves representing the HDL subclasses and to calculate fraction of absorbance in each of the HDL subclasses. We summed the data to generate fractional absorbance for HDL3 particles and for those larger than HDL3. To avoid the problems of singularity for proportional data (ie, the two fractions must sum to one and are not independent), we used a log ratio transformation as suggested by Aitchison,14 in which the fractional absorbance for larger HDLs was divided by the fractional absorbance for HDL3. Thus, each of the size distributions was summarized into a single variable, appropriate for statistical analyses, which was positively correlated with average size of HDLs; these variables are HDL-Csize, HDL-A-Isize, and HDL-A-IIsize.
Covariates
In addition to the seven HDL measures, a number of covariate effects were also included in the analysis. Since all of these traits had previously been screened for a variety of covariate effects in univariate analyses, the multivariate analysis included only covariates that had previously been found to be significant for at least one of the seven traits examined. In addition, all seven HDL measures had previously been screened for shared household effects, but none were detected. The covariate effects that were included in the multivariate analysis were sex, sex-specific age, hormone replacement therapy, diabetic status, and use of diabetic medications. It must be remembered, however, that no model is a true reflection of the natural world, and undoubtedly as yet unidentified covariates exist that could have a significant effect. However, the inclusion of such environ-mental covariates would only serve to strengthen the results (ie, increase the relative genetic contribution) by adding further resolution to the model.
Analytic Methods
Variance Decomposition
By use of established quantitative genetic theory, it is possible to extend a univariate genetic analysis to encompass the multivariate state.15 16 17 18 Variance decomposition techniques16 17 18 using maximum-likelihood methods and implemented in a modified version of the computer program PAP19 were used to simultaneously decompose the phenotypic correlations among seven measures of size and concentration of HDL into their underlying additive
G and random
E. This decomposition was accomplished by conditioning the phenotypic correlations on the degree of relationship among all the individuals in the sample by use of extended pedigree information. Specifically, the parameterization of the additive genetic model used here allowed for the simultaneous estimation of mean effects, six covariate effects, phenotypic SDs, and heritabilities for each trait and estimates of
G and
E for all pairwise combinations of the traits.
The phenotypic correlation (
P) between a pair of traits within an individual can be expressed in terms of the underlying
G and
E, correcting for the use of related individuals, by the equation (Equation 1
) where h12 is the heritability of trait 1 and h22 is the heritability of trait 2. The heritability (h2) represents the portion of the phenotypic variance accounted for by the total additive genetic variance (h2=
2G/
2P), or in other words, the portion of the similarity in a given trait between related individuals that is predicted by degree of relatedness. The phenotypic correlation between a pair of traits is expressed as a function of the shared genetic and environmental effects (expressed as the
G and
E) between a pair of traits weighted by the portion of their respective phenotypic variances accounted for by the effect of genes (ie, their heritabilities).
To estimate the additive
G and the random
E for pairs of traits, the multivariate phenotype of an individual is first modeled as a linear function of the individual's trait measurements, and the matrix of kinship coefficients is generated, expressing relatedness among all pairs of individuals in the pedigree. Using standard quantitative genetic theory, the phenotypic variance-covariance matrix and its additive genetic and random environmental components are then obtained. From these matrices,
G and
E are estimated directly. A significant nonzero
G is a direct measure of the extent of genes shared in common between a pair of traits (ie, pleiotropy). A more detailed explanation of the extension of this methodology to the multivariate state can be found in previously published work.15 20 21
To examine the unique genetic component of each of these seven traits independent of all the others, we calculated the conditional phenotypic, genetic, and environmental covariance matrices using standard matrix formulas employing the observed maximum-likelihood parameter estimates of the full model in which all the traits were included as separate phenotypes, with these conditional estimates themselves being maximum-likelihood estimates. We then obtained estimates of residual heritabilities (ie, heritabilities of each of these traits after the common genetic and shared environmental effects of the remaining six traits were removed). In other words, the heritabilities, along with the variance, for each of the seven traits were reestimated while their common interactions with the six remaining traits were accounted for.
Even though the multivariate model used in this study allows only for additive genetic effects, it is possible that major genes could be involved in determination of the common variation shared among these seven phenotypes, which would lead to deviations from multivariate normality. The maximum-likelihood methods used here, however, are robust to such deviations in the underlying distribution. Therefore, valid maximum-likelihood estimates for the parameters of the genetic model can be obtained.22
Log-Likelihood Ratio Test
The significance of both the
G and
E and the estimates of heritability and covariate effects were tested by comparing the ln likelihoods from restricted models, in which each of these parameters was in turn constrained to equal zero, against the ln likelihood for the general model in which all parameters were estimated. The ln likelihood values of the general and the restricted models were compared by use of a likelihood ratio test. This test yields a statistic that is distributed asymptotically as a
2 with df equal to the difference in the number of parameters estimated in the two models being compared and is calculated as -2(ln likelihoodrestricted model-ln likelihoodgeneral model). In this case, the comparisons of the restricted models to the general model have 1 df. An ln likelihood score for either of the restricted models that was significantly worse than that of the general model (P
.05) was considered to be evidence of a significant, nonzero, correlation, heritability, or covariate effect.
| Results |
|---|
|
|
|---|
.05) covariate effects as simultaneously estimated in the full heptavariate model. These parameters include the male mean (µmale), a sex effect (ßsex), male age effect (ßagemale), female age effect (ßagefemale), effect of exogenous hormone use (ßhormone), effect of diabetes status (ßdiabetes), and the effect of diabetic medication (ßdmeds). In the case of apo A-I, apo A-II, and HDL-UC, there are no significant sex effects, so the male mean is in fact the MLE for the mean in both sexes. With respect to measures of HDL-EC, HDL-A-Isize, HDL-A-IIsize, and HDL-Csize, there are significant sex effects (P
.05), so the female means for these traits are obtained by adding 0.10, 0.20, 0.14, and 0.21, respectively, to the MLE for the male mean for the corresponding trait. Significant effects of age were detected for three measures of HDL concentration: apo A-I, apo A-II, and HDL-EC. Significant effects for hormone replacement therapy were detected for apo A-I, apo A-II, and HDL-UC. A significant effect of diabetic status was detected for HDL-EC, and significant effects of diabetes (ßdiabetes) and use of diabetic medications (ßdmeds) were detected for HDL-A-Isize.
|
SDs and Heritabilities
Table 2
provides the maximum-likelihood estimates and associated standard errors for the phenotypic SDs (
) and heritabilities (h2) for each of the seven HDL measures in the heptavariate analysis. The SD estimates range from 0.17 for apo A-I to 0.47 for HDL-Csize. The heritability estimates range from 0.17 for HDL-A-IIsize to 0.61 for apo A-I, with all seven traits having heritabilities significantly greater than zero at P
.05.
|
Genetic Correlations
Table 3
provides the additive
G among these seven HDL measures along with their associated standard errors. As indicated in the table, the
G are all positive, and all are significantly greater than zero at P
.05. These additive
G range from a moderate value of 0.33 (apo A-II:HDL-A-Isize) to a relatively large value of 0.87 (HDL-UC:HDL-Csize).
|
Environmental Correlations
Table 4
provides the random
E among these seven HDL measures. Only the random
E between apo A-I:HDL-A-Isize and apo A-II:HDL-UC were not significantly different from zero at P
.05. Unlike the pattern of
G, the random
Evary in sign (both positive and negative correlations), and the significant negative random
E were detected only between apo A-II and the three HDL size measures, with values of -0.33, -0.22, and -0.36 between apo A-II:HDL-A-Isize, apo A-II:HDL-A-IIsize, and apo A-II:HDL-Csize, respectively. The significant positive random
E ranged from 0.18 (apo A-I:HDL-A-IIsize) to 0.65 (HDL-A-IIsize:HDL-Csize).
|
SDs and Heritabilities of the Unique Genetic Component
Table 5
provides the maximum-likelihood estimates and associated standard errors for the phenotypic SDs (
) and heritabilities (h2) for each of the seven HDL measures in the heptavariate analysis after removing the common effects of the remaining six traits. The SDs ranged from 0.13 for apo A-I to 0.29 for HDL-Csize. The heritability estimates range from 0.10 for HDL-UC to 0.31 for apo A-I. For all but two of the traits (HDL-UC and HDL-A-IIsize), residual heritabilities were significantly greater than zero at P<.05.
|
Percent Variance Remaining After Corrections for Shared Genetic Effects
Table 6
shows the percentage of the additive genetic, random environmental, and phenotypic variance remaining after correction for the shared effects with the other six traits. Between 35% and 71% of the total phenotypic variance remained after each trait was adjusted for the pleiotropic effects shared with the other six traits. Table 6
also gives the percentage of the additive genetic and random environmental variances remaining after adjustment for the shared genetic effects. As would be expected, a substantial amount of random environmental variance remains after this adjustment. However, there is a marked decrease in the additive genetic variance of each of the seven traits, with percentages of additive genetic variance remaining ranging from 13% (HDL-UC) to 36% (HDL-A-IIsize).
|
| Discussion |
|---|
|
|
|---|
G among these seven HDL measures, we can conclude the existence of a substantial degree of pleiotropy (ie, additive genetic effects common to all seven traits). The finding of significant heritabilities for all of the phenotypes examined, together with the pattern of significant positive additive
G among all pairs of these phenotypes, implies that the expression of these traits to varying degrees is under genetic control and that those genes responsible for elevated values for one of these traits also leads to higher values for the other traits. It should also be emphasized that because a preliminary screening of these seven HDL traits failed to detect any significant common household effect, these heritabilities and the
G reflect solely the action of genes and are not confounded by nongenetic familial effects. The magnitudes of the heritabilities for each of these traits and the
G among them indicate that the phenotypic intercorrelation observed among this group of variables is largely driven by the common influence of additive effects of a set of shared genes (ie, pleiotropy). Despite these findings of a high degree of commonality with respect to genetic effects among this set of seven HDL phenotypes, there is no evidence of complete pleiotropy (ie,
G in each case is significantly less than one). This means that there are additional genes that uniquely exert significant effects on each of the traits (and subsets of them). This point is further supported by the fact that there are relatively pronounced residual genetic variances for each of these traits after adjustment for the effects of pleiotropy. In the case of HDL-UC and HDL-A-IIsize, however, with their nonsignificant heritabilities after correction for pleiotropy, it would appear that perhaps the bulk of their genetic component is derived from the genes common among all seven of the traits.
Although the overall values for the random
E are not quite as large as those for the
G, there is still clear evidence of shared random environmental effects among the majority of this set of phenotypically related HDL traits. The significant negative random
E between apo A-II and the three measures of particle size (HDL-A-Isize, HDL-A-IIsize, and HDL-Csize) indicate that the random environmental factors leading to increased values of apo A-II in turn result in lower values in each of the three size measures. In contrast, all of the significant positive
E among this set of traits indicate that those random environmental factors leading to increased values in one trait also lead to increased values in each of the others.
We found significant random
E for most traits. These indicate the existence of other important covariates that were not included in this model or that are still unmeasured; however, the identification of these environmental covariates would not decrease the genetic effects already detected. Likely candidates for these random environmental covariates include diet and endocrine hormones. Age and sex are well established factors associated with lipoprotein variation. Diabetes has been associated with low concentrations of HDL cholesterol, and in this study we found a significant effect for HDL-EC but not for HDL-UC or for apo A-I, the principal apolipoprotein of HDL. Recently, we reported that diabetic individuals had significantly smaller HDL particles when lipoprotiens were stained for apo A-I distributions but not when stained for cholesterol distributions.23 The present study, including a substantially larger data set, confirms these earlier results and reveals a significant counteracting effect of diabetic medications on HDL size phenotypes.
In general, these findings reveal a strong pattern of
G among all seven of these traits, which is indicative of pronounced pleiotropy and suggests that at some level, each of them is in fact a proxy for a larger and perhaps as yet undefined phenotypic entity of HDL. Given the results of this analysis, it would be ill advised to treat these seven phenotypes as truly independent measures of different aspects of HDL without first accounting for their rather significant shared genetic and environmental components. However, we have also shown that a significant unique genetic component exists for each of these traits after the removal of the shared effects. Therefore, the next task is to identify ways to exploit both the common information inherent in each of these seven measures and those that would allow for the identification of the unique elements that are also inherent in each of these phenotypes. Recognition of such genetic and environmental interactions and their ultimate disentanglement will further the understanding of this group of lipoprotein variables and of the normal metabolic and pathogenic processes in which they are involved.
|
| (1) |
| Selected Abbreviations and Acronyms |
|---|
|
| Acknowledgments |
|---|
Received July 9, 1996; accepted October 2, 1996.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. E. Cassidy-Bushrow, L. F. Bielak, P. F. Sheedy II, S. T. Turner, I. J. Kullo, X. Lin, and P. A. Peyser Coronary Artery Calcification Progression Is Heritable Circulation, July 3, 2007; 116(1): 25 - 31. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Menzaghi, A. Coco, L. Salvemini, R. Thompson, S. De Cosmo, A. Doria, and V. Trischitta Heritability of Serum Resistin and Its Genetic Correlation with Insulin Resistance-Related Features in Nondiabetic Caucasians J. Clin. Endocrinol. Metab., July 1, 2006; 91(7): 2792 - 2795. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. L. Rainwater, L. J. Martin, and A. G. Comuzzie Genetic Control of Coordinated Changes in HDL and LDL Size Phenotypes Arterioscler. Thromb. Vasc. Biol., November 1, 2001; 21(11): 1829 - 1833. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. G. Comuzzie, T. Funahashi, G. Sonnenberg, L. J. Martin, H. J. Jacob, A. E. K. Black, D. Maas, M. Takahashi, S. Kihara, S. Tanaka, et al. The Genetic Basis of Plasma Variation in Adiponectin, a Global Endophenotype for Obesity and the Metabolic Syndrome J. Clin. Endocrinol. Metab., September 1, 2001; 86(9): 4321 - 4325. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. L. Rainwater, B. D. Mitchell, M. C. Mahaney, and S. M. Haffner Genetic Relationship Between Measures of HDL Phenotypes and Insulin Concentrations Arterioscler. Thromb. Vasc. Biol., December 1, 1997; 17(12): 3414 - 3419. [Abstract] [Full Text] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
ATVB Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 1997 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |