Evaluation of Endobronchial Ultrasound-Guided Transbronchial Needle Aspiration (EBUS-TBNA) Samples from Advanced Non-Small Cell Lung Cancer for Whole Genome, Whole Exome and Comprehensive Panel Sequencing.

Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2-3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

The oesophageal adenocarcinoma tumour immune microenvironment dictates outcomes with different modalities of neoadjuvant therapy - results from the AGITG DOCTOR trial and the cancer evolution biobank.

A plateau in treatment effect can be seen for the current 'one-size-fits-all' approach to oesophageal adenocarcinoma (OAC) management using neoadjuvant chemoradiotherapy (nCRT) or chemotherapy (nCT). In OAC, the tumour microenvironment (TME) is largely immunosuppressed, however a subgroup of patients with an immune-inflamed TME exist and show improved outcomes. We aimed to understand the overall immune-based mechanisms underlying treatment responses and patient outcomes in OAC, and in relation to neoadjuvant therapy modality. This study included 107 patients; 68 patients were enrolled in the Australian Gastro-Intestinal Trials Group sponsored DOCTOR Trial, and 38 patients were included from the Cancer Evolution Biobank. Matched pre-treatment and post-treatment tumour biopsies were used to perform multi-modality analysis of the OAC TME including NanoString mRNA expression analysis, multiplex and single colour immunohistochemistry (IHC), and peripheral blood mononuclear cell analysis of tumour-antigen specific T cell responses. Patients with the best clinicopathological outcomes and survival had an immune-inflamed TME enriched with anti-tumour immune cells and pathways. Those with the worst survival showed a myeloid T regulatory cell enriched TME, with decreased CD8+ cell infiltration and increased pro-tumour immune cells. Multiplex IHC analysis identified that high intra-tumoural infiltration of CD8+ cells, and low infiltration with CD163+ cells was associated with improved survival. High tumour core CD8+ T cell infiltration, and a low tumour margin infiltration of CD163+ cells was also associated with improved survival. nCRT showed improved survival compared with nCT for patients with low CD8+, or high CD163+ cell infiltration. Poly-functional T cell responses were seen with tumour-antigen specific T cells. Overall, our study supports the development of personalised therapeutic approaches based on the immune microenvironment in OAC. Patients with an immune-inflamed TME show favourable outcomes regardless of treatment modality. However, in those with an immunosuppressed TME with CD163+ cell infiltration, treatment with nCRT can improve outcomes. Our findings support previous studies into the TME of OAC and with more research, immune based biomarker selection of treatment modality may lead in improved outcomes in this deadly disease.

The clinical utility and costs of whole-genome sequencing to detect cancer susceptibility variants-a multi-site prospective cohort study.

BACKGROUND: Many families and individuals do not meet criteria for a known hereditary cancer syndrome but display unusual clusters of cancers. These families may carry pathogenic variants in cancer predisposition genes and be at higher risk for developing cancer. METHODS: This multi-centre prospective study recruited 195 cancer-affected participants suspected to have a hereditary cancer syndrome for whom previous clinical targeted genetic testing was either not informative or not available. To identify pathogenic disease-causing variants explaining participant presentation, germline whole-genome sequencing (WGS) and a comprehensive cancer virtual gene panel analysis were undertaken. RESULTS: Pathogenic variants consistent with the presenting cancer(s) were identified in 5.1% (10/195) of participants and pathogenic variants considered secondary findings with potential risk management implications were identified in another 9.7% (19/195) of participants. Health economic analysis estimated the marginal cost per case with an actionable variant was significantly lower for upfront WGS with virtual panel ($8744AUD) compared to standard testing followed by WGS ($24,894AUD). Financial analysis suggests that national adoption of diagnostic WGS testing would require a ninefold increase in government annual expenditure compared to conventional testing. CONCLUSIONS: These findings make a case for replacing conventional testing with WGS to deliver clinically important benefits for cancer patients and families. The uptake of such an approach will depend on the perspectives of different payers on affordability.

Performance of tumour microenvironment deconvolution methods in breast cancer using single-cell simulated bulk mixtures.

