- 1Department of Spinal Surgery, The Affiliated Hospital of Qingdao University, Qingdao, Shandong, China
- 2Department of Health Care, The Affiliated Hospital of Qingdao University, Qingdao, Shandong, China
- 3Department of Operation Room, The Affiliated Hospital of Qingdao University, Qingdao, Shandong, China
- 4Department of Anesthesiology, The Affiliated Hospital of Qingdao University, Qingdao, Shandong, China
- 5Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
Background: To investigate the genotypic and family genetic characteristics of neurofibromatosis type 1 (NF-1) patients associated with dystrophic scoliosis and to further evaluate the clinical efficiency of surgical intervention to these patients.
Methods: A total of seven NF-1 patients with dystrophic scoliosis and their 11 immediate relatives were included in this study who visited The Affiliated Hospital of Qingdao University spinal surgery department from January 2020 to December 2022. The outcomes were summarized by investigating the clinical and imaging parameters before and after the treatment. Whole-exome sequencing (WES) was conducted to analyze the genotypic and family genetic characteristics of all patients and their families.
Results: Among the seven patients, six patients underwent surgical treatment after follow-up. Compared to preoperative Cobb angle, the maximum postoperative correction rate was 85.3%. We identified eight pathogenic variants in the NF-1 gene. The variants c.c4084T and c. T4445C were located in the GAP-related domain, and the variant c. G1885A was shared by patient 3 and his diseased sibling. In the family of patient 3, variant c. G1885A was detected in both neurofibromatosis patients. The rs112819846 exhibited the most pronounced frequency disparity, whereas rs2916067 and rs80221306 were found in all patients.
Conclusion: In this cohort study, NF-1 patients with dystrophic scoliosis predominantly presented with upper thoracic involvement, and surgical correction achieved a maximum postoperative Cobb angle correction rate of 85.3% without loss of correction during follow-up. The c. G1885A variant in the NF-1 gene may influence the phenotypic severity of the disease, while rs112819846, rs2916067, and rs80221306 may represent disease-relevant pathogenic variants of NF-1 patients with dystrophic scoliosis.
Introduction
Neurofibromatosis is an autosomal dominant genetic disease with familial inheritance (1), with neurofibromatosis type 1 (NF-1) being the most common form of neurofibromatosis. NF-1primarily affects the peripheral nervous system (2), and often involves bones, skin, and soft tissues with an incidence of 10–60% (3). Approximately 50% of NF-1 patients are susceptible to bone malnutrition, with scoliosis being the highly prevalent bone damage (4). NF-1 scoliosis can be classified as dystrophic and non-dystrophic, with significant differences of natural history and treatment plans. Non-malnourished scoliosis, which represents the more common spinal deformity of NF-1 patients, resembles adolescent idiopathic scoliosis in terms of imaging manifestations, treatment plans, and complications (5). If the Cobb angle of the main curvature of this type of scoliosis exceeds 40°, posterior spinal fusion is required. Dystrophic scoliosis, the more severe subtype, usually presents with early onset, rapid progression, and distinctive radiographic features, including vertebral scalloping and rib penciling. Given its aggressive nature, timely surgical intervention and long-term management are essential to prevent severe deformity and related complications.
Previous studies have shown that the main cause of NF-1 is related to mutations of NF-1 gene, leading to abnormalities in the encoding of neurofibromatosis proteins (6). The major reason for bone dysplasia in NF-1 is attributed to the reduced expression of neurofibromatosis proteins in osteoblasts, osteoclasts, chondrocytes, fibroblasts, and vascular endothelial cells (7). The NF-1 gene is located on chromosome 17 (17q11.2), containing approximately 60 exons and spanning over 300 bases (8). Mutations types include nonsense, missense, insertion or deletion (frameshift), splicing, and complete gene deletion. Neurofibromatosis proteins are the gene product of NF-1, containing functional protein of 2,818 amino acids with molecular weight of 280kD. Neurofibromatosis proteins have homology between the 25–40 kDa fragment of NF-1GRD and the 250–400aa fragment of p120GAP, and act as GTPase-activating proteins (GAPs) on Ras to perform their primary function (9). GAPs are currently the only clear functional fragments of neurofibromatosis proteins. NF-1 mRNA is widely expressed in neural crest-derived tissues, bones, and vascular endothelium, and its protein product, neurofibromin regulates cell proliferation and differentiation mainly through the Ras/MAPK and PI3K/AKT signaling pathways. Moreover, specific NF-1 variants have been correlated with distinct clinical phenotypes, suggesting important genotype–phenotype associations in disease severity and manifestations.
Unlike adolescent idiopathic scoliosis, malnourished scoliosis is often accompanied by significant spinal abnormalities, such as acute angular deformity, vertebral rotation or subluxation, shortened ribs, and slender pedicle (10). The typical site is the upper thoracic spine, which is often accompanied by obvious shoulder imbalance. Early and active surgical intervention is necessary to prevent the progression of scoliosis in patients with neurofibromatosis. While previous studies have conducted comparative analyses of the surgical treatment, imaging manifestations, and other aspects of scoliosis in neurofibromatosis (11–14), there remains limited research on the intrinsic genetic initiating mutations and genetics of scoliosis in neurofibromatosis. Identification of specific pathogenic mutations and genetic characteristics in neurofibromatosis patients is of great significance for the prevention of scoliosis.
Whole-exome sequencing (WES) allows for the sequencing of more than 4,000 genes that currently have known biological functions. Exons, which comprise about 1.5% of the genome, account for approximately 85% of the pathogenic variants (15, 16). WES exhibit several benefits over other genomic techniques, including lower cost, higher throughput, and greater diagnostic potential. This method has wide coverage and can comprehensively screen for pathogenic genes related to skeletal developmental abnormalities. Research has shown that WES technology improves the diagnostic rate of fetal skeletal system diseases (17, 18). With decrease in sequencing costs, its clinical applications have become increasingly widespread in the recent years.
Therefore, this study conducted WES of neurofibromatosis patients with scoliosis, combined with relevant imaging data, to explore specific gene mutations and provide a theoretical basis for the prevention and treatment of scoliosis.
Materials and methods
Clinical specimen acquisition
This study enrolled seven NF-1 patients with dystrophic scoliosis who visited The Affiliated Hospital of Qingdao University spinal surgery department from January 2020 to December 2022. Two independent doctors diagnosed NF-1-associated dystrophic scoliosis based on established criteria. The diagnostic criteria were referenced as follows: (1) six or more café-au-lait spots ≥5 mm in diameter before puberty or ≥1.5 mm in diameter after puberty, (2) axillary or inguinal skinfold freckling, (3) two or more dermal neurofibromas or one plexiform neurofibroma, (4) two or more iris hamartomas (Lisch nodules), (5) an optic pathway glioma, (6) distinctive long bone dysplasia involving the sphenoid wing or thinning of the long bone cortex with or without pseudarthrosis, and (7) a first-degree relative with NF-1. NF-1 can be diagnosed if an individual presents with two or more of these symptoms (19). All patients who met the eligibility criteria were invited to participate in the study, and written consent was obtained from each participant to publish clinical and genetic data. The clinical data of all patients was also recorded. In addition, venous blood samples of patients and family members also collected for whole-genome sequencing to identify specific genes.
DNA extraction
Blood samples were collected in ethylenediaminetetraacetic acid (EDTA) tubes and centrifuge them at 1,600 × g for 10 min at room temperature to separate the plasma. We used the MGIEasy Circulating DNA Isolation Kit (BGI) according to the manufacturer’s instructions. A total of 300 μL plasma samples were extracted by magnetic bead method. The extracted DNA was quantified using the Qubit Fluorometer (Invitrogen) and assessed for degradation and contamination via 0.8% agarose gel electrophoresis.
Library preparation and sequencing
About 0.7 μg of total gDNA of each sample was firstly sonicated to 300–500 bp size range using the Covaris LE220 (Covaris). DNA fragments were then end-repaired, A-tailed, and purified. WES library preparation was performed with the MGIEasy Exome Capture V4 Probe following the manufacturer’s instructions. Pooled libraries were sequenced on the DNBSEQ-T1&T5 (PE 100) platform (BGI-Shenzhen), with a mean sequence coverage depth of 7 × for each tumor and NAT sample, respectively.
Data preprocessing and alignment
FastQC v0.11.3 was conducted to qualify the raw sequencing reads. The low-quality reads were trimmed as clean reads using the fast (v0.19.6) software with default parameters. The remaining high-quality clean reads were mapped to the human reference genome (GRCh38) using BWA-MEM v0.7.17 and duplicate reads were marked using GATK MarkDuplicate (GATK v4.2.2.0). After alignment and trimming, all BAM files were sorted and indexed using SAMtools v1.15 for downstream processing analyses, including somatic copy number alteration (SCNA) and ecDNA analysis, somatic mutation, and structural variation detection.
Somatic mutation calling, filtering, and annotation
Somatic mutations, including single nucleotide variants (SNVs) and insertions/deletions (Indels), were identified via local assembly of haplotypes by Mutect2, with the patient and matched normal BAMs as inputs (20). To ensure data reliability, we applied multiple quality control measures, including sequencing depth ≥100 × for targeted regions, minimum base quality score ≥30, mapping quality ≥50, and variant allele frequency (VAF) threshold ≥0.05. FilterMutectCalls in GATK (v4.2.5.0) was used to filter the raw output of Mutect2, removing sequencing artifacts and low-confidence calls. Candidate variants were further evaluated by excluding common polymorphisms using population databases (gnomAD, 1,000 Genomes, ExAC) and annotated using ANNOVAR. In addition, we manually inspected variants with Integrative Genomics Viewer (IGV) to confirm their authenticity. This systematic pipeline, integrating variant calling, multi-level filtering, and functional annotation, has been widely validated in previous studies (21, 22).
Mutational signature analysis
We excluded non-synonymous variants that were shared by all family members. We then removed variants with a population frequency above 0.01 (23). Lastly, we annotated the variants with their clinical significance and identified the pathogenic variants among all family members using the ClinVar database, which is a public repository that collects and disseminates data on genomic variants, their clinical relevance, and supporting evidence (24). We assigned 18 participants to either the patient group (seven probands and one relative) or the control group (10 non-diseased family members) according to their neurofibromatosis status. Differentially pathogenic variants between the patient and control groups were identified using a one-sided Fisher’s exact test, with statistical significance defined as p < 0.05. This test was selected because of the small sample size and categorical nature of the variant frequency data.
Results
Baseline characteristics of the participants
Three male patients and four female NF-1 patients with dystrophic scoliosis were enrolled in this study (Table 1). The range of age was from 10 years to 17 years. A total of six patients’ immediate family also underwent venous blood sampling for comprehensively analyzing the specific pathogenic genes caused by neurofibromatosis. Abnormalities of thoracic vertebrae were detected in five of seven patients. Among the seven patients, three had thoracic scoliosis (including one with kyphotic deformity), two had lumbar scoliosis, and two had thoracolumbar involvement. The apical vertebrae were mainly located in the upper thoracic region (T4–T12), and two patients presented with shoulder imbalance. One patient (patient 3) showed complex vertebral anomalies, including intravertebral rib head, hemivertebra, and abnormal pedicles. Family history was positive in two patients, with affected first-degree relatives.
Among the seven patients, the shortest and longest follow-up time are 6 months and 2 years, respectively, since the diagnosis of scoliosis. Ultimately, six patients underwent surgical treatment after follow-up and there were no related complications such as increased Cobb angle, failed internal fixation, or neurological dysfunction.
Seven neurofibromatosis patients underwent regular X-ray from the first diagnosis to determine the progression of scoliosis. Among the six surgical patients, the minimum and maximum Cobb angle were 36.9° and 68°, respectively, before surgery. Compared to preoperative Cobb angle, the maximum postoperative correction rate was 85.3%. One patient had abnormal vertebral development, including small pedicle, rib head invasion into the vertebral canal, and hemivertebra. Among the six patients who underwent surgery, the longest fixed segment spanned 10 vertebral levels (T6-L4). Postoperative follow-up of all patients showed no loss of correction (Table 1).
Identification of pathogenic variants of NF-1 gene
We performed WES on seven patients and their 11 immediate relatives, and identified eight pathogenic variants in the NF-1 gene after filtering and annotation. Of these, seven variants were present in neurofibromatosis patients (five patients and sibling of patient 3). Only c.8413A > T was unique to unaffected family members. All of these identified variants were previously reported in public databases and are not novel. We examined the association between amino acid positions altered by NF-1 variants in patients and peptide chain domains (Figure 1). We discovered that variants c.4084C > T and c.4445 T > C are located in the GAP-related domain (GRD), which is a well-studied domain. We also observed that variant c.1885G > A, which was present in the patient exhibiting abnormal vertebral morphology, resided in the cysteine- and serine-rich domain (CSRD) (Table 2).

