Adolescent idiopathic scoliosis (AIS) is the most common spinal disorder in children and is best understood as a complex polygenic condition. Genome-wide association studies (GWAS) have identified several AIS risk loci, including regions near ADGRG6 and SOX9, yet the functional mechanisms underlying AIS heritability remain poorly defined. Here, we use spatial transcriptomics to characterize altered gene expression in the spine of a conditional Adgrg6 mutant mouse model of AIS (Adgrg6-cKO), revealing reduced expression of Sox9 and several components of extracellular matrix organization in the intervertebral disc. We further show that SOX9 occupies regions of open chromatin within the Adgrg6 locus in cells isolated from the mouse intervertebral disc. Finally, we demonstrate a strong genetic interaction between Adgrg6-cKO and a hypomorphic Sox9 allele that increases both the penetrance and severity of AIS-like pathology in mice. Collectively, these findings support a self-reinforcing feedforward regulatory circuit, where Adgrg6 and Sox9 are co-regulated to maintain extracellular matrix gene expression in the annulus fibrosus and paraspinal tissues. These findings provide mechanistic insight into the functional significance of AIS-associated GWAS loci near ADGRG6 and SOX9 and establish combined Adgrg6-Sox9 insufficiency as a tractable model of polygenic scoliosis susceptibility.

Protein-protein interactions underlie core brain functions, including neurotransmitter release, receptor activation, and intracellular signaling. Here, we present a compendium of protein-protein interactions conserved across vertebrate brains, spanning over 300 million years of evolution. From 2,197 biochemical fractions across five vertebrate species, we identify over 81,000 high-confidence interactions among 6,108 conserved proteins. This interaction map (VerteBrain) reveals both regulatory and structural complexes, including extensive synaptonemal protein associations likely involved in inter-neuronal coordination. Conservation across species underscores essential roles in neuronal and glial function, as well as in additional tissues for more widely expressed complexes. The VerteBrain dataset uncovers candidate disease mechanisms, including roles for ARHGEF1 in short stature syndromes, synaptic vesicle trafficking complexes in epilepsy, and RELCH (Rab11-binding and LisH domain, coiled-coil and HEAT repeat-containing) in congenital deafness. VerteBrain provides a resource for investigating brain protein interactions and their relevance to human neurological disorders.

Here, we present a protocol for quantifying zebrafish vertebral column morphology using micro-computed tomography (microCT) datasets for scoliosis researchers. We describe steps for orientation of skeletal and cranial microCT datasets in DataViewer, vertebral position measurements in ImageJ/Fiji, and analysis of vertebral position in R. The approach described here graphically depicts vertebral position while providing a means to quantify scoliosis severity and perform statistical analysis of differences in vertebral column morphology in zebrafish.For complete details on the use and execution of this protocol, please refer to Voigt et al.1

The TRPV4 skeletal dysplasias are characterized by short stature, short limbs with prominent large joints, and progressive scoliosis. They result from dominant missense mutations that activate the TRPV4 calcium permeable ion channel. As a platform to understand the mechanism of disease and to test the hypothesis that channel inhibition could treat these disorders, we developed a knock-in mouse that conditionally expresses the p.R594H Trpv4 mutation. Embryonic, chondrocyte-specific induction of the mutation using Col2a1-Cre resulted in a skeletal dysplasia affecting the long bones, spine, and craniofacial skeletal elements, consistent with the human skeletal dysplasia phenotypes produced by TRPV4 mutations. Cartilage growth plate histological abnormalities included disorganized proliferating chondrocyte columns and reduced hypertrophic chondrocyte development, reflecting abnormal endochondral ossification. In vivo treatment with the TRPV4-specific inhibitor GSK2798745 markedly improved the radiographic skeletal phenotype and rescued the growth plate histological abnormalities. ScRNA-Seq of chondrocyte transcripts from affected mice identified calcium-mediated effects on multiple signaling pathways as potential mechanisms underlying the defects in linear and cartilage appositional growth observed in both mutant mice and patients. These results provide preclinical evidence demonstrating TRPV4 inhibition as a rational, mechanism-based therapeutic strategy to ameliorate disease progression and severity in the TRPV4 skeletal dysplasias.

