RAB1A haploinsufficiency phenocopies the 2p14-p15 microdeletion and is associated with impaired neuronal differentiation.

Hereditary spastic parapareses (HSPs) are clinically heterogeneous motor neuron diseases with variable age of onset and severity. Although variants in dozens of genes are implicated in HSPs, much of the genetic basis for pediatric-onset HSP remains unexplained. Here, we re-analyzed clinical exome-sequencing data from siblings with HSP of unknown genetic etiology and identified an inherited nonsense mutation (c.523C>T [p.Arg175Ter]) in the highly conserved RAB1A. The mutation is predicted to produce a truncated protein with an intact RAB GTPase domain but without two C-terminal cysteine residues required for proper subcellular protein localization. Additional RAB1A mutations, including two frameshift mutations and a mosaic missense mutation (c.83T>C [p.Leu28Pro]), were identified in three individuals with similar neurodevelopmental presentations. In rescue experiments, production of the full-length, but not the truncated, RAB1a rescued Golgi structure and cell proliferation in Rab1-depleted cells. In contrast, the missense-variant RAB1a disrupted Golgi structure despite intact Rab1 expression, suggesting a dominant-negative function of the mosaic missense mutation. Knock-down of RAB1A in cultured human embryonic stem cell-derived neurons resulted in impaired neuronal arborization. Finally, RAB1A is located within the 2p14-p15 microdeletion syndrome locus. The similar clinical presentations of individuals with RAB1A loss-of-function mutations and the 2p14-p15 microdeletion syndrome implicate loss of RAB1A in the pathogenesis of neurodevelopmental manifestations of this microdeletion syndrome. Our study identifies a RAB1A-related neurocognitive disorder with speech and motor delay, demonstrates an essential role for RAB1a in neuronal differentiation, and implicates RAB1A in the etiology of the neurodevelopmental sequelae associated with the 2p14-p15 microdeletion syndrome.

Dominant missense variants in SREBF2 are associated with complex dermatological, neurological, and skeletal abnormalities.

PURPOSE: We identified 2 individuals with de novo variants in SREBF2 that disrupt a conserved site 1 protease (S1P) cleavage motif required for processing SREBP2 into its mature transcription factor. These individuals exhibit complex phenotypic manifestations that partially overlap with sterol regulatory element binding proteins (SREBP) pathway-related disease phenotypes, but SREBF2-related disease has not been previously reported. Thus, we set out to assess the effects of SREBF2 variants on SREBP pathway activation. METHODS: We undertook ultrastructure and gene expression analyses using fibroblasts from an affected individual and utilized a fly model of lipid droplet (LD) formation to investigate the consequences of SREBF2 variants on SREBP pathway function. RESULTS: We observed reduced LD formation, endoplasmic reticulum expansion, accumulation of aberrant lysosomes, and deficits in SREBP2 target gene expression in fibroblasts from an affected individual, indicating that the SREBF2 variant inhibits SREBP pathway activation. Using our fly model, we discovered that SREBF2 variants fail to induce LD production and act in a dominant-negative manner, which can be rescued by overexpression of S1P. CONCLUSION: Taken together, these data reveal a mechanism by which SREBF2 pathogenic variants that disrupt the S1P cleavage motif cause disease via dominant-negative antagonism of S1P, limiting the cleavage of S1P targets, including SREBP1 and SREBP2.

Mutations in genes encoding condensin complex proteins cause microcephaly through decatenation failure at mitosis.

Compaction of chromosomes is essential for accurate segregation of the genome during mitosis. In vertebrates, two condensin complexes ensure timely chromosome condensation, sister chromatid disentanglement, and maintenance of mitotic chromosome structure. Here, we report that biallelic mutations in NCAPD2, NCAPH, or NCAPD3, encoding subunits of these complexes, cause microcephaly. In addition, hypomorphic Ncaph2 mice have significantly reduced brain size, with frequent anaphase chromatin bridge formation observed in apical neural progenitors during neurogenesis. Such DNA bridges also arise in condensin-deficient patient cells, where they are the consequence of failed sister chromatid disentanglement during chromosome compaction. This results in chromosome segregation errors, leading to micronucleus formation and increased aneuploidy in daughter cells. These findings establish "condensinopathies" as microcephalic disorders, with decatenation failure as an additional disease mechanism for microcephaly, implicating mitotic chromosome condensation as a key process ensuring mammalian cerebral cortex size.

