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.

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.

regioneReloaded: evaluating the association of multiple genomic region sets.

MOTIVATION: Next-generation sequencing methods continue improving the annotation of genomes in part by determining the distribution of features such as epigenetic marks. Evaluating and interpreting the association between genomic regions and their features has become a common and challenging analysis in genomic and epigenomic studies. RESULTS: With regioneR we provided an R package allowing to assess the statistical significance of pairwise associations between genomic region sets using permutation tests. We now present the R package regioneReloaded that builds upon regioneR's statistical foundation and extends the functionality for the simultaneous analysis and visualization of the associations between multiple genomic region sets. Thus, we provide a novel discovery tool for the identification of significant associations that warrant to be tested for functional interdependence. AVAILABILITY AND IMPLEMENTATION: regioneReloaded is an R package released under an Artistic-2.0 License. The source code and documentation are freely available through Bioconductor: http://www.bioconductor.org/packages/regioneReloaded.

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.

Deep genomic analysis of malignant peripheral nerve sheath tumor cell lines challenges current malignant peripheral nerve sheath tumor diagnosis.

Malignant peripheral nerve sheath tumors (MPNSTs) are soft-tissue sarcomas of the peripheral nervous system that develop either sporadically or in the context of neurofibromatosis type 1 (NF1). MPNST diagnosis can be challenging and treatment outcomes are poor. We present here a resource consisting of the genomic characterization of 9 widely used human MPNST cell lines for their use in translational research. NF1-related cell lines recapitulated primary MPNST copy number profiles, exhibited NF1, CDKN2A, and SUZ12/EED tumor suppressor gene (TSG) inactivation, and presented no gain-of-function mutations. In contrast, sporadic cell lines collectively displayed different TSG inactivation patterns and presented kinase-activating mutations, fusion genes, altered mutational frequencies and COSMIC signatures, and different methylome-based classifications. Cell lines re-classified as melanomas and other sarcomas exhibited a different drug-treatment response. Deep genomic analysis, methylome-based classification, and cell-identity marker expression, challenged the identity of common MPNST cell lines, opening an opportunity to revise MPNST differential diagnosis.

Modeling iPSC-derived human neurofibroma-like tumors in mice uncovers the heterogeneity of Schwann cells within plexiform neurofibromas.

Plexiform neurofibromas (pNFs) are developmental tumors that appear in neurofibromatosis type 1 individuals, constituting a major source of morbidity and potentially transforming into a highly metastatic sarcoma (MPNST). pNFs arise after NF1 inactivation in a cell of the neural crest (NC)-Schwann cell (SC) lineage. Here, we develop an iPSC-based NC-SC in vitro differentiation system and construct a lineage expression roadmap for the analysis of different 2D and 3D NF models. The best model consists of generating heterotypic spheroids (neurofibromaspheres) composed of iPSC-derived differentiating NF1(-/-) SCs and NF1(+/-) pNF-derived fibroblasts (Fbs). Neurofibromaspheres form by maintaining highly proliferative NF1(-/-) cells committed to the NC-SC axis due to SC-SC and SC-Fb interactions, resulting in SC linage cells at different maturation points. Upon engraftment on the mouse sciatic nerve, neurofibromaspheres consistently generate human NF-like tumors. Analysis of expression roadmap genes in human pNF single-cell RNA-seq data uncovers the presence of SC subpopulations at distinct differentiation states.

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.

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.

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.

Using antisense oligonucleotides for the physiological modulation of the alternative splicing of NF1 exon 23a during PC12 neuronal differentiation.

Neurofibromatosis Type 1 (NF1) is a genetic condition affecting approximately 1:3500 persons worldwide. The NF1 gene codes for neurofibromin protein, a GTPase activating protein (GAP) and a negative regulator of RAS. The NF1 gene undergoes alternative splicing of exon 23a (E23a) that codes for 21 amino acids placed at the center of the GAP related domain (GRD). E23a-containing type II neurofibromin exhibits a weaker Ras-GAP activity compared to E23a-less type I isoform. Exon E23a has been related with the cognitive impairment present in NF1 individuals. We designed antisense Phosphorodiamidate Morpholino Oligomers (PMOs) to modulate E23a alternative splicing at physiological conditions of gene expression and tested their impact during PC12 cell line neuronal differentiation. Results show that any dynamic modification of the natural ratio between type I and type II isoforms disturbed neuronal differentiation, altering the proper formation of neurites and deregulating both the MAPK/ERK and cAMP/PKA signaling pathways. Our results suggest an opposite regulation of these pathways by neurofibromin and the possible existence of a feedback loop sensing neurofibromin-related signaling. The present work illustrates the utility of PMOs to study alternative splicing that could be applied to other alternatively spliced genes in vitro and in vivo.

