Association of SNPs in EGR3 and ARC with Schizophrenia Supports a Biological Pathway for Schizophrenia Risk

We have previously hypothesized a biological pathway of activity-dependent synaptic plasticity proteins that addresses the dual genetic and environmental contributions to schizophrenia. Accordingly, variations in the immediate early gene EGR3, and its target ARC, should influence schizophrenia susceptibility. We used a pooled Next-Generation Sequencing approach to identify variants across these genes in U.S. populations of European (EU) and African (AA) descent. Three EGR3 and one ARC SNP were selected and genotyped for validation, and three SNPs were tested for association in a replication cohort. In the EU group of 386 schizophrenia cases and 150 controls EGR3 SNP rs1877670 and ARC SNP rs35900184 showed significant associations (p = 0.0078 and p = 0.0275, respectively). In the AA group of 185 cases and 50 controls, only the ARC SNP revealed significant association (p = 0.0448). The ARC SNP did not show association in the Han Chinese (CH) population. However, combining the EU, AA, and CH groups revealed a highly significant association of ARC SNP rs35900184 (p = 2.353 x 10⁻⁷; OR [95% CI] = 1.54 [1.310–1.820]). These findings support previously reported associations between EGR3 and schizophrenia. Moreover, this is the first report associating an ARC SNP with schizophrenia and supports recent large-scale GWAS findings implicating the ARC complex in schizophrenia risk. These results support the need for further investigation of the proposed pathway of environmentally responsive, synaptic plasticity-related, schizophrenia genes.     

Analysis of PPARGC1B,RUNX3 and TBKBP1 Polymorphisms in Chinese Han Patients with Ankylosing Spondylitis: A Case-Control Study

Background

Susceptibility to and severity of ankylosing spondylitis (AS) are largely genetically determined. PPARGC1B, RUNX3 and TBKBP1 have recently been found to be associated with AS in patients of western European descent. Our purpose is to examine the influence of PPARGC1B, RUNX3 and TBKBP1 polymorphisms on the susceptibility to and the severity of ankylosing spondylitis in Chinese ethnic majority Han population.

Methods

Blood samples are drawn from 396 AS patients and 404 unrelated healthy controls. All the patients and the controls are Han Chinese and the patients are HLA-B27 positive. The AS patients are classified based on the severity of the disease. Twelve tag single nucleotide polymorphisms (tagSNPs) in PPARGC1B, RUNX3 and TBKBP1 are selected and genotyped. Frequencies of different genotypes and alleles are analyzed among the different severity AS patients and the controls.

Results

After Bonferroni correction, the rs7379457 SNP in PPARGC1B shows significant difference when comparing all AS patients to controls (p = 0.005). This SNP also shows significant difference when comparing normal AS patients to controls (p = 0.002). The rs1395621 SNP in RUNX3 shows significant difference when comparing severe AS patients to controls (p = 0.007). The rs9438876 SNP in RUNX3 shows significant difference when comparing normal AS patients to controls (p = 0.007). The rs8070463 SNP in TBKBP1 shows significant difference in genotype distribution when comparing severe AS patients to controls (p = 0.003).

Conclusions

The rs7379457 SNP in PPARGC1B is related to susceptibility to AS in Chinese Han population. The rs7379457 SNP in PPARGC1B, the rs1395621 and rs9438876 SNPs in RUNX3, and the rs8070463 SNP in TBKBP1 are related to the severity of AS in Chinese Han population.     

WW-domain-containing oxidoreductase is associated with low plasma HDL-C levels

