Polymorphism in the Calsequestrin 1 (CASQ1) Gene on Chromosome 1q21 Is Associated With Type 2 Diabetes in the Old Order Amish Mao Fu, 1,2 Coleen M. Damcott, 1 Mona Sabra, 1 Toni I. Pollin, 1 Sandra H. Ott, 1 Jian Wang, 1 Michael J. Garant, 1 Jeffrey R. O’Connell, 1 Braxton D. Mitchell, 1 and Alan R. Shuldiner 1,3 Calsequestrin (CASQ)1 is involved in intracellular stor- age and release of calcium, a process that has been shown to mediate glucose transport in muscle. Its gene, CASQ1, is encoded on chromosome 1q21, a region that has been linked to type 2 diabetes in the Amish and several other populations. We screened all 11 exons, exon-intron junctions, and the proximal regulatory re- gion of CASQ1 for mutations. We detected four novel single nucleotide polymorphisms (SNPs) (�1470C3T, �1456delG, �1366insG, and 593C3T). Ten informative SNPs within CASQ1 were genotyped in Amish subjects with type 2 diabetes (n � 145), impaired glucose toler- ance (n � 148), and normal glucose tolerance (n � 358). Rs2275703 and rs617698 in introns 4 and 2 were signif- icantly associated with type 2 diabetes (P � 0.008 and 0.04, respectively); three other SNPs showed borderline evidence for association to type 2 diabetes (P � 0.076 – 0.093). Furthermore, in nondiabetic subjects (n � 754), both rs2275703 and rs617698 were significantly associ- ated with glucose area under the curve during an oral glucose tolerance test (P � 0.035 and 0.013, respec- tively). Haplotype analysis suggested that no haplotype could explain these associations better than rs2275703. These findings, coupled with similar findings in Utah Caucasians, suggest that sequence variation in CASQ1 may influence risk of type 2 diabetes. Diabetes 53: 3292–3299, 2004 T ype 2 diabetes is a common heterogeneous dis- order in which genetic and nongenetic factors (i.e., excess caloric intake and physical inactiv- ity) potently influence risk (1,2). Genome-wide scans have led to the chromosomal localization of suscep- tibility loci for type 2 diabetes, with the most promising and replicated findings on chromosomes 1q21-q24 (3–10), 2q37 (11), 12q24 (12–14), and chromosome 20 (15–18). Consistent with current thinking, multiple linked regions suggest that several susceptibility genes contribute to the development of type 2 diabetes (2). To date, positional cloning has identified variants in calpain 10 (CAPN10) on chromosome 2q37 (19) and hepatocyte nuclear factor 4� (HNF4A) on chromosome 20q12-q13.1 (20,21) as type 2 diabetes susceptibility genes in some populations. We previously reported evidence for linkage of type 2 diabetes and impaired glucose homeostasis to a locus on chromosome 1q21-q24 (logarithm of odds 2.35, P � 0.0008) in the Old Order Amish (3). Linkage of this region on chromosome 1 to type 2 diabetes has been reported in at least five other populations, including Pima Indians (4), Utah Caucasians (5), French Caucasians (6), U.K. Cauca- sians (7), and Chinese (8); linkages of diabetes-related traits to this region have also been reported (9,10). Such robust replication provides strong evidence for a shared type 2 diabetes susceptibility locus in this chromosomal interval. However, this region contains at least 450 genes that include several excellent positional candidate genes, introducing the possibility that more than one gene variant in the region may influence diabetes risk. A number of positional candidate genes within this region have been studied to date, with variation in potassium inwardly rectifying channel, subfamily J, members 9 and 10 (KCNJ9/10) (22,23), liver-specific pyruvate kinase (PKLR) (24), C-reactive protein (CRP) (25), lamin A/C (LMNA) (26), and omentin (ITLN) (27) showing modest evidence for association with type 2 diabetes or related traits. To further refine this interval with the aim of position- ally cloning the type 2 diabetes gene(s) in this region, we embarked upon linkage disequilibrium (LD) mapping by performing genotype and association analysis of single nucleotide polymorphisms (SNPs) in Amish type 2 dia- betic cases and nondiabetic control subjects. To date, �250 SNPs over the 28.5-Mb region from 148.5 to 170 Mb (NCBI [National Cancer for Biotechnology] build 34.3) From the 1Division of Endocrinology, Diabetes and Nutrition, University of Maryland School of Medicine, Baltimore, Maryland; the 2Division of Endocri- nology, Second Affiliated Hospital, Sun Yat-sen University, Guangzhou, China; and the 3Geriatric Research and Education Clinical Center, Veterans Admin- istration Medical Center, Baltimore, Maryland. Address correspondence and reprint requests to Alan R. Shuldiner, MD, Division of Endocrinology, Diabetes and Nutrition, University of Maryland School of Medicine, 660 West Redwood St., Room 494, Baltimore, MD 21201. E-mail: ashuldin@medicine.umaryland.edu. Received for publication 25 May 2004 and accepted in revised form 9 September 2004. Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. AFDS, Amish Family Diabetes Study; AUC, area under the curve; CASQ, calsequestrin; HOMA, homeostasis model assessment; HOMA-IR, HOMA of insulin resistance; IGT, impaired glucose tolerance; LD, linkage disequilib- rium; NGT, normal glucose tolerance; OGTT, oral glucose tolerance test; SNP, single nucleotide polymorphism. © 2004 by the American Diabetes Association. 3292 DIABETES, VOL. 53, DECEMBER 2004 have been genotyped. Multiple SNPs associated with type 2 diabetes were identified within the interval between �157 and 158.5 Mb, including rs617698, located within the calsequestrin 1 (CASQ1) gene. CASQ1 is a calcium storage protein localized in the sarcoplasmic reticulum of fast- twitch skeletal muscle cells. It is a constituent of the protein backbone of the luminal Ca2� reservoir and plays a pivotal role in calcium flux through regulation of calcium channel activity and interaction with other proteins present in the sarcoplasmic reticulum (28). Calcium re- lease from the sarcoplasmic reticulum into the cytosol has been shown to regulate GLUT4 expression (29) and glu- cose transport in muscle (30). Furthermore, Howarth et al. (31) reported that CASQ1 expression and calcium binding is increased in streptozotocin-induced diabetic rat skeletal muscle. We hypothesized that variation in CASQ1 may affect insulin action on glucose uptake and glycogen synthesis in skeletal muscle and thus type 2 diabetes susceptibility. To examine this positional candidate gene further, we screened CASQ1 for mutations and deter- mined whether the observed sequence variation was asso- ciated with type 2 diabetes and related traits in the Old Order Amish. RESEARCH DESIGN AND METHODS The Old Order Amish are a genetically well-defined Caucasian founder population who live in a relatively homogeneous environment and have large sibships. Nearly all of these individuals share common ancestors; �95% of the chromosomal material of the Amish community of Lancaster County (now numbering �30,000 individuals) appears to have descended from �100 founders who emigrated to this area in the mid 1700s (12–14 generations) (32,33). The Amish Family Diabetes Study (AFDS) was begun in 1995 to identify susceptibility genes for type 2 diabetes and related traits. Probands with diabetes onset between 35 and 65 years and all willing first- and second-degree family members �18 years of age were invited to participate. DNA from 18 non–first-degree relatives with type 2 diabetes (36 alleles) from families providing evidence for linkage of diabetes to chromosome 1q21-q24 was used for sequence analysis of CASQ1. To define LD and haplotype structure of CASQ1, all SNPs were genotyped in a limited set of 258 Amish subjects chosen to be relatively unrelated (non–first-degree relatives) to each other. Next, to examine the relationship between haplotype-tagging polymor- phisms in CASQ1 and diabetes, we genotyped DNA from subjects with type 2 diabetes (n � 145), impaired glucose tolerance (IGT) (n � 148), and normal glucose tolerance (NGT) (n � 358) for case-control association analysis. NGT subjects included in this analysis were required to be at least 38 years of age at the time of examination. For association analysis of quantitative traits (e.g., glucose and insulin), a larger set of 754 nondiabetic subjects (including the 506 NGT and IGT subjects described above) was studied. Informed consent was obtained from all AFDS participants, and the study protocol was approved by the institutional review board at the University of Maryland School of Medicine. Trait measurements. Details of methods for phenotyping of AFDS subjects have been previously described (3,34). Briefly, body height and weight were measured in all subjects using standard protocols. After a 12-h overnight fast, a 75-g oral glucose tolerance test (OGTT) was administered, with venous blood samples obtained at 0, 30, 60, 90, 120, 150, and 180 min for glucose and insulin measurements. The total glucose and insulin areas under the curve (AUCs) during the 3-h OGTT were calculated using the trapezoid method. Insulin resistance was estimated from fasting blood glucose and insulin levels with the homeostasis model assessment (HOMA) (HOMA of insulin resistance [HOMA-IR] � [fasting insulin {�U/ml} � fasting glucose {mmol/l}]/22.5). Diabetes was defined by fasting plasma glucose level (�7 mmol/l), 2-h OGTT plasma glucose level (�11.1 mmol/l), random plasma glucose level (�11.1 mmol/l), the use of insulin or oral glucose-lowering agents, or a diagnosis of diabetes, with age of onset between 35 and 65 years, documented by a physician. IGT was diagnosed based on OGTT plasma glucose levels (2-h OGTT plasma glucose level between 7.8 and 11.1 mmol/l). NGT was defined based on fasting plasma glucose level (�6.1 mmol/l) and 2-h OGTT plasma glucose level (�7.8 mmol/l). Mutation screening. Primers were designed to specifically amplify CASQ1 exons, splice junctions, and 1,000 bp of the 5 flanking sequence using primer T A B L E 1 F re qu en ci es o f C A S Q 1 S N P s in su b je ct s w it h ty p e 2 d ia b et es , IG T , an d N G T M in o r al le le fr eq u en cy T yp e 2 d ia b et es vs . N G T * T yp e 2 d ia b et es � IG T vs . N G T * S N P n am e L o ca ti o n P o si ti o n † M aj o r/ m in o r al el e T yp e 2 d ia b et es (n � 14 5) IG T (n � 14 8) N G T (n � 35 8) O R (9 5% C I) P O R (9 5% C I) P C as qS N P 1 5’ fl an k 14 70 C /T 0. 47 1 0. 46 3 0. 42 1 1. 33 (0 .8 9– 1. 98 ) 0. 15 1. 24 (0 .9 2– 1. 68 ) 0. 17 R s1 18 66 94 5’ fl an k 14 04 C /T 0. 26 5 0. 25 0 0. 26 3 1. 14 (0 .7 8– 1. 66 ) 0. 47 0. 98 (0 .7 4– 1. 32 ) 0. 92 R s3 83 82 16 In tr o n 1 � 10 74 /g gc at tc ag at ag gc ct 0. 09 2 0. 10 1 0. 13 9 0. 73 (0 .4 1– 1. 31 ) 0. 31 0. 70 (0 .4 4– 1. 08 ) 0. 09 1 R s6 17 69 8 In tr o n 2 � 23 12 A /G 0. 27 9 0. 32 7 0. 35 5 0. 65 (0 .4 2– 0. 99 ) 0 .0 4 0 0. 72 (0 .5 2– 1. 00 ) 0 .0 4 7 R s6 17 59 9 In tr o n 2 � 23 99 G /A 0. 28 1 0. 33 9 0. 36 8 0. 58 (0 .3 2– 1. 07 ) 0. 07 6 0. 77 (0 .4 7– 1. 26 ) 0. 42 R s3 74 76 22 In tr o n 2 � 28 89 G /C 0. 07 1 0. 08 7 0. 12 4 0. 64 (0 .3 6– 1. 15 ) 0. 18 0. 56 (0 .3 3– 0. 94 ) 0 .0 2 6 C as qS N P 4 E x o n 3 � 30 59 C /T 0. 07 1 0. 05 9 0. 03 0 1. 94 (0 .8 7– 4. 31 ) 0. 09 3 1. 82 (0 .9 3– 3. 59 ) 0. 06 4 R s2 27 57 03 In tr o n 4 � 45 35 C /A 0. 31 2 0. 37 1 0. 42 6 0. 57 (0 .4 1– 0. 79 ) 0 .0 0 8 0. 64 (0 .4 7– 0. 85 ) 0 .0 0 5 R s8 22 45 0 In tr o n 7 � 78 08 G /A 0. 28 0 0. 31 1 0. 34 1 0. 70 (0 .4 5– 1. 07 ) 0. 09 2 0. 75 (0 .5 0– 1. 14 ) 0. 10 R s3 82 75 32 In tr o n 10 � 10 17 1 G /A 0. 03 4 0. 02 6 0. 06 7 0. 35 (0 .1 0– 1. 32 ) 0. 13 0. 26 (0 .0 8– 0. 88 ) 0 .0 3 0 *P va lu es ar e b as ed o n ge n o ty p e fr eq u en ci es an d O R s re fl ec t th e o d d s o f d is ea se as so ci at ed w it h h av in g tw o co p ie s o f th e m in o r al le le ve rs u s th e o d d s o f d is ea se as so ci at ed w it h h av in g tw o co p ie s o f th e m aj o r al le le . F o r rs 22 75 70 3 an d rs 61 76 98 , th e O R s fo r ty p e 2 d ia b et es o f th e m aj o r al le le s ar e 1. 75 (9 5% C I 1. 27 –2 .4 3) an d 1. 54 (1 .0 1– 2. 38 ), re sp ec ti ve ly . P va lu es � 0. 05 ar e sh o w n in b o ld . †T h e n u cl eo ti d e p o si ti o n o f ea ch p o ly m o rp h is m is co u n te d fr o m th e A o f th e A T G st ar t co d o n , w h ic h w as d es ig n at ed as p o si ti o n � 1. M. FU AND ASSOCIATES DIABETES, VOL. 53, DECEMBER 2004 3293 3 (available at http//www-genome.wi.mit.edu/cgi-bin/primer/primer3_www. cgi) (see Table 1 of online appendix [available at http://diabetes.diabetesjour- nals.org]). PCR conditions were 95°C for 2 min, followed by 35 cycles of 95°C for 45 s, 56°C for 45 s, and 72°C for 1 min, with a final extension at 72°C for 7 min. Both strands were sequenced from DNA of the 18 subjects on an ABI 3700 DNA sequencer and analyzed with Sequence Analysis 3.2 software (Applied Biosystems Division, PerkinElmer, Foster City, CA). Genotyping. The 593C3T polymorphism in exon 3 (casqSNP4) was ampli- fied and detected by sequence analysis on an ABI 3700 DNA sequencer (A of ATG start codon of Genbank accession no. NM_001231, designated as position �1). CasqSNP1 ( 1470 C3T), rs1186694, rs617698, and rs617599 were typed using the AcycloPrime-FP SNP detection system (PerkinElmer Life Sciences, Boston, MA), which determines the base at a SNP location by template- directed dye-terminator incorporation followed by detection with fluores- cence polarization (35). Rs3747622, rs3747623, rs2275703, rs822450, rs2275705, and rs686015 were typed using the SNPstream UHT Genotyping System (Beckman Coulter, Fullerton, CA) (36). SNP 1456delG (CasqSNP2), rs3838216, and rs3827532 were genotyped by pyrosequencing using a PSQ HS 96A Pyrosequencer according to the manufacturer’s protocol (Pyrosequenc- ing, Uppsala, Sweden). In addition to SNPs within CASQ1, we also genotyped SNPs flanking the gene to determine whether associated SNPs within CASQ1 were in LD with extragenic SNPs. These included rs860294, rs1408664, rs881291, rs625129, rs7545350, rs6427504, and rs2820553, which are located 1.9 –145.7 kb upstream of CASQ1, and rs680083, rs1041066, and rs2147472, which are located 8.9 –169.9 kb downstream of CASQ1, all performed by SNPstream UHT genotyping. All typing included at least 5% duplicate samples to determine mistyping rates, which were 0 –2%. Statistical analysis. Before analysis, all SNPs were checked for Mendel errors using the PedCheck software program (37) in the extended Amish pedigree. A small number of Mendel errors were either resolved or genotypes removed before analysis. Allele frequencies of polymorphisms were deter- mined by gene counting. Because of the extensive relatedness of study subjects in our Amish pedigrees, we tested for associations of SNP genotypes while accounting for the pedigree structure. Briefly, we used a variance component approach in which we modeled the probability that the subject was a case or control as a function of the individual’s age, sex, and genotype, conditional on the correlations in phenotype among relative pairs. For the primary analysis, we considered an additive genetic model in which the genotype was coded as 0, 1, or 2, depending on whether the subject was homozygous for the minor allele (genotype � 2), heterozygous (genotype � 1), or homozygous for the major allele (genotype � 0). Statistical testing was performed by testing the likelihood of the data under a model in which the genotype effect was estimated against the likelihood of a nested model in which the genotype effect was constrained to be zero. Under the likelihood ratio test, two times the difference in minus log10 base likelihoods is distributed as a �2 statistic with degrees of freedom equal to the difference in number of parameters between the models (in this instance, one). Secondary analyses were carried out assuming dominant and recessive genetic models, in which instance the genotype was coded as 0 or 1 depending on whether the subject carried the minor allele (genotype � 1 if minor allele carrier, 0 if not [dominant model]; genotype � 1 if homozygous for minor allele, 0 if not, [recessive model]). Parameter estimates (i.e., �-coefficients) were obtained by maximum likelihood, and odds ratios (ORs) were obtained by taking the inverse log of the �-coefficient. The OR for the additive model was scaled to reflect the odds that a case was homozygous for the minor allele versus the odds that the case was homozygous for the major allele. The variance components analysis was carried out using the SOLAR software program (38). In similar fashion, we also evaluated the effect of genotype on levels of quantitative traits (e.g., glucose and insulin levels) by comparing mean trait levels across genotypes. We compared the likelihood of a model in which the trait values were allowed to vary by genotype (unconstrained model) with that in which the genotype effects were constrained to be zero. As before, the model likelihoods and parameter estimates were computed conditional on the pedigree structure, and likelihoods obtained from competing models were compared by likelihoods ratio test. Again, the additive model was used as our primary analysis, and dominant and recessive models were used as secondary analyses. The quantitative trait analyses were conducted in nondiabetic (NGT and IGT) individuals only. The disequilibrium statistics r2 and �D � between each pair of SNPs were calculated using ZAPLO (39). For SNP pairs with r2 � 0.9, only one of the two SNPs was genotyped in the full sample set. To determine the haplotype structure of SNPs with CASQ1, we used Haploview (40). We used Haploscore (41), which uses efficient score statistics within a general linear model framework, to test for haplotype associations between case and control subjects. Haploscore computes empirical P values by simulation. For all analyses, P values �0.05 were considered statistically significant, while P values between 0.05 and 0.1 were considered of borderline statistical signifi- cance. RESULTS Sequencing of all exons, exon-intron boundaries, and the proximal regulatory region of CASQ1 revealed a total of six polymorphisms. These include three novel SNPs up- stream of the ATG start codon, 1470T3C (CasqSNP1), 1456delG (CasqSNP2), and 1366insG (CasqSNP3), and one novel SNP in exon 3, 593C3T (CasqSNP4), which was silent (Genbank accession no. NM_001231) (Fig. 1). In addition, we detected two SNPs previously reported in the dbSNP database (rs1186694 and rs3747623), located at 1,404 bp upstream of the ATG start codon and in intron 2, respectively. Of the 22 SNPs in and immediately flanking this gene reported in the public databases at the time of this study (including rs1186694 and rs3747623), we con- firmed 11 by genotype analysis; the other 11 were not polymorphic (or had very low allele frequency) in the Amish (Fig. 1). In total, there were 15 SNPs in CASQ1; one occurred in the coding region and did not predict alter- FIG. 1. Schematic of CASQ1 and SNPs. Exons are shown as boxes and horizontal lines as introns, regulatory region or flanking segments. Novel SNPs discovered by sequence analysis are underlined, and dbSNPs are designated by the “Rs ” numbers; those that were informative are above, while those that were not polymorphic in the Amish are below. SNPs with asterisks next to them were designated as informative haplotype- tagging polymorphisms and genotyped in the full sample set (see text for details). Note exon and intron sizes are not strictly drawn to scale. CASQ1 AND TYPE 2 DIABETES IN THE AMISH 3294 DIABETES, VOL. 53, DECEMBER 2004 ations in the amino acid sequence, four were in the 5 flanking region, one was in the 3 flanking region, and all others were in introns and did not predict any obvious alterations in RNA splicing. Based on our initial genotype and LD analysis per- formed in 258 relatively unrelated Amish subjects, CasqSNP1 and CasqSNP2, located only 14 bp apart in the 5 flanking region, were in complete LD with one another; therefore, only CasqSNP1 was genotyped in the full sam- ple. Similarly, strong LD was observed between rs617698 and rs374723 (r2 � 0.94) and rs822450 and rs686015 (r2 � 0.93); thus, only rs617698 and rs822450 are reported in the full sample. In addition, CasqSNP3 and rs2275705 were found to be rare (allele frequencies 0.034 and 0.002, respectively) in the Amish and not pursued further. We typed the remaining 10 informative haplotype-tagging polymorphisms in 145 type 2 diabetic individuals, 148 individuals with IGT, and 358 individuals with NGT and tested for association between CASQ1 polymorphism and diabetes. Using an additive model, genotype frequencies of rs2275703 and rs617698 in introns 4 and 2 of CASQ1 differed significantly between subjects with type 2 diabe- tes and those with NGT (P � 0.008 and 0.04, respectively), with the more common allele being the at-risk allele for type 2 diabetes (OR for common alleles 1.75 [95% CI 1.27–2.43] and 1.54 [1.01–2.38], respectively) (Table 1). Three other SNPs spanning from intron 2 to the immediate 3 flanking region showed borderline evidence for associ- ation to type 2 diabetes (P � 0.076 – 0.093). Rs2275703 was also associated with the combined trait type 2 diabetes/ IGT (OR for common allele 1.56 [1.18 –2.23], P � 0.005), as were the common alleles of rs3747622 (1.78 [1.06 –3.03], P � 0.026) and rs617698 (1.39 [1.00 –1.92], P � 0.047), both in intron 2, and rs3827532 in intron 10 (3.85 [1.14 –12.50], P � 0.030). Three other polymorphisms showed borderline association with type 2 diabetes/IGT (P � 0.064 – 0.10). Genotypic associations estimated under recessive and dominant genetic models did not provide a better fit to the data than that estimated under the additive model (data not shown). Eight SNPs lying 1.9 –217 kb upstream and three SNPs lying 8.9 –169.9 kb downstream of CASQ1 were not significantly associated with type 2 diabetes (data not shown). We assessed the relationship between the SNPs and quantitative traits in an expanded set of 754 nondiabetic members of the AFDS. Three SNPs (rs617698, rs2275703, and rs822450) were significantly associated with glucose AUC during the OGTT (P � 0.05) (Fig. 2). These findings in nondiabetic subjects provide additional support for an effect of CASQ1 polymorphism on glucose homeostasis. FIG. 2. Association of CASQ1 SNPs rs617698 and rs2275703 with glucose and insulin levels during a 3-h OGTT in nondiabetic individuals. Both SNPs were associated with total OGTT glucose AUC (P � 0.013 and 0.035, addidtive model for rs617698 and rs2275703, respectively). The A allele of rs617698 appeared to better fit a recessive genetic model, with AA homozygotes having significantly higher glucose AUC than heterozygotes and GG homozygotes (P � 0.013). By contrast, the C allele of 2275703 appeared to better fit a dominant genetic model, with CC homozygotes and heteroyzgotes having higher glucose AUC than AA homozygotes (P � 0.03). *P < 0.05 between genotypes for a given time point. Rs822450 was also significantly associated with glucose AUC (P � 0.042) (data not shown). There was no evidence for association of rs617698 or rs2275703 with OGTT insulin levels. M. FU AND ASSOCIATES DIABETES, VOL. 53, DECEMBER 2004 3295 There was no evidence for association of any of the SNPs with BMI, insulin levels (Fig. 2), or HOMA-IR. The two SNPs associated with type 2 diabetes, rs617698 and rs2275703, are in introns 2 and 4, and are located only 2,223 bp apart. We thus sought to further define the haplotype structure of CASQ1 to determine whether these two SNPs might be tagging a single haplotype that is associated with type 2 diabetes or whether they may be defining distinct haplotypes, both of which are associated with type 2 diabetes. Using the Hapview program (40), the CASQ1 region formed a complex pattern of haplotype structure. There were four main blocks (or bins): bin 1 (rs881291, rs1408664, and rs860294; promoter), bin 2 (CasqSNP1 and rs1186694; promoter), bin 3 (rs3838216, rs617698, rs617599, rs3747622, and rs3747623; intron 1–2), and bin 4 (rs822450, rs3827532, 686015, rs680083, and rs1041066; intron 7 to 3 flank). Rs2275703 in intron 4 and several SNPs in bin 4 (e.g., rs822450, rs3827532, rs686015, and rs680083 from intron 7 to the 3 flanking region) appeared to be in partial LD with bin 3 SNPs (Fig. 3) (pairwise r2 and �D � for all SNP pairs also reported in Table 2 in the online appendix). Although the two SNPs most significantly associated with type 2 diabetes, rs617698 and rs2275703, are not assigned to the same bin, they are moderately correlated (�D � � 0.78, r2 � 0.49). In fact, the borderline-significant P value of 0.04 for association of rs617698 with type 2 diabetes can be accounted for by this correlation by examining the haplotype analysis results in Table 2. These results indicate that the more common C allele at rs2275703 is associated with increased risk of diabetes (as indicated by a positive Haploscore) and that the A allele at rs617698 is also associated with increased risk primarily because it tracks with the rs2275703 C allele. In other words, the rs617698 at-risk A allele association with type 2 diabetes is derived from its association with the common A–C haplotype containing the rs2275703 at-risk C allele. Furthermore, examination of three-, four-, and five-SNP haplotypes did not provide evidence for association with FIG. 3. Haplotype structure of CASQ1 and flanking regions. Shown is the haplotype structure of CASQ1, as estimated by Haploview (40), on genotypes from 258 Amish subjects. �D�� values �100 are shown in the boxes, with empty boxes being 100 (or �D�� � 1). Not shown are CasqSNP1, which was in complete LD with CasqSNP2, and CasqSNP3 and rs2275705, which were both rare in the Amish (allele frequency 0.034 and 0.002, respectively). SNPs showing significant association to type 2 diabetes are underlined, while SNPs showing borderline significant association to type 2 diabetes are shown in italics. TABLE 2 Haploscore results comparing type 2 diabetic cases with NGT control subjects using SNPs rs617698 and rs2275703 Rs617698 Rs2275703 Frequency Haploscore P* G A 0.286 2.29 0.022 A A 0.105 1.79 0.073 G C 0.0466 0.357 0.72 A C 0.563 2.95 0.0032 *Global P value � 0.016. CASQ1 AND TYPE 2 DIABETES IN THE AMISH 3296 DIABETES, VOL. 53, DECEMBER 2004 type 2 diabetes that was any greater than that of rs2275703 alone (data not shown). DISCUSSION Chromosome 1q21-q24 has perhaps the best replicated linkage to type 2 diabetes in diverse populations, including Caucasians, Native Americans, and Chinese (3–10). Thus, there is a high likelihood that at least one type 2 diabetes susceptibility gene resides in this region. We and others have embarked upon a systematic LD mapping and posi- tional candidate gene effort toward the positional cloning of the gene(s). Our LD mapping in Amish samples, as well as finer linkage analysis in Utah Caucasians (42), suggest at least two distinct type 2 diabetes loci within the 1q21-q24 interval. One of these regions, at �157–158.5 Mb, contains a number of good positional candidate genes, including CASQ1. Calcium release from the sarcoplasmic reticulum into the cytosol has been shown to regulate GLUT4 expression (29) and glucose transport in muscle (30). Hayashi et al. (43) showed that an increase in cytoplasmic calcium is a key step in the activation of intracellular signaling cas- cades that mediate the immediate and prolonged effect of exercise on glucose transport. CASQ is present in muscle cells and located in the lumen of the sarcoplasmic reticu- lum. The major CASQ isoforms in skeletal (CASQ1) and cardiac (CASQ2) muscle exist as predominantly fast and slow isoforms, respectively (44). CASQ1 binds and re- leases large quantities of Ca2� rapidly through its high capacity (40 –50 Ca2� per CASQ molecule) and relatively low-affinity interactions with Ca2� (28,45). Howarth et al. (31) investigated the expression of CASQ and Ca2� bind- ing in cardiac and skeletal muscle from the streptozotocin- induced diabetic rat and found no significant changes in heart but an increase in the relative abundance of CASQ and CASQ-like proteins in skeletal muscle. Together, these data suggest that CASQ1 may play an important role in insulin- and/or contraction-stimulated glucose transport. We thus hypothesized that sequence variants in CASQ1 that alter expression or function may change intracellular Ca2� flux and decrease glucose transport leading to glu- cose intolerance and/or diabetes. We have extensively screened CASQ1 for variants in Amish subjects with type 2 diabetes. We identified four novel SNPs (three in the promoter and one in exon 3), none of which encode variants that would be expected to alter the protein’s structure. Although it is possible that the SNPs in the promoter may alter CASQ1 expression, none were associated with type 2 diabetes. Genotyping and association analysis of 10 informative polymorphisms re- vealed that two were significantly associated with type 2 diabetes and three others showed borderline evidence for association with type 2 diabetes. The two most signifi- cantly associated SNPs, rs2275703 and rs617698, are lo- cated in introns 4 and 2, respectively. Interestingly, the more common alleles of rs2275703 and rs617698 were the at-risk alleles, which is in line with the common disease, common variant hypothesis and with findings for other type 2 diabetes susceptibility genes, e.g., calpain 10 and Pro12Ala persoxisome proliferator–activated receptor 2. Supportive evidence that one or more of these variants influence glucose homeostasis is that allele frequencies of type 2 diabetes–associated SNPs were intermediate be- tween type 2 diabetes and NGT in subjects with IGT (Table 1). Furthermore, rs2275703 and rs617698 were associated with glucose levels in nondiabetic subjects (Fig. 2). Unfortunately, since our original linkage was obtained from large multiplex Amish families, determining how much of the original linkage can be accounted for by the CASQ1 SNPs is not straightforward. However, given the relatively modest evidence for association of CASQ1 SNPs with type 2 diabetes, it is unlikely that these SNPs (or haplotypes) can explain our original linkage result. Indeed, two-point linkage analysis with CASQ1 SNPs showed little evidence for linkage to type 2 diabetes or type 2 diabetes/ IGT. Might these associations represent false-positive results due to multiple comparisons? For association studies such as these, this is always a possibility. Indeed a full Bonfer- roni correction for multiple comparisons would render our associations nonsignificant (at the P � 0.05 level). How- ever, we believe that correction for multiple comparisons is overly conservative, since some of our traits are corre- lated with each other. Furthermore, we have evidence for linkage in this region not only in the Amish but also in several other populations, which led to this hypothesis- driven LD mapping and positional candidate gene effort. Further reducing concern regarding false-positive results are two internal consistencies. First, the same SNPs showed association to type 2 diabetes and glucose traits in nondiabetic subjects. Second, for 8 of the 10 SNPs, the allele frequencies for the IGT cases were intermediate between those for the type 2 diabetes cases and the NGT control subjects. Perhaps most compelling is that our findings are similar to studies conducted in an independent study of Utah Caucasians showing association of CASQ1 polymorphism with type 2 diabetes (46). The region of association with type 2 diabetes that overlaps in the Amish and the more outbred Utah Caucasian population appears to be within intron 2, specifically rs617599 and rs617598 in Utah Cau- casians and rs617598 in the Amish (although this SNP was not the most significantly associated SNP in the Amish). Rs617599 and rs617598 are only 87 bp from one another and in moderate LD (d � 0.82, r2 � 0.62) in the Amish. However, contrary to our findings in the Amish, in which the rs617598 A allele was the at-risk allele for type 2 diabetes, the frequency of the G allele was significantly higher in Utah Caucasian subjects with type 2 diabetes. This discrepancy between the two populations may indi- cate that this SNP is not the functional SNP but is marking an at-risk haplotype that differs between the Amish and Utah Caucasians. Alternatively, these discrepancies be- tween populations could represent false-positive or false- negative results. The findings reported here do not provide insights into possible functional consequences of the variants in CASQ1. It is possible that variants in introns may affect mRNA splicing and/or expression. However, we cannot rule out the possibility that other variants in CASQ1 or other nearby genes in LD with CASQ1 SNPs might be the responsible functional variants. We have observed evidence for association of other 1q21-q24 positional candidate genes with type 2 diabetes M. FU AND ASSOCIATES DIABETES, VOL. 53, DECEMBER 2004 3297 or related traits, namely lamin A (LMNA) (26), guanine exchange factor-11 (ARHGEF11) (M.S., A.R.S., unpub- lished observations), and omentin 2 (ITLN2) (27). While all of these genes reside under our linkage peak, they are 4.0, 3.2, and 0.75 Mb from CASQ1, respectively. There does not appear to be any significant LD for the two CASQ1 type 2 diabetes–associated SNPs, rs617598 and rs2275703, with SNPs in LMNA, ARHGEF11, or ITLN2, suggesting that the associations of variation in these genes with type 2 diabetes and related traits are independent of one another. These findings, coupled with more than one association signal in Pima Indians and Utah Caucasians, suggest that there are likely to be more than one gene (possibly several genes) responsible for the linkage at 1q21-q24. In summary, we have reported a systematic search for polymorphism in CASQ1 on chromosome 1q21 and iden- tified four novel SNPs. None predict alterations in the protein sequence. We describe for the first time a signifi- cant association between CASQ1 polymorphism and type 2 diabetes, possibly suggesting a new susceptibility gene for type 2 diabetes. Further analysis in other populations, as well as functional studies, will be necessary to further elucidate the role of polymorphism in CASQ1 in the pathogenesis of type 2 diabetes. ACKNOWLEDGMENTS This work was supported by research grants R01- DK54261, K24-DK02673, U01 DK58026, and K07-CA67960; the University of Maryland General Clinical Research Center Grant M01 RR 16500; the General Clinical Research Centers Program, the National Center for Research Re- sources; the National Institutes of Health; and the Balti- more Veterans Administration Geriatric Research and Education Clinical Center. We gratefully acknowledge our Amish liaisons and field workers and the extraordinary cooperation and support of the Amish community, without whom these studies would not have been possible. REFERENCES 1. Horenstein RB, Shuldiner AR: Genetics of diabetes. Rev Endocr Metab Disord 5:25–36, 2004 2. McCarthy MI: Growing evidence for diabetes susceptibility genes from genome scan data. Curr Diab Rep 3:159 –167, 2003 3. Hsueh W, Jean PLS, Mitchel BD, Pollin PI, Knowler WC, Ehm MG, Bell CJ, Sakul H, Wagner MJ, Burns DK, Shuldiner AR: Genome-wide and fine- mapping linkage studies of type 2 diabetes and glucose traits in the Old Order Amish. Diabetes 52:550 –557, 2003 4. Hanson RL, Ehm MG, Pettitt DJ, Prochazka M, Thompson BD, Timberlake D, Foroud T, Kobes S, Baier L, Burns DK, Almasy L, Blangero J, Garvey WT, Bennett PH, Knowler WC: An autosomal genomic scan for loci linked to type II diabetes mellitus and body-mass index in Pima Indians. Am J Hum Genet 63:1130 –1138, 1998 5. Elbein S, Hoffman M, Teng K, Leppert M, Hasstedt S: A genome-wide search for type 2 diabetes susceptibility genes in Utah Caucasians. Diabetes 48:1175–1182, 1999 6. Vionnet N, Hani EH, Dupont S, Gallina S, Francke S, Dotte S, Matos FD, Durand E, Leprêtre F, Lecoeur C, Gallina P, Zekiri L, Dina C, Froguel P: Genomewide search for type 2 diabetes-susceptibility genes in French whites: evidence for a novel susceptibility locus for early-onset diabetes on chromosome 3q27-qter and independent replicate of a type 2-diabetes locus on chromosome1q21– q24. Am J Hum Genet 67:1470 –1480, 2000 7. Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes U.K. Warren 2 repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet 69:553–569, 2001 8. Xiang K, Wang Y, Zheng T, Jia W, Li J, Chen L, Shen K, Wu S, Lin X, Zhang G, Wang C, Wang S, Lu H, Fang Q, Shi Y, Zhang R, Xu J, Weng Q: Genome-wide search for type 2 diabetes/impaired glucose homeostasis susceptibility genes in the Chinese: significant linkage to chromosome 6q21– q23 and chromosome 1q21– q24. Diabetes 53:228 –234, 2004 9. Meigs JB, Panhuysen CI, Myers RH, Wilson PW, Cupples LA: A genome- wide scan for loci linked plasma levels of glucose an HbA1c in a community-based of Caucasian pedigrees: the Framingham Offspring Study. Diabetes 51:833– 840, 2002 10. Langefeld CD, Wagenknecht LE, Rotter JI, Williams AH, Hokanson JE, Saad MF, Bowden DW, Haffner S, Norris JM, Rich SS, Mitchell BD: Linkage of the metabolic syndrome to 1q23– q31 in Hispanic families: the Insulin Resistance Atherosclerosis Study Family Study. Diabetes 53:1170 –1174, 2004 11. Hanis Cl, Boerwinkle E, Chakraborty R, Ellsworth DL, Concannon P, Stirling B, Morrison VA, Wapelhorst B, Spielman RS, Gogolin-Ewens KJ, Shepard JM, Williams SR, Risch N, Hinds D, Iwasaki N, Ogata M, Omori Y, Petzold C, Rietzch H, Schroder HE, Schulze J, Cox NJ, Menzel S, Boriraj VV, Chen X: A genome-wide search for human non-insulin-dependent (type 2) diabetes genes reveals a major susceptibility locus on chromosome 2. Nat Genet 13:161–166, 1996 12. Mahtani MM, Widen E, Lehto M, Thomas J, McCarthy M, Brayer J, Bryant B, Chan G, Daly M, Forsblom C, Kanninen T, Kirby A, Kruglyak L, Munnelly K, Parkkonen M, Reeve-Daly MP, Weaver A, Brettin T, Duyk G, Lander ES, Groop LC: Mapping of a gene for type 2 diabetes associated with an insulin secretion defect by a genome scan in Finnish families. Nat Genet 14:90 –94, 1996 13. Ehm MG, Karnoub MC, Sakul H, Gottschalk K, Holt DC, Weber JL, Vaske D, Briley D, Briley L, Kopf J, McMillen P, Nguyen Q, Reisman M, Lai EH, Joslyn G, Shepherd NS, Bell C, Wagner MJ, Burns DK, the American Diabetes Association GENNID Study Group: Genomewide scan for type 2 diabetes susceptibility genes in four American populations. Am J Hum Genet 66:1871–1881, 2000 14. Lindgren CM, Mahtani MM, Widén E, McCarthy MI, Daly MJ, Kirby A, Reeve MP, Kruglyak L, Parker A, Meyer J, Almgren P, Lehto M, Kanninen T, Tuomi T, Groop LC, Lander ES: Genomewide search for type 2 diabetes mellitus susceptibility loci in Finnish families: the Botnia study. Am J Hum Genet 70:509 –516, 2002 15. Ghosh S, Watanabe RM, Valle TT, Hauser ER, Magnuson VL, Langefeld CD, Ally DS, Mohlke KL, Silander K, Kohtamäki K, Chines P, Balow J, Birznieks JG, 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 dia- betes mellitus genetics (FUSION) study. I. An autosomal genome scan for genes that predispose to type 2 diabetes. Am J Hum Genet 67:1174 –1185, 2000 16. Bowden DW, Sale M, Howard TD, Qadri A, Spray BJ, Rothschild CB, Akots G, Rich SS, Freedman BI: Linkage of genetic markers on human chromo- somes 20 and 12 to NIDDM in Caucasian sib pairs with a history of diabetic nephropathy. Diabetes 46:882– 886, 1997 17. Ji L, Malecki M, Warram JH, Yang Y, Rich SS, Krolewski AS: New susceptibility locus for NIDDM is localized to human chromosome 20q. Diabetes 46:876 – 881, 1997 18. Permutt MA, Wasson JC, Suarez BK, Lin J, Thomas J, Meyer J, Lewitzky S, Rennich JS, Parker A, DuPrat L, Maruti S, Chayen S, Glaser B: A genome scan for type 2 diabetes susceptibility loci in a genetically isolated population. Diabetes 50:681– 685, 2001 19. Horikawa Y, Oda N, Cox NJ, Li XQ, Orho-Melander M, Hara M, Hinokio Y, Lindner TH, Mashima H, Schwarz PEH, Bosque-plata L, Horikawa Y, Oda Y: Yoshiuchi I, Colilla S, Polonsky KS, Wei S, Concannon P, Iwasaki N, Schulze J, Baier LJ, Bogardus C, Groop L, Boerwinkle E, Hanis CL, Bell GI: Genetic variation in the gene encoding calpain-10 is associated with type 2 diabetes mellitus. Nat Genet 26:163–175, 2000 20. Love-Gregory LD, Wasson J, Ma J, Jin CH, Glaser B, Suarez BK, Permutt MA: A common polymorphism in the upstream promoter region of the hepatocyte nuclear factor-4� gene on chromosome 20q is associated with type 2 diabetes and appears to contribute to the evidence for linkage in an Ashkenazi Jewish population. Diabetes 53:1134 –1140, 2004 21. Silander K, Mohlke KL, Scott LJ, Peck EC, Hollstein P, Skol AD, Jackson AU, Deloukas P, Hunt S, Stavrides G, Chines PS, Erdos MR, Narisu N, Conneely KN, Li C, Fingerlin TE, Dhanjal SK, Valle TT, Bergman RN, Tuomilehto J, Watanabe RM, Boehnke M, Collins FS: Genetic variation CASQ1 AND TYPE 2 DIABETES IN THE AMISH 3298 DIABETES, VOL. 53, DECEMBER 2004 near the hepatocyte nuclear factor-4� gene predicts susceptibility to type 2 diabetes. Diabetes 53:1141–1149, 2004 22. Wolford JK, Hanson RL, Kobes S, Bogardus C, Prochazka M: Analysis of linkage disequilibrium between polymorphisms in the KCNJ9 gene with type 2 diabetes mellitus in Pima Indians. Mol Genet Metab 73:97–103, 2001 23. Farook VS, Hanson RL, Wolford JK, Bogardus C, Prochazka M: Molecular analysis of KCNJ10 on 1q as a candidate gene for type 2 diabetes in Pima Indians. Diabetes 51:3342–3346, 2002 24. Wang H, Chu W, Das SK, Ren Q, Hasstedt SJ, Elbein SC: Liver pyruvate kinase polymorphisms are associated with type 2 diabetes in northern European Caucasians. Diabetes 51:2861–2865, 2002 25. Wolford JK, Gruber JD, Ossowski VM, Vozarova B, Antonio Tataranni P, Bogardus C, Hanson RL: A C-reactive protein promoter polymorphism is associated with type 2 diabetes mellitus in Pima Indians. Mol Genet Metab 78:136 –144, 2003 26. Steinle NI, Kazlauskaite R, Imumorin IG, Hsueh WC, Pollin TI, O’Connell JR, Mitchell BD, Shuldiner AR: Variation in the lamin A/C gene: associa- tions with metabolic syndrome. Arterioscler Thromb Vasc 24:1708 –1713, 2004 27. Fu M, Gong D-W, Damcott C, Sabra M, Yang R, Pollin TI, Tanner K, Ott S, McLenithan JC, Fried S, O’Connell JR, Mitchell BD, Shuldiner AR: Sys- temic analysis of omentin 1 and omentin 2 on 1q23 as candidate genes for type 2 diabetes in the Old Order Amish (Abstract). Diabetes 53 (Suppl. 2):A59, 2004 28. Shin DW, Pan Z, Kim EK, Lee JM, Bhat MB, Parness J, Kim DH, Ma J: A retrograde signal from calsequestrin for the regulation of store-operated Ca2� entry in skeletal muscle. J Biol Chem 278:3286 –3292, 2003 29. Ojuka EO, Jones TE, Notlte LA, Chen M, Wamhoff BR, Sturek M, Holloszy JO: Regulation of GLUT4 biogenesis in muscle: evidence for involvement of AMPK and Ca. Am J Physiol Endocrinol Metab 283:E1008 –E1013, 2002 30. Young JH, Gulve EA, Holloszy JO: Calcium stimulates glucose transport in skeletal muscle by a pathway independent of contraction. Am J Physiol Cell Physiol 260:C555–C561, 1991 31. Howarth FC, Glover L, Culligan K, Qureshi MA, Ohlendieck K: Caseques- trin expression and calcium binding is increased in streptozotocin-induced diabetic rat skeletal muscle though not in cardiac muscle. Pflugers Arch 444:52–58, 2002 32. Agarwala R, Biesecker LG, Hopkins KA, Francomano CA, Schaffer AA: Software for constructing and verifying pedigrees within large genealogies and an application to the Old Older Amish of Lancaster County. Genome Res 8:211–221, 1998 33. Pollin TI, Agarwala R, Schaffer AA, Shuldiner AR, Mitchell BD, O’Connell JR: Investigations of male founder structure and the Y chromosome in the Old Order Amish (Abstract). Am J Hum Genet 7 (Suppl.):187, 2003 34. Hsueh WC, Mitchell BD, Aburomia R, Pollin T, Sakul H, GelderEhm M, Michelsen BK, Wagner MJ, St Jean PL, Knowler WC, Burns DK, Bell CJ, Shuldiner AR: Diabetes in the Old Older Amish: characterization and heritability analysis of the Amish Family Diabetes Study. Diabetes Care 23:595– 601, 2000 35. Hsu TM, Chen X, Duan S, Miller RD, Kwok P-Y: Universal SNP genotyping assay with fluorescence polarization detection. Biotechniques 31:560 –570, 2001 36. Bell PA, Chaturvedi S, Gelfand CA, Huang CY, Kochersperger M, Kopla R, Modica F, Pohl M, Varde S, Zhao R, Zhao X, Boyce-Jacino MT: SNPstream UHT: ultra-high throughput SNP genotyping for pharmacogenomics and drug discovery. Biotechniques (Suppl.):70 –72, 74, 76 –77, 2002 [erratum in Biotechniques 34:496, 2003] 37. O’Connell JR, Weeks DE: PedCheck: A program for identifying genotype incompatibilities in linkage analysis. Am J Hum Genet 63:259 –266, 1998 38. Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet 62:1198 –1211, 1998 39. O’Connell JR: Zero-recombinant haplotyping: applications to fine mapping using SNPs. Genet Epidemiol 19 (Suppl. 1):S64 –S70, 2000 40. Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES: High-resolution haplotype structure in the human genome. Nat Genet 29:229 –232, 2001 41. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA: Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet 70:425– 434, 2002 42. Das SK, Hasstedt SJ, Zhang Z, Elbein SC: Linkage and association mapping of a chromosome 1q21– q24 type 2 diabetes susceptibility locus in northern European Caucasians. Diabetes 53:492– 499, 2004 43. Hayashi T, Wojtaszewski JFP, Goodyear LJ: Exercise regulation of glucose transport in skeletal muscle. Am J Physiol 273:E1039 –E1051, 1997 44. Fliegel L, Ohnishi M, Carpenter MR, Khanna VK, Reithmeier RA, MacLen- nan DH: Amino acid sequence of rabbit fast-twitch skeletal muscle calsequestrin deduced from cDNA and peptide sequencing. Proc Natl Acad Sci U S A 84:1167–1171, 1987 45. Fliegel L, Leberer E, Green NM, Maclennan DH: The fast-twitch muscle calsequestrin isoform predominates in rabbit slow-twitch soleus muscle. FEBS Lett 242:297–300, 1989 46. Das SK, Chu W, Zhang Z, Hasstedt SJ, Elbein SC: Calsquestrin 1 (CASQ1) gene polymorphisms under chromosome 1q21 linkage peak are associated with type 2 diabetes in Northern European Caucasians. Diabetes 53:3300 – 3306, 2004 M. FU AND ASSOCIATES DIABETES, VOL. 53, DECEMBER 2004 3299