Genome-wide association studies on Northern Italy isolated populations provide further support concerning genetic susceptibility for major depressive disorder.

OBJECTIVES: Major depressive disorder (MDD) is a psychiatric disorder with pathogenesis influenced by both genetic and environmental factors. To date, the molecular-level understanding of its aetiology remains unclear. Thus, we aimed to identify genetic variants and susceptibility genes for MDD with a genome-wide association study (GWAS) approach. METHODS: We performed a meta-analysis of GWASs and a gene-based analysis on two Northern Italy isolated populations (cases/controls n = 166/472 and 33/320), followed by replication and polygenic risk score (PRS) analyses in Italian independent samples (cases n = 464, controls n = 339). RESULTS: We identified two novel MDD-associated genes, KCNQ5 (lead SNP rs867262, p = 3.82 × 10-9) and CTNNA2 (rs6729523, p = 1.25 × 10-8). The gene-based analysis revealed another six genes (p < 2.703 × 10-6): GRM7, CTNT4, SNRK, SRGAP3, TRAPPC9, and FHIT. No replication of the genome-wide significant SNPs was found in the independent cohort, even if 14 SNPs around CTNNA2 showed association with MDD and related phenotypes at the nominal level of p (<0.05). Furthermore, the PRS model developed in the discovery cohort discriminated cases and controls in the replication cohort. CONCLUSIONS: Our work suggests new possible genes associated with MDD, and the PRS analysis confirms the polygenic nature of this disorder. Future studies are required to better understand the role of these findings in MDD.

Case Report: Two cases of apparent discordance between non-invasive prenatal testing (NIPT) and amniocentesis resulting in feto-placental mosaicism of trisomy 21. Issues in diagnosis, investigation and counselling.

The sequencing of cell-free fetal DNA in the maternal plasma through non-invasive prenatal testing (NIPT) is an accurate genetic screening test to detect the most common fetal aneuploidies during pregnancy. The extensive use of NIPT, as a screening method, has highlighted the limits of the technique, including false positive and negative results. Feto-placental mosaicism is a challenging biological issue and is the most frequent cause of false positive and negative results in NIPT screening, and of discrepancy between NIPT and invasive test results. We are reporting on two cases of feto-placental mosaicism of trisomy 21, both with a low-risk NIPT result, identified by ultrasound signs and a subsequent amniocentesis consistent with a trisomy 21. In both cases, after the pregnancy termination, cytogenetic and/or cytogenomic analyses were performed on the placenta and fetal tissues, showing in the first case a mosaicism of trisomy 21 in both the placenta and the fetus, but a mosaicism in the placenta and a complete trisomy 21 in the fetus in the second case. These cases emphasize the need for accurate and complete pre-test NIPT counselling, as well as to identify situations at risk for a possible false negative NIPT result, which may underestimate a potential pathological condition, such as feto-placental mosaicism or fetal trisomy. Post-mortem molecular autopsy may discriminate between placental, fetal and feto-placental mosaicism, and between complete or mosaic fetal chromosomal anomalies. A multidisciplinary approach in counselling, as well as in the interpretation of biological events, is essential for the clarification of complex cases, such as feto-placental mosaicisms.

Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin.

BACKGROUND: Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown. METHODS: We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing. RESULTS: Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion. CONCLUSIONS: Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

Identification of a New Mutation in RSK2, the Gene for Coffin-Lowry Syndrome (CLS), in Two Related Patients with Mild and Atypical Phenotypes.

BACKGROUND: Coffin-Lowry syndrome (CLS) is a syndromic form of X-linked intellectual disability, in which specific associated facial, hand, and skeletal abnormalities are diagnostic features. METHODS: In the present study, an unreported missense genetic variant of the ribosomal S6 kinase 2 (RSK2) gene has been identified, by next-generation sequencing, in two related males with two different phenotypes of intellectual disability (ID) and peculiar facial dysmorphisms. We performed functional studies on this variant and another one, already reported in the literature, involving the same amino acid residue but, to date, without an efficient characterization. RESULTS: Our study demonstrated that the two variants involving residue 189 significantly impaired its kinase activity. CONCLUSIONS: We detected a loss-of-function RSK2 mutation with loss in kinase activity in a three-generation family with an X-linked ID.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction.

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Genetic determinants of circulating VEGF levels in major depressive disorder and electroconvulsive therapy response.

Alterations in peripheral vascular endothelial growth factor (VEGF) levels were observed in major depressive disorder and relative treatments and were shown to be influenced by genetic variants. The study objective was to explore, at a genome-wide level, possible interplaying effects between the genetic background and major depressive disorder in regulating VEGF levels. Moreover, we aimed to investigate the association between these variants and response to electroconvulsive therapy. A genome-wide association study was carried out both on controls and patients with major depressive disorder (n = 145; n = 121) in correlation with serum VEGF levels determined by ELISA. Five SNPs not included in SNP arrays were additionally genotyped. Seventy-one patients with treatment-resistant depression underwent electroconvulsive therapy and were evaluated as responders/nonresponders. An association between VEGF levels and a locus in 6p21.1, downstream the VEGF gene, was evidenced both in controls (best SNP: FDR-corrected p = 2.4 × 10-5 ) and in patients with major depressive disorder (best SNP: FDR-corrected p = 2.6 × 10-3 ). The alleles associated with lower VEGF concentrations in patients were also associated with nonresponse to electroconvulsive therapy (p = .01). These results confirm a role of SNPs in 6p21.1 locus as major influencers of circulating VEGF levels also in patients affected by major depressive disorder and indicate a possible implication in response to electroconvulsive therapy.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

