A Recurrent Rare SOX9 Variant (M469V) is Associated with Congenital Vertebral Malformations

Molecular landscape of congenital vertebral malformations: recent discoveries and future directions

Vertebral malformations (VMs) pose a significant global health problem, causing chronic pain and disability. Vertebral defects occur as isolated conditions or within the spectrum of various congenital disorders, such as Klippel-Feil syndrome, congenital scoliosis, spondylocostal dysostosis, sacral agenesis, and neural tube defects. Although both genetic abnormalities and environmental factors can contribute to abnormal vertebral development, our knowledge on molecular mechanisms of numerous VMs is still limited. Furthermore, there is a lack of resource that consolidates the current knowledge in this field. In this pioneering review, we provide a comprehensive analysis of the latest research on the molecular basis of VMs and the association of the VMs-related causative genes with bone developmental signaling pathways. Our study identifies 118 genes linked to VMs, with 98 genes involved in biological pathways crucial for the formation of the vertebral column. Overall, the review summarizes the current knowledge on VM genetics, and provides new insights into potential involvement of biological pathways in VM pathogenesis. We also present an overview of available data regarding the role of epigenetic and environmental factors in VMs. We identify areas where knowledge is lacking, such as precise molecular mechanisms in which specific genes contribute to the development of VMs. Finally, we propose future research avenues that could address knowledge gaps.

Variants in the SOX9 transactivation middle domain induce axial skeleton dysplasia and scoliosis

SOX9 is an essential transcriptional regulator of cartilage development and homeostasis. In humans, dysregulation of SOX9 is associated with a wide spectrum of skeletal disorders, including campomelic and acampomelic dysplasia, and scoliosis. The mechanism of how SOX9 variants contribute to the spectrum of axial skeletal disorders is not well understood. Here, we report four novel pathogenic variants of SOX9 identified in a large cohort of patients with congenital vertebral malformations. Three of these heterozygous variants are in the HMG and DIM domains, and for the first time, we report a pathogenic variant within the transactivation middle (TAM) domain of SOX9 . Probands with these variants exhibit variable skeletal dysplasia, ranging from isolated vertebral malformation to acampomelic dysplasia. We also generated a Sox9 hypomorphic mutant mouse model bearing a microdeletion within the TAM domain ( Sox9 ^Asp272del ). We demonstrated that disturbance of the TAM domain with missense mutation or microdeletion results in reduced protein stability but does not affect the transcriptional activity of SOX9. Homozygous Sox9 ^Asp272del mice exhibited axial skeletal dysplasia including kinked tails, ribcage anomalies, and scoliosis, recapitulating phenotypes observed in human, while heterozygous mutants display a milder phenotype. Analysis of primary chondrocytes and the intervertebral discs in Sox9 ^Asp272del mutant mice revealed dysregulation of a panel of genes with major contributions of the extracellular matrix, angiogenesis, and ossification-related processes. In summary, our work identified the first pathologic variant of SOX9 within the TAM domain and demonstrated that this variant is associated with reduced SOX9 protein stability. Our finding suggests that reduced SOX9 stability caused by variants in the TAM domain may be responsible for the milder forms of axial skeleton dysplasia in humans.

Environmentally relevant uptake, elimination, and metabolic changes following early embryonic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin in zebrafish

Dioxin and dioxin-like compounds are ubiquitous environmental contaminants that induce toxicity by binding to the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. The zebrafish model has been used to define the developmental toxicity observed following exposure to exogenous AHR ligands such as the potent agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (dioxin, TCDD). While the model has successfully identified cellular targets of TCDD and molecular mechanisms mediating TCDD-induced phenotypes, fundamental information such as the body burden produced by standard exposure models is still unknown. We performed targeted gas chromatography (GC) high-resolution mass spectrometry (HRMS) in tandem with non-targeted liquid chromatography (LC) HRMS to quantify TCDD uptake, model the elimination dynamics of TCDD, and determine how TCDD exposure affects the zebrafish metabolome. We found that 50 ppt, 10 ppb, and 1 ppb waterborne exposures to TCDD during early embryogenesis produced environmentally relevant body burdens: 38 ± 4.34, 26.6 ± 1.2, and 8.53 ± 0.341 pg/embryo, respectively, at 24 hours post fertilization. TCDD exposure was associated with the dysregulation of metabolic pathways that are associated with the AHR signaling pathway as well as pathways shown to be affected in mammals following TCDD exposure. In addition, we discovered that TCDD exposure affected several metabolic pathways that are critical for brain development and function including glutamate metabolism, chondroitin sulfate biosynthesis, and tyrosine metabolism. Together, these data demonstrate that existing exposure methods produce environmentally relevant body burdens of TCDD in zebrafish and provide insight into the biochemical pathways impacted by toxicant-induced AHR activation.

Pseudo-188D: Phage Protein Prediction Based on a Model of Pseudo-188D

Phages have seriously affected the biochemical systems of the world, and not only are phages related to our health, but medical treatments for many cancers and skin infections are related to phages; therefore, this paper sought to identify phage proteins. In this paper, a Pseudo-188D model was established. The digital features of the phage were extracted by PseudoKNC, an appropriate vector was selected by the AdaBoost tool, and features were extracted by 188D. Then, the extracted digital features were combined together, and finally, the viral proteins of the phage were predicted by a stochastic gradient descent algorithm. Our model effect reached 93.4853%. To verify the stability of our model, we randomly selected 80% of the downloaded data to train the model and used the remaining 20% of the data to verify the robustness of our model.

