High SHBG and Low Bioavailable Testosterone are Strongly Causally Associated with Increased Forearm Fracture Risk in Women: An MR Study Leveraging Novel Female-Specific Data.

The effects of androgens on women's bone health are not fully understood. Mendelian randomization (MR) studies using sex-combined data suggest that sex hormone-binding globulin (SHBG) and bioavailable testosterone (BioT) causally affect bone traits. Given significant sex differences in hormone regulation and effects, female-specific MR studies are necessary. In the current study, we explored the causal relationships of SHBG, BioT, and total testosterone (TT) with forearm fracture (FAFx) risk in women using two-sample MR analyses. We utilized a unique female-specific FAFx outcome dataset from three European biobanks (UFO, HUNT, Estonian Biobank) comprising 111,351 women and 8823 FAFx cases, along with female-specific genetic instruments of SHBG, BioT, and TT identified in the UK Biobank. We also assessed bone mineral density (BMD) at the forearm (FA), femoral neck (FN), and lumbar spine (LS) using female-specific GWAS data from the GEFOS consortium. High SHBG (odds ratio per standard deviation increase (OR/SD): 1.53, 95% confidence intervals (CIs): 1.34-1.75), low BioT (OR/SD: 0.77, 0.71-0.84) and low TT (OR/SD 0.90, 0.83-0.98) were causally associated with increased FAFx risk. BioT was positively, and SHBG inversely, causally associated with especially FA-BMD, but also LS-BMD and FN-BMD, while TT was only significantly positively associated with FA-BMD and LS-BMD. We propose that endogenous androgens and SHBG are important for women's bone health at distal trabecular-rich bone sites such as the distal forearm and may serve as predictors for FAFx risk.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Correction: Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

[This corrects the article DOI: 10.1371/journal.pgen.1006866.].

Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

A small percentage of women with cervical HPV infection progress to cervical neoplasia, and the risk factors determining progression are incompletely understood. We sought to define the genetic loci involved in cervical neoplasia and to assess its heritability using unbiased unrelated case/control statistical approaches. We demonstrated strong association of cervical neoplasia with risk and protective HLA haplotypes that are determined by the amino-acids carried at positions 13 and 71 in pocket 4 of HLA-DRB1 and position 156 in HLA-B. Furthermore, 36% (standard error 2.4%) of liability of HPV-associated cervical pre-cancer and cancer is determined by common genetic variants. Women in the highest 10% of genetic risk scores have approximately >7.1% risk, and those in the highest 5% have approximately >21.6% risk, of developing cervical neoplasia. Future studies should examine genetic risk prediction in assessing the risk of cervical neoplasia further, in combination with other screening methods.

HLA and KIR Associations of Cervical Neoplasia.

Background: Cervical cancer is the fourth most common cancer in women, and we recently reported human leukocyte antigen (HLA) alleles showing strong associations with cervical neoplasia risk and protection. HLA ligands are recognized by killer immunoglobulin-like receptors (KIRs) expressed on a range of immune cell subsets, governing their proinflammatory activity. We hypothesized that the inheritance of particular HLA-KIR combinations would increase cervical neoplasia risk. Methods: Here, we used HLA and KIR dosages imputed from single-nucleotide polymorphism genotype data from 2143 cervical neoplasia cases and 13858 healthy controls of European decent. Results: The following 4 novel HLA alleles were identified in association with cervical neoplasia, owing to their linkage disequilibrium with known cervical neoplasia-associated HLA-DRB1 alleles: HLA-DRB3*9901 (odds ratio [OR], 1.24; P = 2.49 × 10-9), HLA-DRB5*0101 (OR, 1.29; P = 2.26 × 10-8), HLA-DRB5*9901 (OR, 0.77; P = 1.90 × 10-9), and HLA-DRB3*0301 (OR, 0.63; P = 4.06 × 10-5). We also found that homozygosity of HLA-C1 group alleles is a protective factor for human papillomavirus type 16 (HPV16)-related cervical neoplasia (C1/C1; OR, 0.79; P = .005). This protective association was restricted to carriers of either KIR2DL2 (OR, 0.67; P = .00045) or KIR2DS2 (OR, 0.69; P = .0006). Conclusions: Our findings suggest that HLA-C1 group alleles play a role in protecting against HPV16-related cervical neoplasia, mainly through a KIR-mediated mechanism.

Pre-diagnostic vitamin D concentrations and cancer risks in older individuals: an analysis of cohorts participating in the CHANCES consortium.

