Functional genomic landscape of cancer-intrinsic evasion of killing by T cells.

The genetic circuits that allow cancer cells to evade destruction by the host immune system remain poorly understood1-3. Here, to identify a phenotypically robust core set of genes and pathways that enable cancer cells to evade killing mediated by cytotoxic T lymphocytes (CTLs), we performed genome-wide CRISPR screens across a panel of genetically diverse mouse cancer cell lines that were cultured in the presence of CTLs. We identify a core set of 182 genes across these mouse cancer models, the individual perturbation of which increases either the sensitivity or the resistance of cancer cells to CTL-mediated toxicity. Systematic exploration of our dataset using genetic co-similarity reveals the hierarchical and coordinated manner in which genes and pathways act in cancer cells to orchestrate their evasion of CTLs, and shows that discrete functional modules that control the interferon response and tumour necrosis factor (TNF)-induced cytotoxicity are dominant sub-phenotypes. Our data establish a central role for genes that were previously identified as negative regulators of the type-II interferon response (for example, Ptpn2, Socs1 and Adar1) in mediating CTL evasion, and show that the lipid-droplet-related gene Fitm2 is required for maintaining cell fitness after exposure to interferon-γ (IFNγ). In addition, we identify the autophagy pathway as a conserved mediator of the evasion of CTLs by cancer cells, and show that this pathway is required to resist cytotoxicity induced by the cytokines IFNγ and TNF. Through the mapping of cytokine- and CTL-based genetic interactions, together with in vivo CRISPR screens, we show how the pleiotropic effects of autophagy control cancer-cell-intrinsic evasion of killing by CTLs and we highlight the importance of these effects within the tumour microenvironment. Collectively, these data expand our knowledge of the genetic circuits that are involved in the evasion of the immune system by cancer cells, and highlight genetic interactions that contribute to phenotypes associated with escape from killing by CTLs.

Chemical genomics with pyrvinium identifies C1orf115 as a regulator of drug efflux.

Pyrvinium is a quinoline-derived cyanine dye and an approved anti-helminthic drug reported to inhibit WNT signaling and have anti-proliferative effects in various cancer cell lines. To further understand the mechanism by which pyrvinium is cytotoxic, we conducted a pooled genome-wide CRISPR loss-of-function screen in the human HAP1 cell model. The top drug-gene sensitizer interactions implicated the malate-aspartate and glycerol-3-phosphate shuttles as mediators of cytotoxicity to mitochondrial complex I inhibition including pyrvinium. By contrast, perturbation of the poorly characterized gene C1orf115/RDD1 resulted in strong resistance to the cytotoxic effects of pyrvinium through dysregulation of the major drug efflux pump ABCB1/MDR1. Interestingly, C1orf115/RDD1 was found to physically associate with ABCB1/MDR1 through proximity-labeling experiments and perturbation of C1orf115 led to mis-localization of ABCB1/MDR1. Our results are consistent with a model whereby C1orf115 modulates drug efflux through regulation of the major drug exporter ABCB1/MDR1.

Recurrent Germline DLST Mutations in Individuals with Multiple Pheochromocytomas and Paragangliomas.

Pheochromocytomas and paragangliomas (PPGLs) provide some of the clearest genetic evidence for the critical role of metabolism in the tumorigenesis process. Approximately 40% of PPGLs are caused by driver germline mutations in 16 known susceptibility genes, and approximately half of these genes encode members of the tricarboxylic acid (TCA) cycle. Taking as a starting point the involvement of the TCA cycle in PPGL development, we aimed to identify unreported mutations that occurred in genes involved in this key metabolic pathway and that could explain the phenotypes of additional individuals who lack mutations in known susceptibility genes. To accomplish this, we applied a targeted sequencing of 37 TCA-cycle-related genes to DNA from 104 PPGL-affected individuals with no mutations in the major known predisposing genes. We also performed omics-based analyses, TCA-related metabolite determination, and 13C5-glutamate labeling assays. We identified five germline variants affecting DLST in eight unrelated individuals (∼7%); all except one were diagnosed with multiple PPGLs. A recurrent variant, c.1121G>A (p.Gly374Glu), found in four of the eight individuals triggered accumulation of 2-hydroxyglutarate, both in tumors and in a heterologous cell-based assay designed to functionally evaluate DLST variants. p.Gly374Glu-DLST tumors exhibited loss of heterozygosity, and their methylation and expression profiles are similar to those of EPAS1-mutated PPGLs; this similarity suggests a link between DLST disruption and pseudohypoxia. Moreover, we found positive DLST immunostaining exclusively in tumors carrying TCA-cycle or EPAS1 mutations. In summary, this study reveals DLST as a PPGL-susceptibility gene and further strengthens the relevance of the TCA cycle in PPGL development.

