A CRISPRi/a screening platform to study cellular nutrient transport in diverse microenvironments

Blocking the import of nutrients essential for cancer cell proliferation represents a therapeutic opportunity, but it is unclear which transporters to target. Here we report a CRISPR interference/activation screening platform to systematically interrogate the contribution of nutrient transporters to support cancer cell proliferation in environments ranging from standard culture media to tumours. We applied this platform to identify the transporters of amino acids in leukaemia cells and found that amino acid transport involves high bidirectional flux dependent on the microenvironment composition. While investigating the role of transporters in cystine starved cells, we uncovered a role for serotonin uptake in preventing ferroptosis. Finally, we identified transporters essential for cell proliferation in subcutaneous tumours and found that levels of glucose and amino acids can restrain proliferation in that environment. This study establishes a framework for systematically identifying critical cellular nutrient transporters, characterizing their function and exploring how the tumour microenvironment impacts cancer metabolism.

Microenvironmental reorganization in brain tumors following radiotherapy and recurrence revealed by hyperplexed immunofluorescence imaging

The tumor microenvironment plays a crucial role in determining response to treatment. This involves a series of interconnected changes in the cellular landscape, spatial organization, and extracellular matrix composition. However, assessing these alterations simultaneously is challenging from a spatial perspective, due to the limitations of current high-dimensional imaging techniques and the extent of intratumoral heterogeneity over large lesion areas. In this study, we introduce a spatial proteomic workflow termed Hyperplexed Immunofluorescence Imaging (HIFI) that overcomes these limitations. HIFI allows for the simultaneous analysis of > 45 markers in fragile tissue sections at high magnification, using a cost-effective high-throughput workflow. We integrate HIFI with machine learning feature detection, graph-based network analysis, and cluster-based neighborhood analysis to analyze the microenvironment response to radiation therapy in a preclinical model of glioblastoma, and compare this response to a mouse model of breast-to-brain metastasis. Here we show that glioblastomas undergo extensive spatial reorganization of immune cell populations and structural architecture in response to treatment, while brain metastases show no comparable reorganization. Our integrated spatial analyses reveal highly divergent responses to radiation therapy between brain tumor models, despite equivalent radiotherapy benefit.

Cell state dependent effects of Bmal1 on melanoma immunity and tumorigenicity

The circadian clock regulator Bmal1 modulates tumorigenesis, but its reported effects are inconsistent. Here, we show that Bmal1 has a context-dependent role in mouse melanoma tumor growth. Loss of Bmal1 in YUMM2.1 or B16-F10 melanoma cells eliminates clock function and diminishes hypoxic gene expression and tumorigenesis, which could be rescued by ectopic expression of HIF1α in YUMM2.1 cells. By contrast, over-expressed wild-type or a transcriptionally inactive mutant Bmal1 non-canonically sequester myosin heavy chain 9 (Myh9) to increase MRTF-SRF activity and AP-1 transcriptional signature, and shift YUMM2.1 cells from a Sox10high to a Sox9high immune resistant, mesenchymal cell state that is found in human melanomas. Our work describes a link between Bmal1, Myh9, mouse melanoma cell plasticity, and tumor immunity. This connection may underlie cancer therapeutic resistance and underpin the link between the circadian clock, MRTF-SRF and the cytoskeleton.

High-plex immunofluorescence imaging and traditional histology of the same tissue section for discovering image-based biomarkers

Precision medicine is critically dependent on better methods for diagnosing and staging disease and predicting drug response. Histopathology using hematoxylin and eosin (H&E)-stained tissue (not genomics) remains the primary diagnostic method in cancer. Recently developed highly multiplexed tissue imaging methods promise to enhance research studies and clinical practice with precise, spatially resolved single-cell data. Here, we describe the 'Orion' platform for collecting H&E and high-plex immunofluorescence images from the same cells in a whole-slide format suitable for diagnosis. Using a retrospective cohort of 74 colorectal cancer resections, we show that immunofluorescence and H&E images provide human experts and machine learning algorithms with complementary information that can be used to generate interpretable, multiplexed image-based models predictive of progression-free survival. Combining models of immune infiltration and tumor-intrinsic features achieves a 10- to 20-fold discrimination between rapid and slow (or no) progression, demonstrating the ability of multimodal tissue imaging to generate high-performance biomarkers.