Genomic characterization of the adolescent idiopathic scoliosis-associated transcriptome and regulome.

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is the most common pediatric musculoskeletal disorder, affecting ~3% of the population worldwide. However, its genetic bases and tissues of origin remain largely unknown. Several genome-wide association studies (GWAS) have implicated nucleotide variants in non-coding sequences that control genes with important roles in cartilage, muscle, bone, connective tissue and intervertebral disks (IVDs) as drivers of AIS susceptibility. Here, we set out to define the expression of AIS-associated genes and active regulatory elements by performing RNA-seq and chromatin immunoprecipitation-sequencing against H3 lysine 27 acetylation in these tissues in mouse and human. Our study highlights genetic pathways involving AIS-associated loci that regulate chondrogenesis, IVD development and connective tissue maintenance and homeostasis. In addition, we identify thousands of putative AIS-associated regulatory elements which may orchestrate tissue-specific expression in musculoskeletal tissues of the spine. Quantification of enhancer activity of several candidate regulatory elements from our study identifies three functional enhancers carrying AIS-associated GWAS SNPs at the ADGRG6 and BNC2 loci. Our findings provide a novel genome-wide catalog of AIS-relevant genes and regulatory elements and aid in the identification of novel targets for AIS causality and treatment.

Do Children With Medicaid Insurance Have Increased Revision Rates 5 Years After Posterior Spinal Fusions?

BACKGROUND: Socioeconomic disparities in musculoskeletal care are increasingly recognized, however, no studies to date have investigated the role of the insurance carrier on outcomes after posterior spinal fusion (PSF) with segmental spinal instrumentation for adolescent idiopathic scoliosis (AIS). METHODS: A US insurance dataset was queried using the PearlDiver Mariner software for all patients aged 10 to 18 undergoing PSF for a primary diagnosis of AIS between 2010 and 2020. Age, sex, geographic region, number of levels fused, and baseline medical comorbidities were queried. Complications occurring within 90 days of the index surgery were queried using the International Classification of Diseases, Ninth Revision (ICD-9) and International Classification of Diseases, 10th Revision (ICD-10) codes. Revision surgery was also queried up to 5 years after the index PSF. Categorical variables were compared using the Fisher χ 2 tests and continuous variables were compared using independent t tests. All-cause revision within 5 years was compared using the Kaplan-Meier analysis and a log-rank test. Significance was set at P -value <0.05. RESULTS: A total of 10,794 patients were identified with 9006 (83.4%) patients with private insurance and 1788 (16.6%) patients insured by Medicaid. The mean follow-up in the database was 5.36±3 years for patients with private insurance and 4.78±2.9 years for patients with Medicaid insurance ( P <0.001). Children with AIS and Medicaid insurance had a significantly higher prevalence of asthma, hypertension, and obesity. A larger percentage of children with Medicaid insurance (41.3%) underwent a ≥13-level PSF compared with privately insured children (34.5%) ( P <0.001). Medicaid patients did not experience higher odds of postoperative complications; in addition, revision surgeries occurred in 1.1% and 1.8% of patients with private insurance and Medicaid insurance, respectively at 5 years postoperatively ( P =0.223). CONCLUSION: Despite worse baseline comorbidities and longer fusion constructs, AIS patients insured with Medicaid did not have higher rates of complications or revisions at 5-year follow-up versus privately insured patients. LEVEL OF EVIDENCE: Level III-retrospective cohort study.

TBX6 missense variants expand the mutational spectrum in a non-Mendelian inheritance disease.

