Search

Abstract Sarcopenia is a common skeletal muscle disease in older people. Lower limb muscle strength is a good predictive value for sarcopenia; however, little is known about its genetic components. Here, we conducted a genome-wide association study (GWAS) for knee extension strength in a total of 3452 Japanese aged 60 years or older from two independent cohorts. We identified a significant locus, rs10749438 which is an intronic variant in TACC2 (transforming acidic coiled-coil-containing 2) (P = 4.2 × 10−8). TACC2, encoding a cytoskeleton-related protein, is highly expressed in skeletal muscle, and is reported as a target of myotonic dystrophy 1-associated splicing alterations. These suggest that changes in TACC2 expression are associated with variations in muscle strength in older people. The association was consistently observed in young and middle-aged subjects. Our findings would shed light on genetic components of lower limb muscle strength and indicate TACC2 as a potential therapeutic target for sarcopenia.

Abstract Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Abstract Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10−10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

Abstract Background Polygenic risk score (PRS) analysis is used to predict disease risk. Although PRS has been shown to have great potential in improving clinical care, PRS accuracy assessment has been mainly focused on European ancestry. This study aimed to develop an accurate genetic risk score for knee osteoarthritis (OA) using a multi-population PRS and leveraging a multi-trait PRS in the Japanese population. Methods We calculated PRS using PRS-CS-auto, derived from genome-wide association study (GWAS) summary statistics for knee OA in the Japanese population (same ancestry) and multi-population. We further identified risk factor traits for which PRS could predict knee OA and subsequently developed an integrated PRS based on multi-trait analysis of GWAS (MTAG), including genetically correlated risk traits. PRS performance was evaluated in participants of the Nagahama cohort study who underwent radiographic evaluation of the knees (n = 3,279). PRSs were incorporated into knee OA integrated risk models along with clinical risk factors. Results A total of 2,852 genotyped individuals were included in the PRS analysis. The PRS based on Japanese knee OA GWAS was not associated with knee OA (p = 0.228). In contrast, PRS based on multi-population knee OA GWAS showed a significant association with knee OA (p = 6.7 × 10−5, odds ratio (OR) per standard deviation = 1.19), whereas PRS based on MTAG of multi-population knee OA, along with risk factor traits such as body mass index GWAS, displayed an even stronger association with knee OA (p = 5.4 × 10−7, OR = 1.24). Incorporating this PRS into traditional risk factors improved the predictive ability of knee OA (area under the curve, 74.4% to 74.7%; p = 0.029). Conclusions This study showed that multi-trait PRS based on MTAG, combined with traditional risk factors, and using large sample size multi-population GWAS, significantly improved predictive accuracy for knee OA in the Japanese population, even when the sample size of GWAS of the same ancestry was small. To the best of our knowledge, this is the first study to show a statistically significant association between the PRS and knee OA in a non-European population. Trial registration No. C278.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E–11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

ABSTRACT Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty‐related exposures and AIS risk in Europeans and Asians. For our two‐sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist‐hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25‐hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome‐wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance‐weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy‐sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Abstract Background Skeletal dysplasia is a heterogeneous group of disorders. Spondyloepiphyseal dysplasias comprise one subgroup. Deficiency of carbohydrate sulfotransferase 3 has been reported in a small number of patients with recessively inherited spondyloepiphyseal dysplasia with joint dislocation, short stature and scoliosis. We report here molecular and clinical findings of affected individuals in three consanguineous Pakistani families. Affected individuals in all three families had a uniform phenotype including severe short stature, multiple dislocated joints, progressive scoliosis and facial dysmorphism. Methods Clinical evaluation was done for three unrelated families. Radiological survey of bones was completed for patients from two of the families. Whole exome sequencing index patients from each family was performed followed by Sanger sequencing for validation of segregation of identified variants in respective families. In-silico analysis for determining pathogenicity of identified variants and conservation was done. Results Whole-exome sequencing revealed biallelic variants c.590 T > C;p.(Leu197Pro), c.603C > A;p.(Tyr201Ter) and c.661C > T;p.(Arg221Cys) in CHST3 (NM_004273.5) in the three families with eight, five and two affected individuals, respectively. Contrary to previous reports, affected individuals in none of the families exhibited a hearing loss. Conclusion We describe genotypic and phenotypic findings of three unrelated families with spondyloepiphyseal dysplasia. Our study confirms phenotypic variability and adds to the genotypic spectrum of spondyloepiphyseal dysplasia.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.

IntroductionAdolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated.Material and methodsMendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese.ResultsSignificant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI.ConclusionsOur Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Osteoarthritis (OA), in which M1 macrophage polarization in the synovium exacerbates disease progression, is a major cause of cartilage degeneration and functional disabilities. Therapeutic strategies of OA designed to interfere with the polarization of macrophages have rarely been reported. Here, we report that SHP099, as an allosteric inhibitor of src-homology 2-containing protein tyrosine phosphatase 2 (SHP2), attenuated osteoarthritis progression by inhibiting M1 macrophage polarization. We demonstrated that M1 macrophage polarization was accompanied by the overexpression of SHP2 in the synovial tissues of OA patients and OA model mice. Compared to wild-type (WT) mice, myeloid lineage conditional Shp2 knockout (cKO) mice showed decreased M1 macrophage polarization and attenuated severity of synovitis, an elevated expression of cartilage phenotype protein collagen II (COL2), and a decreased expression of cartilage degradation markers collagen X (COL10) and matrix metalloproteinase 3 (MMP3) in OA cartilage. Further mechanistic analysis showed thatSHP099 inhibited lipopolysaccharide (LPS)-induced Toll-like receptor (TLR) signaling mediated by nuclear factor kappa B (NF-κB) and PI3K–AKT signaling. Moreover, intra-articular injection of SHP099 also significantly attenuated OA progression, including joint synovitis and cartilage damage. These results indicated that allosteric inhibition of SHP2 might be a promising therapeutic strategy for the treatment of OA.

Abstract Background Hypophosphatemic rickets (HR) is a rare genetic disorder associated with renal phosphate wasting and characterized by bone defects. Inactivating mutations in the phosphate regulating endopeptidase homolog X‑linked gene (PHEX) account for most cases of HR. The aim of this study was to identify causative variants in nine unrelated Chinese families associated with HR, and to determine potential pathogenicity of the identified variants. Methods Genomic DNA was isolated from the peripheral blood of HR patients and their healthy relatives, followed by next-generation sequencing and/or Sanger sequencing. In silico prediction combined with conservation analysis was performed to assess the effects of the variants, and 3D protein modeling was conducted to predict the functional effects on the encoded protein. Results All HR patients recruited in this study displayed bone deformities and tooth agenesis, as well as reduced serum phosphate levels and elevated urine phosphate levels. Nine PHEX variants were identified in eight families, including four novel variants (c.1661_1726del, c.980A > G, c.1078A > T, and c.1017_1051dup). Of the nine identified PHEX variants, five caused a truncated protein, two caused an altered amino acid, and the other two were the canonical splicing variants. Novel variants c.1336G > A and c.1364 T > C in SLC34A3 were also found in one family. Conservation analysis showed that all the amino acids corresponding to the missense variants were highly conserved. In silico analysis and 3D protein structure modeling confirmed the pathogenicity of these variants. Conclusions This study identified four novel variants in PHEX and two novel variants in SLC34A3 in a Chinese cohort with HR. Our findings highlight the dominant role of PHEX in HR, and expand the genotypic and phenotypic spectra of this disorder.