Polymerase θ inhibition activates the cGAS-STING pathway and cooperates with immune checkpoint blockade in models of BRCA-deficient cancer

Recently developed inhibitors of polymerase theta (POLθ) have demonstrated synthetic lethality in BRCA-deficient tumor models. To examine the contribution of the immune microenvironment to antitumor efficacy, we characterized the effects of POLθ inhibition in immunocompetent models of BRCA1-deficient triple-negative breast cancer (TNBC) or BRCA2-deficient pancreatic ductal adenocarcinoma (PDAC). We demonstrate that genetic POLQ depletion or pharmacological POLθ inhibition induces both innate and adaptive immune responses in these models. POLθ inhibition resulted in increased micronuclei, cGAS/STING pathway activation, type I interferon gene expression, CD8⁺ T cell infiltration and activation, local paracrine activation of dendritic cells and upregulation of PD-L1 expression. Depletion of CD8⁺ T cells compromised the efficacy of POLθ inhibition, whereas antitumor effects were augmented in combination with anti-PD-1 immunotherapy. Collectively, our findings demonstrate that POLθ inhibition induces immune responses in a cGAS/STING-dependent manner and provide a rationale for combining POLθ inhibition with immune checkpoint blockade for the treatment of HR-deficient cancers.

Heat-inactivated modified vaccinia virus Ankara boosts Th1 cellular and humoral immunity as a vaccine adjuvant

Protein or peptide-based subunit vaccines have generated excitement and renewed interest in combating human cancer or COVID-19 outbreak. One major concern for subunit vaccine application is the weak immune responses induced by protein or peptides. Developing novel and effective vaccine adjuvants are critical for the success of subunit vaccines. Here we explored the potential of heat-inactivated MVA (heat-iMVA) as a vaccine adjuvant. Heat-iMVA dramatically enhances T cell responses and antibodies responses, mainly toward Th1 immune responses when combined with protein or peptide-based immunogen. The adjuvant effect of Heat-iMVA is stronger than live MVA and is dependent on the cGAS/STING-mediated cytosolic DNA-sensing pathway. In a therapeutic vaccination model based on tumor neoantigen peptide vaccine, Heat-iMVA significantly extended the survival and delayed tumor growth. When combined with SARS-CoV-2 spike protein, Heat-iMVA induced more robust spike-specific antibody production and more potent neutralization antibodies. Our results support that Heat-iMVA can be developed as a safe and potent vaccine adjuvant for subunit vaccines against cancer or SARS-CoV-2.

Isoform-resolved mRNA profiling of ribosome load defines interplay of HIF and mTOR dysregulation in kidney cancer

Hypoxia inducible factor (HIF) and mammalian target of rapamycin (mTOR) pathways orchestrate responses to oxygen and nutrient availability. These pathways are frequently dysregulated in cancer, but their interplay is poorly understood, in part because of difficulties in simultaneous measurement of global and mRNA-specific translation. Here, we describe a workflow for measurement of ribosome load of mRNAs resolved by their transcription start sites (TSSs). Its application to kidney cancer cells reveals extensive translational reprogramming by mTOR, strongly affecting many metabolic enzymes and pathways. By contrast, global effects of HIF on translation are limited, and we do not observe reported translational activation by HIF2A. In contrast, HIF-dependent alterations in TSS usage are associated with robust changes in translational efficiency in a subset of genes. Analyses of the interplay of HIF and mTOR reveal that specific classes of HIF1A and HIF2A transcriptional target gene manifest different sensitivity to mTOR, in a manner that supports combined use of HIF2A and mTOR inhibitors in treatment of kidney cancer.

Comparative optimization of combinatorial CRISPR screens

Combinatorial CRISPR technologies have emerged as a transformative approach to systematically probe genetic interactions and dependencies of redundant gene pairs. However, the performance of different functional genomic tools for multiplexing sgRNAs vary widely. Here, we generate and benchmark ten distinct pooled combinatorial CRISPR libraries targeting paralog pairs to optimize digenic knockout screens. Libraries composed of dual Streptococcus pyogenes Cas9 (spCas9), orthogonal spCas9 and Staphylococcus aureus (saCas9), and enhanced Cas12a from Acidaminococcus were evaluated. We demonstrate a combination of alternative tracrRNA sequences from spCas9 consistently show superior effect size and positional balance between the sgRNAs as a robust combinatorial approach to profile genetic interactions of multiple genes.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Integrative molecular and clinical profiling of acral melanoma links focal amplification of 22q11.21 to metastasis