Cells within the tumour microenvironment (TME) can impact tumour development and influence treatment response. Computational approaches have been developed to deconvolve the TME from bulk RNA-seq. Using scRNA-seq profiling from breast tumours we simulate thousands of bulk mixtures, representing tumour purities and cell lineages, to compare the performance of nine TME deconvolution methods (BayesPrism, Scaden, CIBERSORTx, MuSiC, DWLS, hspe, CPM, Bisque, and EPIC). Some methods are more robust in deconvolving mixtures with high tumour purity levels. Most methods tend to mis-predict normal epithelial for cancer epithelial as tumour purity increases, a finding that is validated in two independent datasets. The breast cancer molecular subtype influences this mis-prediction. BayesPrism and DWLS have the lowest combined numbers of false positives and false negatives, and have the best performance when deconvolving granular immune lineages. Our findings highlight the need for more single-cell characterisation of rarer cell types, and suggest that tumour cell compositions should be considered when deconvolving the TME.

Circulating Tumor DNA: A Promising Biomarker for Predicting Recurrence in Patients with BRAF-Negative Melanoma.

For patients with BRAF wild-type stage III and IV melanoma, there is an urgent clinical need to identify prognostic biomarkers and biomarkers predictive of treatment response. Circulating tumor DNA (ctDNA) is emerging as a blood-based biomarker and has shown promising results for many cancers, including melanoma. The purpose of this study was to identify targetable, tumor-derived mutations in patient blood that may lead to treatment alternatives and improved outcomes for patients with BRAF-negative melanoma. Using a CAncer Personalized Profiling by deep Sequencing (CAPP-seq) pan-cancer gene panel, ctDNA from 150 plasma samples (n = 106 patients) was assessed, including serial blood collections for a subset of patients (n = 16). ctDNA variants were detected in 85% of patients, all in targetable pathways, such as vascular endothelial growth factor receptor, epidermal growth factor receptor, phosphatidylinositol 3-kinase/AKT, Bcl2/mammalian target of rapamycin (mTOR), ALK/MET, and cyclin-dependent kinase 4/6. Patients with stage IV melanoma with low ctDNA concentrations, <10 ng/mL, had significantly better disease-specific survival and progression-free survival. Patients with both a high concentration of ctDNA and any detectable ctDNA variants had the worst prognosis. In addition, these results indicated that longitudinal changes in ctDNA correlated with treatment response and disease progression determined by radiology. This study confirms that ctDNA may be used as a noninvasive liquid biopsy to identify recurrent disease and detect targetable variants in patients with late-stage melanoma.

Generalising uncertainty improves accuracy and safety of deep learning analytics applied to oncology.

Uncertainty estimation is crucial for understanding the reliability of deep learning (DL) predictions, and critical for deploying DL in the clinic. Differences between training and production datasets can lead to incorrect predictions with underestimated uncertainty. To investigate this pitfall, we benchmarked one pointwise and three approximate Bayesian DL models for predicting cancer of unknown primary, using three RNA-seq datasets with 10,968 samples across 57 cancer types. Our results highlight that simple and scalable Bayesian DL significantly improves the generalisation of uncertainty estimation. Moreover, we designed a prototypical metric-the area between development and production curve (ADP), which evaluates the accuracy loss when deploying models from development to production. Using ADP, we demonstrate that Bayesian DL improves accuracy under data distributional shifts when utilising 'uncertainty thresholding'. In summary, Bayesian DL is a promising approach for generalising uncertainty, improving performance, transparency, and safety of DL models for deployment in the real world.

Multi-omic features of oesophageal adenocarcinoma in patients treated with preoperative neoadjuvant therapy.

Oesophageal adenocarcinoma is a poor prognosis cancer and the molecular features underpinning response to treatment remain unclear. We investigate whole genome, transcriptomic and methylation data from 115 oesophageal adenocarcinoma patients mostly from the DOCTOR phase II clinical trial (Australian New Zealand Clinical Trials Registry-ACTRN12609000665235), with exploratory analysis pre-specified in the study protocol of the trial. We report genomic features associated with poorer overall survival, such as the APOBEC mutational and RS3-like rearrangement signatures. We also show that positron emission tomography non-responders have more sub-clonal genomic copy number alterations. Transcriptomic analysis categorises patients into four immune clusters correlated with survival. The immune suppressed cluster is associated with worse survival, enriched with myeloid-derived cells, and an epithelial-mesenchymal transition signature. The immune hot cluster is associated with better survival, enriched with lymphocytes, myeloid-derived cells, and an immune signature including CCL5, CD8A, and NKG7. The immune clusters highlight patients who may respond to immunotherapy and thus may guide future clinical trials.

Intratumoral CD16+ Macrophages Are Associated with Clinical Outcomes of Patients with Metastatic Melanoma Treated with Combination Anti-PD-1 and Anti-CTLA-4 Therapy.

