Identification of 613 new loci associated with heel bone mineral density and a polygenic risk score for bone mineral density, osteoporosis and fracture

Bridging Genomic Research Disparities in Osteoporosis GWAS: Insights for Diverse Populations

PURPOSE OF REVIEW

Genome-wide association studies (GWAS) have significantly advanced osteoporosis research by identifying genetic loci associated with bone mineral density (BMD) and fracture risk. However, disparities persist due to the underrepresentation of non-European populations, limiting the applicability of polygenic risk scores (PRS). This review examines recent advancements in osteoporosis genetics, highlights existing disparities, and explores strategies for more inclusive research.

RECENT FINDINGS

European-focused GWAS have identified key loci for osteoporosis, including WNT signaling (SOST, LRP5) and RUNX2 transcriptional regulation. However, fewer than 40% of these variants can be replicated in Asian and African populations. Emerging studies in non-European groups reveal population-specific loci, sex-specific associations, and gene-environment interactions. Advances in machine learning (ML)-assisted GWAS and multi-omics integration are improving genetic discovery. Expanding GWAS in diverse populations, integrating multi-omics data, refining ML-based risk models, and standardizing biobank data are essential for equitable osteoporosis research. Future efforts must prioritize clinical translation to enhance personalized osteoporosis prevention and treatment.

Mapping the role of macro and micronutrients in bone mineral density: a comprehensive Mendelian randomization study

BACKGROUND

Macro and micronutrients may play an important role in osteoporosis development; however, observational studies have yielded inconsistent results. Clarifying these associations is vital to the development of nutritional recommendations aimed at preventing osteoporosis.

METHODS

Utilizing the largest available genome-wide association study (GWAS) summary statistics to date, we performed a two-sample Mendelian randomization (MR) analysis to investigate the causal effects of energy-adjusted macronutrient intake (fat, protein, carbohydrate, and sugar) and circulating levels of 20 micronutrients (ten each for vitamins and minerals) on heel estimated bone mineral density (eBMD), a promising marker for osteoporosis risk and fracture susceptibility. Sensitivity, sex-specific, and one-sample MR analyses were applied to further validate and annotate the results.

RESULTS

Among all nutrients, four genetically predicted circulating levels of micronutrients were suggestively associated with eBMD (vitamin A: β IVW = - 0.054, PIVW = 3.70 × 10-2; vitamin B12: β IVW = - 0.020, PIVW = 3.71 × 10-6; vitamin E: β IVW = 0.277, PIVW = 3.22 × 10-2; selenium: β IVW = 0.023, PIVW = 5.37 × 10-3; all Pintercept > 0.05). All these results were also considered robust, as sensitivity analyses yielded directionally consistent results. However, only the causal effects of vitamin B12 and selenium on eBMD remained significant after Bonferroni correction and were not confounded by obesity, smoking status, or alcohol consumption. Sex-specific analysis revealed a male-specific causal association between vitamin E and eBMD ( β IVW = 0.275, PIVW = 9.81 × 10-14). Additionally, using individuallevel data from the UK Biobank cohort, one-sample MR analysis found no causal relationships between diet-derived nutrient intake and eBMD in the overall population, as well as in the females or the males.

CONCLUSIONS

Our findings suggest that appropriate levels of plasma vitamin B12 and adequate levels of serum selenium are crucial for delaying bone loss and promoting bone health, emphasizing the need for nutritional recommendations to maintain optimal levels of these nutrients to promote eBMD and prevent osteoporosis.

Regulation of WNT16 in bone may involve upstream enhancers within CPED1

WNT16 stands up as an essential gene for bone homeostasis. Here, we present new evidence of the functional role of a particular region within WNT16. Performing 4 C chromatin conformation analysis in three osteoblast-related cells (the human fetal osteoblast hFOB 1.19 cell line, Saos 2 osteosarcoma cell line and mesenchymal Stem Cells -MSC-), we identify physical interactions between the proximal part of WNT16 intron 2, shown here to be an active promoter in Saos 2 osteosarcoma cells, and several putative regulatory regions within CPED1. Analysis of previously published RNA-seq data from hFOB cells disclosed low expression of a region located downstream of this promoter. Our results suggest a novel regulatory mechanism of WNT16 in bone, mediated by physical interaction with various enhancer regions within CPED1.

Unveiling two distinct osteolineage cell populations linked to age-related osteoporosis in adult mice through integrative single-cell analyses

The bone marrow microenvironment contains heterogeneous stromal cells, which are critical for bone remodeling and provide essential supportive roles for hematopoietic functions. Although the diversity of PDGFRα+β+ mesenchymal stromal/stem cells (MSCs) get consensus, the osteo-lineage cells (OLCs) that constitute the developmental trajectory of osteoblasts are largely remain unclear. Here, we construct a comprehensive atlas of stromal cell via performing integrative single cell analyses for 77 samples from 14 datasets. Besides previously defined Lepr+ BM-MSCs derived OLC1, we present a novel OLC2 with unique molecular signatures and differentiation pathway. Both OLC1 and OLC2 show significant polygenic association with bone mineral density (BMD), while extracellular matrix (ECM) proteins specifically expressed by OLC2 are particularly decreased in bone tissue of aged mice. Notably, OLC1 and OLC2 are consistently reduced in aged mice, which might be induced by the decreased expression of several lineage drivers. Collectively, our study presents a thorough prospect of OLCs in the bone marrow microenvironment and highlights their important roles in age-related osteoporosis.

Early-life tobacco smoke elevating later-life osteoporosis risk: Mediated by telomere length and interplayed with genetic predisposition

INTRODUCTION

The growing prevalence of osteoporosis (OP) in an aging global population presents a significant public health concern. Tobacco smoke negatively affects bone turnover, leading to reduced bone mass and heightened OP and fracture risk. However, the impact of early-life tobacco smoke exposure on later-life OP risk remains unclear.

OBJECTIVES

This study was to explore the effects of early-life tobacco smoke exposure on incident OP risk in later life. The mediating role of telomere length (TL) and the interaction with genetic predisposition were also studied.

METHODS

Data on in utero tobacco smoke exposure (IUTSE) status and age of tobacco use initiation from the UK Biobank were used to estimate early-life tobacco smoke exposure. Incident OP cases were identified according to health-related records. Linear, Cox, and Laplace regression models were mainly used for data analysis.

RESULTS

Individuals with IUTSE showed a higher OP risk [hazard ratio (HR): 1.06, 95 % confidence interval (CI): 1.01, 1.11] and experienced earlier OP onset by 0.30 years [50th percentile difference = -0.30, 95 % CI: -0.51, -0.09] compared to those without. Participants initiating tobacco smoke in childhood, adolescence, and adulthood had 1.41 times (95 % CI: 1.23, 1.61), 1.17 times (95 % CI:1.10, 1.24), and 1.14 times (95 % CI: 1.07, 1.20) the risk of OP, respectively, compared to never smokers. They also experienced earlier OP onset by 2.16, 0.95, and 0.71 years, sequentially. The TL significantly mediated the early-life tobacco exposure and OP association. Significant joint and interactive effects were detected between early-life tobacco smoke exposure and genetic elements.

CONCLUSIONS

Our findings implicate that early-life tobacco smoke exposure elevates the later-life OP risk, mediated by telomere length and interplayed with genetic predisposition. These findings highlight the importance of early-life intervention against tobacco smoke exposure and ageing status for precise OP prevention, especially in individuals with a high genetic risk.

