Search

Study Design: A sex- and age-matched case-control study and a cross-sectional study. Objective: In our previous study, patients with early-onset (Methods: In Study 1, the simple brief-type self-administered diet history questionnaire (BDHQ) was used to compare nutrient intake levels of early-onset OPLL patients (n = 13) with those of sex- and age-matched non-OPLL controls (n = 39) or with those of common OPLL (onset age ≥ 50 years, n = 62). In Study 2, serological validation was conducted for thoracic OPLL patients (n = 77) and non-OPLL controls (n = 101) in a nationwide multicenter study in Japan. Results: The BDHQ showed that the early-onset OPLL patients had significantly lower intakes of vitamins A and B6 than non-OPLL controls. These results were validated by lower serum vitamins A and B6 levels in the early-onset thoracic OPLL patients. The severity of OPLL negatively correlated with serum vitamin A levels in male early-onset OPLL patients. The multiple regression analysis revealed that the severity of thoracic OPLL had an association with onset age and serum vitamin A level. Conclusions: Vitamin A deficiency resulting from unbalanced dietary habits is associated with exacerbation of male early-onset OPLL.

No genome-wide association studies (GWAS) were reported for colorectal polyps and the overlap in polygenic backgrounds conferring risk of colorectal cancer and polyps remains unclear. We performed GWAS on subjects with colorectal polyps using the BioBank Japan data with 4447 cases and 157,226 controls. We evaluated genetic correlations between colorectal polyps and cancer, and effects on colorectal polyps of single nucleotide polymorphisms (SNPs) known to be associated with colorectal cancer. We identified CUX2, a known genetic locus to colorectal cancer, as a susceptibility locus to colorectal polyps (p value = 1.1 × 10−15). Subsequent fine-mapping analysis indicated that rs11065828 in CUX2 is the causal variant for colorectal polyps. We found that known colorectal cancer-susceptible SNPs were also associated with colorectal polyps. The genetic correlation between colorectal cancer and polyps is very high (r = 0.98 and p value = 0.0006). We additionally identified 14 significant loci of colorectal polyps and three significant loci of colorectal cancer by applying the multi-trait analysis of GWAS of colorectal cancer and colorectal polyps. We showed very similar germline polygenic features, which gives us the additional insight into potential cancers at polygenic levels for patients with polyps who are followed up at outpatients’ clinic; thus, close observation and polypectomy is critical to prevent colorectal cancers.

Osteoarthritis (OA) is the most prevalent age-related disorder due to cartilage degradation. Previous studies have identified aberrant chondrocyte homeostasis under extracellular stress as a key pathological mechanism behind cartilage degradation in OA. TDP-43, a DNA/RNA-binding protein has been demonstrated to participate in processing many extracellular stress responses; however, understanding of the role of TDP-43 in OA is limited. This study aims to investigate the role of TDP-43 in chondrocyte homeostasis and cartilage degradation in OA.The role of TDP-43 during degradation of cartilage is examined by experimental posttraumatic OA animal models and human cartilage specimens. Cartilage degradation is assessed by histological analysis, qPCR, and Western blot. The molecular mechanisms are investigated in vitro using human primary chondrocytes.TDP-43 decreases significantly in degenerated cartilage. TDP-43 concentration is positively correlated with IL-1β concentration in synovial fluid derived from OA patients (Pearson r = 0.95, CI (95%) [0.80, 0.99], P 0.0001). Intra-articular injection of recombinant TDP-43 significantly alleviates cartilage degradation and subchondral bone remodeling in vivo. In vitro mechanistic analyses show that TDP-43 maintains chondrocyte homeostasis under oxidative stress through regulating stress granule dynamics via G3BP1.The present study indicates that TDP-43 maintains chondrocyte homeostasis under oxidative stress and alleviates cartilage degeneration in osteoarthritis, identifying TDP-43 as a potential target for the diagnosis and treatment of knee OA.

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).