Acral melanoma, the most common melanoma subtype among non-White individuals, is associated with poor prognosis. However, its key molecular drivers remain obscure. Here, we perform integrative genomic and clinical profiling of acral melanomas from 104 patients treated in North America (n = 37) or China (n = 67). We find that recurrent, late-arising focal amplifications of cytoband 22q11.21 are a leading determinant of inferior survival, strongly associated with metastasis, and linked to downregulation of immunomodulatory genes associated with response to immune checkpoint blockade. Unexpectedly, LZTR1 - a known tumor suppressor in other cancers - is a key candidate oncogene in this cytoband. Silencing of LZTR1 in melanoma cell lines causes apoptotic cell death independent of major hotspot mutations or melanoma subtypes. Conversely, overexpression of LZTR1 in normal human melanocytes initiates processes associated with metastasis, including anchorage-independent growth, formation of spheroids, and an increase in MAPK and SRC activities. Our results provide insights into the etiology of acral melanoma and implicate LZTR1 as a key tumor promoter and therapeutic target.

Isoform specific anti-TGFβ therapy enhances antitumor efficacy in mouse models of cancer

TGFβ is a potential target in cancer treatment due to its dual role in tumorigenesis and homeostasis. However, the expression of TGFβ and its inhibition within the tumor microenvironment has mainly been investigated in stroma-heavy tumors. Using B16 mouse melanoma and CT26 colon carcinoma as models of stroma-poor tumors, we demonstrate that myeloid/dendritic cells are the main sources of TGFβ1 and TGFβ3. Depending on local expression of TGFβ isoforms, isoform specific inhibition of either TGFβ1 or TGFβ3 may be effective. The TGFβ signature of CT26 colon carcinoma is defined by TGFβ1 and TGFβ1 inhibition results in tumor delay; B16 melanoma has equal expression of both isoforms and inhibition of either TGFβ1 or TGFβ3 controls tumor growth. Using T cell functional assays, we show that the mechanism of tumor delay is through and dependent on enhanced CD8+ T cell function. To overcome the local immunosuppressive environment, we found that combining TGFβ inhibition with immune checkpoint blockade results in improved tumor control. Our data suggest that TGFβ inhibition in stroma poor tumors shifts the local immune environment to favor tumor suppression.

Suppression of pancreatic ductal adenocarcinoma growth and metastasis by fibrillar collagens produced selectively by tumor cells

Pancreatic ductal adenocarcinoma (PDAC) has a collagen-rich dense extracellular matrix (ECM) that promotes malignancy of cancer cells and presents a barrier for drug delivery. Data analysis of our published mass spectrometry (MS)-based studies on enriched ECM from samples of progressive PDAC stages reveal that the C-terminal prodomains of fibrillar collagens are partially uncleaved in PDAC ECM, suggesting reduced procollagen C-proteinase activity. We further show that the enzyme responsible for procollagen C-proteinase activity, bone morphogenetic protein1 (BMP1), selectively suppresses tumor growth and metastasis in cells expressing high levels of COL1A1. Although BMP1, as a secreted proteinase, promotes fibrillar collagen deposition from both cancer cells and stromal cells, only cancer-cell-derived procollagen cleavage and deposition suppresses tumor malignancy. These studies reveal a role for cancer-cell-derived fibrillar collagen in selectively restraining tumor growth and suggest stratification of patients based on their tumor epithelial collagen I expression when considering treatments related to perturbation of fibrillar collagens.

Functional significance of U2AF1 S34F mutations in lung adenocarcinomas

The functional role of U2AF1 mutations in lung adenocarcinomas (LUADs) remains incompletely understood. Here, we report a significant co-occurrence of U2AF1 S34F mutations with ROS1 translocations in LUADs. To characterize this interaction, we profiled effects of S34F on the transcriptome-wide distribution of RNA binding and alternative splicing in cells harboring the ROS1 translocation. Compared to its wild-type counterpart, U2AF1 S34F preferentially binds and modulates splicing of introns containing CAG trinucleotides at their 3' splice junctions. The presence of S34F caused a shift in cross-linking at 3' splice sites, which was significantly associated with alternative splicing of skipped exons. U2AF1 S34F induced expression of genes involved in the epithelial-mesenchymal transition (EMT) and increased tumor cell invasion. Finally, S34F increased splicing of the long over the short SLC34A2-ROS1 isoform, which was also associated with enhanced invasiveness. Taken together, our results suggest a mechanistic interaction between mutant U2AF1 and ROS1 in LUAD.