Low serum HDL-cholesterol (HDL-C) is a major risk factor for coronary artery disease. We performed targeted genotyping of a 12.4 Mb linked region on 16q to test for association with low HDL-C by using a regional-tag SNP strategy. We identified one SNP, rs2548861, in the WW-domain-containing oxidoreductase (WWOX) gene with region-wide significance for low HDL-C in dyslipidemic families of Mexican and European descent and in low-HDL-C cases and controls of European descent (p = 6.9 × 10⁻⁷). We extended our investigation to the population level by using two independent unascertained population-based Finnish cohorts, the cross-sectional METSIM cohort of 4,463 males and the prospective Young Finns cohort of 2,265 subjects. The combined analysis provided p = 4 × 10⁻⁴ to 2 × 10⁻⁵. Importantly, in the prospective cohort, we observed a significant longitudinal association of rs2548861 with HDL-C levels obtained at four different time points over 21 years (p = 0.003), and the T risk allele explained 1.5% of the variance in HDL-C levels. The rs2548861 resides in a highly conserved region in intron 8 of WWOX. Results from our in vitro reporter assay and electrophoretic mobility-shift assay demonstrate that this region functions as a cis-regulatory element whose associated rs2548861 SNP has a specific allelic effect and that the region forms an allele-specific DNA-nuclear-factor complex. In conclusion, analyses of 9,798 subjects show significant association between HDL-C and a WWOX variant with an allele-specific cis-regulatory function.     

Identification of functional DNA variants in the constitutive promoter region of MDM2

Although mutations in the oncoprotein murine double minute 2 (MDM2) are rare, MDM2 gene overexpression has been observed in several human tumors. Given that even modest changes in MDM2 levels might influence the p53 tumor suppressor signaling pathway, we postulated that sequence variation in the promoter region of MDM2 could lead to disregulated expression and variation in gene dosage. Two promoters have been reported for MDM2; an internal promoter (P2), which is located near the end of intron 1 and is p53-responsive, and an upstream constitutive promoter (P1), which is p53-independent. Both promoter regions contain DNA variants that could influence the expression levels of MDM2, including the well-studied single nucleotide polymorphism (SNP) SNP309, which is located in the promoter P2; i.e., upstream of exon 2. In this report, we screened the promoter P1 for DNA variants and assessed the functional impact of the corresponding SNPs. Using the dbSNP database and genotyping validation in individuals of European descent, we identified three common SNPs (?1494?G?>?A; indel 40?bp; and ?182?C?>?G). Three major promoter haplotypes were inferred by using these three promoter SNPs together with rs2279744 (SNP309). Following subcloning into a gene reporter system, we found that two of the haplotypes significantly influenced MDM2 promoter activity in a haplotype-specific manner. Site-directed mutagenesis experiments indicated that the 40?bp insertion/deletion variation is causing the observed allelic promoter activity. This study suggests that part of the variability in the MDM2 expression levels could be explained by allelic p53-independent P1 promoter activity.     

LncRNA <Emphasis Type="Italic">ANRIL</Emphasis> Expression and <Emphasis Type="Italic">ANRIL</Emphasis> Gene Polymorphisms Contribute to the Risk of Ischemic Stroke in the Chinese Han Population

The aim of the present study was to explore the role of lncRNA ANRIL in the pathogenesis of ischemic stroke (IS) and coronary artery disease (CAD) and to determine the association between ANRIL variants and the genetic susceptibility of IS and CAD in the Chinese Han population. A genetic association study including 550 IS patients, 550 CAD patients, and 550 healthy controls was conducted. The expression levels of lncRNA ANRIL, CDKN2A, and CDKN2B were detected using qRT-PCR. Genotyping was performed by Sequenom MassARRAY on an Agena platform. Our study showed that IS patients had an increased lncRNA ANRIL expression (P?=?0.002) and a decreased CDKN2A expression (P?<?0.001) compared with normal controls. A significant difference with regard to the genotype distribution of rs2383207 was found between male IS patients and controls (P?=?0.011). The minor allele of rs2383207 significantly increased the IS risk under a recessive model (OR?=?1.52, 95% CI?=?1.05–2.21, P?=?0.027). The minor allele of rs1333049 was significantly associated with the risk of IS among the male patients under a recessive model (OR?=?1.56, 95% CI?=?1.04–2.35, P?=?0.031). However, no significant association was found between the ANRIL variants and the risk of CAD (all P?>?0.050). In addition, we found a decreased lncRNA ANRIL expression in IS patients who carried the GG genotype of rs1333049 compared with IS patients who carried the CC or CG genotype (P?=?0.041). In summary, we found that IS patients had an increased lncRNA ANRIL expression and a decreased CDKN2A expression compared with the controls, which might play an impellent role in pathological processes of IS. The ANRIL variants rs2383207 and rs1333049 were significantly associated with the risk of IS among males but not females in the Chinese Han population.     

