| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Clinical Research |
From the Franz Volhard Clinic and Max Delbrück Center for Molecular Medicine (H.A.-K., S.B., K.H., A.B., G.N., F.C.L.), Medical Faculty of the Charité, Humboldt University of Berlin, Germany; the Department of Biochemistry (H.A.-K., M.J.), College of Pharmacy, University of Damascus, Syria; Institute of Molecular Genetics and Faculty of Medicine (H.A.-K., E.K.F.B., H.A.D.), University of Heidelberg; Institute for Medical Biometry (K.S.), Informatics and Epidemiology, University of Bonn; and Infogen GmbH (A.B.), subsidiary of Valigen NV, Berlin, Germany.
Correspondence to Friedrich C. Luft, Wiltberg Strasse 50, 13125 Berlin, Germany. E-mail luft{at}fvk-berlin.de
| Abstract |
|---|
|
|
|---|
Key Words: cholesterol LDL hyperlipidemia genetics familial hypercholesterolemia
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
An expanded Materials and Methods section can be found in the online data supplement available at http://www.circresaha.org.
| Results |
|---|
|
|
|---|
=0.00). We genotyped the markers D1S2667, D1S186, and 12 markers in between. These markers cover a wide region that contains the tumor necrosis factor receptor superfamily member 1B (TNFRSF1B). This gene has been implicated in familial-combined hyperlipidemia.12,13 To exclude TNFRSF1B, we sequenced the gene and found no mutations.
|
The recombination event in the individual 5.14 at D1S513 established the proximal boundary of the region as shown by the haplotypes in Figure 1. An ancestral recombination event defined the distal boundary of the region at D1S2826. We then genotyped all the family members for the markers in the region. These marker data enabled us to estimate a more accurate allele frequency, independent of that available in the Genome Database. Table 2 shows the calculated 2-point LOD score in the core family depending on the observed allele frequency. The highest 2-point LOD score was 2.61 at the marker D1S2843 (
=0.00).
|
|
Figure 2 (top panel) shows the calculated multipoint LOD scores for all genotyped markers in the region. We examined whether or not this locus contributes quantitatively (QTL) to cholesterol, LDL, HDL, and triglyceride levels. Heterozygous persons had total cholesterol, LDL, HDL, and triglyceride concentrations that were not significantly different from the 44 individuals with the wild-type haplotypes, consistent with an autosomal recessive mode of inheritance.
|
The genome scan yielded an additional LOD score of 2.63 at the markers D13S1267 and D13S1240. Genotyping more markers enabled us to identify a region of homozygosity on chromosome 13. The boundaries of this region are defined by the recombination events that had occurred in the individual 4.5 at the markers D13S1300 and D13S1266. Haplotypes in the core family are shown in Figure 3. The region thus extends over 13q22-q32, a 14-cM interval. We found this result remarkable in light of our earlier finding that a cholesterol-influencing gene resides at this locus on chromosome 13.8 We therefore genotyped chromosome 13 markers in the entire pedigree. We found one additional homozygous person (5.15) who has the same haplotype as the affected offspring. This individual is also a daughter of the probands father, but stems from a different mother. Table 2 gives the 2-point LOD scores. The multipoint analysis yielded a LOD score of 3.08 shown in Figure 2 (bottom panel). We tested this locus as a QTL in the family; however, the locus did not have a QTL effect.
|
We had observed in our earlier study that the chromosome 13 locus was linked to total and LDL cholesterol in a panel of dizygotic twins and their parents from Berlin.8 We therefore also genotyped this group for the markers D1S552, D1S199, and D1S2843. This analysis yielded linkage of total and LDL cholesterol to the chromosome 1 locus. The multipoint LOD score for total cholesterol was 1.82 (P=0.0019) and 1.86 (P=0.0017) for LDL cholesterol.
In the single-locus analysis, the FH trait appeared to be recessive at both the chromosome 1 and 13 loci. All 3 affected children in the core family (individuals 5.9, 5.10, and 5.11) were homozygous by descent for a few markers at both regions. We regarded this remarkable result as a clear incentive to perform a linkage analysis with 2-trait loci. We used GENHUNTER-TWOLOCUS to test the hypothesis that only a simultaneous homozygous-mutant genotype at the trait loci on both chromosome 1 and 13 leads to the disease ("recessive and recessive"), a multiplicative 2-locus model.14 In addition to the nonparametric analysis with 2 disease loci, we performed a parametric 2-locus analysis under the assumption of this multiplicative model. We also analyzed the data under a heterogeneity model and under an additive model assuming a recessive action of both loci. The models are shown in Table 3. Because 2-locus analysis is computationally much more intensive than an analysis with only 1 disease locus, we were only able to include 4 of the 6 children of the core family. We first chose to include the 3 affected children and 1 unaffected child (5.13) in the GENEHUNTER-TWOLOCUS analysis. The nonparametric 2-locus analysis yielded a maximum nonparametric linkage (NPL) score of 8.64 (P=3.8x10-6). The P value reflects a test of the null hypothesis that both disease loci are unlinked to their corresponding marker maps. The parametric analysis yielded a maximum 2-locus LOD score of 5.41 under the multiplicative model. The heterogeneity and additive models yielded maximum LOD scores of 2.82 and 2.75, respectively. For nonparametric analysis, as well as parametric analysis under all 3 models, the maximum score occurred directly at the markers D1S2843 and D13S1267, the same position for the maximal single-locus LOD score. The results remained the same when we replaced individual 5.13 by another unaffected child; the maximal change in the LOD-score was less than 1%. To compare the results of the parametric 2-trait locus analysis to the single-locus results, we next recalculated the single-locus LOD scores for the pedigree without the individuals 5.12 and 5.14. The maximum LOD scores were 2.80 for chromosome 1 and 2.83 for chromosome 13. Hence, the 2-locus LOD score for the multiplicative model was superior to both single-locus LOD scores.
|
The superiority of the multiplicative model, compared with the other models, could be because of the fact that 3 of the 4 children were affected. To test this possibility, we repeated the 2-trait locus analysis using all 3 unaffected children, but only 1 affected child 5.9. The resulting LOD scores were 3.06 for the multiplicative model, 1.86 for the heterogeneity model, and 1.66 for the additive model. Again, the LOD scores did not change by more than 1% when another affected child was selected for the analysis.
We identified the genomic structure of ARH as shown in Table 4. First, we scanned the bacterial artificial chromosome (BAC) library with primers from the cDNA of ARH and found that the 121-O3 clone harbored the ARH gene. ARH-specific primers, designed based on the ARH cDNA sequence, were used for sequencing the cloned genomic DNA (BAC clone) to determine the positions of the introns and the intron-exon boundaries given in Table 4. Intron-exon boundaries were defined by comparing the genomic sequence with the cDNA sequence. ARH consists of 9 exons. The exon size ranges from 35 (exon 8) to 2066 bp (exon 9, the last exon). Exon 1 contains 21 bp of untranslated sequence followed by the initiation codon ATG. All the exon-intron boundary sequences conform to the consensus splice donor (GT) and acceptor (AG) sites, with the exception of the splice donor site of intron 8, which is GG instead of GT.
|
We sequenced the ARH gene in the core family and found a mutation that lies in the acceptor site of intron 1. The mutation converts the original acceptor-splice site from G
C, a transversion mutation. Figure 4A shows the intron 1/exon 2 junction in a normal proband (top) and in an affected child (bottom). The affected offspring are all homozygous for the mutation, whereas the parents are heterozygous. Figure 4B shows the mutated nucleotide sequence of the coding and noncoding strands at the intron 1/exon 2 junction. We found a restriction enzyme PvuII to cut the PCR product of exon 2 at the mutated donor-acceptor site. We used this enzyme to screen the entire family by using the exon 2 PCR product. On recognition, the enzyme cuts the 316-bp PCR product of exon 2 producing 2 fragments, 256 and 60 bp, in the homozygous state and 3 fragments, 316, 256, and 60 bp, in the heterozygous state. Figure 4C shows the status of the core family members with regard to the ARH G transversion mutation. The parents (4.5 and 3.3) are heterozygous for the mutation. The 3 affected children (5.9, 5.10, and 5.11) are homozygous. In agreement with the reported haplotypes, 2 of unaffected children (5.13 and 5.14) are carriers and 1 (5.12) is a noncarrier. We screened the entire family and found 21 individuals who were heterozygous for the mutation. These individuals are the same who share the haplotype with the affected patients for the microsatellite markers genotyped in the region. Only the 3 affected offspring in the core family are homozygous for the mutation. We also screened 31 normal Syrian individuals for the mutation and another 30 normal unrelated German individuals. None had the mutation.
|
We also found 3 single nucleotide polymorphisms (SNPs) in ARH gene. One variant lies in exon 6 at bp 625 of ARH cDNA converting C
T and leads to an amino acid substitution from Pro to Ser. The second SNP lies in the intron 6 at bp 28 converting A
G. The third SNP lies in exon 7 at bp 675 of ARH cDNA converting A
G but without any change in the coding sequence. To explore a possible effect of the variant in exon 6 on the phenotype, we genotyped the whole family for these 3 SNPs. Five persons were homozygous for the variant in exon 6. We also genotyped our twin panel for the 3 SNPs. The data showed no significant association between these variants and LDL levels (P<0.1).
| Discussion |
|---|
|
|
|---|
In addition to the chromosome 1p36.1-p35 locus, we also found compelling evidence for linkage to 13q22-q32. We do not believe that this finding is related to chance. In our earlier study, we identified this locus as the site for a gene that lowers LDL.8 We studied a family with autosomal-dominant familial hypercholesterolemia that had affected heterozygous and homozygous members with lower than expected LDL values. We corroborated this locus with a linkage analysis in dizygotic twins and their parents. Eden et al6 did not report any evidence of linkage to 13q22-q32 in their study. However, the family in our present study and the family we reported earlier are both unrelated Arab families8; their genetic background may be similar.
Garcia et al7 cloned the gene responsible for ARH residing at chromosome 1p35. Their exciting finding was a new protein that is responsible for LDL receptor malfunction. The gene they identified codes for a protein that has a phosphotyrosine binding domain PTD. This domain was found to interact with the cytoplasmic tail of the LDL receptor.15 Malfunction of this interaction could account for a disabled LDL receptor. ARH is expressed in many tissues, even in those that have only a very low LDL receptor expression. Conceivably ARH has other functions; however, similar to the families described by Garcia et al,7 our affected family members have no other discernible phenotypes.
As reported in the present study, we determined the exon-intron structure of ARH. All splice sites conform to the GT/AG rule, except for intron 8, where the donor splice site GT was converted to GG. However, there are a few exceptions to this rule.1620 We identified the intron 1 acceptor-site mutation in the affected persons of the core family. The mutation replaces the G bp of the acceptor-splice site (AG) by C. We corroborated our findings by analysis of the other family members. We also examined 31 healthy Syrians and 30 Germans. All were wild-type for the mutation. A defect in intron 1 splice-donor site will either skip exon 2 on splicing or will retain intron 1 in the processed mRNA. Both cases contain stop codons that would result in a truncated protein missing the amino acids coded by exons 2 to 9. We also found 3 SNPs in exon 6, intron 6, and exon 7. The SNP in exon 6 causes an amino acid exchange Pro
Ser. We excluded this variant when we found other homozygous family members who had normal total cholesterol and LDL levels. Moreover, we examined these variants in our twin panel and found that they were not associated with total cholesterol and LDL levels.
The parametric and nonparametric 2-locus analysis was not only done to verify that ARH is indeed governed by 2 loci, but also to gain insight into how the 2 genes may act together. Due to the restrictions on the pedigree size with GENEHUNTER-TWOLOCUS, we were not able to include all 6 children in the analysis simultaneously. We performed the analysis with all 3 affected children and 1 unaffected child. A recessive-and-recessive multiplicative model gave the strongest results, irrespective of the unaffected child selected. The model reflects that only individuals with a homozygous mutant genotype at both loci will develop the trait. Here, the 2-locus LOD score was 5.41, which was almost twice as high as the result obtained for the heterogeneity and additive models. With the latter models, the disease may be caused by a homozygous-mutant genotype at only 1 locus.
It may be argued that the multiplicative 2-locus model gave the strongest results merely because 2 unaffected individuals were not included in the analysis. We therefore repeated the 2-locus calculations with all 3 unaffected children and 1 affected child. The resulting LOD score was 3.06 for the multiplicative model, 1.86 for the heterogeneity model, and 1.66 for the additive model. Hence, even in this setting the superiority of the multiplicative model remained over the other 2 models. These results support a multiplicative action of both loci, where a homozygous mutant genotype at both chromosome 1 and 13 loci is necessary to express the disease.
Our 2-locus multiplicative model suggests that a mutation at only 1 locus is not sufficient to develop the disease. However, our results suggest that the chromosome 1 mutation has a dramatic effect on the adaptor protein. Our analysis raises the possibility that 2 proteins coded by 2 separate genes at distant loci are able to perform the same biological function. A mutation at 1 locus may disable one of the proteins, whereas the second is still available for normal function. Such a model would require disruption of both genes, disabling both proteins to express the phenotype. If our hypothesis is correct, a mutation on chromosome 1 leading to a dramatically altered protein is still in accord with the multiplicative 2-locus trait model. Our linkage and interaction analysis strongly supports the hypothesis that 2 loci are involved in this form of recessive FH in our family.
We were interested to observe linkage for the 1p36.1-p35 locus to total and LDL cholesterol in normal dizygotic twins and their parents. We verified earlier findings that LDL, HDL, total cholesterol, triglycerides, and body mass index are all strongly influenced by genetic variance.21 In the Syrian family, the chromosome 1 locus did not appear as a QTL. The effect here is autosomal-recessive and qualitative in nature. In the twins, the locus appeared as a QTL, suggesting that other variations in the responsible gene are operative. We believe that these studies are unique because we not only mapped a digenic syndrome, but also suggested its relevance in a completely unrelated, normal population. Identification of ARH and its mutations are clearly important to our understanding of cholesterol metabolism and cardiovascular risk. Elucidating the gene at the chromosome 13q locus and the basis for the interaction will have a high priority.
| Acknowledgments |
|---|
Received October 29, 2001; revision received April 1, 2002; accepted April 1, 2002.
| References |
|---|
|
|
|---|
2. Myant NB. Familial defective apolipoprotein B-100: a review, including some comparisons with familial hypercholesterolaemia. Atherosclerosis. 1993; 104: 118.[CrossRef][Medline] [Order article via Infotrieve]
3. Varret M, Rabes JP, Saint-Jore B, Cenarro A, Marinoni JC, Civeira F, Devillers M, Krempf M, Coulon M, Thiart R, Kotze MJ, Schmidt H, Buzzi JC, Kostner GM, Bertolini S, Pocovi M, Rosa A, Farnier M, Martinez M, Junien C, Boileau C. A third major locus for autosomal dominant hypercholesterolemia maps to 1p34.1-p32. Am J Hum Genet. 1999; 64: 13781387.[CrossRef][Medline] [Order article via Infotrieve]
4.
Hunt SC, Hopkins PN, Bulka K, McDermott MT, Thorne TL, Wardell BB, Bowen BR, Ballinger DG, Skolnick MH, Samuels ME. Genetic localization to chromosome 1p32 of the third locus for familial hypercholesterolemia in a Utah kindred. Arterioscler Thromb Vasc Biol. 2000; 20: 10891093.
5. Ciccarese M, Pacifico A, Tonolo G, Pintus P, Nikoshkov A, Zuliani G, Fellin R, Luthman H, Maioli M. A new locus for autosomal recessive hypercholesterolemia maps to human chromosome 15q25-q26. Am J Hum Genet. 2000; 66: 453460.[CrossRef][Medline] [Order article via Infotrieve]
6. Eden ER, Naoumova RP, Burden JJ, McCarthy MI, Soutar AK. Use of homozygosity mapping to identify a region on chromosome 1 bearing a defective gene that causes autosomal recessive homozygous hypercholesterolemia in two unrelated families. Am J Hum Genet. 2001; 68: 653660.[CrossRef][Medline] [Order article via Infotrieve]
7.
Garcia CK, Wilund K, Arca M, Zuliani G, Fellin R, Maioli M, Calandra S, Bertolini S, Cossu F, Grishin N, Barnes R, Cohen JC, Hobbs HH. Autosomal recessive hypercholesterolemia caused by mutations in a putative LDL receptor adaptor protein. Science. 2001; 292: 13941398.
8. Knoblauch H, Müller-Myhsok B, Busjahn A, Ben Avi L, Bähring S, Baron H, Heath SC, Regina Uhlmannm, Faulhaber H-D, Shpitzen S, Aydin A, Reshef A, Rosenthal M, Eliav O, Mühl A, Lowe A, Schurr D, Harats D, Jeschke E, Friedlander Y, Schuster H, Luft FC, Leitersdorf E. A cholesterol-lowering gene exists on chromosome 13q. Am J Hum Genet. 2000; 66: 157166.[CrossRef][Medline] [Order article via Infotrieve]
9. Nevo S. Genetic blood markers in the Arab Druze of Israel. Am J Phys Anthropol. 1988; 77: 183190.[CrossRef][Medline] [Order article via Infotrieve]
10. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972; 18: 499502.[Abstract]
11. Strauch K, Fimmers R, Kurz T, Deichmann KA, Wienker TF, Baur MP. Parametric and nonparametric multipoint linkage analysis with imprinting and two-locus-trait models: application to mite sensitization. Am J Hum Genet. 2000; 66: 19451957.[CrossRef][Medline] [Order article via Infotrieve]
12.
Glenn CL, Wang WY, Benjafield AV, Morris BJ. Linkage and association of tumor necrosis factor receptor 2 locus with hypertension, hypercholesterolemia and plasma shed receptor. Hum Mol Genet. 2000; 9: 19431949.
13.
Geurts JM, Janssen RG, van Greevenbroek MM, van Der Kallen CJ, Cantor RM, Bu Xd, Aouizerat BE, Allayee H, Rotter JI, de Bruin TW. Identification of TNFRSF1B as a novel modifier gene in familial combined hyperlipidemia. Hum Mol Genet. 2000; 9: 20672074.
14. Risch N. Linkage strategies for genetically complex traits, I: multilocus models. Am J Hum Genet. 1990; 46: 222228.[Medline] [Order article via Infotrieve]
15.
Trommsdorff M, Borg JP, Margolis B, Herz J. Interaction of cytosolic adaptor proteins with neuronal apolipoprotein E receptors and the amyloid precursor protein. J Biol Chem. 1998; 273: 3355633560.
16. Green MR. Biochemical mechanisms of constitutive and regulated pre-mRNA splicing. Annu Rev Cell Biol. 1991; 7: 559599.[CrossRef]
17. Dong J, Kukula AK, Toyoshima M, Nakajima M. Genomic organization and chromosome localization of the newly identified human heparanase gene. Gene. 2000; 253: 171178.[CrossRef][Medline] [Order article via Infotrieve]
18.
Broderick TP, Schaff DA, Bertino AM, Dush MK, Tischfield JA, Stambrook PJ. Comparative anatomy of the human APRT gene and enzyme: nucleotide sequence divergence and conservation of a nonrandom CpG dinucleotide arrangement. Proc Natl Acad Sci U S A. 1987; 84: 33493353.
19.
Schroeder HW Jr, Walter MA, Hofker MH, Ebens A, Willems van Dijk K, Liao LC, Cox DW, Milner EC, Perlmutter RM. Physical linkage of a human immunoglobulin heavy chain variable region gene segment to diversity and joining region elements. Proc Natl Acad Sci U S A. 1988; 85: 81968200.
20.
Alexandre D, Chuchana P, Brockly F, Blancher A, Lefranc G, Lefranc MP. First genomic sequence of a human Ig variable lambda gene belonging to subgroup I: functional genes, pseudogenes and vestigial sequences are interspersed in the IGLV locus. Nucleic Acids Res. 1989; 17: 3975.
21.
Knoblauch H, Busjahn A, Münter S, Nagy Z, Faulhaber H-D, Schuster H, Luft FC. Heritability analysis of lipids and three gene loci in twins link the macrophage scavenger receptor to high-density lipoprotein cholesterol concentrations. Arterioscler Thromb Vasc Biol. 1997; 17: 20542060.
This article has been cited by other articles:
![]() |
P. L.M. Dahia, K. Hao, J. Rogus, C. Colin, M. A.G. Pujana, K. Ross, D. Magoffin, N. Aronin, A. Cascon, C. Y. Hayashida, et al. Novel Pheochromocytoma Susceptibility Loci Identified by Integrative Genomics Cancer Res., November 1, 2005; 65(21): 9651 - 9658. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Bosse, Y. C. Chagnon, J.-P. Despres, T. Rice, D. C. Rao, C. Bouchard, L. Perusse, and M.-C. Vohl Compendium of genome-wide scans of lipid-related phenotypes: adding a new genome-wide search of apolipoprotein levels J. Lipid Res., December 1, 2004; 45(12): 2174 - 2184. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Soutar, R. P. Naoumova, and L. M. Traub Genetics, Clinical Phenotype, and Molecular Cell Biology of Autosomal Recessive Hypercholesterolemia Arterioscler Thromb Vasc Biol, November 1, 2003; 23(11): 1963 - 1970. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Al-Kateb, E. K.F. Bautz, F. C. Luft, and S. Bahring A Splice Mutation in a Syrian Autosomal Recessive Hypercholesterolemia Family Causes a Two-Nucleotide Deletion of mRNA Circ. Res., September 5, 2003; 93 (5): e49 - e50. [Full Text] [PDF] |
||||
![]() |
M. Harada-Shiba, A. Takagi, Y. Miyamoto, M. Tsushima, Y. Ikeda, S. Yokoyama, and A. Yamamoto Clinical Features and Genetic Analysis of Autosomal Recessive Hypercholesterolemia J. Clin. Endocrinol. Metab., June 1, 2003; 88(6): 2541 - 2547. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Mishra, S. C. Watkins, and L. M. Traub The autosomal recessive hypercholesterolemia (ARH) protein interfaces directly with the clathrin-coat machinery PNAS, December 10, 2002; 99(25): 16099 - 16104. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. R. Wilund, M. Yi, F. Campagna, M. Arca, G. Zuliani, R. Fellin, Y.-K. Ho, J. V. Garcia, H. H. Hobbs, and J. C. Cohen Molecular mechanisms of autosomal recessive hypercholesterolemia Hum. Mol. Genet., November 15, 2002; 11(24): 3019 - 3030. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.H. HOBBS, G.A. GRAF, L. YU, K.R. WILUND, and J.C. COHEN Genetic Defenses against Hypercholesterolemia Cold Spring Harb Symp Quant Biol, January 1, 2002; 67(0): 499 - 506. [Abstract] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Research Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2002 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |