A genome-wide association study of female sexual dysfunction

Background

Female sexual dysfunction (FSD) is an important but controversial problem with serious negative impact on women’s quality of life. Data from twin studies have shown a genetic contribution to the development and maintenance of FSD.

Methodology/Principal Findings

We performed a genome-wide association study (GWAS) on 2.5 million single-nucleotide polymorphisms (SNPs) in 1,104 female twins (25–81 years of age) in a population-based register and phenotypic data on lifelong sexual functioning. Although none reached conventional genome-wide level of significance (10×-8), we found strongly suggestive associations with the phenotypic dimension of arousal (rs13202860, P = 1.2×10⁻⁷; rs1876525, P = 1.2×10⁻⁷; and rs13209281 P = 8.3×10⁻⁷) on chromosome 6, around 500kb upstream of the locus HTR1E (5-hydroxytryptamine receptor 1E) locus, related to the serotonin brain pathways. We could not replicate previously reported candidate SNPs associated with FSD in the DRD4, 5HT2A and IL-1B loci.

Conclusions/Significance

We report the first GWAS of FSD symptoms in humans. This has pointed to several “risk alleles” and the implication of the serotonin and GABA pathways. Ultimately, understanding key mechanisms via this research may lead to new FSD treatments and inform clinical practice and developments in psychiatric nosology.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.     

Genome-wide association study identifies genetic variation in neurocan as a susceptibility factor for bipolar disorder

We conducted a genome-wide association study (GWAS) and a follow-up study of bipolar disorder (BD), a common neuropsychiatric disorder. In the GWAS, we investigated 499,494 autosomal and 12,484 X-chromosomal SNPs in 682 patients with BD and in 1300 controls. In the first follow-up step, we tested the most significant 48 SNPs in 1729 patients with BD and in 2313 controls. Eight SNPs showed nominally significant association with BD and were introduced to a meta-analysis of the GWAS and the first follow-up samples. Genetic variation in the neurocan gene (NCAN) showed genome-wide significant association with BD in 2411 patients and 3613 controls (rs1064395, p = 3.02 × 10⁻⁸; odds ratio = 1.31). In a second follow-up step, we replicated this finding in independent samples of BD, totaling 6030 patients and 31,749 controls (p = 2.74 × 10⁻⁴; odds ratio = 1.12). The combined analysis of all study samples yielded a p value of 2.14 × 10⁻⁹ (odds ratio = 1.17). Our results provide evidence that rs1064395 is a common risk factor for BD. NCAN encodes neurocan, an extracellular matrix glycoprotein, which is thought to be involved in cell adhesion and migration. We found that expression in mice is localized within cortical and hippocampal areas. These areas are involved in cognition and emotion regulation and have previously been implicated in BD by neuropsychological, neuroimaging, and postmortem studies.     

Evaluation of Linkage Disequilibrium Pattern and Association Study on Seed Oil Content in Brassica napus Using ddRAD Sequencing

High-density genetic markers are the prerequisite for understanding linkage disequilibrium (LD) and genome-wide association studies (GWASs) of complex traits in crops. To evaluate the LD pattern in oilseed rape, we sequenced a previous association panel containing 189 B. napus inbred lines using double-digested restriction-site associated DNA (ddRAD) and genotyped 19,327 RAD tags. A total of 15,921 RAD tags were assigned to a published genetic linkage map and the majority (71.1%) of these tags was uniquely mapped to the draft reference genome “Darmor-bzh.” The distance of LD decay was 1,214 kb across the genome at the background level (r² = 0.26), with the distances of LD decay being 405 kb and 2,111 kb in the A and C subgenomes, respectively. A total of 361 haplotype blocks with length > 100 kb were identified in the entire genome. The association panel could be classified into two groups, P1 and P2, which are essentially consistent with the geographical origins of varieties. A large number of group-specific haplotypes were identified, reflecting that varieties in the P1 and P2 groups experienced distinct selection in breeding programs to adapt their different growth habitats. GWAS repeatedly detected two loci significantly associated with oil content of seeds based on the developed SNPs, suggesting that the high-density SNPs were useful for understanding the genetic determinants of complex traits in GWAS.     

Genome-wide association study identifies HLA-DP as a susceptibility gene for pediatric asthma in Asian populations

Asthma is a complex phenotype influenced by genetic and environmental factors. We conducted a genome-wide association study (GWAS) with 938 Japanese pediatric asthma patients and 2,376 controls. Single-nucleotide polymorphisms (SNPs) showing strong associations (P<1×10⁻⁸) in GWAS were further genotyped in an independent Japanese samples (818 cases and 1,032 controls) and in Korean samples (835 cases and 421 controls). SNP rs987870, located between HLA-DPA1 and HLA-DPB1, was consistently associated with pediatric asthma in 3 independent populations (P _combined = 2.3×10⁻¹⁰, odds ratio [OR] = 1.40). HLA-DP allele analysis showed that DPA1*0201 and DPB1*0901, which were in strong linkage disequilibrium, were strongly associated with pediatric asthma (DPA1*0201: P = 5.5×10⁻¹⁰, OR = 1.52, and DPB1*0901: P = 2.0×10⁻⁷, OR = 1.49). Our findings show that genetic variants in the HLA-DP locus are associated with the risk of pediatric asthma in Asian populations.     

Genome-wide Association Study Reveals Multiple Nasopharyngeal Carcinoma-Associated Loci within the HLA Region at Chromosome 6p21.3

Nasopharyngeal carcinoma (NPC) is a multifactorial malignancy closely associated with genetic factors and Epstein-Barr virus infection. To identify the common genetic variants linked to NPC susceptibility, we conducted a genome-wide association study (GWAS) in 277 NPC patients and 285 healthy controls within the Taiwanese population, analyzing 480,365 single-nucleotide polymorphisms (SNPs). Twelve statistically significant SNPs were identified and mapped to chromosome 6p21.3. Associations were replicated in two independent sets of case-control samples. Two of the most significant SNPs (rs2517713 and rs2975042; p_combined = 3.9 × 10⁻²⁰ and 1.6 × 10⁻¹⁹, respectively) were located in the HLA-A gene. Moreover, we detected significant associations between NPC and two genes: specifically, gamma aminobutyric acid b receptor 1 (GABBR1) (rs29232; p_combined = 8.97 × 10⁻¹⁷) and HLA-F (rs3129055 and rs9258122; p_combined = 7.36 × 10⁻¹¹ and 3.33 × 10⁻¹⁰, respectively). Notably, the association of rs29232 remained significant (residual p < 5 × 10⁻⁴) after adjustment for age, gender, and HLA-related SNPs. Furthermore, higher GABA_B receptor 1 expression levels can be found in the tumor cells in comparison to the adjacent epithelial cells (p < 0.001) in NPC biopsies, implying a biological role of GABBR1 in NPC carcinogenesis. To our knowledge, it is the first GWAS report of NPC showing that multiple loci (HLA-A, HLA-F, and GABBR1) within chromosome 6p21.3 are associated with NPC. Although some of these relationships may be attributed to linkage disequilibrium between the loci, the findings clearly provide a fresh direction for the study of NPC development.     

A systems genetics approach provides a bridge from discovered genetic variants to biological pathways in rheumatoid arthritis

Genome-wide association studies (GWAS) have yielded novel genetic loci underlying common diseases. We propose a systems genetics approach to utilize these discoveries for better understanding of the genetic architecture of rheumatoid arthritis (RA). Current evidence of genetic associations with RA was sought through PubMed and the NHGRI GWAS catalog. The associations of 15 single nucleotide polymorphisms and HLA-DRB1 alleles were confirmed in 1,287 cases and 1,500 controls of Japanese subjects. Among these, HLA-DRB1 alleles and eight SNPs showed significant associations and all but one of the variants had the same direction of effect as identified in the previous studies, indicating that the genetic risk factors underlying RA are shared across populations. By receiver operating characteristic curve analysis, the area under the curve (AUC) for the genetic risk score based on the selected variants was 68.4%. For seropositive RA patients only, the AUC improved to 70.9%, indicating good but suboptimal predictive ability. A simulation study shows that more than 200 additional loci with similar effect size as recent GWAS findings or 20 rare variants with intermediate effects are needed to achieve AUC = 80.0%. We performed the random walk with restart (RWR) algorithm to prioritize genes for future mapping studies. The performance of the algorithm was confirmed by leave-one-out cross-validation. The RWR algorithm pointed to ZAP70 in the first rank, in which mutation causes RA-like autoimmune arthritis in mice. By applying the hierarchical clustering method to a subnetwork comprising RA-associated genes and top-ranked genes by the RWR, we found three functional modules relevant to RA etiology: “leukocyte activation and differentiation”, “pattern-recognition receptor signaling pathway”, and “chemokines and their receptors”.These results suggest that the systems genetics approach is useful to find directions of future mapping strategies to illuminate biological pathways.     

Network-assisted investigation of combined causal signals from genome-wide association studies in schizophrenia

With the recent success of genome-wide association studies (GWAS), a wealth of association data has been accomplished for more than 200 complex diseases/traits, proposing a strong demand for data integration and interpretation. A combinatory analysis of multiple GWAS datasets, or an integrative analysis of GWAS data and other high-throughput data, has been particularly promising. In this study, we proposed an integrative analysis framework of multiple GWAS datasets by overlaying association signals onto the protein-protein interaction network, and demonstrated it using schizophrenia datasets. Building on a dense module search algorithm, we first searched for significantly enriched subnetworks for schizophrenia in each single GWAS dataset and then implemented a discovery-evaluation strategy to identify module genes with consistent association signals. We validated the module genes in an independent dataset, and also examined them through meta-analysis of the related SNPs using multiple GWAS datasets. As a result, we identified 205 module genes with a joint effect significantly associated with schizophrenia; these module genes included a number of well-studied candidate genes such as DISC1, GNA12, GNA13, GNAI1, GPR17, and GRIN2B. Further functional analysis suggested these genes are involved in neuronal related processes. Additionally, meta-analysis found that 18 SNPs in 9 module genes had P _meta<1×10⁻⁴, including the gene HLA-DQA1 located in the MHC region on chromosome 6, which was reported in previous studies using the largest cohort of schizophrenia patients to date. These results demonstrated our bi-directional network-based strategy is efficient for identifying disease-associated genes with modest signals in GWAS datasets. This approach can be applied to any other complex diseases/traits where multiple GWAS datasets are available.     

Genome-Wide Association Study of Gene by Smoking Interactions in Coronary Artery Calcification

Many GWAS have identified novel loci associated with common diseases, but have focused only on main effects of individual genetic variants rather than interactions with environmental factors (GxE). Identification of GxE interactions is particularly important for coronary heart disease (CHD), a major preventable source of morbidity and mortality with strong non-genetic risk factors. Atherosclerosis is the major cause of CHD, and coronary artery calcification (CAC) is directly correlated with quantity of coronary atherosclerotic plaque. In the current study, we tested for genetic variants influencing extent of CAC via interaction with smoking (GxS), by conducting a GxS discovery GWAS in Genetic Epidemiology Network of Arteriopathy (GENOA) sibships (N = 915 European Americans) followed by replication in Framingham Heart Study (FHS) sibships (N = 1025 European Americans). Generalized estimating equations accounted for the correlation within sibships in strata-specific groups of smokers and nonsmokers, as well as GxS interaction. Primary analysis found SNPs that showed suggestive associations (p≤10⁻⁵) in GENOA GWAS, but these index SNPs did not replicate in FHS. However, secondary analysis was able to replicate candidate gene regions in FHS using other SNPs (+/−250 kb of GENOA index SNP). In smoker and nonsmoker groups, replicated genes included TCF7L2 (p = 6.0×10⁻⁵) and WWOX (p = 4.5×10⁻⁶); and TNFRSF8 (p = 7.8×10⁻⁵), respectively. For GxS interactions, replicated genes included TBC1D4 (p = 6.9×10⁻⁵) and ADAMTS9 (P = 7.1×10⁻⁵). Interestingly, these genes are involved in inflammatory pathways mediated by the NF-κB axis. Since smoking is known to induce chronic and systemic inflammation, association of these genes likely reflects roles in CAC development via inflammatory pathways. Furthermore, the NF-κB axis regulates bone remodeling, a key physiological process in CAC development. In conclusion, GxS GWAS has yielded evidence for novel loci that are associated with CAC via interaction with smoking, providing promising new targets for future population-based and functional studies of CAC development.     

A novel, functional and replicable risk gene region for alcohol dependence identified by genome-wide association study