GWAS-Informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblasts (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from stratified LD-score regression involving 98 cell types grouped into 15 tissues. 23 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed unexpectedly that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues. Extending our CRISPRi screening approach to these tissues could play a key role in fully elucidating the etiology of BMD.

The Genetic Variants Influencing Hypertension Prevalence Based on the Risk of Insulin Resistance as Assessed Using the Metabolic Score for Insulin Resistance (METS-IR)

Insulin resistance is a major indicator of cardiovascular diseases, including hypertension. The Metabolic Score for Insulin Resistance (METS-IR) offers a simplified and cost-effective way to evaluate insulin resistance. This study aimed to identify genetic variants associated with the prevalence of hypertension stratified by METS-IR score levels. Data from the Korean Genome and Epidemiology Study (KoGES) were analyzed. The METS-IR was calculated using the following formula: ln [(2 × fasting blood glucose (FBG) + triglycerides (TG)) × body mass index (BMI)]/ ln [high-density lipoprotein cholesterol (HDL-C)]. The participants were divided into tertiles 1 (T1) and 3 (T3) based on their METS-IR scores. Genome-wide association studies (GWAS) were performed for hypertensive cases and non-hypertensive controls within these tertile groups using logistic regression adjusted for age, sex, and lifestyle factors. Among the METS-IR tertile groups, 3517 of the 19,774 participants (17.8%) at T1 had hypertension, whereas 8653 of the 20,374 participants (42.5%) at T3 had hypertension. A total of 113 single-nucleotide polymorphisms (SNPs) reached the GWAS significance threshold (p < 5 × 10-8) in at least one tertile group, mapping to six distinct genetic loci. Notably, four loci, rs11899121 (chr2p24), rs7556898 (chr2q24.3), rs17249754 (ATP2B1), and rs1980854 (chr20p12.2), were significantly associated with hypertension in the high-METS-score group (T3). rs10857147 (FGF5) was significant in both the T1 and T3 groups, whereas rs671 (ALDH2) was significant only in the T1 group. The GWASs identified six genetic loci significantly associated with hypertension, with distinct patterns across METS-IR tertiles, highlighting the role of metabolic context in genetic susceptibility. These findings underscore critical genetic factors influencing hypertension prevalence and provide insights into the metabolic-genetic interplay underlying this condition.

Mapping rare protein-coding variants on multi-organ imaging traits

Human organ structure and function are important endophenotypes for clinical outcomes. Genome-wide association studies (GWAS) have identified numerous common variants associated with phenotypes derived from magnetic resonance imaging (MRI) of the brain and body. However, the role of rare protein-coding variations affecting organ size and function is largely unknown. Here we present an exome-wide association study that evaluates 596 multi-organ MRI traits across over 50,000 individuals from the UK Biobank. We identified 107 variant-level associations and 224 gene-based burden associations (67 unique gene-trait pairs) across all MRI modalities, including PTEN with total brain volume, TTN with regional peak circumferential strain in the heart left ventricle, and TNFRSF13B with spleen volume. The singleton burden model and AlphaMissense annotations contributed 8 unique gene-trait pairs including the association between an approved drug target gene of KCNA5 and brain functional activity. The identified rare coding signals elucidate some shared genetic regulation across organs, prioritize previously identified GWAS loci, and are enriched for drug targets. Overall, we demonstrate how rare variants enhance our understanding of genetic effects on human organ morphology and function and their connections to complex diseases.

A proteomics approach to study mouse long bones: examining baseline differences and mechanical loading-induced bone formation in young-adult and old mice

With aging, bone mass declines and the anabolic effects of skeletal loading diminish. While much research has focused on gene transcription, how bone ages and loses its mechanoresponsiveness at the protein level remains unclear. We developed a novel proteomics approach and performed a paired mass spectrometry and RNA-seq analysis on tibias from young-adult (5-month) and old (22-month) mice. We report the first correlation estimate between the bone proteome and transcriptome (Spearman ρ = 0.40), which is in line with other tissues but indicates that a relatively low amount of variation in protein levels is explained by the variation in transcript levels. Of 71 shared targets that differed with age, eight were associated with bone mineral density in previous GWAS, including understudied targets Asrgl1 and Timp2. We used complementary RNA in situ hybridization to confirm that Asrgl1 and Timp2 had reduced expression in osteoblasts/osteocytes in old bones. We also found evidence for reduced TGF-beta signaling with aging, in particular Tgfb2. Next, we defined proteomic changes following mechanical loading. At the protein level, bone differed more with age than with loading, and aged bone had fewer loading-induced changes. Overall, our findings underscore the need for complementary protein-level assays in skeletal biology research.

Shared genetic basis connects smoking behaviors and bone health: insights from a genome-wide cross-trait analysis

Although the negative association of tobacco smoking with osteoporosis is well-documented, little is known regarding the shared genetic basis underlying these conditions. In this study, we aim to investigate a shared genetic architecture between smoking and heel estimated bone mineral density (eBMD), a reliable proxy for osteoporosis. We conducted a comprehensive genome-wide cross-trait analysis to identify genetic correlation, pleiotropic loci and causal relationship of smoking with eBMD, leveraging summary statistics of the hitherto largest genome-wide association studies conducted in European ancestry for smoking initiation (Nsmoker = 1 175 108, Nnonsmoker = 1 493 921), heaviness (cigarettes per day, N = 618 489), cessation (Ncurrent smoker = 304 244, Nformer smoker = 843 028), and eBMD (N = 426 824). A significant negative global genetic correlation was found for smoking cessation and eBMD (${r}_g$ = -0.051, P = 0.01), while we failed to identify a significant global genetic correlation of smoking initiation or heaviness with eBMD. Partitioning the whole genome into independent blocks, we observed 6 significant shared local signals for smoking and eBMD, with 22q13.1 showing the strongest regional genetic correlation. Such a genetic overlap was further supported by 71 pleiotropic loci identified in the cross-trait meta-analysis. Mendelian randomization identified no causal effect of smoking initiation (beta = -0.003 g/cm2, 95% CI = -0.033 to 0.027) or heaviness (beta = -0.017 g/cm2, 95% CI = -0.072 to 0.038) on eBMD, but a putative causal effect of genetic predisposition to being a current smoker was associated with a lower eBMD compared to former smokers (beta = -0.100 g/cm2, 95% CI = -0.181 to -0.018). Our study demonstrates a pronounced biological pleiotropy as well as a putative causal link between current smoking status and eBMD, providing novel insights into the primary prevention and modifiable intervention of osteoporosis by advocating individuals to avoid, reduce or quit smoking as early as possible.

Identification of a novel QTL for lean meat percentage using imputed genotypes

Lean meat percentage is a critical production trait in pig breeding systems with direct implications for the sustainability of the industry. In this study, we conducted a genome-wide association study for lean meat percentage using a cohort of 850 Duroc × (Landrace × Yorkshire) crossbred pigs and we identified QTL on SSC3 and SSC18. Based on the predicted effect of imputed variants and using the PigGTEx database of molecular QTL, we prioritized candidate genes and SNPs located within the QTL regions, which may be involved in the regulation of porcine leanness. Our results indicate that a nonsense mutation in ZC3HAV1L on SSC18 has a direct effect on lean meat percentage.

Elucidating causal relationships of diet-derived circulating antioxidants and the risk of osteoporosis: A Mendelian randomization study

Background

Osteoporosis (OP) is typically diagnosed by evaluating bone mineral density (BMD), and it frequently results in fractures. Here, we investigated the causal relationships between diet-derived circulating antioxidants and the risk of OP using Mendelian randomization (MR).

Methods

Published studies were used to identify instrumental variables related to absolute levels of circulating antioxidants like lycopene, retinol, ascorbate, and β-carotene, as well as antioxidant metabolites such as ascorbate, retinol, α-tocopherol, and γ-tocopherol. Outcome variables included BMD (in femoral neck, lumbar spine, forearm, heel, total body, total body (age over 60), total body (age 45-60), total body (age 30-45), total body (age 15-30), and total body (age 0-15)), fractures (in arm, spine, leg, heel, and osteoporotic fractures), and OP. Inverse variance weighted or Wald ratio was chosen as the main method for MR analysis based on the number of single nucleotide polymorphisms (SNPs). Furthermore, we performed sensitivity analyses to confirm the reliability of the findings.

Results

We found a causal relationship between absolute retinol levels and heel BMD (p = 7.6E-05). The results of fixed effects IVW showed a protective effect of absolute retinol levels against heel BMD, with per 0.1 ln-transformed retinol being associated with a 28% increase in heel BMD (OR: 1.28, 95% CI: 1.13-1.44). In addition, a sex-specific effect of the absolute circulating retinol levels on the heel BMD has been observed in men. No other significant causal relationship was found.

Conclusion

There is a positive causal relationship between absolute retinol levels and heel BMD. The implications of our results should be taken into account in future studies and in the creation of public health policies and OP prevention tactics.

Genetic Evaluation for Monogenic Disorders of Low Bone Mass and Increased Bone Fragility: What Clinicians Need to Know

PURPOSE OF REVIEW

The purpose of this review is to outline the principles of clinical genetic testing and to provide practical guidance to clinicians in navigating genetic testing for patients with suspected monogenic forms of osteoporosis.

RECENT FINDINGS

Heritability assessments and genome-wide association studies have clearly shown the significant contributions of genetic variations to the pathogenesis of osteoporosis. Currently, over 50 monogenic disorders that present primarily with low bone mass and increased risk of fractures have been described. The widespread availability of clinical genetic testing offers a valuable opportunity to correctly diagnose individuals with monogenic forms of osteoporosis, thus instituting appropriate surveillance and treatment. Clinical genetic testing may identify the appropriate diagnosis in a subset of patients with low bone mass, multiple or unusual fractures, and severe or early-onset osteoporosis, and thus clinicians should be aware of how to incorporate such testing into their clinical practices.

Reconstructing the ancestral gene pool to uncover the origins and genetic links of Hmong-Mien speakers

BACKGROUND

Hmong-Mien (HM) speakers are linguistically related and live primarily in China, but little is known about their ancestral origins or the evolutionary mechanism shaping their genomic diversity. In particular, the lack of whole-genome sequencing data on the Yao population has prevented a full investigation of the origins and evolutionary history of HM speakers. As such, their origins are debatable.

RESULTS

Here, we made a deep sequencing effort of 80 Yao genomes, and our analysis together with 28 East Asian populations and 968 ancient Asian genomes suggested that there is a strong genetic basis for the formation of the HM language family. We estimated that the most recent common ancestor dates to 5800 years ago, while the genetic divergence between the HM and Tai-Kadai speakers was estimated to be 8200 years ago. We proposed that HM speakers originated from the Yangtze River Basin and spread with agricultural civilization. We identified highly differentiated variants between HM and Han Chinese, in particular, a deafness-related missense variant (rs72474224) in the GJB2 gene is in a higher frequency in HM speakers than in others.

CONCLUSIONS

Our results indicated complex gene flow and medically relevant variants involved in the HM speakers' evolution history.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

Sarm1 knockout prevents type 1 diabetic bone disease in females independent of neuropathy

Patients with diabetes have a high risk of developing skeletal diseases accompanied by diabetic peripheral neuropathy (DPN). In this study, we isolated the role of DPN in skeletal disease with global and conditional knockout models of sterile-α and TIR-motif-containing protein-1 (Sarm1). SARM1, an NADase highly expressed in the nervous system, regulates axon degeneration upon a range of insults, including DPN. Global knockout of Sarm1 prevented DPN, but not skeletal disease, in male mice with type 1 diabetes (T1D). Female wild-type mice also developed diabetic bone disease but without DPN. Unexpectedly, global Sarm1 knockout completely protected female mice from T1D-associated bone suppression and skeletal fragility despite comparable muscle atrophy and hyperglycemia. Global Sarm1 knockout rescued bone health through sustained osteoblast function with abrogation of local oxidative stress responses. This was independent of the neural actions of SARM1, as beneficial effects on bone were lost with neural conditional Sarm1 knockout. This study demonstrates that the onset of skeletal disease occurs rapidly in both male and female mice with T1D completely independently of DPN. In addition, this reveals that clinical SARM1 inhibitors, currently being developed for treatment of neuropathy, may also have benefits for diabetic bone through actions outside of the nervous system.

Genetic Variants of Adrenoceptors

Adrenoceptors are class A G-protein-coupled receptors grouped into three families (α1-, α2-, and β-adrenoceptors), each one including three members. All nine corresponding adrenoceptor genes display genetic variation in their coding and adjacent non-coding genomic region. Coding variants, i.e., nucleotide exchanges within the transcribed and translated receptor sequence, may result in a difference in amino acid sequence thus altering receptor function and signaling. Such variants have been intensely studied in vitro in overexpression systems and addressed in candidate-gene studies for distinct clinical parameters. In recent years, large cohorts were analyzed in genome-wide association studies (GWAS), where variants are detected as significant in context with specific traits. These studies identified two of the in-depth characterized 18 coding variants in adrenoceptors as repeatedly statistically significant genetic risk factors - p.Arg389Gly in the β1- and p.Thr164Ile in the β2-adrenoceptor, along with 56 variants in the non-coding regions adjacent to the adrenoceptor gene loci, the functional role of which is largely unknown at present. This chapter summarizes current knowledge on the two coding variants in adrenoceptors that have been consistently validated in GWAS and provides a prospective overview on the numerous non-coding variants more recently attributed to adrenoceptor gene loci.

Synovial Tissue Heterogeneity in Japanese Patients With Rheumatoid Arthritis Elucidated Using a Cell-Type Deconvolution Approach

OBJECTIVE

Recent advances in single-cell RNA sequencing technology have improved our understanding of the immunological landscape of rheumatoid arthritis (RA). We aimed to stratify the synovium from East Asian patients with RA by immune cell compositions and gain insight into the inflammatory drivers of each synovial phenotype.

METHODS

Synovial tissues were obtained from East Asian patients in Japan with RA (n = 41) undergoing articular surgery. The cellular composition was quantified by a deconvolution approach using a public single-cell-based reference. Inflammatory pathway activity was calculated by gene set variation analysis, and chromatin accessibility was evaluated using assay of transposase accessible chromatin-sequencing.

RESULTS

We stratified RA synovium into three distinct subtypes based on the hierarchical clustering of cellular composition data. One subtype was characterized by abundant HLA-DRAhigh synovial fibroblasts, autoimmune-associated B cells, GZMK+ GZMB+ CD8+ T cells, interleukin (IL)1-β+ monocytes, and plasmablasts. In addition, tumor necrosis factor (TNF)-α, interferons (IFNs), and IL-6 signaling were highly activated in this subtype, and the expression of various chemokines was significantly enhanced. Moreover, we found an open chromatin region overlapping with RA risk locus rs9405192 near the IRF4 gene, suggesting the genetic background influences the development of this inflammatory synovial state. The other two subtypes were characterized by increased IFNs and IL-6 signaling, and expression of molecules associated with degeneration, respectively.

CONCLUSION

This study adds insights into the synovial heterogeneity in East Asian patients and shows a promising link with predominant inflammatory signals. Evaluating the site of inflammation has the potential to lead to appropriate drug selection that matches the individual pathology.

Identification of heel bone mineral density as a risk factor of Alzheimer's disease by analyzing large-scale genome-wide association studies datasets

Introduction: Both low bone mineral density (BMD) and Alzheimer's disease (AD) commonly co_occur in the older adult. Until now, the association between AD and BMD has been widely reported by observational studies. However, Mendelian randomization (MR) studies did not support the causal association between BMD and AD. We think that the lack of significant causal association between AD and BMD identified by recent MR studies may be caused by small number of potential instrumental variables. Methods: We conduct a MR study to evaluate the causal effect of heel BMD on the risk of AD using 1,362 genome-wide significant and independent (p < 5.00E-08) heel BMD genetic variants as the potential instrumental variables, which are identified by a large-scale genome wide association study (GWAS) of heel BMD in 394,929 UK Biobank individuals. Using these 1,362 genome-wide significant and independent heel BMD genetic variants, we extracted their corresponding AD GWAS summary results in IGAP AD GWAS dataset (n = 63,926) and FinnGen AD GWAS dataset (n = 377,277). Five methods including inverse-variance weighted meta-analysis (IVW), weighted median, MR-Egger, MR-PRESSO, and MRlap were selected to perform the MR analysis. 951 of these 1,362 genetic variants are available in AD GWAS dataset. Results: We observed statistically significant causal effect of heel BMD on the risk of AD using IVW in IGAP AD GWAS dataset (OR = 1.048, 95%CI: 1.002-1.095, p = 0.04) and FinnGen AD GWAS dataset (OR = 1.053, 95% CI:1.011-1.098, p = 0.011). Importantly, meta-analysis of IVW estimates from IGAP and FinnGen further supported the causal effect of heel BMD on the risk of AD (OR = 1.051, 95% CI: 1.02-1.083, p = 0.0013). Discussion: Collectively, our current MR study supports heel BMD to be a risk factor of AD by analyzing the large-scale heel BMD and AD GWAS datasets. The potential mechanisms underlying the association between heel BMD and AD should be further evaluated in future.

Genetic Prediction of Lifetime Risk of Fracture

CONTEXT

Fragility fracture is a significant public health problem because it is associated with increased mortality. We want to find out whether the risk of fracture can be predicted from the time of birth.

OBJECTIVE

To examine the association between a polygenic risk score (PRS) and lifetime fracture risk.

METHODS

This population-based prospective study involved 3515 community-dwelling individuals aged 60+ years who have been followed for up to 20 years. Femoral neck bone mineral density (BMD) was measured by dual-energy x-ray absorptiometry. A PRS was created by summing the weighted number of risk alleles for each single nucleotide polymorphism using BMD-associated coefficients. Fragility fractures were radiologically ascertained, whereas mortality was ascertained through a state registry. Residual lifetime risk of fracture (RLRF) was estimated by survival analysis.

RESULTS

The mortality-adjusted RLRF for women and men was 36% (95% CI, 34%-39%) and 21% (18%-24%), respectively. Individuals with PRS > 4.24 (median) had a greater risk (1.2-fold in women and 1.1-fold in men) than the population average risk. For hip fracture, the average RLRF was 10% (95% CI, 8%-12%) for women and ∼5% (3%-7%) for men; however, the risk was significantly increased by 1.5-fold and 1.3-fold for women and men with high PRS, respectively.

CONCLUSION

A genetic profiling of BMD-associated genetic variants is associated with the residual lifetime risk of fracture, suggesting the potential for incorporating the polygenic risk score in personalized fracture risk assessment.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits

BACKGROUND

Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits.

RESULTS

We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19).