Genome-Wide Meta-Analysis Unravels Interactions between Magnesium Homeostasis and Metabolic Phenotypes.

Magnesium (Mg2+) homeostasis is critical for metabolism. However, the genetic determinants of the renal handling of Mg2+, which is crucial for Mg2+ homeostasis, and the potential influence on metabolic traits in the general population are unknown. We obtained plasma and urine parameters from 9099 individuals from seven cohorts, and conducted a genome-wide meta-analysis of Mg2+ homeostasis. We identified two loci associated with urinary magnesium (uMg), rs3824347 (P=4.4×10-13) near TRPM6, which encodes an epithelial Mg2+ channel, and rs35929 (P=2.1×10-11), a variant of ARL15, which encodes a GTP-binding protein. Together, these loci account for 2.3% of the variation in 24-hour uMg excretion. In human kidney cells, ARL15 regulated TRPM6-mediated currents. In zebrafish, dietary Mg2+ regulated the expression of the highly conserved ARL15 ortholog arl15b, and arl15b knockdown resulted in renal Mg2+ wasting and metabolic disturbances. Finally, ARL15 rs35929 modified the association of uMg with fasting insulin and fat mass in a general population. In conclusion, this combined observational and experimental approach uncovered a gene-environment interaction linking Mg2+ deficiency to insulin resistance and obesity.

Common variants in CLDN14 are associated with differential excretion of magnesium over calcium in urine.

The nature and importance of genetic factors regulating the differential handling of Ca2+ and Mg2+ by the renal tubule in the general population are poorly defined. We conducted a genome-wide meta-analysis of urinary magnesium-to-calcium ratio to identify associated common genetic variants. We included 9320 adults of European descent from four genetic isolates and three urban cohorts. Urinary magnesium and calcium concentrations were measured centrally in spot urine, and each study conducted linear regression analysis of urinary magnesium-to-calcium ratio on ~2.5 million single-nucleotide polymorphisms (SNPs) using an additive model. We investigated, in mouse, the renal expression profile of the top candidate gene and its variation upon changes in dietary magnesium. The genome-wide analysis evidenced a top locus (rs172639, p = 1.7 × 10-12), encompassing CLDN14, the gene coding for claudin-14, that was genome-wide significant when using urinary magnesium-to-calcium ratio, but not either one taken separately. In mouse, claudin-14 is expressed in the distal nephron segments specifically handling magnesium, and its expression is regulated by chronic changes in dietary magnesium content. A genome-wide approach identified common variants in the CLDN14 gene exerting a robust influence on the differential excretion of Mg2+ over Ca2+ in urine. These data highlight the power of urinary electrolyte ratios to unravel genetic determinants of renal tubular function. Coupled with mouse experiments, these results support a major role for claudin-14, a gene associated with kidney stones, in the differential paracellular handling of divalent cations by the renal tubule.

52 Genetic Loci Influencing Myocardial Mass.

BACKGROUND: Myocardial mass is a key determinant of cardiac muscle function and hypertrophy. Myocardial depolarization leading to cardiac muscle contraction is reflected by the amplitude and duration of the QRS complex on the electrocardiogram (ECG). Abnormal QRS amplitude or duration reflect changes in myocardial mass and conduction, and are associated with increased risk of heart failure and death. OBJECTIVES: This meta-analysis sought to gain insights into the genetic determinants of myocardial mass. METHODS: We carried out a genome-wide association meta-analysis of 4 QRS traits in up to 73,518 individuals of European ancestry, followed by extensive biological and functional assessment. RESULTS: We identified 52 genomic loci, of which 32 are novel, that are reliably associated with 1 or more QRS phenotypes at p < 1 × 10(-8). These loci are enriched in regions of open chromatin, histone modifications, and transcription factor binding, suggesting that they represent regions of the genome that are actively transcribed in the human heart. Pathway analyses provided evidence that these loci play a role in cardiac hypertrophy. We further highlighted 67 candidate genes at the identified loci that are preferentially expressed in cardiac tissue and associated with cardiac abnormalities in Drosophila melanogaster and Mus musculus. We validated the regulatory function of a novel variant in the SCN5A/SCN10A locus in vitro and in vivo. CONCLUSIONS: Taken together, our findings provide new insights into genes and biological pathways controlling myocardial mass and may help identify novel therapeutic targets.