Preclinical toxicology profile of squalene epoxidase inhibitors.

Small cell lung cancer (SCLC) is a particularly aggressive subset of lung cancer, and identification of new therapeutic options is of significant interest. We recently reported that SCLC cell lines display a specific vulnerability to inhibition of squalene epoxidase (SQLE), an enzyme in the cholesterol biosynthetic pathway that catalyzes the conversion of squalene to 2,3-oxidosqualene. Since it has been reported that SQLE inhibition can result in dermatitis in dogs, we conducted a series of experiments to determine if SQLE inhibitors would be tolerated at exposures predicted to drive maximal efficacy in SCLC tumors. Detailed profiling of the SQLE inhibitor NB-598 showed that dogs did not tolerate predicted efficacious exposures, with dose-limiting toxicity due to gastrointestinal clinical observations, although skin toxicities were also observed. To extend these studies, two SQLE inhibitors, NB-598 and Cmpd-4″, and their structurally inactive analogs, NB-598.ia and Cmpd-4″.ia, were profiled in monkeys. While both active SQLE inhibitors resulted in dose-limiting gastrointestinal toxicity, the structurally similar inactive analogs did not. Collectively, our data demonstrate that significant toxicities arise at exposures well below the predicted levels needed for anti-tumor activity. The on-target nature of the toxicities identified is likely to limit the potential therapeutic utility of SQLE inhibition for the treatment of SCLC.

A mouse model of a human congenital disorder of glycosylation caused by loss of PMM2.

The most common congenital disorder of glycosylation (CDG), phosphomannomutase 2 (PMM2)-CDG, is caused by mutations in PMM2 that limit availability of mannose precursors required for protein N-glycosylation. The disorder has no therapy and there are no models to test new treatments. We generated compound heterozygous mice with the R137H and F115L mutations in Pmm2 that correspond to the most prevalent alleles found in patients with PMM2-CDG. Many Pmm2R137H/F115L mice died prenatally, while survivors had significantly stunted growth. These animals and cells derived from them showed protein glycosylation deficiencies similar to those found in patients with PMM2-CDG. Growth-related glycoproteins insulin-like growth factor (IGF) 1, IGF binding protein-3 and acid-labile subunit, along with antithrombin III, were all deficient in Pmm2R137H/F115L mice, but their levels in heterozygous mice were comparable to wild-type (WT) littermates. These imbalances, resulting from defective glycosylation, are likely the cause of the stunted growth seen both in our model and in PMM2-CDG patients. Both Pmm2R137H/F115L mouse and PMM2-CDG patient-derived fibroblasts displayed reductions in PMM activity, guanosine diphosphate mannose, lipid-linked oligosaccharide precursor and total cellular protein glycosylation, along with hypoglycosylation of a new endogenous biomarker, glycoprotein 130 (gp130). Over-expression of WT-PMM2 in patient-derived fibroblasts rescued all these defects, showing that restoration of mutant PMM2 activity is a viable therapeutic strategy. This functional mouse model of PMM2-CDG, in vitro assays and identification of the novel gp130 biomarker all shed light on the human disease, and moreover, provide the essential tools to test potential therapeutics for this untreatable disease.

A genome-wide RNAi screen identifies multiple RSK-dependent regulators of cell migration.

