Screening of CNVs using NGS data improves mutation detection yield and decreases costs in genetic testing for hereditary cancer.

INTRODUCTION: Germline CNVs are important contributors to hereditary cancer. In genetic diagnostics, multiplex ligation-dependent probe amplification (MLPA) is commonly used to identify them. However, MLPA is time-consuming and expensive if applied to many genes, hence many routine laboratories test only a subset of genes of interest. METHODS AND RESULTS: We evaluated a next-generation sequencing (NGS)-based CNV detection tool (DECoN) as first-tier screening to decrease costs and turnaround time and expand CNV analysis to all genes of clinical interest in our diagnostics routine. We used DECoN in a retrospective cohort of 1860 patients where a limited number of genes were previously analysed by MLPA, and in a prospective cohort of 2041 patients, without MLPA analysis. In the retrospective cohort, 6 new CNVs were identified and confirmed by MLPA. In the prospective cohort, 19 CNVs were identified and confirmed by MLPA, 8 of these would have been lost in our previous MLPA-restricted detection strategy. Also, the number of genes tested by MLPA across all samples decreased by 93.0% in the prospective cohort. CONCLUSION: Including an in silico germline NGS CNV detection tool improved our genetic diagnostics strategy in hereditary cancer, both increasing the number of CNVs detected and reducing turnaround time and costs.

Neurofibromatosis type 1 families with first-degree relatives harbouring distinct NF1 pathogenic variants. Genetic counselling and familial diagnosis: what should be offered?

Neurofibromatosis type 1 (NF1) is an autosomal dominant disorder caused by pathogenic variants in NF1 Recently, NF1 testing has been included as a clinical criterion for NF1 diagnosis. Additionally, preconception genetic counselling in patients with NF1 focuses on a 50% risk of transmitting the familial variant as the risk of having a sporadic NF1 is considered the same as the general population. METHODS: 829 individuals, 583 NF1 sporadic cases and 246 patients with NF1 with documented family history, underwent genetic testing for NF1. Genotyping and segregation analysis of NF1 familial variants was determined by microsatellite analysis and NF1 sequencing. RESULTS: The mutational analysis of NF1 in 154 families with two or more affected cases studied showed the co-occurrence of two different NF1 germline pathogenic variants in four families. The estimated mutation rate in those families was 3.89×10-3, 20 times higher than the NF1 mutation rate (~2×10-4) (p=0.0008). Furthermore, the co-occurrence of two different NF1 germline pathogenic variants in these families was 1:39, 60 times the frequency of sporadic NF1 (1:2500) (p=0.003). In all cases, the de novo NF1 pathogenic variant was present in a descendant of an affected male. In two cases, variants were detected in the inherited paternal wild-type allele. CONCLUSIONS: Our results, together with previous cases reported, suggest that the offspring of male patients with NF1 could have an increased risk of experiencing de novo NF1 pathogenic variants. This observation, if confirmed in additional cohorts, could have relevant implications for NF1 genetic counselling, family planning and NF1 genetic testing.

Revisiting the UK Genetic Severity Score for NF2: a proposal for the addition of a functional genetic component.

BACKGROUND: Neurofibromatosis type 2 (NF2) is an autosomal dominant disorder characterised by the development of multiple schwannomas, especially on vestibular nerves, and meningiomas. The UK NF2 Genetic Severity Score (GSS) is useful to predict the progression of the disease from germline NF2 pathogenic variants, which allows the clinical follow-up and the genetic counselling offered to affected families to be optimised. METHODS: 52 Spanish patients were classified using the GSS, and patients' clinical severity was measured and compared between GSS groups. The GSS was reviewed with the addition of phenotype quantification, genetic variant classification and functional assays of Merlin and its downstream pathways. Principal component analysis and regression models were used to evaluate the differences between severity and the effect of NF2 germline variants. RESULTS: The GSS was validated in the Spanish NF2 cohort. However, for 25% of mosaic patients and patients harbouring variants associated with mild and moderate phenotypes, it did not perform as well for predicting clinical outcomes as it did for pathogenic variants associated with severe phenotypes. We studied the possibility of modifying the mutation classification in the GSS by adding the impact of pathogenic variants on the function of Merlin in 27 cases. This revision helped to reduce variability within NF2 mutation classes and moderately enhanced the correlation between patient phenotype and the different prognosis parameters analysed (R2=0.38 vs R2=0.32, p>0001). CONCLUSIONS: We validated the UK NF2 GSS in a Spanish NF2 cohort, despite the significant phenotypic variability identified within it. The revision of the GSS, named Functional Genetic Severity Score, could add value for the classification of mosaic patients and patients showing mild and moderate phenotypes once it has been validated in other cohorts.