CONCLUSIONS

These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

Are commercial genetic injury tests premature?

INTRODUCTION

Several direct-to-consumer (DTC) genetic testing companies have emerged that claim to be able to test for susceptibility for musculoskeletal injuries. Although there are several publications on the emergence of this industry, none have critically evaluated the evidence for the use of genetic polymorphisms in commercial tests. The aim of this review was to identify, where possible, the polymorphisms and to evaluate the current scientific evidence for their inclusion.

RESULTS

The most common polymorphisms included COL1A1 rs1800012, COL5A1 rs12722, and GDF5 rs143383. The current evidence suggests that it is premature or even not viable to include these three polymorphisms as markers of injury risk. A unique set of injury-specific polymorphisms, which do not include COL1A1, COL5A1, or GDF5, identified from genome-wide association studies (GWAS) is used by one company in their tests for 13 sports injuries. However, of the 39 reviewed polymorphisms, 22 effective alleles are rare and absent in African, American, and/or Asian populations. Even when informative in all populations, the sensitivity of many of the genetic markers was low and/or has not been independently validated in follow-up studies.

CONCLUSIONS

The current evidence suggests it is premature to include any of the reviewed polymorphisms identified by GWAS or candidate gene approaches in commercial genetic tests. The association of MMP7 rs1937810 with Achilles tendon injuries, and SAP30BP rs820218 and GLCCI1 rs4725069 with rotator cuff injuries does warrant further investigation. Based on current evidence, it remains premature to market any commercial genetic test to determine susceptibility to musculoskeletal injuries.

