1 Center for Epidemiology and Biostatistics, Cincinnati Childrens Hospital Medical Center, Cincinnati, Ohio
2 Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas
3 Institute of Biology, Institute Pasteur of Lille, Lille, France
4 Centre National de Génotypage, Evry Cedex, France
5 Barts & The London Genome Centre, Queen Mary & Smithfield College, London, U.K
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
A major cause for this difficulty is disease heterogeneity. Given the complexity of glycemic control mechanisms, this heterogeneity is not surprising. For example, the insulin receptor signals by phosporylating at least six intracellular substrates (5). Thus, a mutation affecting production of any of these substrates could result in an increased susceptibility to diabetes.
Indeed, significant linkages have been identified on several chromosomes, including 2q37 (6), 15q21 (7), and 3p (8) in Mexican-Americans; 3q27 (9) in French Caucasians; 10q26 (10), 12q24 (11), and 1q21-q23 (9,12,13) in Caucasians and Pima Indians; 20q (9,14) in Caucasian whites; and 9q21 (15) in Chinese. These results underscore the potential genetic heterogeneity of type 2 diabetes across diverse ethnic groups and the difficulties for replication of linkage.
Given the high degree of genetic heterogeneity, multiple genes with small to moderate effects may influence susceptibility to diabetes. One strategy to circumvent this limitation would be to focus on individuals that are likely to share the same polymorphism, such as in an extended family. Although a single family may share household effects, the use of all possible relative pairs in an extended family minimizes this bias because a large proportion of relative pairs do not share a common household. The purpose of this paper is to search for quantitative trait loci (QTLs) for susceptibility of type 2 diabetes in a single extended family from France.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Affected individuals were either diabetic or glucose intolerant. Although individuals with impaired glucose tolerance (IGT) have not developed diabetes, these individuals have an increased risk of type 2 diabetes (18). Indeed, in a Dutch population (the Hoorn study), 64.5% of the participants who had both impaired fasting and impaired postload glucose levels at baseline progressed to diabetes during a 6-year follow-up (19). Yet, a substantial number of IGT individuals never progress to diabetes. Nonetheless, as these individuals are from a single family, there may be a common genetic mechanism underlying the diabetic and IGT states.
Genotyping.
Genomic DNA was isolated from leukocytes using the Puregene DNA isolation kit, according to the manufacturers protocol (Gentra Systems). Genotyping was performed using a fluorescently labeled human linkage mapping set (PE-LMSV2) comprising 400 highly informative microsatellite markers, with an average intermarker spacing of 9.7 cM. Multiplex PCR conditions were set up to amplify the markers in 87 PCRs. PCR (95°C for 12 min, followed by 30 cycles of 94°C for 15 s, 55°C for 15 s, 72°C of 30 s, and 72°C for 10 min) was performed on GeneAmp PCR system 9700 (Perkin-Elmer) (10-ml reactions; 40 ng of genomic DNA, 2.5 mmol/l MgCl2, and 0.25 mmol/l each dNTP; Pharmacia), with a variable amount (0.21.5 pmol) of 5' and 3' primers and 0.4 units AmpliTaq Gold DNA polymerase (Perkin Elmer) in 1x PCR Bufffer II (Perkin Elmer) (multiplex PCR conditions are available on request).
Pooled amplification products were electrophoresed through 5% polyacrylamide gels (Long Ranger Singel pack; Perkin Elmer) for 1.5 h at 2,000 V on 24-cm plates on an ABI 377 DNA sequencer. An automated 96-channel pipettor Multimek 96 (Beckman) was used for all pipetting steps. Semiautomated fragment sizing was performed using GENESCAN 3.0 software (ABI), followed by allele calling with GENOTYPER 2.1 software (ABI). To confirm the accuracy of allele calling, two individuals reviewed each genotype independently. Average heterozygosity was 0.79. Incompatibilities were identified with PED-CHECK 1.1 (20), and inconsistencies were resolved by blanking (220 of 62,957 genotypes).
Variance components linkage analysis.
Our variance components approach is an extension of the strategy developed by Amos (21), which is based on specifying expected genetic covariances between relatives as a function of the identity by descent (IBD) relationship at the marker locus. We tested the null hypothesis that the additive genetic variance due to a QTL (q2) equals 0 (no linkage) by comparing the likelihood of this restricted model with that of a model in which
q2 is estimated (22). The difference between the two log10 likelihoods produces a logarithm of odds (LOD) score that is equivalent to the classical LOD score of linkage analysis. Twice the difference in loge likelihoods of these models yields a test statistic that is asymptotically distributed as a 1/2:1/2 mixture of a
2 variable and a point mass at 0 (23). This method has been implemented in the program package SOLAR (24), which determines whether genetic variation at a specific chromosomal location can explain the variation in the phenotype (21,22,25).
A pairwise maximum likelihood-based procedure was used to estimate IBD probabilities (24). To permit multipoint analysis for QTL mapping, an extension (23) of the technique of Fulker et al. (26) was used. Estimates of the IBD probabilities were generated at any point on a chromosome using a constrained linear function of observed IBD probabilities of markers at known locations within the region. Exact multipoint IBDs were calculated using SIMWALK; however, results did not differ from the SOLAR-generated IBDs, therefore results are reported for the SOLAR IBDs. A LOD-score evaluation was performed at 1-cM intervals along the chromosome; the distances between markers having been obtained from the sex-averaged maps compiled by Généthon.
Extension to dichotomous traits.
We extended the basic variance components method to dichotomous traits by assuming that an individual belongs to a specific affection status if an underlying genetically determined risk (i.e., liability) exceeds a certain threshold (27). Although disease status is usually considered qualitative, with individuals either scored as affected or unaffected, it is generally assumed that there is an underlying quantitative liability determining affection status. If an individuals liability score exceeds a specified threshold (determined using age- and sex-specific parameters, to produce the appropriate population prevalence), disease ensues. In contrast, if an individuals score is below the threshold, the individual remains unaffected. An integral over the appropriate region of the curve is used to estimate each persons liability value, and the latent liability is assumed to have an underlying normal distribution.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Before quantitative genetic analysis, effects of age and sex were removed. Increasing age was a risk factor for affecteds. Additionally, men had a greater risk of being affected. These covariates explained 15% of the total phenotypic variation in affection status. No evidence of age by sex interactions was found, thus these were not included. BMI was also included as a covariate, and results were similar to those without BMI. Because BMI is a risk factor for type 2 diabetes, we report results of analyses without BMI.
Figure 2 displays results from the linkage analysis for diabetes status by chromosome. We detected one chromosome with suggestive evidence of linkage (LOD score >1.9) (28). A maximum LOD score of 2.4 (P = 0.00044) near marker D5S428 was obtained on chromosome 5 (99 cM) and was confirmed by simulation analyses. The 1-LOD unit support interval, surrounding the peak LOD score on chromosome 5, ranges from 87114 cM from p-terminal (Fig. 3). Three other chromosomes yielded LOD scores >1: chromosome 10 at 89 cM (LOD = 1.0), chromosome 11 at 147 cM (LOD = 1.6), and chromosome 13 at 4 cM (LOD = 1.6). However, these failed to reach the level for suggestive evidence of linkage (28).
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Our results suggest that a gene located near D5S428 is responsible for variability in affection status. As we focused on a single family, the applicability of these findings to the general population is a concern. However, Lindsay et al. (29) identified tentative evidence of linkage (LOD = 1.5) 15 cM from our peak when conditioning on maternal inheritance. Additionally, Hager and et al. (30) identified a QTL on chromosome 5 for serum leptin levels. Our signal for diabetes susceptibility localizes within the 1-LOD unit support interval of this QTL. As leptin is an indicator of total adiposity and obesity is a risk factor for diabetes, the same mechanism may be causing obesity and diabetes in these populations. Therefore, our findings may be generalizable.
Within the 1-LOD unit support interval, there are several mapped genes with an effect on insulin action. A potential candidate in this region is the gene for prohormone convertase 1 (PC1/PC3, SUBTILISIN/KEXIN-TYPE 1, PCSK1), which is implicated in the processing of several prohormones into mature hormones, such as proinsulin and proopiomelanocortin (POMC). The gene products of CPE, PC1, and PC2 cooperate in prohormone processing. By in situ hybridization, Seidah et al. (31) mapped PC1 to human 5q15-q21. Glucose regulates insulin and PC1 expression, supporting the important role of PC1 in regulating proinsulin processing. Although type 2 diabetes is associated with increased secretion of proinsulin and proinsulin-like molecules, screening for mutations of the entire coding region of the PC1 gene in Japanese type 2 diabetes subjects found no mutation associated with type 2 diabetes (32). However, in a woman with extreme childhood obesity, abnormal glucose homeostasis and clinical manifestations suggestive of defective prohormone processing, Jackson et al. (33) identified compound heterozygosity for two mutations of the PC1 gene. Three of their four clinically unaffected children had the gly483-to-arg missense mutation, while the fourth child had the splice site mutation.
Another positional candidate gene encodes calpastatin (CAST), which is an endogenous modulator of intracellular calpain activity and is localized within the same cytoplasmic region as calpains. Calapastatin is a calcium-dependent cystein protease probably involved in muscle protein degradation in living tissue. CAST has been mapped to 5q15-q21 by in situ hybridization and spot-blot analysis of sorted chromosomes (34). Although the physiological role of calpains remains unclear, Horikawa and et al. (35) described an association between type 2 diabetes and a genetic variant identified in the gene encoding calpain-10. Calpains are ubiquitous nonlysosomal cystein proteases, which have been implicated in a wide variety of cellular functions, such as cell proliferation and differentiation, intracellular signaling, and apoptosis. Calpastatin also downregulates genes encoding several membrane-associated proteins or nuclear proteins and upregulates genes of collagen 2. Hence, in addition to their proteolytic activities on cytoskeletal proteins and other cellular regulatory proteins, calpain-calpastatin systems can affect expression levels of genes encoding structural or regulatory proteins.
Although no other regions reached the level of suggestive linkage, the signal on chromosome 11 may be a replication of the QTL for diabetes on chromosome 11 originally detected in Pima Indians (13). We obtained an LOD of 1.6 on chromosome 11 near marker D11S1320. Based on Marshfield map distances (http://research.marshfieldclinic.org/genetics/), this signal is <20 cM from the QTL for BMI and diabetes in Pima Indians (D11S4464).
It is important to note that we grouped diabetic and IGT subjects in the affected category. However, not all IGT individuals progress to diabetes, therefore by using IGT instead of diabetes status we may be misclassifying individuals. Although it would be interesting to examine the effect of this misclassification (including IGTs who never progress to diabetes status) on our linkage signal, the examination of such an effect would require that we have some predictive measure of who is most likely to convert to type 2 diabetes. As we do not have such a measure, estimation of the effect of misclassification would be random, therefore not informative. Moreover, the identification of genes that are involved in both type 2 diabetes and IGT would clinically relevant because early intervention for high-risk individuals may reduce the number of complications associated with type 2 diabetes. Biologically, IGT is associated with insulin resistance, while diabetes is associated with insulin resistance and dysregulation of insulin secretion. Therefore, our approach would likely identify genes involved with insulin resistance genes. However, this approach would likely lessen our ability to detect a genetic signal associated with genes influencing insulin secretion.
In summary, we obtained suggestive evidence of linkage (LOD = 2.4) for type 2 diabetes susceptibility on chromosome 5q in a single extended family. This signal is in the same region as a QTL for diabetes when conditioning on maternal inheritance (28) and as a QTL for leptin levels (30). Within the 1-LOD unit support interval, there are two strong candidates: PCSK1 and CAST. Furthermore, we have replicated localization of a QTL for diabetes status on chromosome 11 detected by Hanson et al. (13). These findings suggest that both novel and previously identified QTLs may be influencing the development of diabetes.
![]() |
ACKNOWLEDGMENTS |
---|
We are indebted to the family members who participated in this study.
![]() |
FOOTNOTES |
---|
Received for publication 18 April 2002 and accepted in revised form 16 September 2002.
IBD, identity by descent; IGT, impaired glucose tolerance; LOD, logarithm of odds; QTL, quantitative trait locus.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|