Scoliosis is the most common developmental spinal deformity, but its genetic underpinnings remain only partially understood. To identify scoliosis-related loci, we utilized dual energy X-ray absorptiometry (DXA) scans from 57,588 individuals in the UK Biobank (UKB), and quantified spinal curvature using deep learning-based vertebral segmentation and landmarking to measure cumulative horizontal displacement. On a subset of 150 individuals, our automated image-derived curvature measurements showed a correlation of 0.83 with clinical Cobb angle assessments, supporting their validity as a proxy for scoliosis severity. To connect spinal curvature to genetics, we conducted a genome-wide association study (GWAS). Our quantitative imaging phenotype identified 2 novel loci associated with scoliosis in a European population. These loci are in SEM1/SHFM1 and on an lncRNA on chr 3 located between EDEM1 and GRM7. Genetic correlation analysis revealed significant overlap between our image-based GWAS and ICD-10-based GWAS in the UKB and the Biobank of Japan. We show that our quantitative GWAS identifies more genome-wide significant loci than a case-control scoliosis dataset with ten times the sample size. Our results illustrate the potential of quantitative imaging phenotypes to uncover genetic associations that are challenging to capture with medical records alone and identify new loci for functional follow-up.

Motile cilia are evolutionarily conserved organelles performing essential roles in development and tissue homeostasis. Unlike the core scaffold, the distal regions remain relatively less explored and display great diversity across species. Here, we describe a previously uncharacterized ciliary protein Jhc1 (just the head of cilia 1, Loc108698169), localizing at the extreme distal tip of multiciliated cell (MCC) cilia and encoded only in the genomes of non-mammalian vertebrates. Jhc1 is essential for normal cilia structure and function in Xenopus, and this activity is conserved in Jhc1 from reptiles and fish. Phylogenetic analysis and structure modeling suggest that Jhc1 arose by duplication and neofunctionalization of thiamine triphosphatase, and the residues crucial for that enzyme's function have been lost and replaced by residues essential for ciliary localization. These data provide insights into the molecular mechanisms underlying the broad diversification of the structures at the tip of motile cilia during vertebrate evolution.

Cell size is a key contributor to tissue morphogenesis. As a notable example, growth plate hypertrophic chondrocytes use cellular biogenesis and disproportionate fluid uptake to expand 10 to 20 times in size to drive lengthening of endochondral bone. Similarly, notochord vacuolated cells expand to one of the largest cell types in the developing embryo to drive axial extension. In zebrafish, the notochord vacuolated cells undergo vacuole fusion to form a single large, fluid-filled vacuole that fills the cytoplasmic space and contributes to vacuolated cell expansion. When this process goes awry, the notochord lacks sufficient hydrostatic pressure to support vertebral bone deposition, resulting in adult spines with misshapen vertebral bones and scoliosis. However, it remains unclear whether endochondral bone and the notochord share common genetic and cellular mechanisms for regulating cell and tissue expansion. Here, we demonstrate that the 5'-inositol phosphatase gene, inppl1a, regulates notochord expansion independent of vacuole fusion, thereby genetically decoupling these processes. We demonstrate that inppl1a-dependent vacuolated cell expansion is essential to establish normal mechanical properties of the notochord and to facilitate the development of a straight spine. Finally, we find that inppl1a is also important for hypertrophic chondrocyte differentiation and endochondral bone lengthening in fish, as has been shown in the human INPPL1-related endochondral bone disorder, opsismodysplasia. Overall, this work reveals a shared mechanism of cell size regulation that influences disparate tissues critical for skeletal development and short-stature disorders.