A GWAS Study on Liver Function Test Using eMERGE Network Participants

Introduction

Liver enzyme levels and total serum bilirubin are under genetic control and in recent years genome-wide population-based association studies have identified different susceptibility loci for these traits. We conducted a genome-wide association study in European ancestry participants from the Electronic Medical Records and Genomics (eMERGE) Network dataset of patient medical records with available genotyping data in order to identify genetic contributors to variability in serum bilirubin levels and other liver function tests and to compare the effects between adult and pediatric populations.

Methods

The process of whole genome imputation of eMERGE samples with standard quality control measures have been described previously. After removing missing data and outliers based on principal components (PC) analyses, 3294 samples from European ancestry were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and total serum bilirubin and other liver function tests was tested using linear regression, adjusting for age, gender, site, platform and ancestry principal components (PC).

Results

Consistent with previous results, a strong association signal has been detected for UGT1A gene cluster (best SNP rs887829, beta = 0.15, p = 1.30x10^-118) for total serum bilirubin level. Indeed, in this region more than 176 SNPs (or indels) had p<10⁻⁸ spanning 150Kb on the long arm of chromosome 2q37.1. In addition, we found a similar level of magnitude in a pediatric group (p = 8.26x10^-47, beta = 0.17). Further imputation using sequencing data as a reference panel revealed association of other markers including known TA7 repeat indels (rs8175347) (p = 9.78x10^-117) and rs111741722 (p = 5.41x10^-119) which were in proxy (r2 = 0.99) with rs887829. Among rare variants, two Asian subjects homozygous for coding SNP rs4148323 (G71R) were identified. Additional known effects for total serum bilirubin were also confirmed including organic anion transporters SLCO1B1-SLCO1B3, TDRP and ZMYND8 at FDR<0.05 with no gene-gene interaction effects. Phenome-wide association studies (PheWAS) suggest a protective effect of TA7 repeat against cerebrovascular disease in an adult cohort (OR = 0.75, p = 0.0008). Among other liver function tests, we also confirmed the previous effect of the ABO blood group locus for variation in serum alkaline phosphatase (rs579459, p = 9.44x10^-15).

Conclusions

Taken together, our data present interesting findings with strong confirmation of previous effects by simply using the eMERGE electronic health record phenotyping. In addition, our findings indicate that similar to the adult population, the UGT1A1 is the main locus responsible for normal variation of serum bilirubin in pediatric populations.     

Testing allele homogeneity: the problem of nested hypotheses

Background

Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas).

Results

We identified four genes with SNPs that showed the strongest and most consistent associations with LDL-C and CAD: EBF1, PPP2R2B, SPOCK1, and PRELID2. The most significant results for association of SNPs with LDL-C were: EBF1, rs6865969, p = 0.01; PPP2R2B, rs2125443, p = 0.005; SPOCK1, rs17600115, p = 0.003; and PRELID2, rs10074645, p = 0.0002). The most significant results for CAD were EBF1, rs6865969, p = 0.007; PPP2R2B, rs7736604, p = 0.0003; SPOCK1, rs17170899, p = 0.004; and PRELID2, rs7713855, p = 0.003.

Conclusion

Using an intermediate disease-related quantitative trait of LDL-C we have identified four novel CAD genes, EBF1, PRELID2, SPOCK1, and PPP2R2B. These four genes should be further examined in future functional studies as candidate susceptibility loci for cardiovascular disease mediated through LDL-cholesterol pathways.     

Increased prevalence of the CVD-associated ANRIL allele in the Roma/Gypsy population in comparison with the majority Czech population