Summary: Background: Osteoarthritis (OA) is the most common degenerative joint disease primarily characterized by cartilage destruction. The aim of this study was to investigate the role, molecular characteristics and potential therapeutic target of chondrocyte ferroptosis in the pathogenesis of OA. Methods: The expression of ferroptotic hallmarks (iron and lipid peroxidation accumulation, glutathione deletion) were analyzed in paired intact and damaged cartilages from OA patients. Single cell RNA sequencing (scRNA-seq) analysis was performed on 17,638 chondrocytes to verify the presence, investigate the molecular signatures and unveil the potential therapeutic target of ferroptotic chondrocyte cluster in human OA cartilages. Destabilization of medial meniscus (DMM)-induced OA model and tert-butyl hydroperoxide (TBHP)-treated primary mouse chondrocytes and human cartilage explants were used to evaluate the protective effect of pharmacologically activated transient receptor potential vanilloid 1 (TRPV1). The downstream molecular mechanisms of TRPV1 was further investigated in glutathione peroxidase 4 (Gpx4) heterozygous genetic deletion mice (Gpx4+/−). Findings: The concentrations of iron and lipid peroxidation and the expression of ferroptotic drivers in the damaged areas of human OA cartilages were significantly higher than those in the intact cartilage. scRNA-seq analysis revealed a chondrocyte cluster characterized by preferentially expressed ferroptotic hallmarks and genes, namely ferroptotic chondrocyte cluster. Comprehensive gene set variation analysis revealed TRPV1 as an anti-ferroptotic target in human OA cartilage. Pharmacological activation of TRPV1 significantly abrogated cartilage degeneration by protecting chondrocytes from ferroptosis. Mechanistically, TRPV1 promoted the expression of GPX4, and its anti-ferroptotic role was largely mitigated in the OA model of Gpx4+/− mice. Interpretation: TRPV1 activation protects chondrocytes from ferroptosis and ameliorates OA progression by upregulating GPX4. Funding: National Key R&D Program of China (2018YFC1105904), Key Program of NSFC (81730067), National Science Foundation of China (81772335, 81941009, 81802196), Natural Science Foundation of Jiangsu Province, China (BK20180127), Jiangsu Provincial Key Medical Talent Foundation, Six Talent Peaks Project of Jiangsu Province (WSW-079).

Background: Synovial inflammation plays a major role in the pathogenesis of osteoarthritis (OA). This study investigated the effect of andrographolide (Andro) on synovial inflammation mediated by tumor necrosis factor-alpha receptor 2 (TNFR2) trafficking and its utility in attenuating OA progression. Methods: Knee joints were harvested from rats subjected to radial transection of the medial collateral ligament (MCLT) and medial meniscus (MMT) to examine the effect of Andro on synovial inflammation and OA progression. Quantitative real-time polymerase chain reaction was used to evaluate the expression of inflammatory factors in primary fibroblast-like synoviocytes (FLSs) after Andro treatment in vitro. The mechanism underlying Andro-mediated regulation of TNFR2 distribution and nuclear factor-κB (NF-κB) expression was verified using endosome maturation inhibitor hydroxychloroquine (HCQ) through flow cytometry, immunofluorescence, and western blot analysis. Results: Andro treatment was found to reduce synovial inflammation and OA progression in vivo. Furthermore, a decrease in pain hypersensitivity and dorsal horn neuron activation was observed after treatment. Andro also downregulated the expression of inflammatory mediators and TNFR2 in FLSs. TNFR2 is crucial for the activation of the NF-κB signaling pathway, and Andro-induced degradation of TNFR2 was associated with lysosomal function, which in turn, reduced the downstream phosphorylation of p65 in the NF-κB signaling pathway. Conclusions: Andro could suppress synovial inflammation via regulation of TNFR2 trafficking and degradation. This also suggests it could be a potential treatment for the prevention of synovial inflammation and OA progression. The translational potential of this article: This study provides strong evidence that Andro reduces NF-κB activation and inflammatory responses in OA FLSs via regulation of TNFR2 trafficking. The inhibition of TNFR2 and Andro could be a novel therapeutic approach for OA and pain management.

Abstract Osteoarthritis (OA) is the major course of joint deterioration, in which M1 macrophage-driven synovitis exacerbates the pathological process. However, precise therapies for M1 macrophage to decrease synovitis and attenuate OA progression have been scarcely proposed. Transient receptor potential vanilloid 1 (TRPV1) is a cation channel that has been implicated in pain perception and inflammation. In this study, we investigated the role of TRPV1 in the M1 macrophage polarization and pathogenesis of OA. We demonstrated that TRPV1 expression and M1 macrophage infiltration were simultaneously increased in both human and rat OA synovium. More than 90% of the infiltrated M1 macrophages expressed TRPV1. In the rat OA model, intra-articular injection of capsaicin (CPS), a specific TRPV1 agonist, significantly attenuated OA phenotypes, including joint swelling, synovitis, cartilage damage, and osteophyte formation. CPS treatment markedly reduced M1 macrophage infiltration in the synovium. Further mechanistic analyses showed that TRPV1-evoked Ca2+ influx promoted the phosphorylation of calcium/calmodulin-dependent protein kinase II (CaMKII) and facilitated the nuclear localization of nuclear factor-erythroid 2-related factor 2 (Nrf2), which ultimately resulted in the inhibition of M1 macrophage polarization. Taken together, our findings establish that TRPV1 attenuates the progression of OA by inhibiting M1 macrophage polarization in synovium via the Ca2+/CaMKII/Nrf2 signaling pathway. These results highlight the effect of targeting TRPV1 for the development of a promising therapeutic strategy for OA.