Protein-protein interactions underlie core brain functions, including neurotransmitter release, receptor activation, and intracellular signaling essential for learning, memory, and cognition. Here, we systematically map conserved brain protein interactions across five vertebrate species-rabbit, chicken, dolphin, pig, and mouse-using co-fractionation and immunoprecipitation mass spectrometry. From 2,197 biochemical fractions, we identify over 81,000 high-confidence interactions among 6,108 conserved proteins. This interaction map (VerteBrain) reveals both regulatory and structural complexes, including extensive synaptonemal protein associations likely involved in inter-neuronal coordination. Conservation across species underscores essential roles in neuronal and glial function, as well as in additional tissues for more widely expressed complexes. The VerteBrain dataset uncovers candidate disease mechanisms, including roles for ARHGEF1 in short stature syndromes, synaptic vesicle trafficking complexes in epilepsy, and RELCH in congenital deafness. VerteBrain provides a publicly accessible framework for investigating brain protein interactions and their relevance to human neurological disorders.

Motile cilia coordinate fluid flows that are essential for normal tissue physiology and function. Cilia display diverse beating waveforms, and while pronounced defects are strongly associated with motile ciliopathies, more subtle alterations may also influence disease manifestations1-3. Finer quantification of ciliary dynamics is critical for a full understanding of cilia-associated disorders, but the heterogeneity of cilia beating dynamics makes accurate and robust characterization challenging. Because existing tools have proven to be limiting in noisy in vivo environments4-6, we developed Cilia.io, a machine learning (ML)-based quantification tool that uses state-of-the-art vision transformers to segment cilia out of the background based on their biological features. Cilia.io enables fast, accurate, and reproducible quantification of motile cilia morphodynamics and outperforms existing tools. Indeed, using Cilia.io, we discovered distinct regional differences in ciliary waveforms in the zebrafish spinal cord. Cilia.io's ability to capture subtle ciliary defects was further demonstrated by analyzing a novel allele in the ciliopathy gene bbs2 that causes a highly heterogeneous scoliosis phenotype. In these mutants, only dorsal cilia displayed altered beating dynamics, while ventral cilia remained largely unaffected. Our new tool therefore represents a substantial advance on existing methods and suggests that additional fine-scale analyses of ciliary beating will be important for understanding organismal phenotypes and cilia-driven disease.

Retrons can produce multicopy single-stranded DNA in cells through self-primed reverse transcription. However, their potential for inserting genetic cargos in eukaryotes remains largely unexplored. Here we report the discovery and engineering of highly efficient retron-based gene editors for mammalian cells and vertebrates. Through bioinformatic analysis of metagenomic data and functional screening, we identify retron reverse transcriptases that are highly active in mammalian cells. Rational design further improves the editing efficiency to levels comparable with conventional single-stranded oligodeoxynucleotide donors but from a genetically encoded cassette. Retron editors exhibit robust activity with Cas12a nuclease and Cas9 nickase, expanding the genomic target scope and bypassing the need for a DNA double-stranded break. Using a rationally engineered retron editor, we incorporate a split GFP epitope tag for live-cell imaging. Lastly, we develop an all-RNA delivery strategy to enable DNA-free gene editing in cells and vertebrate embryos. This work establishes retron editors as a versatile and efficient tool for precise gene editing. A metagenomic screen identifies retron reverse transcriptases for precise genome-editing applications.

Oxygen deprivation within large or poorly vascularized bone defects remains a key barrier to successful regeneration, especially during the early postimplantation period before vascular ingrowth. Here, the development of COSnPPOD (CaO2-Silica NP Platform for Osteogenic Development) is reported, a visible light digital light processing-printed hydrogel scaffold that integrates oxygen-releasing nanoparticles (NPs) within a Primitive-type triply periodic minimal surface architecture. The scaffold combines a gelatin methacrylate-poly(ethylene glycol) diacrylate matrix with calcium peroxide (CaO2)-loaded hollow silica NPs, enabling localized, short-term oxygen release while preserving structural fidelity. COSnPPOD scaffolds demonstrate favorable degradation kinetics, tunable stiffness, and increased protein adsorption in vitro. In a preosteoblast model, COSnPPOD maintains cell viability and supports osteogenic gene expression without cytotoxic effects. While overall gene expression is comparable to controls, a 16-fold increased expression of phosphoprotein 1 (Spp1) suggests scaffold-driven activation of matrix remodeling pathways. In vivo, COSnPPOD scaffolds enhance bone regeneration in a murine calvarial defect model, with significantly greater bone formation and collagen deposition than untreated defects and hydrogel controls. Additionally, vascular endothelial growth factor immunostaining is increased within the defect, consistent with a proangiogenic response, and no systemic toxicity is observed. These findings establish COSnPPOD as a promising scaffold system that combines sustained oxygenation with biomimetic geometry to support localized bone regeneration.