PURPOSE: This study characterizes intratumoral macrophage populations within baseline melanoma biopsies from patients with advanced melanoma who received either anti-PD-1 monotherapy or a combination with anti-CTLA-4. Particularly, FcγRIIIa (CD16)-expressing macrophage densities were investigated for associations with response and progression-free survival. EXPERIMENTAL DESIGN: Patients with advanced melanoma who received either anti-PD-1 monotherapy or combination anti-PD-1 and anti-CTLA-4 were retrospectively identified. Macrophage populations were analyzed within baseline melanoma biopsies via multiplex IHC in relation to treatment outcomes. RESULTS: Patients who responded to combination immune checkpoint inhibitor contained higher CD16+ macrophage densities than those who did not respond (196 vs. 7 cells/mm2; P = 0.0041). There was no diffidence in CD16+ macrophage densities in the PD-1 monotherapy-treated patients based on response (118 vs. 89 cells/mm2; P = 0.29). A significantly longer 3-year progression-free survival was observed in combination-treated patients with high intratumoral densities of CD16+ macrophages compared with those with low densities (87% vs. 42%, P = 0.0056, n = 40). No association was observed in anti-PD-1 monotherapy-treated patients (50% vs. 47%, P = 0.4636, n = 50). Melanoma biopsies with high densities of CD16+ macrophages contained upregulated gene expression of critical T-cell recruiting chemokines (CXCL9, CXCL10, and CXCL11). CONCLUSIONS: Our data demonstrate that tumor microenvironments enriched with CD16+ macrophages are favorable for response to combination anti-PD-1 and anti-CTLA-4 therapy but not anti-PD-1 monotherapy. These data provides a potential biomarker of response for combination immunotherapies in patients with metastatic melanoma. See related commentary by Smithy and Luke, p. 2345.

Efficient detection and monitoring of pediatric brain malignancies with liquid biopsy based on patient-specific somatic mutation screening.

BACKGROUND: Brain cancer is the leading cause of cancer-related death in children. Early detection and serial monitoring are essential for better therapeutic outcomes. Liquid biopsy has recently emerged as a promising approach for detecting these tumors by screening body fluids for the presence of circulating tumor DNA (ctDNA). Here we tested the limits of liquid biopsy using patient-specific somatic mutations to detect and monitor primary and metastatic pediatric brain cancer. METHODS: Somatic mutations were identified in 3 ependymoma, 1 embryonal tumor with multilayered rosettes, 1 central nervous system neuroblastoma, and 7 medulloblastoma patients. The mutations were used as liquid biomarkers for serial assessment of cerebrospinal fluid (CSF) samples using a droplet digital PCR (ddPCR) system. The findings were correlated to the imaging data and clinical assessment to evaluate the utility of the approach for clinical translation. RESULTS: We developed personalized somatic mutation ddPCR assays which we show are highly specific, sensitive, and efficient in detection and monitoring of ctDNA, with a positive correlation between presence of ctDNA, disease course, and clinical outcomes in the majority of patients. CONCLUSIONS: We demonstrate the feasibility and clinical utility of personalized mutation-based liquid biopsy for the surveillance of brain cancer in children. However, even with this specific and sensitive approach, we identified some potential false negative analyses. Overall, our results indicate that changes in ctDNA profiles over time demonstrate the great potential of our specific approach for predicting tumor progression, burden, and response to treatment.

IL-10-producing Th1 cells possess a distinct molecular signature in malaria

Control of intracellular parasites responsible for malaria requires host IFN-γ+T-bet+CD4+ T cells (Th1 cells) with IL-10 produced by Th1 cells to mitigate the pathology induced by this inflammatory response. However, these IL-10-producing Th1 (induced type I regulatory [Tr1]) cells can also promote parasite persistence or impair immunity to reinfection or vaccination. Here, we identified molecular and phenotypic signatures that distinguished IL-10-Th1 cells from IL-10+Tr1 cells in Plasmodium falciparum-infected people who participated in controlled human malaria infection studies, as well as C57BL/6 mice with experimental malaria caused by P. berghei ANKA. We also identified a conserved Tr1 cell molecular signature shared between patients with malaria, dengue, and graft-versus-host disease. Genetic manipulation of primary human CD4+ T cells showed that the transcription factor cMAF played an important role in the induction of IL-10, while BLIMP-1 promoted the development of human CD4+ T cells expressing multiple coinhibitory receptors. We also describe heterogeneity of Tr1 cell coinhibitory receptor expression that has implications for targeting these molecules for clinical advantage during infection. Overall, this work provides insights into CD4+ T cell development during malaria that offer opportunities for creation of strategies to modulate CD4+ T cell functions and improve antiparasitic immunity.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes.

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes. SIGNIFICANCE: This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

