Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade.

Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.

Zebrafish ApoB-Containing Lipoprotein Metabolism: A Closer Look.

Zebrafish have become a powerful model of mammalian lipoprotein metabolism and lipid cell biology. Most key proteins involved in lipid metabolism, including cholesteryl ester transfer protein, are conserved in zebrafish. Consequently, zebrafish exhibit a human-like lipoprotein profile. Zebrafish with mutations in genes linked to human metabolic diseases often mimic the human phenotype. Zebrafish larvae develop rapidly and externally around the maternally deposited yolk. Recent work revealed that any disturbance of lipoprotein formation leads to the accumulation of cytoplasmic lipid droplets and an opaque yolk, providing a visible phenotype to investigate disturbances of the lipoprotein pathway, already leading to discoveries in MTTP (microsomal triglyceride transfer protein) and ApoB (apolipoprotein B). By 5 days of development, the digestive system is functional, making it possible to study fluorescently labeled lipid uptake in the transparent larvae. These and other approaches enabled the first in vivo description of the STAB (stabilin) receptors, showing lipoprotein uptake in endothelial cells. Various zebrafish models have been developed to mimic human diseases by mutating genes known to influence lipoproteins (eg, ldlra, apoC2). This review aims to discuss the most recent research in the zebrafish ApoB-containing lipoprotein and lipid metabolism field. We also summarize new insights into lipid processing within the yolk cell and how changes in lipid flux alter yolk opacity. This curious new finding, coupled with the development of several techniques, can be deployed to identify new players in lipoprotein research directly relevant to human disease.

Molecular profiling of reticular gigantocellularis neurons indicates that eNOS modulates environmentally dependent levels of arousal.

Neurons of the medullary reticular nucleus gigantocellularis (NGC) and their targets have recently been a focus of research on mechanisms supporting generalized CNS arousal (GA) required for proper cognitive functions. Using the retro-TRAP method, we characterized transcripts enriched in NGC neurons which have projections to the thalamus. The unique expression and activation of the endothelial nitric oxide (eNOS) signaling pathway in these cells and their intimate connections with blood vessels indicate that these neurons exert direct neurovascular coupling. Production of nitric oxide (NO) within eNOS-positive NGC neurons increases after environmental perturbations, indicating a role for eNOS/NO in modulating environmentally appropriate levels of GA. Inhibition of NO production causes dysregulated behavioral arousal after exposure to environmental perturbation. Further, our findings suggest interpretations for associations between psychiatric disorders and mutations in the eNOS locus.

Genome-wide prediction of synthetic rescue mediators of resistance to targeted and immunotherapy.

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

Directly Measuring Atherogenic Lipoprotein Kinetics in Zebrafish with the Photoconvertible LipoTimer Reporter.

Lipoprotein kinetics are a crucial factor in understanding lipoprotein metabolism since a prolonged time in circulation can contribute to the atherogenic character of apolipoprotein-B (ApoB)-containing lipoproteins (B-lps). Here, we report a method to directly measure lipoprotein kinetics in live developing animals. We developed a zebrafish geneticly encoded reporter, LipoTimer, in which endogenous ApoBb.1 is fused to the photoconvertible fluorophore Dendra2 which shift its emission profile from green to red upon UV exposure. By quantifying the red population of ApoB-Dendra2 over time, we found that B-lp turnover in wild-type larvae becomes faster as development proceeds. Mutants with impaired B-lp uptake or lipolysis present with increased B-lp levels and half-life. In contrast, mutants with impaired B-lp triglyceride loading display slightly fewer and smaller-B-lps, which have a significantly shorter B-lp half-life. Further, we showed that chronic high-cholesterol feeding is associated with a longer B-lp half-life in wild-type juveniles but does not lead to changes in B-lp half-life in lipolysis deficient apoC2 mutants. These data support the hypothesis that B-lp lipolysis is suppressed by the flood of intestinal-derived B-lps that follow a high-fat meal.

Pla2g12b drives expansion of triglyceride-rich lipoproteins.

Vertebrates transport hydrophobic triglycerides through the circulatory system by packaging them within amphipathic particles called Triglyceride-Rich Lipoproteins. Yet, it remains largely unknown how triglycerides are loaded onto these particles. Mutations in Phospholipase A2 group 12B (PLA2G12B) are known to disrupt lipoprotein homeostasis, but its mechanistic role in this process remains unclear. Here we report that PLA2G12B channels lipids within the lumen of the endoplasmic reticulum into nascent lipoproteins. This activity promotes efficient lipid secretion while preventing excess accumulation of intracellular lipids. We characterize the functional domains, subcellular localization, and interacting partners of PLA2G12B, demonstrating that PLA2G12B is calcium-dependent and tightly associated with the membrane of the endoplasmic reticulum. We also detect profound resistance to atherosclerosis in PLA2G12B mutant mice, suggesting an evolutionary tradeoff between triglyceride transport and cardiovascular disease risk. Here we identify PLA2G12B as a key driver of triglyceride incorporation into vertebrate lipoproteins.