Identification of Copy Number Variants in a Southern Chinese Cohort of Patients with Congenital Scoliosis

Congenital scoliosis (CS) is a lateral curvature of the spine resulting from congenital vertebral malformations (CVMs) and affects 0.5–1/1000 live births. The copy number variant (CNV) at chromosome 16p11.2 has been implicated in CVMs and recent studies identified a compound heterozygosity of 16p11.2 microdeletion and TBX6 variant/haplotype causing CS in multiple cohorts, which explains about 5–10% of the affected cases. Here, we studied the genetic etiology of CS by analyzing CNVs in a cohort of 67 patients with congenital hemivertebrae and 125 family controls. We employed both candidate gene and family-based approaches to filter CNVs called from whole exome sequencing data. This identified 12 CNVs in four scoliosis-associated genes (TBX6, NOTCH2, DSCAM, and SNTG1) as well as eight recessive and 64 novel rare CNVs in 15 additional genes. Some candidates, such as DHX40, NBPF20, RASA2, and MYSM1, have been found to be associated with syndromes with scoliosis or implicated in bone/spine development. In particular, the MYSM1 mutant mouse showed spinal deformities. Our findings suggest that, in addition to the 16p11.2 microdeletion, other CNVs are potentially important in predisposing to CS.

Genetic variants of TBX6 and TBXT identified in patients with congenital scoliosis in Southern China

Congenital scoliosis (CS) is a spinal deformity present at birth due to underlying congenital vertebral malformation (CVM) that occurs during embryonic development. Hemivertebrae is the most common anomaly that causes CS. Recently, compound heterozygosity in TBX6 has been identified in Northern Chinese, Japanese, and European CS patient cohorts, which explains about 7%-10% of the affected population. In this report, we recruited 67 CS patients characterized with hemivertebrae in the Southern Chinese population and investigated the TBX6 variant and risk haplotype. We found that two patients with hemivertebrae in the thoracic spine and one patient with hemivertebrae in the lumbar spine carry the previously defined pathogenic TBX6 compound heterozygous variants. In addition, whole exome sequencing of patients with CS and their family members identified a de novo missense mutation (c.G47T: p.R16L) in another member of the T-box family, TBXT. This rare mutation compromised the binding of TBXT to its target sequence, leading to reduced transcriptional activity, and exhibited dominant-negative effect on wild-type TBXT. Our findings further highlight the importance of T-box family genes in the development of congenital scoliosis.

Variants Affecting the C-Terminal of CSF1R Cause Congenital Vertebral Malformation Through a Gain-of-Function Mechanism

CSF1R encodes the colony-stimulating factor 1 receptor which regulates the proliferation, differentiation, and biological activity of monocyte/macrophage lineages. Pathogenic variants in CSF1R could lead to autosomal dominant adult-onset leukoencephalopathy with axonal spheroids and pigmented glia or autosomal recessive skeletal dysplasia. In this study, we identified three heterozygous deleterious rare variants in CSF1R from a congenital vertebral malformation (CVM) cohort. All of the three variants are located within the carboxy-terminal region of CSF1R protein and could lead to an increased stability of the protein. Therefore, we established a zebrafish model overexpressing CSF1R. The zebrafish model exhibits CVM phenotypes such as hemivertebral and vertebral fusion. Furthermore, overexpression of the mutated CSF1R mRNA depleted of the carboxy-terminus led to a higher proportion of zebrafish with vertebral malformations than wild-type CSF1R mRNA did (p = 0.03452), implicating a gain-of-function effect of the C-terminal variant. In conclusion, variants affecting the C-terminal of CSF1R could cause CVM though a potential gain-of-function mechanism.

Mutational burden and potential oligogenic model of TBX6-mediated genes in congenital scoliosis

BACKGROUND

Congenital scoliosis (CS) is a spinal deformity due to vertebral malformations. Although insufficiency of TBX6 dosage contributes to a substantial proportion of CS, the molecular etiology for the majority of CS remains largely unknown. TBX6-mediated genes involved in the process of somitogenesis represent promising candidates.

METHODS

Individuals affected with CS and without a positive genetic finding were referred to this study. Proband-only exome sequencing (ES) were performed on the recruited individuals, followed by analysis of TBX6-mediated candidate genes, namely MEOX1, MEOX2, MESP2, MYOD1, MYF5, RIPPLY1, and RIPPLY2.

RESULTS

A total of 584 patients with CS of unknown molecular etiology were recruited. After ES analysis, protein-truncating variants in RIPPLY1 and MYF5 were identified from two individuals, respectively. In addition, we identified five deleterious missense variants (MYOD1, n = 4; RIPPLY2, n = 1) in TBX6-mediated genes. We observed a significant mutational burden of MYOD1 in CS (p = 0.032) compared with the in-house controls (n = 1854). Moreover, a potential oligogenic disease-causing mode was proposed based on the observed mutational co-existence of MYOD1/MEOX1 and MYOD1/RIPPLY1.

CONCLUSION

Our study characterized the mutational spectrum of TBX6-mediated genes, prioritized core candidate genes/variants, and provided insight into a potential oligogenic disease-causing mode in CS.