Single-cell genomics analysis reveals complex genetic interactions in an in vivo model of acquired BRAF inhibitor resistance.

The evolution of therapeutic resistance is a major obstacle to the success of targeted oncology drugs. While both inter- and intratumoral heterogeneity limit our ability to detect resistant subpopulations that pre-exist or emerge during treatment, our ability to analyze tumors with single-cell resolution is limited. Here, we utilized a cell-based transposon mutagenesis method to identify mechanisms of BRAF inhibitor resistance in a model of cutaneous melanoma. This screen identified overexpression of NEDD4L and VGLL3 as significant drivers of BRAF inhibitor resistance in vivo. In addition, we describe a novel single-cell genomics profiling method to genotype thousands of individual cells within tumors driven by transposon mutagenesis. This approach revealed a surprising genetic diversity among xenograft tumors and identified recurrent co-occurring mutations that emerge within distinct tumor subclones. Taken together, these observations reveal an unappreciated genetic complexity that drives BRAF inhibitor resistance.

Understanding cancer drug resistance with Sleeping Beauty functional genomic screens: Application to MAPK inhibition in cutaneous melanoma.

Combined BRAF and MEK inhibition is an effective treatment for BRAF-mutant cutaneous melanoma. However, most patients progress on this treatment due to drug resistance. Here, we applied the Sleeping Beauty transposon system to understand how melanoma evades MAPK inhibition. We found that the specific drug resistance mechanisms differed across melanomas in our genetic screens of five cutaneous melanoma cell lines. While drivers that reactivated MAPK were highly conserved, many others were cell-line specific. One such driver, VAV1, activated a de-differentiated transcriptional program like that of hyperactive RAC1, RAC1P29S. To target this mechanism, we showed that an inhibitor of SRC, saracatinib, blunts the VAV1-induced transcriptional reprogramming. Overall, we highlighted the importance of accounting for melanoma heterogeneity in treating cutaneous melanoma with MAPK inhibitors. Moreover, we demonstrated the utility of the Sleeping Beauty transposon system in understanding cancer drug resistance.

SRC-RAC1 signaling drives drug resistance to BRAF inhibition in de-differentiated cutaneous melanomas.

Rare gain-of-function mutations in RAC1 drive drug resistance to targeted BRAF inhibition in cutaneous melanoma. Here, we show that wildtype RAC1 is a critical driver of growth and drug resistance, but only in a subset of melanomas with elevated markers of de-differentiation. Similarly, SRC inhibition also selectively sensitized de-differentiated melanomas to BRAF inhibition. One possible mechanism may be the suppression of the de-differentiated state, as SRC and RAC1 maintained markers of de-differentiation in human melanoma cells. The functional differences between melanoma subtypes suggest that the clinical management of cutaneous melanoma can be enhanced by the knowledge of differentiation status. To simplify the task of classification, we developed a binary classification strategy based on a small set of ten genes. Using this gene set, we reliably determined the differentiation status previously defined by hundreds of genes. Overall, our study informs strategies that enhance the precision of BRAFi by discovering unique vulnerabilities of the de-differentiated cutaneous melanoma subtype and creating a practical method to resolve differentiation status.

WDR26 and MTF2 are therapeutic targets in multiple myeloma.

Unbiased genetic forward screening using retroviral insertional mutagenesis in a genetically engineered mouse model of human multiple myeloma may further our understanding of the genetic pathways that govern neoplastic plasma cell development. To evaluate this hypothesis, we performed a tumor induction study in MYC-transgenic mice infected as neonates with the Moloney-derived murine leukemia virus, MOL4070LTR. Next-generation DNA sequencing of proviral genomic integration sites yielded rank-ordered candidate tumor progression genes that accelerated plasma cell neoplasia in mice. Rigorous clinical and biological validation of these genes led to the discovery of two novel myeloma genes: WDR26 (WD repeat-containing protein 26) and MTF2 (metal response element binding transcription factor 2). WDR26, a core component of the carboxy-terminal to LisH (CTLH) complex, is overexpressed or mutated in solid cancers. MTF2, an ancillary subunit of the polycomb repressive complex 2 (PRC2), is a close functional relative of PHD finger protein 19 (PHF19) which is currently emerging as an important driver of myeloma. These findings underline the utility of genetic forward screens in mice for uncovering novel blood cancer genes and suggest that WDR26-CTLH and MTF2-PRC2 are promising molecular targets for new approaches to myeloma treatment and prevention.