Increased amphiregulin expression by CD4+ T cells from individuals with asymptomatic Leishmania donovani infection.

Objectives: There is an urgent need to be able to identify individuals with asymptomatic Leishmania donovani infection, so their risk of progressing to VL and transmitting parasites can be managed. This study examined transcriptional markers expressed by CD4+ T cells that could distinguish asymptomatic individuals from endemic controls and visceral leishmaniasis (VL) patients. Methods: CD4+ T cells were isolated from individuals with asymptomatic L. donovani infection, endemic controls and VL patients. RNA was extracted and RNAseq employed to identify differentially expressed genes. The expression of one gene and its protein product during asymptomatic infection were evaluated. Results: Amphiregulin (AREG) was identified as a distinguishing gene product in CD4+ T cells from individuals with asymptomatic L. donovani infection, compared to VL patients and healthy endemic control individuals. AREG levels in plasma and antigen-stimulated whole-blood assay cell culture supernatants were significantly elevated in asymptomatic individuals, compared to endemic controls and VL patients. Regulatory T (Treg) cells were identified as an important source of AREG amongst CD4+ T-cell subsets in asymptomatic individuals. Conclusion: Increased Treg cell AREG expression was identified in individuals with asymptomatic L. donovani infection, suggesting the presence of an ongoing inflammatory response in these individuals required for controlling infection and that AREG may play an important role in preventing inflammation-induced tissue damage and subsequent disease in asymptomatic individuals.

Comprehensive genomic and tumour immune profiling reveals potential therapeutic targets in malignant pleural mesothelioma.

BACKGROUND: Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials. METHODS: We analysed somatic mutations from 229 MPM samples, including previously published data and 58 samples that had undergone WGS within this study. This was combined with RNA-seq analysis to characterize the tumour immune environment. RESULTS: The comprehensive genome analysis identified 12 driver genes, including new candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed SBS5/40 were dominant in 93% of samples, and defects in homologous recombination repair were infrequent in our cohort. The tumour immune environment contained high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immune suppressive environment. The expression of TGFB1 was associated with overall survival. A small subset of samples (less than 10%) had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a 'hot' immune environment independent of the somatic mutations. CONCLUSIONS: We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.

Patient-derived xenograft models capture genomic heterogeneity in endometrial cancer.

BACKGROUND: Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprises four molecular subtypes with differing etiology, prognoses, and responses to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. As pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models. METHODS: Here, we report whole exome or whole genome sequencing of 11 PDX models and their matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC. RESULTS: PDX models were successfully generated from grade 2/3 tumors, including three uterine carcinosarcomas. The models showed similar histomorphology to the primary tumors and represented all four molecular subtypes of EC, including five mismatch-repair deficient models. The different PDX lineages showed a wide range of inter-tumor and intra-tumor heterogeneity. However, for most PDX models, one arm recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease; however, all lacked the HR deficiency genomic signature. CONCLUSIONS: EC PDX models represent the four molecular subtypes of disease and can capture intra-tumor heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi; however, deeper and more durable responses will likely require combination of PARPi with other agents.

Donor bone marrow-derived macrophage MHC II drives neuroinflammation and altered behavior during chronic GVHD in mice.

Graft-versus-host disease (GVHD) remains the leading cause of nonrelapse mortality after allogeneic stem cell transplantation for hematological malignancies. Manifestations of GVHD in the central nervous system (CNS) present as neurocognitive dysfunction in up to 60% of patients; however, the mechanisms driving chronic GVHD (cGVHD) in the CNS are yet to be elucidated. Our studies of murine cGVHD revealed behavioral deficits associated with broad neuroinflammation and persistent Ifng upregulation. By flow cytometry, we observed a proportional shift in the donor-derived T-cell population in the cGVHD brain from early CD8 dominance to later CD4 sequestration. RNA sequencing of the hippocampus identified perturbations to structural and functional synapse-related gene expression, together with the upregulation of genes associated with interferon-γ responses and antigen presentation. Neuroinflammation in the cortex of mice and humans during acute GVHD was recently shown to be mediated by resident microglia-derived tumor necrosis factor. In contrast, infiltration of proinflammatory major histocompatibility complex (MHC) class II+ donor bone marrow (BM)-derived macrophages (BMDMs) was identified as a distinguishing feature of CNS cGVHD. Donor BMDMs, which composed up to 50% of the CNS myeloid population, exhibited a transcriptional signature distinct from resident microglia. Recipients of MHC class II knockout BM grafts exhibited attenuated neuroinflammation and behavior comparable to controls, suggestive of a critical role of donor BMDM MHC class II expression in CNS cGVHD. Our identification of disease mediators distinct from those in the acute phase indicates the necessity to pursue alternative therapeutic targets for late-stage neurological manifestations.

Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance.