Sclerosing bone disorder (SBD) includes a broad spectrum of monogenic diseases characterised by increased bone density. Here, the authors describe a previously unknown SBD in four families caused by mutations in TMEM53 and demonstrate the role this protein plays in BMP signalling during bone formation.

Dysosteosclerosis (DOS) is a rare sclerosing bone dysplasia characterized by unique osteosclerosis of the long tubular bones and platyspondyly. DOS is inherited in an autosomal recessive manner and is genetically and clinically heterogeneous. To date, four individuals with DOS who have five different TNFRSF11A mutations have been reported. Based on their data, it is hypothesized that mutations producing aberrant mutant RANK proteins (missense or truncated or elongated) cause DOS, while null mutations lead to osteopetrosis, autosomal recessive 7 (OPTB7). Herein, we present the fifth case of TNFRSF11A-associated DOS with a novel homozygous frame-shift mutation (c.19_31del; p.[Arg7CysfsTer172]). The mutation is predicted to cause nonsense mutation-mediated mRNA decay (NMD) in all RANK isoform transcripts, resulting in totally null allele. Our findings suggest genotype-phenotype relationship in TNFRSF11A-associated OPTB7 and DOS remains unclear, and that the deficiency of TNFRSF11A functions might cause DOS, rather than osteopetrosis. More data are necessary to understand the phenotypic spectrum caused by TNFRSF11A mutations.

Desbuquois dysplasia (DD) type 1 is a rare skeletal dysplasia characterized by a short stature, round face, progressive scoliosis, and joint laxity. The causative gene has been identified as calcium‐activated nucleotidase 1 (CANT1), which encodes a nucleotidase that preferentially hydrolyzes UDP to UMP and phosphate. In this study, we generated Cant1 KO mice using CRISPR/Cas9‐mediated genome editing. All F0 mice possessing frameshift deletions at both Cant1 alleles exhibited a dwarf phenotype. Germline transmission of the edited allele was confirmed in an F0 heterozygous mouse, and KO mice were generated by crossing of the heterozygous breeding pairs. Cant1 KO mice exhibited skeletal defects, including short stature, thoracic kyphosis, and delta phalanx, all of which are observed in DD type 1 patients. The glycosaminoglycan (GAG) content and extracellular matrix space were reduced in the growth plate cartilage of mutants, and proliferating chondrocytes lost their typical flat shape and became round. Chondrocyte differentiation, especially terminal differentiation to hypertrophic chondrocytes, was impaired in Cant1 KO mice. These findings indicate that CANT1 is involved in the synthesis of GAG and regulation of chondrocyte differentiation in the cartilage and contribute to a better understanding of the pathogenesis of DD type 1.

Mesenchymal stem cells (MSCs) are well known for their multi-directional differentiation potential and are widely applied in cartilage and bone disease. Synovial mesenchymal stem cells (SMSCs) exhibit a high proliferation rate, low immunogenicity, and greater chondrogenic differentiation potential. Microtubule (MT) plays a key role in various cellular processes. Perturbation of MT stability and their associated proteins is an underlying cause for diseases. Little is known about the role of MT stabilization in the differentiation and homeostasis of SMSCs. In this study, we demonstrated that MT stabilization via docetaxel treatment had a significant effect on enhancing the chondrogenic differentiation of SMSCs. MT stabilization inhibited the expression of Yes-associated proteins (YAP) and the formation of primary cilia in SMSCs to drive chondrogenesis. This finding suggested that MT stabilization might be a promising therapeutic target of cartilage regeneration.