Colony-stimulating factor 1 receptor (CSF1R) plays key roles in the development and function of the cells in the monocyte/macrophage lineage, including microglia and osteoclasts. It is well known that mono-allelic mutations of CSF1R cause hereditary diffuse leukoencephalopathy with spheroids (HDLS, OMIM # 221820), an adult-onset progressive neurodegenerative disorder. Recently, a more severe phenotypic spectrum has been identified in individuals with bi-allelic mutations of CSF1R. In addition to leukoencephalopathy of earlier onset than HDLS, the new disease shows brain malformations and skeletal dysplasia compatible with dysosteosclerosis (DOS), thus named "brain abnormalities, neurodegeneration, and dysosteosclerosis" (BANDDOS, OMIM # 618476). In addition, some individuals with bi-allelic missense mutations of CSF1R have been found to present with incomplete BANDDOS where skeletal dysplasia is absent. In this review, we summarize the monogenic disorders caused by mutations in CSF1R and their mutational spectra, and propose a dose-dependent model to explain the complex genotype-phenotype association.

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.

Summary Chondrodysplasias are hereditary diseases caused by mutations in the components of growth cartilage. Although the unfolded protein response (UPR) has been identified as a key disease mechanism in mouse models, no suitable in vitro system has been reported to analyze the pathology in humans. Here, we developed a three-dimensional culture protocol to differentiate hypertrophic chondrocytes from induced pluripotent stem cells (iPSCs) and examine the phenotype caused by MATN3 and COL10A1 mutations. Intracellular MATN3 or COL10 retention resulted in increased ER stress markers and ER size in most mutants, but activation of the UPR was dependent on the mutation. Transcriptome analysis confirmed a UPR with wide-ranging changes in bone homeostasis, extracellular matrix composition, and lipid metabolism in the MATN3 T120M mutant, which further showed altered cellular morphology in iPSC-derived growth-plate-like structures in vivo. We then applied our in vitro model to drug testing, whereby trimethylamine N-oxide led to a reduction of ER stress and intracellular MATN3., Highlights • A new induction method enables hypertrophic chondrocyte differentiation from iPSCs • Chondrodysplasia mutants show intracellular MATN3 or COL10 accumulation • Some, but not all, mutants have ER stress with an unfolded protein response • This induction system is applicable to transcriptomic analysis and drug development, Chondrodysplasias are highly heterogeneous genetic cartilage diseases. Here, Toguchida and colleagues establish a new induction system to differentiate hypertrophic chondrocytes from iPSCs and analyze these disorders in vitro. They found that different chondrodysplasia mutants, despite mutations being in the same gene, showed varying phenotypes and transcriptomic changes. This system provides an initial platform for further investigation and drug development.

ObjectiveSystemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations.MethodsWe newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations.ResultsWe identified 113 genetic regions including 46 novel loci at genome-wide significance (p−8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=−0.242) and non-albumin protein (rg=0.238).ConclusionThis study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Dysosteosclerosis (DOS) is a distinct form of sclerosing bone disease characterized by platyspondyly and progressive osteosclerosis. DOS is genetically heterogeneous. Three causal genes, SLC29A3, CSF1R, and TNFRSF11A are reported. TNFRSF11A-associated DOS has been identified in two patients; however, TNFRSF11A is also a causal gene for osteopetrosis, autosomal recessive 7 (OP-AR7). Whole-exome sequencing in a patient with sclerosing bone disease identified novel compound heterozygous variants (c.414_427 + 7del, c.1664del) in TNFRSF11A. We examined the impact of the two variants on five splicing isoforms of TNFRSF11A by RT-PCR. We found that c.1664del resulted in elongated proteins (p.S555Cfs*121, etc.), while c.414_427 + 7del generated two aberrant splicing products (p.A139Wfs*19 and p.E132Dfs*19) that lead to nonsense mediated mRNA decay (NMD). In the previous two cases of TNFRSF11A-associated DOS, their mutations produced truncated TNFRSF11A protein isoforms. The mutations in all three cases thus contrast with TNFRSF11A mutations reported in OP-AR7, which does not generated truncated or elongated TNFRSF11A proteins. Thus, we identified the third case of TNFRSF11A-associated DOS and reinforced the genotype-phenotype correlation that aberrant protein-producing TNFRSF11A mutations cause DOS.

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., Mutations in the gene encoding for CANT1 (calcium‐activated nucleotidase 1) are responsible for Desbuquois dysplasia (DD) type 1. Cant1 KO mice exhibit skeletal defects that recapitulate DD type 1. Reduced glycosaminoglycan (GAG) content and impaired chondrocyte differentiation are observed in the growth plate cartilage of these mutants. These findings indicate that CANT1 is involved in the synthesis of GAG and regulation of chondrocyte differentiation.

Osteonecrosis of the femoral head (ONFH) involves necrosis of bone and bone marrow of the femoral head caused by ischemia with unknown etiology. Previous genetic studies on ONFH failed to produce consistent results, presumably because ONFH has various causes with different genetic backgrounds and the underlying diseases confounded the associations. Steroid-associated ONFH (S-ONFH) accounts for one-half of all ONFH, and systemic lupus erythematosus (SLE) is a representative disease underlying S-ONFH. We performed a genome-wide association study (GWAS) to identify genetic risk factors for S-ONFH in patients with SLE. We conducted a two-staged GWAS on 636 SLE patients with S-ONFH and 95 588 non-SLE controls. Among the novel loci identified, we determined S-ONFH-specific loci by comparing allele frequencies between SLE patients without S-ONFH and non-SLE controls. We also used Korean datasets comprising 148 S-ONFH cases and 37 015 controls to assess overall significance. We evaluated the functional annotations of significant variants by in silico analyses. The Japanese GWAS identified 4 significant loci together with 12 known SLE susceptibility loci. The four significant variants showed comparable effect sizes on S-ONFH compared with SLE controls and non-SLE controls. Three of the four loci, MIR4293/MIR1265 [odds ratio (OR) = 1.99, P-value = 1.1 × 10−9)], TRIM49/NAALAD2 (OR = 1.65, P-value = 4.8 × 10−8) and MYO16 (OR = 3.91, P-value = 4.9 × 10−10), showed significant associations in the meta-analysis with Korean datasets. Bioinformatics analyses identified MIR4293, NAALAD2 and MYO16 as candidate causal genes. MIR4293 regulates a PPARG-related adipogenesis pathway relevant to S-ONFH. We identified three novel susceptibility loci for S-ONFH in SLE.

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.

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 HCO3 - . Bi-allelic 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 bi-allelic 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 importance of SLC4A2 in human osteoclast differentiation. This article is protected by copyright. All rights reserved.

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.

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.