Several genome-wide association studies (GWASs) reported tens of risk genes for alcohol dependence, but most of them have not been replicated or confirmed by functional studies. The present study used a GWAS to search for novel, functional and replicable risk gene regions for alcohol dependence. Associations of all top-ranked SNPs identified in a discovery sample of 681 African-American (AA) cases with alcohol dependence and 508 AA controls were retested in a primary replication sample of 1,409 European-American (EA) cases and 1,518 EA controls. The replicable associations were then subjected to secondary replication in a sample of 6,438 Australian family subjects. A functional expression quantitative trait locus (eQTL) analysis of these replicable risk SNPs was followed-up in order to explore their cis-acting regulatory effects on gene expression. We found that within a 90 Mb region around PHF3-PTP4A1 locus in AAs, a linkage disequilibrium (LD) block in PHF3-PTP4A1 formed the only peak associated with alcohol dependence at p<10⁻⁴. Within this block, 30 SNPs associated with alcohol dependence in AAs (1.6×10⁻⁵≤p≤0.050) were replicated in EAs (1.3×10⁻³≤p≤0.038), and 18 of them were also replicated in Australians (1.8×10⁻³≤p≤0.048). Most of these risk SNPs had strong cis-acting regulatory effects on PHF3-PTP4A1 mRNA expression across three HapMap samples. The distributions of −log(p) values for association and functional signals throughout this LD block were highly consistent across AAs, EAs, Australians and three HapMap samples. We conclude that the PHF3-PTP4A1 region appears to harbor a causal locus for alcohol dependence, and proteins encoded by PHF3 and/or PTP4A1 might play a functional role in the disorder.     

Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-β activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10⁻¹² in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85–0.94) and 0.80 (0.75–0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10⁻⁶ from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10⁻⁷), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively.     

Extensive natural variation for cellular hydrogen peroxide release is genetically controlled

Natural variation in DNA sequence contributes to individual differences in quantitative traits. While multiple studies have shown genetic control over gene expression variation, few additional cellular traits have been investigated. Here, we investigated the natural variation of NADPH oxidase-dependent hydrogen peroxide (H₂O₂ release), which is the joint effect of reactive oxygen species (ROS) production, superoxide metabolism and degradation, and is related to a number of human disorders. We assessed the normal variation of H₂O₂ release in lymphoblastoid cell lines (LCL) in a family-based 3-generation cohort (CEPH-HapMap), and in 3 population-based cohorts (KORA, GenCord, HapMap). Substantial individual variation was observed, 45% of which were associated with heritability in the CEPH-HapMap cohort. We identified 2 genome-wide significant loci of Hsa12 and Hsa15 in genome-wide linkage analysis. Next, we performed genome-wide association study (GWAS) for the combined KORA-GenCord cohorts (n = 279) using enhanced marker resolution by imputation (>1.4 million SNPs). We found 5 significant associations (p<5.00×10−8) and 54 suggestive associations (p<1.00×10−5), one of which confirmed the linked region on Hsa15. To replicate our findings, we performed GWAS using 58 HapMap individuals and ∼2.1 million SNPs. We identified 40 genome-wide significant and 302 suggestive SNPs, and confirmed genome signals on Hsa1, Hsa12, and Hsa15. Genetic loci within 900 kb from the known candidate gene p67phox on Hsa1 were identified in GWAS in both cohorts. We did not find replication of SNPs across all cohorts, but replication within the same genomic region. Finally, a highly significant decrease in H₂O₂ release was observed in Down Syndrome (DS) individuals (p<2.88×10−12). Taken together, our results show strong evidence of genetic control of H₂O₂ in LCL of healthy and DS cohorts and suggest that cellular phenotypes, which themselves are also complex, may be used as proxies for dissection of complex disorders.     

Genome wide association (GWA) study for early onset extreme obesity supports the role of fat mass and obesity associated gene (FTO) variants

Background

Obesity is a major health problem. Although heritability is substantial, genetic mechanisms predisposing to obesity are not very well understood. We have performed a genome wide association study (GWA) for early onset (extreme) obesity.

Methodology/Principal Findings

