Bone mineral density (BMD) is the most important predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and East Asian ancestry. We tested the top-associated BMD markers for replication in 50,933 independent subjects and for risk of low-trauma fracture in 31,016 cases and 102,444 controls. We identified 56 loci (32 novel)associated with BMD atgenome-wide significant level (P<5×10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal-stem-cell differentiation, endochondral ossification and the Wnt signalling pathways. However, we also discovered loci containing genes not known to play a role in bone biology. Fourteen BMD loci were also associated with fracture risk (P<5×10−4, Bonferroni corrected), of which six reached P<5×10−8 including: 18p11.21 (C18orf19), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

SUMMARY 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 largely unknown. Bone mineral density (BMD) is highly heritable, is a major predictor of osteoporotic fractures and has been previously associated with common genetic variants1–8, and rare, population-specific, coding variants9. 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), whole-exome sequencing (n= 3,549), deep imputation of genotyped samples using a combined UK10K/1000Genomes 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 4-fold larger than the mean of previously reported common variants for lumbar spine BMD8 (rs11692564[T], MAF = 1.7%, replication effect size = +0.20 standard deviations [SD], Pmeta = 2×10−14), which was also associated with a decreased risk of fracture (OR = 0.85; P = 2×10−11; ncases = 98,742 and ncontrols = 409,511). Using an En1Cre/flox mouse model, we observed that conditional loss of En1 results in low bone mass, likely as a consequence of high bone turn-over. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817[T], MAF = 1.1%, replication effect size = +0.39 SD, 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.

Objective. To determine genome-wide methylation profiles of bone from patients with hip osteoarthritis (OA) and those with osteoporotic (OP) hip fractures.Methods. Trabecular bone pieces were obtained from the central part of the femoral head of 27 patients with hip fractures and 26 patients with hip OA. DNA was isolated, and methylation was explored with Illumina methylation arrays. RNA was extracted, pooled, and deep-sequenced to obtain the whole transcriptome. Differentially methylated regions were identified, and connections between genes with differentially methylated regions were explored by pathway and text-mining analyses.Results. After quality control, methylation of 23,367 CpG sites (13,463 genes) was analyzed. There was a genome-wide inverse relationship between methylation and gene expression in both patient groups. Comparison of OP and OA bones revealed 241 CpG sites, located in 228 genes, with significant differences in methylation (false discovery rate <0.05). Of them, 217 were less methylated in OP than in OA. The absolute methylation differences were >5% in 128 CpG sites and >10% in 45 CpG sites. The differentially methylated genes were enriched for association with bone traits in the genome-wide association study catalog. Pathway analysis and text-mining analysis with Gene Relationships Across Implicated Loci software revealed enrichment in genes participating in glycoprotein metabolism or cell differentiation, and particularly in the homeobox superfamily of transcription factors.Conclusion. Genome-wide methylation profiling of bone samples revealed differentially methylated regions in OP and OA. These regions were enriched in genes associated with cell differentiation and skeletal embryogenesis, such as those in the homeobox superfamily, suggesting the existence of a developmental component in the predisposition to these disorders.Bone increases in size during the growth period by a modeling process driven by the formation of new bone. Thereafter, it is constantly remodeled by the concerted action of bone-resorbing osteoclasts and bone-forming osteoblasts, originating from hematopoietic and mesenchymal precursors, respectively. When remodeling is to start at a certain site, osteoclast precursors are attracted and committed to differentiate into mature osteoclasts. When the resorption phase ends, surrounding osteoblast precursors proliferate and differentiate into mature osteoblasts that synthesize bone matrix that eventually mineralizes and replaces the old bone resorbed by osteoclasts. Thus, bone remodeling requires the cyclic and sequential proliferation and differentiation of osteoclast and osteoblast precursors. Any disturbance of this process will result in abnormal bone mass. This is the case in osteoporosis (OP), which is characterized by a decrease in bone mass, due to high bone resorption and/or low bone formation, and consequently propensity to fracture. In contrast to OP, in osteoarthritis (OA) bone formation is increased in the

A collaborative work was carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG) to estimate Y-STR mutation rates. Seventeen Y chromosome STR loci (DYS19, DYS385, DYS389I and II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS460, DYS461, DYS635 [GATA C4], GATA H4, and GATA A10) were analyzed in a sample of 3,026 father/son pairs. Among 27,029 allele transfers, 54 mutations were observed, with an overall mutation rate across the 17 loci of 1.998 x 10(-3) (95% CI, 1.501 x 10(-3) to 2.606 x 10(-3)). With just one exception, all of the mutations were single-step, and they were observed only once per gametogenesis. Repeat gains were more frequent than losses, longer alleles were found to be more mutable, and the mutation rate seemed to increase with the father's age. Hum Mutat 26(6), 520-528, 2005. (c) 2005 Wiley-Liss, Inc.

In a collaborative work carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG), a polymerase chain reaction multiplex was optimized in order to type ten X-chromosome short tandem repeats (STRs) in a single reaction, including: DXS8378, DXS9902, DXS7132, DXS9898, DXS6809, DXS6789, DXS7133, GATA172D05, GATA31E08, and DXS7423. Using this X-decaplex, each 17 of the participating laboratories typed a population sample of approximately 200 unrelated individuals (100 males and 100 females). In this work, we report the allele frequencies for the ten X-STRs in 15 samples from Argentina (Buenos Aires, Córdoba, Río Negro, Entre Ríos, and Misiones), Brazil (São Paulo, Rio de Janeiro, Paraná, and Mato Grosso do Sul), Colombia (Antioquia), Costa Rica, Portugal (Northern and Central regions), and Spain (Galicia and Cantabria). Gene diversities were calculated for the ten markers in each population and all values were above 56%. The average diversity per locus varied between 66%, for DXS7133, and 82%, for DXS6809. For this set of STRs, a high discrimination power was obtained in all populations, both in males (> or =1 in 5 x 10(5)) and females (> or =1 in 3 x 10(9)), as well as high mean exclusion chance in father/daughter duos (> or =99.953%) and in father/mother/daughter trios (> or =99.999%). Genetic distance analysis showed no significant differences between northern and central Portugal or between the two Spanish samples from Galicia and Cantabria. Inside Brazil, significant differences were found between Rio de Janeiro and the other three populations, as well as between São Paulo and Paraná. For the five Argentinean samples, significant distances were only observed when comparing Misiones with Entre Ríos and with Río Negro, the only two samples that do not differ significantly from Costa Rica. Antioquia differed from all other samples, except the one from Río Negro.