Congenital scoliosis (CS) is a form of scoliosis caused by congenital vertebral malformations. Genetic predisposition has been demonstrated in CS. We previously reported that TBX6 loss-of-function causes CS in a compound heterozygous model; however, this model can explain only 10% of CS. Many monogenic and polygenic CS genes remain to be elucidated. In this study, we analyzed exome sequencing (ES) data of 615 Chinese CS from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) project. Cosegregation studies for 103 familial CS identified a novel heterozygous nonsense variant, c.2649G>A (p.Trp883Ter) in FBN1. The association between FBN1 and CS was then analyzed by extracting FBN1 variants from ES data of 574 sporadic CS and 828 controls; 30 novel variants were identified and prioritized for further analyses. A mutational burden test showed that the deleterious FBN1 variants were significantly enriched in CS subjects (OR = 3.9, P = 0.03 by Fisher’s exact test). One missense variant, c.2613A>C (p.Leu871Phe) was recurrent in two unrelated CS subjects, and in vitro functional experiments for the variant suggest that FBN1 may contribute to CS by upregulating the transforming growth factor beta (TGF-β) signaling. Our study expanded the phenotypic spectrum of FBN1, and provided nove insights into the genetic etiology of CS.

Cysteine-rich transmembrane bone morphogenetic protein regulator 1 (CRIM1) is a type I transmembrane protein involved in the organogenesis of many tissues via its interactions with growth factors including BMP, TGF-β, and VEGF. In this study, we used whole-exome sequencing and linkage analysis to identify a novel Crim1 mutant allele generated by ENU mutagenesis in mice. This allele is a missense mutation that causes a cysteine-to-serine substitution at position 140, and is referred to as Crim1C140S. In addition to the previously reported phenotypes in Crim1 mutants, Crim1C140S homozygous mice exhibited several novel phenotypes, including dwarfism, enlarged seminal vesicles, and rectal prolapse. In vitro analyses showed that Crim1C140S mutation affected the formation of CRIM1 complexes and decreased the amount of the overexpressed CRIM1 proteins in the cell culture supernatants. Cys140 is located in the internal region 1 (IR1) of the N-terminal extracellular region of CRIM1 and resides outside any identified functional domains. Inference of the domain architecture suggested that the Crim1C140S mutation disturbs an intramolecular disulfide bond in IR1, leading to the protein instability and the functional defects of CRIM1. Crim1C140S highlights the functional importance of the IR1, and Crim1C140S mice should serve as a valuable model for investigating the functions of CRIM1 that are unidentified as yet.