Fish Playpens - Method for raising individual juvenile zebrafish on a recirculating system for studies requiring repeated measures.

Even though many experimental approaches benefit from tracking individual larval animals, there is yet to be a commercial zebrafish rack system designed to accomplish this task. Thus, we invented playpens, an acrylic and screen container, to raise 12 individual zebrafish juveniles per standard 10 L tank on an existing recirculating fish system. During a week-long experiment, fish raised in playpens grow to the same size as conventionally raised juveniles.

Fish Playpens: Method for Raising Individual Juvenile Zebrafish on a Recirculating System for Studies Requiring Repeated Measures.

Even though many experimental approaches benefit from tracking individual juvenile animals, there is yet to be a commercial zebrafish rack system designed to accomplish this task. Thus, we invented playpens, an acrylic, and screen container, to raise 12 individual zebrafish juveniles per standard 10 L tank on an existing recirculating fish system. During a week-long experiment, fish raised in playpens grow to the same size as conventionally raised juveniles.

A high-cholesterol zebrafish diet promotes hypercholesterolemia and fasting-associated liver triglycerides accumulation.

Zebrafish are an ideal model organism to study lipid metabolism and to elucidate the molecular underpinnings of human lipid-associated disorders. In this study, we provide an improved protocol to assay the impact of a high-cholesterol diet (HCD) on zebrafish lipid deposition and lipoprotein regulation. Fish fed HCD developed hypercholesterolemia as indicated by significantly elevated ApoB-containing lipoproteins (ApoB-LP) and increased plasma levels of cholesterol and cholesterol esters. Feeding of the HCD to larvae (8 days followed by a 1 day fast) and adult female fish (2 weeks, followed by 3 days of fasting) was also associated with a fatty liver phenotype that presented as severe hepatic steatosis. The HCD feeding paradigm doubled the levels of liver triacylglycerol (TG), which was striking because our HCD was only supplemented with cholesterol. The accumulated liver TG was unlikely due to increased de novo lipogenesis or inhibited ß-oxidation since no differentially expressed genes in these pathways were found between the livers of fish fed the HCD versus control diets. However, fasted HCD fish had significantly increased lipogenesis gene fasn in adipose tissue and higher free fatty acids (FFA) in plasma. This suggested that elevated dietary cholesterol resulted in lipid accumulation in adipocytes, which supplied more FFA during fasting, promoting hepatic steatosis. In conclusion, our HCD zebrafish protocol represents an effective and reliable approach for studying the temporal characteristics of the physiological and biochemical responses to high levels of dietary cholesterol and provides insights into the mechanisms that may underlie fatty liver disease.

Pathway signatures derived from on-treatment tumor specimens predict response to anti-PD1 blockade in metastatic melanoma.

Both genomic and transcriptomic signatures have been developed to predict responses of metastatic melanoma to immune checkpoint blockade (ICB) therapies; however, most of these signatures are derived from pre-treatment biopsy samples. Here, we build pathway-based super signatures in pre-treatment (PASS-PRE) and on-treatment (PASS-ON) tumor specimens based on transcriptomic data and clinical information from a large dataset of metastatic melanoma treated with anti-PD1-based therapies as the training set. Both PASS-PRE and PASS-ON signatures are validated in three independent datasets of metastatic melanoma as the validation set, achieving area under the curve (AUC) values of 0.45-0.69 and 0.85-0.89, respectively. We also combine all test samples and obtain AUCs of 0.65 and 0.88 for PASS-PRE and PASS-ON signatures, respectively. When compared with existing signatures, the PASS-ON signature demonstrates more robust and superior predictive performance across all four datasets. Overall, we provide a framework for building pathway-based signatures that is highly and accurately predictive of response to anti-PD1 therapies based on on-treatment tumor specimens. This work would provide a rationale for applying pathway-based signatures derived from on-treatment tumor samples to predict patients' therapeutic response to ICB therapies.

Evolution of delayed resistance to immunotherapy in a melanoma responder.