Congenital scoliosis (CS) is a birth defect with variable clinical and anatomical manifestations due to spinal malformation. The genetic etiology underlying about 10% of CS cases in the Chinese population is compound inheritance by which the gene dosage is reduced below that of haploinsufficiency. In this genetic model, the trait manifests as a result of the combined effect of a rare variant and common pathogenic variant allele at a locus. From exome sequencing (ES) data of 523 patients in Asia and two patients in Texas, we identified six TBX6 gene-disruptive variants from 11 unrelated CS patients via ES and in vitro functional testing. The in trans mild hypomorphic allele was identified in 10 of the 11 subjects; as anticipated these 10 shared a similar spinal deformity of hemivertebrae. The remaining case has a homozygous variant in TBX6 (c.418C>T) and presents a more severe spinal deformity phenotype. We found decreased transcriptional activity and abnormal cellular localization as the molecular mechanisms for TBX6 missense loss-of-function alleles. Expanding the mutational spectrum of TBX6 pathogenic alleles enabled an increased molecular diagnostic detection rate, provided further evidence for the gene dosage-dependent genetic model underlying CS, and refined clinical classification.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Distal chromosome 16p11.2 duplications containing SH2B1 in patients with scoliosis.

INTRODUCTION: Adolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. CNVs play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored. METHODS: Exome sequence data from 1197 AIS cases and 1664 in-house controls was analysed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with >10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics. RESULTS: Distal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harbouring a chromosome 16p11.2 microduplication, compared with 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584) (p=2.28×10-11, OR=16.15). Furthermore, examination of electronic health records of 92 455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared with 7.6% (10/132) of controls (p=5.6×10-4, OR=3.9). CONCLUSIONS: Recurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

Mutations Preventing Regulated Exon Skipping in MET Cause Osteofibrous Dysplasia.

The periosteum contributes to bone repair and maintenance of cortical bone mass. In contrast to the understanding of bone development within the epiphyseal growth plate, factors that regulate periosteal osteogenesis have not been studied as intensively. Osteofibrous dysplasia (OFD) is a congenital disorder of osteogenesis and is typically sporadic and characterized by radiolucent lesions affecting the cortical bone immediately under the periosteum of the tibia and fibula. We identified germline mutations in MET, encoding a receptor tyrosine kinase, that segregate with an autosomal-dominant form of OFD in three families and a mutation in a fourth affected subject from a simplex family and with bilateral disease. Mutations identified in all families with dominant inheritance and in the one simplex subject with bilateral disease abolished the splice inclusion of exon 14 in MET transcripts, which resulted in a MET receptor (MET(Δ14)) lacking a cytoplasmic juxtamembrane domain. Splice exclusion of this domain occurs during normal embryonic development, and forced induction of this exon-exclusion event retarded osteoblastic differentiation in vitro and inhibited bone-matrix mineralization. In an additional subject with unilateral OFD, we identified a somatic MET mutation, also affecting exon 14, that substituted a tyrosine residue critical for MET receptor turnover and, as in the case of the MET(Δ14) mutations, had a stabilizing effect on the mature protein. Taken together, these data show that aberrant MET regulation via the juxtamembrane domain subverts core MET receptor functions that regulate osteogenesis within cortical diaphyseal bone.

Summary of the first inaugural joint meeting of the International Consortium for scoliosis genetics and the International Consortium for vertebral anomalies and scoliosis, March 16-18, 2017, Dallas, Texas.

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.

Somatic gain-of-function mutations in PIK3CA in patients with macrodactyly.