KIF11 and KIF15 mitotic kinesins are potential therapeutic vulnerabilities for malignant peripheral nerve sheath tumors.

BACKGROUND: Malignant peripheral nerve sheath tumor (MPNST) constitutes the leading cause of neurofibromatosis type 1-related mortality. MPNSTs contain highly rearranged hyperploid genomes and exhibit a high division rate and aggressiveness. We have studied in vitro whether the mitotic kinesins KIF11, KIF15, and KIF23 have a functional role in maintaining MPNST cell survival and can represent potential therapeutic vulnerabilities. METHODS: We studied the expression of kinesin mRNAs and proteins in tumors and cell lines and used several in vitro functional assays to analyze the impact of kinesin genetic suppression (KIF15, KIF23) and drug inhibition (KIF11) in MPNST cells. We also performed in vitro combined treatments targeting KIF11 together with other described MPNST targets. RESULTS: The studied kinesins were overexpressed in MPNST samples. KIF15 and KIF23 were required for the survival of MPNST cell lines, which were also more sensitive than benign control fibroblasts to the KIF11 inhibitors ispinesib and ARRY-520. Co-targeting KIF11 and BRD4 with ARRY-520 and JQ1 reduced MPNST cell viability, synergistically killing a much higher proportion of MPNST cells than control fibroblasts. In addition, genetic suppression of KIF15 conferred an increased sensitivity to KIF11 inhibitors alone or in combination with JQ1. CONCLUSIONS: The mitotic spindle kinesins KIF11 and KIF15 and the cytokinetic kinesin KIF23 play a clear role in maintaining MPNST cell survival and may represent potential therapeutic vulnerabilities. Although further in vivo evidences are still mandatory, we propose a simultaneous suppression of KIF11, KIF15, and BRD4 as a potential therapy for MPNSTs.

Use of patient derived orthotopic xenograft models for real-time therapy guidance in a pediatric sporadic malignant peripheral nerve sheath tumor.

BACKGROUND: The aim of this study was to test the feasibility and utility of developing patient-derived orthotopic xenograft (PDOX) models for patients with malignant peripheral nerve sheath tumors (MPNSTs) to aid therapeutic interventions in real time. PATIENT & METHODS: A sporadic relapsed MPNST developed in a 14-year-old boy was engrafted in mice, generating a PDOX model for use in co-clinical trials after informed consent. SNP-array and exome sequencing was performed on the relapsed tumor. Genomics, drug availability, and published literature guided PDOX treatments. RESULTS: A MPNST PDOX model was generated and expanded. Analysis of the patient's relapsed tumor revealed mutations in the MAPK1, EED, and CDK2NA/B genes. First, the PDOX model was treated with the same therapeutic regimen as received by the patient (everolimus and trametinib); after observing partial response, tumors were left to regrow. Regrown tumors were treated based on mutations (palbociclib and JQ1), drug availability, and published literature (nab-paclitaxel; bevacizumab; sorafenib plus doxorubicin; and gemcitabine plus docetaxel). The patient had a lung metastatic relapse and was treated according to PDOX results, first with nab-paclitaxel, second with sorafenib plus doxorubicin after progression, although a complete response was not achieved and multiple metastasectomies were performed. The patient is currently disease free 46 months after first relapse. CONCLUSION: Our results indicate the feasibility of generating MPNST-PDOX and genomic characterization to guide treatment in real time. Although the treatment responses observed in our model did not fully recapitulate the patient's response, this pilot study identify key aspects to improve our co-clinical testing approach in real time.

Evaluation of CNV detection tools for NGS panel data in genetic diagnostics.