a) GWA (Genome-Wide Human SNP Array 5.0 comprising 440,794 single nucleotide polymorphisms) for early onset extreme obesity based on 487 extremely obese young German individuals and 442 healthy lean German controls; b) confirmatory analyses on 644 independent families with at least one obese offspring and both parents. We aimed to identify and subsequently confirm the 15 SNPs (minor allele frequency ≥10%) with the lowest p-values of the GWA by four genetic models: additive, recessive, dominant and allelic. Six single nucleotide polymorphisms (SNPs) in FTO (fat mass and obesity associated gene) within one linkage disequilibrium (LD) block including the GWA SNP rendering the lowest p-value (rs1121980; log-additive model: nominal p = 1.13×10⁻⁷, corrected p = 0.0494; odds ratio (OR)_CT 1.67, 95% confidence interval (CI) 1.22–2.27; OR_TT 2.76, 95% CI 1.88–4.03) belonged to the 15 SNPs showing the strongest evidence for association with obesity. For confirmation we genotyped 11 of these in the 644 independent families (of the six FTO SNPs we chose only two representing the LD bock). For both FTO SNPs the initial association was confirmed (both Bonferroni corrected p<0.01). However, none of the nine non-FTO SNPs revealed significant transmission disequilibrium.

Conclusions/Significance

Our GWA for extreme early onset obesity substantiates that variation in FTO strongly contributes to early onset obesity. This is a further proof of concept for GWA to detect genes relevant for highly complex phenotypes. We concurrently show that nine additional SNPs with initially low p-values in the GWA were not confirmed in our family study, thus suggesting that of the best 15 SNPs in the GWA only the FTO SNPs represent true positive findings.     

Fine-Mapping the HOXB Region Detects Common Variants Tagging a Rare Coding Allele: Evidence for Synthetic Association in Prostate Cancer

The HOXB13 gene has been implicated in prostate cancer (PrCa) susceptibility. We performed a high resolution fine-mapping analysis to comprehensively evaluate the association between common genetic variation across the HOXB genetic locus at 17q21 and PrCa risk. This involved genotyping 700 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of 3195 SNPs in 20,440 PrCa cases and 21,469 controls in The PRACTICAL consortium. We identified a cluster of highly correlated common variants situated within or closely upstream of HOXB13 that were significantly associated with PrCa risk, described by rs117576373 (OR 1.30, P = 2.62×10⁻¹⁴). Additional genotyping, conditional regression and haplotype analyses indicated that the newly identified common variants tag a rare, partially correlated coding variant in the HOXB13 gene (G84E, rs138213197), which has been identified recently as a moderate penetrance PrCa susceptibility allele. The potential for GWAS associations detected through common SNPs to be driven by rare causal variants with higher relative risks has long been proposed; however, to our knowledge this is the first experimental evidence for this phenomenon of synthetic association contributing to cancer susceptibility.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

Association of JAG1 with Bone Mineral Density and Osteoporotic Fractures: A Genome-wide Association Study and Follow-up Replication Studies

Bone mineral density (BMD), a diagnostic parameter for osteoporosis and a clinical predictor of fracture, is a polygenic trait with high heritability. To identify genetic variants that influence BMD in different ethnic groups, we performed a genome-wide association study (GWAS) on 800 unrelated Southern Chinese women with extreme BMD and carried out follow-up replication studies in six independent study populations of European descent and Asian populations including 18,098 subjects. In the meta-analysis, rs2273061 of the Jagged1 (JAG1) gene was associated with high BMD (p = 5.27 × 10⁻⁸ for lumbar spine [LS] and p = 4.15 × 10⁻⁵ for femoral neck [FN], n = 18,898). This SNP was further found to be associated with the low risk of osteoporotic fracture (p = 0.009, OR = 0.7, 95% CI 0.57–0.93, n = 1881). Region-wide and haplotype analysis showed that the strongest association evidence was from the linkage disequilibrium block 5, which included rs2273061 of the JAG1 gene (p = 8.52 × 10⁻⁹ for LS and 3.47 × 10⁻⁵ at FN). To assess the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of c-Myc to the “G” but not “A” allele of rs2273061. A mRNA expression study in both human bone-derived cells and peripheral blood mononuclear cells confirmed association of the high BMD-related allele G of rs2273061 with higher JAG1 expression. Our results identify the JAG1 gene as a candidate for BMD regulation in different ethnic groups, and it is a potential key factor for fracture pathogenesis.