Identification of an intronic enhancer regulating RANKL expression in osteocytic cells

The bony skeleton is continuously renewed throughout adult life by the bone remodeling process, in which old or damaged bone is removed by osteoclasts via largely unknown mechanisms. Osteocytes regulate bone remodeling by producing the osteoclast differentiation factor RANKL (encoded by the TNFSF11 gene). However, the precise mechanisms underlying RANKL expression in osteocytes are still elusive. Here, we explored the epigenomic landscape of osteocytic cells and identified a hitherto-undescribed osteocytic cell-specific intronic enhancer in the TNFSF11 gene locus. Bioinformatics analyses showed that transcription factors involved in cell death and senescence act on this intronic enhancer region. Single-cell transcriptomic data analysis demonstrated that cell death signaling increased RANKL expression in osteocytic cells. Genetic deletion of the intronic enhancer led to a high-bone-mass phenotype with decreased levels of RANKL in osteocytic cells and osteoclastogenesis in the adult stage, while RANKL expression was not affected in osteoblasts or lymphocytes. These data suggest that osteocytes may utilize a specialized regulatory element to facilitate osteoclast formation at the bone surface to be resorbed by linking signals from cellular senescence/death and RANKL expression.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Preliminary study on genetic factors related to Demirjian's tooth age estimation method based on genome-wide association analysis

The age determination of individuals, especially minors, is critical in forensic research. In forensic practice, dental age estimation is one of the most commonly used methods for determining age as teeth are easy to preserve and relatively resistant to environmental factors. Tooth development is affected and regulated by genetic factors; however, these are not incorporated into current commonly used tooth age inference methods, leading to unreliable results. Here, we established a Demirjian and a Cameriere tooth age estimation-based methods suitable for use in children in southern China. By using the difference between the inferred age and the actual age (MD) as the phenotype, we identified 65 and 49 SNPs related to tooth age estimation from 743,722 loci among 171 children in southern China through a genome-wide association analysis (p<0.0001). We also conducted a genome-wide association study on dental development stage (DD) using the Demirjian tooth age estimation method and screened two sets of SNP sites (52 and 26) based on whether age difference was considered. The gene function enrichment analysis of these SNPs found that they were related to bone development and mineralization. Although SNP sites screened based on MD seem to improve the accuracy of tooth age estimation, there is little correlation between these SNPs and an individual's Demirjian morphological stage. In conclusion, we found that individual genotypes can affect tooth age estimation, and based on different phenotypic analysis models, we have identified some novel SNP sites related to tooth age inference and Demirjian's tooth development stage. These studies provide a reference for subsequent phenotypic selection based on tooth age inference analysis, and the results could possibly be used in the future to make forensic age estimation more accurate.

The clinical utility of the BMD-related comprehensive genome-wide polygenic score in identifying individuals with a high risk of osteoporotic fractures

The potential of bone mineral density (BMD)-related genome-wide polygenic score (PGS) in identifying individuals with a high risk of fractures remains unclear. This study suggests that an efficient PGS enables the identification of strata with up to a 1.5-fold difference in fracture incidence. Incorporating PGS into clinical diagnosis is anticipated to increase the population-level screening benefits.

PURPOSE

This study sought to construct genome-wide polygenic scores for femoral neck and total body BMD and to estimate their potential in identifying individuals with a high risk of osteoporotic fractures.

METHODS

Genome-wide polygenic scores were developed and validated for femoral neck and total body BMD. We externally tested the PGSs, both by themselves and in combination with available clinical risk factors, in 455,663 European ancestry individuals from the UK Biobank. The predictive accuracy of the developed genome-wide PGS was also compared with previously published restricted PGS employed in fracture risk assessment.

RESULTS

For each unit decrease in PGSs, the genome-wide PGSs were associated with up to 1.17-fold increased fracture risk. Out of four studied PGSs, [Formula: see text] (HR: 1.03; 95%CI 1.01-1.05, p = 0.001) had the weakest and the [Formula: see text] (HR: 1.17; 95%CI 1.15-1.19, p < 0.0001) had the strongest association with an incident fracture. In the reclassification analysis, compared to the FRAX base model, the models with [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] improved the reclassification of fracture by 1.2% (95% CI, 1.0 to 1.3%), 0.2% (95% CI, 0.1 to 0.3%), 1.4% (95% CI, 1.3 to 1.5%), and 2.2% (95% CI, 2.1 to 2.4%), respectively.

CONCLUSIONS

Our findings suggested that an efficient PGS estimate enables the identification of strata with up to a 1.7-fold difference in fracture incidence. Incorporating PGS information into clinical diagnosis is anticipated to increase the benefits of screening programs at the population level.