To clarify the genetic mechanisms underlying intervertebral disc degeneration (IDD), we examined the associations between single-nucleotide polymorphisms (SNPs) and indicated as coefficient of interaction term (IDD) in a general population in Japan.This was a cross-sectional study. In 1,605 participants, C2-3 to L5/S1 in the total spine magnetic resonance imaging (MRI) were evaluated using the Pfirrmann's scoring system. Disc scores of 4 and 5 were defined as IDD. Eight SNPs in eight genes associated with IDD were examined at each disc level, considering the non-genetic risk factors of age, sex, and body mass index (BMI).The highest odds ratio was found for rs9406328 in the THBS2 gene at disc level T12-L1 (OR 1.27, 95%CI 1.05 to 1.53), and this association was strengthened after adjustment for age using logistic regression (OR 1.37, 95%CI 1.12 to 1.67). Among participants aged50 years and 50-59, the average IDD score in those with 2 risk alleles of rs9406328 was markedly higher than in those with 0 or 1 risk allele, and the difference is much wider than the elderly participants. It indicates the genetic effect of rs9406328 is stronger in the younger age groups. Finally, multiple linear regression analyses of the association between rs9406328 and IDD, adjusted for age, sex, and BMI at each disc level, showed a statistical interaction between age and the number of risk alleles at C7-T1, T3-4 and T4-T5 as well as T12-L1.CONCLUSION: The association between rs9406328 in THBS2 and IDD was replicated. The contributions of genetic and environmental factors to IDD differed by disc level.

To characterize clinically measurable endophenotypes, implicating the TBX6 compound inheritance model. Patients with congenital scoliosis (CS) from China(N = 345, cohort 1), Japan (N = 142, cohort 2), and the United States (N = 10, cohort 3) were studied. Clinically measurable endophenotypes were compared according to the TBX6 genotypes. A mouse model for Tbx6 compound inheritance (N = 52) was investigated by micro computed tomography (micro-CT). A clinical diagnostic algorithm (TACScore) was developed to assist in clinical recognition of TBX6-associated CS (TACS). In cohort 1, TACS patients (N = 33) were significantly younger at onset than the remaining CS patients (P = 0.02), presented with one or more hemivertebrae/butterfly vertebrae (P = 4.9 × 10‒8), and exhibited vertebral malformations involving the lower part of the spine (T8–S5, P = 4.4 × 10‒3); observations were confirmed in two replication cohorts. Simple rib anomalies were prevalent in TACS patients (P = 3.1 × 10‒7), while intraspinal anomalies were uncommon (P = 7.0 × 10‒7). A clinically usable TACScore was developed with an area under the curve (AUC) of 0.9 (P = 1.6 × 10‒15). A Tbx6-/mh (mild-hypomorphic) mouse model supported that a gene dosage effect underlies the TACS phenotype. TACS is a clinically distinguishable entity with consistent clinically measurable endophenotypes. The type and distribution of vertebral column abnormalities in TBX6/Tbx6 compound inheritance implicate subtle perturbations in gene dosage as a cause of spine developmental birth defects responsible for about 10% of CS.