Macrodactyly is a discrete congenital anomaly consisting of enlargement of all tissues localized to the terminal portions of a limb, typically within a 'nerve territory'. The classic terminology for this condition is 'lipofibromatous hamartoma of nerve' or Type I macrodactyly. The peripheral nerve, itself, is enlarged both in circumference and in length. It is not related to neurofibromatosis (NF1), nor is it associated with vascular malformations, such as in the recently reported CLOVES syndrome. The specific nerve pathophysiology in this form of macrodactyly has not been well described and a genetic etiology for this specific form of enlargement is unknown. To identify the genetic cause of macrodactyly, we used whole-exome sequencing to identify somatic mutations present in the affected nerve of a single patient. We confirmed a novel mutation in PIK3CA (R115P) present in the patient's affected nerve tissue but not in blood DNA. Sequencing PIK3CA exons identified gain-of-function mutations (E542K, H1047L or H1047R) in the affected tissue of five additional unrelated patients; mutations were absent in blood DNA available from three patients. Immunocytochemistry confirmed AKT activation in cultured cells from the nerve of a macrodactyly patient. Additionally, we found that the most abnormal structure within the involved nerve in a macrodactylous digit is the perineurium, with additional secondary effects on the axon number and size. Thus, isolated congenital macrodactyly is caused by somatic activation of the PI3K/AKT cell-signaling pathway and is genetically and biochemically related to other overgrowth syndromes.

Rare and common variants in CARD14, encoding an epidermal regulator of NF-kappaB, in psoriasis.

Psoriasis is a common inflammatory disorder of the skin and other organs. We have determined that mutations in CARD14, encoding a nuclear factor of kappa light chain enhancer in B cells (NF-kB) activator within skin epidermis, account for PSORS2. Here, we describe fifteen additional rare missense variants in CARD14, their distribution in seven psoriasis cohorts (>6,000 cases and >4,000 controls), and their effects on NF-kB activation and the transcriptome of keratinocytes. There were more CARD14 rare variants in cases than in controls (burden test p value = 0.0015). Some variants were only seen in a single case, and these included putative pathogenic mutations (c.424G>A [p.Glu142Lys] and c.425A>G [p.Glu142Gly]) and the generalized-pustular-psoriasis mutation, c.413A>C (p.Glu138Ala); these three mutations lie within the coiled-coil domain of CARD14. The c.349G>A (p.Gly117Ser) familial-psoriasis mutation was present at a frequency of 0.0005 in cases of European ancestry. CARD14 variants led to a range of NF-kB activities; in particular, putative pathogenic variants led to levels >2.5× higher than did wild-type CARD14. Two variants (c.511C>A [p.His171Asn] and c.536G>A [p.Arg179His]) required stimulation with tumor necrosis factor alpha (TNF-α) to achieve significant increases in NF-kB levels. Transcriptome profiling of wild-type and variant CARD14 transfectants in keratinocytes differentiated probably pathogenic mutations from neutral variants such as polymorphisms. Over 20 CARD14 polymorphisms were also genotyped, and meta-analysis revealed an association between psoriasis and rs11652075 (c.2458C>T [p.Arg820Trp]; p value = 2.1 × 10(-6)). In the two largest psoriasis cohorts, evidence for association increased when rs11652075 was conditioned on HLA-Cw*0602 (PSORS1). These studies contribute to our understanding of the genetic basis of psoriasis and illustrate the challenges faced in identifying pathogenic variants in common disease.

A meta-analysis identifies adolescent idiopathic scoliosis association with LBX1 locus in multiple ethnic groups.