Knee osteoarthritis (KOA) is the most common form of joint degeneration with increasing prevalence and incidence in recent decades. KOA is a molecular disorder characterized by the interplay of numerous molecules, a considerable number of which can be detected in body fluids, including synovial fluid, urine, and blood. However, the current diagnosis and treatment of KOA mainly rely on clinical and imaging manifestations, neglecting its molecular pathophysiology. The mismatch between participants’ molecular characteristics and drug therapeutic mechanisms might explain the failure of some disease-modifying drugs in clinical trials. Hence, according to the temporal alteration of representative molecules, we propose a novel molecular classification of KOA divided into pre-KOA, early KOA, progressive KOA, and end-stage KOA. Then, progressive KOA is furtherly divided into four subtypes as cartilage degradation-driven, bone remodeling-driven, inflammation-driven, and pain-driven subtype, based on the major pathophysiology in patient clusters. Multiple clinical findings of representatively investigated molecules in recent years will be reviewed and categorized. This molecular classification allows for the prediction of high-risk KOA individuals, the diagnosis of early KOA patients, the assessment of therapeutic efficacy, and in particular, the selection of homogenous patients who may benefit most from the appropriate therapeutic agents.

Abstract Genetic white matter disorders have heterogeneous etiologies and overlapping clinical presentations. We performed a study of the diagnostic efficacy of genome sequencing in 41 unsolved cases with prior exome sequencing, resolving an additional 14 from an historical cohort (n = 191). Reanalysis in the context of novel disease‐associated genes and improved variant curation and annotation resolved 64% of cases. The remaining diagnoses were directly attributable to genome sequencing, including cases with small and large copy number variants (CNVs) and variants in deep intronic and technically difficult regions. Genome sequencing, in combination with other methodologies, achieved a diagnostic yield of 85% in this retrospective cohort.

Thousands of genetic loci are known to associate with human height, but these are mainly based on studies in European ancestry populations. Here, Akiyama et al. construct a genotype reference panel for the Japanese population followed by GWAS and report 573 height associated variants in 191,787 Japanese.

Adolescent idiopathic scoliosis (AIS) is a common pediatric disease leading to spinal deformities. Here, the authors report GWAS followed by genome-wide meta-analysis in up to 79,211 Japanese individuals, identifying 20 genetic loci for AIS, 14 of which were previously unreported, and perform in vitro validation for rs1978060.

Abstract Background Patients with ankylosing spines are susceptible to developing spinal fractures even with minor trauma and can develop early or late neurological injuries. These fractures require early and aggressive surgical management to enable spinal stability and/or neural decompression. Being highly unstable by nature, they require relatively long segment instrumentation and fusion, which can increase paravertebral soft tissue damage and perioperative bleeding. The purpose of this report is to describe a rare case of traumatic double fractures at the cervico-thoracic and thoraco-lumbar transition zones in ankylosing spine with spondylo-epiphyseal dysplasia (SED) of unknown cause, which were successfully treated with a combined open and percutaneous spinal fusion procedure. Case presentation A 46-year-old woman who was diagnosed with non-contiguous fractures in cervico-thoracic and thoraco-lumbar junction zones among multiple injuries sustained in a traffic accident was treated with hybrid techniques for posterior instrumentation with an open approach using a computed tomography (CT)-based navigation system and percutaneous pedicle-screwing method. She regained mobility to pre-admission levels and started walking on crutches 3 months postoperatively. Genetic testing for the cause of SED revealed no mutation in the COL2A1 or TRPVR4 genes. The union of fractured spine was confirmed on CT scan 1 year postoperatively. Conclusion This is the first report of double spinal fractures in an ankylosing spine with genetically undetermined spondyloepiphyseal dysplasia. A long-segment posterior instrumentation procedure incorporating the invasive treatment of spinal fractures in ankylosing spondylitis or diffuse idiopathic hyperostosis was effective.