The associations of circulating 25-hydroxyvitamin D [25(OH)D] concentrations with total and site-specific cancer incidence have been examined in several epidemiological studies with overall inconclusive findings. Very little is known about the association of vitamin D with cancer incidence in older populations. We assessed the association of pre-diagnostic serum 25(OH)D levels with incidence of all cancers combined and incidence of lung, colorectal, breast, prostate and lymphoid malignancies among older adults. Pre-diagnostic 25(OH)D concentrations and cancer incidence were available in total for 15,486 older adults (mean age 63, range 50-84 years) participating in two cohort studies: ESTHER (Germany) and TROMSØ (Norway); and a subset of previously published nested-case control data from a another cohort study: EPIC-Elderly (Greece, Denmark, Netherlands, Spain and Sweden) from the CHANCES consortium on health and aging. Cox proportional hazards or logistic regression were used to derive multivariable adjusted hazard and odds ratios, respectively, and their 95% confidence intervals across 25(OH)D categories. Meta-analyses with random effects models were used to pool study-specific risk estimates. Overall, lower 25(OH)D concentrations were not significantly associated with increased incidence of most of the cancers assessed. However, there was some evidence of increased breast cancer and decreased lymphoma risk with higher 25(OH)D concentrations. Our meta-analyses with individual participant data from three large European population-based cohort studies provide at best limited support for the hypothesis that vitamin D may have a major role in cancer development and prevention among European older adults.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Hip fracture risk and cadmium in erythrocytes: a nested case-control study with prospectively collected samples.

Several studies have investigated the relation between bone mass density and cadmium exposure, but only few studies have been performed on fractures and biomarkers of cadmium. This study analyzed the association between hip fracture risk and cadmium in erythrocytes (Ery-Cd). Prospective samples from the Northern Sweden Health and Disease Study's biobank were used for 109 individuals who later in life had sustained a low-trauma hip fracture, matched with two controls of the same age and gender. The mean concentration of Ery-Cd (±SD) in case samples was 1.3 ± 1.4 versus 0.9 ± 1.0 µg/L in controls. The odds ratio (OR) was 1.63 [95% confidence interval (CI) 1.10-2.42] for suffering a hip fracture for each microgram per liter increase in Ery-Cd. However, when taking smoking into consideration (never, former, or current), neither Ery-Cd nor smoking showed a statistically significant increase in fracture risk. Using multiple conditional logistic regression with BMI, height, and smoking, the estimated OR for a 1-µg/L increase in Ery-Cd was 1.52 (95% CI 0.77-2.97). Subgroup analysis showed an increased fracture risk among women (OR = 1.94, 95% CI 1.18-3.20, for a 1 µg/L increase), which also remained in the multiple analysis (OR = 3.33, 95% CI 1.29-8.56). This study shows that fracture risk is associated with Ery-Cd. It is, however, not possible to draw firm conclusions on whether cadmium is the causal factor or whether other smoking-related factors cause this association. Subgroup analysis shows that cadmium is a risk factor for hip fracture among women.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm.

