Genome Alert!: A standardized procedure for genomic variant reinterpretation and automated gene-phenotype reassessment in clinical routine.

PURPOSE: Retrospective interpretation of sequenced data in light of the current literature is a major concern of the field. Such reinterpretation is manual and both human resources and variable operating procedures are the main bottlenecks. METHODS: Genome Alert! method automatically reports changes with potential clinical significance in variant classification between releases of the ClinVar database. Using ClinVar submissions across time, this method assigns validity category to gene-disease associations. RESULTS: Between July 2017 and December 2019, the retrospective analysis of ClinVar submissions revealed a monthly median of 1247 changes in variant classification with potential clinical significance and 23 new gene-disease associations. Re-examination of 4929 targeted sequencing files highlighted 45 changes in variant classification, and of these classifications, 89% were expert validated, leading to 4 additional diagnoses. Genome Alert! gene-disease association catalog provided 75 high-confidence associations not available in the OMIM morbid list; of which, 20% became available in OMIM morbid list For more than 356 negative exome sequencing data that were reannotated for variants in these 75 genes, this elective approach led to a new diagnosis. CONCLUSION: Genome Alert! (https://genomealert.univ-grenoble-alpes.fr/) enables systematic and reproducible reinterpretation of acquired sequencing data in a clinical routine with limited human resource effect.

EZH1/2 function mostly within canonical PRC2 and exhibit proliferation-dependent redundancy that shapes mutational signatures in cancer.

Genetic mutations affecting chromatin modifiers are widespread in cancers. In malignant peripheral nerve sheath tumors (MPNSTs), Polycomb repressive complex 2 (PRC2), which plays a crucial role in gene silencing, is inactivated through recurrent mutations in core subunits embryonic ectoderm development (EED) and suppressor of zeste 12 homolog (SUZ12), but mutations in PRC2's main catalytic subunit enhancer of zeste homolog 2 (EZH2) have never been found. This is in contrast to myeloid and lymphoid malignancies, which harbor frequent loss-of-function mutations in EZH2. Here, we investigated whether the absence of EZH2 mutations in MPNST is due to a PRC2-independent (i.e., noncanonical) function of the enzyme or to redundancy with EZH1. We show that, in the absence of SUZ12, EZH2 remains bound to EED but loses its interaction with all other core and accessory PRC2 subunits. Through genetic and pharmacological analyses, we unambiguously establish that EZH2 is functionally inert in this context, thereby excluding a PRC2-independent function. Instead, we show that EZH1 and EZH2 are functionally redundant in the slowly proliferating MPNST precursors. We provide evidence that the compensatory function of EZH1 is alleviated upon higher proliferation. This work reveals how context-dependent redundancies can shape tumor-type specific mutation patterns in chromatin regulators.

PRC2 loss amplifies Ras-driven transcription and confers sensitivity to BRD4-based therapies.

The polycomb repressive complex 2 (PRC2) exerts oncogenic effects in many tumour types. However, loss-of-function mutations in PRC2 components occur in a subset of haematopoietic malignancies, suggesting that this complex plays a dichotomous and poorly understood role in cancer. Here we provide genomic, cellular, and mouse modelling data demonstrating that the polycomb group gene SUZ12 functions as tumour suppressor in PNS tumours, high-grade gliomas and melanomas by cooperating with mutations in NF1. NF1 encodes a Ras GTPase-activating protein (RasGAP) and its loss drives cancer by activating Ras. We show that SUZ12 loss potentiates the effects of NF1 mutations by amplifying Ras-driven transcription through effects on chromatin. Importantly, however, SUZ12 inactivation also triggers an epigenetic switch that sensitizes these cancers to bromodomain inhibitors. Collectively, these studies not only reveal an unexpected connection between the PRC2 complex, NF1 and Ras, but also identify a promising epigenetic-based therapeutic strategy that may be exploited for a variety of cancers.

Confirmation of mutation landscape of NF1-associated malignant peripheral nerve sheath tumors.

The commonest tumors associated with neurofibromatosis type 1 (NF1) are benign peripheral nerve sheath tumors, called neurofibromas. Malignant transformation of neurofibromas into aggressive MPNSTs may occur with a poor patient prognosis. A cooperative role of SUZ12 or EED inactivation, along with NF1, TP53, and CDKN2A loss-of-function, has been proposed to drive progression to MPNSTs. An exome sequencing analysis of eight MPNSTs, one plexiform neurofibroma, and seven cutaneous neurofibromas was undertaken. Biallelic inactivation of the NF1 gene was observed in the plexiform neurofibroma and the MPNSTs, underlining that somatic biallelic NF1 inactivation is likely to be the initiating event for plexiform neurofibroma genesis, although it is unlikely to be sufficient for the subsequent MPNST development. The majority (5/8) of MPNSTs in our analyses demonstrated homozygous or heterozygous deletions of CDKN2A, which may represent an early event following NF1 LOH in the malignant transformation of Schwann cells from plexiform neurofibroma to MPNST. Biallelic somatic alterations of SUZ12 was also found in 4/8 MPNSTs. EED biallelic alterations were detected in 2 of the other four MPNSTs, with one tumor having a homozygous EED deletion. A missense mutation in the chromatin regulator KDM2B was also identified in one MPNST. No TP53 point mutations were found in this study, confirming previous data that TP53 mutations may be relatively rare in NF1-associated MPNSTs. Our study confirms the frequent biallelic inactivation of PRC2 subunits SUZ12 and EED in MPNSTs, and suggests the implication of KDM2B.

SETD2 and DNMT3A screen in the Sotos-like syndrome French cohort.

BACKGROUND: Heterozygous NSD1 mutations were identified in 60%-90% of patients with Sotos syndrome. Recently, mutations of the SETD2 and DNMT3A genes were identified in patients exhibiting only some Sotos syndrome features. Both NSD1 and SETD2 genes encode epigenetic 'writer' proteins that catalyse methylation of histone 3 lysine 36 (H3K36me). The DNMT3A gene encodes an epigenetic 'reader' protein of the H3K36me chromatin mark. METHODS: We aimed at confirming the implication of DNMT3A and SETD2 mutations in an overgrowth phenotype, through a comprehensive targeted-next generation sequencing (NGS) screening in 210 well-phenotyped index cases with a Sotos-like phenotype and no NSD1 mutation, from a French cohort. RESULTS: Six unreported heterozygous likely pathogenic variants in DNMT3A were identified in seven patients: two nonsense variants and four de novo missense variants. One de novo unreported heterozygous frameshift variant was identified in SETD2 in one patient. All the four DNMT3A missense variants affected DNMT3A functional domains, suggesting a potential deleterious impact. DNMT3A-mutated index cases shared similar clinical features including overgrowth phenotype characterised by postnatal tall stature (≥+2SD), macrocephaly (≥+2SD), overweight or obesity at older age, intellectual deficiency and minor facial features. The phenotype associated with SETD2 mutations remains to be described more precisely. The p.Arg882Cys missense de novo constitutional DNMT3A variant found in two patients is the most frequent DNMT3A somatic mutation in acute leukaemia. CONCLUSIONS: Our results illustrate the power of targeted NGS to identify rare disease-causing variants. These observations provided evidence for a unifying mechanism (disruption of apposition and reading of the epigenetic chromatin mark H3K36me) that causes an overgrowth syndrome phenotype. Further studies are needed in order to assess the role of SETD2 and DNMT3A in intellectual deficiency without overgrowth.

NF1 single and multi-exons copy number variations in neurofibromatosis type 1.

Neurofibromatosis type 1 (NF1) is caused by dominant loss-of-function mutations of the tumor suppressor NF1 containing 57 constitutive coding exons. A huge number of different pathogenic NF1 alterations has been reported. The aim of the present study was to evaluate the usefulness of a multiplex ligation-dependent probe amplification (MLPA) approach in NF1 patients to detect single and multi-exon NF1 gene copy number variations. A genotype-phenotype correlation was then performed in NF1 patients carrying these types of genetic alterations. Among 565 NF1 index cases from the French NF1 cohort, single and multi-exon deletions/duplications screening identified NF1 partial deletions/duplications in 22 patients (~4%) using MLPA analysis. Eight single exon deletions, 11 multiple exons deletions, 1 complex rearrangement and 2 duplications were identified. All results were confirmed using a custom array-CGH. MLPA and custom array-CGH allowed the identification of rearrangements that were missed by cDNA/DNA sequencing or microsatellite analysis. We then performed a targeted next-generation sequencing of NF1 that allowed confirmation of all 22 rearrangements. No clear genotype-phenotype correlations were found for the most clinically significant disease features of NF1 in patients with single and multi-exons NF1 gene copy number changes.

Mutations in SETD2 cause a novel overgrowth condition.

BACKGROUND: Overgrowth conditions are a heterogeneous group of disorders characterised by increased growth and variable features, including macrocephaly, distinctive facial appearance and various degrees of learning difficulties and intellectual disability. Among them, Sotos and Weaver syndromes are clinically well defined and due to heterozygous mutations in NSD1 and EZH2, respectively. NSD1 and EZH2 are both histone-modifying enzymes. These two epigenetic writers catalyse two specific post-translational modifications of histones: methylation of histone 3 lysine 36 (H3K36) and lysine 27 (H3K27). We postulated that mutations in writers of these two chromatin marks could cause overgrowth conditions, resembling Sotos or Weaver syndromes, in patients with no NSD1 or EZH2 abnormalities. METHODS: We analysed the coding sequences of 14 H3K27 methylation-related genes and eight H3K36 methylation-related genes using a targeted next-generation sequencing approach in three Sotos, 11 'Sotos-like' and two Weaver syndrome patients. RESULTS: We identified two heterozygous mutations in the SETD2 gene in two patients with 'Sotos-like' syndrome: one missense p.Leu1815Trp de novo mutation in a boy and one nonsense p.Gln274* mutation in an adopted girl. SETD2 is non-redundantly responsible for H3K36 trimethylation. The two probands shared similar clinical features, including postnatal overgrowth, macrocephaly, obesity, speech delay and advanced carpal ossification. CONCLUSIONS: Our results illustrate the power of targeted next-generation sequencing to identify rare disease-causing variants. We provide a compelling argument for Sotos and Sotos-like syndromes as epigenetic diseases caused by loss-of-function mutations of epigenetic writers of the H3K36 histone mark.

EZH2 mutations in follicular lymphoma distort H3K27me3 profiles and alter transcriptional responses to PRC2 inhibition.

Mutations in chromatin regulators are widespread in cancer. Among them, the histone H3 lysine 27 methyltransferase Polycomb Repressive Complex 2 (PRC2) shows distinct alterations according to tumor type. This specificity is poorly understood. Here, we model several PRC2 alterations in one isogenic system to reveal their comparative effects. Focusing then on lymphoma-associated EZH2 mutations, we show that Ezh2Y641F induces aberrant H3K27 methylation patterns even without wild-type Ezh2, which are alleviated by partial PRC2 inhibition. Remarkably, Ezh2Y641F rewires the response to PRC2 inhibition, leading to induction of antigen presentation genes. Using a unique longitudinal follicular lymphoma cohort, we further link EZH2 status to abnormal H3K27 methylation. We also uncover unexpected variability in the mutational landscape of successive biopsies, pointing to frequent co-existence of different clones and cautioning against stratifying patients based on single sampling. Our results clarify how oncogenic PRC2 mutations disrupt chromatin and transcription, and the therapeutic vulnerabilities this creates.

NF1 molecular characterization and neurofibromatosis type I genotype-phenotype correlation: the French experience.

Neurofibromatosis type 1 (NF1) affects about one in 3,500 people in all ethnic groups. Most NF1 patients have private loss-of-function mutations scattered along the NF1 gene. Here, we present an original NF1 investigation strategy and report a comprehensive mutation analysis of 565 unrelated patients from the NF-France Network. A NF1 mutation was identified in 546 of the 565 patients, giving a mutation detection rate of 97%. The combined cDNA/DNA approach showed that a significant proportion of NF1 missense mutations (30%) were deleterious by affecting pre-mRNA splicing. Multiplex ligation-dependent probe amplification allowed the identification of restricted rearrangements that would have been missed if only sequencing or microsatellite analysis had been performed. In four unrelated families, we identified two distinct NF1 mutations within the same family. This fortuitous association points out the need to perform an exhaustive NF1 screening in the case of molecular discordant-related patients. A genotype-phenotype study was performed in patients harboring a truncating (N = 368), in-frame splicing (N = 36), or missense (N = 35) mutation. The association analysis of these mutation types with 12 common NF1 clinical features confirmed a weak contribution of the allelic heterogeneity of the NF1 mutation to the NF1 variable expressivity.

MicroRNAome profiling in benign and malignant neurofibromatosis type 1-associated nerve sheath tumors: evidences of PTEN pathway alterations in early NF1 tumorigenesis.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a common dominant tumor predisposition syndrome affecting 1 in 3,500 individuals. The hallmarks of NF1 are the development of peripheral nerve sheath tumors either benign (dermal and plexiform neurofibromas) or malignant (MPNSTs). RESULTS: To comprehensively characterize the role of microRNAs in NF1 tumorigenesis, we analyzed 377 miRNAs expression in a large panel of dermal and plexiform neurofibromas, and MPNSTs. The most significantly upregulated miRNA in plexiform neurofibromas was miR-486-3p that targets the major tumor suppressor gene, PTEN. We confirmed PTEN downregulation at mRNA level. In plexiform neurofibromas, we also report aberrant expression of four miRNAs involved in the RAS-MAPK pathway (miR-370, miR-143, miR-181a, and miR-145). In MPNSTs, significant deregulated miRNAs were involved in PTEN repression (miR-301a, miR-19a, and miR-106b), RAS-MAPK pathway regulation (Let-7b, miR-195, and miR-10b), mesenchymal transition (miR-200c, let-7b, miR-135a, miR-135b, and miR-9), HOX genes expression (miR-210, miR-196b, miR-10a, miR-10b, and miR-9), and cell cycle progression (miR-195, let-7b, miR-20a, miR-210, miR-129-3p, miR-449a, and miR-106b). CONCLUSION: We confirmed the implication of PTEN in genesis of plexiform neurofibromas and MPNSTs in NF1. Markedly deregulated miRNAs might have potential diagnostic or prognostic value and could represent novel strategies for effective pharmacological therapies of NF1 tumors.

Competitive allele specific TaqMan PCR for KRAS, BRAF and EGFR mutation detection in clinical formalin fixed paraffin embedded samples.

BACKGROUND: The development of targeted therapies has created a need for robust molecular characterization of cancer and it has become a challenge to validate methods to ensure accuracy in tumor mutation testing. METHODS: The current study was designed to evaluate KRAS, BRAF and EGFR genotyping by Competitive Allele Specific hydrolysis probes (TaqMan) PCR technology (CAST), on suboptimal formalin fixed paraffin embedded (FFPE) tumor samples. Assays were calibrated on FFPE samples and a minimal quantification cycle (Cq) cut-off was determined to standardize analyses and avoid over-interpretation of degraded material. Sensibility, specificity and blinded clinical sample screenings (n=63) were evaluated. RESULTS: CAST PCR allowed efficient amplification of FFPE samples, probes were highly specific and all assays had a sensibility inferior to 1% except for the EGFR p.T790M assay. 60/63 samples were correctly typed. The three missed mutations were EGFR exon 19 deletions that were not recognized by the DEL19 assays that were used. CONCLUSIONS: This technology is less laborious and prevent crossover of PCR products as compared to multistep methods. TaqMan® Mutation Detection assay is an important technology to consider in the field of mutation detection for KRAS, BRAF and EGFR point mutation screening. Assay calibration on FFPE samples may prevent erroneous interpretations that will ultimately harm clinical oncology practice.

Case Report: How whole-genome sequencing-based cell-free DNA prenatal testing can help identify a marker mhromosome.

A supernumerary marker chromosome (SMC) is a structurally abnormal chromosome that cannot be characterized by conventional banding cytogenetics. Marker chromosomes are present in 0.075% of prenatal cases. They are associated with variable phenotypes, ranging from normal to severely abnormal, and the prognosis is largely dependent on the results of further cytogenomic analysis. Here, we report the identification and characterization of a marker chromosome following prenatal screening in a 39-year-old pregnant patient. The patient had a normal first trimester ultrasound but was high-risk for fetal chromosome anomalies based on the results of maternal serum parameters. Chorionic villus sampling was performed, and analysis of chorionic villi revealed the presence of two identical marker chromosomes. In the interest of a rapid identification of the markers, we performed noninvasive prenatal testing (NIPT) together with chorionic villus sampling. A pericentromeric 29 Mb duplication of chromosome 20: dup (20) (p13q11.21) was identified and thereafter confirmed by targeted metaphasic FISH. Whole-genome sequencing-based NIPT was instrumental in rapid characterization of the SMCs and allowed us to obviate the need for multiple expensive and time-consuming FISH analyses.

Noninvasive Prenatal Screening for Trisomy 21 in Patients with a Vanishing Twin.

A vanishing twin (VT) occurs in up to 30% of early diagnosed twin pregnancies and is associated with an increased risk of fetal aneuploidy. Here, we describe our experience in a large VT population of 847 patients that underwent noninvasive prenatal testing (NIPT) for common fetal trisomies over a three-year period. All patients underwent an ultrasound examination prior to NIPT. Two comparison populations were included, namely, the singleton (n = 105,560) and the viable multiple gestation pregnancy samples (n = 9691) collected over the same period. All NIPT samples in the VT population received a result, of which 14 were high-risk for trisomy 21 (1.6%), nine for trisomy 18 (1.1%), and six for trisomy 13 (0.7%). Diagnostic testing confirmed the presence of trisomy 21 in 6/12 samples, giving a positive predictive value of 50%. One trisomy 18 case and no trisomy 13 cases were confirmed. The time between fetal demise and NIPT sampling did not appear to affect the number of true- or false-positive cases. In conclusion, NIPT is an effective screening method for trisomy 21 in the surviving fetus(es) in VT pregnancies. For trisomies 18 and 13, a positive NIPT should be interpreted carefully and ultrasound monitoring is preferrable over invasive diagnostic testing.

Strategy for Use of Genome-Wide Non-Invasive Prenatal Testing for Rare Autosomal Aneuploidies and Unbalanced Structural Chromosomal Anomalies.

Atypical fetal chromosomal anomalies are more frequent than previously recognized and can affect fetal development. We propose a screening strategy for a genome-wide non-invasive prenatal test (NIPT) to detect these atypical chromosomal anomalies (ACAs). Two sample cohorts were tested. Assay performances were determined using Cohort A, which consisted of 192 biobanked plasma samples-42 with ACAs, and 150 without. The rate of additional invasive diagnostic procedures was determined using Cohort B, which consisted of 3097 pregnant women referred for routine NIPT. Of the 192 samples in Cohort A, there were four initial test failures and six discordant calls; overall sensitivity was 88.1% (37/42; CI 75.00-94.81) and specificity was 99.3% (145/146; CI 96.22-99.88). In Cohort B, there were 90 first-pass failures (2.9%). The rate of positive results indicating an anomaly was 1.2% (36/3007) and 0.57% (17/3007) when limited to significant unbalanced chromosomal anomalies and trisomies 8, 9, 12, 14, 15, 16, and 22. These results show that genome-wide NIPT can screen for ACAs with an acceptable sensitivity and a small increase in invasive testing, particularly for women with increased risk following maternal serum screening and by limiting screening to structural anomalies and the most clinically meaningful trisomies.

NF1 mutations identify molecular and clinical subtypes of lung adenocarcinomas.

The tumor suppressor gene neurofibromin 1 (NF1) is a major regulator of the RAS-MAPK pathway. NF1 mutations occur in lung cancer but were not extensively explored. We hypothesized that NF1-mutated tumors could define a specific population with a distinct clinical and molecular profile. We performed NF1 sequencing using next generation sequencing (NGS) in 154 lung adenocarcinoma surgical specimens with known KRAS, EGFR, TP53, BRAF, HER2, and PIK3CA status, to evaluate the molecular and clinical specificities of NF1-mutated lung cancers. Clinical data were retrospectively collected, and their associations with molecular profiles assessed. In this series, 24 tumors were NF1 mutated (17.5%) and 11 were NF1 deleted (8%). There was no mutation hotspot. NF1 mutations were rarely associated with other RAS-MAPK pathway mutations. Most of patients with NF1 alterations were males (74.3%) and smokers (74.3%). Overall survival and disease-free survival were statistically better in patients with NF1 alterations (N = 34) than in patients with KRAS mutations (N = 30) in univariate analysis. Our results confirm that NF1 is frequently mutated and represents a distinct molecular and clinical subtype of lung adenocarcinoma.

Could Digital PCR Be an Alternative as a Non-Invasive Prenatal Test for Trisomy 21: A Proof of Concept Study.

OBJECTIVE: NIPT for fetal aneuploidy by digital PCR has been hampered by the large number of PCR reactions needed to meet statistical requirements, preventing clinical application. Here, we designed an octoplex droplet digital PCR (ddPCR) assay which allows increasing the number of available targets and thus overcomes statistical obstacles. METHOD: After technical optimization of the multiplex PCR on mixtures of trisomic and euploid DNA, we performed a validation study on samples of plasma DNA from 213 pregnant women. Molecular counting of circulating cell-free DNA was performed using a mix of hydrolysis probes targeting chromosome 21 and a reference chromosome. RESULTS: The results of our validation experiments showed that ddPCR detected trisomy 21 even when the sample's trisomic DNA content is as low as 5%. In a validation study of plasma samples from 213 pregnant women, ddPCR discriminated clearly between the trisomy 21 and the euploidy groups. CONCLUSION: Our results demonstrate that digital PCR can meet the requirements for non-invasive prenatal testing of trisomy 21. This approach is technically simple, relatively cheap, easy to implement in a diagnostic setting and compatible with ethical concerns regarding access to nucleotide sequence information. These advantages make it a potential technique of choice for population-wide screening for trisomy 21 in pregnant women.

Neurofibromatosis type 1 molecular diagnosis: what can NGS do for you when you have a large gene with loss of function mutations?

Molecular diagnosis of neurofibromatosis type 1 (NF1) is challenging owing to the large size of the tumour suppressor gene NF1, and the lack of mutation hotspots. A somatic alteration of the wild-type NF1 allele is observed in NF1-associated tumours. Genetic heterogeneity in NF1 was confirmed in patients with SPRED1 mutations. Here, we present a targeted next-generation sequencing (NGS) of NF1 and SPRED1 using a multiplex PCR approach (230 amplicons of ∼150 bp) on a PGM sequencer. The chip capacity allowed mixing 48 bar-coded samples in a 4-day workflow. We validated the NGS approach by retrospectively testing 30 NF1-mutated samples, and then prospectively analysed 279 patients in routine diagnosis. On average, 98.5% of all targeted bases were covered by at least 20X and 96% by at least 100X. An NF1 or SPRED1 alteration was found in 246/279 (88%) and 10/279 (4%) patients, respectively. Genotyping throughput was increased over 10 times, as compared with Sanger, with ∼90[euro ] for consumables per sample. Interestingly, our targeted NGS approach also provided quantitative information based on sequencing depth allowing identification of multiexons deletion or duplication. We then addressed the NF1 somatic mutation detection sensitivity in mosaic NF1 patients and tumours.

The activation of the WNT signaling pathway is a Hallmark in neurofibromatosis type 1 tumorigenesis.

PURPOSE: The hallmark of neurofibromatosis type 1 (NF1) is the onset of dermal or plexiform neurofibromas, mainly composed of Schwann cells. Plexiform neurofibromas can transform into malignant peripheral nerve sheath tumors (MPNST) that are resistant to therapies. EXPERIMENTAL DESIGN: The aim of this study was to identify an additional pathway in the NF1 tumorigenesis. We focused our work on Wnt signaling that is highly implicated in cancer, mainly in regulating the proliferation of cancer stem cells. We quantified mRNAs of 89 Wnt pathway genes in 57 NF1-associated tumors including dermal and plexiform neurofibromas and MPNSTs. Expression of two major stem cell marker genes and five major epithelial-mesenchymal transition marker genes was also assessed. The expression of significantly deregulated Wnt genes was then studied in normal human Schwann cells, fibroblasts, endothelial cells, and mast cells and in seven MPNST cell lines. RESULTS: The expression of nine Wnt genes was significantly deregulated in plexiform neurofibromas in comparison with dermal neurofibromas. Twenty Wnt genes showed altered expression in MPNST biopsies and cell lines. Immunohistochemical studies confirmed the Wnt pathway activation in NF1-associated MPNSTs. We then confirmed that the knockdown of NF1 in Schwann cells but not in epithelial cells provoked the activation of Wnt pathway by functional transfection assays. Furthermore, we showed that the protein expression of active ß-catenin was increased in NF1-silenced cell lines. Wnt pathway activation was strongly associated to both cancer stem cell reservoir and Schwann-mesenchymal transition. CONCLUSION: We highlighted the implication of Wnt pathway in NF1-associated tumorigenesis.