To define the functional pathways regulating epithelial cell migration, we performed a genome-wide RNAi screen using 55,000 pooled lentiviral shRNAs targeting ∼11,000 genes, selecting for transduced cells with increased motility. A stringent validation protocol generated a set of 31 genes representing diverse pathways whose knockdown dramatically enhances cellular migration. Some of these pathways share features of epithelial-to-mesenchymal transition (EMT), and together they implicate key regulators of transcription, cellular signaling, and metabolism, as well as novel modulators of cellular trafficking, such as DLG5. In delineating downstream pathways mediating these migration phenotypes, we observed universal activation of ERKs and a profound dependence on their RSK effectors. Pharmacological inhibition of RSK dramatically suppresses epithelial cell migration induced by knockdown of all 31 genes, suggesting that convergence of diverse migratory pathways on this kinase may provide a therapeutic opportunity in disorders of cell migration, including cancer metastasis.

A developmentally regulated inducer of EMT, LBX1, contributes to breast cancer progression.

Epithelial-to-mesenchymal transition (EMT) plays an important role during normal embryogenesis, and it has been implicated in cancer invasion and metastasis. Here, we report that Ladybird homeobox 1 (LBX1), a developmentally regulated homeobox gene, directs expression of the known EMT inducers ZEB1, ZEB2, Snail1, and transforming growth factor beta2 (TGFB2). In mammary epithelial cells, overexpression of LBX1 leads to morphological transformation, expression of mesenchymal markers, enhanced cell migration, increased CD44(high)/CD24(low) progenitor cell population, and tumorigenic cooperation with known oncogenes. In human breast cancer, LBX1 is up-regulated in the unfavorable estrogen receptor (ER)/progesterone (PR)/HER2 triple-negative basal-like subtype. Thus, aberrant expression of LBX1 may lead to the activation of a developmentally regulated EMT pathway in human breast cancer.

Rational design of a ligand-controlled protein conformational switch.

Design of a regulatable multistate protein is a challenge for protein engineering. Here we design a protein with a unique topology, called uniRapR, whose conformation is controlled by the binding of a small molecule. We confirm switching and control ability of uniRapR in silico, in vitro, and in vivo. As a proof of concept, uniRapR is used as an artificial regulatory domain to control activity of kinases. By activating Src kinase using uniRapR in single cells and whole organism, we observe two unique phenotypes consistent with its role in metastasis. Activation of Src kinase leads to rapid induction of protrusion with polarized spreading in HeLa cells, and morphological changes with loss of cell-cell contacts in the epidermal tissue of zebrafish. The rational creation of uniRapR exemplifies the strength of computational protein design, and offers a powerful means for targeted activation of many pathways to study signaling in living organisms.

SUMO modification is required for in vivo Hox gene regulation by the Caenorhabditis elegans Polycomb group protein SOP-2.

Post-translational modification of proteins by the ubiquitin-like molecule SUMO (sumoylation) regulates their subcellular localization and affects their functional properties in vitro, but the physiological function of sumoylation in multicellular organisms is largely unknown. Here, we show that the C. elegans Polycomb group (PcG) protein SOP-2 interacts with the SUMO-conjugating enzyme UBC-9 through its evolutionarily conserved SAM domain. Sumoylation of SOP-2 is required for its localization to nuclear bodies in vivo and for its physiological repression of Hox genes. Global disruption of sumoylation phenocopies a sop-2 mutation by causing ectopic Hox gene expression and homeotic transformations. Chimeric constructs in which the SOP-2 SAM domain is replaced with that derived from fruit fly or mammalian PcG proteins, but not those in which the SOP-2 SAM domain is replaced with the SAM domains of non-PcG proteins, confer appropriate in vivo nuclear localization and Hox gene repression. These observations indicate that sumoylation of PcG proteins, modulated by their evolutionarily conserved SAM domain, is essential to their physiological repression of Hox genes.

Tumor-intrinsic LKB1-LIF signaling axis establishes a myeloid niche to promote immune evasion and tumor growth.