Brain-derived neurotrophic factor serum levels in genetically isolated populations: gender-specific association with anxiety disorder subtypes but not with anxiety levels or Val66Met polymorphism.

Anxiety disorders (ADs) are disabling chronic disorders with exaggerated behavioral response to threats. This study was aimed at testing the hypothesis that ADs may be associated with reduced neurotrophic activity, particularly of Brain-derived neurotrophic factor (BDNF), and determining possible effects of genetics on serum BDNF concentrations. In 672 adult subjects from six isolated villages in North-Eastern Italy with high inbreeding, we determined serum BDNF levels and identified subjects with different ADs subtypes such as Social and Specific Phobias (PHSOC, PHSP), Generalized Anxiety Disorder (GAD), and Panic Disorder (PAD). Analysis of the population as a whole or individual village showed no significant correlation between serum BDNF levels and Val66Met polymorphism and no association with anxiety levels. Stratification of subjects highlighted a significant decrease in serum BDNF in females with GAD and males with PHSP. This study indicates low heritability and absence of any impact of the Val66Met polymorphism on circulating concentrations of BDNF. Our results show that BDNF is not a general biomarker of anxiety but serum BDNF levels correlate in a gender-specific manner with ADs subtypes.

Analysis of functional variants reveals new candidate genes associated with alexithymia.

In this study we explored the possible association between 36,915 functional variants and alexithymia, a personality trait characterized by the inability to identify and describe emotions and feelings. From our analysis, variants in the genes ABCB4, TP53AIP1, ARHGAP32 and TMEM88B were identified linked to the alexithymia phenotype.

Modulation of genetic associations with serum urate levels by body-mass-index in humans.

We tested for interactions between body mass index (BMI) and common genetic variants affecting serum urate levels, genome-wide, in up to 42569 participants. Both stratified genome-wide association (GWAS) analyses, in lean, overweight and obese individuals, and regression-type analyses in a non BMI-stratified overall sample were performed. The former did not uncover any novel locus with a major main effect, but supported modulation of effects for some known and potentially new urate loci. The latter highlighted a SNP at RBFOX3 reaching genome-wide significant level (effect size 0.014, 95% CI 0.008-0.02, Pinter= 2.6 x 10-8). Two top loci in interaction term analyses, RBFOX3 and ERO1LB-EDARADD, also displayed suggestive differences in main effect size between the lean and obese strata. All top ranking loci for urate effect differences between BMI categories were novel and most had small magnitude but opposite direction effects between strata. They include the locus RBMS1-TANK (men, Pdifflean-overweight= 4.7 x 10-8), a region that has been associated with several obesity related traits, and TSPYL5 (men, Pdifflean-overweight= 9.1 x 10-8), regulating adipocytes-produced estradiol. The top-ranking known urate loci was ABCG2, the strongest known gout risk locus, with an effect halved in obese compared to lean men (Pdifflean-obese= 2 x 10-4). Finally, pathway analysis suggested a role for N-glycan biosynthesis as a prominent urate-associated pathway in the lean stratum. These results illustrate a potentially powerful way to monitor changes occurring in obesogenic environment.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization.

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal mendelian long-QT syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals, we identified 35 common variant loci associated with QT interval that collectively explain ∼8-10% of QT-interval variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 new QT interval-associated loci in 298 unrelated probands with LQTS identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies new candidate genes for ventricular arrhythmias, LQTS and SCD.

A general approach for haplotype phasing across the full spectrum of relatedness.

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Insight into genetic determinants of resting heart rate.

BACKGROUND: Recent studies suggested that resting heart rate (RHR) might be an independent predictor of cardiovascular mortality and morbidity. Nonetheless, the interrelation between RHR and cardiovascular diseases is not clear. In order to resolve this puzzle, the importance of genetic determinants of RHR has been recently suggested, but it needs to be further investigated. OBJECTIVE: The aim of this study was to estimate the contribution of common genetic variations on RHR using Genome Wide Association Study. METHODS: We performed a Genome Wide Association Study in an isolated population cohort of 1737 individuals, the Italian Network on Genetic Isolates - Friuli Venezia Giulia (INGI-FVG). Moreover, a haplotype analysis was performed. A regression tree analysis was run to highlight the effect of each haplotype combination on the phenotype. RESULTS: A significant level of association (p<5 × 10(-8)) was detected for Single Nucleotide Polymorphisms (SNPs) in two genes expressed in the heart: MAML1 and CANX. Founding that the three different variants of the haplotype, which encompass both genes, yielded a phenotypic correlation. Indeed, a haplotype in homozygosity is significantly associated with the lower quartile of RHR (RHR ≤ 58 bpm). Moreover no significant association was found between cardiovascular risk factors and the different haplotype combinations. CONCLUSION: Mastermind-like 1 and Calnexin were found to be associated with RHR. We demonstrated a relation between a haplotype and the lower quartile of RHR in our populations. Our findings highlight that genetic determinants of RHR may be implicated in determining cardiovascular diseases and could allow a better risk stratification.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 µg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Common variants in Mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations.

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤ 21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.