BACKGROUND: Adolescent idiopathic scoliosis (AIS) is a common rotational deformity of the spine that presents in children worldwide, yet its etiology is poorly understood. Recent genome-wide association studies (GWAS) have identified a few candidate risk loci. One locus near the chromosome 10q24.31 LBX1 gene (OMIM #604255) was originally identified by a GWAS of Japanese subjects and replicated in additional Asian populations. To extend this result, and to create larger AIS cohorts for the purpose of large-scale meta-analyses in multiple ethnicities, we formed a collaborative group called the International Consortium for Scoliosis Genetics (ICSG). METHODS: Here, we report the first ICSG study, a meta-analysis of the LBX1 locus in six Asian and three non-Asian cohorts. RESULTS: We find significant evidence for association of this locus with AIS susceptibility in all nine cohorts. Results for seven cohorts containing both genders yielded P=1.22×10-43 for rs11190870, and P=2.94×10-48 for females in all nine cohorts. Comparing the regional haplotype structures for three populations, we refined the boundaries of association to a ∼25 kb block encompassing the LBX1 gene. The LBX1 protein, a homeobox transcription factor that is orthologous to the Drosophila ladybird late gene, is involved in proper migration of muscle precursor cells, specification of cardiac neural crest cells, and neuronal determination in developing neural tubes. CONCLUSIONS: Our results firmly establish the LBX1 region as the first major susceptibility locus for AIS in Asian and non-Hispanic white groups, and provide a platform for larger studies in additional ancestral groups.

Extreme growth failure is a common presentation of ligase IV deficiency.

Ligase IV syndrome is a rare differential diagnosis for Nijmegen breakage syndrome owing to a shared predisposition to lympho-reticular malignancies, significant microcephaly, and radiation hypersensitivity. Only 16 cases with mutations in LIG4 have been described to date with phenotypes varying from malignancy in developmentally normal individuals, to severe combined immunodeficiency and early mortality. Here, we report the identification of biallelic truncating LIG4 mutations in 11 patients with microcephalic primordial dwarfism presenting with restricted prenatal growth and extreme postnatal global growth failure (average OFC -10.1 s.d., height -5.1 s.d.). Subsequently, most patients developed thrombocytopenia and leucopenia later in childhood and many were found to have previously unrecognized immunodeficiency following molecular diagnosis. None have yet developed malignancy, though all patients tested had cellular radiosensitivity. A genotype-phenotype correlation was also noted with position of truncating mutations corresponding to disease severity. This work extends the phenotypic spectrum associated with LIG4 mutations, establishing that extreme growth retardation with microcephaly is a common presentation of bilallelic truncating mutations. Such growth failure is therefore sufficient to consider a diagnosis of LIG4 deficiency and early recognition of such cases is important as bone marrow failure, immunodeficiency, and sometimes malignancy are long term sequelae of this disorder.

A novel method for analyzing genetic association with longitudinal phenotypes.

Knowledge of genes influencing longitudinal patterns may offer information about predicting disease progression. We developed a systematic procedure for testing association between SNP genotypes and longitudinal phenotypes. We evaluated false positive rates and statistical power to localize genes for disease progression. We used genome-wide SNP data from the Framingham Heart Study. With longitudinal data from two real studies unrelated to Framingham, we estimated three trajectory curves from each study. We performed simulations by randomly selecting 500 individuals. In each simulation replicate, we assigned each individual to one of the three trajectory groups based on the underlying hypothesis (null or alternative), and generated corresponding longitudinal data. Individual Bayesian posterior probabilities (BPPs) for belonging to a specific trajectory curve were estimated. These BPPs were treated as a quantitative trait and tested (using the Wald test) for genome-wide association. Empirical false positive rates and power were calculated. Our method maintained the expected false positive rate for all simulation models. Also, our method achieved high empirical power for most simulations. Our work presents a method for disease progression gene mapping. This method is potentially clinically significant as it may allow doctors to predict disease progression based on genotype and determine treatment accordingly.

Susceptibility to childhood-onset rheumatoid arthritis: investigation of a weighted genetic risk score that integrates cumulative effects of variants at five genetic loci.

OBJECTIVE: Children with childhood-onset rheumatoid arthritis (RA) include those with rheumatoid factor or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. To test the hypothesis that adult-onset RA-associated variants are also associated with childhood-onset RA, we investigated RA-associated variants at 5 loci in a cohort of patients with childhood-onset RA. We also assessed the cumulative association of these variants in susceptibility to childhood-onset RA using a weighted genetic risk score (wGRS). METHODS: A total of 155 children with childhood-onset RA and 684 healthy controls were genotyped for 5 variants in the PTPN22, TRAF1/C5, STAT4, and TNFAIP3 loci. High-resolution HLA-DRB1 genotypes were available for 149 cases and 373 controls. We tested each locus for association with childhood-onset RA via logistic regression. We also computed a wGRS for each subject, with weights based on the natural log of the published odds ratios (ORs) for the alleles investigated, and used logistic regression to test the wGRS for association with childhood-onset RA. RESULTS: Childhood-onset RA was associated with TNFAIP3 rs10499194 (OR 0.60 [95% confidence interval 0.44-0.83]), PTPN22 rs2476601 (OR 1.61 [95% confidence interval 1.11-2.31]), and STAT4 rs7574865 (OR 1.41 [95% confidence interval 1.06-1.87]) variants. The wGRS was significantly different between cases and controls (P < 2 × 10(-16) ). Individuals in the third to fifth quintiles of wGRS had a significantly increased disease risk compared to baseline (individuals in the first quintile). Higher wGRS was associated with increased risk of childhood-onset RA, especially among males. CONCLUSION: The magnitude and direction of the association between TNFAIP3, STAT4, and PTPN22 variants and childhood-onset RA are similar to those observed in RA, suggesting that adult-onset RA and childhood-onset RA share common genetic risk factors. Using a wGRS, we have demonstrated the cumulative association of RA-associated variants with susceptibility to childhood-onset RA.

Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality.

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A genome-wide association study identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we use mouse enhancer assays, three mouse enhancer knockouts, and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterize a sequence, PEC7, which overlaps the AIS-associated variant, and find it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, or deletion of both sequences lead to a kinky tail phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicates several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Spatial transcriptomics implicates impaired BMP signaling in NF1 fracture pseudarthrosis in murine and patient tissues.

The neurofibromatosis type 1 (NF1) RASopathy is associated with persistent fibrotic nonunions (pseudarthrosis) in human and mouse skeletal tissue. Here, we performed spatial transcriptomics to define the molecular signatures occurring during normal endochondral healing following fracture in mice. Within the control fracture callus, we observed spatially restricted activation of morphogenetic pathways, such as TGF-ß, WNT, and BMP. To investigate the molecular mechanisms contributing to Nf1-deficient delayed fracture healing, we performed spatial transcriptomic analysis on a Postn-cre;Nf1fl/- (Nf1Postn) fracture callus. Transcriptional analyses, subsequently confirmed through phospho-SMAD1/5/8 immunohistochemistry, demonstrated a lack of BMP pathway induction in Nf1Postn mice. To gain further insight into the human condition, we performed spatial transcriptomic analysis of fracture pseudarthrosis tissue from a patient with NF1. Analyses detected increased MAPK signaling at the fibrocartilaginous-osseus junction. Similar to that in the Nf1Postn fracture, BMP pathway activation was absent within the pseudarthrosis tissue. Our results demonstrate the feasibility of delineating the molecular and tissue-specific heterogeneity inherent in complex regenerative processes, such as fracture healing, and reconstructing phase transitions representing endochondral bone formation in vivo. Furthermore, our results provide in situ molecular evidence of impaired BMP signaling underlying NF1 pseudarthrosis, potentially informing the clinical relevance of off-label BMP2 as a therapeutic intervention.

Deep Learning-Based Automated Measurement of Murine Bone Length in Radiographs.

Genetic mouse models of skeletal abnormalities have demonstrated promise in the identification of phenotypes relevant to human skeletal diseases. Traditionally, phenotypes are assessed by manually examining radiographs, a tedious and potentially error-prone process. In response, this study developed a deep learning-based model that streamlines the measurement of murine bone lengths from radiographs in an accurate and reproducible manner. A bone detection and measurement pipeline utilizing the Keypoint R-CNN algorithm with an EfficientNet-B3 feature extraction backbone was developed to detect murine bone positions and measure their lengths. The pipeline was developed utilizing 94 X-ray images with expert annotations on the start and end position of each murine bone. The accuracy of our pipeline was evaluated on an independent dataset test with 592 images, and further validated on a previously published dataset of 21,300 mouse radiographs. The results showed that our model performed comparably to humans in measuring tibia and femur lengths (R2 > 0.92, p-value = 0) and significantly outperformed humans in measuring pelvic lengths in terms of precision and consistency. Furthermore, the model improved the precision and consistency of genetic association mapping results, identifying significant associations between genetic mutations and skeletal phenotypes with reduced variability. This study demonstrates the feasibility and efficiency of automated murine bone length measurement in the identification of mouse models of abnormal skeletal phenotypes.