Despite initial responses1-3, most melanoma patients develop resistance4 to immune checkpoint blockade (ICB). To understand the evolution of resistance, we studied 37 tumor samples over 9 years from a patient with metastatic melanoma with complete clinical response to ICB followed by delayed recurrence and death. Phylogenetic analysis revealed co-evolution of seven lineages with multiple convergent, but independent resistance-associated alterations. All recurrent tumors emerged from a lineage characterized by loss of chromosome 15q, with post-treatment clones acquiring additional genomic driver events. Deconvolution of bulk RNA sequencing and highly multiplexed immunofluorescence (t-CyCIF) revealed differences in immune composition among different lineages. Imaging revealed a vasculogenic mimicry phenotype in NGFRhi tumor cells with high PD-L1 expression in close proximity to immune cells. Rapid autopsy demonstrated two distinct NGFR spatial patterns with high polarity and proximity to immune cells in subcutaneous tumors versus a diffuse spatial pattern in lung tumors, suggesting different roles of this neural-crest-like program in different tumor microenvironments. Broadly, this study establishes a high-resolution map of the evolutionary dynamics of resistance to ICB, characterizes a de-differentiated neural-crest tumor population in melanoma immunotherapy resistance and describes site-specific differences in tumor-immune interactions via longitudinal analysis of a patient with melanoma with an unusual clinical course.

Targeting Extracellular Matrix Remodeling Restores BRAF Inhibitor Sensitivity in BRAFi-resistant Melanoma.

PURPOSE: The extracellular matrix (ECM) is an intriguing, yet understudied component of therapy resistance. Here, we investigated the role of ECM remodeling by the collagenase, MT1-MMP, in conferring resistance of v-Raf murine sarcoma viral oncogene homolog B1 (BRAF)-mutant melanoma to BRAF inhibitor (BRAFi) therapy. EXPERIMENTAL DESIGN: Publicly available RNA-sequencing data and reverse phase protein array were used to determine the relevance of MT1-MMP upregulation in BRAFi-resistant melanoma in patients, patient-derived xenografts, and cell line-derived tumors. Short hairpin RNA (shRNA)-mediated knockdown of MT1-MMP, inhibition via the selective MT1-MMP/MMP2 inhibitor, ND322, or overexpression of MT1-MMP was used to assess the role of MT1-MMP in mediating resistance to BRAFi. RESULTS: MT1-MMP was consistently upregulated in posttreatment tumor samples derived from patients upon disease progression and in melanoma xenografts and cell lines that acquired resistance to BRAFi. shRNA- or ND322-mediated inhibition of MT1-MMP synergized with BRAFi leading to resensitization of resistant cells and tumors to BRAFi. The resistant phenotype depends on the ability of cells to cleave the ECM. Resistant cells seeded in MT1-MMP uncleavable matrixes were resensitized to BRAFi similarly to MT1-MMP inhibition. This is due to the inability of cells to activate integrinß1 (ITGB1)/FAK signaling, as restoration of ITGB1 activity is sufficient to maintain resistance to BRAFi in the context of MT1-MMP inhibition. Finally, the increase in MT1-MMP in BRAFi-resistant cells is TGFß dependent, as inhibition of TGFß receptors I/II dampens MT1-MMP overexpression and restores sensitivity to BRAF inhibition. CONCLUSIONS: BRAF inhibition results in a selective pressure toward higher expression of MT1-MMP. MT1-MMP is pivotal to an ECM-based signaling pathway that confers resistance to BRAFi therapy.

Publisher Correction: Robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma.

In the version of this article originally published, there was an error in the URL linked to by an accession code in the data availability section of the methods. The erroneous URL was: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE100351 . The correct URL is: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE115821 . The error has been corrected in the HTML and PDF versions of this article.

Robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma.

Immune checkpoint blockade (ICB) therapy provides remarkable clinical gains and has been very successful in treatment of melanoma. However, only a subset of patients with advanced tumors currently benefit from ICB therapies, which at times incur considerable side effects and costs. Constructing predictors of patient response has remained a serious challenge because of the complexity of the immune response and the shortage of large cohorts of ICB-treated patients that include both 'omics' and response data. Here we build immuno-predictive score (IMPRES), a predictor of ICB response in melanoma which encompasses 15 pairwise transcriptomics relations between immune checkpoint genes. It is based on two key conjectures: (i) immune mechanisms underlying spontaneous regression in neuroblastoma can predict melanoma response to ICB, and (ii) key immune interactions can be captured via specific pairwise relations of the expression of immune checkpoint genes. IMPRES is validated on nine published datasets1-6 and on a newly generated dataset with 31 patients treated with anti-PD-1 and 10 with anti-CTLA-4, spanning 297 samples in total. It achieves an overall accuracy of AUC = 0.83, outperforming existing predictors and capturing almost all true responders while misclassifying less than half of the nonresponders. Future studies are warranted to determine the value of the approach presented here in other cancer types.