SOX9 is a crucial transcriptional regulator of cartilage development and homeostasis. Dysregulation of SOX9 is associated with a wide spectrum of skeletal disorders, including campomelic dysplasia, acampomelic campomelic dysplasia, and scoliosis. Yet how SOX9 variants contribute to the spectrum of axial skeletal disorders is not well understood. Here, we report four pathogenic variants of SOX9 identified in a cohort of patients with congenital vertebral malformations. We report a pathogenic missense variant in the transactivation middle (TAM) domain of SOX9 associated with mild skeletal dysplasia and scoliosis. We isolated a Sox9 mutant mouse with an in-frame microdeletion in the TAM domain ( Sox9 Asp272del ), which exhibits skeletal dysplasia including kinked tails, rib cage anomalies, and scoliosis in homozygous mutants. We find that both the human missense and the mouse microdeletion mutations resulted in reduced SOX9 protein stability in cell culture, while Sox9 Asp272del mutant mice show decreased SOX9 expression in the growth plate and annulus fibrosus tissues of the spine. This reduction in SOX9 expression was correlated with the reduction of extracellular matrix components, such as tenascin-X and the Adhesion G-protein coupled receptor ADGRG6. In summary, our work identified and modeled a pathologic variant of SOX9 within the TAM domain and demonstrated its importance for SOX9 protein stability. Our work demonstrates that SOX9 stability is important for the regulation of ADGRG6 expression, which is a known regulator of postnatal spine homeostasis, underscoring the essential role of SOX9 dosage in a spectrum of axial skeleton dysplasia in humans.

Adolescent idiopathic scoliosis (AIS) is a twisting spinal deformity that occurs in otherwise healthy children at the time of rapid pre-pubescent growth. AIS affects ∼3% of children worldwide and is the most common musculoskeletal diagnosis in pediatric populations, posing a significant physiological, psychosocial and financial burden to patients. Genetic predisposition is a clear and major contributor to AIS, and insights from genomic discoveries are inspiring translational studies ultimately aimed at developing novel diagnostics and therapies. Pre-clinical animal models of AIS are now essential to validate human genetic findings, understand gene-by-environment interactions, and speed etiologic and therapeutic discovery. In this Perspective, we highlight the current status of pre-clinical models of AIS and discuss the challenges posed by the nature of the disorder combined with the limitations of standard approaches. Current research suggests that straightforward genetic targeting of orthologous AIS disease genes in vertebrates may not necessarily yield equivalent physiological phenotypes but nevertheless can be utilized to understand disease mechanisms. Longer-term, appropriately complex models are needed to fully recapitulate the human AIS phenotype arising from genetic, physiological and mechanical interactions.

