Introduction: Neurofibromatosis type 1 (NF-1) is caused by mutations in the NF1 gene that encodes neurofibromin, a negative regulator of RAS proto-oncogene. Approximately one-third of the reported pathogenic mutations in NF1 are splicing mutations, but most consequences are unclear. The objective of this study was to identify the pathogenicity of splicing mutation in a Chinese family with NF-1 and determine the effects of the pre-mRNA splicing mutation by in vitro functional analysis. Methods: Next-generation sequencing was used to screen candidate mutations. We performed a minigene splicing assay to determine the effect of the splicing mutation on NF1 expression, and three-dimensional structure models of neurofibromin were generated using SWISS-MODEL and PROCHECK methods, respectively. Results: A pathogenic splicing mutation c.479+1G>C in NF1 was found in the proband characterized by childhood-onset refractory hypertension. In vitro analysis demonstrated that c.479+1G>C mutation caused the skipping of exon 4, leading to a glutamine-to-valine substitution at position 97 in neurofibromin and an open reading frame shift terminating at codon 108. Protein modeling showed that several major domains were missing in the truncated neurofibromin protein. Conclusion: The splicing mutation c.479+1G>C identified in a Chinese patient with NF-1 and childhood-onset refractory hypertension caused the skipping of exon 4 and a truncated protein. Our findings offer new evidence for the molecular diagnosis of NF-1.

