Transposon Mutagenesis-Guided CRISPR/Cas9 Screening Strongly Implicates Dysregulation of Hippo/YAP Signaling in Malignant Peripheral Nerve Sheath Tumor Development

Simple Summary

Malignant peripheral nerve sheath tumors (MPNSTs) are highly aggressive tumors with a complex genetic landscape. Patients with neurofibromatosis type 1 syndrome (NF1) are at a high risk of MPNSTs, which usually develop from pre-existing benign Schwann cell tumors called plexiform neurofibromas. In this study, we aimed to find genes that, when altered, resulted in MPNST development. Our results suggest that the functional genetic landscape of human MPNST is complex and implicates the hippo/Yes Activated Protein (YAP) pathway in the transformation of neurofibromas.

Abstract

Malignant peripheral nerve sheath tumors (MPNSTs) are highly aggressive, genomically complex, have soft tissue sarcomas, and are derived from the Schwann cell lineage. Patients with neurofibromatosis type 1 syndrome (NF1), an autosomal dominant tumor predisposition syndrome, are at a high risk for MPNSTs, which usually develop from pre-existing benign Schwann cell tumors called plexiform neurofibromas. NF1 is characterized by loss-of-function mutations in the NF1 gene, which encode neurofibromin, a Ras GTPase activating protein (GAP) and negative regulator of RasGTP-dependent signaling. In addition to bi-allelic loss of NF1, other known tumor suppressor genes include TP53, CDKN2A, SUZ12, and EED, all of which are often inactivated in the process of MPNST growth. A sleeping beauty (SB) transposon-based genetic screen for high-grade Schwann cell tumors in mice, and comparative genomics, implicated Wnt/β-catenin, PI3K-AKT-mTOR, and other pathways in MPNST development and progression. We endeavored to more systematically test genes and pathways implicated by our SB screen in mice, i.e., in a human immortalized Schwann cell-based model and a human MPNST cell line, using CRISPR/Cas9 technology. We individually induced loss-of-function mutations in 103 tumor suppressor genes (TSG) and oncogene candidates. We assessed anchorage-independent growth, transwell migration, and for a subset of genes, tumor formation in vivo. When tested in a loss-of-function fashion, about 60% of all TSG candidates resulted in the transformation of immortalized human Schwann cells, whereas 30% of oncogene candidates resulted in growth arrest in a MPNST cell line. Individual loss-of-function mutations in the TAOK1, GDI2, NF1, and APC genes resulted in transformation of immortalized human Schwann cells and tumor formation in a xenograft model. Moreover, the loss of all four of these genes resulted in activation of Hippo/Yes Activated Protein (YAP) signaling. By combining SB transposon mutagenesis and CRISPR/Cas9 screening, we established a useful pipeline for the validation of MPNST pathways and genes. Our results suggest that the functional genetic landscape of human MPNST is complex and implicate the Hippo/YAP pathway in the transformation of neurofibromas. It is thus imperative to functionally validate individual cancer genes and pathways using human cell-based models, to determinate their role in different stages of MPNST development, growth, and/or metastasis.

Highly multiplexed genome engineering using CRISPR/Cas9 gRNA arrays

The CRISPR/Cas9 system is an RNA guided nuclease system that evolved as a mechanism of adaptive immunity in bacteria. This system has been adopted for numerous genome engineering applications in research and recently, therapeutics. The CRISPR/Cas9 system has been largely implemented by delivery of Cas9 as protein, RNA, or plasmid along with a chimeric crRNA-tracrRNA guide RNA (gRNA) under the expression of a pol III promoter, such as U6. Using this approach, multiplex genome engineering has been achieved by delivering several U6-gRNA plasmids targeting multiple loci. However, this approach is limited due to the efficiently of delivering multiple plasmids to a single cell at one time. To augment the capability and accessibility of multiplexed genome engineering, we developed an efficient golden gate based method to assemble gRNAs linked by optimal Csy4 ribonuclease sequences to deliver up to 10 gRNAs as a single gRNA array transcript. Here we report the optimal expression of our guide RNA array under a strong pol II promoter. This system can be implemented alongside the myriad of CRISPR applications, allowing users to model complex biological processes requiring numerous gRNAs.

A Genome-Wide Scan Identifies Variants in NFIB Associated with Metastasis in Patients with Osteosarcoma

Metastasis is the leading cause of death in patients with osteosarcoma, the most common pediatric bone malignancy. We conducted a multistage genome-wide association study of osteosarcoma metastasis at diagnosis in 935 osteosarcoma patients to determine whether germline genetic variation contributes to risk of metastasis. We identified an SNP, rs7034162, in NFIB significantly associated with metastasis in European osteosarcoma cases, as well as in cases of African and Brazilian ancestry (meta-analysis of all cases: P = 1.2 × 10(-9); OR, 2.43; 95% confidence interval, 1.83-3.24). The risk allele was significantly associated with lowered NFIB expression, which led to increased osteosarcoma cell migration, proliferation, and colony formation. In addition, a transposon screen in mice identified a significant proportion of osteosarcomas harboring inactivating insertions in Nfib and with lowered NFIB expression. These data suggest that germline genetic variation at rs7034162 is important in osteosarcoma metastasis and that NFIB is an osteosarcoma metastasis susceptibility gene.

SIGNIFICANCE

Metastasis at diagnosis in osteosarcoma is the leading cause of death in these patients. Here we show data that are supportive for the NFIB locus as associated with metastatic potential in osteosarcoma.

A Sleeping Beauty forward genetic screen identifies new genes and pathways driving osteosarcoma development and metastasis

Osteosarcomas are sarcomas of the bone, derived from osteoblasts or their precursors, with a high propensity to metastasize. Osteosarcoma is associated with massive genomic instability, making it problematic to identify driver genes using human tumors or prototypical mouse models, many of which involve loss of Trp53 function. To identify the genes driving osteosarcoma development and metastasis, we performed a Sleeping Beauty (SB) transposon-based forward genetic screen in mice with and without somatic loss of Trp53. Common insertion site (CIS) analysis of 119 primary tumors and 134 metastatic nodules identified 232 sites associated with osteosarcoma development and 43 sites associated with metastasis, respectively. Analysis of CIS-associated genes identified numerous known and new osteosarcoma-associated genes enriched in the ErbB, PI3K-AKT-mTOR and MAPK signaling pathways. Lastly, we identified several oncogenes involved in axon guidance, including Sema4d and Sema6d, which we functionally validated as oncogenes in human osteosarcoma.