BackgroundCongenital scoliosis (CS) is a common vertebral malformation. Spondylocostal dysostosis (SCD) is a rare skeletal dysplasia characterised by multiple vertebral malformations and rib anomalies. In a previous study, a compound heterozygosity for a null mutation and a risk haplotype composed by three single-nucleotide polymorphisms inTBX6have been reported as a disease-causing model of CS. Another study identified bi-allelic missense variants in a SCD patient. The purpose of our study is to identifyTBX6variants in CS and SCD and examine their pathogenicity.MethodsWe recruited 200 patients with CS or SCD and investigatedTBX6variants. We evaluated the pathogenicity of the variants by in silico prediction and in vitro experiments.ResultsWe identified five 16p11.2 deletions, one splice-site variant and five missense variants in 10 patients. In vitro functional assays for missense variants identified in the previous and present studies demonstrated that most of the variants caused abnormal localisation of TBX6 proteins. We confirmed mislocalisation of TBX6 proteins in presomitic mesoderm cells induced from SCD patient-derived iPS cells. In induced cells, we found decreased mRNA expressions ofTBX6and its downstream genes were involved in somite formation. All CS patients with missense variants had the risk haplotype in the opposite allele, while a SCD patient with bi-allelic missense variants did not have the haplotype.ConclusionsOur study suggests that bi-allelic loss of function variants ofTBX6cause a spectrum of phenotypes including CS and SCD, depending on the severity of the loss ofTBX6function.

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.

Summary Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation., Graphical abstract, Highlights • A multicohort study identifies 52 previously unknown osteoarthritis genetic risk variants • Similarities and differences in osteoarthritis genetic risk depend on joint sites • Osteoarthritis genetic components are associated with pain-related phenotypes • High-confidence effector genes highlight potential targets for drug intervention, A multicohort genome-wide association meta-analysis of osteoarthritis highlights the impact of joint site types on the features of genetic risk variants and the link between osteoarthritis genetics and pain-related phenotypes, pointing toward potential targets for therapeutic intervention.

Dysosteosclerosis (DOS) is a rare sclerosing bone dysplasia characterized by osteosclerosis and platyspondyly. DOS is genetically heterogeneous and causally associated with mutations in three genes, SLC29A3, CSF1R, and TNFRSF11A. TNFRSF11A has been known as the causal gene for osteopetrosis, autosomal recessive 7, and is recently reported to cause DOS in three cases, which show a complex genotype-phenotype relationship. The phenotypic spectrum of TNFRSF11A-associated sclerosing bone dysplasia remains unclear and needs to be characterized further in more cases with molecular genetic diagnosis. Here, we report another TNFRSF11A-associated DOS case with a homozygous missense mutation (p.R129C). The mutation effect is different from the previous three cases, in which truncated or elongated RANK proteins were generated in isoform specific manner, thus enriching our understanding of the genotype-phenotype association in TNFRSF11A-associated sclerosing bone dysplasia. Besides DOS, our case presented with intracranial extramedullary hematopoiesis, which is an extremely rare condition and has not been identified in any other sclerosing bone dysplasias with molecular genetic diagnosis. Our findings provide the fourth case of TNFRSF11A-associated DOS and further expand its phenotypic spectrum.