Tumor mutations can influence the surrounding microenvironment leading to suppression of anti-tumor immune responses and thereby contributing to tumor progression and failure of cancer therapies. Here we use genetically engineered lung cancer mouse models and patient samples to dissect how LKB1 mutations accelerate tumor growth by reshaping the immune microenvironment. Comprehensive immune profiling of LKB1 -mutant vs wildtype tumors revealed dramatic changes in myeloid cells, specifically enrichment of Arg1 + interstitial macrophages and SiglecF Hi neutrophils. We discovered a novel mechanism whereby autocrine LIF signaling in Lkb1 -mutant tumors drives tumorigenesis by reprogramming myeloid cells in the immune microenvironment. Inhibiting LIF signaling in Lkb1 -mutant tumors, via gene targeting or with a neutralizing antibody, resulted in a striking reduction in Arg1 + interstitial macrophages and SiglecF Hi neutrophils, expansion of antigen specific T cells, and inhibition of tumor progression. Thus, targeting LIF signaling provides a new therapeutic approach to reverse the immunosuppressive microenvironment of LKB1 -mutant tumors.

Inflammation in the tumor-adjacent lung as a predictor of clinical outcome in lung adenocarcinoma.

Approximately 30% of early-stage lung adenocarcinoma patients present with disease progression after successful surgical resection. Despite efforts of mapping the genetic landscape, there has been limited success in discovering predictive biomarkers of disease outcomes. Here we performed a systematic multi-omic assessment of 143 tumors and matched tumor-adjacent, histologically-normal lung tissue with long-term patient follow-up. Through histologic, mutational, and transcriptomic profiling of tumor and adjacent-normal tissue, we identified an inflammatory gene signature in tumor-adjacent tissue as the strongest clinical predictor of disease progression. Single-cell transcriptomic analysis demonstrated the progression-associated inflammatory signature was expressed in both immune and non-immune cells, and cell type-specific profiling in monocytes further improved outcome predictions. Additional analyses of tumor-adjacent transcriptomic data from The Cancer Genome Atlas validated the association of the inflammatory signature with worse outcomes across cancers. Collectively, our study suggests that molecular profiling of tumor-adjacent tissue can identify patients at high risk for disease progression.

Molecular characterization of colorectal cancer related peritoneal metastatic disease.

A significant proportion of colorectal cancer (CRC) patients develop peritoneal metastases (PM) in the course of their disease. PMs are associated with a poor quality of life, significant morbidity and dismal disease outcome. To improve care for this patient group, a better understanding of the molecular characteristics of CRC-PM is required. Here we present a comprehensive molecular characterization of a cohort of 52 patients. This reveals that CRC-PM represent a distinct CRC molecular subtype, CMS4, but can be further divided in three separate categories, each presenting with unique features. We uncover that the CMS4-associated structural protein Moesin plays a key role in peritoneal dissemination. Finally, we define specific evolutionary features of CRC-PM which indicate that polyclonal metastatic seeding underlies these lesions. Together our results suggest that CRC-PM should be perceived as a distinct disease entity.

Coexpression of Argonaute-2 enhances RNA interference toward perfect match binding sites.

RNAi is widely applied to inhibit expression of specific genes, but it is limited by variable efficiency and specificity of empirically designed siRNA or shRNA constructs. This complicates studies targeting individual genes and significantly impairs large-scale screens using genome-wide knockdown libraries. Here, we show that ectopic expression of the RISC slicer Argonaute-2 (Ago2, eIF2C2) dramatically enhances RNAi specifically for mRNA targets with perfectly matched binding sites. This effect depends on its endonuclease activity and is uncoupled from its regulation of microRNA expression. To model the application of Ago2 coexpression with shRNA knockdown, we targeted the EGF receptor (EGFR) in lung cancer cells exhibiting oncogene addiction to EGFR. Whereas multiple empirically designed shRNA constructs exhibited highly divergent efficiencies in mediating EGFR knockdown and cell killing, coexpression of Ago2 resulted in uniform and highly specific target gene suppression and apoptosis in EGFR-dependent cells. Codelivery of Ago2 with shRNA constructs or siRNA duplexes thus provides a strategy to enhance the efficacy and the specificity of RNAi in experimental and potentially therapeutic settings.

Structure and inhibition mechanism of the catalytic domain of human squalene epoxidase.

