# Power of maximum HLOD tests to detect linkage to obesity genes

- Yun Joo Yoo
^{1}Email author, - Yanling Huo
^{1}, - Yuming Ning
^{2}, - Derek Gordon
^{3}, - Stephen Finch
^{1}and - Nancy R Mendell
^{1}

**4(Suppl 1)**:S16

https://doi.org/10.1186/1471-2156-4-S1-S16

© Yoo et al; licensee BioMed Central Ltd 2003

**Published: **31 December 2003

## Abstract

### Background

We investigate the power of heterogeneity LOD test to detect linkage when a trait is determined by several major genes using Genetic Analysis Workshop 13 simulated data. We consider three traits, two of which are disease-causing traits: 1) the rate of change in body mass index (BMI); and 2) the maximum BMI; and 3) the disease itself (hypertension). Of interest is the power of "HLOD2", the maximum heterogeneity LOD obtained upon maximizing over the two genetic models.

### Results

Using a trait phenotype Obesity Slope, we observe that the power to detect the two markers closest to the two genes (S1, S2) at the 0.05 level using HLOD2 is 13% and 10%. The power of HLOD2 for Max BMI phenotype is 12% and 9%. The corresponding values for the Hypertension phenotype are 8% and 6%.

### Conclusion

The power to detect linkage to the slope genes is quite low. But the power using disease-related traits as a phenotype is greater than the power using the disease (hypertension) phenotype.

## Keywords

## Background

One of the issues in analysis of complex diseases is the power of considering an disease related trait as a phenotype than the disease itself [1]. Another issue considered here is the method of analysis, which allows for heterogeneity when one or more major genes determine a trait [2]. The Genetic Analysis Workshop 13 (GAW13) simulated longitudinal data on obesity provide an example of a disease-related phenotype. According to the simulation model, the Obesity Slope, the rate of increase in BMI (body mass index) over time is determined by only two major genes (S1, S2). Also, Obesity Slope in turn is a risk factor for hypertension, and it is also implied in the simulation model that these two genes for obesity are in fact two of many hypertension-causing genes.

In this work, we want to investigate the power of HLOD maximized over two genetic models (dominant and recessive) to detect linkage with the Obesity Slope phenotype and also with the hypertension phenotype, i.e., the disease itself. We also consider another trait, the individual's maximum observed BMI (Max BMI), which is related to both obesity and hypertension.

## Methods

### Definition of disease phenotypes

The first phenotype considered is Obesity Slope. For each individual, we regressed BMI values on log of age values at each time where the BMI was recorded, using the following model: E(BMI _{ij}) = α_{j} + β_{j}log_{10}(X_{ij}), for i = 1,2,..., t_{j}.

Here X_{ij} denotes the j^{th} individual's age at the i^{th} measurement; BMI_{ij} denotes the j^{th} individual's BMI at age X_{ij}; and t_{j} denotes the number of time points for which we have BMI values on the j^{th} individual. The value of BMI_{ij} is computed, again following the model described in the simulation as BMI_{ij} = W_{ij}/H_{ij}^{2}, where W_{ij} denotes j^{th} individual's weight in kilograms and H_{ij} denotes the j^{th} individual's height in meters at the i^{th} time point observed. Ordinary least-squares regression resulted in estimates of β_{j} for each of our subjects. We treated an individual as "missing" slope information if he had too few measurements. Individuals in Cohort 1 with less than three measurements were treated as missing. Individuals in Cohort 2 with less than two measurements were treated as missing. Cohort 1 has 21 time points with two-year intervals and Cohort 2 has five time points with eight- or four-year intervals. We then defined individuals as "Obese" if β_{j} ≥ b_{90} and "Normal" if β_{j} < b_{90}, where b_{90} denotes the 90^{th} percentile of the slopes observed in each replicate being analyzed.

The second disease-related phenotype considered here is the maximum of j^{th} individual's BMI over all ages (Max BMI). We defined individuals as "High" if Max BMI ≥ m_{90} and "Normal" if Max BMI < m_{90}, where m_{90} denotes the 90th percentile of the Max BMI values observed in each replicate being analyzed. The third phenotype is "Hypertension." Individuals defined as "affected" are those with diastolic BP ≥ 90 and/or systolic BP ≥ 140, and "Normal" are those having diastolic BP < 90 and systolic BP < 140.

### Methods of linkage analysis

Single-point linkage analysis was done for the two markers closest to the two genes S1 and S2, and for five unlinked markers. One marker linked to S1 is 11g6 and the other, linked to S2, is 7g7. We also considered five unlinked markers on chromosomes 2, 4, 6, 8, and 10 (2g5, 4g4, 6g3, 8g2, 10g1), on which none of the disease genes were located.

We used the linkage test specified on Abreu et al. [3]. We obtained a value "HLOD-D" by maximizing over all values of the recombination fraction and heterogeneity parameter assuming disease allele frequency 0.01 and dominance with penetrance 0.5. Then we obtained a value "HLOD-R" by maximizing over all values of the recombination fraction and heterogeneity parameter assuming disease allele frequency 0.01 and recessive model with penetrance 0.5. The test statistic HLOD2 was obtained as the maximum of HLOD-D and HLOD-R. Also, for the comparison, we considered LOD2, the maximum of LOD-D and LOD-R using the same two analysis models without heterogeneity parameter. To calculate HLOD values, we used HOMOG program combining with LINKAGE package.

## Results and Discussion

*N*= 500. Then we obtained power values for the three phenotypes by computing the percentage of replicates with higher HLOD2 (or LOD2) values than the critical value in all 100 replicates. Highest power is 13% for Obesity Slope phenotype at the 11g6 marker. Power values for three phenotypes are presented in the Table 1.

Power of HLOD2 and LOD2 test for three phenotypes Critical values of 0.947 for HLOD2 and 0.730 for LOD2 are used for both of markers.

Marker | 11g6 (S1) | 7g7 (S2) | ||
---|---|---|---|---|

Phenotype | HLOD2 | LOD2 | HLOD2 | LOD2 |

Obesity Slope | 13% | 10% | 10% | 10% |

Max BMI | 12% | 11% | 9% | 9% |

Hypertension | 8% | 8% | 6% | 4% |

^{2}= 19.5). On the other hand, there is no association in this sample between being in the top 10% for the Slope Phenotype and having hypertension (χ

^{2}= 2.66).

There are several possible explanations for the low power to detect linkage observed here. The most likely one is that we lose power by dichotomizing the values rather than analyzing these traits as quantitative phenotypes using model-free methods. Second, even though we are focusing on the slope itself, which is the direct product of the gene, we do have error in our estimate of an individuals' slope due to the combination of random noises in the simulations models.

## Conclusions

The overall power values for HLOD2 analysis are quiet low in the case of the analysis of dichotomized BMI slope values and dichotomized BMI values. However, they are higher than those obtained using the hypertension phenotype.

## Authors’ Affiliations

## References

- Rice JP, Saccone NL, Rasmussen E: Definition of the phenotype. Adv Genet. 2001, 42: 69-76.View ArticlePubMedGoogle Scholar
- Greenberg DA, Abreu PC: Determining trait locus position from multipoint analysis: accuracy and power of three different statistics. Genet Epidemiol. 2001, 21: 299-314. 10.1002/gepi.1036.View ArticlePubMedGoogle Scholar
- Abreu PC, Hodge SE, Greenberg DA: Quantification of type I error probabilities for heterogeneity LOD scores. Genet Epidemiol. 2002, 22: 156-169. 10.1002/gepi.0155.View ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.