Although germline copy-number variants (CNVs) are the genetic cause of multiple hereditary diseases, detecting them from targeted next-generation sequencing data (NGS) remains a challenge. Existing tools perform well for large CNVs but struggle with single and multi-exon alterations. The aim of this work is to evaluate CNV calling tools working on gene panel NGS data and their suitability as a screening step before orthogonal confirmation in genetic diagnostics strategies. Five tools (DECoN, CoNVaDING, panelcn.MOPS, ExomeDepth, and CODEX2) were tested against four genetic diagnostics datasets (two in-house and two external) for a total of 495 samples with 231 single and multi-exon validated CNVs. The evaluation was performed using the default and sensitivity-optimized parameters. Results showed that most tools were highly sensitive and specific, but the performance was dataset dependant. When evaluating them in our diagnostics scenario, DECoN and panelcn.MOPS detected all CNVs with the exception of one mosaic CNV missed by DECoN. However, DECoN outperformed panelcn.MOPS specificity achieving values greater than 0.90 when using the optimized parameters. In our in-house datasets, DECoN and panelcn.MOPS showed the highest performance for CNV screening before orthogonal confirmation. Benchmarking and optimization code is freely available at https://github.com/TranslationalBioinformaticsIGTP/CNVbenchmarkeR .

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.

Reprogramming Captures the Genetic and Tumorigenic Properties of Neurofibromatosis Type 1 Plexiform Neurofibromas.

Neurofibromatosis type 1 (NF1) is a tumor predisposition genetic disease caused by mutations in the NF1 tumor suppressor gene. Plexiform neurofibromas (PNFs) are benign Schwann cell (SC) tumors of the peripheral nerve sheath that develop through NF1 inactivation and can progress toward a malignant soft tissue sarcoma. There is a lack of non-perishable model systems to investigate PNF development. We reprogrammed PNF-derived NF1(-/-) cells, descendants from the tumor originating cell. These NF1(-/-)-induced pluripotent stem cells (iPSCs) captured the genomic status of PNFs and were able to differentiate toward neural crest stem cells and further to SCs. iPSC-derived NF1(-/-) SCs exhibited a continuous high proliferation rate, poor myelination ability, and a tendency to form 3D spheres that expressed the same markers as their PNF-derived primary SC counterparts. They represent a valuable model to study and treat PNFs. PNF-derived iPSC lines were banked for making them available.

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.

A 4-gene expression prognostic signature might guide post-remission therapy in patients with intermediate-risk cytogenetic acute myeloid leukemia.

In intermediate-risk cytogenetic acute myeloid leukemia (IRC-AML) patients, novel biomarkers to guide post-remission therapy are needed. We analyzed with high-density arrays 40 IRC-AML patients who received a non-allogeneic hematopoietic stem-cell transplantation-based post-remission therapy, and identified a signature that correlated with early relapse. Subsequently, we analyzed selected 187 genes in 49 additional IRC-AML patients by RT-PCR. BAALC, MN1, SPARC and HOPX overexpression correlated to refractoriness. BAALC or ALDH2 overexpression correlated to shorter overall survival (OS) (5-year OS: 33 ± 8.6% vs. 73.7 ± 10.1%, p = .006; 32 ± 9.3% vs. 66.4 ± 9.7%, p = .016), whereas GPR44 or TP53INP1 overexpression correlated to longer survival (5-year OS: 66.7 ± 10.3% vs. 35.4 ± 9.1%, p = .04; 58.3 ± 8.2% vs. 23.1 ± 11.7%, p = .029). A risk-score combining these four genes expression distinguished low-risk and high-risk patients (5-year OS: 79 ± 9% vs. 30 ± 8%, respectively; p = .001) in our cohort and in an independent set of patients from a public repository. Our 4-gene signature may add prognostic information and guide post-remission treatment in IRC-AML patients.

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.

Benchmarking of Whole Exome Sequencing and Ad Hoc Designed Panels for Genetic Testing of Hereditary Cancer.

Next generation sequencing panels have been developed for hereditary cancer, although there is some debate about their cost-effectiveness compared to exome sequencing. The performance of two panels is compared to exome sequencing. Twenty-four patients were selected: ten with identified mutations (control set) and fourteen suspicious of hereditary cancer but with no mutation (discovery set). TruSight Cancer (94 genes) and a custom panel (122 genes) were assessed alongside exome sequencing. Eighty-three genes were targeted by the two panels and exome sequencing. More than 99% of bases had a read depth of over 30x in the panels, whereas exome sequencing covered 94%. Variant calling with standard settings identified the 10 mutations in the control set, with the exception of MSH6 c.255dupC using TruSight Cancer. In the discovery set, 240 unique non-silent coding and canonic splice-site variants were identified in the panel genes, 7 of them putatively pathogenic (in ATM, BARD1, CHEK2, ERCC3, FANCL, FANCM, MSH2). The three approaches identified a similar number of variants in the shared genes. Exomes were more expensive than panels but provided additional data. In terms of cost and depth, panels are a suitable option for genetic diagnostics, although exomes also identify variants in non-targeted genes.