Squalene epoxidase (SQLE), also known as squalene monooxygenase, catalyzes the stereospecific conversion of squalene to 2,3(S)-oxidosqualene, a key step in cholesterol biosynthesis. SQLE inhibition is targeted for the treatment of hypercholesteremia, cancer, and fungal infections. However, lack of structure-function understanding has hindered further progression of its inhibitors. We have determined the first three-dimensional high-resolution crystal structures of human SQLE catalytic domain with small molecule inhibitors (2.3 Å and 2.5 Å). Comparison with its unliganded state (3.0 Å) reveals conformational rearrangements upon inhibitor binding, thus allowing deeper interpretation of known structure-activity relationships. We use the human SQLE structure to further understand the specificity of terbinafine, an approved agent targeting fungal SQLE, and to provide the structural insights into terbinafine-resistant mutants encountered in the clinic. Collectively, these findings elucidate the structural basis for the specificity of the epoxidation reaction catalyzed by SQLE and enable further rational development of next-generation inhibitors.

A chemical biology screen identifies a vulnerability of neuroendocrine cancer cells to SQLE inhibition.

Aberrant metabolism of cancer cells is well appreciated, but the identification of cancer subsets with specific metabolic vulnerabilities remains challenging. We conducted a chemical biology screen and identified a subset of neuroendocrine tumors displaying a striking pattern of sensitivity to inhibition of the cholesterol biosynthetic pathway enzyme squalene epoxidase (SQLE). Using a variety of orthogonal approaches, we demonstrate that sensitivity to SQLE inhibition results not from cholesterol biosynthesis pathway inhibition, but rather surprisingly from the specific and toxic accumulation of the SQLE substrate, squalene. These findings highlight SQLE as a potential therapeutic target in a subset of neuroendocrine tumors, particularly small cell lung cancers.

Frequent met oncogene amplification in a Brca1/Trp53 mouse model of mammary tumorigenesis.

In a screen for gene copy number alterations in mouse mammary tumors initiated by loss of the Brca1 and Trp53 genes, we observed that the majority (11 of 15; 73%) had high-level amplification of wild-type Met, encoding a growth factor receptor implicated in tumor progression. Met amplification was localized to unstable double minute chromosomes and was uniquely found in mouse breast tumors driven by loss of Brca1 and Trp53. Whereas analogous MET amplification was not found in human breast cancers, the identification of a dominant somatic genetic lesion in the Brca1/Trp53 mouse model suggests that recurrent secondary hits may also exist in BRCA1-initiated human breast cancer.

SUMO-1 modification of the Wilms' tumor suppressor WT1.

SUMO-1 conjugation modulates numerous cellular functions, including the subnuclear localization of its target proteins. The WT1 tumor suppressor encodes a four-zinc finger protein with distinct splicing isoforms. WT1(-KTS), encoding uninterrupted zinc fingers, functions as a transcription factor and has a diffusely nuclear distribution; WT1(+KTS), with an insertion of three amino acids (KTS) between zinc fingers three and four, localizes to discrete nuclear speckles, the function of which is unknown. Because the SUMO-1 E2-conjugating enzyme, Ubc9, interacts with WT1, we tested whether sumoylation modulates the cellular localization of WT1. We find here that both WT1 isoforms are directly sumoylated on lysine residues 73 and 177. Although RNA interference-mediated Ubc9 depletion effectively suppresses WT1 nuclear speckles, a SUMO-1-deficient WT1(+KTS)(K73, 177R) double mutant retains localization to speckles. Thus, direct sumoylation of WT1 is not responsible for its cellular localization, and other sumoylated proteins may target WT1 to these nuclear structures. Identification of other components of WT1-associated speckles is likely to provide clues to their function.

Differential Aspartate Usage Identifies a Subset of Cancer Cells Particularly Dependent on OGDH.

Although aberrant metabolism in tumors has been well described, the identification of cancer subsets with particular metabolic vulnerabilities has remained challenging. Here, we conducted an siRNA screen focusing on enzymes involved in the tricarboxylic acid (TCA) cycle and uncovered a striking range of cancer cell dependencies on OGDH, the E1 subunit of the alpha-ketoglutarate dehydrogenase complex. Using an integrative metabolomics approach, we identified differential aspartate utilization, via the malate-aspartate shuttle, as a predictor of whether OGDH is required for proliferation in 3D culture assays and for the growth of xenograft tumors. These findings highlight an anaplerotic role of aspartate and, more broadly, suggest that differential nutrient utilization patterns can identify subsets of cancers with distinct metabolic dependencies for potential pharmacological intervention.