Association of plant-based dietary patterns with the risk of osteoporosis in community-dwelling adults over 60 years: a cross-sectional study

Plant foods are rich in many important micronutrients that are beneficial for bone health. This cross-sectional study of 9613 community-dwelling older adults found that more consumption of healthy plant foods and less consumption of animal foods and unhealthy plant foods were associated with a lower risk of osteoporosis.

INTRODUCTION

Osteoporosis,a common chronic disease in older adults, threatens their health. Many nutrients in plant foods are important for preventing osteoporosis. However, the relationship between plant-based dietary patterns and osteoporosis remains unclear. This study aimed to investigate the cross-sectional association between plant-based dietary patterns and osteoporosis in older adults.

METHODS

This study was conducted among 9613 community-dwelling older adults in Liaoning Province, China. The effective food frequency questionnaire (FFQ) and plant-based diet index (PDI) were used to evaluate compliance with plant-based dietary patterns. Osteoporosis was defined based on heel ultrasound. We analyzed the association between healthy plant-based diet index (HPDI) and unhealthful plant-based diet index (UPDI) and the risk of osteoporosis.

RESULTS

A higher PDI was associated with higher bone mineral density (BMD) in older adults. In logistic regression models, the highest quartile of PDI and HPDI had a significantly lower risk of osteoporosis than the lowest quartile, whereas UPDI in the highest quartile was associated with a higher risk of osteoporosis. There was a dose-response relationship between the three indices and the risk of osteoporosis. Subgroup analysis revealed differences in the relationship between HPDI and the risk of osteoporosis according to gender and age.

CONCLUSIONS

Older adults, especially women, consume more healthy plant foods and reduce the consumption of animal foods and unhealthy plant foods, which was associated with a lower risk of osteoporosis.

Exome-wide screening identifies novel rare risk variants for bone mineral density

Bone mineral density (BMD) is an independent risk factor of osteoporosis-related fractures. We performed gene-based burden tests to assess the association between rare variants and BMD, and identified several BMD candidate genes.

PURPOSE

BMD is highly heritable and a major predictor of osteoporotic fractures, but its genetic basis remains unclear. We aimed to identify rare risk variants contributing to BMD.

METHODS

Utilizing the newly released UK Biobank 200,643 exome dataset, we conducted a gene-based exome-wide association study in males and females, respectively. First, 100,639 males and 117,338 females with BMD values were included in the polygenic risk scores (PRS) analysis. Among individuals with lower 30% PRS, cases were individuals with top 10% BMD, and individuals with bottom 10% BMD were the controls. Considering the effects of vitamin D (VD), individuals with the highest 30% VD concentration were selected for VD-BMD analysis. After quality control, 741 males and 697 females were included in the BMD analysis, and 717 males and 708 females were included in the VD-BMD analysis. The variants were annotated by ANNOVAR software, then BMD and VD-BMD qualified variants were imported into the SKAT R-package to perform gene-based burden tests, respectively.

RESULTS

The gene-based burden test of the exonic variants identified genome-wide candidate associations in ANKRD18A (P = 1.60 × 10^-5, P_{Bonferroni adjust} = 2.11 × 10^-3), C22orf31 (P = 3.49 × 10^-4, P_{Bonferroni adjust} = 3.17 × 10^-2), and SPATC1L (P = 1.09 × 10^-5, P_{Bonferroni adjust} = 8.80 × 10^-3). For VD-BMD analysis, three genes were associated with BMD, such as NIPAL1 (P = 1.06 × 10^-3, P_{Bonferroni adjust} = 3.91 × 10^-2).

CONCLUSIONS

Our study suggested that rare variants contribute to BMD, providing new sights for broadening the genetic structure of BMD.

Genome-wide polygenic risk score for major osteoporotic fractures in postmenopausal women using associated single nucleotide polymorphisms

BACKGROUND

Osteoporosis is highly polygenic and heritable, with heritability ranging from 50 to 80%; most inherited susceptibility is associated with the cumulative effect of many common genetic variants. However, existing genetic risk scores (GRS) only provide a few percent predictive power for osteoporotic fracture.

METHODS

We derived and validated a novel genome-wide polygenic score (GPS) comprised of 103,155 common genetic variants to quantify this susceptibility and tested this GPS prediction ability in an independent dataset (n = 15,776).

RESULTS

Among postmenopausal women, we found a fivefold gradient in the risk of major osteoporotic fracture (MOF) (p < 0.001) and a 15.25-fold increased risk of severe osteoporosis (p < 0.001) across the GPS deciles. Compared with the remainder of the GPS distribution, the top GPS decile was associated with a 3.59-, 2.48-, 1.92-, and 1.58-fold increased risk of any fracture, MOF, hip fracture, and spine fracture, respectively. The top GPS decile also identified nearly twofold more high-risk osteoporotic patients than the top decile of conventional GRS based on 1103 conditionally independent genome-wide significant SNPs. Although the relative risk of severe osteoporosis for postmenopausal women at around 50 is relatively similar, the cumulative incident at 20-year follow-up is significantly different between the top GPS decile (13.7%) and the bottom decile (< 1%). In the subgroup analysis, the GPS transferability in non-Hispanic White is better than in other racial/ethnic groups.

CONCLUSIONS

This new method to quantify inherited susceptibility to osteoporosis and osteoporotic fracture affords new opportunities for clinical prevention and risk assessment.

Multiple solutions at the genomic level in response to selective breeding for high locomotor activity

Replicate lines under uniform selection often evolve in different ways. Previously, analyses using whole-genome sequence data for individual mice (Mus musculus) from 4 replicate High Runner lines and 4 nonselected control lines demonstrated genomic regions that have responded consistently to selection for voluntary wheel-running behavior. Here, we ask whether the High Runner lines have evolved differently from each other, even though they reached selection limits at similar levels. We focus on 1 High Runner line (HR3) that became fixed for a mutation at a gene of major effect (Myh4Minimsc) that, in the homozygous condition, causes a 50% reduction in hindlimb muscle mass and many pleiotropic effects. We excluded HR3 from SNP analyses and identified 19 regions not consistently identified in analyses with all 4 lines. Repeating analyses while dropping each of the other High Runner lines identified 12, 8, and 6 such regions. (Of these 45 regions, 37 were unique.) These results suggest that each High Runner line indeed responded to selection somewhat uniquely, but also that HR3 is the most distinct. We then applied 2 additional analytical approaches when dropping HR3 only (based on haplotypes and nonstatistical tests involving fixation patterns). All 3 approaches identified 7 new regions (as compared with analyses using all 4 High Runner lines) that include genes associated with activity levels, dopamine signaling, hippocampus morphology, heart size, and body size, all of which differ between High Runner and control lines. Our results illustrate how multiple solutions and "private" alleles can obscure general signatures of selection involving "public" alleles.

Biallelic frameshift variants inPHLDB1cause mild-type osteogenesis imperfecta with regressive spondylometaphyseal changes

Background

Osteogenesis imperfecta (OI) is a heterogeneous group of inherited disorders characterised by susceptibility to fractures, primarily due to defects in type 1 collagen. The aim of this study is to present a novel OI phenotype and its causative candidate gene.

Methods

Whole-exome sequencing and clinical evaluation were performed in five patients from two unrelated families.PHLDB1mRNA expression in blood and fibroblasts was investigated by real-time PCR, and western blot analysis was further performed on skin fibroblasts.

Results

The common findings among the five affected children were recurrent fractures and/or osteopaenia, platyspondyly, short and bowed long bones, and widened metaphyses. Metaphyseal and vertebral changes regressed after early childhood, and no fractures occurred under bisphosphonate treatment. We identified biallelic NM_001144758.3:c.2392dup and NM_001144758.3:c.2690_2693del pathogenic variants inPHLDB1in the affected patients, respectively, in the families; parents were heterozygous for these variants.PHLDB1encodes pleckstrin homology-like domain family B member-1 (PHLDB1) protein, which has a role in insulin-dependent Akt phosphorylation. Compared with controls, a decrease in the expression levels ofPHLDB1in the blood and skin fibroblast samples was detected. Western blot analysis of cultured fibroblasts further confirmed the loss of PHLDB1.