Cardiovascular disease (CVD) is a major cause of death around the world, with highest prevalence reported in minority Roma/Gypsy populations living in developed countries. Whether these differences are caused by unhealthy lifestyles or genetic factors remain unknown. The aim of our study was to examine the genotype frequencies of the rs10757274 polymorphism in the 9p.21 locus within ANRIL (antisense non-coding RNA in the INK4 locus), a long non-coding RNA located in the vicinity of the CDKN2A/2B inhibitors loci. ANRIL is understood to be the strongest genetic determinant of CVD in Caucasians. Using PCR-RFLP, we analysed the ANRIL rs10757274 polymorphism in 298 non-Roma (50% male) and 302 Roma/Gypsy (50% male) adult (39.5 ± 15.1 years and 39.2 ± 12.8 years, respectively) subjects. We found that frequencies of the ANRIL GG, GA and AA genotypes were 20.1%, 52.4% and 27.5% in the majority population and 32.9%, 47.9% and 19.2% in Roma/Gypsy subjects, respectively. The distribution of genotypes was deemed significantly different at P < 0.001. Within the Roma/Gypsy population, we detected increased prevalence of the CVD-associated GG genotype. Increased prevalence of CVD among Roma/Gypsies subjects may be significantly linked to genetic background.     

A common variant highly associated with plasma VEGFA levels also contributes to the variation of both LDL-C and HDL-C

Vascular endothelial growth factor A (VEGFA) is among the most-significant stimulators of angiogenesis. Its effect on cardiovascular diseases and on the variation of related risk factors such as lipid parameters is considered important, although as yet unclear. Recently, we identified four common variants (rs6921438, rs4416670, rs6993770, and rs10738760) that explain up to 50% of the heritability of plasma VEGFA levels. In the present study, we aimed at assessing the contribution of these variants to the variation of blood lipid levels (including apoE, triglycerides, total cholesterol, low- and high-density lipoprotein cholesterol levels (LDL-C and HDL-C)] in healthy subjects. The effect of these single-nucleotide polymorphisms (SNPs) on lipid levels was assessed using linear regression in discovery and replication samples (n = 1,006 and n = 1,145; respectively), followed by a meta-analysis. Their gene×gene and gene×environment interactions were also assessed. SNP rs6921438 was associated with HDL-C (β = −0.08 mmol/l, P_overall = 1.2 × 10⁻⁷) and LDL-C (β = 0.13 mmol/l, P_overall = 1.5 × 10⁻⁴). We also identified a significant association between the interaction rs4416670×hypertension and apoE variation (P_overall = 1.7 × 10⁻⁵). Therefore, our present study shows a common genetic regulation between VEGFA and cholesterol homeostasis molecules. The SNP rs6921438 is in linkage disequilibrium with variants located in an enhancer- and promoter-associated histone mark region and could have a regulatory effect in the expression of surrounding genes, including VEGFA.     

Examining ERBB2 as a candidate gene for susceptibility to leprosy (Hansen’s disease) in Brazil

Leprosy remains prevalent in Brazil. ErbB2 is a receptor for leprosy bacilli entering Schwann cells, which mediates Mycobacterium leprae-induced demyelination and the ERBB2 gene lies within a leprosy susceptibility locus on chromosome 17q11-q21. To determine whether polymorphisms at the ERBB2 locus contribute to this linkage peak, three haplotype tagging single nucleotide polymorphisms (tag-SNPs) (rs2517956, rs2952156, rs1058808) were genotyped in 72 families (208 cases; 372 individuals) from the state of Pará (PA). All three tag-SNPs were associated with leprosy per se [best SNP rs2517959 odds ratio (OR) = 2.22; 95% confidence interval (CI) 1.37-3.59; p = 0.001]. Lepromatous (LL) (OR = 3.25; 95% CI 1.37-7.70; p = 0.007) and tuberculoid (TT) (OR = 1.79; 95% CI 1.04-3.05; p = 0.034) leprosy both contributed to the association, which is consistent with the previous linkage to chromosome 17q11-q21 in the population from PA and supports the functional role of ErbB2 in disease pathogenesis. To attempt to replicate these findings, six SNPs (rs2517955, rs2517956, rs1810132, rs2952156, rs1801200, rs1058808) were genotyped in a population-based sample of 570 leprosy cases and 370 controls from the state of Rio Grande do Norte (RN) and the results were analysed using logistic regression analysis. However, none of the associations were replicated in the RN sample, whether analysed for leprosy per se, LL leprosy, TT leprosy, erythema nodosum leprosum or reversal reaction conditions. The role of polymorphisms at ERBB2 in controlling susceptibility to leprosy in Brazil therefore remains unclear.     