Amyotrophic lateral sclerosis (ALS) is a devastating progressive motor neuron disease that affects people of all ethnicities. Approximately 90% of ALS cases are sporadic and thought to have multifactorial pathogenesis. To understand the genetics of sporadic ALS, we conducted a genome-wide association study using 1,173 sporadic ALS cases and 8,925 controls in a Japanese population. A combined meta-analysis of our Japanese cohort with individuals of European ancestry revealed a significant association at the ACSL5 locus (top SNP p = 2.97 × 10−8). We validated the association with ACSL5 in a replication study with a Chinese population and an independent Japanese population (1941 ALS cases, 3821 controls; top SNP p = 1.82 × 10−4). In the combined meta-analysis, the intronic ACSL5 SNP rs3736947 showed the strongest association (p = 7.81 × 10−11). Using a gene-based analysis of the full multi-ethnic dataset, we uncovered additional genes significantly associated with ALS: ERGIC1, RAPGEF5, FNBP1, and ATXN3. These results advance our understanding of the genetic basis of sporadic ALS., Gen Sobue, Masao Nagasaki and colleagues report a genome-wide association study for amyotrophic lateral sclerosis (ALS) in a large, multi-ethnic cohort comprising Japanese, Chinese, and European ancestry populations. They find a significant association to variants within the ACSL5 gene and identify novel associations with 4 additional genes using a gene-based approach.

Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci1-8. Nevertheless, these loci only partially explain SLE heritability and provide limited biological insight. We report the largest study of SLE in East Asians (13,377 cases and 194,993 controls), identifying 233 association signals within 113 (46 novel) genetic loci. We detect six new lead missense variants and prioritize ten most likely putative causal variants, one of which we demonstrate exhibits allele-specific regulatory effect on ACAP1 in vitro. We suggest 677 effector genes with potential for drug repurposing, and provide evidence that two distinct association signals at a single locus act on different genes (NCF2 and SMG7). We demonstrate that SLE-risk variants overlap with cell-specific active regulatory elements, notably EBNA2-mediated super-enhancers in Epstein-Barr Virus-transformed B cells, and implicate the role for non-immune cells in SLE biology. These findings shed light on genetic and biological understandings of SLE.

The RANKL/OPG/RANK signalling pathway is a major regulatory system for osteoclast formation and activity. Mutations in TNFSF11, TNFRSF11B and TNFRSF11A cause defects in bone metabolism and development, thereby leading to skeletal disorders with changes in bone density and/or morphology. To date, nine kinds of monogenic skeletal diseases have been found to be causally associated with TNFSF11, TNFRSF11B and TNFRSF11A mutations. These diseases can be divided into two types according to the mutation effects and the resultant pathogenesis. One is caused by the mutations inducing constitutional RANK activation or OPG deficiency, which increase osteoclastogenesis and accelerate bone turnover, resulting in juvenile Paget's disease 2, Paget disease of bone 2, familial expansile osteolysis, expansile skeletal hyperphosphatasia, panostotic expansile bone disease, and Paget disease of bone 5. The other is caused by the de-activating mutations in TNFRSF11A or TNFSF11, which decrease osteoclastogenesis and elevate bone density, resulting in osteopetrosis, autosomal recessive 2 and 7, and dysosteosclerosis. Here we reviewed the current knowledge about these genetic disorders with paying particular attention to the updating genotype-phenotype association in the TNFRSF11A-caused diseases.

SummaryChondrodysplasias are hereditary diseases caused by mutations in the components of growth cartilage. Although the unfolded protein response (UPR) has been identified as a key disease mechanism in mouse models, no suitablein vitrosystem has been reported to analyze the pathology in humans. Here, utilizing human chondrodysplasia-specific iPSCs, we examined the UPR caused by mutations inMATN3orCOL10A1. In growth plate-like structures formed from iPSC-derived sclerotomein vivo, the hypertrophic zone was disrupted, and induced hypertrophic chondrocytesin vitroshowed varying levels of ER stress depending on the mutation. Autophagy inducers and chemical chaperones succeeded in reducing ER stress only in some mutants, while transcriptome analysis revealed many mutation-specific changes in genes involved in apoptosis, metabolism, and protein trafficking. In this way, our system has allowed the precise evaluation of the UPR caused by each mutation, opening up new avenues for treatment of individual chondrodysplasia patients.

