|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Journal of Lipid Research, Vol. 45, 1876-1884, October 2004
Copyright © 2004 by American Society for Biochemistry and Molecular Biology





,**
* Department of Molecular Medicine, National Public Health Institute, Helsinki, Finland
Department of Medicine, Helsinki University Central Hospital, Helsinki, Finland
Department of Human Genetics, David Geffen School of Medicine at the University of California Los Angeles, Los Angeles, CA
** Department of Medical Genetics, University of Helsinki, Helsinki, Finland
Published, JLR Papers in Press, July 16, 2004. DOI 10.1194/jlr.M400141-JLR200
1 To whom correspondence should be addressed. e-mail: leena.peltonen{at}ktl.fi (L.P.); ppajukanta{at}mednet.ucla.edu (P.P.)
| ABSTRACT |
|---|
|
|
|---|
These analyses provided evidence for linkage (a logarithm of odds score of 3.2) on 10q11 using a quantitative HDL-C trait. Importantly, this region, previously linked to TGs, body mass index, and obesity, provided evidence for association for quantitative TGs (P = 0.0006) and for a combined trait of HDL-C and TGs (P = 0.008) with marker D10S546. Suggestive evidence for linkage also emerged for HDL-C on 2q31 and for TGs on 20q13.32. Finnish families ascertained for dyslipidemias thus suggest that 10q11, 2q31, and 20q13.32 harbor loci for HDL-C and TGs.
Abbreviations: apoA-I, apolipoprotein A-I; BMI, body mass index; CHD, coronary heart disease; cM, centimorgan; FCHL, familial combined hyperlipidemia; HDL-C, high density lipoprotein-cholesterol; HDL2 and HDL3, high density lipoprotein cholesterol particles 2 and 3; LDLps, low density lipoprotein cholesterol particle size; QTL quantitative trait locus; TC, total cholesterol; TG, triglyceride
Supplementary key words low high density lipoprotein-cholesterol familial combined hyperlipidemia coronary heart disease atherosclerosis quantitative traits variance component analyses quantitative trait locus triglycerides type 2 diabetes metabolic syndrome
| INTRODUCTION |
|---|
|
|
|---|
50% of familial dyslipidemias predisposing to premature coronary heart disease (CHD) (1). Both are typical complex disorders, influenced by several environmental and genetic factors. FCHL is characterized by increased levels of serum total cholesterol (TC) and triglycerides (TGs), or both, and also expresses low HDL-C as a component trait (24). It has been estimated that the genetic component accounts at least 50% of the variation in HDL-C levels (5). Although several genes have been identified for rare forms of the low HDL-C trait, the major loci affecting plasma HDL-C levels and explaining most of the HDL-C variation have not been identified (6). We previously performed a genome-wide scan in Finnish low HDL-C families (3). The identified loci for low qualitative HDL-C were also tested in an independent study sample of Finnish FCHL families. In the pooled data analysis of these two study samples, evidence for linkage emerged for low HDL-C on 8q23.1, 16q23.3, and 20q13.32, the markers on 8q23.1 showing the most significant signal with a two-point logarithm of odds (LOD) score of 4.7 (3). All three of these loci have also been linked to HDL-C or to other CHD risk factors in other populations (714). The loci on 8q23.1 and 16q have been linked to quantitative HDL-C in Mexican Americans (7, 12). The locus on 20q12-13 seems to have a role in the development of other component traits of the metabolic syndrome, and the contribution of the same region to the development of type 2 diabetes mellitus (MIM 125853) has been shown in several studies (811, 13, 14). The locus on 20q13 has also been linked to obesity (MIM 601665) (11). In addition to these three loci (8q23.1, 16q23.3, and 20q13), we have earlier identified loci for TGs on 10q11 (15) and for premature CHD on 2q31 and Xq24 (16). In these previous genome scans (3, 15, 16), our analysis strategy focused on the information obtained from the affected family members using dichotomized affection criteria, whereas in the present study, the aim was to use quantitative lipid phenotypes and to analyze these six chromosomal regions as potential quantitative trait loci (QTLs) for HDL-C and other lipid traits characteristic of CHD and the metabolic syndrome. For this purpose, we more than tripled the number of genotyped individuals and analyzed an extended study sample consisting of 1,109 genotyped family members from 92 multigenerational dyslipidemic families for a dense set of markers in each chromosomal region.
| MATERIALS AND METHODS |
|---|
|
|
|---|
In the Finnish HDL-C families, there were 168 individuals affected by low HDL-C and 49 individuals affected by high TGs using the 10th and 90th age/sex-specific Finnish population percentiles for HDL-C and TGs. Inclusion criteria for the low HDL-C probands were age of 3060 years for men and women, HDL-C levels below the 10th age/sex-specific Finnish population percentile, and CHD verified by either coronary angiography (>50% stenosis in one or more coronary arteries) or myocardial infarction (3). Additional lipid criteria for the probands were TC of
6.3 mmol/l in men and
6.0 mmol/l in women and TGs of
2.3 mmol/l in both genders. Exclusion criteria for the low HDL-C probands were type 1 and type 2 diabetes mellitus, hepatic or renal disease, and body mass index (BMI) > 30. If the criteria mentioned above were fulfilled, families with at least two affected members were included in the study and all accessible family members were examined. The affected family members were ascertained for low HDL-C using the Finnish age/sex-specific 10th population percentiles (19, 20). The age/sex-specific HDL-C and TG percentiles are publicly available for the Finnish population (18) (http://www.ktl.fi/molbio/wwwpub/fchl/genomescan/link4.html).
In the Finnish FCHL families, 214 individuals were affected by low HDL-C and 198 by high TGs using the 10th and 90th age/sex-specific Finnish population percentiles for HDL-C and TGs. The inclusion criteria for the probands of FCHL families were as follows: 1) FCHL lipid phenotype (serum TC and/or TGs > 90th age/sex-specific Finnish population percentiles); 2) age of >30 years and <55 years for males and <65 years for females; and 3) at least 50% stenosis in one or more coronary arteries by coronary angiography. Exclusion criteria for the FCHL probands were type 1 diabetes mellitus, hypothyroidism, and hepatic or renal disease (15, 17).
Biochemical analyses and classification of the affection status
Serum TC and TGs were determined with an automated Cobas Mira analyzer (Hoffman-La Roche, Basel, Switzerland) by enzymatic methods (Hoffman-La Roche kits 0722138 and 0715166, respectively). Serum HDL-C was quantified by phosphotungstic acid/magnesium chloride precipitation procedures (Hoffman-La Roche kit 0720674). Concentrations of apolipoprotein A-I (apoA-I), apoA-II, and apoB were measured by immunoturbidimetric methods with commercial kits (Boehringer-Mannheim, Mannheim, Germany). HDL-C particles 2 and 3 (HDL2 and HDL3) were separated by ultracentrifugation, as described (21). LDL particle size (LDLps) was determined using native gradient gel electrophoresis, as described in detail by Vakkilainen et al. (22). BMI was calculated as weight/height2. The clinical and biochemical data in the Finnish low HDL-C and FCHL families are fully compatible, because the data collection, laboratory measurements, and phenotype determinations for these two study samples were performed in the same laboratory at the Helsinki University Central Hospital. Patients who used lipid-lowering drugs were studied after their lipid-lowering treatment was withheld for 4 weeks.
Genotyping of the markers
A total of 67 microsatellite markers were genotyped for six regions (2q31, 8q23.1, 10q11, 16q23.3, 20q13.32, and Xq24) in 1,109 family members: 426 individuals from 39 low HDL-C families and 683 individuals from 53 FCHL families. On average, 11 markers were genotyped per locus. When selecting the markers, we aimed to cover
10 Mb around peak markers in every selected chromosomal region to obtain an average intermarker distance of 1 Mb. The selection of the markers, as well as the order and distances of the markers for the six regions, were based on the genetic maps of the Decode map (23), Celera (http://www.celera.com), the Marshfield comprehensive human genetic map (http://research.marshfieldclinic.org/genetics/), and the Human Genome Browser Gateway (http://www.genome.ucsc.edu). The April 2003 version of the Human Genome Browser Gateway was used for chromosomal band determination. Therefore, the chromosomal bands of the loci analyzed in the current study differ slightly from the original bands reported in the previous publications (3, 15, 16).
Genotyping was performed at the Finnish Genome Center (Helsinki, Finland). The fluorescence-labeled PCR products were pooled into panels and the pooled samples were dialyzed before electrophoresis, which was performed on a MegaBace 1,000-capillary electrophoresis instrument (Molecular Dynamics). Alleles were defined by using Genetic Profiler version 1.1 software (Molecular Dynamics). The data were checked for Mendelian errors with Pedmanager (M. P. Reeve, unpublished data) (http://www.genome.wi.mit.edu/ftp/distribution/software/pedmanager/) and Pedcheck software (24).
Statistical analyses
Quantitative linkage analyses
To investigate the earlier identified regions using quantitative measures by variance component methods in extended families, QTL analyses were conducted using maximum likelihood-based approaches implemented in the computer program SOLAR version 1.7.4 (25). The significance of age, sex, and BMI as covariates was tested in our analyses of HDL-C and TGs. For HDL-C, sex and BMI were significant covariates. For TGs, age, sex, and BMI were significant. These significant covariates were included in the linkage analyses at the first stage. We also further investigated the relationship between TGs and HDL-C by including TGs as a covariate for the quantitative trait analyses of HDL-C, and vice versa, for 10q11. To correct for the study ascertainment, the proband status was scored as a dichotomized variable and included in the SOLAR analyses. Individual TG and apoB levels were log transformed to reduce skewness and kurtosis. No other transformations were needed before statistical analyses.
Qualitative linkage analyses Qualitative analyses were performed using both parametric linkage and nonparametric affected sibpair analyses. Nonparametric analyses were used to better model the unknown mode of inheritance of the traits. For the parametric linkage analyses, the MLINK program of the linkage package (26) version FASTLINK 4.1P (27, 28) was selected. The identical-by-descent status of affected sibpairs was assessed with the help of the SIBPAIR program (29) of the ANALYZE package (30). The parametric linkage analyses were performed with dominant and recessive modes of inheritance. Gene frequencies of 0.4% and 8%, reflecting an estimated population prevalence of 1%, were used for the dominant and recessive modes of inheritance of the HDL-C trait, as described previously (3, 15, 17, 18). For TGs, gene frequencies of 0.6% and 11% were adopted for the dominant and recessive modes of inheritance (15, 17). Linkage analyses were also performed on the component nuclear families. To circumvent problems of incomplete penetrance and genetic ambiguity of the phenotype "unaffected," we used the affecteds-only strategy and coded the subjects either affected or unknown based on their age/sex-specific 10th percentile of measured HDL-C and 90th percentile of TGs (3, 15). The unaffected individuals were genotyped to increase phase information and were treated as if their phenotype was unknown in statistical analyses. For each marker, the allele frequencies were estimated from all individuals with the DOWNFREQ program (31). Genetic heterogeneity between the families was examined using the HOMOG program (32).
Data for the 10q11 region were also analyzed with the QTL program MERLIN (33) to compare the results obtained with SOLAR with another QTL-based program. Furthermore, the most promising region, 10q11, was analyzed with a novel QTL association method included in the software package Mendel (34). Although the method and its implementation will be described at length elsewhere (K. Lange and E. Sobel, unpublished data), we can state the main ideas briefly here. These revolve around the "measured genotype" approach, which treats the genotypes at a candidate locus as linear predictors modifying the mean of an associated quantitative trait (3537). The variance component context of measured genotype models allows one to control for random environment and polygenic background while testing for significant mean effects by maximum likelihood estimation and likelihood ratio tests. Mendel implements this strategy for a) multivariate traits, b) both random and ascertained pedigrees, and c) noncodominant markers. Ascertainment is handled by conditioning on proband trait values. Noncodominant markers are dealt with by imputing marker allele counts for each person conditional on the marker phenotypes observed throughout his or her pedigree. These counts may be fractional in the presence of missing data. In our view, the measured genotype approach to association testing is one of the most powerful approaches available. It does require normally distributed traits and can be confused by ethnic stratification. Fortunately, the latter problem can be largely ignored in the Finnish data. Another advantage in the use of variance component methodology as a linkage or association test is that it makes it possible to combine traits straightforwardly in a vector representation. It is possible to detect whether the putative trait locus is significantly linked or associated with a combined phenotype (i.e., jointly affecting the individual traits). If it is, this might indicate that the biological pathways leading to each of the individual phenotypes overlapped at the function of the gene at this location. This is compared with getting significant results only when analyzing the phenotypes separately, in which case the gene might be involved in each pathway in a different way.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
7 cM from D10S1772 (P = 0.0006) (Table 7). For HDL-C, the association analysis produced a P value of 0.02 with marker D10S1790, located in the close vicinity of D10S546 (0.1 cM, 889 kb) (Table 7). Furthermore, a combined analysis of quantitative HDL-C and TGs produced a P value of 0.008, again with marker D10S546 in the association analysis (Table 7).
Chromosome 20q13.32
The locus on 20q13.32 was originally identified for a binary low HDL-C (LOD 1.9) in the combined analysis of low HDL-C and FCHL families (3) (Tables 46). Here, in the extended study sample of 92 low HDL-C and FCHL families, both the quantitative and qualitative TG traits produced slightly higher two-point LOD scores of 2.2 and 2.4 for this region, with markers D20S173 and D20S102, 11 cM apart. The quantitative and qualitative HDL-C trait also showed some evidence for linkage (Tables 36). When the two study samples were analyzed separately, the 53 FCHL families produced LOD scores of 1.5 (with marker D20S164) and 2.1 (D20S102) for the qualitative HDL-C and TG traits, and the 39 low HDL-C families produced LOD scores of 1.1 (D20S94) and 0.9 (D20S100), respectively. For the quantitative HDL-C and TG traits, similar results were obtained (data not shown). When analyzing the original and new families separately, we observed that although the new families did not show evidence for linkage to this region using the dichotomized approach (Tables 4 6), in the quantitative analyses these new families resulted in a LOD score of 1.4 with the TG trait and in a LOD score of 1.8 with the HDL-C trait (Table 3).
Chromosome 2q31
The locus identified earlier in the Finnish premature CHD families on 2q31 (LOD 3.7) (16) showed evidence for linkage in the 39 low HDL-C families of the present study with the dichotomized low HDL-C trait: a two-point LOD score of 2.7 was obtained with D2S129. The 53 FCHL families produced a LOD score of 1.3 for the low HDL-C trait with marker D2S2324 (Tables 46). The combined study sample of the 92 families resulted in a LOD score of 1.5 with D2S129. Quantitative traits did not provide any significant evidence for linkage to 2q31.
Chromosomes 8q23.1, 16q23.3, and Xq24
The regions on 8q23.1 and 16q23.3 did not provide additional evidence for linkage in the extended study sample of the 92 low HDL-C and FCHL families compared with the original scan results (Tables 46). When we analyzed the original families (3) in their extended form, including all available family members, a LOD score of 3.1 was obtained with the same peak marker, D8S1132, as in the original study (Tables 46), whereas the 38 new families provided no evidence of linkage (Tables 46). Quantitative traits did not provide any significant evidence for linkage (Table 3).
The 16q region, identified previously for low HDL-C in the combined analysis of low HDL-C and FCHL families (a LOD score of 2.2) (3) (Tables 46), provided a LOD score of 1.5 for the dichotomized HDL-C trait and a LOD score of 1.1 for the quantitative HDL-C trait in the extended study sample of 92 low HDL-C and FCHL families (Tables 3 and 4). The region on Xq24, detected originally for the CHD trait in the CHD families (16), resulted in a LOD score of 1.6 for high TGs in the 92 low HDL-C and FCHL families (Tables 46). With quantitative traits, no evidence for linkage to Xq24 was detected.
| DISCUSSION |
|---|
|
|
|---|
We found significant evidence of linkage for quantitative HDL-C to 10q11. In our earlier genome-wide scan, this locus was linked to high TGs (15), a trait highly inversely correlated to HDL-C. The TG trait was also implicated here, especially in the association analysis. Evidence for an association on 10q11 was observed with marker D10S546 (P = 0.0006) using a novel quantitative association program (K. Lange and E. Sobel, unpublished data). This program also allowed us to analyze a combined phenotype of quantitative HDL-C and TGs, which provided some evidence for an association with the same marker (P = 0.008), suggesting that biological pathways leading to each of the individual phenotypes may overlap regarding the function of the underlying gene. Previously, this particular 10q11 region has been linked to obesity and BMI (3840), implying that the underlying gene(s) contributes to various component traits of the metabolic syndrome. The 20q13.32 region was linked to TGs and HDL-C. Several previous studies have reported evidence for linkage between the 20q13 region and type 2 diabetes mellitus, obesity, and HDL-C in different populations (8, 10, 14, 41). However, the TG trait has not, to the best of our knowledge, been linked to 20q13.32 previously. The locus on 2q represents an independent replication: families ascertained for CHD produced the original linkage signal (16), whereas an independent sample of dyslipidemic families was analyzed here. This region has also been linked to type 2 diabetes mellitus (42) and, interestingly, to HDL-C (43, 44). Similarly, in the present study, 2q31 provided evidence for linkage to low HDL-C in the low HDL-C families. This replication thus confirms that the 2q31 region likely harbors a gene for low HDL-C.
The loci on 8q23.1, 16q23.3, and Xq24 did not obtain additional support in the extended study sample. This is not atypical for complex traits. Both FCHL and type 2 diabetes mellitus provide similar examples in which the second study cohort originating from the same population as the original one fails to show evidence for linkage for some of the identified regions (46, 47). It has been suggested that differences in family structures or in the frequency of a disease-causing allele within families could contribute to the power to detect linkage, especially in the case of small study samples (4750). Furthermore, additional families and significant extension of existing families probably also pose a risk of introducing increased genetic heterogeneity within and between the families.
This study does not represent a typical replication attempt with an independent study sample. Rather, we aimed to apply a powerful quantitative approach by extending the genotyping to all families and family members. However, these families were ascertained for multiple dyslipidemic family members, which results in a somewhat restricted variation of the investigated traits and may have reduced the power of QTL analyses. Furthermore, the two inherently different approaches, dichotomized and quantitative, are expected to extract different genetic information from these families and thus can result in different linkage peaks, even in different genomic regions.
Low HDL-C and high TG are highly inversely correlated traits, potentially exhibiting at least a partially shared genetic background. We investigated this relationship by including TGs as a covariate for HDL-C, and vice versa, in the linkage analysis for the most significant signal, 10q11. The LOD score decreased from 3.2 to 2.5, suggesting that the incorporation of TGs as a covariate did not improve the signal-to-noise ratio in our study sample. The 10q11 locus, however, seems to have a pleiotropic effect on the HDL-C and TG levels.
The linked region on 10q11 between markers D10S1233 and D10S546 (9.4 cM, 11.5 Mb) contains three positional candidate genes. First, marker D10S546, which resulted in a P value of 0.0006 in the association analysis, is an intragenic marker of the protocadherin 15 gene (PCDH15), encoding an integral membrane protein that mediates calcium-dependent cell-cell adhesion (51), but is not known to be involved in dyslipidemias or CHD. Second, marker D10S1772, which produced a LOD score of 3.2, is an intragenic marker of gene KIAA1607, a lysosomal trafficking regulator. Although the exact function of this gene is unknown, it may be required for sorting of endosomal resident proteins into late multivesicular endosomes (52). Furthermore, this gene contains an ATP/GTP binding site motif A (P-loop), a domain that also exists in the ABCA1 gene, making KIAA1607 a highly potential candidate gene for HDL-C regulation. The third positional candidate gene is an RNA-binding protein apobec-1 complementation factor that is involved in the editing of the APOB mRNA (53). These three genes provide a starting point to identify the underlying gene(s) involved in the regulation of HDL-C and TG levels.
| ACKNOWLEDGMENTS |
|---|
Manuscript received April 13, 2004 and in revised form July 7, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Mar-Heyming, M. Miyazaki, D. Weissglas-Volkov, N. A. Kolaitis, N. Sadaat, C. Plaisier, P. Pajukanta, R. M. Cantor, T. W.A. de Bruin, J. M. Ntambi, et al. Association of Stearoyl-CoA Desaturase 1 Activity With Familial Combined Hyperlipidemia Arterioscler. Thromb. Vasc. Biol., June 1, 2008; 28(6): 1193 - 1199. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. K. Arnett, A. E. Baird, R. A. Barkley, C. T. Basson, E. Boerwinkle, S. K. Ganesh, D. M. Herrington, Y. Hong, C. Jaquish, D. A. McDermott, et al. Relevance of Genetics and Genomics for Prevention and Treatment of Cardiovascular Disease: A Scientific Statement From the American Heart Association Council on Epidemiology and Prevention, the Stroke Council, and the Functional Genomics and Translational Biology Interdisciplinary Working Group Circulation, June 5, 2007; 115(22): 2878 - 2901. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Wittenburg, M. A. Lyons, R. Li, U. Kurtz, X. Wang, J. Mossner, G. A. Churchill, M. C. Carey, and B. Paigen QTL mapping for genetic determinants of lipoprotein cholesterol levels in combined crosses of inbred mouse strains, J. Lipid Res., August 1, 2006; 47(8): 1780 - 1790. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Weissglas-Volkov, A. Huertas-Vazquez, E. Suviolahti, J. Lee, C. Plaisier, S. Canizales-Quinteros, T. Tusie-Luna, C. Aguilar-Salinas, M.-R. Taskinen, and P. Pajukanta Common Hepatic Nuclear Factor-4{alpha} Variants Are Associated With High Serum Lipid Levels and the Metabolic Syndrome. Diabetes, July 1, 2006; 55(7): 1970 - 1977. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Suviolahti, K. Reue, R. M. Cantor, J. Phan, M. Gentile, J. Naukkarinen, A. Soro-Paavonen, L. Oksanen, J. Kaprio, A. Rissanen, et al. Cross-species analyses implicate Lipin 1 involvement in human glucose metabolism Hum. Mol. Genet., February 1, 2006; 15(3): 377 - 386. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Dastani, L. Quiogue, C. Plaisier, J. C. Engert, M. Marcil, J. Genest, and P. Pajukanta Evidence for a Gene Influencing High-Density Lipoprotein Cholesterol on Chromosome 4q31.21 Arterioscler. Thromb. Vasc. Biol., February 1, 2006; 26(2): 392 - 397. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Yu, D. F. Wyszynski, D. M. Waterworth, S. D. Wilton, P. J. Barter, Y. A. Kesaniemi, R. W. Mahley, R. McPherson, G. Waeber, T. P. Bersot, et al. Multiple QTLs influencing triglyceride and HDL and total cholesterol levels identified in families with atherogenic dyslipidemia J. Lipid Res., October 1, 2005; 46(10): 2202 - 2213. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. L. Pollex and R. A. Hegele Complex Trait Locus Linkage Mapping in Atherosclerosis: Time to Take a Step Back Before Moving Forward? Arterioscler. Thromb. Vasc. Biol., August 1, 2005; 25(8): 1541 - 1544. [Full Text] [PDF] |
||||
![]() |
X. Wang and B. Paigen Genetics of Variation in HDL Cholesterol in Humans and Mice Circ. Res., January 7, 2005; 96(1): 27 - 42. [Abstract] [Full Text] [PDF] |
||||
| ||||||