Neurofibromatosis type 1 (NF-1; OMIM #162200) is one of the most common autosomal dominant heritable cancer predisposition syndromes, affecting about 1 in 2,000–3,000 individuals worldwide [1‒3]. NF-1 results from inactivating mutations in the NF1 gene that encodes neurofibromin, a large protein (approximately 250 kDa) with 2,818 amino acids that plays a role in upstream negative regulation of RAS-guanosine triphosphatases. The pathophysiology of NF1 is attributed to the loss of the negative RAS regulator neurofibromin, activating downstream cascade-signaling pathways including the MAPK/extracellular-signal-regulated kinase/mitogenic extracellular signal-regulated kinase (MEK/ERK) pathway and the phosphoinositide 3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/Akt/mTOR) pathway [4].

NF1, located at 17q11.2, is a large gene that spans approximately 350 kb of genomic DNA and contains 60 exons. Complex and alternative splicing patterns in NF1 have resulted in many isoforms of neurofibromin [5]. To date, nearly 3,890 DNA variants of NF1 are listed in the Human Gene Mutation Database (HGMD, Institute of Medical Genetics, Cardiff, UK; http://www.hgmd.org/). Approximately 50% of the reported mutations are de novo mutations, and postzygotic NF1 mutations that can lead to a mosaic phenotype have also been reported [6]. Most of the pathogenic mutations have produced truncated neurofibromin proteins, and 30% of them are splicing mutations that affect mRNA processing [7, 8]. Destruction of the 3′ or 5′ splice site caused by splicing variants results in abnormal neurofibromins including exon skipping and intron retention [9]. Although in silico tools have been used to predict the pathogenicity of splicing mutations, the results need to be confirmed by, for example, in vivo splicing analysis or in vitro minigene assay experiments.

The clinical phenotype of NF1 is progressive with complete penetrance. Café-au-lait macules (CALMs), axillary and inguinal freckling, and subcutaneous neurofibromas are common features of NF-1, and brain tumors (optic pathway gliomas and glioblastoma), skeletal dysplasia, and malignant peripheral nerve sheath tumors are less common but more serious complications [10]. Cardiovascular abnormalities, including arterial stenosis and aneurysms, are under-recognized complications that are associated with increased mortality and deserve close attention in patients with NF-1 [11]. Early diagnosis of NF-1 and regular screening of high-risk complications can help improve the prognosis and decrease the heavy disease burden. Conventionally, the clinical diagnosis of NF-1 has been based on diagnostic criteria outlined in a National Institutes of Health (NIH) Consensus Development Conference [12]. However, the variable and complex phenotypes of NF-1 make the clinical diagnosis of this disease challenging. More recently, next-generation sequencing methods have facilitated the development of molecular diagnosis approaches that have been applied to the diagnosis and management of NF-1.

Here, we identified a splicing mutation c.479+1G>C in NF1 in a Chinese pedigree that we detected using next-generation sequencing combined with Sanger sequencing. We used a minigene splicing assay to verify the consequence of the mutation, further contributing to the understanding of the pathogenesis of NF-1.

Subjects and Clinical Evaluation

The proband, an 18-year-old male, was referred to Fuwai Hospital (Beijing, China) for refractory hypertension. Previously, he was diagnosed with hypertension, left renal artery, and superior mesenteric artery stenosis and underwent a percutaneous transluminal renal angioplasty (PTRA) procedure. However, his blood pressure (BP) still fluctuated in the range 150–180/100–110 mm Hg. During hospitalization, a detailed physical examination, imaging, and laboratory tests were performed. Six of the proband’s family members were also recruited for genetic and clinical characteristics analyses (Fig. 1a).

Fig. 1.

A splicing mutation c.479+1G>C in NF1 identified in a Chinese family. a Pedigree of the family. b Sanger sequencing results of NF1 gene in the patient with neurofibromatosis type 1 and a normal subject. The black arrow indicates the proband.

Fig. 1.

A splicing mutation c.479+1G>C in NF1 identified in a Chinese family. a Pedigree of the family. b Sanger sequencing results of NF1 gene in the patient with neurofibromatosis type 1 and a normal subject. The black arrow indicates the proband.

Close modal

This study was approved by the Ethics Committee of Fuwai Hospital and carried out in accordance with the Declaration of Helsinki. Each participant or authorized representative signed written informed consent forms in advance.

Sample Collection and DNA Sequencing

Peripheral blood samples were retrieved from the enrolled subjects and genomic DNA was extracted using a QIAamp DNA Mini Kit (Qiagen, Valencia, CA, USA). Next-generation sequencing of the proband’s sample was performed using a monogenic hypertension-related-gene panel (online suppl. Table 1; for all online suppl. material, see https://doi.org/10.1159/000533144). The sequenced data were aligned to the reference human genome (GRCh37/hg19, https://www.ncbi.nlm.nih.gov/assembly/GCF_000001405.13/), then single-nucleotide variants, insertion/deletion variants, and copy number variants were identified and analyzed as described previously [13, 14]. SpliceSiteFinder-like (v3.1; http://umd.be/Redirect.html), NNSPLICE (v0.9; https://www.fruitfly.org/seq_tools/splice.html), MaxEntScan [15], Human Splicing Finder [16], and NetGene2 (https://services.healthtech.dtu.dk/service.php?NetGene2-2.42) were used to predict the effect of the genetic variant on splicing. Sanger sequencing was used to verify the possible pathogenic mutations in the family members’ samples.

Construction of NF1 Expression Plasmids

Briefly, the amplified target sequences contained exons 3–6 of human NF1 (exon 3-intron 3-exon 4-intron 4-exon 5-intron 5-exon 6) and the mutant fragment of NF1 was induced by site-directed mutagenesis (JiKai Gene, Shanghai, China). These fragments were inserted into GV219 reporter vectors with restriction endonucleases Xhol/Kpnl. Sanger sequencing was performed to verify the successful construction of the recombinant wild-type (WT) and mutant-type (MT) expression vectors (Fig. 2).

Fig. 2.

Successful construction of the WT and MT NF1 expression plasmids verified by Sanger sequencing.

Fig. 2.

Successful construction of the WT and MT NF1 expression plasmids verified by Sanger sequencing.

Close modal

Cell Culture, Transfection, and Reverse Transcription PCR

Human HEK293T cells (Cell Resource from National Center for Cardiovascular Diseases, Beijing, China) were routinely cultured in Dulbecco’s Modified Eagle’s Medium (Invitrogen, Waltham, MA, USA) containing 10% fetal bovine serum (Gibco, Carlsbad, CA, USA). Lipofectamine 3000 reagent (Invitrogen, Waltham, MA, USA) was used for transient plasmid transfection in accordance with the manufacturer’s protocol. Cells were incubated for 24 h after transfection, then harvested. Total RNA was extracted from the cells using TRIzol reagent (Invitrogen) and reverse-transcribed using PrimeScriptRT Master Mix (Takara, RR036A). Complementary DNA (cDNA) sequences were amplified using 2× PCR Master Mix (Beyotime, Shanghai, China). The products were purified by 2.4% agarose DNA gel electrophoresis and sequenced. The primers for the PCR amplification were as follows: forward, 5′-AGA​ATA​TTT​GGA​GAA​GCT​GCT​G-3′; reverse, 5′-ATA​ACT​GCT​AAC​TGC​GCA​ACC-3′.

Molecular Modeling of Neurofibromin

We conducted a homology search in NCBI database using the WT sequence of NF1 as the query sequence and found a 100% homologous template of NF1 (Protein Data Bank ID 7R03) [17]. A three-dimensional structure model of WT neurofibromin was generated using SWISS-MODEL [18], and RoseTTAFold was used to predict the three-dimensional structure of MT neurofibromin [19]. We used PROCHECK to evaluate the optimized WT and MT neurofibromin models [20].

Clinical Characteristics

The proband was born with ten CALMs that were distributed in the armpit, groin, neck; the largest spot was 4 mm in diameter after puberty. A few freckles were scattered around the armpits, groin, and neck area. At 6 years of age, the proband complained of intermittent headache, had nausea and vomiting symptoms, and was diagnosed with hypertension (BP 200/160 mm Hg). After screening the etiology of early onset hypertension, he was found to have renal artery and superior mesenteric artery stenosis. At 12 years of age, the proband underwent a PTRA procedure to reduce his BP. More recently, despite of taking metoprolol and amlodipine, his BP fluctuated in the range 150–180/100–110 mm Hg. During hospitalization, an aortic computerized tomography scan showed moderate stenosis in the proximal segment of the left renal artery and severe stenosis in the initial segment of the superior mesenteric artery (Fig. 3). The echocardiographic findings suggested left ventricular wall hypertrophy, positive motor stimulation test, left ventricular outflow tract, and left ventricular midsection obstruction. The ophthalmologic examination revealed bilateral Lisch nodules and hypertensive retinal artery stenosis. BP monitoring of extremities found BP readings of 178/97 mm Hg (right upper limb), 170/85 mm Hg (left upper limb), 181/87 mm Hg (right ankle), and 186/81 mm Hg (left ankle). The results of the biochemical examinations of the proband’s sample are list in Table 1. Mental retardation, hearing impairment, skeletal abnormalities, or subcutaneous neurofibromas were not observed in the proband. Essentially, the patient met the NIH clinical diagnostic criteria for NF-1 [12].

Fig. 3.

Left renal artery stenosis detected in the proband. The red arrow indicates the vascular lesion. a Sagittal view. b Coronal view.

Fig. 3.

Left renal artery stenosis detected in the proband. The red arrow indicates the vascular lesion. a Sagittal view. b Coronal view.

Close modal
Table 1.

Biochemical characteristics of the proband during hospitalization

VariablesProbandNormal range
Renal function test 
 Creatinine, μmmol/L 73.36 44–133 
 BUN, mmol/L 5.19 2.86–7.90 
 mALB, mg/L 100.30 <30 
 URIC, μmmol/L 419.68 148.8–416.5 
Serum electrolytes 
 Potassium, mmol/L 3.70 3.5–5.3 
 Sodium, mmol/L 138.30 137–147 
 Chloride, mmol/L 102.90 99–110 
Serum hormone concentration 
 ACTH, pg/mL 9.63 0–46 
 TSH, μIU/mL 0.96 0.55–4.78 
 Thyroxine, μg/mL 7.40 4.29–12.47 
 Aldosterone, ng/dL 7.8 3.0–23.6 
 PRC, μIU/mL 98.2 2.8–39.9 
 ADRR, (ng/dL)/(μIU/mL) 0.079 <3.7 
 NMN, μg/24 h 469 0–1,464 
 MN, μg/24 h 390 0–394 
 Cortisol, μg/dL 6.1 3.09–16.66 (3–5 pm) 
 Testosterone, ng/dL 428.06 241–827 
 Estradiol, pg/mL 39.65 0–39.8 
VariablesProbandNormal range
Renal function test 
 Creatinine, μmmol/L 73.36 44–133 
 BUN, mmol/L 5.19 2.86–7.90 
 mALB, mg/L 100.30 <30 
 URIC, μmmol/L 419.68 148.8–416.5 
Serum electrolytes 
 Potassium, mmol/L 3.70 3.5–5.3 
 Sodium, mmol/L 138.30 137–147 
 Chloride, mmol/L 102.90 99–110 
Serum hormone concentration 
 ACTH, pg/mL 9.63 0–46 
 TSH, μIU/mL 0.96 0.55–4.78 
 Thyroxine, μg/mL 7.40 4.29–12.47 
 Aldosterone, ng/dL 7.8 3.0–23.6 
 PRC, μIU/mL 98.2 2.8–39.9 
 ADRR, (ng/dL)/(μIU/mL) 0.079 <3.7 
 NMN, μg/24 h 469 0–1,464 
 MN, μg/24 h 390 0–394 
 Cortisol, μg/dL 6.1 3.09–16.66 (3–5 pm) 
 Testosterone, ng/dL 428.06 241–827 
 Estradiol, pg/mL 39.65 0–39.8 

BUN, blood urea nitrogen; URIC, uric acid; mALB, microalbuminuria; ACTH, adrenocorticotropic hormone; TSH, thyroid stimulating hormone; PRC, plasma renin concentration; ADRR, aldosterone/renin ratio; NMN, normetanephrine; MN, metanephrine.

No abnormal clinical features were detected among the proband’s family members, except for his mother who had freckles on both cheeks and had developed glioma at 22 years of age. She died of glioma when she was 25 years old.

Identification of a Splicing Mutation in NF1

A heterozygous splice site mutation c.479+1G>C in NF1 was found in the proband but not in any of the other family members. The Sanger sequencing results are shown in Figure 1b. This mutation resulted in a nucleotide G>C transversion in the donor splice site of intron 4. The SpliceSiteFinder-like, NNSPLICE, MaxEntScan, Human Splicing Finder, and NetGene2 in silico analysis predicted that c.479+1G>C may affect splicing products of pre-mRNA.

The Splicing Mutation Leading to the Skipping of Exon 4

We performed a minigene splicing assay after the splice site mutation was confirmed by Sanger sequencing. After 2.4% gel extraction, we sequenced various PCR products of WT and MT minigenes, and the results showed that exon 4 was missing in the MT minigene product (Fig. 4a, b). The mutation disrupted the normal splice donor site, resulting in the skipping of exon 4 (c.289_479del), a shift in the open reading frame, a premature termination code, and finally a truncated protein (p.Q97Vfs13).

Fig. 4.

Pathogenic splicing mutation c.479+1G>C in NF1. a Schematic diagram of the splicing process (c.479+1G>C). b Skipping of exon 4 in the MT demonstrated by Sanger sequencing. c Predicted three-dimensional models of WT and MT neurofibromin. CSRD, cysteine-serine-rich domain; SEC14-PH, Sec14-homologous domain and pleckstrin homology domain; GRD, GAP-related domain.

Fig. 4.

Pathogenic splicing mutation c.479+1G>C in NF1. a Schematic diagram of the splicing process (c.479+1G>C). b Skipping of exon 4 in the MT demonstrated by Sanger sequencing. c Predicted three-dimensional models of WT and MT neurofibromin. CSRD, cysteine-serine-rich domain; SEC14-PH, Sec14-homologous domain and pleckstrin homology domain; GRD, GAP-related domain.

Close modal

Domain Structure of WT and MT Neurofibromin by Molecular Modeling

The predicted three-dimensional models of WT and MT neurofibromin are shown in Figure 4c. The mutation (p.Q97Vfs13) resulted in the absence of most of the domains that were present in WT neurofibromin, namely the cysteine/serine-rich, GAP-related, Sec14-homologous, and pleckstrin homology domains. This finding provides new evidence for the loss-of-function variant of neurofibromin.

We identified a splicing mutation c.479+1G>C in the NF1 gene of a Chinese male patient with NF-1 who presented with multiple CALMs, scattered freckles, and refractory renovascular hypertension by next-generation sequencing. We performed in silico analysis to predict the potential harm of the mutation, then conducted a minigene splicing assay together with three-dimensional structural modeling to elucidate the pathogenesis of NF-1. We found that the c.479+1G>C mutation disrupted a normal splice donor site and led to the skipping of exon 4, which resulted in a truncated protein product (p.Q97Vfs13) in which most neurofibromin domains were missing.

Cutaneous presentations such as dermal neurofibromas and pigmentation are the most common features of NF-1 with early penetrance, whereas cardiovascular complications are non-neoplastic and easily overlooked but have significant manifestations that can potentially lead to catastrophic outcomes and reduced average life expectancy [21‒23]. Hypertension is frequently detected in patients with NF-1 and is associated with essential hypertension, renal or aortic vasculopathy, and pheochromocytoma/paraganglioma [24]. Renal vasculopathy is the primary cause of refractory renovascular hypertension in patients with NF-1, accounting for approximately 1% of patients [25, 26]. Oderich et al. [27] studied vascular abnormalities in 31 patients with clinical NF-1 and found that aneurysms or stenoses lesions in the aortic, renal, and mesenteric arteries were the most common of these abnormalities. Although the relative success rate of endovascular intervention was lower in patients with NF-1 than in those without NF-1, PTRA is still the first-line therapy for patients with NF-1 and refractory hypertension caused by renal artery vasculopathy whose BP was not controlled by oral antihypertensives, and repeated procedure is advised if initially unsuccessful [22, 28]. In this study, the proband initially presented with CALMs and refractory hypertension secondary from right renal artery stenosis. Superior mesenteric artery and hypertrophic cardiomyopathy were also detected; the latter may have resulted from long-term uncontrolled hypertension. However, despite the restenosis of renal artery appeared after PTRA procedure, the proband refused to accept a repeat PTRA because of the high cost. For NF-1 children or adolescents presenting with refractory hypertension, we recommend considering renal artery Doppler ultrasound or computed tomography angiography to detect renal vasculopathy. Timely medication or endovascular intervention can alleviate the disease, and regular follow-ups are necessary to monitor progression.

Neurofibromin is a cytoplasmic guanosine triphosphatase-activating protein that acts as a negative regulator of the RAS proto-oncogene. The correlation between vessel wall homeostasis and downstream signaling pathway activity can partially explain the pathogenesis of vasculopathy induced by NF-1. Neurofibromin-deficient vascular smooth muscle cells (VSMCs) or monocytes/macrophages were shown to increase cell proliferation and migration in response to growth factor or inflammatory states, and this response was regulated by over-activation of the RAS-MEK/ERK signaling axis, consistent with the histology of vascular disorders in NF-1 [29, 30]. Xu et al. [31] generated mice with homozygous inactivation of NF1 in VSMCs and identified increased activity of mitogen-activated protein kinase after vascular injury followed by enhanced VSMC proliferation and marked intimal hyperproliferation. In 2010, Lasater et al. [32] demonstrated that specific NF1 heterozygous bone marrow-derived hematopoietic cells were major drivers of the amplification of neointima formation by accelerating the accumulation of macrophages in NF1+/− mice. NF1+/− monocytes/macrophages produced increased numbers of reactive oxygen species, which induced the proliferation of VSMCs [33]. Moreover, NF1+/− monocytes/macrophages were found to mobilize the enhanced release of the chemokine receptor CCR2 into peripheral blood [34]. Chemokine ligand (CCL2; also known as monocyte chemotactic protein-1) is the ligand of CCR2, and elevated CCL2 was shown to promote the proliferation and migration of macrophages and activate RAS kinases [34]. By targeting these mechanisms using MEK inhibition [35], NADPH oxidase-2 inhibitor [33], or CCR2 antagonist [34], the neointima was alleviated with satisfactory effects in animal models. Understanding the exact mechanisms of vasculopathy in NF1 will contribute to the development of precision prevention and therapy of this disease; however, such an understanding requires further research.

In this study, we identified a splicing mutation c.479+1G>C in NF1 in the proband. Although his mother died of glioma without genetic testing, her freckle feature and malignant tumor history make it reasonable to speculate that the pathogenic mutation was inherited from his mother. Although this splicing mutation has been reported in a French cohort with NF-1 [36], we carried out an in vitro functional analysis for the first time and found that the mutation resulted in the skipping of exon 4. This frame shift mutation led to the production of a truncated neurofibromin protein that was 108 amino acids long and was absent of major domains, including the cysteine/serine-rich and GAP-related domains. Splicing mutations in NF1 are common, accounting for approximately 30% of the pathogenic variants identified in patients with NF-1 [7, 37, 38]. All splicing mutations that led to exon 4 skipping are summarized in Table 2. Alkindy et al. [38] summarized genotype-phenotype associations in 149 patients with NF1 and suggested there was an increased risk of brain gliomas and malignant peripheral nerve sheath tumors in the patients who carried splice site mutations, which confirms our recommendation that the progression of tumorigenesis should be closely monitored in the proband.

Table 2.

Summary of splicing mutation leading to exon 4 skipping in NF1

LocationMutationAcceptor site (end of intron 2)Exon 4Donor site (start of intron 4)Reference
Wild type …ttttag CAA…CAG gttagt… 
IVS3, E4 c.288+230_480-1812delins∼323 270bp deletion Deletion 1812bp deletion ins 323bp [39
E4 c.479G>C …ttttag CAA…CAC gttagt… [40
E4 c.479G>A …ttttag CAA…CAA gttagt… [41
IVS4 c.479+1G>A …ttttag CAA…CAG attagt… [41, 42
IVS4 c.479+1G>C …ttttag CAA…CAG cttagt… This study, [36
IVS4 c.479+2T>G …ttttag CAA…CAG ggtagt… [43
IVS4 c.479+4delA …ttttag CAA…CAG gttgtg… [42
IVS4 c.479+5G>C …ttttag CAA…CAG gttact… [36, 40
IVS4 c.479+5G>A …ttttag CAA…CAG gttaat… [44, 45
LocationMutationAcceptor site (end of intron 2)Exon 4Donor site (start of intron 4)Reference
Wild type …ttttag CAA…CAG gttagt… 
IVS3, E4 c.288+230_480-1812delins∼323 270bp deletion Deletion 1812bp deletion ins 323bp [39
E4 c.479G>C …ttttag CAA…CAC gttagt… [40
E4 c.479G>A …ttttag CAA…CAA gttagt… [41
IVS4 c.479+1G>A …ttttag CAA…CAG attagt… [41, 42
IVS4 c.479+1G>C …ttttag CAA…CAG cttagt… This study, [36
IVS4 c.479+2T>G …ttttag CAA…CAG ggtagt… [43
IVS4 c.479+4delA …ttttag CAA…CAG gttgtg… [42
IVS4 c.479+5G>C …ttttag CAA…CAG gttact… [36, 40
IVS4 c.479+5G>A …ttttag CAA…CAG gttaat… [44, 45

Molecular diagnosis has important implications in the diagnosis of NF-1, especially in patients who are paucisymptomatic or do not satisfy the clinical diagnostic criteria [41, 46]. Precise diagnosis is needed to improve the prognosis and for genetic counseling. If either parent carries a pathogenic NF1 mutation, prenatal diagnosis is recommended to prevent the birth of infants with the same mutation [10]. Moreover, identification of germline NF1 mutations could be valuable for genotype-phenotype correlation analysis to better guide the management of NF-1. Cost-effective methods and protocols for genetic diagnosis of NF1 are now available, including whole-genome sequencing [47], NF1-targeted panel sequencing [48], cDNA analysis combined with multiplex ligation-dependent probe amplification of genomic DNA [49], and NF1 transcriptome analysis [50]. However, because of the complexity and large size of the NF1 gene, sensitive localization of pathogenic variants is challenging and will probably require the combination of multiple detection technologies [47].

In conclusion, we identified a splicing mutation c.479+1G>C in a Chinese patient with NF-1 by next-generation sequencing. A minigene splicing assay showed that this mutation led to the skipping of exon 4, an open reading frame shift, and termination codon at 108, resulting in a truncated neurofibromin protein (p.Q97Vfs13). Taken together, our results highlight the importance of genetic testing in the diagnosis of NF-1 and confirm that the minigene splicing assay is a useful strategy for analyzing mRNA expression outcomes caused by donor splicing site variation.

We thank all participants for taking part in this study.

This study was approved by the Ethics Committee of Fuwai Hospital, approval number 2020-1321. All of the participants provided written informed consent.

The authors have no conflicts of interest to declare.

This work was supported by CAMS Innovation Fund for Medical Sciences (CIFMS, 2022-I2M-C&T-A-010 and 2022-I2M-C&T-B-041).

Conception and design: Yi-Ting Lu, Xian-Liang Zhou, and Peng Fan; administrative support: Xian-Liang Zhou and Peng Fan; provision of study materials or patients: Di Zhang, Kun-Qi Yang, and Tao Tian; collection and assembly of data: Yi-Ting Lu, Di Zhang, Buweimairemu·Rejiepu, and Dong-Cheng Cai; data analysis and interpretation: Yi-Ting Lu, Buweimairemu·Rejiepu, Dong-Cheng Cai, and Kun-Qi Yang; manuscript writing and final approval of the manuscript: all authors.

All data generated or analyzed during this study are included in this article. Further inquiries can be directed to the corresponding author.

1.
Huson
SM
,
Compston
DA
,
Clark
P
,
Harper
PS
.
A genetic study of von Recklinghausen neurofibromatosis in south east Wales. I. Prevalence, fitness, mutation rate, and effect of parental transmission on severity
.
J Med Genet
.
1989 Nov
26
11
704
11
.
2.
Evans
DG
,
Howard
E
,
Giblin
C
,
Clancy
T
,
Spencer
H
,
Huson
SM
.
Birth incidence and prevalence of tumor-prone syndromes: estimates from a UK family genetic register service
.
Am J Med Genet A
.
2010 Feb
152A
2
327
32
.
3.
Uusitalo
E
,
Leppavirta
J
,
Koffert
A
,
Suominen
S
,
Vahtera
J
,
Vahlberg
T
.
Incidence and mortality of neurofibromatosis: a total population study in Finland
.
J Invest Dermatol
.
2015 Mar
135
3
904
6
.
4.
Solares
I
,
Vinal
D
,
Morales-Conejo
M
,
Rodriguez-Salas
N
,
Feliu
J
.
Novel molecular targeted therapies for patients with neurofibromatosis type 1 with inoperable plexiform neurofibromas: a comprehensive review
.
ESMO Open
.
2021 Aug
6
4
100223
.
5.
Bergoug
M
,
Doudeau
M
,
Godin
F
,
Mosrin
C
,
Vallee
B
,
Benedetti
H
.
Neurofibromin structure, functions and regulation
.
Cells
.
2020 Oct 27
9
11
2365
.
6.
Kehrer-Sawatzki
H
,
Cooper
DN
.
Mosaicism in sporadic neurofibromatosis type 1: variations on a theme common to other hereditary cancer syndromes
.
J Med Genet
.
2008 Oct
45
10
622
31
.
7.
Koster
R
,
Brandao
RD
,
Tserpelis
D
,
van Roozendaal
CEP
,
van Oosterhoud
CN
,
Claes
KBM
.
Pathogenic neurofibromatosis type 1 (NF1) RNA splicing resolved by targeted RNAseq
.
NPJ Genom Med
.
2021 Nov 15
6
1
95
.
8.
Xu
J
,
Zhang
Y
,
Zhu
K
,
Li
J
,
Guan
Y
,
He
X
.
Clinical characteristics and in silico analysis of congenital pseudarthrosis of the tibia combined with neurofibromatosis type 1 caused by a novel NF1 mutation
.
Front Genet
.
2022
;
13
:
991314
.
9.
Ward
AJ
,
Cooper
TA
.
The pathobiology of splicing
.
J Pathol
.
2010 Jan
220
2
152
63
.
10.
Jett
K
,
Friedman
JM
.
Clinical and genetic aspects of neurofibromatosis 1
.
Genet Med
.
2010 Jan
12
1
1
11
.
11.
Li
F
,
Downing
BD
,
Smiley
LC
,
Mund
JA
,
Distasi
MR
,
Bessler
WK
.
Neurofibromin-deficient myeloid cells are critical mediators of aneurysm formation in vivo
.
Circulation
.
2014 Mar 18
129
11
1213
24
.
12.
Neurofibromatosis. Conference statement. National Institutes of Health Consensus development conference
.
Arch Neurol
.
1988 May
45
5
575
8
.
13.
Li
H
,
Handsaker
B
,
Wysoker
A
,
Fennell
T
,
Ruan
J
,
Homer
N
.
The sequence alignment/map format and SAMtools
.
Bioinformatics
.
2009 Aug 15
25
16
2078
9
.
14.
Krumm
N
,
Sudmant
PH
,
Ko
A
,
O’Roak
BJ
,
Malig
M
,
Coe
BP
.
Copy number variation detection and genotyping from exome sequence data
.
Genome Res
.
2012 Aug
22
8
1525
32
.
15.
Yeo
G
,
Burge
CB
.
Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals
.
J Comput Biol
.
2004
11
2–3
377
94
.
16.
Desmet
FO
,
Hamroun
D
,
Lalande
M
,
Collod-Beroud
G
,
Claustres
M
,
Beroud
C
.
Human Splicing Finder: an online bioinformatics tool to predict splicing signals
.
Nucleic Acids Res
.
2009 May
37
9
e67
.
17.
Heuts
DP
,
Weissenborn
MJ
,
Olkhov
RV
,
Shaw
AM
,
Gummadova
J
,
Levy
C
.
Crystal structure of a soluble form of human CD73 with ecto-5′-nucleotidase activity
.
Chembiochem
.
2012 Nov 5
13
16
2384
91
.
18.
Waterhouse
A
,
Bertoni
M
,
Bienert
S
,
Studer
G
,
Tauriello
G
,
Gumienny
R
.
SWISS-MODEL: homology modelling of protein structures and complexes
.
Nucleic Acids Res
.
2018 Jul 2
46
W1
W296
303
.
19.
Baek
M
,
DiMaio
F
,
Anishchenko
I
,
Dauparas
J
,
Ovchinnikov
S
,
Lee
GR
.
Accurate prediction of protein structures and interactions using a three-track neural network
.
Science
.
2021 Aug 20
373
6557
871
6
.
20.
Laskowski
RA
,
MacArthur
MW
,
Moss
DS
,
Thornton
JM
.
PROCHECK: a program to check the stereochemical quality of protein structures
.
J Appl Crystallogr
.
1993
;
26
(
2
):
283
91
.
21.
Rasmussen
SA
,
Yang
Q
,
Friedman
JM
.
Mortality in neurofibromatosis 1: an analysis using U.S. death certificates
.
Am J Hum Genet
.
2001 May
68
5
1110
8
.
22.
Friedman
JM
,
Arbiser
J
,
Epstein
JA
,
Gutmann
DH
,
Huot
SJ
,
Lin
AE
.
Cardiovascular disease in neurofibromatosis 1: report of the NF1 cardiovascular task force
.
Genet Med
.
2002 May-Jun
4
3
105
11
.
23.
Persu
A
,
Canning
C
,
Prejbisz
A
,
Dobrowolski
P
,
Amar
L
,
Chrysochou
C
.
Beyond atherosclerosis and fibromuscular dysplasia: rare causes of renovascular hypertension
.
Hypertension
.
2021 Sep
78
4
898
911
.
24.
Lu
YT
,
Zhang
D
,
Liu
XC
,
Zhang
QY
,
Dong
XQ
,
Fan
P
.
Identification of NF1 frameshift variants in two Chinese families with neurofibromatosis type 1 and early-onset hypertension
.
Front Pediatr
.
2021
;
9
:
785982
.
25.
Williams
VC
,
Lucas
J
,
Babcock
MA
,
Gutmann
DH
,
Korf
B
,
Maria
BL
.
Neurofibromatosis type 1 revisited
.
Pediatrics
.
2009 Jan
123
1
124
33
.
26.
de Oliveira Campos
JL
,
Bitencourt
L
,
Pedrosa
AL
,
Silva
DF
,
Lin
FJJ
,
de Oliveira Dias
LT
.
Renovascular hypertension in pediatric patients: update on diagnosis and management
.
Pediatr Nephrol
.
2021 Dec
36
12
3853
68
.
27.
Oderich
GS
,
Sullivan
TM
,
Bower
TC
,
Gloviczki
P
,
Miller
DV
,
Babovic-Vuksanovic
D
.
Vascular abnormalities in patients with neurofibromatosis syndrome type I: clinical spectrum, management, and results
.
J Vasc Surg
.
2007 Sep
46
3
475
84
.
28.
Gutmann
DH
,
Ferner
RE
,
Listernick
RH
,
Korf
BR
,
Wolters
PL
,
Johnson
KJ
.
Neurofibromatosis type 1
.
Nat Rev Dis Primers
.
2017 Feb 23
3
17004
.
29.
Li
F
,
Munchhof
AM
,
White
HA
,
Mead
LE
,
Krier
TR
,
Fenoglio
A
.
Neurofibromin is a novel regulator of RAS-induced signals in primary vascular smooth muscle cells
.
Hum Mol Genet
.
2006 Jun 1
15
11
1921
30
.
30.
Stansfield
BK
,
Bessler
WK
,
Mali
R
,
Mund
JA
,
Downing
BD
,
Kapur
R
.
Ras-Mek-Erk signaling regulates Nf1 heterozygous neointima formation
.
Am J Pathol
.
2014 Jan
184
1
79
85
.
31.
Xu
J
,
Ismat
FA
,
Wang
T
,
Yang
J
,
Epstein
JA
.
NF1 regulates a Ras-dependent vascular smooth muscle proliferative injury response
.
Circulation
.
2007 Nov 6
116
19
2148
56
.
32.
Lasater
EA
,
Li
F
,
Bessler
WK
,
Estes
ML
,
Vemula
S
,
Hingtgen
CM
.
Genetic and cellular evidence of vascular inflammation in neurofibromin-deficient mice and humans
.
J Clin Invest
.
2010 Mar
120
3
859
70
.
33.
Bessler
WK
,
Hudson
FZ
,
Zhang
H
,
Harris
V
,
Wang
Y
,
Mund
JA
.
Neurofibromin is a novel regulator of Ras-induced reactive oxygen species production in mice and humans
.
Free Radic Biol Med
.
2016 Aug
97
212
22
.
34.
Bessler
WK
,
Kim
G
,
Hudson
FZ
,
Mund
JA
,
Mali
R
,
Menon
K
.
Nf1+/- monocytes/macrophages induce neointima formation via CCR2 activation
.
Hum Mol Genet
.
2016 Mar 15
25
6
1129
39
.
35.
Tritz
R
,
Hudson
FZ
,
Harris
V
,
Ghoshal
P
,
Singla
B
,
Lin
H
.
MEK inhibition exerts temporal and myeloid cell-specific effects in the pathogenesis of neurofibromatosis type 1 arteriopathy
.
Sci Rep
.
2021 Dec 21
11
1
24345
.
36.
Sabbagh
A
,
Pasmant
E
,
Imbard
A
,
Luscan
A
,
Soares
M
,
Blanche
H
.
NF1 molecular characterization and neurofibromatosis type I genotype-phenotype correlation: the French experience
.
Hum Mutat
.
2013 Nov
34
11
1510
8
.
37.
Valero
MC
,
Martin
Y
,
Hernandez-Imaz
E
,
Marina Hernandez
A
,
Melean
G
,
Valero
AM
.
A highly sensitive genetic protocol to detect NF1 mutations
.
J Mol Diagn
.
2011 Mar
13
2
113
22
.
38.
Alkindy
A
,
Chuzhanova
N
,
Kini
U
,
Cooper
DN
,
Upadhyaya
M
.
Genotype-phenotype associations in neurofibromatosis type 1 (NF1): an increased risk of tumor complications in patients with NF1 splice-site mutations
.
Hum Genomics
.
2012 Aug 13
6
12
.
39.
Hsiao
MC
,
Piotrowski
A
,
Callens
T
,
Fu
C
,
Wimmer
K
,
Claes
KB
.
Decoding NF1 intragenic copy-number variations
.
Am J Hum Genet
.
2015 Aug 6
97
2
238
49
.
40.
Wimmer
K
,
Roca
X
,
Beiglbock
H
,
Callens
T
,
Etzler
J
,
Rao
AR
.
Extensive in silico analysis of NF1 splicing defects uncovers determinants for splicing outcome upon 5’ splice-site disruption
.
Hum Mutat
.
2007 Jun
28
6
599
612
.
41.
Pasmant
E
,
Parfait
B
,
Luscan
A
,
Goussard
P
,
Briand-Suleau
A
,
Laurendeau
I
.
Neurofibromatosis type 1 molecular diagnosis: what can NGS do for you when you have a large gene with loss of function mutations
.
Eur J Hum Genet
.
2015 May
23
5
596
601
.
42.
Pros
E
,
Gomez
C
,
Martin
T
,
Fabregas
P
,
Serra
E
,
Lazaro
C
.
Nature and mRNA effect of 282 different NF1 point mutations: focus on splicing alterations
.
Hum Mutat
.
2008 Sep
29
9
E173
93
.
43.
Brinckmann
A
,
Mischung
C
,
Bassmann
I
,
Kuhnisch
J
,
Schuelke
M
,
Tinschert
S
.
Detection of novel NF1 mutations and rapid mutation prescreening with Pyrosequencing
.
Electrophoresis
.
2007 Dec
28
23
4295
301
.
44.
Bolcekova
A
,
Nemethova
M
,
Zatkova
A
,
Hlinkova
K
,
Pozgayova
S
,
Hlavata
A
.
Clustering of mutations in the 5’ tertile of the NF1 gene in Slovakia patients with optic pathway glioma
.
Neoplasma
.
2013
;
60
(
6
):
655
65
.
45.
Nemethova
M
,
Bolcekova
A
,
Ilencikova
D
,
Durovcikova
D
,
Hlinkova
K
,
Hlavata
A
.
Thirty-nine novel neurofibromatosis 1 (NF1) gene mutations identified in Slovak patients
.
Ann Hum Genet
.
2013 Sep
77
5
364
79
.
46.
Kokkinou
E
,
Roka
K
,
Alexopoulos
A
,
Tsina
E
,
Nikas
I
,
Krallis
P
.
Development of a multidisciplinary clinic of neurofibromatosis type 1 and other neurocutaneous disorders in Greece. A 3-year experience
.
Postgrad Med
.
2019 Sep
131
7
445
52
.
47.
Pacot
L
,
Pelletier
V
,
Chansavang
A
,
Briand-Suleau
A
,
Burin des Roziers
C
,
Coustier
A
.
Contribution of whole genome sequencing in the molecular diagnosis of mosaic partial deletion of the NF1 gene in neurofibromatosis type 1
.
Hum Genet
.
2023
;
142
:
1
9
.
48.
Gieldon
L
,
Masjkur
JR
,
Richter
S
,
Darr
R
,
Lahera
M
,
Aust
D
.
Next-generation panel sequencing identifies NF1 germline mutations in three patients with pheochromocytoma but no clinical diagnosis of neurofibromatosis type 1
.
Eur J Endocrinol
.
2018 Feb
178
2
K1
9
.
49.
Alesi
V
,
Lepri
FR
,
Dentici
ML
,
Genovese
S
,
Sallicandro
E
,
Bejo
K
.
Intragenic inversions in NF1 gene as pathogenic mechanism in neurofibromatosis type 1
.
Eur J Hum Genet
.
2022 Nov
30
11
1239
43
.
50.
Douben
HCW
,
Nellist
M
,
van Unen
L
,
Elfferich
P
,
Kasteleijn
E
,
Hoogeveen-Westerveld
M
.
High-yield identification of pathogenic NF1 variants by skin fibroblast transcriptome screening after apparently normal diagnostic DNA testing
.
Hum Mutat
.
2022 Oct 17
43
12
2130
40
.