Background Congenital scoliosis (CS) is defined as a lateral curvature of the spine due to the vertebral malformations and has an incidence of 0.5–1/1,000 births. We previously examined TBX6 in Japanese CS patients and revealed that approximately 10% of CS was caused by TBX6 mutations. However, the genetic cause of remaining CS is unknown. Methods We recruited 78 CS patients without TBX6 mutations and major comorbidities, and investigated the genes previously reported to be associated with CS and congenital vertebral malformations by whole‐exome sequencing. Results We identified the compound heterozygous missense variants in LFNG in one patient. No likely disease‐causing variants were identified in other patients, however. LFNG encodes a GlcNAc‐transferase. The LFNG variants showed loss of their enzyme function. Conclusions A LFNG mutation is reported in a case of spondylocostal dysostosis (SCD), a skeletal dysplasia with severe malformations of vertebra and rib. The CS patient with LFNG mutations had multiple vertebral malformations including hemivertebrae, butterfly vertebrae, and block vertebrae, and rib malformations. LFNG mutations may cause a spectrum of phenotypes including CS and SCD. The current list of known disease genes could explain only a small fraction of genetic cause of CS.

Background Spinal extradural arachnoid cyst (SEDAC) is a cystic lesion that protrudes into the epidural space from a small dural defect. Early diagnosis of SEDAC is important because its expansion causes neurological damage. Two types of SEDAC, syndromic and sporadic, are present. Syndromic SEDAC is inherited as a part of lymphedema-distichiasis syndrome caused by mutations in the FOXC2 gene; however, it is often mistaken as sporadic because of low penetrance. It is not reasonable to conduct a genetic testing for all SEDAC patients and their family members. The aim of this study is to establish an effective screening method to distinguish syndromic SEDAC from sporadic SEDAC. Methods We performed a retrospective review of medical records and imaging studies of 29 subjects who were diagnosed with SEDAC. Clinical features, family history and magnetic resonance imaging (MRI) were analyzed. Mutations in FOXC2 were examined by Sanger-sequencing of the entire coding region of the genes. SEDAC having a mutation in FOXC2 gene was defined with syndromic SEDAC. Results Eleven subjects had a heterozygous mutation in FOXC2 . They were all familial and hence syndromic SEDAC. Only one proband had known family history of SEDAC at diagnosis. MRI findings and physical examinations, especially eye and leg examinations, were quite useful to screen syndromic SEDAC. Physical examination often showed accompanying lymphedema and distichiasis in syndromic SEDAC. Syndromic SEDAC tended to have multiple cysts out of the thoracolumbar area. Conclusions We established an effective screening method based on physical examinations and MRI findings.

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.

Ossification of the posterior longitudinal ligament (OPLL) of the spine is a common pathological condition that causes intractable myelopathy and radiculopathy, mainly the result of an endochondral ossification-like process. Our previous genome-wide association study identified six susceptibility loci for OPLL, including the cell division cycle 5-like (CDC5L) gene region. Here, we found CDC5L to be expressed in type II collagen-producing chondrocyte-like fibroblasts in human OPLL specimens, as well as in differentiating ATDC5 chondrocytes. Cdc5l siRNA transfection in murine chondrocytes decreased the expression of the early chondrogenic genes Sox9 and Col2a1, diminished the cartilage matrix production, and enhanced the expression of parathyroid-hormone-related protein (a resting chondrocyte marker). We also showed that Cdc5l shRNA suppressed the growth of cultured murine embryonal metatarsal cartilage rudiments and that Cdc5l knockdown suppressed the growth of ATDC5 cells. Fluorescence-activated cell sorting analysis revealed that the G2/M cell cycle transition was blocked; our data showed that Cdc5l siRNA transfection enhanced expression of Wee1, an inhibitor of the G2/M transition. Cdc5l siRNA also decreased the pre-mRNA splicing efficiency of Sox9 and Col2a1 genes in both ATDC5 cells and primary chondrocytes; conversely, loss of Cdc5l resulted in enhanced splicing of Wee1 pre-mRNA. Finally, an RNA-binding protein immunoprecipitation assay revealed that Cdc5l bound directly to these target gene transcripts. Overall, we conclude that Cdc5l promotes both early chondrogenesis and cartilage growth and may play a role in the etiology of OPLL, at least in part by fine-tuning the pre-mRNA splicing of chondrogenic genes and Wee1, thus initiating the endochondral ossification process.

