The identification of PAX7 variants and a potential role of muscle development dysfunction in congenital scoliosis

© The Author(s) 2022. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp:// creat iveco mmons. org/ licen ses/ by/4. 0/. The Creative Commons Public Domain Dedication waiver (http:// creat iveco mmons. org/ publi cdoma in/ zero/1. 0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Dear Editor,


Dear Editor,
Congenital scoliosis (CS) is a spinal malformation characterized by failure of vertebral formation or segmentation, or a mix of these deformities, resulting in longitudinal and rotational imbalance, and affects 0.05-0.1% of newborns (Wu et al. 2015). It is generally understood that the development of CS has an underlying genetic basis. Specifically, genes related to somite regulation or osteogenesis during embryonic development are believed to be responsible for the vertebral malformations observed in CS patients (Pourquié 2011).
PAX7 is a pioneer transcription factor that controls cell fate specification through chromatin remodeling. It is critical to regulate satellite cell expansion and differentiation during both neonatal and adult myogenesis (Mayran et al. 2018). Bi-allelic variants in PAX7 have been shown to cause progressive congenital myopathy with scoliosis (MYOSCO). Patients develop muscle weakness and atrophy in the early postnatal period and deteriorate with age (Feichtinger et al. 2019). Previous studies have reported the impact of PAX7 loss-of-function in human muscle development and related diseases. However, the role of PAX7 variants in congenital vertebral malformations, particularly CS, remains understudied.
Considering that scoliosis is within the phenotypic spectrum of the PAX7-related myopathy described above, gene variants of PAX7 are likely to play an essential role in isolated CS. Understanding this relationship has the potential to improve the diagnostic rate of such diseases and further expand our understanding of how muscle development influences vertebral development.
We identified nine individuals with CS who carried rare variants of PAX7 from a Chinese cohort of vertebral malformations in the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study (http:// disco study. org/) at Peking Union Medical College Hospital (PUMCH) (total cohort size n = 583). Exome sequencing and bioinformatics data analyses were performed, and variants in genes associated with CS were not detected (materials and methods are shown in supplementary material). We identified six different disease-causing candidate variants in the PAX7 gene (NM_002584.2) based on each variant's predicted deleterious functional impact, including 5 heterozygous missense variants and a heterozygous stop-gain variant locating on exons 6 and 8 of PAX7 isoforms 1 and 2 (    Figure 2 in supplementary material). The most frequent location of vertebral malformation was the thoracic spine, followed by the lumbar and cervical spine. Failure of segmentation led to rib fusion and pigeon chest in one patient, while hemivertebrae of thoracic and lower lumbar spine resulted in the absence of ribs and sacrum deformity, respectively, in two patients. This study provides evidence that heterozygous missense and stop-gain variants in the PAX7 gene may be a novel genetic cause of CS or congenital vertebral malformation. These variants have high pathogenicity prediction scores and are previously unreported.
In previously published reports, phenotypes associated with PAX7 mutations varied from cleft lip to hypotonia, muscle weakness, MYOSCO, and postnatal growth retardation (Summarized in Table 2 in supplementary material). The predominant phenotype of our patients was congenital scoliosis, characterized by either failure of vertebral formation, segmentation, or mixed deformity, expanding the phenotypic spectrum of PAX7-related diseases. Intraspinal anomalies including diastematomyelia and syringomyelia were found in 5 individuals through spinal MRI. There was no evidence of significant muscle tonus decrease or craniofacial malformation in our cohort. Possible mechanisms underlying the phenotypic differences between our group of patients and those described previously include 1) variants in previous studies are mainly located near or within the paired box domain, which is highly conserved across species; 2) during early stages of embryo development, PAX7 expression in neural crest cells (NCCs) displays a migratory pattern. Therefore its functions in the formation of cartilage and bone of the craniofacial skeleton and neural system development may result from the variation of phenotypes from cleft lip/palate to intraspinal anomalies; 3) different variants may affect the level of PAX7 expression to different extents. The PAX7 variants in our cohort were generally localized on exon 6 and 8, surrounding the paired box protein domain (Shown in Figure 1 in supplementary material), while the phenotypes varied. There may not be a strong correlation between PAX7 variants and CS classification since the patients sharing the same mutation sites or location within exons presented differently in vertebral abnormities regarding formation or segmentation failure. Further, we did not characterize the associations between distinct PAX7 gene variants and specific sub-phenotypes, including the location of vertebral anomalies, intraspinal anomalies, and extra-spinal skeletal malformations.
Paraspinal muscle plays an important role in the posture and movement of the spine. Asymmetry of the paraspinal musculature and imbalance of its strengths are likely to result in scoliosis. Differences in transcriptional profiling or dysregulated expression of genes related to the TGF-β pathway between the concave and convex side of the spinal curve were previously identified in the paravertebral muscles of adolescent idiopathic scoliosis (AIS) patients (Buchan et al. 2014). Scoliosis is a common manifestation in myopathy diseases, including Duchenne muscular dystrophy (DMD) and dystonia diseases such as Parkinson's Disease (PD). Scoliosis involving vertebral malformations has been found in patients with variants in genes related to muscle function or myogenesis but with low frequency. Zieba et al. identified an autosomal dominant form of spondylocarpotarsal synostosis (SCT) due to heterozygous missense or nonsense variants in MYH3 (Zieba et al. 2017). These patients had comorbidities of progressive scoliosis ranging from mild disc space narrowing to severe vertebral fusion. MYH3 encodes embryonic myosin heavy chain 3 and functions in muscle contraction by binding to and exerting force on actin filaments through the hydrolysis of ATP. In wild-type mice, its expression can be seen postnatally in the bone and muscle between the neural arches of the cervical and thoracic spine. MYH3 variants can also cause CS in Freeman-Sheldon Syndrome. More recently, biallelic homozygous loss-of-function variants in MYF5, a key regulator for myogenesis's early stage, were found to be associated with congenital external ophthalmoplegia, vertebral, and rib anomalies in three families (Di Gioia et al. 2018). Additionally, researchers reported kyphoscoliosis in 6 of 15 pediatric-onset cardiomyopathy patients with variants of ALPK3. These variants were predicted to have "muscle contraction" and "myogenesis" phenotypes by GeneNetwork Assisted Diagnostic Optimization (Herkert et al. 2020). Previous studies of NCC demonstrated that PAX7 staining occurs throughout the presumptive musculature surrounding the ribs and vertebrae in Carnegie Stage 18 embryos, which may interfere with the process of segmentation or vertebrae development (Betters et al. 2010). It is speculated that muscle function or myogenesis-related genes may affect vertebrae segmentation and spine growth by affecting the forces surrounding the intervertebral disc. Indeed, the paraspinal muscles, especially small multifidus muscles, connecting to the neural arches in early life stages have been implicated in causing impropriate vertebral segmentation or formation. Axial hypotonia may accelerate the progression of scoliosis with the growth of the patient's age.
In conclusion, we find that heterozygous pathogenic variants in PAX7 are associated with CS without notable dysfunction in the appendicular skeleton muscle. These findings provide important implications for diagnosis and a better understanding of the etiology of CS.
Additional file 1: Figure 1. Distribution of variant on PAX7 protein. Figure 2. Radiological characteristics of the CS patients carrying PAX7 variants. Table 2. PAX7 mutations and their associated phenotypes in previously published reports.