We concurrently examine the whole genome, transcriptome, methylome, and immune cell infiltrates in baseline tumors from 77 patients with advanced cutaneous melanoma treated with anti-PD-1 with or without anti-CTLA-4. We show that high tumor mutation burden (TMB), neoantigen load, expression of IFNγ-related genes, programmed death ligand expression, low PSMB8 methylation (therefore high expression), and T cells in the tumor microenvironment are associated with response to immunotherapy. No specific mutation correlates with therapy response. A multivariable model combining the TMB and IFNγ-related gene expression robustly predicts response (89% sensitivity, 53% specificity, area under the curve [AUC], 0.84); tumors with high TMB and a high IFNγ signature show the best response to immunotherapy. This model validates in an independent cohort (80% sensitivity, 59% specificity, AUC, 0.79). Except for a JAK3 loss-of-function mutation, for patients who did not respond as predicted there is no obvious biological mechanism that clearly explained their outlier status, consistent with intratumor and intertumor heterogeneity in response to immunotherapy.

Targeting novel LSD1-dependent ACE2 demethylation domains inhibits SARS-CoV-2 replication.

Treatment options for COVID-19 remain limited, especially during the early or asymptomatic phase. Here, we report a novel SARS-CoV-2 viral replication mechanism mediated by interactions between ACE2 and the epigenetic eraser enzyme LSD1, and its interplay with the nuclear shuttling importin pathway. Recent studies have shown a critical role for the importin pathway in SARS-CoV-2 infection, and many RNA viruses hijack this axis to re-direct host cell transcription. LSD1 colocalized with ACE2 at the cell surface to maintain demethylated SARS-CoV-2 spike receptor-binding domain lysine 31 to promote virus-ACE2 interactions. Two newly developed peptide inhibitors competitively inhibited virus-ACE2 interactions, and demethylase access to significantly inhibit viral replication. Similar to some other predominantly plasma membrane proteins, ACE2 had a novel nuclear function: its cytoplasmic domain harbors a nuclear shuttling domain, which when demethylated by LSD1 promoted importin-α-dependent nuclear ACE2 entry following infection to regulate active transcription. A novel, cell permeable ACE2 peptide inhibitor prevented ACE2 nuclear entry, significantly inhibiting viral replication in SARS-CoV-2-infected cell lines, outperforming other LSD1 inhibitors. These data raise the prospect of post-exposure prophylaxis for SARS-CoV-2, either through repurposed LSD1 inhibitors or new, nuclear-specific ACE2 inhibitors.

Radiomics Biomarkers Correlate with CD8 Expression and Predict Immune Signatures in Melanoma Patients.

Treatment for metastatic melanoma includes targeted and/or immunotherapy. Although many patients respond, only a subset has complete response. As late-stage patients often have multiple tumors in difficult access sites, non-invasive techniques are necessary for the development of predictive/prognostic biomarkers. PET/CT scans from 52 patients with stage III/IV melanoma were assessed and CT image parameters were evaluated as prognostic biomarkers. Analysis indicated patients with high standard deviation or high mean of positive pixels (MPP) had worse progression-free survival (P = 0.00047 and P = 0.0014, respectively) and worse overall survival (P = 0.0223 and P = 0.0465, respectively). Whole-exome sequencing showed high MPP was associated with BRAF mutation status (P = 0.0389). RNA-sequencing indicated patients with immune "cold" signatures had worse survival, which was associated with CT biomarker, MPP4 (P = 0.0284). Multiplex immunofluorescence confirmed a correlation between CD8 expression and image biomarkers (P = 0.0028). IMPLICATIONS: CT parameters have the potential to be cost-effective biomarkers of survival in melanoma, and reflect the tumor immune-microenvironment. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/19/6/950/F1.large.jpg.