Genome-wide association meta-analyses identified 1q43 and 2q32.2 for hip Ward's triangle areal bone mineral density

Genome-wide association study of osteoporosis identifies genetic risk and interactions with Dietary Approaches to Stop Hypertension diet and sugar-sweetened beverages in a Hispanic cohort of older adults

Osteoporosis (OP) and low bone mass can be debilitating and costly conditions if not acted on quickly. This disease is also difficult to diagnose as the symptoms develop unnoticed until fracture occurs. Therefore, gaining understanding of the genetic risk associated with these conditions could be beneficial for health-care professionals in early detection and prevention. The Boston Puerto Rican Osteoporosis (BPROS) study, an ancillary study to the Boston Puerto Rican Health Study (BPRHS), collected information regarding bone and bone health. All bone measurements were taken during regular BPROS visits using dual-energy X-ray absorptiometry. The OP was defined as T-score ≤ -2.5 (≥2.5 SDs below peak bone mass). Dietary variables were collected at the second wave of the BPRHS via a food frequency questionnaire. We conducted genome-wide associations with bone outcomes, including BMD and OP for 978 participants. We also examined the interactions with dietary quality on the relationships between genotype and bone outcomes. We further tested if candidate genetic variants described in previous GWAS on OP and BMD contribute to OP risk in this population. Four variants were associated with OP: rs114829316 (IQ motif containing J gene), rs76603051, rs12214684 (melanin-concentrating hormone receptor 2 gene), and rs77303493 (Ras and Rab interactor 2 gene), and 2 variants were associated with BMD of lumbar spine (rs11855618, cingulin-like 1 gene) and hip (rs73480593, NTRK2), reaching the genome-wide significance threshold of P ≤ 5E-08. In a gene-diet interaction analysis, we found that 1 SNP showed a significant interaction with the overall Dietary Approaches to Stop Hypertension (DASH) score, and 7 SNPs with sugar-sweetened beverages (SSBs), a major contributor to the DASH score. This study identifies new genetic markers related to OP and BMD in older Hispanic adults. Additionally, we uncovered unique genetic markers that interact with dietary quality, specifically SSBs, in relation to bone health. These findings may be useful to guide early detection and preventative care.

The roles and regulatory mechanisms of TGF-β and BMP signaling in bone and cartilage development, homeostasis and disease

Transforming growth factor-βs (TGF-βs) and bone morphometric proteins (BMPs) belong to the TGF-β superfamily and perform essential functions during osteoblast and chondrocyte lineage commitment and differentiation, skeletal development, and homeostasis. TGF-βs and BMPs transduce signals through SMAD-dependent and -independent pathways; specifically, they recruit different receptor heterotetramers and R-Smad complexes, resulting in unique biological readouts. BMPs promote osteogenesis, osteoclastogenesis, and chondrogenesis at all differentiation stages, while TGF-βs play different roles in a stage-dependent manner. BMPs and TGF-β have opposite functions in articular cartilage homeostasis. Moreover, TGF-β has a specific role in maintaining the osteocyte network. The precise activation of BMP and TGF-β signaling requires regulatory machinery at multiple levels, including latency control in the matrix, extracellular antagonists, ubiquitination and phosphorylation in the cytoplasm, nucleus-cytoplasm transportation, and transcriptional co-regulation in the nuclei. This review weaves the background information with the latest advances in the signaling facilitated by TGF-βs and BMPs, and the advanced understanding of their diverse physiological functions and regulations. This review also summarizes the human diseases and mouse models associated with disordered TGF-β and BMP signaling. A more precise understanding of the BMP and TGF-β signaling could facilitate the development of bona fide clinical applications in treating bone and cartilage disorders.

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Uniparental isodisomy of chromosome 2 causing MRPL44-related multisystem mitochondrial disease

Mutations in nuclear-encoded protein subunits of the mitochondrial ribosome are an increasingly recognised cause of oxidative phosphorylation system (OXPHOS) disorders. Among them, mutations in the MRPL44 gene, encoding a structural protein of the large subunit of the mitochondrial ribosome, have been identified in four patients with OXPHOS defects and early-onset hypertrophic cardiomyopathy with or without additional clinical features. A 23-year-old individual with cardiac and skeletal myopathy, neurological involvement, and combined deficiency of OXPHOS complexes in skeletal muscle was clinically and genetically investigated. Analysis of whole-exome sequencing data revealed a homozygous mutation in MRPL44 (c.467 T > G), which was not present in the biological father, and a region of homozygosity involving most of chromosome 2, raising the possibility of uniparental disomy. Short-tandem repeat and genome-wide SNP microarray analyses of the family trio confirmed complete maternal uniparental isodisomy of chromosome 2. Mitochondrial ribosome assembly and mitochondrial translation were assessed in patient derived-fibroblasts. These studies confirmed that c.467 T > G affects the stability or assembly of the large subunit of the mitochondrial ribosome, leading to impaired mitochondrial protein synthesis and decreased levels of multiple OXPHOS components. This study provides evidence of complete maternal uniparental isodisomy of chromosome 2 in a patient with MRPL44-related disease, and confirms that MRLP44 mutations cause a mitochondrial translation defect that may present as a multisystem disorder with neurological involvement.

Assessing causal relationship from gut microbiota to heel bone mineral density

