Calcium levels in blood serum play an important role in many biological processes. The regulation of serum calcium is under strong genetic control. This study describes the first meta-analysis of a genome-wide association study from four cohorts totaling 12,865 participants of European and Indian Asian descent. Confirming previous results in some candidate gene studies, we find that common polymorphisms at the calcium-sensing receptor (CASR) gene locus are associated with serum calcium concentrations. We show that CASR variants give rise to the strongest signals associated with serum calcium levels in both European and Indian Asian populations, while no other locus reaches genome-wide significance. Our results show that CASR is a key player in genetic regulation of serum calcium in the adult general population.
Calcium has a pivotal role in biological functions, and serum calcium levels have been associated with numerous disorders of bone and mineral metabolism, as well as with cardiovascular mortality. Here we report results from a genome-wide association study of serum calcium, integrating data from four independent cohorts including a total of 12,865 individuals of European and Indian Asian descent. Our meta-analysis shows that serum calcium is associated with SNPs in or near the calcium-sensing receptor (CASR) gene on 3q13. The top hit with a p-value of 6.3×10-37 is rs1801725, a missense variant, explaining 1.26% of the variance in serum calcium. This SNP had the strongest association in individuals of European descent, while for individuals of Indian Asian descent the top hit was rs17251221 (p = 1.1×10-21), a SNP in strong linkage disequilibrium with rs1801725. The strongest locus in CASR was shown to replicate in an independent Icelandic cohort of 4,126 individuals (p = 1.02×10-4). This genome-wide meta-analysis shows that common CASR variants modulate serum calcium levels in the adult general population, which confirms previous results in some candidate gene studies of the CASR locus. This study highlights the key role of CASR in calcium regulation.
IntroductionCalcium is the most abundant mineral in the human body contributing approximately one kilogram to the average adult human body mass. Whereas 99% of calcium is stored in the skeleton and teeth, the remaining 1% circulates in the bloodstream and is involved in many physiological processes including its function as a universal cellular signaling molecule [1]–[2]. Calcium plays a key role in membrane potential, which is important for muscle contraction, heart rate regulation and generation of nerve impulses. Calcium also influences bone metabolism, ion transport and many other cellular processes [3]. Approximately 2/5 of calcium in the extracellular fluid is found in blood serum. The level of serum calcium is under tight hormonal control with a normal concentration of 2.15–2.55 mmol/L. Serum calcium is under strong genetic control, with twin studies showing that the variance in total calcium due to genetic effects is between 50% and 78% [4]–[5].While skeletal calcium is important in numerous clinical disorders, in particular bone and mineral disorders, the clinical role of serum calcium is less clear. Several [6]–[7] (but not all [8]) studies indicated that elevated serum calcium levels are associated with an increased risk of cardiovascular disease. Patients with hyperparathyroidism, who suffer from chronic hypercalcemia, have a high prevalence of hypertension and increased cardiovascular mortality [9]–[11]. However, the mechanisms underlying the putative association of serum calcium with increased cardiovascular morbidity and mortality remain unclear.Rare monogenic forms of hypo- or hypercalcemia have been described, including disorders involving the calcium-sensing receptor (CASR, locus 3q13) gene. Heterozygous and homozygous CASR mutations that inactivate CASR are responsible, respectively, for familial hypocalciuric hypercalcemia, type 1 (also known as familial benign hypercalcemia) (OMIM #145980) [12]–[13] and neonatal severe hyperparathyroidism (OMIM #239200) [13]. On the other hand, mutations that result in CASR activation lead to autosomal dominant hypocalcemia (OMIM #146200) [14]. Mutations in many other genes have also been found to lead to disturbances of serum calcium levels (Table 1).10.1371/journal.pgen.1001035.t001Table 1Serum calcium candidate genes.GeneGene regionDiseaseOMIM# SNPs*
Top SNP†
Top SNP p-value† AIRE
chr21: 44,530,191-44,542,530Autoimmune polyendocrine syndrome, type I240300870rs2838473, rs13052277, rs7171772.48E-04, 2.78E-05, 4.55E-03
ALPL
chr1: 21,581,175-21,650,208Hypophosphatasia, infantile241500786rs6426723, rs1256348, rs46549736.88E-03, 1.14E-03, 7.30E-04
BSND
chr1: 55,176,638-55,186,485Bartter syndrome, type 46025221007rs17111592, rs11584093, rs65885284.61E-03, 8.40E-05, 1.70E-02
CASR
chr3:123,385,220-123,488,032Familial hypocalciuric hypercalcemia, type 1; neonatal severe hyperparathyroidism; autosomal dominant hypocalcemia145980; 239200; 146200; 241200921rs1801725, rs1801725, rs172512216.29E-37, 2.58E-18, 1.07E-21
CDKN1B
chr12: 12,761,576-12,766,569Bartter syndrome, type 4131100855rs3825271, rs888200, rs110552257.05E-04, 3.56E-03, 8.66E-04
CLCNKA
chr1: 16,093,672-16,105,850Bartter syndrome, type 4602522545rs12405694, rs16852052, rs66610121.34E-02, 2.51E-02, 1.00E-02
CLCNKB
chr1: 16,115,658-16,128,782Bartter syndrome, type 3; Bartter syndrome, type 4607364; 602522520rs12405694, rs16852052, rs66610121.34E-02, 2.51E-02, 1.00E-02
CLDN16
chr3: 191,588,543-191,611,035Hypomagnesemia 3, renal2482501381rs11714779, rs11714779, rs96825993.15E-04, 2.86E-04, 3.18E-02
CYP27B1
chr12: 56,442,384-56,447,145Vitamin D-dependant rickets type I264700567rs11172284, rs810204, rs7159302.18E-03, 2.36E-03, 3.26E-02
GATA3
chr10: 8,136,673-8,157,170Hypoparathyroidism; sensorineural deafness; renal disease1462551492rs11812109, rs12359361, rs27653991.30E-03, 2.86E-03, 4.35E-03
GCM2
chr6: 10,981,450-10,990,084Familial hyperparathyroidism146200908rs16870899, rs16870899, rs64571604.88E-03, 4.59E-03, 3.36E-03
GNAS
chr20: 56,900,130-56,919,640Pseudohypoparathyroidism, type IA; Pseudohypoparathyroidism, type IB103580; 603233949rs2145477, rs911297, rs60153752.38E-04, 3.95E-03, 1.12E-02
HRPT2
chr1: 189,822,81-189,952,713Hyperparathyroidism (familial isolated hyperparathyroidism); parathyroid carcinoma145000; 608266836rs10737627, rs913478, rs28876131.15E-02, 6.50E-03, 4.71E-02
KCNJ1
chr11: 128,213,12-128,242,478Bartter syndrome, antenatal, type 22412001059rs948215, rs3897566, rs71166061.08E-03, 1.50E-03, 3.46E-05
MEN1
chr11: 64,327,572-64,335,342Hyperparathyroidism (familial isolated hyperparathyroidism); multiple endocrine neoplasia, type I145000; 131100498rs7947143, rs7947143, rs118203221.49E-02, 6.83E-03, 2.92E-03
PHEX
chrX: 21,810,216-22,025,985Hypophosphatemic rickets, X-linked dominant307800NANANA
PTH
chr11: 13,470,178-13,474,143Familial hyperparathyroidism1462001292rs10832087, rs10500780, rs15022429.61E-04, 1.71E-03, 1.96E-03
PTH1R
chr3: 46,894,240-46,920,291Jansen's metaphyseal chondrodysplasia156400511rs1402151, rs883739, rs64420371.03E-02, 7.54E-03, 2.18E-02
RET
chr10: 42,892,533-42,944,955Multiple endocrine neoplasia, type I131100856rs3026762, rs3026762, rs122657922.10E-04, 1.90E-04, 4.10E-02
SLC12A1
chr15: 46,285,790-46,383,568Bartter syndrome, antenatal, type I601678925rs1025759, rs596942, rs9191292.25E-03, 1.29E-02, 2.30E-04
SLC4A1
chr17: 39,682,566-39,700,993Renal tubular acidosis, distal, autosomal dominant179800494rs12602991, rs12602991, rs7083842.18E-03, 7.20E-03, 4.05E-03
TBCE
chr1: 231,856,81-231,938,321Hypoparathyroidism-retardation-dysmorphism syndrome241410591rs12133603, rs12133603, rs2913532.34E-03, 2.18E-03, 5.38E-03
TRPM6
chr9: 74,566,965-74,732,564Hypomagnesemia with secondary hypocalcemia6020141123rs877809, rs877809, rs125509039.31E-04, 1.20E-03, 6.14E-03
VDR
chr12: 46,521,589-46,585,081Vitamin D-resistant rickets type II; vitamin D-dependent rickets, type II277440; 2597001066rs1859441, rs1859441, rs111683541.22E-03, 9.59E-04, 3.34E-04Genes which have been shown to lead to disturbances of serum calcium levels. For each gene, we report the top SNP for the meta-analyses of all cohorts, European cohorts only and Indian Asian cohorts only.*500 kb upstream and downstream the gene region.†Combined, European, and Indian Asian.In the present study, we report results obtained from meta-analysis of genome-wide associations of serum calcium levels from four cohorts with a total of 12,865 participants. We first describe the design of the study and its main finding, that variants in CASR give rise to the strongest signals associated with serum calcium levels in both European and Indian Asian populations. Our results confirm previous studies showing that mutations in CASR are associated with serum calcium levels in young healthy women [15]–[16] and extend this observation to men and women across a large spectrum of age. We show that CASR is a key player in the genetic regulation of serum calcium in men and women from the general adult population.ResultsWe performed a meta-analysis for genome-wide associations of serum calcium, determined by subtracting the estimated amount of calcium bound to albumin from the total serum calcium, to infer the amount of ionized calcium (see Materials and Methods). Our study included four cohorts: (i) 5404 European individuals from the Cohorte Lausanne (CoLaus) [17]–[18], (ii) 5548 European and Indian Asian individuals from the London Life Sciences Population (LOLIPOP) Study from West London UK [19]–[20], (iii) 1196 European individuals from the InCHIANTI Study (Tuscany, Italy) [21], and (iv) 717 individuals of European descent from the Baltimore Longitudinal Study of Aging (BLSA) study based in the Baltimore-Washington DC area [22], totaling 12,865 participants (see Table 2 for more detailed characteristics of each cohort).10.1371/journal.pgen.1001035.t002Table 2Characteristics of participants, by study.CoLausLOLIPOP European WhitesLOLIPOP Indian AsiansInCHIANTIBLSAdeCODESample size54041601394711967174126Gender (males/females)2542/28621397/2043832/115533/663390/3271313/2813Age(years)*
53.43(34.9,75.4)54.5(22.6,75.0)50.7(35.0,74.9)68.22(21,102)70.4(22,98)60.1(7,103)Pre/Post menopause1210/1652100/10456/5979/58438/289872/1836 (105 with pre and post measurements)Serum calcium(mmol/L)*
2.29(0.094)2.41(0.12)2.37(0.11)2.36(0.10)2.31(0.11)2.38(0.14)Corrected serumCalcium (mmol/L)*
2.18(0.09)2.31(0.09)2.29(0.09)2.30(0.09)2.30(0.1)2.31(0.14)Serum albumin[g/L]*
44.2(2.5)43.7(2.9)43.4(2.9)42.3(3.1)40.4(3.5)42.9(3.9)Characteristics are shown for CoLaus, LOLIPOP European, LOLIPOP Indian Asian, InCHIANTI, BLSA and deCODE. Corrected serum calcium, designed to estimate the amount of biologically active serum calcium, is defined as Ca_corrected = total serum calcium [mmol/L] + (40 - albumin [g/L])/40.*Values represent mean (range) or mean (sd).Genome-wide association scans were performed first independently for each cohort using linear regression and then the effect sizes from each cohort were meta-analyzed (see Materials and Methods). Due to the possibility of population substructure obscuring effects of genetic variants, meta-analysis was performed separately for (i) combined European and Indian Asian cohorts (N = 12,865) and restricted to cohorts of (ii) European (N = 8,919), and (iii) Indian Asian descent (N = 3,947). The meta-analyses yielded 100 SNPs from the combined cohorts, 70 SNPs when restricting to European cohorts and 22 SNPs restricting to Indian Asian cohorts that exceeded the genome-wide significance threshold of 5×10-8 (Figure 1A–1C) (the full list is provided in Table S2A, S2B, S2C). All SNPs reaching statistical significance clustered around the CASR locus at 3q13. The most significant SNP in the (i) combined and (ii) European meta-analyses was rs1801725 (p = 6.29×10-37, p = 2.58×10-18, respectively) and in the (iii) Indian Asian meta-analysis was rs17251221 (p = 1.07×10-21). These two SNPs are less than 11 kb apart and are in high linkage disequilibrium with each other (r2 = 0.946, 0.494, 1.0, 1.0 in HapMap CEU, CHB, JPT, YRI, respectively), and therefore most likely derive from the same association signal. We find that rs1801725 explains 1.26% of the variance in serum calcium, with the effect sizes and standard errors of the serum calcium increasing T allele in individual cohorts shown in Figure 2 and Table S3. According to our additive model, each rs1801725 T allele increases log10 serum calcium (in units mmol/L) by 3.61×10-3, equivalent to a multiplicative effect of 1.008 on serum calcium (see also Table S2). At an average serum calcium level of 2.25 mmol/L, each rs1801725 T allele yields an increase of 0.01874 mmol/L, or 21% of one standard deviation of serum calcium levels in a normal population. The regional pattern of association of SNPs around the CASR locus, and their linkage disequilibrium with rs1801725, are shown in Figure 3. Of note, rs1042636, which has been associated with decreased serum calcium [23], also achieved genome-wide significance with the G minor allele associated with decreased serum calcuim (p = 4.96×10-9). However, conditional on the rs1801725 locus, located 12 bps upstream, the rs1042636 p-value became 3.32×10−4, indicating that the two loci share contributions to serum calcium levels.10.1371/journal.pgen.1001035.g001Figure 1Genome-wide association results.Manhattan plots showing significance of association of all SNPs in the meta-analysis for (A) combined European and Indian Asian cohorts, (B) European cohorts only and (C) Indian Asian cohorts only. SNPs are plotted on the x-axis according to their position on each chromosome against association with serum calcium concentrations on the y-axis (shown as −log10 p-values).10.1371/journal.pgen.1001035.g002Figure 2Comparison of rs1801725 significance across cohorts.The effect size and 95% confidence intervals of the serum calcium increasing T allele of SNP rs1801725 are shown separately for each cohort (CoLaus, LOLIPOP_EWA, LOLIPOP_EWP, LOLIPOP_IAA, LOLIPOP_IAI, LOLIPOP_IAP, BLSA, InCHIANTI) and for the replication cohort deCODE. European cohorts are shown in blue and Indian Asian cohorts are drawn in green. The size of the box is proportional to the precision 1/se2 and the meta-analysis estimate and 95% confidence interval across all cohorts is given by a diamond.10.1371/journal.pgen.1001035.g003Figure 3Regional association plot of the CASR locus.Plots show genomic position on the x-axis and −log10 p-values on the y-axis for SNPs in the CASR locus on chromosome 3. The sentinel hit is shown as a red square. Patterns of linkage disequilibrium between the sentinel SNP and all other SNPs are color-coded. Red circles indicate high correlation (r2>0.8), green circles indicate moderate correlation (r2>0.5), blue circles indicate low correlation (r2>0.2) and black circles indicate no correlation (r2<0.2). The fine-scale recombination rates from [72]–[73] are plotted in light blue.To confirm the rs1801725 signal, we analyzed the association pattern with serum calcium in a separate cohort. We used a subset of 4,126 Icelandic individuals from the deCODE study [24]–[25] with serum calcium measurements. We found the rs1801725 T allele to be strongly associated with increased serum calcium (p = 1.02×10-4), replicating the key meta-analysis result.While only the CASR locus reached nominal genome-wide significance for association with serum calcium, the top regions with p<10-5 are shown in Table 3. These SNPs cover 12 regions, the significance of which is displayed across cohorts in Figure S3. There were no SNPs in other candidate genes (which have previously been shown to be involved in disorders associated with disturbed serum calcium levels) that were associated with serum calcium at genome-wide significance. The most significant SNPs within 500 kb of the gene transcripts are shown in Table 1. Considering the set of 18,611 distinct SNPs mapping to the set of serum calcium candidate genes excluding CASR, we find no significant association (at significance level 0.05 and applying the Bonferroni correction for multiple testing, giving a cut-off p-value of 2.69×10-6, see also Figure S4). Indeed, fixing the sample size and genome-wide significance threshold our study is well-powered (≥0.80) to detect SNPs explaining at least 0.31% of the variance. Therefore the common SNPs within the candidate genes (excluding CASR) likely play at best a small role in serum calcium regulation.10.1371/journal.pgen.1001035.t003Table 3Significance of top SNPs.Combined European and Indian Asian Cohortsdb SNPChrPosition (Build 35)Coded AlleleNon- Coded AlleleFrequency Coded AlleleBetaSeLower 95% CIUpper 95% CIP-valueR2
rs18017253123486447TG17.753.61E-032.81E-043.06E-034.16E-03
6.29E-37
0.0126rs17120351814731368TC96.244.03E-037.67E-042.52E-035.53E-032.06E-070.0021rs74480175117800594TG17.83-1.73E-033.57E-04-2.43E-03-1.03E-031.65E-060.0018rs742393125271187GA78.862.11E-034.44E-041.24E-032.98E-032.41E-060.0018rs16827695140956147GA94.523.46E-037.41E-042.01E-034.92E-033.74E-060.0017rs15505322234046848GC74.24-1.15E-032.48E-04-1.64E-03-6.66E-044.28E-060.0017rs2615031381305121GA91.352.45E-035.29E-041.41E-033.49E-034.60E-060.0017rs101191950098513GA73.82-1.63E-033.53E-04-2.33E-03-9.41E-044.76E-060.0017rs176664606149298132GA11.161.93E-034.18E-041.11E-032.75E-034.81E-060.0017rs169024868129024247GC4.26-2.66E-035.75E-04-3.78E-03-1.53E-034.83E-060.0017rs17005914270721826TC21.4-1.47E-033.24E-04-2.11E-03-8.38E-046.76E-060.0016rs10455097674550153CA47.75-9.58E-042.15E-04-1.38E-03-5.37E-049.94E-060.0015European Cohortsdb SNPChrPosition (Build 35)Coded AlleleNon-Coded AlleleFrequency Coded AlleleBetaSeLower 95% CIUpper 95% CIP-valueR2
rs18017253123486447TG16.763.11E-033.55E-042.41E-033.80E-03
2.58E-18
0.0085rs17120351814731368TC96.244.03E-037.67E-042.52E-035.53E-031.69E-070.0031rs74480175117800594TG17.83-1.73E-033.57E-04-2.43E-03-1.03E-031.40E-060.0026rs16827695140956147GA94.453.46E-037.41E-042.01E-034.92E-033.20E-060.0024rs2615031381305121GA91.352.45E-035.29E-041.41E-033.49E-033.95E-060.0024rs169024868129024247GC4.7-2.92E-036.30E-04-4.15E-03-1.68E-034.04E-060.0024rs176664606149298132GA11.161.93E-034.18E-041.11E-032.75E-034.13E-060.0024rs123251141672861392TC76.14-1.40E-033.05E-04-2.00E-03-8.01E-044.86E-060.0024rs64273101150628843TC16.29-1.65E-033.60E-04-2.35E-03-9.41E-045.09E-060.0023rs17005914270721826TC21.4-1.47E-033.24E-04-2.11E-03-8.38E-045.83E-060.0023rs124166681083645028TC94.94-2.77E-036.18E-04-3.98E-03-1.56E-037.60E-060.0023rs61110212016221853CA41.751.20E-032.67E-046.71E-041.72E-038.36E-060.002Indian Asian Cohortsdb SNPChrPosition (Build 35)Coded AlleleNon-Coded AlleleFrequency Coded AlleleBetaSeLower 95% CIUpper 95% CIP-valueR2
rs172512213123475937GA19.244.72E-034.86E-043.76E-035.67E-03
1.07E-21
0.0233rs132033356166194438TC28.04-2.91E-035.70E-04-4.03E-03-1.79E-034.70E-070.0066rs4695355448030638GA40.261.72E-033.75E-049.89E-042.46E-035.70E-060.0053rs1084691712124568003TC59.691.74E-033.78E-049.98E-042.48E-035.79E-060.0053We report SNPs with p-values <1E-05 filtered by distinct regions, determined by merging SNPs within 1 Mb of each other. Results are given for different cohort subsets of combined European and Indian Asian cohorts, European cohorts, and Indian Asian cohorts. The Beta and SE values represent the effect of the coded allele on log10 corrected serum calcium levels; R2 represents the fraction of variation explained by the SNP.We analyzed the association of the top SNP with several calcium-related outcomes (coronary heart disease, myocardial infarction, hypertension, stroke, osteoarthritis, osteoporosis and kidney stones). The number of cases and controls for each outcome and each cohort is given in Table S4. Logistic regression including age and pseudosex (see Materials and Methods) as covariates did not find any significant association between rs1801725 and the calcium-related outcomes, after correcting for multiple testing (effect sizes and standard errors for the T allele are listed in Table S5). Power calculations show that given the sample sizes for the clinical traits above, our study has good power (≥0.80) to detect odds ratios of 1.20, 1.13, 1.77, 1.27, 1.27, 1.24 and 1.75, respectively. As the smallest p-values from calcium-related traits were for osteoarthritis and osteoporosis (bonferroni-corrected p = 0.21, 0.44, respectively), we further investigated bone density traits. None of deCODE hip bone mineral density or spine bone mineral density (N = 6657 and 6838, respectively) nor InCHIANTI total bone density, trabecular bone density, cortical bone density, cortical bone thickness or cortical bone area (N = 1196) bonferroni-adjusted p-values for eight traits were significant.DiscussionThis genome-wide scan of 12,865 individuals revealed CASR as the most significant (and only genome-wide significant) locus influencing serum calcium levels. Specifically, we found evidence for a strong association of SNPs located in the CASR locus with serum calcium levels in both Europeans and Indian Asians. The strongest locus in CASR was further shown to replicate in an independent Icelandic cohort of 4,126 individuals.The top signal (rs1801725, 2956G>T) explains 1.26% of the variance in serum calcium. Indeed, this is similar to results from other GWAS of human height [26]–[29], body mass index [30]–[31], serum urate [32]–[34] and serum lipid concentrations [34]–[36], for which the genome-wide significant loci uncovered thus far collectively explain only a small fraction of the phenotypic variance (usually at least one order of magnitude less than the total additive genetic variance estimated from heritability studies [37]–[38]). The rs1801725 T allele (A986S) was associated with higher serum calcium, consistent with previous findings (see Table S6). The rs1801725 polymorphism (with T allele frequencies of 16.76%, 19.98% in European and Indian Asian cohorts, respectively) affects serum calcium levels of a substantial proportion of the population.The rs1801725 polymorphism encodes a missense variant in exon 7 of the CASR gene, which leads to a non-conservative amino-acid change (serine substitution for alanine-986, A986S corresponding to nucleotides 2956G>T) in the cytoplasmic tail of CASR. In vitro studies showed that mutations within the C-terminal tail may influence several aspects of CASR function, such as signal transduction, intracellular trafficking and cell surface expression [39]–[41]. However, PolyPhen predicts rs1801725 to be a benign substitution. It is presently unclear whether this substitution gives rise to functional variants, as functional studies have yielded conflicting results [42]–[43]. Deep sequencing of this region may help identifying the causal variants. While it is still not possible to infer a direct causal role, it is of interest to note that the SNP gives rise to an amino acid change in the C-terminal tail of CASR, a domain which plays a key role in the receptor function and may potentially influence intracellular trafficking following CASR activation by extracellular calcium.Several studies have reported associations of A986S and nearby CASR mutations with various phenotypes. The A986S CASR polymorphism has been associated with variations in circulating calcium levels in healthy adults in some studies [15], [23], [44]–[45], but not in others [46]–[47]. The fact that the latter studies were underpowered (sample size ranging from 148 to 1252) to detect a small effect size likely explains these inconsistent results. The rs1042636 (R990G) polymorphism has been associated with the magnitude of parathyroid hormone (PTH) secretion in patients with primary hyperparathyroidism [48], and preliminary results suggest that it could influence response to cinacalcet, a calcimimetic used to treat secondary hyperparathyroidism in patients with end-stage renal disease [49]. In a meta-analysis, 986S was associated with a 49% increased risk (P = 0.002) of primary hyperparathyroidism [47], [50]–[51]. Among patients with primary hyperparathyroidism, the AGQ haplotype (i.e. 986A, 990G, 1011Q, which is associated with lower serum calcium levels and hypercalciuria [52]) was associated with increased risk, and the SRQ haplotype with decreased risk, of kidney stones [50].
CASR has been previously considered as a candidate gene for osteoporosis [53] and coronary heart disease as well as increased total and cardiovascular mortality [54]. In our meta-analysis, we found no significant association of rs1801725 with these calcium-related phenotypes. A recent meta-analysis focusing on effects of candidate genes on osteoporosis also reports negative results for CASR
[55]. Furthermore, results on the association of elevated serum calcium with increased cardiovascular risk in the general population are controversial [6]–[8]. It is therefore not clear to what extent serum calcium might predict cardiovascular risk. The SNPs identified in this meta-analysis could serve as genetic instruments in future studies, such as Mendelian randomization analysis in longitudinal cohorts, to further investigate the causal effect of serum calcium on osteoporosis and on cardiovascular disease risk (see Table S5 for rs1801725 effects and standard errors).Our meta-analysis suffers from some limitations. First, we used corrected serum calcium and not directly measured ionized serum calcium. The correlation between corrected serum calcium and ionized serum calcium varies between 0.66 and 0.87 [56]–[58]. We can hypothesize that the association of ionized serum calcium with CASR variants would be stronger than the one with corrected serum calcium because ionized calcium is the form physiologically active on CASR. Second, data on serum phosphate, PTH or vitamin D are not available, so that we cannot explore further these relationships. Third, sample sizes for calcium-related clinical traits were limited, many clinical traits in CoLaus were self-reported instead of clinically diagnosed, and we incur a multiple testing penalty due to the number of clinical traits posited to be associated with serum calcium. However, the major strengths of the study are the hypothesis-free nature of GWAS studies, the large sample meta-analysis and the inclusion of multiple populations.Materials and MethodsCohortsCoLaus is a population-based sample from Lausanne, Switzerland, consisting of 5435 individuals between 35 and 75 years old (after QC) of which a subset of 5404 had available serum calcium measurements. The study design and protocols have been described previously [17]–[18]. The CoLaus study was approved by the Institutional Ethic's Committee of the University of Lausanne. The London Life Sciences Prospective Population Study (LOLIPOP) is an ongoing population-based cohort study of ∼30,000 Indian Asian and European white men and women, aged 35–75 years living in West London, United Kingdom [59]. All study participants gave written consent including for genetic studies. The LOLIPOP study is approved by the local Research Ethics Committee. The participants included in the present study are a subset of 3947 Indian Asians and 1601 Europeans from the LOLIPOP cohort study. LOLIPOP individuals are separated by origin as well as the genotyping platform, with IAA, IAI or IAP denoting Indian Asians genotyped on Affymetrix, Illumina or Perlegen platforms, respectively, and EWA and EWI denoting Europeans genotyped on Affymetrix or Illumina platforms, respectively (see Table S1). The InCHIANTI study is a population-based epidemiological study aimed at evaluating the factors that influence mobility in the older population living in the Chianti region in Tuscany, Italy. The details of the study have been previously reported [60]. Overnight fasted blood samples were taken for genomic DNA extraction, and measurement of serum calcium. For this study, 1196 subjects with serum calcium and GWAS data were analyzed. The study protocol was approved by the Italian National Institute of Research and Care of Aging Institutional Review and Medstar Research Institute (Baltimore, MD). The Baltimore longitudinal study on Aging (BLSA) study is a population-based study aimed to evaluate contributors of healthy aging in the older population residing predominantly in the Baltimore-Washington DC area [61]. Starting in 1958, participants are examined every one to four years depending on their age. Blood samples were collected for DNA extraction. This analysis focused on a subset of the participants (N = 717) of European ancestry. The BLSA has continuing approval from the Institutional Review Board (IRB) of Medstar Research Institute. Approval was obtained from local ethic committees for all studies and all participants signed informed written consent. The deCODE study consists of individuals who visited a private outpatient laboratory, the Laboratory in Mjodd, Reykjavik, Iceland between 1997 and 2008. The main referral center for this laboratory is a multispecialty medical clinic in Reykjavik (Laeknasetrid). For the serum calcium analysis we used information on 4,126 individuals with both genome-wide SNP data and measured serum calcium and serum albumin. The samples for bone density analysis have previously been described in detail [24]–[25]. For this study 6,657 individuals with total hip bone mineral density (BMD) and 6,838 individuals with lumbar spine BMD and SNP data were available for analysis. All participants gave informed consent and the study was approved by the Data Protection Commission of Iceland (DPC) and the National Bioethics Committee of Iceland.Clinical dataFor each CoLaus participant a venous blood sample was collected under fasting conditions. Measurements were conducted using a Modular P apparatus (Roche Diagnostics, Switzerland). Total serum calcium was measured by O-cresolphtalein (2.1% – 1.5% maximum inter and intra-batch CVs); albumin was measured by bromocresol green (2.5% – 0.4%). To further characterize the identified genetic variants, we analyzed the association with several outcomes postulated to be correlated with serum calcium. Within the CoLaus study, we have questionnaire responses to queries about personal histories of osteoporosis, osteoarthritis, myocardial infarction and stroke in addition to clinical data determining hypertension status, defined as previously described [17]. The assessment of LOLIPOP study participants was carried out by a trained research nurse, during a 45 minute appointment according to a standardized protocol and with regular QC audits. An interviewer-administered questionnaire was used to collect data on medical history, family history, current prescribed medication, and cardiovascular risk factors. Physical assessment included anthropometric measurements (height, weight, waist, hip) and blood pressure. Blood was collected after an 8 hour fast for biochemical analysis, including glucose, insulin, total and HDL cholesterol and triglycerides, and whole blood was taken for DNA extraction [59]. InCHIANTI serum albumin concentrations were determined as percentage of total protein using agarose electrophoretic technique (Hydragel Protein (E) 15/30, Sebia, Issy-les-Moulineaux, France). Serum calcium was measured using calorimetric assay (Roch Diagnostic, GmbH, Mannheim, Germany) by a Roche-Hitachi autoanalyzer (The intra-assay CV and 0.9% and the inter-assay CV was 1.5%). Measures of bone density, bone dimensions and osteoporosis diagnosis were assessed by peripheral quantitative computed tomography (pQCT) using the XCT 2000 device (Stratec Medizintechnik, Pforzheim, Germany) [62]. BLSA albumin concentrations were measured by a calorimetric assay using bromocresol green (Ortho-Clinical Diagnostics). Calcium concentrations were measured by a calorimetric assay (Vitros 5,1,FS).Genome-wide genotyping and imputationCoLaus participants were genotyped using Affymetrix Human Mapping 500 K Array. For the genome-wide association stage, genotyping in LOLIPOP participants was carried out using the Illumina 317 K mapping array, Affymetrix Human Mapping 500 K array, and Perlegen, 284 K platforms (Table S1). Participants of the InCHIANTI and BLSA studies were genotyped using Illumina Infinium HumanHap 550 K SNP arrays were used for genotyping [21]. Imputation of allele dosage of SNPs was performed using either MACH [63] or IMPUTE [64] with parameters and quality control filters as described in . All European cohorts imputed SNPs typed in the HapMap CEU population; LOLIPOP Indian Asian cohort imputed SNPs using mixed HapMap populations, given that this showed greater concordance with real genotypes compared with use of any one HapMap population. SNPs were excluded if cohort-specific imputation quality as assessed by r2.hat (MACH) or .info (IMPUTE) metrics were <0.30. In total, 2,557,252 genotyped or imputed SNPs had data from one or more cohorts and were analyzed. Genotypes in deCODE were measured using either humanHap300, humanHap300-duo or humanCNV370.Statistical analysisIndividual genome-wide association analysisBiologically active serum calcium is estimated by the correction, Ca_corrected = total serum calcium [mmol/L] + (40 - albumin [g/L])/40. Individuals with values <1.9 or >3 were removed as these were extreme outliers. Linear-regression analyses were carried out using an additive genetic model on log10-transformed corrected calcium levels adjusted for age and pseudosex (a factor variable with three values: males, pre-menopausal females and post-menopausal females). BLSA also included the first two and LOLIPOP included the first four ancestry principal components in the regression, respectively. Regression analyses were performed with QUICKTEST [65] (CoLaus), MACH2qtl (LOLIPOP) [63] or MERLIN [66] (InCHIANTI, BLSA).Meta-analysisThe results from all cohorts were combined into a fixed-effects meta-analysis using inverse variance weighting. Tests for heterogeneity were assessed using Cochran's Q statistic and the log of the related H statistic [67] after grouping LOLIPOP subsets into European and Indian subsets. For rs1801725 and rs1042636 the p-values were (0.07657, 0.1432) and (0.3450, 0.8876), respectively, indicating limited between-study variability. The analysis was implemented in R and run on a quad-core Linux machine. SNPs were reported provided they had effect size estimates in at least 2 of the 5 European cohorts, in at least 2 of the 3 Indian Asian cohorts, or in at least 3 of the 8 total cohorts. For the overall meta-analysis, residual inflation of the test statistic was corrected using genomic control [68]. The inflation factor was 1.0207 for the all combined cohorts, 1.0068 for European cohorts and 1.0286 for Indian Asian cohorts. Where reported, individual study p-values are corrected for inflation using genomic control methods for genotyped and imputed SNPs combined (inflation factors for individual studies were 1.0139 (CoLaus), 0.9891 (LOLIPOP EWA), 0.9994 (LOLIPOP EWP), 0.9967 (LOLIPOP IAA), 1.0131 (LOLIPOP IAI), 0.9985 (LOLIPOP IAP), 0.9842 (InCHIANTI), 1.0019 (BLSA)). The regional association plot (Figure 3) was created modifying a publically available R script [69]. The map of fine-scale recombination rates was downloaded from the HapMap website http://www.hapmap.org/downloads/recombination/ using Phase II HapMap data (release 21). Quantile-quantile plots of the association results are shown in Figure S1A, S1B, S1C, study-specific quantile-quantile plots are shown in Figure S2. Associations below p = 5×10−8 were considered genome-wide significant, which corresponds to a Bonferroni correction for the estimated one million independent common variant tests in the human genome of European individuals [70]. The analysis of osteoporosis status in CoLaus and InCHIANTI was performed using logistic regression including age and pseudosex as covariates in QUICKTEST [65]. Linkage disequilibrium was estimated from HapMap CEU (2007-01, build 35 non-redundant) genotypes. LD r2 statistics were estimated for SNPs within 500 kb using Haploview [71].Association of rs1801725 with calcium-related outcomesFor each related trait (coronary heart disease, hypertension, kidney stones, myocardial infarction, osteoarthritis, osteoporosis and stroke) we performed a fixed-effects meta-analysis of the logistic regression coefficients. We applied the bonferroni correction to adjust for multiple testing. We performed Wald-based power calculations using a type I error of 0.05/7 and meta-analysis coefficient estimates and standard errors to estimate the sample size for each trait giving power 0.80.Supporting InformationFigure S1Quantile-Quantile plots of genome-wide association results. Observed -log10 p-values on the y-axis are plotted against theoretical -log10 p-values on the x-axis resulting from each meta-analysis. Results are color-coded by chromosome. The top results largely derive from the CASR locus on chromosome 3. Results are shown separately for (A) all cohorts, (B) European cohorts and (C) Indian Asian cohorts.(6.48 MB TIF)Click here for additional data file.Figure S2Study-specific quantile-quantile plots. Shown are observed -log10 p-values plotted against expected -log10 p-values resulting from each single study after applying genomic control correction. The study-specific λ-values were λ = 1.0139 (CoLaus), λ = 0.9891 (LOLIPOP_EWA), λ = 0.9994 (LOLIPOP_EWP), λ = 0.9967 (LOLIPOP_IAA), λ = 1.0131 (LOLIPOP_IAI), λ = 0.9985 (LOLIPOP_IAP), λ = 0.9842 (InCHIANTI), λ = 1.0019 (BLSA). For the combined European and Indian Asian, European only and Indian Asian only meta-analyses the inflation factors were 1.0207, 1.0068, and 1.0286, respectively.(4.36 MB TIF)Click here for additional data file.Figure S3Comparison of significance across cohorts. The effect size and 95% confidence intervals of SNPs which do not reach genome-wide significance in the combined European and Indian Asian meta-analysis are shown separately for each cohort (CoLaus, LOLIPOP_EWA, LOLIPOP_EWP, LOLIPOP_IAA, LOLIPOP_IAI, LOLIPOP_IAP, BLSA, InCHIANTI). European cohorts are drawn in blue and Indian Asian cohorts are drawn in green. The size of the box is proportional to the precision 1/se2 and the meta-analysis estimate and 95% confidence interval across all cohorts is given by a diamond.(0.38 MB TIF)Click here for additional data file.Figure S4Candidate gene QQ-plots. For 18611 SNPs mapping to candidate genes (excluding CASR), we compare observed -log10 p-values to the mean quantiles of the uniform distribution. As a comparison, we randomly choose a set of genes from which we select the same number of SNPs. From 1,000 random draws we calculate the 95th percentile of -log10 p-values (in blue). Results are shown separately for all cohorts, European only and Indian Asian only (A–C). CoLaus permuted phenotype results comparing observed p-values to the 95th percentile of -log10 p-values from 100 permutations are shown in (D).(0.41 MB TIF)Click here for additional data file.Table S1Genotyping, imputation, and analysis procedures by study. The genotyping platforms, quality control (QC) filters applied before imputation, imputation software, number of SNPs, and genotype-phenotype association software are shown for each study.(0.05 MB DOC)Click here for additional data file.Table S2Complete list of genome-wide significant SNPs. Below is the list of all SNPs that exceeded the threshold of genome-wide significance (p<5×10e-8). Position is given for NCBI Build 35. Meta-analysis is performed by inverse variance weighted fixed effect regression. The coded allele is the allele to which the beta (effect) estimate refers. Results are shown separately for (A) European and Indian Asian cohorts, (B) European cohorts only, and (C) Indian Asian cohorts only.(0.24 MB DOC)Click here for additional data file.Table S3Significance of top SNPs by cohort. Shown are study-specific results of the SNPs with genomic control (GC) p-values <1E-05 filtered by distinct regions, determined by merging SNPs within 1 Mb of each other. Results are shown separately for (A) European and Indian Asian cohorts, (B) European cohorts, and (C) Indian Asian cohorts.(0.08 MB DOC)Click here for additional data file.Table S4Number of cases and controls for calcium-related outcomes. For several related phenotypes, we test the association of rs1801725 with these binary responses. Shown here are the number of cases and controls for each phenotype in each cohort and the total across cohorts.(0.04 MB DOC)Click here for additional data file.Table S5Logistic regression of clinical phenotypes on rs1801725. We report the effect size and standard error of the rs1801725 T allele from logistic regressions of each clinical phenotype.(0.04 MB DOC)Click here for additional data file.Table S6Studies of CASR mutations and serum calcium. A survey of previous studies which investigate the relationship between CASR mutations and levels of serum calcium.(0.05 MB DOC)Click here for additional data file.
Computational Biology/Genomics Computational Biology/Population Genetics Genetics and Genomics Public Health and Epidemiology