1 Diabetes Center, Tokyo Womens Medical University, Tokyo, Japan
2 Departments of Human Genetics, Medicine, and Biochemistry and Molecular Biology, the University of Chicago, Chicago, Illinois
3 Department of Diabetes, Shiseikai Daini Hospital, Tokyo, Japan
4 Department of Endocrinology, Yaizu Citizens Hospital, Shizuoka, Japan
5 Division of Statistical Genetics, Institute of Rheumatology, Tokyo Womens Medical University, Tokyo, Japan
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() |
---|
Type 2 diabetes is one of the most common diseases of middle-age, and its prevalence is increasing in many countries, including Japan (1,2). The increase in prevalence of type 2 diabetes is due, at least in part, to the presence of latent genetic risk factors that are being unmasked by changing patterns of diet, exercise, and other lifestyle/environmental factors. Clinical studies indicate that type 2 diabetes is a phenotypically heterogeneous disorder, and it is likely that it is not one disease but many, with hyperglycemia resulting from different combinations of susceptibility genes superimposed on different nongenetic factors. This heterogeneity is evident on comparing the clinical characteristics of diabetic individuals of different racial groups. For example, Japanese patients are characterized by a lower BMI (25.84 ± 4.01 kg/m2) and lower fasting insulin levels (8.74 ± 6.91 µU/ml) than individuals of European descent (30.43 ± 6.72 kg/m2 and 14.01 ± 13.26 µU/ml), Mexican-Americans (31.61 ± 8.18 kg/m2 and 15.71 ± 11.55 µU/ml), or African-Americans (33.04 ± 8.69 kg/m2 and 17.39 ± 21.42 µU/ml) (3). It has been suggested that ß-cell dysfunction is the primary physiological defect that leads to type 2 diabetes in Japanese subjects, whereas it is insulin resistance in European subjects (4). Both genetic and nongenetic factors are likely to contribute to these clinical differences, and genetic studies in diverse populations may identify the molecular bases for the phenotypic differences between and within populations.
We carried out an autosomal genome scan for type 2 diabetes genes in Japanese subjects. Since obesity is a risk factor for type 2 diabetes (1), we also tested for linkage to log-transformed BMI (lnBMI). The study population (164 families, 368 subjects, and 256 affected sib-pairs [ASPs]), including sibship structure and baseline clinical and biochemical features is summarized in Tables A1 and A2 in the online supplement at http://diabetes.diabetesjournals.org. This sample size has >80% power for detecting a locus with a s >1.6 (suggestive evidence for linkage, P < 0.00074) (5) and 50% power to detect a locus with a
s >1.4.
The study population was genotyped using a panel of 414 autosomal microsatellite markers with an intermarker interval of 8.6 ± 5.2 cM and average heterozygosity of 0.70. A total of 19 markers showed nominal evidence for linkage with type 2 diabetes in two-point analyses (maximum likelihood score [MLS] >0.74, P < 0.05) (Table 1) and 12 markers showed nominal evidence for linkage (logarithm of odds [LOD] >0.59, P < 0.05) in multipoint analyses (Table 1 and Fig. 1). The strongest evidence for linkage with type 2 diabetes was with markers on chromosome 9q (140.0 cM from pter, LOD = 1.40, P = 0.006) and chromosome 21q (48.0 cM, LOD = 1.92, P = 0.001). None of the regions showing evidence for linkage reached the threshold for suggestive (LOD = 2.2, P = 0.0007) or significant evidence for linkage (LOD = 3.6, P = 0.00002) (5).
|
|
We carried out ordered subset analyses using age at diagnosis of diabetes and current BMI as variables to facilitate comparisons with other studies (8) using a threshold of LOD >2.32 (nominal P < 0.001). Three regions exceeded this threshold: 1) chromosome 15q (45.8 cM) with baseline LOD score 0.50 increasing to 2.41 when we used the 55 families with BMI <22; 2) chromosome 21q (48.0 cM) with baseline LOD score 1.92 increasing to 2.42 and 2.59 when we used the leanest 116 families (BMI <24); and 146 families with age at diagnosis <56 years, respectively; and 3) chromosome 9q (140 cM) with baseline LOD score 1.40 increasing to 2.57 in the 136 families with BMI >21.
We also carried out quantitative trait locus (QTL) analyses in 348 subjects for markers linked to current BMI after log transformation using age and sex as covariates. A total of 12 regions showed nominally significant evidence for linkage with lnBMI (Table 2 and Fig. 1). The region on chromosome 2q at 210.5 cM overlaps a region showing evidence of linkage with type 2 diabetes (Fig. 1), suggesting it may harbor a diabetes/obesity gene.
|
We assessed the significance of the linkage results by simulation studies. The results indicate that P values estimated by simulation are substantially different from those calculated analytically assuming an infinitely dense marker map (e.g., Tables 1 and 2). This is due to the fact that there is substantial missing data in families ascertained through two siblings with type 2 diabetes because parents are rarely available for study. In addition, the markers are neither perfectly informative nor infinitely dense. For type 2 diabetes status, a genome-wide P value of 0.05 corresponds to an LOD score of 2.77, a genome-wide P value of 0.10 corresponds to an LOD score of 2.38, and a genome-wide P value of 0.20 corresponds to an LOD score of 2.09. For the QTL BMI, a genome-wide P value of 0.05 corresponds to an LOD score of 3.07, a genome-wide P value of 0.10 corresponds to an LOD score of 2.76, and a genome-wide P value of 0.20 corresponds to an LOD score of 2.40. Thus, even the regions providing the strongest pointwise evidence for linkage with either type 2 diabetes (chromosome 21, 48 cM, LOD = 1.92, empiric genome-wide P value = 0.25) or lnBMI (chromosome 2, 210.5 cM, LOD = 1.45, empiric genome-wide P value = 0.83) do not approach genome-wide levels of significance even allowing for lower information content of our sample.
We also used simulation studies to estimate the number of nominally significant signals that might be expected in a genome-wide screen of these data if no genes were linked to type 2 diabetes or BMI. We found 135 of the 1,000 replicates had 12 or more nominally significant LOD scores (LOD >0.59) for type 2 diabetes. In variance components analyses, 402 of the 1,000 replicates had at least 12 nominally significant LOD scores, but only 71 of 1,000 had at least 9 LOD scores >1.0. These studies suggest that some of the regions providing evidence for linkage to type 2 diabetes or BMI may contain genetic variation affecting these phenotypes. However, there is little in the size of the linkage signal obtained in the primary analyses that would differentiate the true from the false positive signals.
Of the 12 regions, 8 (chromosome 2 169.9 and 236.8 cM; chromosome 4 104.9 cM; chromosome 5 114.8 cM; chromosome 6 42.3 cM; chromosome 8 15.3 cM; chromosome 17 36.1 cM; and chromosome 21 48.0 cM) showing linkage with type 2 diabetes are candidate regions based on other linkage studies (Tables A3 and A4 in the online supplement). Seven (chromosome 2 167.9 and 210.5 cM; chromosome 3 185.7 cM; chromosome 5 131.9 cM; chromosome 7 7.4 cM; chromosome 16 30.0 cM; and chromosome 17 47.8 cM) of the 12 regions linked with lnBMI have shown evidence of linkage with BMI or BMI-related traits in other studies (Table A5 in the online supplement).
This is the third genome scan for type 2 diabetes genes in the Japanese population. The first involved a group of 18 Japanese-American families (45 ASPs) from the Seattle, Washington area (3) and the second was an independent group of 159 Japanese families (224 ASPs) from the Tokyo, Japan, area (8). The 164 families (256 ASPs) that we studied were residents of Tokyo and surrounding areas. (We did not collect information on the ancestral homes of these individuals and it is likely that they include individuals from many different areas of Japan.) No region showing nominal evidence for linkage to type 2 diabetes was common to all three studies (Table A4 in the online supplement). However, six regions showed linkage to type 2 diabetes in two of the three studies. Two regions (chromosome 2 236.8 cM and chromosome 6 42.2 cM; Table A4) overlap between this study and that of Mori et al. (8), and both used relatively large numbers of ASPs ascertained from the Kanto area of eastern Japan. These may be promising regions for follow-up studies in the Japanese population.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() |
---|
Fasting plasma insulin levels were obtained from some subjects not treated with insulin. In the subjects treated with insulin, fasting plasma C-peptide levels were determined. Serum insulin levels were measured by enzyme-linked immunosorbent assay using a COBAS CORE II Insulin-EIA Kit (Roche Diagnostics, Basel, Switzerland), and C-peptide levels were measured by radioimmunoassay.
Genotyping.
Genomic DNA was prepared from peripheral blood. We used simple tandem-repeat polymorphisms (CHLC Human Screening Set 6A and Map Pairs Human Version 9/9aRG) of known map location on the 22 autosomes for the primary genome scan. Primers were obtained from Research Genetics (Huntsville, AL).
Linkage analyses.
The data were screened for the presence of misspecified family relationships using the RELPAIR software (11) and genotyping errors using the SIBMED software (12) before linkage analyses. Two-point linkage analyses with type 2 diabetes were carried out using SPLINK (13). Multipoint analyses to test for linkage with type 2 diabetes were conducted using GENEHUNTER PLUS (14,15). We used the score (pairs) function and the exponential model for all analyses. We contrasted results of analyses in which families were weighted equally (regardless of size, structure, or number of affected individuals), with results of analyses in which families were weighted to the number of (pairwise) independent pairs. These alternative weighting functions will yield different results only for families with four or more affected individuals. We report results for families weighted equally.
Multipoint QTL analyses were conducted using the variance components analysis option of GENEHUNTER PLUS (16) on current BMI after log transformation. Age and sex were included as covariates in these analyses, and only additive variance components were estimated. Simulation studies were conducted to assess genome-wide significance of the LOD scores obtained for both disease (diabetes affection status) and lnBMI. The simulations were conducted using MERLIN (17). The same maps, marker allele frequencies, and pedigree structures (including missing data patterns) used for the actual analyses were used to simulate data for the entire genome. We simulated 1,000 replicates for diabetes status and 1,000 replicates for lnBMI (separate simulations were conducted because of differences in the missing data patterns) and then conducted the same analyses as described above for diabetes and for lnBMI.
We tested for interaction among regions with nominally significant multipoint evidence for linkage with type 2 diabetes by calculating correlations in the nonparametric linkage (NPL) scores across families for the chromosome 9 and 21 regions versus each of the other regions with nominally significant evidence for linkage (6). We also performed ordered subset analyses to see if any region provided stronger evidence for linkage in a clinically defined subgroup. For all analyses, the family was classified according to the average for the variable (age at diagnosis or BMI) in that family. Families were then ordered (e.g., from those with the earliest age at diagnosis to those with the latest age at diagnosis), and subgroups of families were constructed. The subgroups are inclusive. For example, for age at diagnosis, we report results for all families with average age at diagnosis <30, <36, <40, <42, <44, <48, <50, <52, <56, and <60 years. All of those with age at diagnosis <30 are included in the <36 group as well, etc. Although we considered age at onset subset going in the direction from early to late, we considered BMI subsets from both lean and obese directions (>20, >21, >22, >23, >24, >25, and >27 and <20, <21, <22, <23, <24, <25, and <27 kg/m2). The groupings are arbitrary and reflect convenient choices for the sample sizes falling into each group. The maximum number of families for the baseline analysis is 164, and a total of 24 subsets were used for these analyses.
![]() |
ACKNOWLEDGMENTS |
---|
We thank all the families for their participation, without which this study would not have been possible. We also thank Drs. C. Gragnoli and L. del Bosque-Plata for their assistance with genotyping, Dr. A. Pluzhnikov for her assistance with the statistical analyses, C. Roe for her help in the preparation of Fig. 1, and Y. Sagisaka and A. Nogami for their help with sample preparation and genotyping. G.I.B. is an Investigator of the Howard Hughes Medical Institute.
![]() |
FOOTNOTES |
---|
Received for publication 10 August 2002 and accepted in revised form 8 October 2002.
Additional information can be found in an online appendix at http://diabetes.diabetesjournals.org.
ASP, affected sib-pair; lnBMI, log-transformed BMI; LOD, logarithm of odds; MLS, maximum likelihood score; MODY, maturity-onset diabetes of the young; NPL, nonparametric linkage; QTL, quantitative trait locus.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() |
---|