Apobec1 complementation factor overexpression promotes hepatic steatosis, fibrosis, and hepatocellular cancer.

The RNA-binding protein Apobec1 complementation factor (A1CF) regulates posttranscriptional ApoB mRNA editing, but the range of RNA targets and the long-term effect of altered A1CF expression on liver function are unknown. Here we studied hepatocyte-specific A1cf-transgenic (A1cf+/Tg), A1cf+/Tg Apobec1-/-, and A1cf-/- mice fed chow or high-fat/high-fructose diets using RNA-Seq, RNA CLIP-Seq, and tissue microarrays from human hepatocellular cancer (HCC). A1cf+/Tg mice exhibited increased hepatic proliferation and steatosis, with increased lipogenic gene expression (Mogat1, Mogat2, Cidea, Cd36) associated with shifts in polysomal RNA distribution. Aged A1cf+/Tg mice developed spontaneous fibrosis, dysplasia, and HCC, and this development was accelerated on a high-fat/high-fructose diet and was independent of Apobec1. RNA-Seq revealed increased expression of mRNAs involved in oxidative stress (Gstm3, Gpx3, Cbr3), inflammatory response (Il19, Cxcl14, Tnfα, Ly6c), extracellular matrix organization (Mmp2, Col1a1, Col4a1), and proliferation (Kif20a, Mcm2, Mcm4, Mcm6), and a subset of mRNAs (including Sox4, Sox9, Cdh1) were identified in RNA CLIP-Seq. Increased A1CF expression in human HCC correlated with advanced fibrosis and with reduced survival in a subset with nonalcoholic fatty liver disease. In conclusion, we show that hepatic A1CF overexpression selectively alters polysomal distribution and mRNA expression, promoting lipogenic, proliferative, and inflammatory pathways leading to HCC.

Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma.

The use of selective BRAF inhibitors (BRAFi) has produced remarkable outcomes for patients with advanced cutaneous melanoma harboring a BRAFV600E mutation. Unfortunately, the majority of patients eventually develop drug-resistant disease. We employed a genetic screening approach to identify gain-of-function mechanisms of BRAFi resistance in two independent melanoma cell lines. Our screens identified both known and unappreciated drivers of BRAFi resistance, including multiple members of the DBL family. Mechanistic studies identified a DBL/RAC1/PAK signaling axis capable of driving resistance to both current and next-generation BRAFis. However, we show that the SRC inhibitor, saracatinib, can block the DBL-driven resistance. Our work highlights the utility of our straightforward genetic screening method in identifying new drug combinations to combat acquired BRAFi resistance. SIGNIFICANCE: A simple, rapid, and flexible genetic screening approach identifies genes that drive resistance to MAPK inhibitors when overexpressed in human melanoma cells.

A simplified transposon mutagenesis method to perform phenotypic forward genetic screens in cultured cells.

BACKGROUND: The introduction of genome-wide shRNA and CRISPR libraries has facilitated cell-based screens to identify loss-of-function mutations associated with a phenotype of interest. Approaches to perform analogous gain-of-function screens are less common, although some reports have utilized arrayed viral expression libraries or the CRISPR activation system. However, a variety of technical and logistical challenges make these approaches difficult for many labs to execute. In addition, genome-wide shRNA or CRISPR libraries typically contain of hundreds of thousands of individual engineered elements, and the associated complexity creates issues with replication and reproducibility for these methods. RESULTS: Here we describe a simple, reproducible approach using the SB transposon system to perform phenotypic cell-based genetic screens. This approach employs only three plasmids to perform unbiased, whole-genome transposon mutagenesis. We also describe a ligation-mediated PCR method that can be used in conjunction with the included software tools to map raw sequence data, identify candidate genes associated with phenotypes of interest, and predict the impact of recurrent transposon insertions on candidate gene function. Finally, we demonstrate the high reproducibility of our approach by having three individuals perform independent replicates of a mutagenesis screen to identify drivers of vemurafenib resistance in cultured melanoma cells. CONCLUSIONS: Collectively, our work establishes a facile, adaptable method that can be performed by labs of any size to perform robust, genome-wide screens to identify genes that influence phenotypes of interest.

