Genome-wide Scan for Type 2 Diabetes Loci in Hong Kong Chinese and Confirmation of a Susceptibility Locus on Chromosome 1q21-q25

Maggie C.Y. Ng1,2, Wing-Yee So1, Nancy J. Cox3,4, Vincent K.L. Lam1, Clive S. Cockram1, Julian A.J.H. Critchley1,{dagger}, Graeme I. Bell2,3,4, and Juliana C.N. Chan1

1 Department of Medicine and Therapeutics, The Chinese University of Hong Kong, Prince of Wales Hospital, Shatin, Hong Kong
2 Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois
3 Department of Human Genetics, University of Chicago, Chicago, Illinois
4 Department of Medicine, University of Chicago, Chicago, Illinois


    ABSTRACT
 TOP
 ABSTRACT
 RESEARCH DESIGN AND METHODS
 REFERENCES
 
We conducted an autosomal genome scan to map loci for type 2 diabetes in a Hong Kong Chinese population. We studied 64 families, segregating type 2 diabetes, of which 57 had at least one member with an age at diagnosis of <=40 years. These families included a total of 126 affected sibpairs and 4 other affected relative pairs. Nonparametric linkage analysis revealed seven regions showing nominal evidence for linkage with type 2 diabetes (logarithm of odds [LOD] >0.59, Ppointwise < 0.05): chromosome 1 at 173.9 cM (LOD = 3.09), chromosome 3 at 26.3 cM (LOD = 1.27), chromosome 4 at 135.3 cM (LOD = 2.63), chromosome 5 at 139.3 cM (LOD = 0.84), chromosome 6 at 178.9 cM (LOD = 1.91), chromosome 12 at 48.7 cM (LOD = 1.99), and chromosome 18 at 28.1 cM (LOD = 1.00). Simulation studies showed genome-wide significant evidence for linkage of the chromosome 1 region (Pgenome-wide = 0.036). We have confirmed the results of previous studies for the presence of a susceptibility locus on chromosome 1q21-q25 (173.9 cM) and suggest the locations of other loci that may contribute to the development of type 2 diabetes in Hong Kong Chinese.

Type 2 diabetes is a heterogeneous disease characterized by insulin resistance and pancreatic ß-cell dysfunction (1). Genetic factors play an important role in the development of type 2 diabetes. Despite considerable effort, there has been relatively little progress in identifying the genes that affect risk. This may be due, at least in part, to phenotypic heterogeneity; i.e., type 2 diabetes is not one disease but many, each characterized by hyperglycemia. Genome scans to map type 2 diabetes susceptibility loci have been conducted in many different populations (24). Some of the mapped loci have been observed in multiple populations, including those on chromosomes 1q21-q24, 12q, and 20. Other regions, however, may be unique to specific populations, e.g., Mexican Americans and chromosome 2q37.3 and the Amish and chromosome 14q11. It is unclear if this reflects underlying phenotypic heterogeneity, racial/ethnic differences in susceptibility allele frequencies, or differences in sample size, study design, and analytical methods.

The prevalence of diabetes in Hong Kong is among the highest in Asia, with an age-adjusted rate of 8.9% in 1995 (5). This high prevalence is thought to be due to the affluent and westernized lifestyle of Hong Kong Chinese, resulting from rapid socioeconomic development over the last few decades. This is compatible with the thrifty genotype hypothesis, in which favorable genetic traits that facilitate energy storage and mobilization in times of famine are detrimental during time of plenty, leading to the development of modern diseases such as diabetes (6).