Cell size is a key contributor to tissue morphogenesis 1 . As a notable example, growth plate hypertrophic chondrocytes use cellular biogenesis and disproportionate fluid uptake to expand 10-20 times in size to drive lengthening of endochondral bone 2,3 . Similarly, notochordal cells expand to one of the largest cell types in the developing embryo to drive axial extension 4-6 . In zebrafish, the notochord vacuolated cells undergo vacuole fusion to form a single large, fluid-filled vacuole that fills the cytoplasmic space and contributes to vacuolated cell expansion 7 . When this process goes awry, the notochord lacks sufficient hydrostatic pressure to support vertebral bone deposition resulting in adult spines with misshapen vertebral bones and scoliosis 8 . However, it remains unclear whether endochondral bone and the notochord share common genetic and cellular mechanisms for regulating cell and tissue expansion. Here, we demonstrate that the 5'-inositol phosphatase gene, inppl1a , regulates notochord expansion, spine morphogenesis, and endochondral bone lengthening in zebrafish. Furthermore, we show that inppl1a regulates notochord expansion independent of vacuole fusion, thereby genetically decoupling these processes. We demonstrate that inppl1a -dependent notochord expansion is essential to establish normal mechanical properties of the notochord to facilitate the development of a straight spine. Finally, we find that inppl1a is also important for endochondral bone lengthening in fish, as has been shown in the human INPPL1 -related endochondral bone disorder, Opsismodysplasia 9 . Overall, this work reveals a conserved mechanism of cell size regulation that influences disparate tissues critical for skeletal development and short-stature disorders.

Ependymal cells, lining brain ventricular walls, display tufts of cilia that beat in concert promoting laminar Cerebrospinal fluid (CSF) flow within brain ventricles. The ciliary axonemes of multiciliated ependymal cells display a 9+2 microtubule array common to motile cilia. Dyneins and kinesins are ATPase microtubule motor proteins that promote the rhythmic beating of cilia axonemes. Despite common consensus about the importance of axonemal dynein motor proteins, little is known about how Kinesin motors contribute to cilia motility. Here, we define the function of Kinesin family member 6 (Kif6) using a mutation that lacks a highly conserved C-terminal tail domain ( Kif6 p.G555fs ) and which displays progressive hydrocephalus in mice. An analogous mutation was isolated in a proband displaying macrocephaly, hypotonia, and seizures implicating an evolutionarily conserved function for Kif6 in neurodevelopment. We find that loss of Kif6 function caused decreased ependymal cilia motility and subsequently decreased fluid flow on the surface of brain ventricular walls. Kif6 protein was localized at ependymal cilia and displayed processive motor movement (676 nm/s) on microtubules in vitro . Loss of the Kif6 C-terminal tail domain did not affect the initial ciliogenesis in vivo , but did result in defects in cilia orientation, the formation of robust apical actin networks, and stabilization of basal bodies at the apical surface. This suggests a novel role for the Kif6 motor in maintenance of ciliary homeostasis of ependymal cells. Summary statement We found that Kif6 is localized to the axonemes of ependymal cells. In vitro analysis shows that Kif6 moves on microtubules and that its loss mice decrease cilia motility and cilia-driven flow, resulting in hydrocephalus.

The cartilage growth plate is essential for maintaining skeletal growth; however, the mechanisms governing postnatal growth plate homeostasis are still poorly understood. Using approaches of molecular mouse genetics and spatial transcriptomics applied to formalin-fixed, paraffin-embedded tissues, we show that ADGRG6/GPR126, a cartilage-enriched adhesion G protein-coupled receptor (GPCR), is essential for maintaining slow-cycling resting zone cells, appropriate chondrocyte proliferation and differentiation, and growth plate homeostasis in mice. Constitutive ablation of Adgrg6 in osteochondral progenitor cells with Col2a1Cre leads to a shortened resting zone, formation of cell clusters within the proliferative zone, and an elongated hypertrophic growth plate, marked by limited expression of parathyroid hormone-related protein (PTHrP) but increased Indian Hedgehog (IHH) signaling throughout the growth plate. Attenuation of smoothened-dependent hedgehog signaling restored the Adgrg6 deficiency-induced expansion of hypertrophic chondrocytes, confirming that IHH signaling can promote chondrocyte hypertrophy in a PTHrP-independent manner. In contrast, postnatal ablation of Adgrg6 in mature chondrocytes with AcanCreERT2, induced after the formation of the resting zone, does not affect PTHrP expression but causes an overall reduction of growth plate thickness marked by increased cell death specifically in the resting zone cells and a general reduction of chondrocyte proliferation and differentiation. Spatial transcriptomics reveals that ADGRG6 is essential for maintaining chondrocyte homeostasis by regulating osteogenic and catabolic genes in all the zones of the postnatal growth plates, potentially through positive regulation of SOX9 expression. Our findings elucidate the essential role of a cartilage-enriched adhesion GPCR in regulating cell proliferation and hypertrophic differentiation by regulation of PTHrP/IHH signaling, maintenance of slow-cycle resting zone chondrocytes, and safeguarding chondrocyte homeostasis in postnatal mouse growth plates.

