Genetic Mapping of Disposition Index and Acute Insulin Response Loci on Chromosome 11q
http://www.100md.com
糖尿病学杂志 2006年第4期
1 Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, North Carolina
2 Center for Human Genomics, Wake Forest University School of Medicine, Winston-Salem, North Carolina
3 Department of Public Health Sciences, Wake Forest University School of Medicine, Winston-Salem, North Carolina
4 Department of Medicine, University of California at Los Angeles, Los Angeles, California
5 Department of Preventive Medicine and Biometrics, University of Colorado Health Sciences Center, Denver, Colorado
6 Department of Medicine, University of Texas Health Sciences Center at San Antonio, San Antonio, Texas
7 Medical Genetics Institute, Burns and Allen Cedars-Sinai Research Institute, Los Angeles, California
8 Department of Physiology and Biophysics, Keck School of Medicine, University of Southern California, Los Angeles, California
9 Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, North Carolina
AIR, acute insulin response; DI, disposition index; GEE, generalized estimating equation; IBD, identity-by-descent; IRASFS, Insulin Resistance Atherosclerosis Study Family Study; LOD, logarithm of odds; MAF, minor allele frequency; MLS, maximum LOD score; PC, pyruvate carboxylase; PYGM, muscle glycogen phophorylase; QTL, quantitative trait locus; SNP, single nucleotide polymorphism
ABSTRACT
Glucose homeostasis, a defining characteristic of physiological glucose metabolism, is the result of complex feedback relationships with both genetic and environmental determinants that influence insulin sensitivity and -cell function. Relatively little is known about the genetic basis of glucose homeostasis phenotypes or their relationship to risk of diabetes. Our group previously published a genome scan for glucose homeostasis traits in 284 African-American subjects from 21 pedigrees in the Insulin Resistance Atherosclerosis Study Family Study (IRASFS) and presented evidence for linkage to disposition index (DI) on chromosome 11q with a logarithm of odds (LOD) of 3.21 at 81 cM flanked by markers D11S2371 and D11S2002 (support interval from 71 to 96 cM). In this study, genotyping and analysis of an additional 214 African-American subjects in 21 pedigrees from the IRASFS yielded independent evidence of linkage to DI. When these two datasets were combined, a DI linkage peak was observed with an LOD of 3.89 at 78 cM (support interval from 67 to 89 cM). Fine mapping with 15 additional microsatellite markers in this 11q region for the entire 42 pedigrees resulted in an LOD score of 4.80 at 80 cM near marker D11S937 (support interval from 76 to 84 cM). In these 42 pedigrees, there was also suggestive evidence for linkage to acute insulin response (AIR) at two separate locations flanking the DI peak (64 cM, LOD 2.77, flanked by markers D11S4076 and D11S981; and 85 cM, LOD 2.54, flanked by markers D11S4172 and D11S2002). No evidence of linkage to the insulin sensitivity index (Si) was observed. Nine positional candidate genes were evaluated for association to DI and AIR. Among these candidates, single nucleotide polymorphisms (SNPs) in muscle glycogen phosphorylase showed evidence of association with DI (P < 0.011). In addition, SNPs in the pyruvate carboxylase gene showed evidence of association (P < 0.002) with AIR. Further analysis of these candidate genes, however, did not provide evidence that these SNPs accounted for the evidence of linkage to either DI or AIR. These detailed genetic analyses provide strong evidence of a DI locus on 11q in African-American pedigrees, with additional suggestive evidence of independent AIR loci in the same region.
It is widely accepted that type 2 diabetes is due to a combination of defects in insulin sensitivity and insulin secretion that result in impaired glucose homeostasis. Disposition index (DI) is a quantitative measure that describes the relationship between -cell sensitivity and insulin sensitivity (1). In metabolically normal individuals, changes in insulin sensitivity are accompanied by compensatory alterations in the -cell’s sensitivity to glucose. In practice, DI is measured as the product of the insulin sensitivity index (Si) and -cell function as measured by the acute insulin response (AIR). Despite being key features of glucose homeostasis, the molecular mechanism(s) by which insulin sensitivity and -cell function are coregulated to create a homeostatic environment are not well understood (1eC3).
Genetic analysis in appropriately phenotyped families is a potentially powerful approach for gaining insights into complex biological traits such as DI. To date, few mapping studies have incorporated the analysis of DI as a quantitative trait (4,5). Our group recently reported evidence for a DI locus between D11S2371 and D11S2002 on chromosome 11q at 81 cM with a logarithm of odds (LOD) score of 3.21 (support interval 71eC96 cM) in 284 African-American subjects from 21 pedigrees in the Insulin Resistance Atherosclerosis Study Family Study (IRASFS) (6). In addition to DI, other measures of glucose homeostasis have been previously reported to be linked to this region (7eC10).
The purpose of the current study is threefold. First, we report confirmatory evidence that DI and other measures of glucose homeostasis link to chromosome 11q in an independent and similarly phenotyped and ascertained sample of African-American pedigrees from the IRASFS. Second, we report fine-mapping quantitative trait locus (QTL) linkage analysis of glucose homeostasis measures for this region on 11q. Third, we report an evaluation of selected candidate genes within 11q and explore their contribution to the variation in measures of glucose homeostasis.
RESEARCH DESIGN AND METHODS
The study design, recruitment, and phenotyping in the IRASFS have been described in detail (11). Briefly, the IRASFS is a multicenter study designed to identify the genetic determinants of insulin resistance through quantitative, intermediate measures of glucose homeostasis. For this component of the study, members of large families of self-reported African-American descent from the Los Angeles, California, area were recruited and examined. Ascertainment of families was not based on phenotypes but instead on large family size. A clinical examination was performed that included interviews, a frequently sampled intravenous glucose tolerance test, anthropometric measurements, abdominal computed tomography scanning, resting blood pressure, blood collection, and spot urine. Specific to this report, measures of glucose homeostasis included those from the frequently sampled intravenous glucose tolerance test: Si, AIR, and DI.
Genotyping
Genome-wide scan and fine mapping.
The genotyping procedure has been previously described (6). Briefly, 383 microsatellite markers (average one marker/9.3 cM) were genotyped at the National Heart, Lung, and Blood Institute Mammalian Genotyping Service (Marshfield, WI; Marshfield set 10) in two independent sets of pedigrees, denoted as set 1 and set 2. Based on the evidence of linkage to 11q in set 1 (6) and the replication of the linkage in set 2 reported herein, an additional 15 microsatellite markers (fine-mapping markers) were genotyped on all African-American pedigrees (set 1 and set 2) within a 62-cM region centered on the linkage peak. Genotyping data from these 15 microsatellite markers were merged with the genome scan data for the fine-mapping analysis.
Candidate genes.
Nine positional candidate genes were chosen within the linked region based on their hypothesized functional role in glucose homeostasis. Initially, an average of five approximately equally spaced single nucleotide polymorphisms (SNPs) per positional candidate gene were selected from dbSNP (http://www.ncbi.nlm.nih.gov/projects/SNP/). SNPs with available frequency information were preferentially selected from the database. Following preliminary analysis, additional SNPs were chosen in positional candidate genes that showed association to the glucose homeostasis phenotypes DI, AIR, or Si. These additional SNPs were chosen to provide greater coverage of the linkage disequilibrium block and generally had a minor allele frequency (MAF) >0.20 in the IRASFS African-American pedigrees. Genotyping was performed on the Sequenom MassArray Genotyping System using methods previously described (12). Twenty-two individuals from two families served as blind duplicates (duplicated samples given new unique identifiers) to evaluate genotyping accuracy and exhibited a 0.48% discordance percentage.
Statistical analysis.
The statistical analyses reported are for the three primary purposes of this study: 1) confirmatory test that DI links to 11q in an independent and similarly phenotyped and ascertained replication sample of African-American pedigrees from the IRASFS, 2) fine mapping of both African-American samples on 11q, and 3) association analysis of selected candidate genes within the 11q region.
Genome-wide scan QTL linkage analysis.
The IRASFS genome scan was completed in two phases, denoted set 1 (284 African-American subjects from 21 pedigrees) and set 2 (214 African-American subjects from 21 pedigrees). The analysis of the first set of 21 African-American pedigrees showed evidence of linkage to 11q and has been previously reported (6). Since that publication, the locations of markers D11S2365 and ATA9B04 on chromosome 11 that were reported at 58 cM on the Marshfield map were changed to 55.96 and 56.27 cM, respectively, in the recent release of the human genome physical map (NCBI build 35, hg assembly 17). These shifts in location resulted in changes in marker order within the region, requiring these data to be reanalyzed due to the sensitivity of linkage analyses to marker order. The reanalysis for linkage of the first set of 21 pedigrees and the linkage analysis of the second set of 21 pedigrees estimated the identity-by-descent (IBD) statistics using maximum likelihood estimates of allele frequencies from the African-American pedigrees and the Monte Carlo approach implemented in the computer program LOKI (13,14).
For the analysis of set 1, set 2, and the combined dataset, quantitative traits were transformed to best approximate the distributional assumptions of the variance-component QTL linkage analysis (i.e., conditional normality after adjusting for covariates, homogeneity of variance). Thus, the results reported represent analyses on the square root of DI and AIR and the natural logarithm of Si after adding 1. For individuals with a self-reported diagnosis of diabetes or fasting glucose >126 mg/dl, the glucose homeostasis parameters of insulin sensitivity (Si), DI, and AIR were converted to missing. However, all available genotypic data were used to compute the IBD statistics using the Bayesian Markov Chain Monte Carlo approach implemented in the software LOKI (13,14). Previously, self-reported familial relationships were examined and modified using the genome scan data from both set 1 and set 2 pedigrees and the software PREST (15). Testing for evidence of linkage to a QTL was done using the variance-component approach implemented in the SOLAR software package (16). All LOD scores reported herein are empirical LOD scores determined within SOLAR by simulation (lodadj procedure). The entropy-based measure of linkage information content was computed using Merlin (17).
Fine mapping.
Map order and genetic distance for each of the fine-mapping markers were taken from the deCODE genetic map (18). Markers for fine mapping were analyzed using PedCheck (19) to identify any Mendelian inconsistencies. Ninety-nine total genotypes (from 79 individuals across 26 families) were converted to missing. The genotyping data from these 15 microsatellite markers were then merged with the combined genome scan genotyping data from set 1 and set 2, and QTL linkage analyses were computed as described above.
Candidate genes.
The genetic location of each positional candidate gene was calculated by extracting all markers from the deCODE genetic map (18) centered around the linkage region (D11S1357:D11S4104) and mapped to the chromosomal sex-averaged genetic positions to the latest human genome physical map positions (NCBI build 35, hg assembly 17). Using R version 1.9.1 (http://www.r-project.org/index.html), a smoothing spline was fit through these points and the spline used to interpolate genetic map positions corresponding to physical positions of gene loci. For each gene, we computed a single representative physical map position as the mean of the 5' end of the most telomeric exon in any known annotated gene splice variant and the 3' end of the most centromeric, assuming the genetic interval of the entire locus was small relative to that of the linkage region.
Initially, each of the SNPs evaluated in this work were examined for Mendelian inconsistencies in their genotypes using PedCheck (19). Any genotypes inconsistent with Mendelian inheritance were converted to missing (0.07%). Subsequently, maximum likelihood estimates of allele frequencies were computed using the largest set of unrelated individuals from the combined dataset and then tested for departures from Hardy-Weinberg proportions. To test for an association between individual SNPs and each quantitative phenotype (Si, AIR, and DI), a series of generalized estimating equations (GEEs; GEE1) was computed (20). The correlation between subjects within a pedigree was adjusted in the analyses by assuming an exchangeable correlation structure among siblings within a pedigree and computing the sandwich estimator of the variance (21). The sandwich estimator is also denoted the robust or empirical estimator of the variance, as it is robust to misspecification of the correlation matrix because it estimates the within-pedigree correlation matrix from moments of the data. For each phenotype (Si, AIR, and DI) and each SNP, the 2-df (degrees of freedom) test of genotypic association was performed. In addition, three individual contrasts defined by the a priori genetic models (dominant, additive, and recessive) were computed. The contrasts are defined by the genetic models (i.e., dominant model contrasts those with versus without the polymorphism, additive model tests for a dose effect in the number of alleles, and recessive model contrasts individuals homozygous for the polymorphisms versus those who are not) and are tested using a Wald test. If the overall genotypic association is significant, the a priori contrasts are examined directly. If the overall genotypic association is not significant, the a priori contrasts are examined after adjusting for the three comparisons using a Bonferroni adjustment. Tests reported here were computed adjusting for age, sex, and BMI. Subjects with type 2 diabetes were excluded from the analysis of glucose homeostasis traits.
Statistical significance for all tests was evaluated by a simulation study to obtain empiric P values. Specifically, for each SNP and each trait, the gene-dropping algorithm implemented in Mendel (22) was used and the analysis repeated 10,000 times. The empiric P value is the proportion of time that the simulated data had a P value greater than or equal to the observed P value from the real data. This approach allows the retention of the original correlation among pedigree members in the phenotypic and genotypic data.
The pairwise linkage disequilibrium (LD) statistics D' and r2 were computed, and a block of LD was defined based on marker-marker D' > 0.70 for all markers in the maximal set of unrelated individuals. Haplotypes were constructed within the LD block from SNPs with an MAF >10% in positional candidate genes that showed evidence of association to glucose homeostasis phenotypes. Analysis of these haplotypes was completed using a GEE1 analysis as described above, except the quasi-likelihood was weighted by the probability for each possible haplo-genotype for an individual. Specifically, we computed the expectation-maximization algorithm estimates for haplotype frequencies from the family data using the software ZAPLO (23). To reduce the complexity of the problem, ZAPLO assumes zero recombination between markers and generates all possible haplo-genotypes (i.e., haplotype combinations) for each individual. These haplo-genotypes are used as input into PROFILER (24) to compute the joint probability distribution for the haplotypes of individuals within the pedigree. This joint probability distribution is then used to estimate the probability of each haplotype pair combination (haplo-genotype) possible for each individual, conditional on the family data. Each individual enters into the GEE1 analysis once for each haplo-genotype possibility, weighted by the haplo-genotype probability. Thus, the weight for each individual sums to 1. The weighted GEE1 analyses were completed as above using the same transformations and sandwich estimator of the variance to account for the within-cluster correlation for each of the haplotypes observed with >5% frequency.
To test whether a particular SNP, either directly or through LD with the QTL, contributes to the evidence for linkage, each SNP under a dominant model was entered into the QTL linkage analysis and the change in the magnitude of the LOD score calculated. If the polymorphism directly or indirectly contributes to the evidence for linkage, the initial LOD score will be reduced in a model that includes the polymorphism as a covariate.
RESULTS
The descriptive characteristics of set 1 and set 2 participants are shown in Table 1. Diabetes-affected subjects (12.3%) were excluded from analysis of glucose homeostasis phenotypes. On average, the pedigrees in set 2 had fewer participants than pedigrees in set 1. In addition, set 2 was younger (41.5 ± 12.9), had a greater mean BMI (31.1 ± 6.9 kg/m2), and included more women (63.5%) than men.
Mapping of DI.
A summary of the linkage results for measures of glucose homeostasis on chromosome 11 is presented in online appendix Table 1 (available at http://diabetes.diabetesjournals.org). The variance-component linkage analysis for DI from SOLAR in set 1 was reanalyzed using the new marker order from the latest release of the human genome (NCBI build 35, hg assembly 17) and the LOKI (13,14) IBD calculation software. This reanalysis yielded a maximum LOD score (MLS) of 3.10 at 84 cM, flanked by markers D11S2371 and D11S2002, and a support interval from 72 to 93 cM. In set 2, the MLS was 1.26 at 102 cM, flanked by markers D11S2000 and D11S1986, and a support interval from 80 to 120 cM. In the combined set (set 1 + set 2), the MLS was 3.89 at 78 cM, flanked by markers D11S2371 and D11S2002, and a support interval from 67 to 89 cM.
With evidence of linkage of DI to 11q, fine mapping of the region was carried out by genotyping an additional 15 microsatellite markers across the region to give an average density of one marker/3 cM, with the largest gap being 6.94 cM. A QTL linkage analysis combining the genome scan markers and the 15 additional fine-mapping markers (within the 63.8-cM interval centered at the linkage peak) was completed in the entire African-American sample (Fig. 1A and online appendix Table 1). In set 1, the MLS for DI decreased slightly to 3.03 at 81 cM, flanked by markers D11S937 and D11S4172, and a support interval from 75 to 86 cM. In set 2, the MLS increased to 2.29 at 79 cM, flanked by markers D11S2371 and D11S937, and a support interval from 69 to 84 cM. Thus, the MLS maximized within 2 cM in these two independent samples, with set 2 replicating the QTL linkage effect identified in set 1. In the combined set, the MLS was increased to 4.80 at 80 cM at marker D11S937 with a support interval from 76 to 84 cM. Information content of the markers in this region was assessed using the program Merlin (17). Multipoint information content ranged from 0.83 to 0.95 across the fine-mapping region (59.2eC123.0 cM) (data not shown).
Mapping of AIR.
Variance-component linkage analyses of AIR and Si, components of DI, were also carried out with the entire marker set including the fine-mapping markers (Figs. 1B and online appendix Table 1). Linkage analysis of AIR showed suggestive evidence of linkage with two peaks flanking the DI peak in the combined data (set 1, set 2, and fine-mapping markers) analysis. The most significant linkage evidence for AIR was an MLS of 2.77 at 64 cM, flanked by markers D11S4076 and D11S981, and a support interval from 60 to 73 cM. In addition, a second MLS of 2.54 at 85 cM, flanked by markers D11S4172 and D11S2002, and a support interval from 59 to 91 cM was also observed. Si showed no significant evidence of linkage to this region.
Positional candidate gene analysis.
A total of nine positional candidate genes were evaluated within the linked region for association with glucose homeostasis traits. The location and possible function of these genes in modulation of glucose homeostasis is shown in Table 2. In the association analysis of SNPs genotyped within these genes (online appendix Table 2), the most significant evidence of association was observed with muscle glycogen phophorylase (PYGM) and pyruvate carboxylase (PC) precursor. SNPs in PYGM and PC were consistent with Hardy-Weinberg proportions with the exception of two SNPs in PYGM (rs565688, P = 0.009; rs477549, P = 0.025).
As seen in Table 3, multiple PYGM SNPs were significantly associated with DI from rs625172 through rs2073798 (all P < 0.032), with the exception of SNP rs686171, which trended toward association (P = 0.097). Simulations were carried out to calculate empiric P values for the individual SNPs. SNPs with an MAF >10% remained significantly associated with DI (empiric P < 0.041), with the exception of SNP rs565688 (empiric P = 0.068). Examination of LD in this region suggests one LD block encompassing SNPs rs625172 through rs2073798 (inter-SNP D' 0.96 over 10.7 kb) (online appendix Table 3), which is consistent with the association results. Results of the haplotypic association analysis of PYGM are summarized in Table 4. Three common haplotypes (>5%) were estimated within the LD block. Of these haplotypes, haplotype TAA, with a frequency of 29%, was significantly associated with DI (P = 0.020) following a dominant model of inheritance.
The association analysis results for PC with AIR are presented in Table 5. PC was significantly associated with AIR for SNPs rs896438 through rs7926787 (P < 0.002), with the exception of rs7932018 (P = 0.602) and rs2077432 (P = 0.235). Simulations provided evidence of association with AIR (empiric P < 0.014). Examination of LD in this region revealed one LD block encompassing SNPs rs7119426 through rs6591223 (inter-SNP D' 0.74 over 40.3 kb) (online appendix Table 4) and containing the trait-associated SNPs. Results of the haplotypic association analysis of PC are presented in Table 6. Five haplotypes with frequencies >5% were estimated within the LD block. Of these haplotypes, haplotype CACAGA, with a frequency of 20%, was significantly associated with AIR in all three genetic models, with the recessive model providing the strongest evidence of association (P < 0.001).
To assess the contribution to linkage, individual SNPs from all positional candidate genes were used as covariates in variance-component linkage analysis. For the 17 SNPs analyzed in PC and PYGM, the largest change in MLS occurred with SNP rs2167457 in PC with AIR. This SNP had an MAF of 0.24 and was significantly associated with AIR (empiric P = 0.004). Initially, the peak for AIR had an MLS of 2.76 at 64 cM. In our subset of participants genotyped for this SNP, the peak for AIR had an MLS of 2.29 at 64 cM. After adjustment for SNP rs216457, the largest LOD score change was a decrease to an MLS of 1.90 following a recessive model (data not shown). While this decrease in LOD suggests modest contribution of this SNP to the evidence of linkage, it does not reach statistical significance. Overall, when used as covariates in the QTL linkage analysis, single SNPs in the positional candidate genes did not significantly alter evidence of linkage.
DISCUSSION
We have carried out a detailed genetic analysis of glucose homeostasis traits on chromosome 11q in African-American families from the IRASFS. These QTL linkage analyses included genetic mapping in an independent and identically phenotyped and ascertained replication sample of pedigrees, fine mapping, and analysis of a set of relevant candidate genes for association with glucose homeostasis traits and their influence on the evidence for linkage. In two independent samples (set 1 and set 2), we observed evidence for linkage to DI maximizing within 2 cM near 80 cM in both samples. QTL linkage analyses of the entire 42 IRASFS African-American sample, including both the genome scan markers and 15 additional microsatellite markers, further delineated this region. QTL linkage analyses provide evidence of a DI locus on chromosome 11q at 80 cM (LOD 4.80) and suggestive evidence of linkage to AIR with peaks that appear to flank the DI linkage signal. The results raise the possibility that multiple loci within the region may be contributing to evidence of linkage with measures of glucose homeostasis.
The biologic and mathematical relationship among DI, AIR, and Si (DI = AIR x Si) and the coincident linkages of DI and AIR in this region suggested that we explore the possibility that pleiotrophy may exist within this region. We performed a series of bivariate QTL linkage analyses (AIR x Si, AIR x DI, and DI x Si) to test this hypothesis and did not observe either a meaningful increase or decrease in the respective LOD scores. It is important to recognize that the bivariate linkage and pleiotrophy tests do not account for multiple loci in linkage with each other. Thus, the results of the linkage test are inconclusive relative to pleiotrophy, and a more definitive test will be performed through detailed SNP mapping and association tests.
These results represent the first study to identify strong evidence of a QTL for DI. This region of chromosome 11 has, however, been linked to alternative measures of glucose homeostasis not measured in this study. Pratley et al. (10) reported evidence of linkage to 2-h postload insulin with an MLS of 1.31 at marker D11S2371 (82.6 cM) in Pima Indians. This marker is one that flanks our linkage peak for DI. In addition, the FUSION (Finland-United States Investigation of NIDDM Genetics) study reported preliminary findings of linkage to type 2 diabetes (MLS of 1.75 at 84 cM near D11S1314), fasting insulin (MLS of 2.07 at 89.5 cM), and fasting insulin-to-glucose ratio (MLS of 2.32 at 91.0 cM) in diabetic subjects. Evidence for replication of this linkage to type 2 diabetes with an MLS of 2.98 at 82 cM near D11S4172 was subsequently reported by these investigators (7). Most recently, the Framingham Offspring Study reported evidence of linkage to this region for fasting insulin (MLS of 2.43 at 85 cM near D11S2002) and homeostasis model assessment (MLS 2.01 at 92 cM near D11S2002 and D11S2000) (8).
Nine positional candidate genes were evaluated for contributions to the evidence of linkage. Analysis of PYGM and PC resulted in the most significant consistent evidence of association with DI and AIR, respectively. PYGM catalyzes and regulates the breakdown of glycogen to glucose-1-phosphate by removing the glucose residues from -(1,4)-linkages within glycogen molecules. PYGM was significantly associated with DI within a single LD block. Additionally, SNPs rs565688 and rs625172, which were significantly associated with DI, were also significantly associated with Si (rs565688, P = 0.001; rs625172, P = 0.017) (online appendix Table 2). Therefore, the evidence of association seen with DI could be driven by Si, as the two measures are mathematically and metabolically related (DI = Si x AIR). This may also explain the absence of a contribution to the DI linkage peak when these SNPs were used as covariates in the variance-component linkage analysis. Haplotype analysis of PYGM SNPs revealed a 29% frequency haplotype (TAA) significantly associated with DI (P = 0.020). This result is nominally significant, but there appears to be no apparent trend in the genotypic trait means associated with this haplotype.
PC is a nuclear-encoded mitochondrial enzyme that catalyzes the conversion of pyruvate to oxaloacetate in a tissue-specific manner; the initial reaction of glucose (liver, kidney) and lipid (adipose tissue, liver, brain) synthesis from pyruvate. PC was significantly associated with AIR in a single LD block. In haplotype analysis within the LD block, a 20% frequency haplotype (CACAGA) was significantly associated with AIR (P < 0.001). The trend in the genotypic trait means associated with this haplotype would suggest that its presence is associated with "risk" (i.e., lower AIR). A second less frequent (5%) haplotype (TGCACA) was also significantly associated with AIR (P = 0.001). It is interesting to note that Antinozzi and Rothman (25) have recently observed evidence that RNAi inhibition of PC gene expression substantially suppresses secretion of an insulin surrogate in an in vivo cell culture system, a role that is consistent with our evidence of genetic association with AIR. Unfortunately, as with PYGM, there appears to be no consistent trend in the genotypic trait means associated with these haplotypes.
To examine their contribution to the evidence of linkage, each SNP in the candidate genes was analyzed as a covariate in variance-component linkage analysis for traits Si, AIR, and DI. Individually, no single SNP significantly altered the linkage peaks. This result would preclude a direct contribution of these individual SNPs to the observed evidence of linkage. However, from these results alone, we cannot confidently rule out the possibility that these two positional candidate genes (PYGM and PC) are contributors to the overall evidence of linkage to this region on chromosome 11. Causal variants associated with the variance in DI and AIR could be located within the PYGM and PC gene regions, respectively. If the allele frequencies of the SNPs tested in these genes and the actual (unknown) causal variants differ substantially or are not in LD (e.g., low inter-SNP r2 values between the loci examined), it may not be possible to detect the effects of these loci (26). Additionally, the results from this study are consistent with the possibility that multiple genes, each with a potentially limited contribution to the trait, result in the observed linkage peak. It is important to also note that 100 known and hypothetical genes are located in the support interval (76eC84 cM) for DI. Further insights into the genes contributing to DI and AIR in this region may result from detailed SNP fine mapping currently in progress for this region.
ACKNOWLEDGMENTS
This research was supported in part by National Institutes of Health Grants HL060894 (to D.W.B.), HL060919 (to S.M.H.), HL060931 (to R.N.B.), HL060944 (L.E.W.), HL061019 (to J.M.N.), and HL061210 (to M.S.).
FOOTNOTES
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Bergman RN, Ader M, Huecking K, Van Citters G: Accurate assessment of -cell function: the hyperbolic correction. Diabetes 51 (Suppl. 1):S212eC220, 2002
Kahn SE, Prigeon RL, McCulloch DK, Boyko EJ, Bergman RN, Schwartz MW, Neifing JL, Ward WK, Beard JC, Palmer JP, Porte D Jr: Quantification of the relationship between insulin sensitivity and -cell function in human subjects: evidence for a hyperbolic function. Diabetes 42:1663eC1672, 1993
Bergman RN, Finegood DT, Kahn SE: The evolution of beta-cell dysfunction and insulin resistance in type 2 diabetes. Eur J Clin Invest 32 (Suppl. 3):35eC45, 2002
Watanabe RM, Ghosh S, Langefeld CD, Valle TT, Hauser ER, Magnuson VL, Mohlke KL, Silander K, Ally DS, Chines P, Blaschak-Harvan J, Douglas JA, Duren WL, Epstein MP, Fingerlin TE, Kaleta HS, Lange EM, Li C, McEachin RC, Stringham HM, Trager E, White PP, Balow J Jr, Birznieks G, Chang J, Eldridge W: The Finland-United States investigation of non-insulin-dependent diabetes mellitus genetics (FUSION) study. II. An autosomal genome scan for diabetes-related quantitative-trait loci. Am J Hum Genet 67:1186eC1200, 2000
An P, Hong Y, Weisnagel SJ, Rice T, Rankinen T, Leon AS, Skinner JS, Wilmore JH, Chagnon YC, Bergman RN, Bouchard C, Rao DC: Genomic scan of glucose and insulin metabolism phenotypes: the HERITAGE Family Study. Metabolism 52:246eC253, 2003
Rich SS, Bowden DW, Haffner SM, Norris JM, Saad MF, Mitchell BD, Rotter JI, Langefeld CD, Wagenknecht LE, Bergman RN: Identification of quantitative trait loci for glucose homeostasis: the Insulin Resistance Atherosclerosis Study (IRAS) Family Study. Diabetes 53:1866eC1875, 2004
Silander K, Scott LJ, Valle TT, Mohlke KL, Stringham HM, Wiles KR, Duren WL, Doheny KF, Pugh EW, Chines P, Narisu N, White PP, Fingerlin TE, Jackson AU, Li C, Ghosh S, Magnuson VL, Colby K, Erdos MR, Hill JE, Hollstein P, Humphreys KM, Kasad RA, Lambert J, Lazaridis KN, Lin G, Morales-Mena A, Patzkowski K, Pfahl C, Porter R, Rha D, Segal L, Suh YD, Tovar J, Unni A, Welch C, Douglas JA, Epstein MP, Hauser ER, Hagopian W, Buchanan TA, Watanabe RM, Bergman RN, Tuomilehto J, Collins FS, Boehnke M: A large set of Finnish affected sibling pair families with type 2 diabetes suggests susceptibility loci on chromosomes 6, 11, and 14. Diabetes 53:821eC829, 2004
Panhuysen CI, Cupples LA, Wilson PW, Herbert AG, Myers RH, Meigs JB: A genome scan for loci linked to quantitative insulin traits in persons without diabetes: the Framingham Offspring Study. Diabetologia 46:579eC587, 2003
Ghosh S, Watanabe RM, Valle TT, Hauser ER, Magnuson VL, Langefeld CD, Ally DS, Mohlke KL, Silander K, Kohtamaki K, Chines P, Balow J Jr, Birznieks G, Chang J, Eldridge W, Erdos MR, Karanjawala ZE, Knapp JI, Kudelko K, Martin C, Morales-Mena A, Musick A, Musick T, Pfahl C, Porter R, Rayman JB: The Finland-United States investigation of non-insulin-dependent diabetes mellitus genetics (FUSION) study. I. An autosomal genome scan for genes that predispose to type 2 diabetes. Am J Hum Genet 67:1174eC1185, 2000
Pratley RE, Thompson DB, Prochazka M, Baier L, Mott D, Ravussin E, Sakul H, Ehm MG, Burns DK, Foroud T, Garvey WT, Hanson RL, Knowler WC, Bennett PH, Bogardus C: An autosomal genomic scan for loci linked to prediabetic phenotypes in Pima Indians. J Clin Invest 101:1757eC1764, 1998
Henkin L, Bergman RN, Bowden DW, Ellsworth DL, Haffner SM, Langefeld CD, Mitchell BD, Norris JM, Rewers M, Saad MF, Stamm E, Wagenknecht LE, Rich SS: Genetic epidemiology of insulin resistance and visceral adiposity: the IRAS Family Study design and methods. Ann Epidemiol 13:211eC217, 2003
Buetow KH, Edmonson M, MacDonald R, Clifford R, Yip P, Kelley J, Little DP, Strausberg R, Koester H, Cantor CR, Braun A: High-throughput development and characterization of a genomewide collection of gene-based single nucleotide polymorphism markers by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Proc Natl Acad Sci U S A 98:581eC584, 2001
Heath SC, Snow GL, Thompson EA, Tseng C, Wijsman EM: MCMC segregation and linkage analysis. Genet Epidemiol 14:1011eC1016, 1997
Heath SC: Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet 61:748eC760, 1997
McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet 66:1076eC1094, 2000
Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet 62:1198eC1211, 1998
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30:97eC101, 2002
Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, Richardsson B, Sigurdardottir S, Barnard J, Hallbeck B, Masson G, Shlien A, Palsson ST, Frigge ML, Thorgeirsson TE, Gulcher JR, Stefansson K: A high-resolution recombination map of the human genome. Nat Genet 31:241eC247, 2002
O’Connell JR, Weeks DE: PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 63:259eC266, 1998
Zeger SL, Liang KY: Longitudinal data analysis for discrete and continuous outcomes. Biometrics 42:121eC130, 1986
Hardin J, Hilbe J: Generalized Estimating Equations. New York, Chapman & Hall/CRC, 2003
Lange KCR, Horvath S, Perola M, Sabatti C, Sinsheimer J, Sobel E: Mendel version 4.0: a complete package for the exact genetic analysis of discrete traits in pedigree and population data sets (Abstract). Am J Hum Genet 69:A1886, 2001
O’Connell JR: Zero-recombinant haplotyping: applications to fine mapping using SNPs. Genet Epidemiol 19 (Suppl. 1):S64eCS70, 2000
O’Connell J: PROFILER: a program to compute cumulative probability profiles with applications to pedigree analysis (Abstract). Am J Hum Genet 67:A328, 2000
Antinozzi PA, Rothman JE: Using RNAi-based genetic screens to identify essential genes involved in pancreatic -cell function (Abstract). Diabetes 53:A31, 2004
Almasy L, Blangero J: Exploring positional candidate genes: linkage conditional on measured genotype. Behav Genet 34:173eC177, 2004
Rock MT, Brooks WH, Roszman TL: Calcium-dependent signaling pathways in T cells: potential role of calpain, protein tyrosine phosphatase 1b, and p130Cas in integrin-mediated signaling events. J Biol Chem 272:33377eC33383, 1997
Frangioni JV, Oda A, Smith M, Salzman EW, Neel BG: Calpain-catalyzed cleavage and subcellular relocation of protein phosphotyrosine phosphatase 1B (PTP-1B) in human platelets. EMBO J 12:4843eC4856, 1993
Dunning BE, Taborsky GJ Jr: Galanin: sympathetic neurotransmitter in endocrine pancreas Diabetes 37:1157eC1162, 1988
Clement S, Krause U, Desmedt F, Tanti JF, Behrends J, Pesesse X, Sasaki T, Penninger J, Doherty M, Malaisse W, Dumont JE, Le Marchand-Brustel Y, Erneux C, Hue L, Schurmans S: The lipid phosphatase SHIP2 controls insulin sensitivity. Nature 409:92eC97, 2001
Chan CB, De Leo D, Joseph JW, McQuaid TS, Ha XF, Xu F, Tsushima RG, Pennefather PS, Salapatek AM, Wheeler MB: Increased uncoupling protein-2 levels in -cells are associated with impaired glucose-stimulated insulin secretion: mechanism of action. Diabetes 50:1302eC1310, 2001(Nicholette D. Palmer, Car)
2 Center for Human Genomics, Wake Forest University School of Medicine, Winston-Salem, North Carolina
3 Department of Public Health Sciences, Wake Forest University School of Medicine, Winston-Salem, North Carolina
4 Department of Medicine, University of California at Los Angeles, Los Angeles, California
5 Department of Preventive Medicine and Biometrics, University of Colorado Health Sciences Center, Denver, Colorado
6 Department of Medicine, University of Texas Health Sciences Center at San Antonio, San Antonio, Texas
7 Medical Genetics Institute, Burns and Allen Cedars-Sinai Research Institute, Los Angeles, California
8 Department of Physiology and Biophysics, Keck School of Medicine, University of Southern California, Los Angeles, California
9 Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, North Carolina
AIR, acute insulin response; DI, disposition index; GEE, generalized estimating equation; IBD, identity-by-descent; IRASFS, Insulin Resistance Atherosclerosis Study Family Study; LOD, logarithm of odds; MAF, minor allele frequency; MLS, maximum LOD score; PC, pyruvate carboxylase; PYGM, muscle glycogen phophorylase; QTL, quantitative trait locus; SNP, single nucleotide polymorphism
ABSTRACT
Glucose homeostasis, a defining characteristic of physiological glucose metabolism, is the result of complex feedback relationships with both genetic and environmental determinants that influence insulin sensitivity and -cell function. Relatively little is known about the genetic basis of glucose homeostasis phenotypes or their relationship to risk of diabetes. Our group previously published a genome scan for glucose homeostasis traits in 284 African-American subjects from 21 pedigrees in the Insulin Resistance Atherosclerosis Study Family Study (IRASFS) and presented evidence for linkage to disposition index (DI) on chromosome 11q with a logarithm of odds (LOD) of 3.21 at 81 cM flanked by markers D11S2371 and D11S2002 (support interval from 71 to 96 cM). In this study, genotyping and analysis of an additional 214 African-American subjects in 21 pedigrees from the IRASFS yielded independent evidence of linkage to DI. When these two datasets were combined, a DI linkage peak was observed with an LOD of 3.89 at 78 cM (support interval from 67 to 89 cM). Fine mapping with 15 additional microsatellite markers in this 11q region for the entire 42 pedigrees resulted in an LOD score of 4.80 at 80 cM near marker D11S937 (support interval from 76 to 84 cM). In these 42 pedigrees, there was also suggestive evidence for linkage to acute insulin response (AIR) at two separate locations flanking the DI peak (64 cM, LOD 2.77, flanked by markers D11S4076 and D11S981; and 85 cM, LOD 2.54, flanked by markers D11S4172 and D11S2002). No evidence of linkage to the insulin sensitivity index (Si) was observed. Nine positional candidate genes were evaluated for association to DI and AIR. Among these candidates, single nucleotide polymorphisms (SNPs) in muscle glycogen phosphorylase showed evidence of association with DI (P < 0.011). In addition, SNPs in the pyruvate carboxylase gene showed evidence of association (P < 0.002) with AIR. Further analysis of these candidate genes, however, did not provide evidence that these SNPs accounted for the evidence of linkage to either DI or AIR. These detailed genetic analyses provide strong evidence of a DI locus on 11q in African-American pedigrees, with additional suggestive evidence of independent AIR loci in the same region.
It is widely accepted that type 2 diabetes is due to a combination of defects in insulin sensitivity and insulin secretion that result in impaired glucose homeostasis. Disposition index (DI) is a quantitative measure that describes the relationship between -cell sensitivity and insulin sensitivity (1). In metabolically normal individuals, changes in insulin sensitivity are accompanied by compensatory alterations in the -cell’s sensitivity to glucose. In practice, DI is measured as the product of the insulin sensitivity index (Si) and -cell function as measured by the acute insulin response (AIR). Despite being key features of glucose homeostasis, the molecular mechanism(s) by which insulin sensitivity and -cell function are coregulated to create a homeostatic environment are not well understood (1eC3).
Genetic analysis in appropriately phenotyped families is a potentially powerful approach for gaining insights into complex biological traits such as DI. To date, few mapping studies have incorporated the analysis of DI as a quantitative trait (4,5). Our group recently reported evidence for a DI locus between D11S2371 and D11S2002 on chromosome 11q at 81 cM with a logarithm of odds (LOD) score of 3.21 (support interval 71eC96 cM) in 284 African-American subjects from 21 pedigrees in the Insulin Resistance Atherosclerosis Study Family Study (IRASFS) (6). In addition to DI, other measures of glucose homeostasis have been previously reported to be linked to this region (7eC10).
The purpose of the current study is threefold. First, we report confirmatory evidence that DI and other measures of glucose homeostasis link to chromosome 11q in an independent and similarly phenotyped and ascertained sample of African-American pedigrees from the IRASFS. Second, we report fine-mapping quantitative trait locus (QTL) linkage analysis of glucose homeostasis measures for this region on 11q. Third, we report an evaluation of selected candidate genes within 11q and explore their contribution to the variation in measures of glucose homeostasis.
RESEARCH DESIGN AND METHODS
The study design, recruitment, and phenotyping in the IRASFS have been described in detail (11). Briefly, the IRASFS is a multicenter study designed to identify the genetic determinants of insulin resistance through quantitative, intermediate measures of glucose homeostasis. For this component of the study, members of large families of self-reported African-American descent from the Los Angeles, California, area were recruited and examined. Ascertainment of families was not based on phenotypes but instead on large family size. A clinical examination was performed that included interviews, a frequently sampled intravenous glucose tolerance test, anthropometric measurements, abdominal computed tomography scanning, resting blood pressure, blood collection, and spot urine. Specific to this report, measures of glucose homeostasis included those from the frequently sampled intravenous glucose tolerance test: Si, AIR, and DI.
Genotyping
Genome-wide scan and fine mapping.
The genotyping procedure has been previously described (6). Briefly, 383 microsatellite markers (average one marker/9.3 cM) were genotyped at the National Heart, Lung, and Blood Institute Mammalian Genotyping Service (Marshfield, WI; Marshfield set 10) in two independent sets of pedigrees, denoted as set 1 and set 2. Based on the evidence of linkage to 11q in set 1 (6) and the replication of the linkage in set 2 reported herein, an additional 15 microsatellite markers (fine-mapping markers) were genotyped on all African-American pedigrees (set 1 and set 2) within a 62-cM region centered on the linkage peak. Genotyping data from these 15 microsatellite markers were merged with the genome scan data for the fine-mapping analysis.
Candidate genes.
Nine positional candidate genes were chosen within the linked region based on their hypothesized functional role in glucose homeostasis. Initially, an average of five approximately equally spaced single nucleotide polymorphisms (SNPs) per positional candidate gene were selected from dbSNP (http://www.ncbi.nlm.nih.gov/projects/SNP/). SNPs with available frequency information were preferentially selected from the database. Following preliminary analysis, additional SNPs were chosen in positional candidate genes that showed association to the glucose homeostasis phenotypes DI, AIR, or Si. These additional SNPs were chosen to provide greater coverage of the linkage disequilibrium block and generally had a minor allele frequency (MAF) >0.20 in the IRASFS African-American pedigrees. Genotyping was performed on the Sequenom MassArray Genotyping System using methods previously described (12). Twenty-two individuals from two families served as blind duplicates (duplicated samples given new unique identifiers) to evaluate genotyping accuracy and exhibited a 0.48% discordance percentage.
Statistical analysis.
The statistical analyses reported are for the three primary purposes of this study: 1) confirmatory test that DI links to 11q in an independent and similarly phenotyped and ascertained replication sample of African-American pedigrees from the IRASFS, 2) fine mapping of both African-American samples on 11q, and 3) association analysis of selected candidate genes within the 11q region.
Genome-wide scan QTL linkage analysis.
The IRASFS genome scan was completed in two phases, denoted set 1 (284 African-American subjects from 21 pedigrees) and set 2 (214 African-American subjects from 21 pedigrees). The analysis of the first set of 21 African-American pedigrees showed evidence of linkage to 11q and has been previously reported (6). Since that publication, the locations of markers D11S2365 and ATA9B04 on chromosome 11 that were reported at 58 cM on the Marshfield map were changed to 55.96 and 56.27 cM, respectively, in the recent release of the human genome physical map (NCBI build 35, hg assembly 17). These shifts in location resulted in changes in marker order within the region, requiring these data to be reanalyzed due to the sensitivity of linkage analyses to marker order. The reanalysis for linkage of the first set of 21 pedigrees and the linkage analysis of the second set of 21 pedigrees estimated the identity-by-descent (IBD) statistics using maximum likelihood estimates of allele frequencies from the African-American pedigrees and the Monte Carlo approach implemented in the computer program LOKI (13,14).
For the analysis of set 1, set 2, and the combined dataset, quantitative traits were transformed to best approximate the distributional assumptions of the variance-component QTL linkage analysis (i.e., conditional normality after adjusting for covariates, homogeneity of variance). Thus, the results reported represent analyses on the square root of DI and AIR and the natural logarithm of Si after adding 1. For individuals with a self-reported diagnosis of diabetes or fasting glucose >126 mg/dl, the glucose homeostasis parameters of insulin sensitivity (Si), DI, and AIR were converted to missing. However, all available genotypic data were used to compute the IBD statistics using the Bayesian Markov Chain Monte Carlo approach implemented in the software LOKI (13,14). Previously, self-reported familial relationships were examined and modified using the genome scan data from both set 1 and set 2 pedigrees and the software PREST (15). Testing for evidence of linkage to a QTL was done using the variance-component approach implemented in the SOLAR software package (16). All LOD scores reported herein are empirical LOD scores determined within SOLAR by simulation (lodadj procedure). The entropy-based measure of linkage information content was computed using Merlin (17).
Fine mapping.
Map order and genetic distance for each of the fine-mapping markers were taken from the deCODE genetic map (18). Markers for fine mapping were analyzed using PedCheck (19) to identify any Mendelian inconsistencies. Ninety-nine total genotypes (from 79 individuals across 26 families) were converted to missing. The genotyping data from these 15 microsatellite markers were then merged with the combined genome scan genotyping data from set 1 and set 2, and QTL linkage analyses were computed as described above.
Candidate genes.
The genetic location of each positional candidate gene was calculated by extracting all markers from the deCODE genetic map (18) centered around the linkage region (D11S1357:D11S4104) and mapped to the chromosomal sex-averaged genetic positions to the latest human genome physical map positions (NCBI build 35, hg assembly 17). Using R version 1.9.1 (http://www.r-project.org/index.html), a smoothing spline was fit through these points and the spline used to interpolate genetic map positions corresponding to physical positions of gene loci. For each gene, we computed a single representative physical map position as the mean of the 5' end of the most telomeric exon in any known annotated gene splice variant and the 3' end of the most centromeric, assuming the genetic interval of the entire locus was small relative to that of the linkage region.
Initially, each of the SNPs evaluated in this work were examined for Mendelian inconsistencies in their genotypes using PedCheck (19). Any genotypes inconsistent with Mendelian inheritance were converted to missing (0.07%). Subsequently, maximum likelihood estimates of allele frequencies were computed using the largest set of unrelated individuals from the combined dataset and then tested for departures from Hardy-Weinberg proportions. To test for an association between individual SNPs and each quantitative phenotype (Si, AIR, and DI), a series of generalized estimating equations (GEEs; GEE1) was computed (20). The correlation between subjects within a pedigree was adjusted in the analyses by assuming an exchangeable correlation structure among siblings within a pedigree and computing the sandwich estimator of the variance (21). The sandwich estimator is also denoted the robust or empirical estimator of the variance, as it is robust to misspecification of the correlation matrix because it estimates the within-pedigree correlation matrix from moments of the data. For each phenotype (Si, AIR, and DI) and each SNP, the 2-df (degrees of freedom) test of genotypic association was performed. In addition, three individual contrasts defined by the a priori genetic models (dominant, additive, and recessive) were computed. The contrasts are defined by the genetic models (i.e., dominant model contrasts those with versus without the polymorphism, additive model tests for a dose effect in the number of alleles, and recessive model contrasts individuals homozygous for the polymorphisms versus those who are not) and are tested using a Wald test. If the overall genotypic association is significant, the a priori contrasts are examined directly. If the overall genotypic association is not significant, the a priori contrasts are examined after adjusting for the three comparisons using a Bonferroni adjustment. Tests reported here were computed adjusting for age, sex, and BMI. Subjects with type 2 diabetes were excluded from the analysis of glucose homeostasis traits.
Statistical significance for all tests was evaluated by a simulation study to obtain empiric P values. Specifically, for each SNP and each trait, the gene-dropping algorithm implemented in Mendel (22) was used and the analysis repeated 10,000 times. The empiric P value is the proportion of time that the simulated data had a P value greater than or equal to the observed P value from the real data. This approach allows the retention of the original correlation among pedigree members in the phenotypic and genotypic data.
The pairwise linkage disequilibrium (LD) statistics D' and r2 were computed, and a block of LD was defined based on marker-marker D' > 0.70 for all markers in the maximal set of unrelated individuals. Haplotypes were constructed within the LD block from SNPs with an MAF >10% in positional candidate genes that showed evidence of association to glucose homeostasis phenotypes. Analysis of these haplotypes was completed using a GEE1 analysis as described above, except the quasi-likelihood was weighted by the probability for each possible haplo-genotype for an individual. Specifically, we computed the expectation-maximization algorithm estimates for haplotype frequencies from the family data using the software ZAPLO (23). To reduce the complexity of the problem, ZAPLO assumes zero recombination between markers and generates all possible haplo-genotypes (i.e., haplotype combinations) for each individual. These haplo-genotypes are used as input into PROFILER (24) to compute the joint probability distribution for the haplotypes of individuals within the pedigree. This joint probability distribution is then used to estimate the probability of each haplotype pair combination (haplo-genotype) possible for each individual, conditional on the family data. Each individual enters into the GEE1 analysis once for each haplo-genotype possibility, weighted by the haplo-genotype probability. Thus, the weight for each individual sums to 1. The weighted GEE1 analyses were completed as above using the same transformations and sandwich estimator of the variance to account for the within-cluster correlation for each of the haplotypes observed with >5% frequency.
To test whether a particular SNP, either directly or through LD with the QTL, contributes to the evidence for linkage, each SNP under a dominant model was entered into the QTL linkage analysis and the change in the magnitude of the LOD score calculated. If the polymorphism directly or indirectly contributes to the evidence for linkage, the initial LOD score will be reduced in a model that includes the polymorphism as a covariate.
RESULTS
The descriptive characteristics of set 1 and set 2 participants are shown in Table 1. Diabetes-affected subjects (12.3%) were excluded from analysis of glucose homeostasis phenotypes. On average, the pedigrees in set 2 had fewer participants than pedigrees in set 1. In addition, set 2 was younger (41.5 ± 12.9), had a greater mean BMI (31.1 ± 6.9 kg/m2), and included more women (63.5%) than men.
Mapping of DI.
A summary of the linkage results for measures of glucose homeostasis on chromosome 11 is presented in online appendix Table 1 (available at http://diabetes.diabetesjournals.org). The variance-component linkage analysis for DI from SOLAR in set 1 was reanalyzed using the new marker order from the latest release of the human genome (NCBI build 35, hg assembly 17) and the LOKI (13,14) IBD calculation software. This reanalysis yielded a maximum LOD score (MLS) of 3.10 at 84 cM, flanked by markers D11S2371 and D11S2002, and a support interval from 72 to 93 cM. In set 2, the MLS was 1.26 at 102 cM, flanked by markers D11S2000 and D11S1986, and a support interval from 80 to 120 cM. In the combined set (set 1 + set 2), the MLS was 3.89 at 78 cM, flanked by markers D11S2371 and D11S2002, and a support interval from 67 to 89 cM.
With evidence of linkage of DI to 11q, fine mapping of the region was carried out by genotyping an additional 15 microsatellite markers across the region to give an average density of one marker/3 cM, with the largest gap being 6.94 cM. A QTL linkage analysis combining the genome scan markers and the 15 additional fine-mapping markers (within the 63.8-cM interval centered at the linkage peak) was completed in the entire African-American sample (Fig. 1A and online appendix Table 1). In set 1, the MLS for DI decreased slightly to 3.03 at 81 cM, flanked by markers D11S937 and D11S4172, and a support interval from 75 to 86 cM. In set 2, the MLS increased to 2.29 at 79 cM, flanked by markers D11S2371 and D11S937, and a support interval from 69 to 84 cM. Thus, the MLS maximized within 2 cM in these two independent samples, with set 2 replicating the QTL linkage effect identified in set 1. In the combined set, the MLS was increased to 4.80 at 80 cM at marker D11S937 with a support interval from 76 to 84 cM. Information content of the markers in this region was assessed using the program Merlin (17). Multipoint information content ranged from 0.83 to 0.95 across the fine-mapping region (59.2eC123.0 cM) (data not shown).
Mapping of AIR.
Variance-component linkage analyses of AIR and Si, components of DI, were also carried out with the entire marker set including the fine-mapping markers (Figs. 1B and online appendix Table 1). Linkage analysis of AIR showed suggestive evidence of linkage with two peaks flanking the DI peak in the combined data (set 1, set 2, and fine-mapping markers) analysis. The most significant linkage evidence for AIR was an MLS of 2.77 at 64 cM, flanked by markers D11S4076 and D11S981, and a support interval from 60 to 73 cM. In addition, a second MLS of 2.54 at 85 cM, flanked by markers D11S4172 and D11S2002, and a support interval from 59 to 91 cM was also observed. Si showed no significant evidence of linkage to this region.
Positional candidate gene analysis.
A total of nine positional candidate genes were evaluated within the linked region for association with glucose homeostasis traits. The location and possible function of these genes in modulation of glucose homeostasis is shown in Table 2. In the association analysis of SNPs genotyped within these genes (online appendix Table 2), the most significant evidence of association was observed with muscle glycogen phophorylase (PYGM) and pyruvate carboxylase (PC) precursor. SNPs in PYGM and PC were consistent with Hardy-Weinberg proportions with the exception of two SNPs in PYGM (rs565688, P = 0.009; rs477549, P = 0.025).
As seen in Table 3, multiple PYGM SNPs were significantly associated with DI from rs625172 through rs2073798 (all P < 0.032), with the exception of SNP rs686171, which trended toward association (P = 0.097). Simulations were carried out to calculate empiric P values for the individual SNPs. SNPs with an MAF >10% remained significantly associated with DI (empiric P < 0.041), with the exception of SNP rs565688 (empiric P = 0.068). Examination of LD in this region suggests one LD block encompassing SNPs rs625172 through rs2073798 (inter-SNP D' 0.96 over 10.7 kb) (online appendix Table 3), which is consistent with the association results. Results of the haplotypic association analysis of PYGM are summarized in Table 4. Three common haplotypes (>5%) were estimated within the LD block. Of these haplotypes, haplotype TAA, with a frequency of 29%, was significantly associated with DI (P = 0.020) following a dominant model of inheritance.
The association analysis results for PC with AIR are presented in Table 5. PC was significantly associated with AIR for SNPs rs896438 through rs7926787 (P < 0.002), with the exception of rs7932018 (P = 0.602) and rs2077432 (P = 0.235). Simulations provided evidence of association with AIR (empiric P < 0.014). Examination of LD in this region revealed one LD block encompassing SNPs rs7119426 through rs6591223 (inter-SNP D' 0.74 over 40.3 kb) (online appendix Table 4) and containing the trait-associated SNPs. Results of the haplotypic association analysis of PC are presented in Table 6. Five haplotypes with frequencies >5% were estimated within the LD block. Of these haplotypes, haplotype CACAGA, with a frequency of 20%, was significantly associated with AIR in all three genetic models, with the recessive model providing the strongest evidence of association (P < 0.001).
To assess the contribution to linkage, individual SNPs from all positional candidate genes were used as covariates in variance-component linkage analysis. For the 17 SNPs analyzed in PC and PYGM, the largest change in MLS occurred with SNP rs2167457 in PC with AIR. This SNP had an MAF of 0.24 and was significantly associated with AIR (empiric P = 0.004). Initially, the peak for AIR had an MLS of 2.76 at 64 cM. In our subset of participants genotyped for this SNP, the peak for AIR had an MLS of 2.29 at 64 cM. After adjustment for SNP rs216457, the largest LOD score change was a decrease to an MLS of 1.90 following a recessive model (data not shown). While this decrease in LOD suggests modest contribution of this SNP to the evidence of linkage, it does not reach statistical significance. Overall, when used as covariates in the QTL linkage analysis, single SNPs in the positional candidate genes did not significantly alter evidence of linkage.
DISCUSSION
We have carried out a detailed genetic analysis of glucose homeostasis traits on chromosome 11q in African-American families from the IRASFS. These QTL linkage analyses included genetic mapping in an independent and identically phenotyped and ascertained replication sample of pedigrees, fine mapping, and analysis of a set of relevant candidate genes for association with glucose homeostasis traits and their influence on the evidence for linkage. In two independent samples (set 1 and set 2), we observed evidence for linkage to DI maximizing within 2 cM near 80 cM in both samples. QTL linkage analyses of the entire 42 IRASFS African-American sample, including both the genome scan markers and 15 additional microsatellite markers, further delineated this region. QTL linkage analyses provide evidence of a DI locus on chromosome 11q at 80 cM (LOD 4.80) and suggestive evidence of linkage to AIR with peaks that appear to flank the DI linkage signal. The results raise the possibility that multiple loci within the region may be contributing to evidence of linkage with measures of glucose homeostasis.
The biologic and mathematical relationship among DI, AIR, and Si (DI = AIR x Si) and the coincident linkages of DI and AIR in this region suggested that we explore the possibility that pleiotrophy may exist within this region. We performed a series of bivariate QTL linkage analyses (AIR x Si, AIR x DI, and DI x Si) to test this hypothesis and did not observe either a meaningful increase or decrease in the respective LOD scores. It is important to recognize that the bivariate linkage and pleiotrophy tests do not account for multiple loci in linkage with each other. Thus, the results of the linkage test are inconclusive relative to pleiotrophy, and a more definitive test will be performed through detailed SNP mapping and association tests.
These results represent the first study to identify strong evidence of a QTL for DI. This region of chromosome 11 has, however, been linked to alternative measures of glucose homeostasis not measured in this study. Pratley et al. (10) reported evidence of linkage to 2-h postload insulin with an MLS of 1.31 at marker D11S2371 (82.6 cM) in Pima Indians. This marker is one that flanks our linkage peak for DI. In addition, the FUSION (Finland-United States Investigation of NIDDM Genetics) study reported preliminary findings of linkage to type 2 diabetes (MLS of 1.75 at 84 cM near D11S1314), fasting insulin (MLS of 2.07 at 89.5 cM), and fasting insulin-to-glucose ratio (MLS of 2.32 at 91.0 cM) in diabetic subjects. Evidence for replication of this linkage to type 2 diabetes with an MLS of 2.98 at 82 cM near D11S4172 was subsequently reported by these investigators (7). Most recently, the Framingham Offspring Study reported evidence of linkage to this region for fasting insulin (MLS of 2.43 at 85 cM near D11S2002) and homeostasis model assessment (MLS 2.01 at 92 cM near D11S2002 and D11S2000) (8).
Nine positional candidate genes were evaluated for contributions to the evidence of linkage. Analysis of PYGM and PC resulted in the most significant consistent evidence of association with DI and AIR, respectively. PYGM catalyzes and regulates the breakdown of glycogen to glucose-1-phosphate by removing the glucose residues from -(1,4)-linkages within glycogen molecules. PYGM was significantly associated with DI within a single LD block. Additionally, SNPs rs565688 and rs625172, which were significantly associated with DI, were also significantly associated with Si (rs565688, P = 0.001; rs625172, P = 0.017) (online appendix Table 2). Therefore, the evidence of association seen with DI could be driven by Si, as the two measures are mathematically and metabolically related (DI = Si x AIR). This may also explain the absence of a contribution to the DI linkage peak when these SNPs were used as covariates in the variance-component linkage analysis. Haplotype analysis of PYGM SNPs revealed a 29% frequency haplotype (TAA) significantly associated with DI (P = 0.020). This result is nominally significant, but there appears to be no apparent trend in the genotypic trait means associated with this haplotype.
PC is a nuclear-encoded mitochondrial enzyme that catalyzes the conversion of pyruvate to oxaloacetate in a tissue-specific manner; the initial reaction of glucose (liver, kidney) and lipid (adipose tissue, liver, brain) synthesis from pyruvate. PC was significantly associated with AIR in a single LD block. In haplotype analysis within the LD block, a 20% frequency haplotype (CACAGA) was significantly associated with AIR (P < 0.001). The trend in the genotypic trait means associated with this haplotype would suggest that its presence is associated with "risk" (i.e., lower AIR). A second less frequent (5%) haplotype (TGCACA) was also significantly associated with AIR (P = 0.001). It is interesting to note that Antinozzi and Rothman (25) have recently observed evidence that RNAi inhibition of PC gene expression substantially suppresses secretion of an insulin surrogate in an in vivo cell culture system, a role that is consistent with our evidence of genetic association with AIR. Unfortunately, as with PYGM, there appears to be no consistent trend in the genotypic trait means associated with these haplotypes.
To examine their contribution to the evidence of linkage, each SNP in the candidate genes was analyzed as a covariate in variance-component linkage analysis for traits Si, AIR, and DI. Individually, no single SNP significantly altered the linkage peaks. This result would preclude a direct contribution of these individual SNPs to the observed evidence of linkage. However, from these results alone, we cannot confidently rule out the possibility that these two positional candidate genes (PYGM and PC) are contributors to the overall evidence of linkage to this region on chromosome 11. Causal variants associated with the variance in DI and AIR could be located within the PYGM and PC gene regions, respectively. If the allele frequencies of the SNPs tested in these genes and the actual (unknown) causal variants differ substantially or are not in LD (e.g., low inter-SNP r2 values between the loci examined), it may not be possible to detect the effects of these loci (26). Additionally, the results from this study are consistent with the possibility that multiple genes, each with a potentially limited contribution to the trait, result in the observed linkage peak. It is important to also note that 100 known and hypothetical genes are located in the support interval (76eC84 cM) for DI. Further insights into the genes contributing to DI and AIR in this region may result from detailed SNP fine mapping currently in progress for this region.
ACKNOWLEDGMENTS
This research was supported in part by National Institutes of Health Grants HL060894 (to D.W.B.), HL060919 (to S.M.H.), HL060931 (to R.N.B.), HL060944 (L.E.W.), HL061019 (to J.M.N.), and HL061210 (to M.S.).
FOOTNOTES
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Bergman RN, Ader M, Huecking K, Van Citters G: Accurate assessment of -cell function: the hyperbolic correction. Diabetes 51 (Suppl. 1):S212eC220, 2002
Kahn SE, Prigeon RL, McCulloch DK, Boyko EJ, Bergman RN, Schwartz MW, Neifing JL, Ward WK, Beard JC, Palmer JP, Porte D Jr: Quantification of the relationship between insulin sensitivity and -cell function in human subjects: evidence for a hyperbolic function. Diabetes 42:1663eC1672, 1993
Bergman RN, Finegood DT, Kahn SE: The evolution of beta-cell dysfunction and insulin resistance in type 2 diabetes. Eur J Clin Invest 32 (Suppl. 3):35eC45, 2002
Watanabe RM, Ghosh S, Langefeld CD, Valle TT, Hauser ER, Magnuson VL, Mohlke KL, Silander K, Ally DS, Chines P, Blaschak-Harvan J, Douglas JA, Duren WL, Epstein MP, Fingerlin TE, Kaleta HS, Lange EM, Li C, McEachin RC, Stringham HM, Trager E, White PP, Balow J Jr, Birznieks G, Chang J, Eldridge W: The Finland-United States investigation of non-insulin-dependent diabetes mellitus genetics (FUSION) study. II. An autosomal genome scan for diabetes-related quantitative-trait loci. Am J Hum Genet 67:1186eC1200, 2000
An P, Hong Y, Weisnagel SJ, Rice T, Rankinen T, Leon AS, Skinner JS, Wilmore JH, Chagnon YC, Bergman RN, Bouchard C, Rao DC: Genomic scan of glucose and insulin metabolism phenotypes: the HERITAGE Family Study. Metabolism 52:246eC253, 2003
Rich SS, Bowden DW, Haffner SM, Norris JM, Saad MF, Mitchell BD, Rotter JI, Langefeld CD, Wagenknecht LE, Bergman RN: Identification of quantitative trait loci for glucose homeostasis: the Insulin Resistance Atherosclerosis Study (IRAS) Family Study. Diabetes 53:1866eC1875, 2004
Silander K, Scott LJ, Valle TT, Mohlke KL, Stringham HM, Wiles KR, Duren WL, Doheny KF, Pugh EW, Chines P, Narisu N, White PP, Fingerlin TE, Jackson AU, Li C, Ghosh S, Magnuson VL, Colby K, Erdos MR, Hill JE, Hollstein P, Humphreys KM, Kasad RA, Lambert J, Lazaridis KN, Lin G, Morales-Mena A, Patzkowski K, Pfahl C, Porter R, Rha D, Segal L, Suh YD, Tovar J, Unni A, Welch C, Douglas JA, Epstein MP, Hauser ER, Hagopian W, Buchanan TA, Watanabe RM, Bergman RN, Tuomilehto J, Collins FS, Boehnke M: A large set of Finnish affected sibling pair families with type 2 diabetes suggests susceptibility loci on chromosomes 6, 11, and 14. Diabetes 53:821eC829, 2004
Panhuysen CI, Cupples LA, Wilson PW, Herbert AG, Myers RH, Meigs JB: A genome scan for loci linked to quantitative insulin traits in persons without diabetes: the Framingham Offspring Study. Diabetologia 46:579eC587, 2003
Ghosh S, Watanabe RM, Valle TT, Hauser ER, Magnuson VL, Langefeld CD, Ally DS, Mohlke KL, Silander K, Kohtamaki K, Chines P, Balow J Jr, Birznieks G, Chang J, Eldridge W, Erdos MR, Karanjawala ZE, Knapp JI, Kudelko K, Martin C, Morales-Mena A, Musick A, Musick T, Pfahl C, Porter R, Rayman JB: The Finland-United States investigation of non-insulin-dependent diabetes mellitus genetics (FUSION) study. I. An autosomal genome scan for genes that predispose to type 2 diabetes. Am J Hum Genet 67:1174eC1185, 2000
Pratley RE, Thompson DB, Prochazka M, Baier L, Mott D, Ravussin E, Sakul H, Ehm MG, Burns DK, Foroud T, Garvey WT, Hanson RL, Knowler WC, Bennett PH, Bogardus C: An autosomal genomic scan for loci linked to prediabetic phenotypes in Pima Indians. J Clin Invest 101:1757eC1764, 1998
Henkin L, Bergman RN, Bowden DW, Ellsworth DL, Haffner SM, Langefeld CD, Mitchell BD, Norris JM, Rewers M, Saad MF, Stamm E, Wagenknecht LE, Rich SS: Genetic epidemiology of insulin resistance and visceral adiposity: the IRAS Family Study design and methods. Ann Epidemiol 13:211eC217, 2003
Buetow KH, Edmonson M, MacDonald R, Clifford R, Yip P, Kelley J, Little DP, Strausberg R, Koester H, Cantor CR, Braun A: High-throughput development and characterization of a genomewide collection of gene-based single nucleotide polymorphism markers by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Proc Natl Acad Sci U S A 98:581eC584, 2001
Heath SC, Snow GL, Thompson EA, Tseng C, Wijsman EM: MCMC segregation and linkage analysis. Genet Epidemiol 14:1011eC1016, 1997
Heath SC: Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet 61:748eC760, 1997
McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet 66:1076eC1094, 2000
Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet 62:1198eC1211, 1998
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30:97eC101, 2002
Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, Richardsson B, Sigurdardottir S, Barnard J, Hallbeck B, Masson G, Shlien A, Palsson ST, Frigge ML, Thorgeirsson TE, Gulcher JR, Stefansson K: A high-resolution recombination map of the human genome. Nat Genet 31:241eC247, 2002
O’Connell JR, Weeks DE: PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 63:259eC266, 1998
Zeger SL, Liang KY: Longitudinal data analysis for discrete and continuous outcomes. Biometrics 42:121eC130, 1986
Hardin J, Hilbe J: Generalized Estimating Equations. New York, Chapman & Hall/CRC, 2003
Lange KCR, Horvath S, Perola M, Sabatti C, Sinsheimer J, Sobel E: Mendel version 4.0: a complete package for the exact genetic analysis of discrete traits in pedigree and population data sets (Abstract). Am J Hum Genet 69:A1886, 2001
O’Connell JR: Zero-recombinant haplotyping: applications to fine mapping using SNPs. Genet Epidemiol 19 (Suppl. 1):S64eCS70, 2000
O’Connell J: PROFILER: a program to compute cumulative probability profiles with applications to pedigree analysis (Abstract). Am J Hum Genet 67:A328, 2000
Antinozzi PA, Rothman JE: Using RNAi-based genetic screens to identify essential genes involved in pancreatic -cell function (Abstract). Diabetes 53:A31, 2004
Almasy L, Blangero J: Exploring positional candidate genes: linkage conditional on measured genotype. Behav Genet 34:173eC177, 2004
Rock MT, Brooks WH, Roszman TL: Calcium-dependent signaling pathways in T cells: potential role of calpain, protein tyrosine phosphatase 1b, and p130Cas in integrin-mediated signaling events. J Biol Chem 272:33377eC33383, 1997
Frangioni JV, Oda A, Smith M, Salzman EW, Neel BG: Calpain-catalyzed cleavage and subcellular relocation of protein phosphotyrosine phosphatase 1B (PTP-1B) in human platelets. EMBO J 12:4843eC4856, 1993
Dunning BE, Taborsky GJ Jr: Galanin: sympathetic neurotransmitter in endocrine pancreas Diabetes 37:1157eC1162, 1988
Clement S, Krause U, Desmedt F, Tanti JF, Behrends J, Pesesse X, Sasaki T, Penninger J, Doherty M, Malaisse W, Dumont JE, Le Marchand-Brustel Y, Erneux C, Hue L, Schurmans S: The lipid phosphatase SHIP2 controls insulin sensitivity. Nature 409:92eC97, 2001
Chan CB, De Leo D, Joseph JW, McQuaid TS, Ha XF, Xu F, Tsushima RG, Pennefather PS, Salapatek AM, Wheeler MB: Increased uncoupling protein-2 levels in -cells are associated with impaired glucose-stimulated insulin secretion: mechanism of action. Diabetes 50:1302eC1310, 2001(Nicholette D. Palmer, Car)