We carried out an autosomal genome scan in 64 families with type 2 diabetes (126 affected sibpairs and 4 other affected relative pairs). We also considered glucose intolerance (GIT) as a trait (defined as type 2 diabetes, impaired fasting glucose, or impaired glucose tolerance) because impaired fasting glucose and impaired glucose tolerance may be precursors to type 2 diabetes. We carried out an autosomal scan for GIT in 102 families (306 affected sibpairs and 46 affected relative pairs). The clinical characteristics of affected subjects involved in linkage analyses were similar except for the higher BMI and fasting glucose levels in the type 2 diabetes group (Table A1 in the online appendix [available at http://diabetes.diabetesjournals.org]).

Seven regions showed nominal evidence for linkage with type 2 diabetes (logarithm of odds [LOD] >0.59, Ppointwise < 0.05) (Table 1, Fig. 1, and online appendix Table A2). The region showing the strongest evidence was on chromosome 1q21-q25 at 173.9 cM (LOD = 3.09, Ppointwise = 0.00008, 1-LOD CI = 164–180 cM). The LOD scores in these seven regions were not increased when we used GIT as the trait.


View this table:
[in this window]
[in a new window]
 
TABLE 1 Regions showing nominal evidence of linkage to type 2 diabetes in Hong Kong Chinese*

 


View larger version (36K):
[in this window]
[in a new window]
 
FIG. 1. Multipoint linkage analyses for diabetes and glucose intolerance. The horizontal axis is centimorgans from the p-terminus.

 
We used simulation to assess the significance of the results and found the LOD score threshold of 2.93 for >=1 independent linked regions (ILRs) for type 2 diabetes, yielding an estimated Pgenome-wide < 0.05. The chromosome 1q region demonstrated genome-wide significant linkage for type 2 diabetes (LOD = 3.09, Pgenome-wide = 0.036). Chromosomes 4q (LOD = 2.63), 6q (LOD = 1.91), and 12p (LOD = 1.99) showed suggestive evidence for linkage based on observation in our simulations of an average of one ILR per genome scan (7) at a LOD score of 1.59 (Table 1).

We also used locus-counting methods to evaluate the linkage results (8). We observed seven, five, and two loci with LOD >0.59 (Ppointwise < 0.05), >1.18 (Ppointwise < 0.01), and >2.08 (Ppointwise < 0.001), respectively. The simulations showed that 392 of the 500 replicates (78%) had seven or more loci at LOD >0.59, whereas 56 (11%) had five or more loci at LOD >1.18, and 21 (4%) had two or more loci at LOD >2.08. Overall, the results suggest that there may be one or more susceptibility loci in the other six regions with nominal evidence for linkage.

Could our results be due to the presence of genotyping error or missing data? In general, these problems are believed to reduce the power to detect linkage (9). Our genotyping error rate as assessed by blind duplicates is comparable with that of previous genome screens for type 2 diabetes (10). LOD scores obtained before and after removal of Merlin-derived potential genotyping errors were similar in the seven regions with nominal evidence for linkage (chromosome 1, 3.07 vs. 3.09; chromosome 3, 1.27 vs. 1.27; chromosome 4, 2.42 vs. 2.63; chromosome 5, 0.88 vs. 0.84; chromosome 6, 2.08 vs. 1.91; chromosome 12, 1.99 vs. 1.99; and chromosome 18, 1.00 vs. 1.00). Our missing data rate is comparable with that observed in another type 2 diabetes genome scan (11), where it was demonstrated that 15% missing data did not substantially decrease the power to detect linkage. It should be noted that our moderately sized pedigrees allow recovery of some missing information through haplotype reconstruction. Finally, our study assessed significance through simulation, and thus the consequences of the observed missing data patterns are appropriately taken into account.

In this study, we used both strict (type 2 diabetes) and loose (GIT) criteria of glucose homeostasis to define affection status. The former should yield a more homogeneous sample with respect to disease, whereas the latter provides a larger sample size. Our results generally indicate overlapping linkage signals for type 2 diabetes and GIT (chromosomes 1, 3–5, and 12) but the linkage signals for type 2 diabetes tended to be stronger than those for GIT (chromosomes 1, 4, and 12) (Fig. 1), despite the fact that the sample size for type 2 diabetes was smaller. The higher signals for type 2 diabetes suggest that the reduced heterogeneity is important for detecting the larger signals.

The region on chromosome 1q (173.9 cM from p-terminus) has also been linked to type 2 diabetes and related traits in Chinese from Shanghai, Pima Indians, and Europeans (24,1014) (Table 2). There are >100 known genes in this region. Results from association studies on genes encoding for potassium inwardly rectifying channel J9 (KCNJ9) (15), C-reactive protein (CRP) (16), and phospholipase A2 4A (PLA2G4A) (17) suggest there are associations that require confirmation. This region is the focus of an international collaborative project to positionally clone the gene(s) that affect type 2 diabetes risk. At least one or more of the six other regions showing nominal evidence for type 2 diabetes may harbor true susceptibility loci based on the results of locus-counting methods (8) and linkage to these regions in other studies (Table 2). The identification of the genes and alleles that affect risk to type 2 diabetes may lead to a better understanding of the genetic basis of type 2 diabetes in Hong Kong Chinese and other populations.


View this table:
[in this window]
[in a new window]
 
TABLE 2 Regions showing nominal evidence for linkage with type 2 diabetes in Hong Kong Chinese and their replications in other studies

 

    RESEARCH DESIGN AND METHODS
 TOP
 ABSTRACT
 RESEARCH DESIGN AND METHODS
 REFERENCES
 
Since 1995, all patients attending the diabetes clinic of the Prince of Wales Hospital have been documented for detailed phenotypes and family history (parents, siblings, and offspring), which formed the Prince of Wales Hospital Diabetes Registry (18). Through this registry, Hong Kong Family Diabetes Study ascertained families through a diabetic proband with available first-degree relatives for screening at two stages. At the first stage, 161 families were recruited through a diabetic proband with age at diagnosis <=40 years (early-onset group). At the second stage, an additional 30 families were recruited through a diabetic proband with age at diagnosis >40 years and with positive family history of diabetes (late-onset group). The second cohort aimed to increase the sample size for gene-mapping study. A total of 191 families were recruited. Among the young-onset families, 100 of the families were screened for antibodies to GAD. Individuals with clinical type 1 diabetes (diabetic ketoacidosis or heavy ketonuria [>3] or required continuous insulin treatment within 1 year of diagnosis) or with antibodies to GAD were excluded. Subsequently, families without any GIT or diabetic members remaining were excluded. In addition, 48 of the probands with positive family history had participated in a previous study to screen for mutations in the MODY2 and MODY3 genes and in the mitochondrial DNA nucleotide 3243 A->G mutation (19). Families with these gene mutations were also excluded. A total of 179 families remained. From these 179 families, studies on the dichotomous traits of type 2 diabetes and GIT included 64 and 102 families, respectively, due to the requirement for at least two affected (nonparent offspring) relatives. The probands, all of their first- and second-degree relatives, and the probands’ spouses were recruited between January 1998 and March 2002. The clinical characteristics, family structure, and distribution of the affected relative pairs of these families are summarized in Tables A1, A3, and A4, respectively, in the online appendix. Of special note, 57 of the 64 families used in the type 2 diabetes studies and 86 of the 102 in the GIT studies had at least one affected member diagnosed before age 40 years. Informed consent was obtained for each participating subject. This study was approved by the clinical research ethics committee of the Chinese University of Hong Kong.

Clinical studies.
All available family members were assessed for blood pressure and standard anthropometric parameters and completed a questionnaire on demographic data, family and medical histories, and lifestyle. Fasting blood samples were collected for measurement of plasma glucose, insulin, C-peptide, lipid profile, liver and renal function, and DNA extraction. Family members with no known history of type 2 diabetes were tested with a 75-g oral glucose tolerance test. Type 2 diabetes, impaired glucose tolerance, and impaired fasting glucose were diagnosed using the 1998 World Health Organization criteria (20). Plasma glucose was measured by a glucose oxidase method (Diagnostic Chemicals, Charlottetown, PEI, Canada). Plasma insulin and C-peptide were measured by enzyme-linked immunosorbent assay (Dako Diagnostics, Glostrup, Denmark).

Genotyping.
A total of 355 autosomal microsatellite markers (Research Genetics, Huntsville, AL) were typed. These included 320 markers from the Human Screening Set, version 10, and an additional 35 markers that replaced 58 discarded markers from the screening set due to lack of amplification, monomorphism, or significant departure from Hardy-Weinberg equilibrium (P < 0.05). The average intermarker distance was 10 cM, and average heterozygosity was 71%. The largest gap was 27.3 cM between D18S1357 and D18S1371 due to removal of two markers not in Hardy-Weinberg equilibrium. Microsatellite markers were amplified with fluorescent-labeled primers using a universal multiplex PCR protocol. Genotypes were performed by capillary electrophoresis (ABI 3100 DNA analyzer). Alleles were called and scored using the Genescan 3.51 and Genotyper 3.6 software (ABI) followed by manual checking. Sixteen percent of genotype data were missing. Blind duplicates were included, with an allele-wise genotyping error rate of 1.2%. Genetic relationships among family members were checked by Relpair (21) and Prest (22). After correction of family relationships and removal of unrelated individuals, Mendelian incompatibilities for all markers were identified by PedCheck (version 1.1) (23). The genotypes that were most likely to have errors within a family were removed, if possible. Otherwise, the genotypes of the incompatible marker were removed in all subjects within the family. Potential genotyping errors observed in 0.17% of the data using Merlin (version 0.9.12b) were removed (24).

Linkage analyses.
Marker allele frequencies were estimated from all subjects in 179 families, including those involved in type 2 diabetes and GIT studies. Marker position was based on the sex-average maps from the Marshfield Medical Research Foundation (http://research.marshfieldclinic.org/genetics/Map_Markers/maps/IndexMapFrames.html). Nonparametric multipoint linkage analyzes were performed with Merlin (version 0.9.12b) using the score-pairs statistic (24). LOD scores with associated pointwise P values are reported (25). Two linked regions were defined as independent if their maximum LOD score positions differed by >=40 cM. Simulation studies were conducted to assess the significance of the observed linkage signals. Five hundred simulated genotype datasets were generated by gene dropping with Merlin under the null hypothesis of no linkage. The maps, marker allele frequencies, pedigree structures, and missing genotype patterns used in these simulations were identical to those in the actual data. The linkage analyzes were then conducted on the simulated datasets. In each simulation, the number of ILRs greater than a particular LOD score threshold was tallied. Pgenome-wide for n ILRs at a LOD score threshold was calculated as the proportion of simulations with >=n ILRs at that particular LOD score. A significant excess of ILRs in the observed dataset was defined as Pgenome-wide < 0.05.


    ACKNOWLEDGMENTS
 
This work was supported by grants from the Hong Kong Research Grants Committee, the Chinese University of Hong Kong Strategic Grant Program, the Innovation and Technology Support Fund, Hong Kong (ITS/33/00), and U.S. Public Health Service Grants DK-20595, DK-47486, and DK-55889. G.I.B. is an Investigator of the Howard Hughes Medical Institute.

We thank Emily Poon and Patty Tse for their technical support and Dr. June Li and Cherry Chiu for their help with the recruitment of patients and their family members. We also thank William Wen for computational support and Cheryl Roe for simulation data analysis. We are also indebted to the late Professor Robert Turner, Oxford University, for his collaboration in our family study. We thank all of the nursing and medical staff at the Prince of Wales Hospital Diabetes and Endocrine Centre for their dedication and professionalism.

This work is dedicated to the memory of Professor Julian Critchley.


    FOOTNOTES
 
{dagger} J.A.J.H.C. passed away in 2001. Back

Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.

Address correspondence and reprint requests to Dr. Maggie C.Y. Ng, Howard Hughes Medical Institute, University of Chicago, 5841 S. Maryland Ave., MC1028, Chicago, IL 60637. E-mail: maggieng{at}uchicago.edu

Received for publication March 5, 2004 and accepted in revised form March 8, 2004

GIT, glucose intolerance; ILR, independent linked region; LOD, logarithm of odds


    REFERENCES
 TOP
 ABSTRACT
 RESEARCH DESIGN AND METHODS
 REFERENCES
 

  1. Zimmet P, Alberti KG, Shaw J: Global and societal implications of the diabetes epidemic. Nature414 :782 –787,2001[Medline]
  2. Gloyn AL, McCarthy MI: The genetics of type 2 diabetes. Best Pract Res Clin Endocrinol Metab15 :293 –308,2001[Medline]
  3. Florez JC, Hirschhorn J, Altshuler D: The inherited basis of diabetes mellitus: implications for the genetic analysis of complex traits. Annu Rev Genomics Hum Genet4 :257 –291,2003
  4. McCarthy MI: Growing evidence for diabetes susceptibility genes from genome scan data. Curr Diab Rep3 :159 –167,2003[Medline]
  5. Cockram CS, Chan JCN: The epidemiology of diabetes in the western pacific region (excluding Japan). In Diabetes in the New Millennium. Turtle JR, Ed. Sydney, Australia, Pot Still Press,1999 , p.11 –22
  6. Neel JV, Weder AB, Julius S: Type II diabetes, essential hypertension, and obesity as "syndromes of impaired genetic homeostasis": the "thrifty genotype" hypothesis enters the 21st century. Perspect Biol Med42 :44 –74,1998[Medline]
  7. Lander ES, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet11 :241 –247,1995[Medline]
  8. Wiltshire S, Cardon LR, McCarthy MI: Evaluating the results of genomewide linkage scans of complex traits by locus counting. Am J Hum Genet71 :1175 –1182,2002[Medline]
  9. Sullivan PF, Neale BM, Neale MC, van den Oord E, Kendler KS: Multipoint and single point non-parametric linkage analysis with imperfect data. Am J Med Genet121B :89 –94,2003[Medline]
  10. Hanson RL, Ehm MG, Pettitt DJ, Prochazka M, Thompson DB, Timberlake D, Foroud T, Kobes S, Baier L, Burns DK, Almasy L, Blangero J, Garvey WT, Bennett PH, Knowler WC: An autosomal genomic scan for loci linked to type II diabetes mellitus and body-mass index in Pima Indians. Am J Hum Genet63 :1130 –1138,1998[Medline]
  11. Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes U.K. Warren 2 Repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet69 :553 –569,2001[Medline]
  12. Xiang K, Wang Y, Zheng T, Jia W, Li J, Chen L, Shen K, Wu S, Lin X, Zhang G, Wang C, Wang S, Lu H, Fang Q, Shi Y, Zhang R, Xu J, Weng Q: Genome-wide search for type 2 diabetes/impaired glucose homeostasis susceptibility genes in the Chinese: significant linkage to chromosome 6q21-q23 and chromosome 1q21-q24. Diabetes53 :228 –234,2004[Abstract/Free Full Text]
  13. Hsueh WC, St. Jean PL, Mitchell BD, Pollin TI, Knowler WC, Ehm MG, Bell CJ, Sakul H, Wagner MJ, Burns DK, Shuldiner AR: Genome-wide and fine-mapping linkage studies of type 2 diabetes and glucose traits in the Old Order Amish: evidence for a new diabetes locus on chromosome 14q11 and confirmation of a locus on chromosome 1q21-q24. Diabetes52 :550 –557,2003[Abstract/Free Full Text]
  14. Vionnet N, Hani EH, Dupont S, Gallina S, Francke S, Dotte S, De Matos F, Durand E, Lepretre F, Lecoeur C, Gallina P, Zekiri L, Dina C, Froguel P: Genomewide search for type 2 diabetes-susceptibility genes in French whites: evidence for a novel susceptibility locus for early-onset diabetes on chromosome 3q27-qter and independent replication of a type 2-diabetes locus on chromosome 1q21-q24. Am J Hum Genet67 :1470 –1480,2000[Medline]
  15. Wolford JK, Hanson RL, Kobes S, Bogardus C, Prochazka M: Analysis of linkage disequilibrium between polymorphisms in the KCNJ9 gene with type 2 diabetes mellitus in Pima Indians. Mol Genet Metab73 :97 –103,2001[Medline]
  16. Wolford JK, Gruber JD, Ossowski VM, Vozarova B, Antonio Tataranni P, Bogardus C, Hanson RL: A C-reactive protein promoter polymorphism is associated with type 2 diabetes mellitus in Pima Indians. Mol Genet Metab78 :136 –144,2003[Medline]
  17. Wolford JK, Konheim YL, Colligan PB, Bogardus C: Association of a F479L variant in the cytosolic phospholipase A2 gene (PLA2G4A) with decreased glucose turnover and oxidation rates in Pima Indians. Mol Genet Metab79 :61 –66,2003[Medline]
  18. Chan JCN, Lau M, Wong R, Chow CC, Yeung V, Loo KM, Mong M, Yeung T, Ko GTC, Li KY, So WY, Chan WB, Lu K, Cockram CS: Delivery of diabetes care: the experience at the Prince of Wales Hospital. Hospital Authority Quality Bulletin2 :3 –21,1997
  19. Ng MCY, Lee SC, Ko GTC, Li JKY, So WY, Hashim Y, Barnett AH, Mackay IR, Critchley JAJH, Cockram CS, Chan JCN: Familial early-onset type 2 diabetes in Chinese patients: obesity and genetics have more significant roles than autoimmunity. Diabetes Care24 :663 –671,2001[Abstract/Free Full Text]
  20. Alberti KGMM, Zimmet PZ: Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus, provisional report of a WHO consultation. Diabet Med15 :539 –553,1998[Medline]
  21. Boehnke M, Cox NJ: Accurate inference of relationships in sib-pair linkage studies. Am J Hum Genet61 :423 –429,1997[Medline]
  22. McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet66 :1076 –1094,2000[Medline]
  23. O’Connell JR, Weeks DE: PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet63 :259 –266,1998[Medline]
  24. Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet30 :97 –101,2002[Medline]
  25. Kong A, Cox NJ: Allele-sharing models: LOD scores and accurate linkage tests. Am J Hum Genet61 :1179 –1188,1997[Medline]
  26. Lindgren CM, Widen E, Tuomi T, Li H, Almgren P, Kanninen T, Melander O, Weng J, Lehto M, Groop LC: Contribution of known and unknown susceptibility genes to early-onset diabetes in Scandinavia: evidence for heterogeneity. Diabetes51 :1609 –1617,2002[Abstract/Free Full Text]
  27. Parker A, Meyer J, Lewitzky S, Rennich JS, Chan G, Thomas JD, Orho-Melander M, Lehtovirta M, Forsblom C, Hyrkko A, Carlsson M, Lindgren C, Groop LC: A gene conferring susceptibility to type 2 diabetes in conjunction with obesity is located on chromosome 18p11. Diabetes50 :675 –680,2001[Abstract/Free Full Text]
  28. Bowden DW, Sale MM, Williams AH, Langefeld CD, Rich SS, Freedman BI: Genome-wide search for type 2 diabetes susceptibility genes in African Americans (Abstract). In Proceedings of the 18th Congress of the International Diabetes Federation, Paris, France, 24–29 August 2003. Diabetologia46 (Suppl.2) :A36 –A37,2003
  29. van Tilburg JH, Sandkuijl LA, Strengman E, van Someren H, Rigters-Aris CA, Pearson PL, van Haeften TW, Wijmenga C: A genome-wide scan in type 2 diabetes mellitus provides independent replication of a susceptibility locus on 18p11 and suggests the existence of novel loci on 2q12 and 19q13. J Clin Endocrinol Metab88 :2223 –2230,2003[Abstract/Free Full Text]