Co-incidence of RCC-susceptibility polymorphisms with HIF cis-acting sequences supports a pathway tuning model of cancer

Emerging evidence suggests that dysregulation of oncogenic pathways requires precise tuning in order for cancer to develop. To test this, we examined the overlap between cis-acting elements of the hypoxia-inducible factor (HIF) pathway and cancer-susceptibility polymorphisms as defined in genome-wide association studies (GWAS). In renal cancer, where HIF is constitutively and un-physiologically activated by mutation of the von Hippel-Lindau tumour suppressor, we observed marked excess overlap, which extended to potential susceptibility polymorphisms that are below the conventional threshold applied in GWAS. In contrast, in other cancers where HIF is upregulated by different mechanisms, including micro-environmental hypoxia, we observed no excess in overlap. Our findings support a 'pathway tuning' model of cancer, whereby precise modulation of multiple outputs of specific, activated pathways is important in oncogenesis. This implies that selective pressures to modulate such pathways operate during cancer development and should focus attempts to identify their nature and consequences.

The polymorphic variant rs1800734 influences methylation acquisition and allele-specific TFAP4 binding in the MLH1 promoter leading to differential mRNA expression

Expression of the mismatch repair gene MutL homolog 1 (MLH1) is silenced in a clinically important subgroup of sporadic colorectal cancers. These cancers exhibit hypermutability with microsatellite instability (MSI) and differ from microsatellite-stable (MSS) colorectal cancers in both prognosis and response to therapies. Loss of MLH1 is usually due to epigenetic silencing with associated promoter methylation; coding somatic mutations rarely occur. Here we use the presence of a colorectal cancer (CRC) risk variant (rs1800734) within the MLH1 promoter to investigate the poorly understood mechanisms of MLH1 promoter methylation and loss of expression. We confirm the association of rs1800734 with MSI+ but not MSS cancer risk in our own data and by meta-analysis. Using sensitive allele-specific detection methods, we demonstrate that MLH1 is the target gene for rs1800734 mediated cancer risk. In normal colon tissue, small allele-specific differences exist only in MLH1 promoter methylation, but not gene expression. In contrast, allele-specific differences in both MLH1 methylation and expression are present in MSI+ cancers. We show that MLH1 transcriptional repression is dependent on DNA methylation and can be reversed by a methylation inhibitor. The rs1800734 allele influences the rate of methylation loss and amount of re-expression. The transcription factor TFAP4 binds to the rs1800734 region but with much weaker binding to the risk than the protective allele. TFAP4 binding is absent on both alleles when promoter methylation is present. Thus we propose that TFAP4 binding shields the protective rs1800734 allele of the MLH1 promoter from BRAF induced DNA methylation more effectively than the risk allele.

NDRG4 promoter hypermethylation is a mechanistic biomarker associated with metastatic progression in breast cancer patients

The risk of developing metastatic disease in breast cancer patients is traditionally predictable based on the number of positive axillary lymph nodes, complemented with additional clinicopathological factors. However, since lymph node-negative patients have a 20-30% probability of developing metastatic disease, lymph node information alone is insufficient to accurately assess individual risk. Molecular approaches, such as multigene expression panels, analyze a set of cancer-related genes that more accurately predict the early risk of metastasis and the treatment response. Here, we present N-Myc downstream-regulated gene 4 (NDRG4) epigenetic silencing as a mechanistic biomarker of metastasis in ductal invasive breast tumors. While aberrant NDRG4 DNA hypermethylation is significantly associated with the development of metastatic disease, downregulation of NDRG4 transcription and protein expression is functionally associated with enhanced lymph node adhesion and cell mobility. Here, we show that epigenetic silencing of NDRG4 modulates integrin signaling by assembling β1-integrins into large punctate clusters at the leading edge of tumor cells to promote an "adhesive switch," decreasing cell adhesion to fibronectin and increasing cell adhesion and migration towards vitronectin, an important component of human lymph nodes. Taken together, our functional and clinical observations suggest that NDRG4 is a potential mechanistic biomarker in breast cancer that is functionally associated with metastatic disease.