Targeted massively parallel sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas: An analysis of 51 individuals in a single-center experience.

Constitutional LZTR1 or SMARCB1 pathogenic variants (PVs) have been found in ∼86% of familial and ∼40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1, and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit PVs in the LZTR1/SMARCB1/NF2 coding sequences; however, with 22q deletion and/or a different NF2 PV in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) PVs could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3'-untranslated region (3'-UTR) c.*17C>T, c.*70C>T, or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3'-UTR variants on the transcript. The 3'-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting the SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat-masked genomics sequences of chromosome 22q in individuals negative for PVs in the SMARCB1/LZTR1/NF2 coding- and noncoding regions revealed five potential schwannomatosis-predisposing candidate genes, that is, MYO18B, NEFH, SGSM1, SGSM3, and SBF1, pending further verification.

CNVfilteR: an R/Bioconductor package to identify false positives produced by germline NGS CNV detection tools.

SUMMARY: Germline copy-number variants (CNVs) are relevant mutations for multiple genetics fields, such as the study of hereditary diseases. However, available benchmarks show that all next-generation sequencing (NGS) CNV calling tools produce false positives. We developed CNVfilteR, an R package that uses the single-nucleotide variant calls usually obtained in germline NGS pipelines to identify those false positives. The package can detect both false deletions and false duplications. We evaluated CNVfilteR performance on callsets generated by 13 CNV calling tools on three whole-genome sequencing and 541 panel samples, showing a decrease of up to 44.8% in false positives and consistent F1-score increase. Using CNVfilteR to detect false-positive calls can improve the overall performance of existing CNV calling pipelines. AVAILABILITY AND IMPLEMENTATION: CNVfilteR is released under Artistic-2.0 License. Source code and documentation are freely available at Bioconductor (http://www.bioconductor.org/packages/CNVfilteR). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromosomal translocations inactivating CDKN2A support a single path for malignant peripheral nerve sheath tumor initiation.

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas with poor prognosis, developing either sporadically or in persons with neurofibromatosis type 1 (NF1). Loss of CDKN2A/B is an important early event in MPNST progression. However, many reported MPNSTs exhibit partial or no inactivation of CDKN2A/B, raising the question of whether there is more than one molecular path for MPNST initiation. We present here a comprehensive genomic analysis of MPNST cell lines and tumors to explore in depth the status of CDKN2A. After accounting for CDKN2A deletions and point mutations, we uncovered a previously unnoticed high frequency of chromosomal translocations involving CDKN2A in both MPNST cell lines and primary tumors. Most identified translocation breakpoints were validated by PCR amplification and Sanger sequencing. Many breakpoints clustered in an intronic 500 bp hotspot region adjacent to CDKN2A exon 2. We demonstrate the bi-allelic inactivation of CDKN2A in all tumors (n = 15) and cell lines (n = 8) analyzed, supporting a single molecular path for MPNST initiation in both sporadic and NF1-related MPNSTs. This general CDKN2A inactivation in MPNSTs has implications for MPNST diagnostics and treatment. Our findings might be relevant for other tumor types with high frequencies of CDKN2A inactivation.

Mutational spectrum by phenotype: panel-based NGS testing of patients with clinical suspicion of RASopathy and children with multiple café-au-lait macules.

Children with neurofibromatosis type 1 (NF1) may exhibit an incomplete clinical presentation, making difficult to reach a clinical diagnosis. A phenotypic overlap may exist in children with other RASopathies or with other genetic conditions if only multiple café-au-lait macules (CALMs) are present. The syndromes that can converge in these inconclusive phenotypes have different clinical courses. In this context, an early genetic testing has been proposed to be clinically useful to manage these patients. We present the validation and implementation into diagnostics of a custom NGS panel (I2HCP, ICO-IMPPC Hereditary Cancer Panel) for testing patients with a clinical suspicion of a RASopathy (n = 48) and children presenting multiple CALMs (n = 102). We describe the mutational spectrum and the detection rates identified in these two groups of individuals. We identified pathogenic variants in 21 out of 48 patients with clinical suspicion of RASopathy, with mutations in NF1 accounting for 10% of cases. Furthermore, we identified pathogenic mutations mainly in the NF1 gene, but also in SPRED1, in more than 50% of children with multiple CALMs, exhibiting an NF1 mutational spectrum different from a group of clinically diagnosed NF1 patients (n = 80). An NGS panel strategy for the genetic testing of these two phenotype-defined groups outperforms previous strategies.

Analysis of intratumor heterogeneity in Neurofibromatosis type 1 plexiform neurofibromas and neurofibromas with atypical features: Correlating histological and genomic findings.

Plexiform neurofibromas (PNFs) are benign peripheral nerve sheath tumors involving large nerves present in 30%-50% Neurofibromatosis type 1 (NF1) patients. Atypical neurofibromas (ANF) are distinct nodular lesions with atypical features on histology that arise from PNFs. The risk and timeline of malignant transformation in ANF is difficult to assess. A recent NIH workshop has stratified ANFs and separated a subgroup with multiple atypical features and higher risk of malignant transformation termed atypical neurofibromatous neoplasms with uncertain biological potential (ANNUBP). We performed an analysis of intratumor heterogeneity on eight PNFs to link histological and genomic findings. Tumors were homogeneous although histological and molecular heterogeneity was identified. All tumors were 2n, almost mutation-free and had a clonal NF1(-/-) origin. Two ANFs from the same patient showed atypical features on histology and deletions of CDKN2A/B. One of the ANFs exhibited different areas in which the degree of histological atypia correlated with the heterozygous or homozygous loss of the CDKN2A/B loci. CDKN2A/B deletions in different areas originated independently. Results may indicate that loss of a single CDKN2A/B copy in NF1(-/-) cells is sufficient to start ANF development and that total inactivation of both copies of CDKN2A/B is necessary to form an ANNUBP.

karyoploteR: an R/Bioconductor package to plot customizable genomes displaying arbitrary data.

MOTIVATION: Data visualization is a crucial tool for data exploration, analysis and interpretation. For the visualization of genomic data there lacks a tool to create customizable non-circular plots of whole genomes from any species. RESULTS: We have developed karyoploteR, an R/Bioconductor package to create linear chromosomal representations of any genome with genomic annotations and experimental data plotted along them. Plot creation process is inspired in R base graphics, with a main function creating karyoplots with no data and multiple additional functions, including custom functions written by the end-user, adding data and other graphical elements. This approach allows the creation of highly customizable plots from arbitrary data with complete freedom on data positioning and representation. AVAILABILITY AND IMPLEMENTATION: karyoploteR is released under Artistic-2.0 License. Source code and documentation are freely available through Bioconductor (http://www.bioconductor.org/packages/karyoploteR) and at the examples and tutorial page at https://bernatgel.github.io/karyoploter_tutorial. CONTACT: bgel@igtp.cat.

2016 Children's Tumor Foundation conference on neurofibromatosis type 1, neurofibromatosis type 2, and schwannomatosis.

Organized and hosted by the Children's Tumor Foundation (CTF), the Neurofibromatosis (NF) conference is the premier annual gathering for clinicians and researchers interested in neurofibromatosis type 1 (NF1), neurofibromatosis type 2 (NF2), and schwannomatosis (SWN). The 2016 edition constituted a blend of clinical and basic aspects of NF research that helped in clarifying different advances in the field. The incorporation of next generation sequencing is changing the way genetic diagnostics is performed for NF and related disorders, providing solutions to problems like genetic heterogeneity, overlapping clinical manifestations, or the presence of mosaicism. The transformation from plexiform neurofibroma (PNF) to malignant peripheral nerve sheath tumor (MPNST) is being clarified, along with new management and treatments for benign and premalignant tumors. Promising new cellular and in vivo models for understanding the musculoskeletal abnormalities in NF1, the development of NF2 or SWN associated schwannomas, and clarifying the cells that give rise to NF1-associated optic pathway glioma were presented. The interaction of neurofibromin and SPRED1 was described comprehensively, providing functional insight that will help in the interpretation of pathogenicity of certain missense variants identified in NF1 and Legius syndrome patients. Novel promising imaging techniques are being developed, as well as new integrative and holistic management models for patients that take into account psychological, social, and biological factors. Importantly, new therapeutic approaches for schwannomas, meningiomas, ependymomas, PNF, and MPNST are being pursued. This report highlights the major advances that were presented at the 2016 CTF NF conference.

regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests.

MOTIVATION: Statistically assessing the relation between a set of genomic regions and other genomic features is a common challenging task in genomic and epigenomic analyses. Randomization based approaches implicitly take into account the complexity of the genome without the need of assuming an underlying statistical model. SUMMARY: regioneR is an R package that implements a permutation test framework specifically designed to work with genomic regions. In addition to the predefined randomization and evaluation strategies, regioneR is fully customizable allowing the use of custom strategies to adapt it to specific questions. Finally, it also implements a novel function to evaluate the local specificity of the detected association. AVAILABILITY AND IMPLEMENTATION: regioneR is an R package released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor (http://www.bioconductor.org/packages/regioneR). CONTACT: rmalinverni@carrerasresearch.org.

Epigenetic induction of the Ink4a/Arf locus prevents Schwann cell overproliferation during nerve regeneration and after tumorigenic challenge.

The number of Schwann cells is fitted to axonal length in peripheral nerves. This relationship is lost when tumorigenic stimuli induce uncontrolled Schwann cell proliferation, generating tumours such us neurofibromas and schwannomas. Schwann cells also re-enter the cell cycle following nerve injury during the process of Wallerian degeneration. In both cases proliferation is finally arrested. We show that in neurofibroma, the induction of Jmjd3 (jumonji domain containing 3, histone lysine demethylase) removes trimethyl groups on lysine-27 of histone-H3 and epigenetically activates the Ink4a/Arf-locus, forcing Schwann cells towards replicative senescence. Remarkably, blocking this mechanism allows unrestricted proliferation, inducing malignant transformation of neurofibromas. Interestingly, our data suggest that in injured nerves, Schwann cells epigenetically activate the same locus to switch off proliferation and enter the senescence programme. Indeed, when this pathway is genetically blocked, Schwann cells fail to drop out of the cell cycle and continue to proliferate. We postulate that the Ink4a/Arf-locus is expressed as part of a physiological response that prevents uncontrolled proliferation of the de-differentiated Schwann cell generated during nerve regeneration, a response that is also activated to avoid overproliferation after tumorigenic stimuli in the peripheral nervous system.

Negative feedback that improves information transmission in yeast signalling.

Haploid Saccharomyces cerevisiae yeast cells use a prototypic cell signalling system to transmit information about the extracellular concentration of mating pheromone secreted by potential mating partners. The ability of cells to respond distinguishably to different pheromone concentrations depends on how much information about pheromone concentration the system can transmit. Here we show that the mitogen-activated protein kinase Fus3 mediates fast-acting negative feedback that adjusts the dose response of the downstream system response to match the dose response of receptor-ligand binding. This 'dose-response alignment', defined by a linear relationship between receptor occupancy and downstream response, can improve the fidelity of information transmission by making downstream responses corresponding to different receptor occupancies more distinguishable and reducing amplification of stochastic noise during signal transmission. We also show that one target of the feedback is a previously uncharacterized signal-promoting function of the regulator of G-protein signalling protein Sst2. Our work suggests that negative feedback is a general mechanism used in signalling systems to align dose responses and thereby increase the fidelity of information transmission.

Unbalancing cAMP and Ras/MAPK pathways as a therapeutic strategy for cutaneous neurofibromas.

Cutaneous neurofibromas (cNFs) are benign Schwann cell (SC) tumors arising from subepidermal glia. Individuals with neurofibromatosis type 1 (NF1) may develop thousands of cNFs, which greatly affect their quality of life. cNF growth is driven by the proliferation of NF1-/- SCs and their interaction with the NF1+/- microenvironment. We analyzed the crosstalk between human cNF-derived SCs and fibroblasts (FBs), identifying an expression signature specific to the SC-FB interaction. We validated the secretion of proteins involved in immune cell migration, suggesting a role of SC-FB crosstalk in immune cell recruitment. The signature also captured components of developmental signaling pathways, including the cAMP elevator G protein-coupled receptor 68 (GPR68). Activation of Gpr68 by ogerin in combination with the MEK inhibitor (MEKi) selumetinib reduced viability and induced differentiation and death of human cNF-derived primary SCs, a result corroborated using an induced pluripotent stem cell-derived 3D neurofibromasphere model. Similar results were obtained using other Gpr68 activators or cAMP analogs/adenylyl cyclase activators in combination with selumetinib. Interestingly, whereas primary SC cultures restarted their proliferation after treatment with selumetinib alone was stopped, the combination of ogerin-selumetinib elicited a permanent halt on SC expansion that persisted after drug removal. These results indicate that unbalancing the Ras and cAMP pathways by combining MEKi and cAMP elevators could be used as a potential treatment for cNFs.

When Unsuspected Crystallinity Ruins Biological Testing in Early Discovery: A Case Study.

The impact of the crystalline or amorphous structure of a solid on the solubility and pharmacokinetic properties of a drug candidate is always considered by the pharmaceutical industry during the development of a new drug; however, it is not so frequently considered during the early drug discovery process by organic and medicinal chemists, particularly those working in academia. We want to share, as an example, the false negative obtained in the biological testing of a solid sample of a tyrosine kinase inhibitor due to its unexpected crystallinity and lower solubility with respect to a solid amorphous batch of the same compound and the experimentation carried out to establish the origin of such a discrepancy.

Expanding a precision medicine platform for malignant peripheral nerve sheath tumors: New patient-derived orthotopic xenografts, cell lines and tumor entities.

Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive soft-tissue sarcomas with a poor survival rate, presenting either sporadically or in the context of neurofibromatosis type 1 (NF1). The histological diagnosis of MPNSTs can be challenging, with different tumors exhibiting great histological and marker expression overlap. This heterogeneity could be partly responsible for the observed disparity in treatment response due to the inherent diversity of the preclinical models used. For several years, our group has been generating a large patient-derived orthotopic xenograft (PDOX) MPNST platform for identifying new precision medicine treatments. Herein, we describe the expansion of this platform using six primary tumors clinically diagnosed as MPNSTs, from which we obtained six additional PDOX mouse models and three cell lines, thus generating three pairs of in vitro-in vivo models. We extensively characterized these tumors and derived preclinical models, including genomic, epigenomic, and histological analyses. Tumors were reclassified after these analyses: three remained as MPNSTs (two being classic MPNSTs), one was a melanoma, another was a neurotrophic tyrosine receptor kinase (NTRK)-rearranged spindle cell neoplasm, and, finally, the last was an unclassifiable tumor bearing neurofibromin-2 (NF2) inactivation, a neuroblastoma RAS viral oncogene homolog (NRAS) oncogenic mutation, and a SWI/SNF-related matrix-associated actin-dependent regulator of chromatin (SMARCA4) heterozygous truncated variant. New cell lines and PDOXs faithfully recapitulated histology, marker expression, and genomic characteristics of the primary tumors. The diversity in tumor identity and their specific associated genomic alterations impacted treatment responses obtained when we used the new cell lines for testing compounds against known altered pathways in MPNSTs. In summary, we present here an extension of our MPNST precision medicine platform, with new PDOXs and cell lines, including tumor entities confounded as MPNSTs in a real clinical scenario. This platform may constitute a useful tool for obtaining correct preclinical information to guide MPNST clinical trials.

Probe-based quantitative PCR assay for detecting constitutional and somatic deletions in the NF1 gene: application to genetic testing and tumor analysis.

BACKGROUND: About 5% of patients with neurofibromatosis type 1 (NF1) bear constitutional microdeletions that encompass NF1 (neurofibromin 1) and neighboring genes. These patients are characterized by the development of a high number of dermal neurofibromas (dNFs), mental retardation, and an increased risk of developing a malignant peripheral nerve sheath tumor (MPNST). Additionally, 10% of somatic second hits identified in dNFs are caused by deletions involving the NF1 gene. To detect constitutional and somatic deletions, we developed a probe-based quantitative PCR (qPCR) assay for interrogating the copy number status of 11 loci distributed along a 2.8-Mb region around the NF1 gene. METHODS: We developed the qPCR assay with Universal ProbeLibrary technology (Roche) and designed a Microsoft Excel spreadsheet to analyze qPCR data for copy number calculations. The assay fulfilled the essential aspects of the MIQE (minimum information for publication of quantitative real-time PCR experiments) guidelines and used the qBase relative quantification framework for calculations. RESULTS: The assay was validated with a set of DNA samples with known constitutional or somatic NF1 deletions. The assay showed high diagnostic sensitivity and specificity and distinguished between Type-1, Type-2, and atypical constitutional microdeletions in 14 different samples. It also identified 16 different somatic deletions in dNFs. These results were confirmed by multiplex ligation-dependent probe amplification. CONCLUSIONS: The qPCR assay provides a methodology for detecting constitutional NF1 microdeletions that could be incorporated as an additional technique in a genetic-testing setting. It also permits the identification of somatic NF1 deletions in tissues with a high percentage of cells bearing 2 copies of the NF1 gene.

Generation of human iPSC-derived neurofibromaspheres for in vitro and in vivo uses.

Neurofibromas are benign peripheral nervous system tumors associated with neurofibromatosis type 1, which originate from NF1(-/-) Schwann cell precursors. We describe a protocol to generate neurofibromaspheres by differentiating NF1(-/-) Schwann cells from induced pluripotent stem cells and combining them with neurofibroma primary fibroblasts. We also describe the development of neurofibroma-like tumors when neurofibromaspheres are engrafted in the sciatic nerve of nude mice. This model constitutes a versatile platform for drug screening and the study of neurofibroma biology. For complete details on the use and execution of this protocol, please refer to Mazuelas et al. (2022).1.

Existing and Developing Preclinical Models for Neurofibromatosis Type 1-Related Cutaneous Neurofibromas.

Neurofibromatosis type 1 (NF1) is caused by a nonfunctional copy of the NF1 tumor suppressor gene that predisposes patients to the development of cutaneous neurofibromas (cNFs), the skin tumor that is the hallmark of this condition. Innumerable benign cNFs, each appearing by an independent somatic inactivation of the remaining functional NF1 allele, form in nearly all patients with NF1. One of the limitations in developing a treatment for cNFs is an incomplete understanding of the underlying pathophysiology and limitations in experimental modeling. Recent advances in preclinical in vitro and in vivo modeling have substantially enhanced our understanding of cNF biology and created unprecedented opportunities for therapeutic discovery. We discuss the current state of cNF preclinical in vitro and in vivo model systems, including two- and three-dimensional cell cultures, organoids, genetically engineered mice, patient-derived xenografts, and porcine models. We highlight the models' relationship to human cNFs and how they can be used to gain insight into cNF development and therapeutic discovery.

Cutaneous Neurofibroma Heterogeneity: Factors that Influence Tumor Burden in Neurofibromatosis Type 1.

Neurofibromatosis type 1 is one of the most common genetic disorders of the nervous system and predisposes patients to develop benign and malignant tumors. Cutaneous neurofibromas (cNFs) are NF1-associated benign tumors that affect nearly 100% of patients with NF1. cNFs dramatically reduce patients' QOL owing to their unaesthetic appearance, physical discomfort, and corresponding psychological burden. There is currently no effective drug therapy option, and treatment is restricted to surgical removal. One of the greatest hurdles for cNF management is the variability of clinical expressivity in NF1, resulting in intrapatient and interpatient cNF tumor burden heterogeneity, that is, the variability in the presentation and evolution of these tumors. There is growing evidence that a wide array of factors are involved in the regulation of cNF heterogeneity. Understanding the mechanisms underlying this heterogeneity of cNF at the molecular, cellular, and environmental levels can facilitate the development of innovative and personalized treatment regimens.