Apobec1 complementation factor (A1CF) and RBM47 interact in tissue-specific regulation of C to U RNA editing in mouse intestine and liver.

Mammalian C to U RNA is mediated by APOBEC1, the catalytic deaminase, together with RNA binding cofactors (including A1CF and RBM47) whose relative physiological requirements are unresolved. Although A1CF complements APOBEC1 for in vitro RNA editing, A1cf-/- mice exhibited no change in apolipoproteinB (apoB) RNA editing, while Rbm47 mutant mice exhibited impaired intestinal RNA editing of apoB as well as other targets. Here we examined the role of A1CF and RBM47 in adult mouse liver and intestine, following deletion of either one or both gene products and also following forced (liver or intestinal) transgenic A1CF expression. There were minimal changes in hepatic and intestinal apoB RNA editing in A1cf-/- mice and no changes in either liver- or intestine-specific A1CF transgenic mice. Rbm47 liver-specific knockout (Rbm47LKO ) mice demonstrated reduced editing in a subset (11 of 20) of RNA targets, including apoB. By contrast, apoB RNA editing was virtually eliminated (<6% activity) in intestine-specific (Rbm47IKO ) mice with only five of 53 targets exhibiting C-to-U RNA editing. Double knockout of A1cf and Rbm47 in liver (ARLKO ) eliminated apoB RNA editing and reduced editing in the majority of other targets, with no changes following adenoviral APOBEC1 administration. Intestinal double knockout mice (ARIKO ) demonstrated further reduced editing (<10% activity) in four of five of the residual APOBEC1 targets identified in ARIKO mice. These data suggest that A1CF and RBM47 each function independently, yet interact in a tissue-specific manner, to regulate the activity and site selection of APOBEC1 dependent C-to-U RNA editing.

A Forward Genetic Screen Targeting the Endothelium Reveals a Regulatory Role for the Lipid Kinase Pi4ka in Myelo- and Erythropoiesis.

Given its role as the source of definitive hematopoietic cells, we sought to determine whether mutations initiated in the hemogenic endothelium would yield hematopoietic abnormalities or malignancies. Here, we find that endothelium-specific transposon mutagenesis in mice promotes hematopoietic pathologies that are both myeloid and lymphoid in nature. Frequently mutated genes included previously recognized cancer drivers and additional candidates, such as Pi4ka, a lipid kinase whose mutation was found to promote myeloid and erythroid dysfunction. Subsequent validation experiments showed that targeted inactivation of the Pi4ka catalytic domain or reduction in mRNA expression inhibited myeloid and erythroid cell differentiation in vitro and promoted anemia in vivo through a mechanism involving deregulation of AKT, MAPK, SRC, and JAK-STAT signaling. Finally, we provide evidence linking PI4KAP2, previously considered a pseudogene, to human myeloid and erythroid leukemia.

Chronic liver injury alters driver mutation profiles in hepatocellular carcinoma in mice.

Most hepatocellular carcinomas (HCCs) develop in a chronically injured liver, yet the extent to which this microenvironment promotes neoplastic transformation or influences selective pressures for genetic drivers of HCC remains unclear. We sought to determine the impact of hepatic injury in an established mouse model of HCC induced by Sleeping Beauty transposon mutagenesis. Chemically induced chronic liver injury dramatically increased tumor penetrance and significantly altered driver mutation profiles, likely reflecting distinct selective pressures. In addition to established human HCC genes and pathways, we identified several injury-associated candidates that represent promising loci for further study. Among them, we found that FIGN is overexpressed in human HCC and promotes hepatocyte invasion. We also validated Gli2's oncogenic potential in vivo, providing direct evidence that Hedgehog signaling can drive liver tumorigenesis in the context of chronic injury. Finally, we show that a subset of injury-associated candidate genes identifies two distinct classes of human HCCs. Further analysis of these two subclasses revealed significant trends among common molecular classification schemes of HCC. The genes and mechanisms identified here provide functional insights into the origin of HCC in a chronic liver damage environment. CONCLUSION: A chronically damaged liver microenvironment influences the genetic mechanisms that drive hepatocarcinogenesis. (Hepatology 2018;67:924-939).

