Linkage analysis of adult height with parent-of-origin effects in the Framingham Heart Study
© Mukhopadhyay and Weeks; licensee BioMed Central Ltd 2003
Published: 31 December 2003
Current linkage analysis methods for quantitative traits do not usually incorporate imprinting effects. Here, we carried out genome-wide linkage analysis for loci influencing adult height in the Framingham Heart Study subjects using variance components while allowing for imprinting effects. We used a sex-averaged map for the 22 autosomes, while chromosomes 6, 14, 18, and 19 were also analyzed using sex-specific maps. We compared results from these four analyses: 1) non-imprinted with sex-averaged maps, 2) imprinted with sex-averaged maps, 3) non-imprinted with sex-specific maps, and 4) imprinted with sex-specific maps. We found four regions on three chromosomes (14q32, 18p11-q21, 18q21-22, and 19q13) with LOD scores above 2.0, with a maximum LOD score of 3.12, allowing for imprinting and sex-specific maps, at D18S1364 on 18q21. While we obtained significant evidence of imprinting effects in both the 18p11-q21 and 19q13 regions when using sex-averaged maps, there were no significant differences between the imprinted and non-imprinted LOD scores when we used sex-specific maps. Our results illustrate the importance of allowing for gender-specific effects in linkage analyses, whether these are in the form of gender-specific recombination frequencies, or in the form of imprinting effects.
Genomic imprinting is the inactivation of a genomic segment depending on the sex of the parent from which it was inherited. Imprinting may play an important role in the inheritance of complex traits in humans. It has been shown by previous studies that incorporation of imprinting into linkage analysis may improve the power of variance component methods for detecting quantitative trait loci [1, 2]. However, Shete and Amos  advise caution in using imprinted analyses, since they require larger sample sizes than non-imprinted analyses. This difference increases as the distance from the quantitative trait locus increases, and as the imprinting effect decreases . On the other hand, where there is a large imprinting effect, using non-imprinted analysis may reduce power to detect linkage.
We decided to analyze linkage of normal adult height because imprinted genes are involved in growth and development . Further, height is easily and reliably measured, and stable within a certain age range. Height is also highly heritable, with h2 values ranging from 0.5 to 0.8 [4–8]. Although our analysis cannot prove that incorporation of imprinting improves linkage analysis, it has the prospect of identifying new regions linked to height. Since sex-specific maps may be more appropriate in this type of analysis , we selected the four chromosomes that produced the highest imprinting LOD scores with sex-averaged maps, and reanalyzed them using sex-specific maps.
The Framingham Heart Study, which started in 1948, is a longitudinal study designed to evaluate the multifactorial nature of cardiovascular disease. The original cohort, collected in 1948, contained adult members from about two-thirds of the homes in Framingham, Massachusetts. In 1971, a second offspring cohort was included in the study, about three-fifths of which are the biological offspring of the original cohort. A more detailed description of these two cohorts is available in previous publications [9, 10]. We analyze here the Framingham Heart Study data as made available to the Genetic Analysis Workshop 13. Our research was approved by the University of Pittsburgh Institutional Review Board (IRB protocol number 020481).
For both cohorts, observers measured and recorded multiple anthropomorphic measurements, at periodic examinations over the duration of the study, at two-year intervals for the original cohort, and at four-year intervals for the offspring cohort. Height measurements were taken with participants standing erect, with heads in the Frankfort plane. For the original cohort, height was measured at biennial examination 1 (in 1948), 5, 10, 13, and then at every subsequent examination onwards. For the offspring cohort, there was approximately an 8-year gap between the first exam (in 1971–1974) and the second exam; subsequent exams were given every 4 years.
Data cleaning and phenotype definition
We examined the 2885 individuals who had multiple height measurements. Each individual's height measurements were compared in a pair-wise fashion for deviations greater than 2.0 inches. Where a single unique measurement was deviant from all others, this measurement was removed, along with its corresponding examination age. There were 22 such individuals. If we were unable to identify a unique deviant, we set that individual's height and age to unknown for the rest of the analysis. Only 12 individuals were zeroed out in this fashion.
For each individual with known height measurements, we first computed a mean adult height measurement between the ages of 20–60 years, and a mean adult age of examination. Linear regression was performed to eliminate gender and age effects. Since there is evidence of change in the population distribution of height in the last few decades, we also regressed mean height against decade-of-birth for each individual. The residuals resulting from the linear regression of height against gender, mean-age, and decade-of-birth were converted into z-scores, which were used as the trait phenotype for linkage analysis.
Variance components with parent-of-origin effects
In order to incorporate imprinting, we have to include separate major-gene components for each parent. If we assume purely additive variance, the variance-covariance matrix Ω is:
Ω = ΠMOσ2aMO + ΠFAσ2aFA + 2Φσ2G + Iσ2e,
where ΠMO is a matrix containing the estimated proportion of maternal alleles shared identically by descent (IBD), ΠFA is a matrix containing the estimated proportion of paternal alleles shared IBD, Φ is the kinship matrix, and I is an identity matrix. Here, we are partitioning the variance into four components: an additive component (σ2aMO) due to the effect of a maternally derived allele linked to the region of interest, an additive component (σ2aFA) due to the effect of a paternally derived allele, a polygenic component (σ2G), and an environmental effect (σ2e). As a simplifying assumption, we ignore the effects of shared environment. We computed five different likelihoods on our pedigree data:
L(σ2aMO = 0, σ2aFA = 0, σ2G, σ2e): No major gene effect. (1)
L(σ2aMO = σ2aFA, σ2G, σ2e): Non-imprinted, where σ2aMO and σ2aFA are constrained to be equal. (2)
L(σ2aMO = 0, σ2aFA, σ2G, σ2e): Complete maternal imprinting, σ2aMO = 0 (3)
L(σ2aMO, σ2aFA = 0, σ2G, σ2e): Complete paternal imprinting, σ2aFA = 0 (4)
L(σ2aMO, σ2aFA, σ2G, σ2e): Imprinted, i.e. σ2aMO, σ2aFA≥ 0 (5)
We tested for linkage at each marker for equations (2–5) by comparing its log-likelihood against that of the polygenic equation (1) using the likelihood ratio test. LOD scores were also computed by dividing the likelihood ratio by 2ln(10). We also examined evidence for imprinting by comparing the log-likelihood of the non-imprinted model to that of the imprinted model.
Linkage analysis was performed on 330 extended pedigrees containing 4692 individuals. SIMWALK2 version 2.83  was used for multipoint IBD estimation  on the complete extended pedigree structures. SIMWALK2 returned sharing estimates for the 15 possible detailed identity states, from which we derived estimates of maternal and paternal IBD sharing for all the parent-offspring pairs and sibling pairs. The variance-component imprinting models described above are currently only applicable to nuclear families; therefore, the extended pedigrees were broken down into their nuclear components for the variance components analyses. Variance components linkage analysis was performed using FISHER  driven by our own imprinting module, which computes the likelihoods described above. All map distances are in Kosambi centimorgans.
Regression and phenotype distribution
Height is dependent on gender, decade-of-birth, and age as follows: βgender = -5.40, βdecade-of-birth = 0.045, and βage = -0.01, with R2 = 0.5521. Gender and decade-of-birth effects are significant at the 0.01 level. The residuals are distributed normally with a mean of 0.00, standard deviation of 2.58 inches, skewness 0.094, and kurtosis = -0.174. The Shapiro-Wilkes test statistic is 0.99 with a p-value of 0.08. After conversion to z-scores, there are only five outliers, two individuals with z-scores less than 3.0, and three with scores larger than 3.0. These individuals were excluded from the linkage analysis.
Linkage analysis results
Regions with LOD scores ≥ 2.0A
Genetic distances (cM)
Sex-Averaged Map LOD Scores
Sex-Specific Map LOD Scores
Marker, Female:Male Ratio
Several genome-wide linkage analysis studies have attempted to pinpoint the genetic factors that control height. A study on Pima Indians showed evidence of a locus on chromosome 20 . Other studies found evidence of linkage of height on chromosomes 6, 7, 12, and 13 , and on chromosomes 7 and 9 . More recently, evidence for linkage was found on chromosome 3 ; and evidence for linkage was observed on chromosomes 6q25, 9p1, and 12q1 ; and on chromosomes 5q31, Xp22, and Xq25 . None of these studies considered imprinting effects or used sex-specific genetic maps. Our study was unable to replicate any of these regions, however, we obtained a LOD score of 1.94 at marker D6S503 in the 6q27 region with imprinting on sex-specific maps, close to regions linked to height in two other studies [15, 18]. It is interesting to note that the estrogen receptor 1 (ESR1) gene is close to this region.
We explore here the impact of carrying out linkage analyses while allowing for sex-specific map differences and imprinting. If the male recombination fraction (θm) is not equal to the female recombination fraction (θf), then a test that assumes θm = θf will be invalid. Furthermore, in a region of true linkage, such map misspecification results in loss of power to detect linkage . Hanson et al.  found by simulation (as confirmed by analytical analyses by Shete and Amos ) that testing for imprinting effects while erroneously assuming θm = θf can increase the type I error rate of the imprinting LOD score, but that these increases are modest if the female:male map-distance ratio is ≤ 5:1. In our study, using sex-averaged maps, we obtained significant evidence of imprinting in precisely those regions with large differences between the male and female genetic distances (Table 1). After we accounted for the sex-specific map differences, then we no longer had significant evidence for imprinting effects, and the most parsimonious explanation of our data is linkage without any imprinting. The linkage signals are strongest when we properly model the male and female recombination fractions by using sex-specific maps. However, it is important to point out that detection of imprinting effects requires large sample sizes, and that region 18p11-q21, where our highest LOD scores were obtained, is known to contain imprinted genes .
It may not be appropriate to ignore dominance effects as we have in our variance components models. In addition, Spencer  shows that the covariance should include an interaction between the dominance and additive terms in the presence of imprinting.
Most previous variance components-based analyses [23, 24] testing for linkage while allowing for imprinting have been limited in that sex-specific maps have not been used in the analyses, mainly due to computational difficulties. The combined use of SIMWALK2  and FISHER  permits incorporation of sex-specific maps into the analyses. As Shete and Amos  point out, "in the imprinting-testing model, it is important to include this difference between the male and female recombination fractions."
The Framingham Heart Study is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with Boston University. This research was supported by the University of Pittsburgh and grant 9R01MH064205-04 from the National Institutes of Health.
- Shete S, Amos CI: Testing for genetic linkage in families by a variance-components approach in the presence of genomic imprinting. Am J Hum Genet. 2002, 70: 751-757. 10.1086/338931.PubMed CentralView ArticlePubMedGoogle Scholar
- Hanson RL, Kobes S, Lindsay RS, Knowler WC: Assessment of parent-of-origin effects in linkage analysis of quantitative traits. Am J Hum Genet. 2001, 68: 951-962. 10.1086/319508.PubMed CentralView ArticlePubMedGoogle Scholar
- Reik W, Walter J: Genomic imprinting: parental influence on the genome. Nat Rev Genet. 2001, 2: 21-32. 10.1038/35047554.View ArticlePubMedGoogle Scholar
- Pearson K, Lee A: On the laws of inheritance in man. I. Inheritance of physical characters. Biometrika. 1903, 2: 357-462. 10.2307/2331507.View ArticleGoogle Scholar
- Roberts DF, Billewicz WZ, McGregor IA: Heritability of stature in a West African population. Ann Hum Genet. 1978, 42: 15-24. 10.1007/BF00291618.View ArticlePubMedGoogle Scholar
- Solomon PJ, Thompson EA, Rissanen A: The inheritance of height in a Finnish population. Ann Hum Biol. 1983, 10: 247-256. 10.1080/03014468300006411.View ArticlePubMedGoogle Scholar
- Furusho T: On the manifestation of genotypes responsible for stature. Hum Biol. 1968, 40: 437-455.PubMedGoogle Scholar
- Kaur DP, Singh R: Parent-adult offspring correlations and heritability of body measurements in a rural Indian population. Ann Hum Biol. 1981, 8: 333-339. 10.1080/03014468100005131.View ArticlePubMedGoogle Scholar
- Kannel WB, Feinleib M, McNamara PM, Garrison RJ, Castelli WP: An investigation of coronary heart disease in families. The Framingham Offspring Study. Am J Epidemiol. 1979, 110: 281-290.PubMedGoogle Scholar
- Dawber TR, Meadors GF, Moore FE: Epidemiological approaches to heart disease: the Framingham study. Am J Public Health. 1951, 41: 279-286.View ArticleGoogle Scholar
- Sobel E, Lange K: Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet. 1996, 58: 1323-1337.PubMed CentralPubMedGoogle Scholar
- Sobel E, Sengul H, Weeks DE: Multipoint estimation of IBD probabilities at arbitrary positions among marker loci on general pedigrees. Hum Hered. 2001, 52: 121-131. 10.1159/000053366.View ArticlePubMedGoogle Scholar
- Lange K, Weeks D, Boehnke M: Programs for pedigree analysis: MENDEL, FISHER, and dGENE. Genet Epidemiol. 1988, 5: 471-472. 10.1002/gepi.1370050611.View ArticlePubMedGoogle Scholar
- Thompson DB, Ossowski V, Janssen RC, Knowler WC, Bogardus C: Linkage between stature and a region on chromosome 20 and analysis of a candidate gene, bone morphogenetic protein 2. Am J Med Genet. 1995, 59: 495-500. 10.1002/ajmg.1320590417.View ArticlePubMedGoogle Scholar
- Hirschhorn JN, Lindgren CM, Daly MJ, Kirby A, Schaffner SF, Burtt NP, Altshuler D, Parker A, Rioux JD, Platko J, Gaudet D, Hudson TJ, Groop LC, Lander ES: Genomewide linkage analysis of stature in multiple populations reveals several regions with evidence of linkage to adult height. Am J Hum Genet. 2001, 69: 106-116. 10.1086/321287.PubMed CentralView ArticlePubMedGoogle Scholar
- Perola M, Ohman M, Hiekkalinna T, Leppavuori J, Pajukanta P, Wessman M, Koskenvuo M, Palotie A, Lange K, Kaprio J, Peltonen L: Quantitative-trait-locus analysis of body-mass index and of stature, by combined analysis of genome scans of five Finnish study groups. Am J Hum Genet. 2001, 69: 117-123. 10.1086/321286.PubMed CentralView ArticlePubMedGoogle Scholar
- Wiltshire S, Frayling TM, Hattersley AT, Hitman GA, Walker M, Levy JC, O'Rahilly S, Groves CJ, Menzel S, Cardon LR, McCarthy MI: Evidence for linkage of stature to chromosome 3p26 in a large U.K. family data set ascertained for type 2 diabetes. Am J Hum Genet. 2002, 70: 543-546. 10.1086/338760.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu J, Bleecker ER, Jongepier H, Howard TD, Koppelman GH, Postma DS, Meyers DA: Major recessive gene(s) with considerable residual polygenic effect regulating adult height: confirmation of genomewide scan results for chromosomes 6, 9, and 12. Am J Hum Genet. 2002, 71: 646-650. 10.1086/342216.PubMed CentralView ArticlePubMedGoogle Scholar
- Deng HW, Xu FH, Liu YZ, Shen H, Deng H, Huang QY, Liu YJ, Conway T, Li JL, Davies KM, Recker RR: A whole-genome linkage scan suggests several genomic regions potentially containing QTLs underlying the variation of stature. Am J Med Genet. 2002, 113: 29-39. 10.1002/ajmg.10742.View ArticlePubMedGoogle Scholar
- Daw EW, Thompson EA, Wijsman EM: Bias in multipoint linkage analysis arising from map misspecification. Genet Epidemiol. 2000, 19: 366-380. 10.1002/1098-2272(200012)19:4<366::AID-GEPI8>3.0.CO;2-F.View ArticlePubMedGoogle Scholar
- Strichman-Almashanu LZ, Lee RS, Onyango PO, Perlman E, Flam F, Frieman MB, Feinberg AP: A genome-wide screen for normally methylated human CpG islands that can identify novel imprinted genes. Genome Res. 2002, 12: 543-554. 10.1101/gr.224102. Article published online before print in March 2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Spencer HG: The correlation between relatives on the supposition of genomic imprinting. Genetics. 2002, 161: 411-417.PubMed CentralPubMedGoogle Scholar
- Lindsay RS, Kobes S, Knowler WC, Bennett PH, Hanson RL: Genome-wide linkage analysis assessing parent-of-origin effects in the inheritance of type 2 diabetes and BMI in Pima Indians. Diabetes. 2001, 50: 2850-2857. 10.2337/diabetes.50.12.2850.View ArticlePubMedGoogle Scholar
- Lindsay RS, Kobes S, Knowler WC, Hanson RL: Genome-wide linkage analysis assessing parent-of-origin effects in the inheritance of birth weight. Hum Genet. 2002, 110: 503-509. 10.1007/s00439-002-0718-2.View ArticlePubMedGoogle Scholar
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.