Conclusion

Two biallelic frameshift variants in the candidate genePHLDB1were identified in independent families with a novel, mild-type, autosomal recessive OI. The demonstration of decreasedPHLDB1mRNA expression levels in blood and fibroblast samples supports the hypothesis thatPHLDB1pathogenic variants are causative for the observed phenotype.

Association of air pollutants and osteoporosis risk: The modifying effect of genetic predisposition

BACKGROUND

Limited studies have examined the association between air pollutants and osteoporosis incidence; however, the results are conflicting. We aimed to quantify the effects of selected air pollutants on osteoporosis risk and explore the modifying effect of genetic predisposition.

METHODS

A total of 422,955 subjects who did not have osteoporosis at baseline in the UK Biobank were included from 2006 to 2010. We conducted a Cox proportional hazards model with adjustment for covariates to examine the association between air pollutant scores and individual air pollutants and incident osteoporosis. Furthermore, a polygenic risk score (PRS) of osteoporosis was built and examined to determine whether genetic susceptibility modified the effect of air pollutants on osteoporosis. The relationship between air pollutants and osteoporosis was examined by using a restricted cubic spline (RCS) method.

RESULTS

After confounder adjustment, the results showed a remarkable increase in the risk of osteoporosis with each 10 unit increase in exposure to air pollution (hazard ratio: 1.06, 95 % confidence interval: 1.03-1.08), PM_2.5 (1.94, 1.52-2.48), NO₂ (1.06, 1.02-1.10), and NO_X (1.03, 1.01-1.04). However, no significant association was observed between PM₁₀ or PM_2.5-10 exposure and osteoporosis. Subjects with high air pollutant exposure levels and a high PRS had a noteworthy increase in osteoporosis risk compared to those with low air pollutant exposure levels and a low PRS. Air pollutants and genetic variants exerted additive effects on the risk of osteoporosis. Positive correlations were observed between osteoporosis and PM_2.5 (P < 0.001), NO₂ (P = 0.001), and NOx (P = 0.002) exposure.

CONCLUSIONS

Exposure to PM_2.5, NO₂ and NOx was associated with an increase in osteoporosis risk, and this effect was more pronounced in populations with high genetic risk. The association between PM_2.5, NO₂ and NOx exposure and osteoporosis is modified by genetic variations.

Targeted Resequencing of Otosclerosis Patients from Different Populations Replicates Results from a Previous Genome-Wide Association Study

Otosclerosis is one of the most common causes of hearing loss in young adults. It has a prevalence of 0.3-0.4% in the European population. Clinical symptoms usually occur between the second and fifth decade of life. Different studies have been performed to unravel the genetic architecture of the disease. Recently, a genome-wide association study (GWAS) identified 15 novel risk loci and replicated the regions of three previously reported candidate genes. In this study, seven candidate genes from the GWAS were resequenced using single molecule molecular inversion probes (smMIPs). smMIPs were used to capture the exonic regions and the 3' and 5' untranslated regions (UTR). Discovered variants were tested for association with the disease using single variant and gene-based association analysis. The single variant results showed that 13 significant variants were associated with otosclerosis. Associated variants were found in five of the seven genes studied here, including AHSG, LINC01482, MARK3, SUPT3H and RELN. Conversely, burden testing did not show a major role of rare variants in the disease. In conclusion, this study was able to replicate five out of seven candidate genes reported in the previous GWAS. This association is likely mainly driven by common variants.

Sex Hormone Candidate Gene Polymorphisms Are Associated with Endometriosis

The present study was designed to examine whether sex hormone polymorphisms proven by GWAS are associated with endometriosis risk. Unrelated female participants totaling 1376 in number (395 endometriosis patients and 981 controls) were recruited into the study. Nine single-nucleotide polymorphisms (SNPs) which GWAS correlated with circulating levels of sex hormones were genotyped using a TaqMan allelic discrimination assay. FSH-lowering, and LH- and testosterone-heightening polymorphisms of the FSHB promoter (allelic variants A rs11031002 and C rs11031005) exhibit a protective effect for endometriosis (OR = 0.60-0.68). By contrast, the TT haplotype loci that were GWAS correlated with higher FSH levels and lower LH and testosterone concentrations determined an increased risk for endometriosis (OR = 2.03). Endometriosis-involved epistatic interactions were found between eight loci of sex hormone genes (without rs148982377 ZNF789) within twelve genetic simulation models. In silico examination established that 8 disorder-related loci and 80 proxy SNPs are genome variants affecting the expression, splicing, epigenetic and amino acid conformation of the 34 genes which enrich the organic anion transport and secondary carrier transporter pathways. In conclusion, the present study showed that sex hormone polymorphisms proven by GWAS are associated with endometriosis risk and involved in the molecular pathophysiology of the disease due to their functionality.

Associations of bone mass and polygenic risk of osteoporosis with indicators of arterial wall condition

Background: The identification of genetic factors that are simultaneously responsible for the predisposition to the development of cardiovascular diseases (CVD) and osteoporosis (OP) is important for the prevention of both conditions.

Aim: The aim of this study is to evaluate three genetic risk scales (GRS) that previously showed an association with bone mineral density (BMD) and fracture risk, as well as to study the associations of these GRS with vascular wall pathology.

Materials and methods: 250 female outpatients (aged 45 to 69) were enrolled into a cross-sectional study. The intima-media thickness (IMT), the presence and number of atherosclerotic plaques (AP) were studied using duplex scanning. Pulse wave velocity (PWV), augmentation index (AI) were measured by applanation tonometry. Coronary vessels calcium deposits were registered by multispiral computed tomography (MSCT) using the Agatston calcium index (CI). The BMD of the spine, hip neck (HN) and proximal hip (PH) was measured using double energy x-ray absorptiometry. Bone resorption marker type-1 collagen C-terminal telopeptide (CTx) was assessed solid-phase enzyme immunoassay. The genetic study included DNA extraction from whole blood samples. Targeted sequencing was performed on the Nextseq550 sequencer (Illumina, USA). Statistical analysis was carried out using the SAS software package for Windows, version 9.0 (SAS Institute Inc., USA).

Results: The chance of detecting low bone mass increased more than 4 times at values of IMT ≥0.9 mm (OR=4.17; 95%CI [1.2–14.4], p<0.02), 2.4 times in the presence of AP in the carotid arteries (OR=2.45; 95%CI [1.12–4.88], p><0.05), by 6.7 times with an Agatstone CI ≥ 100 units (OR=6.68; 95%CI [1.56–28.7], p><0.001), 1.4 times (OR=1.43; 95%CI [0.56–3.68], p><0.438) with a PWV ≥10 m/s, 1.2 times (OR=1.2; 95%CI [0.601–2.43], p><0.60) with increased AI ≥ 27%. According to multivariate linear regression analysis (adjusted for age, duration of postmenopause, marker of bone resorption CTx), a significant association of all GRS with BMD in all parts of the skeleton was revealed. Both univariate and multivariate regression models adjusted for several covariants (age, total cholesterol, systolic blood pressure) showed a reliable association of GRS62 with the presence of plaques and GRS63 — with coronary artery CI. Conclusion: The results of the study demonstrated the association of polygenic genetic risk of GRS-based OP with BMD and vascular wall status indicators in women in the peri and postmenopausal periods.>< 0.02), 2.4 times in the presence of AP in the carotid arteries (OR=2.45; 95%CI [1.12–4.88], p< 0.05), by 6.7 times with an Agatstone CI ≥ 100 units (OR=6.68; 95%CI [1.56–28.7], p< 0.001), 1.4 times (OR=1.43; 95%CI [0.56–3.68], p< 0.438) with a PWV ≥10 m/s, 1.2 times (OR=1.2; 95%CI [0.601–2.43], p<0.60) with increased AI ≥ 27%. According to multivariate linear regression analysis (adjusted for age, duration of postmenopause, marker of bone resorption CTx), a significant association of all GRS with BMD in all parts of the skeleton was revealed. Both univariate and multivariate regression models adjusted for several covariants (age, total cholesterol, systolic blood pressure) showed a reliable association of GRS62 with the presence of plaques and GRS63 — with coronary artery CI.>< 0.60) with increased AI ≥ 27%. According to multivariate linear regression analysis (adjusted for age, duration of postmenopause, marker of bone resorption CTx), a significant association of all GRS with BMD in all parts of the skeleton was revealed. Both univariate and multivariate regression models adjusted for several covariants (age, total cholesterol, systolic blood pressure) showed a reliable association of GRS62 with the presence of plaques and GRS63 — with coronary artery CI.