Sleeping Beauty Insertional Mutagenesis in Mice Identifies Drivers of Steatosis-Associated Hepatic Tumors.

Hepatic steatosis is a strong risk factor for the development of hepatocellular carcinoma (HCC), yet little is known about the molecular pathology associated with this factor. In this study, we performed a forward genetic screen using Sleeping Beauty (SB) transposon insertional mutagenesis in mice treated to induce hepatic steatosis and compared the results to human HCC data. In humans, we determined that steatosis increased the proportion of female HCC patients, a pattern also reflected in mice. Our genetic screen identified 203 candidate steatosis-associated HCC genes, many of which are altered in human HCC and are members of established HCC-driving signaling pathways. The protein kinase A/cyclic AMP signaling pathway was altered frequently in mouse and human steatosis-associated HCC. We found that activated PKA expression drove steatosis-specific liver tumorigenesis in a mouse model. Another candidate HCC driver, the N-acetyltransferase NAT10, which we found to be overexpressed in human steatosis-associated HCC and associated with decreased survival in human HCC, also drove liver tumorigenesis in a steatotic mouse model. This study identifies genes and pathways promoting HCC that may represent novel targets for prevention and treatment in the context of hepatic steatosis, an area of rapidly growing clinical significance. Cancer Res; 77(23); 6576-88. ©2017 AACR.

From Peas to Disease: Modifier Genes, Network Resilience, and the Genetics of Health.

Phenotypes are rarely consistent across genetic backgrounds and environments, but instead vary in many ways depending on allelic variants, unlinked genes, epigenetic factors, and environmental exposures. In the extreme, individuals carrying the same causal DNA sequence variant but on different backgrounds can be classified as having distinct conditions. Similarly, some individuals that carry disease alleles are nevertheless healthy despite affected family members in the same environment. These genetic background effects often result from the action of so-called "modifier genes" that modulate the phenotypic manifestation of target genes in an epistatic manner. While complicating the prospects for gene discovery and the feasibility of mechanistic studies, such effects are opportunities to gain a deeper understanding of gene interaction networks that provide organismal form and function as well as resilience to perturbation. Here, we review the principles of modifier genetics and assess progress in studies of modifier genes and their targets in both simple and complex traits. We propose that modifier effects emerge from gene interaction networks whose structure and function vary with genetic background and argue that these effects can be exploited as safe and effective ways to prevent, stabilize, and reverse disease and dysfunction.

Transposon mutagenesis identifies chromatin modifiers cooperating with Ras in thyroid tumorigenesis and detects ATXN7 as a cancer gene.

Oncogenic RAS mutations are present in 15-30% of thyroid carcinomas. Endogenous expression of mutant Ras is insufficient to initiate thyroid tumorigenesis in murine models, indicating that additional genetic alterations are required. We used Sleeping Beauty (SB) transposon mutagenesis to identify events that cooperate with HrasG12V in thyroid tumor development. Random genomic integration of SB transposons primarily generated loss-of-function events that significantly increased thyroid tumor penetrance in Tpo-Cre/homozygous FR-HrasG12V mice. The thyroid tumors closely phenocopied the histological features of human RAS-driven, poorly differentiated thyroid cancers. Characterization of transposon insertion sites in the SB-induced tumors identified 45 recurrently mutated candidate cancer genes. These mutation profiles were remarkably concordant with mutated cancer genes identified in a large series of human poorly differentiated and anaplastic thyroid cancers screened by next-generation sequencing using the MSK-IMPACT panel of cancer genes, which we modified to include all SB candidates. The disrupted genes primarily clustered in chromatin remodeling functional nodes and in the PI3K pathway. ATXN7, a component of a multiprotein complex with histone acetylase activity, scored as a significant SB hit. It was recurrently mutated in advanced human cancers and significantly co-occurred with RAS or NF1 mutations. Expression of ATXN7 mutants cooperated with oncogenic RAS to induce thyroid cell proliferation, pointing to ATXN7 as a previously unrecognized cancer gene.