A Regulatory MDM4 Genetic Variant Locating in the Binding Sequence of Multiple MicroRNAs Contributes to Susceptibility of Small Cell Lung Cancer

A functional rs4245739 A>C single nucleotide polymorphism (SNP) locating in the MDM43’-untranslated (3’-UTR) region creates a miR-191-5p or miR-887-3p targeting sites. This change results in decreased expression of oncogene MDM4. Therefore, we examined the association between this SNP and small cell lung cancer (SCLC) risk as well as its regulatory function in SCLC cells. Genotypes were determined in two independent case-control sets consisted of 520SCLC cases and 1040 controls from two regions of China. Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated by logistic regression. The impact of the rs4245739 SNP on miR-191-5p/miR-887-3p mediated MDM4 expression regulation was investigated using luciferase reporter gene assays. We found that the MDM4 rs4245739AC and CC genotypes were significantly associated with decreased SCLC susceptibility compared with the AA genotype in both case-control sets (Shandong set: OR = 0.53, 95% CI = 0.32–0.89, P = 0.014; Jiangsu set: OR = 0.47, 95% CI = 0.26–0.879, P = 0.017). Stratified analyses indicated that there was a significantly multiplicative interaction between rs4245739 and smoking (P _interactioin = 0.048). After co-tranfection of miRNAs and different allelic-MDM4 reporter constructs into SCLC cells, we found that the both miR-191-5p and miR-887-3p can lead to significantly decreased MDM4 expression activities in the construct with C-allelic 3’-UTR but not A-allelic 3’-UTR, suggesting a consistent genotype-phenotype correlation. Our data illuminate that the MDM4rs4245739SNP contributes to SCLC risk and support the notion that gene 3’-UTR genetic variants, impacting miRNA-binding, might modify SCLC susceptibility.     

The NOD2 p.Leu1007fsX1008 Mutation (rs2066847) Is a Stronger Predictor of the Clinical Course of Crohn's Disease than the FOXO3A Intron Variant rs12212067

Background

Very recently, a sub-analysis of genome-wide association scans revealed that the non-coding single nucleotide polymorphism (SNP) rs12212067 in the FOXO3A gene is associated with a milder course of Crohn''s disease (CD) (Cell 2013;155:57–69). The aim of our study was to evaluate the clinical value of the SNP rs12212067 in predicting the severity of CD by correlating CD patient genotype status with the most relevant complications of CD such as stenoses, fistulas, and CD-related surgery.

Methodology/Principal Findings

We genotyped 550 CD patients for rs12212067 (FOXO3A) and the three common CD-associated NOD2 mutations rs2066844, rs2066847, and rs2066847 and performed genotype-phenotype analyses.

Results

No significant phenotypic differences were found between the wild-type genotype TT of the FOXO3A SNP rs12212067 and the minor genotypes TG and GG independently from NOD2 variants. The allele frequency of the minor G allele was 12.7%. Age at diagnosis, disease duration, body mass index, surgery rate, stenoses, fistula, need for immunosuppressive therapy, and disease course were not significantly different. In contrast, the NOD2 mutant p.Leu1007fsX1008 (rs2066847) was highly associated with penetrating CD (p = 0.01), the development of fistulas (p = 0.01) and stenoses (p = 0.01), and ileal disease localization (p = 0.03). Importantly, the NOD2 SNP rs2066847 was a strong separator between an aggressive and a mild course of CD (p = 2.99×10⁻⁵), while the FOXO3A SNP rs12212067 did not separate between mild and aggressive CD behavior in our cohort (p = 0.35). 96.2% of the homozygous NOD2 p.Leu1007fsX1008 carriers had an aggressive disease behavior compared to 69.3% of the patients with the NOD2 wild-type genotype (p = 0.007).