Conclusion: The results of the study demonstrated the association of polygenic genetic risk of GRS-based OP with BMD and vascular wall status indicators in women in the peri and postmenopausal periods.

Using a Polygenic Score to Predict the Risk of Developing Primary Osteoporosis

Osteoporosis (OP) is a multifactorial bone disease belonging to the metabolic osteopathies group. Using the polygenic score (PGS) approach, we combined the effects of bone mineral density (BMD) DNA loci, affecting osteoporosis pathogenesis, based on GEFOS/GENOMOS consortium GWAS meta-analysis. We developed models to predict the risk of low fractures in women from the Volga-Ural region of Russia with efficacy of 74% (AUC = 0.740; OR (95% CI) = 2.9 (2.353–3.536)), as well as the formation of low BMD with efficacy of 79% (AUC = 0.790; OR (95% CI) = 3.94 (2.993–5.337)). In addition, we propose a model that predicts fracture risk and low BMD in a comorbid condition with 85% accuracy (AUC = 0.850; OR (95% CI) = 6.6 (4.411–10.608)) in postmenopausal women.

Multiple polygenic scores improve bone mineral density prediction in an independent sample of Caucasian women

PURPOSE OF THE STUDY

To determine if multiple Genetic Risk Scores (GRSs) improve bone mineral density (BMD) prediction over single GRS in an independent sample of Caucasian women.

STUDY DESIGN

Based on summary statistics of four genome-wide association studies related to two osteoporosis-associated traits, namely BMD and heel quantitative ultrasound derived estimated BMD (eBMD), four GRSs were derived for 1205 individuals in the Genome-Wide Scan for Female Osteoporosis Gene Study. The effect of each GRS on BMD variation was assessed using multivariable linear regression, with conventional risk factors adjusted for. Next, the eBMD-related GRS that explained the most variance in BMD was selected to be entered into a multi-score model, along with the BMD-related GRS. Elastic net regularised regression was used to develop the multiscore model, which estimated the joint effect of two GRSs (GRS_BMD and GRS_eBMD) on BMD variation, after being adjusted for conventional risk factors.

RESULTS

With the same clinical risk factors having been adjusted for, the model that included GRS_BMD performed best by explaining 32.53% of the variance in BMD; the single-score model that included GRS_eBMD explained 34.03% of BMD variance. The model that includes both GRS_BMD and GRS_ eBMD, as well as the clinical risk factors, aggregately explained 35.05% in BMD variation. Compared with the single GRS models, the multiscore model explained significantly more variance in BMD.

CONCLUSIONS

The multipolygenic score model explained a considerable amount of BMD variation. Compared with single score models, multipolygenic score model provided significant improvement in explaining BMD variation.

FGFRL1 and FGF genes are associated with height, hypertension, and osteoporosis

Hypertension and osteoporosis are two major disorders, which interact with each other. Specific genetic signals involving the fibroblast growth factor receptor-like 1 (FGFRL1) gene are related to high blood pressure and bone growth in giraffes. FGFRL1 is associated with cardiovascular system and bone formation. We performed an association study to investigate the role of FGFRL1 in hypertension, osteoporosis, and height determination in humans. In addition, we identified three kinds of phenotypes in fibroblast growth factor (FGF) genes and examined their association with the FGFRL1 gene. We identified 42 SNPs in the FGFRL1 gene associated with each trait. We then analyzed the potential functional annotation of each SNP. The FGFRL1 gene was found to be associated with height, hypertension, and osteoporosis, consistent with the results of a previous study. In addition, the FGF2, FGF4, FGF10, FGF18, and FGF22 genes were found to interact with the FGFRL1 gene. Our study suggests that both FGFRL1 and FGFRL1-related genes may determine the height and the prevalence of osteoporosis and hypertension in the Korean population.

Screening of Mineralogenic and Osteogenic Compounds in Zebrafish—Tools to Improve Assay Throughput and Data Accuracy

Bone disorders affect millions of people worldwide and treatments currently available often produce undesirable secondary effects or have limited efficacy. It is therefore of the utmost interest for patients to develop more efficient drugs with reduced off-target activities. In the long process of drug development, screening and preclinical validation have recently gained momentum with the increased use of zebrafish as a model organism to study pathological processes related to human bone disorders, and the development of zebrafish high-throughput screening assays to identify bone anabolic compounds. In this review, we provided a comprehensive overview of the literature on zebrafish bone-related assays and evaluated their performance towards an integration into screening pipelines for the discovery of mineralogenic/osteogenic compounds. Tools available to standardize fish housing and feeding procedures, synchronize embryo production, and automatize specimen sorting and image acquisition/analysis toward faster and more accurate screening outputs were also presented.

Oxytocin receptor expression patterns in the human brain across development

Oxytocin plays a vital role in social behavior and homeostatic processes, with animal models indicating that oxytocin receptor (OXTR) expression patterns in the brain influence behavior and physiology. However, the developmental trajectory of OXTR gene expression is unclear. By analyzing gene expression data in human post-mortem brain samples, from the prenatal period to late adulthood, we demonstrate distinct patterns of OXTR gene expression in the developing brain, with increasing OXTR expression along the course of the prenatal period culminating in a peak during early childhood. This early life OXTR expression peak pattern appears slightly earlier in a comparative macaque sample, which is consistent with the relative immaturity of the human brain during early life compared to macaques. We also show that a network of genes with strong spatiotemporal couplings with OXTR is enriched in several psychiatric illness and body composition phenotypes. Taken together, these results demonstrate that oxytocin signaling plays an important role in a diverse set of psychological and somatic processes across the lifespan.

Absence of causal association between Vitamin D and bone mineral density across the lifespan: a Mendelian randomization study

Vitamin D deficiency is a candidate risk factor for osteoporosis, characterized by decreased bone mineral density (BMD). We performed this two-sample Mendelian randomization (MR) analysis to investigate the causal effect of vitamin D on BMD. We extracted 143 single-nucleotide polymorphisms from a recent GWAS on 417,580 participants of European ancestry as instrumental variables, and used summary statistics for BMD at forearm (n = 10,805), femoral neck (n = 49,988), lumbar spine (n = 44,731) and total-body of different age-stages (< 15, 15-30, 30-45, 45-60, > 60) (n = 67,358). We explored the direct effect of vitamin D on BMD with an adjusted body mass index (BMI) in a multivariable MR analysis. We found no support for causality of 25-hydroxyvitamin D on BMD at forearm, femoral neck, lumbar spine, and total-body BMD across the lifespan. There was no obvious difference between the total and direct effect of vitamin D on BMD after adjusting for BMI. Our MR analysis provided evidence that genetically determined vitamin D was not causally associated with BMD in the general population. Large-scale randomized controlled trials are warranted to investigate the role of vitamin D supplementation in preventing osteoporosis in the high-risk population.

Genetic Diagnostics in Routine Osteological Assessment of Adult Low Bone Mass Disorders

CONTEXT

Many different inherited and acquired conditions can result in premature bone fragility/low bone mass disorders (LBMDs).

OBJECTIVE

We aimed to elucidate the impact of genetic testing on differential diagnosis of adult LBMDs and at defining clinical criteria for predicting monogenic forms.

METHODS

Four clinical centers broadly recruited a cohort of 394 unrelated adult women before menopause and men younger than 55 years with a bone mineral density (BMD) Z-score < -2.0 and/or pathological fractures. After exclusion of secondary causes or unequivocal clinical/biochemical hallmarks of monogenic LBMDs, all participants were genotyped by targeted next-generation sequencing.