Tumor-induced rickets/osteomalacia (TIO) is a rare paraneoplastic syndrome caused by tumors that ectopically express fibroblast growth factor 23 (FGF23). FGF23 is a bone-derived hormone that regulates serum phosphate concentrations. Patients with TIO develop hypophosphatemic rickets/osteomalacia due to FGF23 excess and suffer from symptoms such as leg deformities, bone pain, skeletal muscle myopathy, and multiple fractures/pseudofractures. Usually, successful surgical removal of the causative tumors normalizes serum FGF23 and phosphate concentrations in patients with TIO. Most FGF23-producing tumors associated with TIO are histologically called phosphaturic mesenchymal tumor, mixed connective tissue variant (PMTMCT). The precise mechanism by which these tumors ectopically overproduce FGF23 outside of bone is yet to be clarified. Therefore, we performed an RNA sequencing analysis of a PMTMCT that was found in the left parotid gland of a patient with TIO. Among the upregulated genes, we focused on Klotho, the protein product of which is a single pass transmembrane protein that works along with an FGF receptor 1c as a receptor complex for FGF23. Subsequent histological analysis confirmed the ectopic expression of Klotho in other PMTMCTs. From these results, we assume that the ectopic expression of Klotho in PTMMCTs enables a positive feedback loop in FGF23 production via the activation of FGF receptor 1c and exacerbates disease manifestations in TIO. Keywords: Tumor-induced osteomalacia (TIO), Klotho, Fibroblast growth factor 23 (FGF23), FGF receptor (FGFR), Hypophosphatemia

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Lumbar disc degeneration (LDD) is age-related break-down in the fibrocartilaginous joints between lumbar vertebrae. It is a major cause of low back pain and is conventionally assessed by magnetic resonance imaging (MRI). Like most other complex traits, LDD is likely polygenic and influenced by both genetic and environmental factors. However, genome-wide association studies (GWASs) of LDD have uncovered few susceptibility loci due to the limited sample size. Previous epidemiology studies of LDD also reported multiple heritable risk factors, including height, body mass index (BMI), bone mineral density (BMD), lipid levels, etc. Genetics can help elucidate causality between traits and suggest loci with pleiotropic effects. One such approach is polygenic score (PGS) which summarizes the effect of multiple variants by the summation of alleles weighted by estimated effects from GWAS. To investigate genetic overlaps of LDD and related heritable risk factors, we calculated the PGS of height, BMI, BMD and lipid levels in a Chinese population-based cohort with spine MRI examination and a Japanese case-control cohort of lumbar disc herniation (LDH) requiring surgery. Because most large-scale GWASs were done in European populations, PGS of corresponding traits were created using weights from European GWASs. We calibrated their prediction performance in independent Chinese samples, then tested associations with MRI-derived LDD scores and LDH affection status. The PGS of height, BMI, BMD and lipid levels were strongly associated with respective phenotypes in Chinese, but phenotype variances explained were lower than in Europeans which would reduce the power to detect genetic overlaps. Despite of this, the PGS of BMI and lumbar spine BMD were significantly associated with LDD scores; and the PGS of height was associated with the increased the liability of LDH. Furthermore, linkage disequilibrium score regression suggested that, osteoarthritis, another degenerative disorder that shares common features with LDD, also showed genetic correlations with height, BMI and BMD. The findings suggest a common key contribution of biomechanical stress to the pathogenesis of LDD and will direct the future search for pleiotropic genes.

Abstract The zinc transporter protein ZIP13 plays critical roles in bone, tooth, and connective tissue development, and its dysfunction is responsible for the spondylocheirodysplastic form of Ehlers‐Danlos syndrome (SCD‐EDS, OMIM 612350). Here, we report the molecular pathogenic mechanism of SCD‐EDS caused by two different mutant ZIP13 proteins found in human patients: ZIP13G64D, in which Gly at amino acid position 64 is replaced by Asp, and ZIP13ΔFLA, which contains a deletion of Phe‐Leu‐Ala. We demonstrated that both the ZIP13G64D and ZIP13ΔFLA protein levels are decreased by degradation via the valosin‐containing protein (VCP)‐linked ubiquitin proteasome pathway. The inhibition of degradation pathways rescued the protein expression levels, resulting in improved intracellular Zn homeostasis. Our findings uncover the pathogenic mechanisms elicited by mutant ZIP13 proteins. Further elucidation of these degradation processes may lead to novel therapeutic targets for SCD‐EDS.