Figure 1. Schematic representation of distribution of variants in NF-1 and neurofibromin domains. The amino acid number is indicated below. CSRD, cysteine- and serine-rich domain; TBD, tubulin-binding domain; GRD, GAP-related domain; Sec-PH, Sec14 homologous and pleckstrin homologous domain; CTD, C-terminal domain. The amino acid number is indicated below.
Family-based analysis
Family-based analysis was performed in the families of four patients with NF-1 gene variants (Figure 2). Patient 1 inherited c.4382 T > C from her father. Moreover, due to the unavailability of the mother’s sample, the origin of variant c.1278G > A remains undetermined. In the family of patient 3, variant c.1885G > A was detected in both neurofibromatosis patients, suggesting that it may influence the penetrance of neurofibromatosis. In the family of patient 4, variant c.7043G > A was present in the patient and her non-diseased mother. Finally, as reported in patient 5’s family, c.5291delC was associated with the penetrance of NF-1.

Figure 2. Family-based analysis of the genetic variants of the NF-1 gene in four probands and their families. Variants are represented by blocks of different colors. The pedigree symbols are as follows: square for male, circle for female, arrow for proband, “?” for specimen unavailable, and shaded for NF-1 diseased.
Identification of disease-related variants
We performed a one-sided analysis of the pathogenic variants that are associated with the disease phenotype in the patient and control groups, and we detected five disease-related variants with statistically significant frequency disparities. All these variants were loci with high mutation rates that are cataloged in the dbSNP database. Among them, rs112819846 exhibited the most pronounced frequency disparity (p = 0.00653), whereas rs2916067 (p = 0.01282), and rs80221306 (p = 0.0359) were found in all the patients with NF-1 (Table 3).
Discussion
In 1987, the National Research Institute of the United States proposed seven diagnostic criteria for NF-1, including skin milky coffee spots, neurofibromas, inguinal freckles, optic gliomas, Lisch nodules, characteristic skeletal change, and family history (25). Due to some diagnostic features changing with age, the international NF-1 diagnostic standard has some limitations, which may lead to delayed diagnosis in some patients. Spinal scoliosis, as the most common bone system lesion in NF-1, can be diagnosed on changes in patient’s physical appearance clinical examination, typically presenting during the first peak period between 5 and 10 years of age (26). According to previous studies, the incidence of NF spinal deformities is 10–64% with the faster progression rate of deformities than those of AIS (27). Therefore, scoliosis screening should be performed in all pre-adolescent children with NF-1. Due to the significant differences in treatment and prognosis between the two types of scoliosis associated with NF-1, it is vital to make timely diagnosis on the types of scoliosis.
Patients with NF-1 accompanied by scoliosis can be divided into malnourished and non-malnourished based on imaging features. Non-malnourished scoliosis typically presents with long, smooth curves. In contrast, malnourished scoliosis, compared with AIS, typically involves four to six vertebral bodies and is characterized by short and sharp lateral curvature, with or without angular kyphosis. In addition, the undernourished scoliosis, which often develops gradually in young age, is accompanied by pencil-shaped changes of ribs, vertebral rotation, vertebral wedge-shaped changes, and lateral curvature of atypical positions (28). Li et al. measured 2,652 pedicle abnormalities of 56 NF-1 patients with dystrophic scoliosis and 22 NF-1 patients with non-dystrophic scoliosis (29). The results showed that the overall incidence of pedicle abnormalities in NF-1 patients was as high as 67%, and the upper thoracic spine (87%) was the most concentrated area of abnormal pedicle abnormalities. In addition to scoliosis, patients with neurofibromatosis often have abnormalities in the skin, blood vessels, eyes, and nervous system. Approximately 20–50% of NF-1 patients have osteoporosis, which is associated with alterations in bone transformation biochemical markers, including decreased serum 25-hydroxyvitamin D and elevated serum PTH levels (30). In this study, one patients also experienced changes in multiple imaging indicators, such as intraspinal rib head(T9, T10),a hemivertebra(T9,)and abnormal pedicles(T9, T10). The accurate recognition of these radiological features is beneficial for early diagnosis of patients with malnourished scoliosis to develop timely and appropriate treatment plans.
With the development of modern diagnostic technology, NF-1 can be diagnosed early through genetic testing. In 2021, the International Consensus Group of Diagnostic Standards for Neurofibromatosis proposed revisions of the NF-1 diagnostic standard (1987) mainly incorporating genetic diagnosis (31). Genetic testing was able to confirm diagnosis before clinical diagnosis (or before the appearance of the second clinical feature), which could help screen out children with atypical features. The NF-1 gene has been recognized as the pathogenic factor for NF-1 in the medical community. However, the correlation between different mutation types and sites of NF-1 and the clinical phenotype of NF-1 patients has always been a hot research topic in the medical community. Although scoliosis is the most common type of skeletal changes in NF-1, there is limited research on the correlation between specific mutation types and mutation sites in NF-1 patients with scoliosis.
In this study, we conducted WES on 18 individuals from seven probands and their immediate relatives to explore gene variants associated with NF-1. As previous studies have shown that mutations in the NF-1 gene lead to the occurrence of NF-1, we performed family-based analysis to identify the variants of NF-1 gene associated with phenotypes and characterized their effects on the amino acid sequence and structure (32). The NF-1 gene codes for neurofibromin, a protein of 2,818 amino acids that negatively regulates the RAS pathway, a signaling cascade involved in cell growth and differentiation (33). Neurofibromin acts as a GAP for Ras, a family of small GTP-binding proteins that are key regulators of the RAS pathway. A 360 amino-acid segment of neurofibromin, called the GAP-related domain (GRD), enhances the GTP hydrolysis of Ras, thereby inactivating it (34). NF-1 mutations impair the function of neurofibromin, resulting in the accumulation of active Ras and the hyperactivation of the RAS pathway, which underlies the tumorigenesis and other phenotypes of NF-1. We identified two mutations located in the GRD region in patients 1(c.4445 T > C, p. Ile1482Thr) and 6(c.4084C > T, p. Arg1362Ter). c.4445 T > C is a single nucleotide variation located at chr17:31260383 with a minor allele frequency of 0.000024, resulting in a missense mutation of isoleucine to threonine at codon 1,482. This SNV is associated with not only NF-1 (NF-1) but also hereditary cancer predisposition syndrome. c.4084C > T is located at chr17:31249093 with a minor allele frequency of 0.000008, resulting in a premature stop codon in the mRNA transcript and leading to a truncated protein product. Eisenbarth et al. (35) identified c.4084C > T as a tumor-specific point mutation in seven neurofibromas from four different NF-1 patients and demonstrated that loss of function of NF-1 gene is a common mechanism in tumorigenesis. Hence, variants in the GRD may impair the function of neurofibromin and contribute to the genetic predisposition for neurofibromatosis.
In the family-based analysis, we identified c.1885G > A in the NF-1 gene in patient 3 and his sibling, both of whom exhibited NF-1 phenotypes, and this variant was absent in their normal parents. Furthermore, patient 3 uniquely presented with vertebral body dysplasia, indicating that c.1885G > A is a genetic factor that contributes to the development of NF-1. c.1885G > A is an SNV at chr17:31249093 with a minor allele frequency of 0.000024, resulting in a nonconservative amino acid substitution in the CSRD of neurofibromin. The CSRD of neurofibromin undergoes phosphorylation by protein kinase A (PKA) and protein kinase C (PKC). The PKC-mediated phosphorylation of the CSRD modulates the conformation and function of neurofibromin, a Ras-GAP, and promotes its interaction with actin. In our cohort study, none of the studied families had consanguineous marriage, suggesting that the identified variants were not attributable to consanguinity.
Through a one-sided analysis between the patient and control group, we detected five NF-1 related pathogenic variants with statistically significant frequency. We observed that two variants, rs2916067 and rs80221306, were present in all the individuals with NF-1 in our study. Variant rs2916067 is an SNV in the coding region of MAGEF1, a ubiquitin ligase enhancer and a member of the melanoma-associated antigen (MAGE) super family (36). MAGEA2 is overexpressed in various cancers and suppresses the RAS pathway. Correspondingly, loss-of-function mutations in the NF-1 gene inhibit the RAS pathway, impair bone development, and lead to skeletal deformities such as scoliosis and pseudarthrosis. Variant rs80221306 is located at gene filaggrin (FLG), which affects differentiation of keratinocytes (37). The dysregulation of the keratinocyte–melanocyte axis results in the formation of café-au-lait macules and malignant melanoma in NF-1, but the underlying molecular basis of this dysregulation remains elusive. The FLG gene might play a role in facilitating phenotypic expression.
There are still some limitations that need to be addressed. First, the small sample size of this study reduces the statistical power. However, the NF-1 is quite rare condition, and only a few of patients present with the dystrophic scoliosis. In addition, WES is an expensive test for patients, which limit the sample size. Second, all participants in the present study were from a single ethnic group, which limits the generalizability of the findings to the broader human population. Third, the study would have been more robust if genetically healthy humans were considered as the control group. Another limitation of this study is that we did not perform functional assays to validate the biological effects of the identified NF-1 variants, and further experimental studies are required to confirm their pathogenicity.
Conclusion
In our cohort study, NF-1 patients with dystrophic scoliosis predominantly exhibited upper thoracic involvement, and surgical correction achieved satisfactory outcomes, with a maximum postoperative Cobb angle correction rate of 85.3% and no loss of correction during follow-up. The c. G1885A variant in the NF-1 gene may influence the phenotypic severity of the disease, and rs112819846, rs2916067, and rs80221306 may be the disease-related pathogenic variants of NF-1 patients with dystrophic scoliosis.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the ethics committee of the Affiliated Hospital of Qingdao University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YD: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. TB: Formal analysis, Investigation, Writing – original draft. JS: Investigation, Methodology, Writing – review & editing. LZ: Data curation, Software, Writing – original draft. HH: Methodology, Project administration, Writing – original draft. PH: Data curation, Methodology, Writing – original draft. FG: Data curation, Methodology, Writing – original draft. JG: Formal analysis, Investigation, Methodology, Writing – original draft. JL: Formal analysis, Methodology, Writing – original draft. CL: Data curation, Formal analysis, Investigation, Writing – original draft. JS: Data curation, Writing – review & editing. GZ: Data curation, Formal analysis, Methodology, Writing – review & editing. HT: Conceptualization, Formal analysis, Methodology, Supervision, Writing – review & editing. YX: Conceptualization, Formal analysis, Funding acquisition, Investigation, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was financially supported by the Taishan Scholar Program of Shandong Province, China (tstp20250511); Qingdao Technology-Benefiting-People Program (No. 23-2-8-amjk-9-nsh) and Medicine Plus Joint Research Program of Qingdao University (YX2024303).
Acknowledgments
We thank Jiale Shao, Shuqing Chen, and Chenxu Li for the data collection of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Tamura, R. Current understanding of Neurofibromatosis type 1, 2, and Schwannomatosis. Int J Mol Sci. (2021) 22:850. doi: 10.3390/ijms22115850
2. Fisher, MJ, Blakeley, JO, Weiss, BD, Dombi, E, Ahlawat, S, Akshintala, S, et al. Management of neurofibromatosis type 1-associated plexiform neurofibromas. Neuro-Oncology. (2022) 24:1827–44. doi: 10.1093/neuonc/noac146
3. Wilson, BN, John, AM, Handler, MZ, and Schwartz, RA. Neurofibromatosis type 1: new developments in genetics and treatment. J Am Acad Dermatol. (2021) 84:1667–76. doi: 10.1016/j.jaad.2020.07.105
4. Neifert, SN, Khan, HA, Kurland, DB, Kim, NC, Yohay, K, Segal, D, et al. Management and surgical outcomes of dystrophic scoliosis in neurofibromatosis type 1: a systematic review. Neurosurg Focus. (2022) 52:E7. doi: 10.3171/2022.2.FOCUS21790
5. Xu, AL, Suresh, KV, Gomez, JA, Emans, JB, Larson, AN, Cahill, PJ, et al. Consensus-based best practice guidelines for the Management of Spinal Deformity and Associated Tumors in pediatric Neurofibromatosis type 1: screening and surveillance, surgical intervention, and medical therapy. J Pediatr Orthop. (2023) 43:e531–7. doi: 10.1097/BPO.0000000000002431
6. Li, H, Zhang, W, Yao, Z, Guo, R, Hao, C, and Zhang, X. Genotypes and clinical intervention of patients with neurofibromatosis type 1 associated dystrophic scoliosis. Front Pediatr. (2022) 10:918136. doi: 10.3389/fped.2022.918136
7. Cui, XW, Ren, JY, Gu, YH, Li, QF, and Wang, ZC. NF1, Neurofibromin and gene therapy: prospects of next-generation therapy. Curr Gene Ther. (2020) 20:100–8. doi: 10.2174/1566523220666200806111451
8. Gutmann, DH, Ferner, RE, Listernick, RH, Korf, BR, Wolters, PL, and Johnson, KJ. Neurofibromatosis type 1. Nat Rev Dis Primers. (2017) 3:17004. doi: 10.1038/nrdp.2017.4
9. de Blank, PMK, Gross, AM, Akshintala, S, Blakeley, JO, Bollag, G, Cannon, A, et al. MEK inhibitors for neurofibromatosis type 1 manifestations: clinical evidence and consensus. Neuro-Oncology. (2022) 24:1845–56. doi: 10.1093/neuonc/noac165
10. Haider, S, Le, LQ, Cho, G, Xi, Y, and Chhabra, A. Scoliosis in Neurofibromatosis type 1 on whole-body magnetic resonance imaging: frequency and association with Intraspinal and Paraspinal tumors. J Comput Assist Tomogr. (2022) 46:231–5. doi: 10.1097/RCT.0000000000001263
11. Cai, S, Tian, Y, Qiu, G, Zhang, J, Shen, J, Zhao, H, et al. Neurofibromatosis type 1 with severe dystrophic kyphosis: surgical treatment and prognostic analysis of 27 patients. Orthop Surg. (2020) 12:1923–40. doi: 10.1111/os.12848
12. Mladenov, KV, Spiro, AS, Krajewski, KL, Stücker, R, and Kunkel, P. Management of spinal deformities and tibial pseudarthrosis in children with neurofibromatosis type 1 (NF-1). Childs Nerv Syst. (2020) 36:2409–25. doi: 10.1007/s00381-020-04775-4
13. Yao, Z, Li, H, Zhang, X, Li, C, and Qi, X. Incidence and risk factors for instrumentation-related complications after scoliosis surgery in pediatric patients with NF-1. Spine. (2018) 43:1719–24. doi: 10.1097/BRS.0000000000002720
14. Zhao, J, Meng, Y, Ma, J, Zhou, X, and Jiang, H. Sectional correction technique in dystrophic scoliosis secondary to Neurofibromatosis type 1: a comparison with traditional 2-rod correction technique. World Neurosurg. (2022) 167:e507–14. doi: 10.1016/j.wneu.2022.08.050
15. Jelin, AC, and Vora, N. Whole exome sequencing: applications in prenatal genetics. Obstet Gynecol Clin N Am. (2018) 45:69–81. doi: 10.1016/j.ogc.2017.10.003
16. Marwaha, S, Knowles, JW, and Ashley, EA. A guide for the diagnosis of rare and undiagnosed disease: beyond the exome. Genome Med. (2022) 14:23. doi: 10.1186/s13073-022-01026-w
17. Peng, Y, Yang, S, Huang, X, Pang, J, Liu, J, Hu, J, et al. Whole exome sequencing analysis in fetal skeletal dysplasia detected by ultrasonography: an analysis of 38 cases. Front Genet. (2021) 12:728544. doi: 10.3389/fgene.2021.728544
18. Tang, H, Zhang, Q, Xiang, J, Yin, L, Wang, J, and Wang, T. Whole exome sequencing aids the diagnosis of fetal skeletal dysplasia. Front Genet. (2021) 12:599863. doi: 10.3389/fgene.2021.599863
19. Margraf, RL, VanSant-Webb, C, Mao, R, Viskochil, DH, Carey, J, Hanson, H, et al. NF1 somatic mutation in dystrophic scoliosis. J Mol Neurosci. (2019) 68:11–8. doi: 10.1007/s12031-019-01277-0
20. Han, P, Wei, G, Cai, K, Xiang, X, Deng, WP, Li, YB, et al. Identification and functional characterization of mutations in LPL gene causing severe hypertriglyceridaemia and acute pancreatitis. J Cell Mol Med. (2020) 24:1286–99. doi: 10.1111/jcmm.14768
21. Zhang, R, Chen, S, Han, P, Chen, F, Kuang, S, Meng, Z, et al. Whole exome sequencing identified a homozygous novel variant in CEP290 gene causes Meckel syndrome. J Cell Mol Med. (2020) 24:1906–16. doi: 10.1111/jcmm.14887
22. Dai, Y, Liang, S, Dong, X, Zhao, Y, Ren, H, Guan, Y, et al. Whole exome sequencing identified a novel DAG1 mutation in a patient with rare, mild and late age of onset muscular dystrophy-dystroglycanopathy. J Cell Mol Med. (2019) 23:811–8. doi: 10.1111/jcmm.13979
23. Wu, Z, Jiang, Z, Li, T, Xie, C, Zhao, L, Yang, J, et al. Structural variants in the Chinese population and their impact on phenotypes, diseases and population adaptation. Nat Commun. (2021) 12:6501. doi: 10.1038/s41467-021-26856-x
24. Landrum, MJ, Lee, JM, Benson, M, Brown, GR, Chao, C, Chitipiralla, S, et al. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. (2018) 46:D1062–d1067. doi: 10.1093/nar/gkx1153
25. Korf BR. Neurofibromatosis. Handb Clin Neurol. (2013) 111:333–40. doi: 10.1016/B978-0-444-52891-9.00039-7
26. Pierce, KE, Passias, PG, Alas, H, Brown, AE, Bortz, CA, Lafage, R, et al. Does patient frailty status influence recovery following spinal fusion for adult spinal deformity?: an analysis of patients with 3-year follow-up. Spine. (2020) 45:E397–e405. doi: 10.1097/BRS.0000000000003288
27. Cai, S, Cui, L, Qiu, G, Shen, J, and Zhang, J. Comparison between surgical fusion and the growing-rod technique for early-onset neurofibromatosis type-1 dystrophic scoliosis. BMC Musculoskelet Disord. (2020) 21:455. doi: 10.1186/s12891-020-03460-6
28. Tauchi, R, Kawakami, N, Suzuki, T, Uno, K, Yamaguchi, T, Yanagida, H, et al. Comparison of early definitive fusion and traditional growing rods in early-onset dystrophic scoliosis in Neurofibromatosis type 1: a preliminary report. J Pediatr Orthop. (2020) 40:569–74. doi: 10.1097/BPO.0000000000001579
29. Li, Y, Luo, M, Wang, W, Shen, M, Xu, G, Gao, J, et al. A computed tomography-based comparison of abnormal vertebrae pedicles between dystrophic and nondystrophic scoliosis in Neurofibromatosis type 1. World Neurosurg. (2017) 106:898–904. doi: 10.1016/j.wneu.2017.07.064
30. Gonca Kaçar, A, Kılınc Oktay, B, Çınar Özel, S, Ocak, S, Güneş, N, Uludağ Alkaya, D, et al. Neurofibromatosis type 1 in children: a single-center experience. Turk Arch Pediatr. (2021) 56:339–43. doi: 10.5152/TurkArchPediatr.2021.20165
31. Legius, E, Messiaen, L, Wolkenstein, P, Pancza, P, Avery, RA, Berman, Y, et al. Revised diagnostic criteria for neurofibromatosis type 1 and Legius syndrome: an international consensus recommendation. Genet Med. (2021) 23:1506–13. doi: 10.1038/s41436-021-01170-5
32. Monroe, CL, Dahiya, S, and Gutmann, DH. Dissecting clinical heterogeneity in Neurofibromatosis type 1. Annu Rev Pathol. (2017) 12:53–74. doi: 10.1146/annurev-pathol-052016-100228
33. Bergoug, M, Doudeau, M, Godin, F, Mosrin, C, Vallée, B, and Bénédetti, H. Neurofibromin structure, functions and regulation. Cells. (2020) 9:9. doi: 10.3390/cells9112365
34. Sherekar, M, Han, SW, Ghirlando, R, Messing, S, Drew, M, Rabara, D, et al. Biochemical and structural analyses reveal that the tumor suppressor neurofibromin (NF1) forms a high-affinity dimer. J Biol Chem. (2020) 295:1105–19. doi: 10.1016/S0021-9258(17)49919-4
35. Eisenbarth, I, Beyer, K, Krone, W, and Assum, G. Toward a survey of somatic mutation of the NF1 gene in benign neurofibromas of patients with neurofibromatosis type 1. Am J Hum Genet. (2000) 66:393–401. doi: 10.1086/302747
36. Weon, JL, Yang, SW, and Potts, PR. Cytosolic Iron-sulfur assembly is evolutionarily tuned by a Cancer-amplified ubiquitin ligase. Mol Cell. (2018) 69:113–125.e6. doi: 10.1016/j.molcel.2017.11.010
Keywords: neurofibromatosis type 1, dystrophic scoliosis, whole-exome sequencing, gene variants, surgical correction, family genetic characteristics
Citation: Du Y, Bai T, Song J, Zhang L, Huang H, Han P, Gai F, Guo J, Li J, Lv C, Shao J, Zhang G, Tao H and Xi Y (2025) The genotypic and family characteristics and clinical intervention of neurofibromatosis type 1 gene are associated with dystrophic scoliosis by whole-exome sequencing. Front. Neurol. 16:1641665. doi: 10.3389/fneur.2025.1641665
Edited by:
Félix Javier Jiménez-Jiménez, Hospital Universitario del Sureste, SpainReviewed by:
Seyed Abbas Rafat, University of Tabriz, IranSantasree Banerjee, Jilin University, China
Copyright © 2025 Du, Bai, Song, Zhang, Huang, Han, Gai, Guo, Li, Lv, Shao, Zhang, Tao and Xi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guodong Zhang, ZHJ6aGFuZzg5MDYyNUAxNjMuY29t; Hao Tao, dGFvaGFvYjIwMDhAMTYzLmNvbQ==; Yongming Xi, eHltNzAwMTE4QDE2My5jb20=
†These authors have contributed equally to this work