Conclusion/Significance

In clinical practice, the NOD2 variant p.Leu1007fsX1008 (rs2066847), in particular in homozygous form, is a much stronger marker for a severe clinical phenotype than the FOXO3A rs12212067 SNP for a mild disease course on an individual patient level despite its important impact on the inflammatory response of monocytes.     

A promoter SNP rs4073T>A in the common allele of the interleukin 8 gene is associated with the development of idiopathic pulmonary fibrosis via the IL-8 protein enhancing mode

Background

Interleukin-8 (IL-8) is a potent chemo-attractant cytokine responsible for neutrophil infiltration in lungs with idiopathic pulmonary fibrosis (IPF). The IL-8 protein and mRNA expression are increased in the lung with IPF. We evaluated the effect of single nucleotide polymorphisms (SNPs) of the IL-8 gene on the risk of IPF.

Methods

One promoter (rs4073T>A) and two intronic SNPs (rs2227307T>G and rs2227306C>T) of the IL-8 genes were genotyped in 237 subjects with IPF and 456 normal controls. Logistic regression analysis was applied to evaluate the association of these SNPs with IPF. IL-8 in BAL fluids was measured using a quantitative sandwich enzyme immunoassay, and promoter activity was assessed using the luciferase reporter assay.

Results

The minor allele frequencies of rs4073T>A and rs2227307T>G were significantly lower in the 162 subjects with surgical biopsy-proven IPF and 75 subjects with clinical IPF compared with normal controls in the recessive model (OR = 0.46 and 0.48, p = 0.006 and 0.007, respectively). The IL-8 protein concentration in BAL fluids significantly increased in 24 subjects with IPF compared with 14 controls (p = 0.009). Nine IPF subjects homozygous for the rs4073 T>A common allele exhibited higher levels of the IL-8 protein compared with six subjects homozygous for the minor allele (p = 0.024). The luciferase activity of the rs4073T>A common allele was significantly higher than that of the rs4073T>A minor allele (p = 0.002).

Conclusion

The common allele of a promoter SNP, rs4073T>A, may increase susceptibility to the development of IPF via up-regulation of IL-8.     

Gender-dependent association of HSD11B1 single nucleotide polymorphisms with glucose and HDL-C levels

In this study, we investigated the influence of two SNPs (rs846910 and rs12086634) of the HSD11B1 gene that encodes 11β-hydroxysteroid dehydrogenase type 1(11β-HSD1), the enzyme that catalyzes the conversion of cortisol to cortisone, on variables associated with obesity and metabolic syndrome in 215 individuals of both sexes from southern Brazil. The HSD11B1 gene variants were genotyped using the TaqMan SNP genotyping assay. Glucose, triglycerides, total cholesterol, HDL-cholesterol and LDL-cholesterol were measured by standard automated methods. Significant results were found in women, with carriers of the G allele of SNP rs12086634 having higher glucose levels than non-carriers. Carriers of the A allele of SNP rs846910 had higher levels of HDL-cholesterol. The involvement of both polymorphisms as independent factors in determining the levels of glucose and HDL-cholesterol was confirmed by multiple regression analysis (β = 0.19 ±0.09, p = 0.03 and β= 0.22 ± 0.10, p = 0.03, respectively). Our findings suggest that the HSD11B1SNPs studied may indirectly influence glucose and HDL-cholesterol metabolism in women, possibly through down-regulation of the HSD11B1 gene by estrogen.     

Genome-Wide Meta-Analysis for Serum Calcium Identifies Significantly Associated SNPs near the Calcium-Sensing Receptor (CASR) Gene

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.