RESULTS

In total, 20.8% of the participants carried rare disease-causing variants (DCVs) in genes known to cause osteogenesis imperfecta (COL1A1, COL1A2), hypophosphatasia (ALPL), and early-onset osteoporosis (LRP5, PLS3, and WNT1). In addition, we identified rare DCVs in ENPP1, LMNA, NOTCH2, and ZNF469. Three individuals had autosomal recessive, 75 autosomal dominant, and 4 X-linked disorders. A total of 9.7% of the participants harbored variants of unknown significance. A regression analysis revealed that the likelihood of detecting a DCV correlated with a positive family history of osteoporosis, peripheral fractures (> 2), and a high normal body mass index (BMI). In contrast, mutation frequencies did not correlate with age, prevalent vertebral fractures, BMD, or biochemical parameters. In individuals without monogenic disease-causing rare variants, common variants predisposing for low BMD (eg, in LRP5) were overrepresented.

CONCLUSION

The overlapping spectra of monogenic adult LBMD can be easily disentangled by genetic testing and the proposed clinical criteria can help to maximize the diagnostic yield.

Comprehensive characterization of posttranscriptional impairment-related 3'-UTR mutations in 2413 whole genomes of cancer patients

The 3' untranslated region (3'-UTR) is the vital element regulating gene expression, but most studies have focused on variations in RNA-binding proteins (RBPs), miRNAs, alternative polyadenylation (APA) and RNA modifications. To explore the posttranscriptional function of 3'-UTR somatic mutations in tumorigenesis, we collected whole-genome data from 2413 patients across 18 cancer types. Our updated algorithm, PIVar, revealed 25,216 3'-UTR posttranscriptional impairment-related SNVs (3'-UTR piSNVs) spanning 2930 genes; 24 related RBPs were significantly enriched. The somatic 3'-UTR piSNV ratio was markedly increased across all 18 cancer types, which was associated with worse survival for four cancer types. Several cancer-related genes appeared to facilitate tumorigenesis at the protein and posttranscriptional regulation levels, whereas some 3'-UTR piSNV-affected genes functioned mainly via posttranscriptional mechanisms. Moreover, we assessed immune cell and checkpoint characteristics between the high/low 3'-UTR piSNV ratio groups and predicted 80 compounds associated with the 3'-UTR piSNV-affected gene expression signature. In summary, our study revealed the prevalence and clinical relevance of 3'-UTR piSNVs in cancers, and also demonstrates that in addition to affecting miRNAs, 3'-UTR piSNVs perturb RBPs binding, APA and m6A RNA modification, which emphasized the importance of considering 3'-UTR piSNVs in cancer biology.

Genetic variation in WNT16 and its association with bone mineral density, fractures and osteoporosis in children with bone fragility

Several genome-wide association studies (GWAS), GWAS meta-analyses, and mouse studies have demonstrated that wingless-related integration site 16 (WNT16) gene is associated with bone mineral density (BMD), cortical bone thickness, bone strength and fracture risk. Practically no data exist regarding the significance of WNT16 in childhood-onset osteoporosis and related fractures. We hypothesized that pathogenic variants and genetic variations in WNT16 could explain skeletal fragility in affected children. We screened the WNT16 gene by Sanger sequencing in three pediatric cohorts: 35 with primary osteoporosis, 59 with multiple fractures, and in 95 healthy controls. Altogether, we identified 12 variants in WNT16. Of them one was a rare 5′UTR variant rs1386898215 in genome aggregate and medical trans-omic databases (GnomAD, TOPMED; minor allele frequency (MAF) 0.00 and 0.000008, respectively). One variant rs1554366753, overrepresented in children with osteoporosis (MAF = 0.06 vs healthy controls MAF = 0.01), was significantly associated with lower BMD. This variant was found associated with increased WNT16 gene expression at mRNA level in fibroblast cultures. None of the other identified variants were rare (MAF < 0.001) or deemed pathogenic by predictor programs. WNT16 may play a role in childhood osteoporosis but genetic WNT16 variation is not a common cause of skeletal fragility in childhood.

•

No pathogenic WNT16 variants were found associated with pediatric osteoporosis or fracture-prone patients

•

Altogether, twelve WNT16 variants were found in pediatric osteoporosis or fracture-prone patients

•

The genetic variation rs1554366753 in the WNT16 gene is associated with bone mineral density and primary osteoporosis

Detection of common copy number of variation underlying selection pressure in Middle Eastern horse breeds using whole-genome sequence data

Dareshouri, Arabian, and Akhal-Teke are three Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole-genome of four Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n=3) and Akhal-Teke (n=3) horses . Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size=100kb, step size=50kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8,837,950 single nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and three quantitative trait loci (QTLs), including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted p-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs which were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Shared Genetic Architecture Between Rheumatoid Arthritis and Varying Osteoporotic Phenotypes

Rheumatoid arthritis (RA) and low bone mineral density (BMD), an indicator of osteoporosis (OP), appear epidemiologically associated. Shared genetic factors may explain this association. This study aimed to investigate the presence of pleiotropy to clarify the potential genetic association between RA and OP. We examined BMDs at varying skeletal sites reported in UK Biobank as well as OP fracture acquired from the Genetic Factors for Osteoporosis (GEFOS) Consortium and the TwinsUK study. PRSice-2 was used to assess the potential shared genetic overlap between RA and OP. The presence of pleiotropy was examined using colocalization analysis. PRSice-2 revealed that RA was significantly associated with OP fracture (β = 351.6 ± 83.9, p value = 2.76E-05), total BMD (β = -1763.5 ± 612.8, p = 4.00E-03), spine BMD (β = -919.8 ± 264.6, p value = 5.09E-04), and forearm BMD (β = -66.09 ± 31.40, p value = 3.53E-02). Through colocalization analysis, the same causal genetic variants, associated with both RA and OP, were apparent in 12 genes: PLCL1, BOLL, AC011997.1, TNFAIP3, RP11-158I9.1, CDK6, CHCHD4P2, RP11-505C13.1, PHF19, TRAF1, C5, and C11orf49 with moderate posterior probabilities (>50%). Pleiotropy is involved in the association between RA and OP phenotypes. These findings contribute to the understanding of disease mechanisms and provide insight into possible therapeutic advancements and enhanced screening measures. © 2021 American Society for Bone and Mineral Research (ASBMR).

SLC4A2 Deficiency Causes a New Type of Osteopetrosis

Osteopetrosis is a group of rare inherited skeletal disorders characterized by a marked increase in bone density due to deficient bone resorption. Pathogenic variants in several genes involved in osteoclast differentiation and/or function have been reported to cause osteopetrosis. Solute carrier family 4 member 2 (SLC4A2, encoding anion exchanger 2) plays an important role in osteoclast differentiation and function by exchange of Cl^- with HCO₃^- . Biallelic Slc4a2 loss-of-function mutations in mice and cattle lead to osteopetrosis with osteoclast deficiency; however, pathogenic SLC4A2 variants in humans have not been reported. In this study, we describe a patient with autosomal recessive osteopetrosis due to biallelic pathogenic variants in SLC4A2. We identified novel compound heterozygous variants in SLC4A2 (NM_003040.4: c.556G>A [p.A186T] and c.1658T>C [p.V553A]) by exome sequencing. The measurement of intracellular Cl^- showed that the variants decrease the anion exchange activity of SLC4A2. The impact of the variants on osteoclast differentiation was assessed by a gene knockout-rescue system using a mouse macrophage cell line, RAW 264.7. The Slc4a2-knockout cells show impaired osteoclastogenesis, which was rescued by the wild-type SLC4A2, but not by the mutant SLC4A2s. Immunofluorescence and pit assay revealed that the mutant SLC4A2s leads to abnormal podosome belt formation with impaired bone absorption. This is the first report on an individual affected by SLC4A2-associated osteopetrosis (osteopetrosis, Ikegawa type). With functional studies, we prove that the variants lead to SLC4A2 dysfunction, which altogether supports the importance of SLC4A2 in human osteoclast differentiation. © 2021 American Society for Bone and Mineral Research (ASBMR).

Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations

Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.