RNA sequencing of Sleeping Beauty transposon-induced tumors detects transposon-RNA fusions in forward genetic cancer screens.

Forward genetic screens using Sleeping Beauty (SB)-mobilized T2/Onc transposons have been used to identify common insertion sites (CISs) associated with tumor formation. Recurrent sites of transposon insertion are commonly identified using ligation-mediated PCR (LM-PCR). Here, we use RNA sequencing (RNA-seq) data to directly identify transcriptional events mediated by T2/Onc. Surprisingly, the majority (∼80%) of LM-PCR identified junction fragments do not lead to observable changes in RNA transcripts. However, in CIS regions, direct transcriptional effects of transposon insertions are observed. We developed an automated method to systematically identify T2/Onc-genome RNA fusion sequences in RNA-seq data. RNA fusion-based CISs were identified corresponding to both DNA-based CISs (Cdkn2a, Mycl1, Nf2, Pten, Sema6d, and Rere) and additional regions strongly associated with cancer that were not observed by LM-PCR (Myc, Akt1, Pth, Csf1r, Fgfr2, Wisp1, Map3k5, and Map4k3). In addition to calculating recurrent CISs, we also present complementary methods to identify potential driver events via determination of strongly supported fusions and fusions with large transcript level changes in the absence of multitumor recurrence. These methods independently identify CIS regions and also point to cancer-associated genes like Braf. We anticipate RNA-seq analyses of tumors from forward genetic screens will become an efficient tool to identify causal events.

A Hybrid Adenoviral Vector System Achieves Efficient Long-Term Gene Expression in the Liver via piggyBac Transposition.

Much research has gone into the development of hybrid gene delivery systems that combine the broad tropism and efficient transduction of adenoviral vectors with the ability to achieve stable expression of cargo genes. In addition to gene therapy applications, such a system has considerable advantages for studies of gene function in vivo, permitting fine-tuned genetic manipulation with higher throughput than can be achieved using standard transgenic and DNA targeting techniques. Existing strategies are limited, however, by low integration efficiencies, small cargo capacity, and/or a dependence on target cell division. The utility of this approach could be enhanced by a system that provides all of the following: (1) efficient delivery, (2) stable expression in a high percentage of target cells (whether mitotic or not), (3) large cargo capacity, (4) flexibility to use with a wide range of additional experimental conditions, and (5) simple experimental technique. Here we report the initial characterization of a hybrid system that meets these criteria by utilizing piggyBac (PB) transposition to achieve genomic integration from adenoviral vectors. We demonstrate stable expression of an adenovirus (Ad)-PB-delivered reporter gene in ∼20-40% of hepatocytes following standard tail vein injection. Its high efficiency and flexibility relative to existing hybrid adenoviral gene delivery approaches indicate a considerable potential utility of the Ad-PB system for therapeutic gene delivery and in vivo studies of gene function.

Sequencing methods and datasets to improve functional interpretation of sleeping beauty mutagenesis screens.