Scoliosis represents the most common musculoskeletal disorder in children and affects approximately 3% of the world population. Scoliosis is separated into two major phenotypic classifications: congenital and idiopathic. Idiopathic scoliosis is defined as a curvature of the spine of 10° or greater visualized on plane radiograph and does not have associated vertebral malformations (VM). “Congenital” scoliosis (CS) due to malformations in vertebrae is frequently associated with other birth defects. Recently, significant advances have been made in understanding the genetic basis of both conditions. There is evidence that both conditions are etiologically related. A 2-day conference entitled “Genomic Approaches to Understanding and Treating Scoliosis” was held at Scottish Rite Hospital for Children in Dallas, Texas, to synergize research in this field. This first combined, multidisciplinary conference featured international scoliosis researchers in basic and clinical sciences. A major outcome of the conference advancing scoliosis research was the proposal and subsequent vote in favor of merging the International Consortium for Vertebral Anomalies and Scoliosis (ICVAS) and International Consortium for Scoliosis Genetics (ICSG) into a single entity called International Consortium for Spinal Genetics, Development, and Disease (ICSGDD). The ICSGDD is proposed to meet annually as a forum to synergize multidisciplinary spine deformity research.

Adolescent idiopathic scoliosis (AIS) is a common spinal deformity affecting millions of children. Since treatment and prognosis of AIS depend on curve progression, identifying factors related to AIS curve progression is important in its management. Although several genetic loci for AIS occurrence are reported, no locus for curve progression has been identified. To identify genes associated with AIS progression, we conducted a genome-wide association study followed by a replication study using a total of 2,543 AIS subjects who were evaluated for the curve progression. We identified a significantly associated locus on chromosome 11q14.1 (P = 1.98 × 10-9, odds ratio = 1.56). In silico and in vitro analyses identified a functional variant, rs35333564 in MIR4300HG, the host gene of a microRNA, MIR4300. The genomic region containing rs35333564 had enhancer activity, which was decreased in its risk allele. Our data suggest that decrease of MIR4300 is related to AIS progression.

Background Lumbar disc degeneration (LDD) is known to be the main cause of low back pain, although it is unclear how or when LDD progresses. The purpose of this study was to clarify the process and features of LDD progression in young women. Methods We enrolled 480 nursing students and carried out a prospective cohort study in 84 nursing students. MRI of the lumbar spine was taken once during their time as a student and again at 9.8 years (7–14) after the first MRI when they were working as nurses. The grade of LDD was determined according to Schneiderman's grade. The progression of disc degeneration was evaluated by the change in the degenerative disc disease (DDD) score (the summation of DDD score at each disc level). The subjects were divided into two groups based on MRI findings: Group A, those without disc degeneration at the first MRI (n = 58) and Group B, those showing disc degeneration on the first MRI (n = 26). We evaluated the change in DDD score and assessed the particular disc levels that showed Schneiderman's grade worsening in each group. Results At the time of the 2nd MRI, the DDD score increased from 5.6 (5–9) to 6.3 (5–11). The L5/S disc was the most frequent level showing the progression of LDD. In Group A, 18 subjects (31.0%) and in Group B, 15 subjects (57.7%) had degeneration progression (p = 0.02). Conclusions This study revealed that 31% of the young adult subjects already had disc degeneration in 20's (time of first MRI) and the disc degeneration rapidly progressed in these subjects.