The Reissner fiber (RF) is an acellular thread positioned in the midline of the central canal that aggregates thanks to the beating of numerous cilia from ependymal radial glial cells (ERGs) generating flow in the central canal of the spinal cord. RF together with cerebrospinal fluid (CSF)-contacting neurons (CSF-cNs) form an axial sensory system detecting curvature. How RF, CSF-cNs and the multitude of motile cilia from ERGs interact in vivo appears critical for maintenance of RF and sensory functions of CSF-cNs to keep a straight body axis, but is not well-understood. Using in vivo imaging in larval zebrafish, we show that RF is under tension and resonates dorsoventrally. Focal RF ablations trigger retraction and relaxation of the fiber’s cut ends, with larger retraction speeds for rostral ablations. We built a mechanical model that estimates RF stress diffusion coefficient D at 5 mm2/s and reveals that tension builds up rostrally along the fiber. After RF ablation, spontaneous CSF-cN activity decreased and ciliary motility changed, suggesting physical interactions between RF and cilia projecting into the central canal. We observed that motile cilia were caudally-tilted and frequently interacted with RF. We propose that the numerous ependymal motile monocilia contribute to RF’s heterogenous tension via weak interactions. Our work demonstrates that under tension, the Reissner fiber dynamically interacts with motile cilia generating CSF flow and spinal sensory neurons.

The advent of targeted CRISPR-Cas nuclease technologies has revolutionized the ability to perform precise genome editing in both established and emerging model systems. CRISPR-Cas genome editing systems use a synthetic guide RNA (sgRNA) to target a CRISPR-associated (Cas) endonuclease to specific genomic DNA loci, where the Cas endonuclease generates a double-strand break. The repair of double-strand breaks by intrinsic error-prone mechanisms leads to insertions and/or deletions, disrupting the locus. Alternatively, the inclusion of double-stranded DNA donors or single-stranded DNA oligonucleotides in this process can elicit the inclusion of precise genome edits ranging from single nucleotide polymorphisms to small immunological tags or even large fluorescent protein constructs. However, a major bottleneck in this procedure can be finding and isolating the desired edit in the germline. This protocol outlines a robust method for screening and isolating germline mutations at specific loci in Danio rerio (zebrafish); however, these principles may be adaptable in any model where in vivo sperm collection is possible.

ABSTRACT Intellectual disability is a neurodevelopmental disorder that affects 2-3% of the general population. Syndromic forms of intellectual disability frequently have a genetic basis and are often accompanied by additional developmental anomalies. Pathogenic variants in components of TATA-binding protein associated factors (TAFs) have recently been identified in a subset of patients with intellectual disability, craniofacial hypoplasia, and congenital heart disease. This syndrome has been termed as a TAFopathy and includes mutations in TATA binding protein (TBP), TAF1, TAF2, and TAF6. The underlying mechanism by which TAFopathies give rise to neurodevelopmental, craniofacial, and cardiac abnormalities remains to be defined. Through a forward genetic screen in zebrafish, we have recovered a recessive mutant phenotype characterized by craniofacial hypoplasia, ventricular hypoplasia, heart failure at 96 h post-fertilization and lethality, and show it is caused by a nonsense mutation in taf5. CRISPR/CAS9 mediated gene editing revealed that these defects where phenocopied by mutations in taf1 and taf5. Mechanistically, taf5-/- zebrafish displayed misregulation in metabolic gene expression and metabolism as evidenced by RNA sequencing, respiration assays, and metabolite studies. Collectively, these findings suggest that the TAF complex may contribute to neurologic, craniofacial, and cardiac development through regulation of metabolism.