BACKGROUND: Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. AIM: To identify genetic variants associated with forearm BMD and forearm fractures. METHODS: BMD at distal radius, measured by dual-energy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a meta-analysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. RESULTS: We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10(-8)) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10(-9)). The gene-based association test suggested an association between MEF2C and forearm BMD (p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. CONCLUSIONS: These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

Haplotype analysis identifies functional elements in monoclonal gammopathy of unknown significance.

Genome-wide association studies (GWASs) based on common single nucleotide polymorphisms (SNPs) have identified several loci associated with the risk of monoclonal gammopathy of unknown significance (MGUS), a precursor condition for multiple myeloma (MM). We hypothesized that analyzing haplotypes might be more useful than analyzing individual SNPs, as it could identify functional chromosomal units that collectively contribute to MGUS risk. To test this hypothesis, we used data from our previous GWAS on 992 MGUS cases and 2910 controls from three European populations. We identified 23 haplotypes that were associated with the risk of MGUS at the genome-wide significance level (p < 5 × 10-8) and showed consistent results among all three populations. In 10 genomic regions, strong promoter, enhancer and regulatory element-related histone marks and their connections to target genes as well as genome segmentation data supported the importance of these regions in MGUS susceptibility. Several associated haplotypes affected pathways important for MM cell survival such as ubiquitin-proteasome system (RNF186, OTUD3), PI3K/AKT/mTOR (HINT3), innate immunity (SEC14L1, ZBP1), cell death regulation (BID) and NOTCH signaling (RBPJ). These pathways are important current therapeutic targets for MM, which may highlight the advantage of the haplotype approach homing to functional units.

Deciphering the genetics and mechanisms of predisposition to multiple myeloma.

Multiple myeloma (MM) is an incurable malignancy of plasma cells. Epidemiological studies indicate a substantial heritable component, but the underlying mechanisms remain unclear. Here, in a genome-wide association study totaling 10,906 cases and 366,221 controls, we identify 35 MM risk loci, 12 of which are novel. Through functional fine-mapping and Mendelian randomization, we uncover two causal mechanisms for inherited MM risk: longer telomeres; and elevated levels of B-cell maturation antigen (BCMA) and interleukin-5 receptor alpha (IL5RA) in plasma. The largest increase in BCMA and IL5RA levels is mediated by the risk variant rs34562254-A at TNFRSF13B. While individuals with loss-of-function variants in TNFRSF13B develop B-cell immunodeficiency, rs34562254-A exerts a gain-of-function effect, increasing MM risk through amplified B-cell responses. Our results represent an analysis of genetic MM predisposition, highlighting causal mechanisms contributing to MM development.

Low Circulating Valine Associate With High Risk of Hip Fractures.

CONTEXT: Hip fractures constitute a major health concern. An adequate supply of amino acids is crucial to ensure optimal acquisition and remodeling of bone. Circulating amino acid levels have been proposed as markers of bone mineral density, but data on their ability to predict incident fractures are scarce. OBJECTIVES: To investigate the associations between circulating amino acids and incident fractures. METHODS: We used UK Biobank (n = 111 257; 901 hip fracture cases) as a discovery cohort and the Umeå Fracture and Osteoporosis (UFO) hip fracture study (hip fracture cases n = 2225; controls n = 2225) for replication. Associations with bone microstructure parameters were tested in a subsample of Osteoporotic Fractures in Men Sweden (n = 449). RESULTS: Circulating valine was robustly associated with hip fractures in the UK Biobank (HR per SD increase 0.79, 95% CI 0.73-0.84), and this finding was replicated in the UFO study (combined meta-analysis including 3126 incident hip fracture cases, odds ratio per SD increase 0.84, 95% CI 0.80-0.88). Detailed bone microstructure analyses showed that high circulating valine was associated with high cortical bone area and trabecular thickness. CONCLUSION: Low circulating valine is a robust predictor of incident hip fractures. We propose that circulating valine may add information for hip fracture prediction. Future studies are warranted to determine whether low valine is causally associated with hip fractures.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases.

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

An atlas of genetic determinants of forearm fracture.

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

Assessment of the genetic and clinical determinants of hip fracture risk: Genome-wide association and Mendelian randomization study.

Hip fracture is the clinically most important fracture, but the genetic architecture of hip fracture is unclear. Here, we perform a large-scale hip fracture genome-wide association study meta-analysis and Mendelian randomization study using five cohorts from European biobanks. The results show that five genetic signals associate with hip fractures. Among these, one signal associates with falls, but not with bone mineral density (BMD), while four signals are in loci known to be involved in bone biology. Mendelian randomization analyses demonstrate a strong causal effect of decreased femoral neck BMD and moderate causal effects of Alzheimer's disease and having ever smoked regularly on risk of hip fractures. The substantial causal effect of decreased femoral neck BMD on hip fractures in both young and old subjects and in both men and women supports the use of change in femoral neck BMD as a surrogate outcome for hip fractures in clinical trials.

The influence of adult hip shape genetic variants on adolescent hip shape: Findings from a population-based DXA study.

OBJECTIVE: Hip shape is a well-recognized risk factor for hip osteoarthritis (OA) and hip fracture. We aimed to investigate whether the genetic variants known to be associated with adult hip shape were also associated with adolescent hip shape. METHODS: Hip DXA scans, obtained in offspring from the Avon Longitudinal Study of Parents and Children (ALSPAC) at two time points (mean ages 13.8 and 17.8 years), were used to quantify hip morphology using a 53-point Statistical Shape Model (SSM). Principal component analysis was used to generate hip shape modes (HSMs). Genetic variants which had previously shown genome-wide significant association with specific HSMs in adults were tested for association with the same HSMs in adolescents (at each timepoint separately) using SNPTEST v2. RESULTS: Complete genotypic and phenotypic data were available for 3550 and 3175 individuals at 14 and 18 years, respectively. The strongest evidence for association with adolescent hip shape was for a variant located near SOX9 (rs2158915) with consistent effects across both time points for HSM1 (age 14: beta -0.05, p = 9.9 × 10-8; age 18: -0.05, p = 3.3 × 10-6) and HSM5 (age 14: beta -0.07, p = 1.6 × 10-4; age 18: -0.1, p = 2.7 × 10-6). There was also strong evidence of association between rs10743612 (near PTHLH) and HSM1 (age 14: 0.05, p = 1.1 × 10-5; age 18: 0.04, p = 0.003) and between rs6537291 (near HHIP) and HSM2 (age 14: -0.06, p = 0.001; age 18: -0.07, p = 0.001) across both time points. The genes with the strongest associations with hip shape in adolescents, (SOX9, PTHLH and HHIP) are known to be involved in endochondral bone formation. HSM1 indicates narrower aspect ratio of the upper femur, whereas both HSM2 and HSM5 reflect variation in the femoral head size and femoral neck width, features previously found to be related to the risk of OA in later life. The SOX9 locus has previously been found to associate with increased risk of hip fracture. CONCLUSION: In conclusion, variants implicated in endochondral bone formation appear to consistently influence hip shape between adolescence and adulthood, including those aspects related to risk of hip OA and/or fracture in later life.

BMD-Related Genetic Risk Scores Predict Site-Specific Fractures as Well as Trabecular and Cortical Bone Microstructure.

CONTEXT: It is important to identify patients at highest risk of fractures. OBJECTIVE: To compare the separate and combined performances of bone-related genetic risk scores (GRSs) for prediction of forearm, hip and vertebral fractures separately, as well as of trabecular and cortical bone microstructure parameters separately. DESIGN, SETTING, AND PARTICIPANTS: Using 1103 single nucleotide polymorphisms (SNPs) independently associated with estimated bone mineral density of the heel (eBMD), we developed a weighted GRS for eBMD and determined its contribution to fracture prediction beyond 2 previously developed GRSs for femur neck BMD (49 SNPs) and lumbar spine BMD (48 SNPs). Associations between these GRSs and forearm (ncases = 1020; ncontrols = 2838), hip (ncases = 1123; ncontrols = 2630) and vertebral (ncases = 288; ncontrols = 1187) fractures were evaluated in 3 Swedish cohorts. Associations between the GRSs and trabecular and cortical bone microstructure parameters (n = 426) were evaluated in the MrOS Sweden cohort. RESULTS: We found that eBMDGRS was the only significant independent predictor of forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS were significant independent predictors of hip fractures. The eBMDGRS was the major GRS contributing to prediction of trabecular bone microstructure parameters while both FN-BMDGRS and eBMDGRS contributed information for prediction of cortical bone microstructure parameters. CONCLUSIONS: The eBMDGRS independently predicts forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS contribute independent information for prediction of hip fractures. We propose that eBMDGRS captures unique information about trabecular bone microstructure useful for prediction of forearm and vertebral fractures. These findings may facilitate personalized medicine to predict site-specific fractures as well as cortical and trabecular bone microstructure separately.

[Preparing for the licence to prescribe in medical school - a questionnaire study on medical students professional confidence in the art of prescribing]. / Läkemedelsarbete behöver vara integrerat i klinisk utbildning - Att göra många läkemedelsgenomgångar ökar trygghet och reflektion över patientens behandling, visar enkätstudie.

A prerequisite for rational use of medicines is adequate prescribing skills; drug treatment is a complex task requiring diagnostic competence combined with pharmacologic knowledge and patient communication skills. Acquiring professional confidence in the art of prescribing is essential during medical training. The results of this questionnaire study, conducted in four medical schools in Sweden after the course in internal medicine (252 respondents; response rate: 74%; median age: 24 years, 61% female), show that 45% and 62% were confident in performing medication reviews and writing medication summary reports, respectively, i.e. the basics of prescribing. The confidence increased by the number of reviews and reports performed, i.e. the extent of practice (correlation coefficients: 0.41 and 0.38, respectively, both p<0.0001), as did the extent of the students' reflection on important aspects of drug treatment such as adherence, adverse reactions, renal function, dosing, and drug interactions. In multivariate regression analyses, major predictors for confidence in performing medication reviews were extent of practice and extent of clinical supervision. The results suggest that these factors are keys to acquiring professional confidence in the art of prescribing.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies.

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2 > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.