Recent studies have demonstrated the important role played by gut microbiota in regulating bone development, but the evidence of such causal relationship is still sparse in human population. The aim of this study is to assess the causal relationship from gut microbiota to bone development and to identify specific causal bacteria taxa via a Mendelian randomization (MR) approach. A genome-wide association study (GWAS) summary statistic based two-sample MR analysis was performed. Summary statistics of microbiome GWAS (MGWAS) in 1126 twin pairs of the TwinsUK study was used as discovery sample, and the MGWAS in 984 Dutch participants from the LifeLines-DEEP cohort was used as replication sample. Estimated heel bone mineral density (eBMD) GWAS in 426,824 participants from the UK biobank (UKB) cohort was used as outcome. Bacteria were grouped into taxa features at both order and family levels. In the discovery sample, a total of 25 bacteria features including 9 orders and 16 families were analyzed. Fourteen features (5 orders + 9 families) were nominally significant, including 5 orders (Bacteroidales, Clostridiales, Lactobacillales, Pasteurellales and Verrucomicrobiales) and 9 families (Bacteroidaceae, Clostridiaceae, Lachnospiraceae, Mogibacteriaceae, Pasteurellaceae, Porphyromonadaceae, Streptococcaceae, Verrucomicrobiaceae and Veillonellaceae). One order Clostridiales and its child taxon, family Lachnospiraceae, were successfully replicated in the replication sample (Clostridiales P_discovery = 3.32 × 10^-3P_replication = 7.29 × 10^-3; Lachnospiraceae P_discovery = 0.03 P_replication = 7.29 × 10^-3). Our findings provided evidence of causal relationship from microbiota to bone development, as well as identified specific bacteria taxa that regulated bone mass variation, thus providing new insights into the microbiota mediated bone development mechanism.

Association of Gremlin-2 gene polymorphisms with osteoporosis risk in Chinese postmenopausal women

The Gremlin-2 (GREM2) plays crucial roles in modulating bone homeostasis through the bone morphogenetic protein-2 pathway. However, GREM2 gene variants in osteoporosis were less frequent in a Chinese population. Therefore, the present study recruited 310 patients with osteoporosis and 339 healthy postmenopausal women to assess the correlation of GREM2 gene polymorphisms with the risk of osteoporosis. Polymerase chain reaction (PCR) and Sanger sequencing were utilized to genotype samples. The results showed that GREM2 gene rs4454537, not rs11588607, polymorphism was significantly associated with an increased risk of osteoporosis in postmenopausal women. Moreover, stratified analyses indicated a significant association between rs4454537 polymorphisms and body mass index of <25 kg/m². Additionally, the association between GREM2 rs4454537 polymorphism and clinical characteristics was assessed, which showed that this locus decreased the bone mineral density (BMD) in postmenopausal osteoporotic individuals. Furthermore, individuals with CC genotype appeared to have a higher GREM2 expression compared with those bearing the TT genotype of rs4454537 polymorphism. However, the genotype distribution of rs4454537 polymorphism showed no statistical difference between osteoporotic patients as a function of fracture status. In summary, GREM2 rs4454537 polymorphism decreases BMD and increases osteoporotic risk in postmenopausal women.

Twelve New Genomic Loci Associated With Bone Mineral Density

Aiming to identify more genomic loci associated with bone mineral density (BMD), we conducted a joint association analysis of 2 genome-wide association study (GWAS) by the integrative association method multi-trait analysis of GWAS (MTAG). The first one is the single GWAS of estimated heel BMD (eBMD) in the UK biobank (UKB) cohort (N = 426,824), and the second one is the GWAS meta-analysis of total body BMD (TB-BMD) in 66,628 participants from 30 studies. Approximate conditional association analysis was performed in the identified novel loci to identify secondary association signal. Statistical fine-mapping was conducted to prioritize plausible credible risk variants (CRVs). Candidate genes were prioritized based on the analyses of cis- expression quantitative trait locus (cis-eQTL) and cis-protein QTL (cis-pQTL) information as well as the functional category of the SNP. By integrating the information carried in over 490,000 participants, this largest joint analysis of BMD GWAS identified 12 novel genomic loci at the genome-wide significance level (GWS, p = 5.0 × 10^-8), nine of which were for eBMD and four were for TB-BMD, explaining an additional 0.11 and 0.23% heritability for the two traits, respectively. These loci include 1p33 (lead SNP rs10493130, p_eBMD = 3.19 × 10^-8), 5q13.2 (rs4703589, p_eBMD = 4.78 × 10^-8), 5q31.3 (rs9324887, p_TB-BMD = 1.36 × 10^-9), 6p21.32 (rs6905837, p_eBMD = 3.32 × 10^-8), 6q14.1 (rs10806234, p_eBMD = 2.63 × 10^-8), 7q21.11 (rs10806234, p_TB-BMD = 3.37 × 10^-8), 8q24.12 (rs11995866, p_eBMD = 6.72 × 10^-9), 12p13.31 (rs1639122, p_eBMD = 4.43 × 10^-8), 12p12.1 (rs58489179, p_eBMD = 4.74 × 10^-8), 12q24.23 (rs75499226, p_eBMD = 1.44 × 10^-8), 19q13.31 (rs7255083, p_TB-BMD = 2.18 × 10^-8) and 22q11.23 (rs13056137, p_TB-BMD = 2.54 × 10^-8). All lead SNPs in these 12 loci are nominally significant in both original studies as well as consistent in effect direction between them, providing solid evidence of replication. Approximate conditional analysis identified one secondary signal in 5q13.2 (rs11738874, p_conditional = 5.06 × 10^-9). Statistical fine-mapping analysis prioritized 269 CRVs. A total of 65 candidate genes were prioritized, including those with known biological function to bone development (such as FGF1, COL11A2 and DEPTOR). Our findings provide novel insights into a better understanding of the genetic mechanism underlying bone development as well as candidate genes for future functional investigation.