BACKGROUND: Animal models of cancer are useful to generate complementary datasets for comparison to human tumor data. Insertional mutagenesis screens, such as those utilizing the Sleeping Beauty (SB) transposon system, provide a model that recapitulates the spontaneous development and progression of human disease. This approach has been widely used to model a variety of cancers in mice. Comprehensive mutation profiles are generated for individual tumors through amplification of transposon insertion sites followed by high-throughput sequencing. Subsequent statistical analyses identify common insertion sites (CISs), which are predicted to be functionally involved in tumorigenesis. Current methods utilized for SB insertion site analysis have some significant limitations. For one, they do not account for transposon footprints - a class of mutation generated following transposon remobilization. Existing methods also discard quantitative sequence data due to uncertainty regarding the extent to which it accurately reflects mutation abundance within a heterogeneous tumor. Additionally, computational analyses generally assume that all potential insertion sites have an equal probability of being detected under non-selective conditions, an assumption without sufficient relevant data. The goal of our study was to address these potential confounding factors in order to enhance functional interpretation of insertion site data from tumors. RESULTS: We describe here a novel method to detect footprints generated by transposon remobilization, which revealed minimal evidence of positive selection in tumors. We also present extensive characterization data demonstrating an ability to reproducibly assign semi-quantitative information to individual insertion sites within a tumor sample. Finally, we identify apparent biases for detection of inserted transposons in several genomic regions that may lead to the identification of false positive CISs. CONCLUSION: The information we provide can be used to refine analyses of data from insertional mutagenesis screens, improving functional interpretation of results and facilitating the identification of genes important in cancer development and progression.

Modeling progressive non-alcoholic fatty liver disease in the laboratory mouse.

Non-alcoholic fatty liver disease (NAFLD) is the most common liver disease in the world and its prevalence is rising. In the absence of disease progression, fatty liver poses minimal risk of detrimental health outcomes. However, advancement to non-alcoholic steatohepatitis (NASH) confers a markedly increased likelihood of developing severe liver pathologies, including fibrosis, cirrhosis, organ failure, and cancer. Although a substantial percentage of NAFLD patients develop NASH, the genetic and molecular mechanisms driving this progression are poorly understood, making it difficult to predict which patients will ultimately develop advanced liver disease. Deficiencies in mechanistic understanding preclude the identification of beneficial prognostic indicators and the development of effective therapies. Mouse models of progressive NAFLD serve as a complementary approach to the direct analysis of human patients. By providing an easily manipulated experimental system that can be rigorously controlled, they facilitate an improved understanding of disease development and progression. In this review, we discuss genetically- and chemically-induced models of NAFLD that progress to NASH, fibrosis, and liver cancer in the context of the major signaling pathways whose disruption has been implicated as a driving force for their development. Additionally, an overview of nutritional models of progressive NAFLD is provided.

The stress-regulated transcription factor CHOP promotes hepatic inflammatory gene expression, fibrosis, and oncogenesis.

Viral hepatitis, obesity, and alcoholism all represent major risk factors for hepatocellular carcinoma (HCC). Although these conditions also lead to integrated stress response (ISR) or unfolded protein response (UPR) activation, the extent to which these stress pathways influence the pathogenesis of HCC has not been tested. Here we provide multiple lines of evidence demonstrating that the ISR-regulated transcription factor CHOP promotes liver cancer. We show that CHOP expression is up-regulated in liver tumors in human HCC and two mouse models thereof. Chop-null mice are resistant to chemical hepatocarcinogenesis, and these mice exhibit attenuation of both apoptosis and cellular proliferation. Chop-null mice are also resistant to fibrosis, which is a key risk factor for HCC. Global gene expression profiling suggests that deletion of CHOP reduces the levels of basal inflammatory signaling in the liver. Our results are consistent with a model whereby CHOP contributes to hepatic carcinogenesis by promoting inflammation, fibrosis, cell death, and compensatory proliferation. They implicate CHOP as a common contributing factor in the development of HCC in a variety of chronic liver diseases.

Domesticated transposable element gene products in human cancer.

The adaptation of transposable elements inserted within the genome to serve novel functions in a host cell, a process known as molecular domestication, is a widespread phenomenon in nature. Around fifty protein-coding genes in humans have arisen through this mechanism. Functional characterization of these domesticated genes has revealed involvement in a multitude of diverse cellular processes. Some of these functions are related to cellular activities and pathways known to be involved in cancer development. In this mini-review we discuss such roles of domesticated genes that may be aberrantly regulated in human cancer, as well as studies that have identified disrupted expression in tumors. We also describe studies that have provided definitive experimental evidence for transposable element-derived gene products in promoting tumorigenesis.