Pre-Treatment Mutational and Transcriptomic Landscape of Responding Metastatic Melanoma Patients to Anti-PD1 Immunotherapy

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy (ACT), have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA sequencing (RNA-seq) datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

Androgen receptor activity inversely correlates with immune cell infiltration and immunotherapy response across multiple cancer lineages

There is now increasing recognition of the important role of androgen receptor (AR) in modulating immune function. To gain a comprehensive understanding of the effects of AR activity on cancer immunity, we employed a computational approach to profile AR activity in 33 human tumor types using RNA-Seq datasets from The Cancer Genome Atlas. Our pan-cancer analysis revealed that the genes most negatively correlated with AR activity across cancers are involved in active immune system processes. Importantly, we observed a significant negative correlation between AR activity and IFNγ pathway activity at the pan-cancer level. Indeed, using a matched biopsy dataset from subjects with prostate cancer before and after AR-targeted treatment, we verified that inhibiting AR enriches immune cell abundances and is associated with higher IFNγ pathway activity. Furthermore, by analyzing immunotherapy datasets in multiple cancers, our results demonstrate that low AR activity was significantly associated with a favorable response to immunotherapy. Together, our data provide a comprehensive assessment of the relationship between AR signaling and tumor immunity.

A method for predicting drugs that can boost the efficacy of immune checkpoint blockade

Combination therapy is a promising therapeutic strategy to enhance the efficacy of immune checkpoint blockade (ICB); however, predicting drugs for effective combination is challenging. Here we developed a general data-driven method called CM-Drug for screening compounds that can boost ICB treatment efficacy based on core and minor gene sets identified between responsive and nonresponsive samples in ICB therapy. The CM-Drug method was validated using melanoma and lung cancer mouse models, with combined therapeutic efficacy demonstrated in eight of nine predicted compounds. Among these compounds, taltirelin had the strongest synergistic effect. Mechanistic analysis and experimental verification demonstrated that taltirelin can stimulate CD8+ T cells and is mediated by the induction of thyroid-stimulating hormone. This study provides an effective and general method for predicting and evaluating drugs for combination therapy and identifies candidate compounds for future ICB combination therapy.

Single-cell RNA sequencing integrated with bulk RNA sequencing analysis identifies a tumor immune microenvironment-related lncRNA signature in lung adenocarcinoma

BACKGROUND

Recently, long non-coding RNAs (lncRNAs) have been demonstrated as essential roles in tumor immune microenvironments (TIME). Nevertheless, researches on the clinical significance of TIME-related lncRNAs are limited in lung adenocarcinoma (LUAD).

METHODS

Single-cell RNA sequencing and bulk RNA sequencing data are integrated to identify TIME-related lncRNAs. A total of 1368 LUAD patients are enrolled from 6 independent datasets. An integrative machine learning framework is introduced to develop a TIME-related lncRNA signature (TRLS).

RESULTS

This study identified TIME-related lncRNAs from integrated analysis of single‑cell and bulk RNA sequencing data. According to these lncRNAs, a TIME-related lncRNA signature was developed and validated from an integrative procedure in six independent cohorts. TRLS exhibited a robust and reliable performance in predicting overall survival. Superior prediction performance barged TRLS to the forefront from comparison with general clinical features, molecular characters, and published signatures. Moreover, patients with low TRLS displayed abundant immune cell infiltration and active lipid metabolism, while patients with high TRLS harbored significant genomic alterations, high PD-L1 expression, and elevated DNA damage repair (DDR) relevance. Notably, subclass mapping analysis of nine immunotherapeutic cohorts demonstrated that patients with high TRLS were more sensitive to immunotherapy.

CONCLUSIONS

This study developed a promising tool based on TIME-related lncRNAs, which might contribute to tailored treatment and prognosis management of LUAD patients.

Single-cell histone chaperones patterns guide intercellular communication of tumor microenvironment that contribute to breast cancer metastases

BACKGROUND

Histone chaperones (HCs) are crucial for governing genome stability and gene expression in multiple cancers. However, the functioning of HCs in the tumor microenvironment (TME) is still not clearly understood.

METHODS

Self-tested single-cell RNA-seq data derived from 6 breast cancer (BC) patients with brain and liver metastases were reanalyzed by nonnegative matrix factorization (NMF) algorithm for 36 HCs. TME subclusters were observed with BC and immunotherapy public cohorts to assess their prognosis and immune response. The biological effect of HSPA8, one of the HCs, was verified by transwell assay and wound-healing assays.

RESULTS

Cells including fibroblasts, macrophages, B cells, and T cells, were classified into various subclusters based on marker genes. Additionally, it showed that HCs might be strongly associated with biological and clinical features of BC metastases, along with the pseudotime trajectory of each TME cell type. Besides, the results of bulk-seq analysis revealed that TME cell subclusters mediated by HCs distinguished significant prognostic value for BC patients and were relevant to patients' immunotherapy responses, especially for B cells and macrophages. In particular, CellChat analysis exhibited that HCs-related TME cell subclusters revealed extensive and diverse interactions with malignant cells. Finally, transwell and wound-healing assays exhibited that HSPA8 deficiency inhibited BC cell migration and invasion.

CONCLUSIONS

Collectively, our study first dissected HCs-guided intercellular communication of TME that contribute to BC metastases.

SASH1 interacts with TNKS2 and promotes human melanocyte stem cell maintenance

Both aging spots (hyperpigmentation) and hair graying (lack of pigmentation) are associated with aging, two seemingly opposite pigmentation phenotypes. It is not clear how they are mechanistically connected. This study investigated the underlying mechanism in a family with an inherited pigmentation disorder. Clinical examinations identified accelerated hair graying and skin dyspigmentation (intermixed hyper and hypopigmentation) in the family members carrying the SASH1 S519N variant. Cell assays indicated that SASH1 promoted stem-like characteristics in human melanocytes, and SASH1 S519N was defective in this function. Multiple assays showed that SASH1 binds to tankyrase 2 (TNKS2), which is required for SASH1's promotion of stem-like function. Further, the SASH1 S519N variant is in a bona fide Tankyrase-binding motif, and SASH1 S519N alters the binding kinetics and affinity. Results here indicate SASH1 as a novel protein regulating the appropriate balance between melanocyte stem cells (McSC) and mature melanocytes (MCs), with S519N variant causing defects. We propose that dysfunction of McSC maintenance connects multiple aging-associated pigmentation phenotypes in the general population.

PDE4D/cAMP/IL-23 axis determines the immunotherapy efficacy of lung adenocarcinoma via activating the IL-9 autocrine loop of cytotoxic T lymphocytes

Although immunotherapy has changed the prognosis of many advanced malignancies including lung adenocarcinoma (LUAD), many patients are insensitive to the drugs, with the mechanisms yet to be elucidated. Herein, we identified PDE4D as an immunotherapy efficacy-related gene through bioinformatics screening. By using a co-culture system of LUAD cells and tumor-cell-specific CD8⁺ T cells, a functional PDE4D/cAMP/IL-23 axis was further revealed in LUAD cells. Fluorescent multiplex immunohistochemistry analysis of patient-derived samples and the in vivo mouse LUAD xenograft tumors revealed not only the colocalization of IL-23 and CD8⁺ T cells but also the immune potentiating effect of IL-23 on cytotoxic T lymphocytes (CTLs) in LUAD tissues. Through transcriptome sequencing and functional validations, IL-23 was proven to up-regulate IL-9 expression in CTLs via activating the NF-κB signaling, leading to elevated productions of immune effector molecules and enhanced efficacy of antitumor immunotherapy. Interestingly, an autocrine loop of IL-9 was also uncovered during this process. In conclusion, PDE4D/cAMP/IL-23 axis determines the immunotherapy efficacy of human LUAD. This effect is mediated by the activation of an NF-κB-dependent IL-9 autocrine loop in CTLs.

Patterns of Somatic Variants in Colorectal Adenoma and Carcinoma Tissue and Matched Plasma Samples from the Hungarian Oncogenome Program

Analysis of circulating cell-free DNA (cfDNA) of colorectal adenoma (AD) and cancer (CRC) patients provides a minimally invasive approach that is able to explore genetic alterations. It is unknown whether there are specific genetic variants that could explain the high prevalence of CRC in Hungary. Whole-exome sequencing (WES) was performed on colon tissues (27 AD, 51 CRC) and matched cfDNAs (17 AD, 33 CRC); furthermore, targeted panel sequencing was performed on a subset of cfDNA samples. The most frequently mutated genes were APC, KRAS, and FBN3 in AD, while APC, TP53, TTN, and KRAS were the most frequently mutated in CRC tissue. Variants in KRAS codons 12 (AD: 8/27, CRC: 11/51 (0.216)) and 13 (CRC: 3/51 (0.06)) were the most frequent in our sample set, with G12V (5/27) dominance in ADs and G12D (5/51 (0.098)) in CRCs. In terms of the cfDNA WES results, tumor somatic variants were found in 6/33 of CRC cases. Panel sequencing revealed somatic variants in 8 out of the 12 enrolled patients, identifying 12/20 tumor somatic variants falling on its targeted regions, while WES recovered only 20% in the respective regions in cfDNA of the same patients. In liquid biopsy analyses, WES is less efficient compared to the targeted panel sequencing with a higher coverage depth that can hold a relevant clinical potential to be applied in everyday practice in the future.

Genome-wide gain-of-function screening characterized lncRNA regulators for tumor immune response

The majority of lncRNAs' roles in tumor immunology remain elusive. This project performed a CRISPR activation screening of 9744 lncRNAs in melanoma cells cocultured with human CD8+ T cells. We identified 16 lncRNAs potentially regulating tumor immune response. Further integrative analysis using tumor immunogenomics data revealed that IL10RB-DT and LINC01198 are significantly correlated with tumor immune response and survival in melanoma and breast cancer. Specifically, IL10RB-DT suppresses CD8+ T cells activation via inhibiting IFN-γ-JAK-STAT1 signaling and antigen presentation in melanoma and breast cancer cells. On the other hand, LINC01198's up-regulation sensitizes the killing of tumor cells by CD8+ T cells. Mechanistically, LINC01198 interacts and activates NF-κB component p65 to trigger the type I and type II interferon responses in melanoma and breast cancer cells. Our study systematically characterized novel lncRNAs involved in tumor immune response.

Integrated analysis from multi-center studies identities m7G-derived modification pattern and risk stratification system in skin cutaneous melanoma

The m7G modification has been proven to play an important role in RNA post-transcriptional modification and protein translation. However, the potential role of m7G modification patterns in assessing the prognosis of Skin cutaneous melanoma (SKCM) and tumor microenvironment (TME) has not been well studied. In this study, we investigated and finally identified 21 available m7G-related genes. We used hierarchical clustering (K-means) to classify 743 SKCM patients into three m7G-modified subtypes named m7G/gene cluster-A, B, C. We found that both m7G cluster B and gene cluster B exhibited higher prognosis and higher immune cell infiltration in TME compared to other subtypes. EIF4E3 and IFIT5, two m7G related genes, were both markedly elevated in Cluster B. Then, we constructed an m7G score system utilizing principal component analysis (PCA) in order to evaluate the patients' prognosis. High m7G score subtype was associated with better survival prognosis and active immune response. Overall, this article revealed that m7G modification patterns were involved in the development of the tumor microenvironment. Evaluating patients' m7G modification patterns will enhance our understanding of TME characteristics and help to guide personal treatment in clinics in the future.

Resistance to anti-PD1 therapies in patients with advanced melanoma: systematic literature review and application of the Society for Immunotherapy of Cancer Immunotherapy Resistance Taskforce anti-PD1 resistance definitions

Nearly half of advanced melanoma patients do not achieve a clinical response with anti-programmed cell death 1 protein (PD1) therapy (i.e. primary resistance) or initially achieve a clinical response but eventually progress during or following further treatment (i.e. secondary resistance). A consensus definition for tumor resistance to anti-PD1 monotherapy was published by Society for Immunotherapy of Cancer Immunotherapy Resistance Taskforce (SITC) in 2020. A systematic literature review (SLR) of clinical trials and observational studies was conducted to characterize the proportions of advanced melanoma patients who have progressed on anti-PD1 therapies. The SLR included 55 unique studies and the SITC definition of primary resistance was applied to 37 studies that specified disease progression by best overall response. Median and range of patients with primary resistance in studies that specified first-line and second-line or higher anti-PD1 monotherapy was 35.50% (21.19-39.13%; n = 4 studies) and 41.54% (30.00-56.41%, n = 3 studies); median and range of patients with primary resistance in studies that specified first-line and second-line or higher combination therapy was 30.23% (15.79-33.33%; n = 6 studies), and 70.00% (61.10-73.33%; n = 3 studies). Primary resistance to anti-PD1 monotherapies and when in combination with ipilimumab are higher in patients receiving second-line or higher therapies, in patients with acral, mucosal, and uveal melanoma, and in patients with active brain metastases. The percentage of patients with primary resistance was generally consistent across clinical trials, with variability in resistance noted for observational studies. Limitations include applying the SITC definitions to combination therapies, where consensus definitions are not yet available. Future studies should highly consider utilizing the SITC definitions to harmonize how resistance is classified and facilitate meaningful context for clinical activity.

Organ Specific Copy Number Variations in Visceral Metastases of Human Melanoma

Malignant melanoma is one of the most aggressive skin cancers with high potential of visceral dissemination. Since the information about melanoma genomics is mainly based on primary tumors and lymphatic or skin metastases, an autopsy-based visceral metastasis biobank was established. We used copy number variation arrays (N = 38 samples) to reveal organ specific alterations. Results were partly completed by proteomic analysis. A significant increase of high-copy number gains was found in an organ-specific manner, whereas copy number losses were predominant in brain metastases, including the loss of numerous DNA damage response genes. Amplification of many immune genes was also observed, several of them are novel in melanoma, suggesting that their ectopic expression is possibly underestimated. This "immunogenic mimicry" was exclusive for lung metastasis. We also provided evidence for the possible autocrine activation of c-MET, especially in brain and lung metastases. Furthermore, frequent loss of 9p21 locus in brain metastases may predict higher metastatic potential to this organ. Finally, a significant correlation was observed between BRAF gene copy number and mutant allele frequency, mainly in lung metastases. All of these events may influence therapy efficacy in an organ specific manner, which knowledge may help in alleviating difficulties caused by resistance.

Melanoma Immunotherapy and Precision Medicine in the Era of Tumor Micro-Tissue Engineering: Where Are We Now and Where Are We Going?

Malignant melanoma still remains a cancer with very poor survival rates, although it is at the forefront of personalized medicine. Most patients show partial responses and disease progressed due to adaptative resistance mechanisms, preventing long-lasting clinical benefits to the current treatments. The response to therapies can be shaped by not only taking into account cancer cell heterogeneity and plasticity, but also by its structural context as well as the cellular component of the tumor microenvironment (TME). Here, we review the recent development in the field of immunotherapy and target-based therapy and how, in the era of tumor micro-tissue engineering, ex-vivo assays could help to enhance our melanoma biology knowledge in its complexity, translating it in the development of successful therapeutic strategies, as well as in the prediction of therapeutic benefits.

CTR-DB, an omnibus for patient-derived gene expression signatures correlated with cancer drug response

To date, only some cancer patients can benefit from chemotherapy and targeted therapy. Drug resistance continues to be a major and challenging problem facing current cancer research. Rapidly accumulated patient-derived clinical transcriptomic data with cancer drug response bring opportunities for exploring molecular determinants of drug response, but meanwhile pose challenges for data management, integration, and reuse. Here we present the Cancer Treatment Response gene signature DataBase (CTR-DB, http://ctrdb.ncpsb.org.cn/), a unique database for basic and clinical researchers to access, integrate, and reuse clinical transcriptomes with cancer drug response. CTR-DB has collected and uniformly reprocessed 83 patient-derived pre-treatment transcriptomic source datasets with manually curated cancer drug response information, involving 28 histological cancer types, 123 drugs, and 5139 patient samples. These data are browsable, searchable, and downloadable. Moreover, CTR-DB supports single-dataset exploration (including differential gene expression, receiver operating characteristic curve, functional enrichment, sensitizing drug search, and tumor microenvironment analyses), and multiple-dataset combination and comparison, as well as biomarker validation function, which provide insights into the drug resistance mechanism, predictive biomarker discovery and validation, drug combination, and resistance mechanism heterogeneity.

Targeting tumor-derived NLRP3 reduces melanoma progression by limiting MDSCs expansion

The nucleotide-binding domain, leucine-rich containing family, pyrin domain-containing-3 (NLRP3) inflammasome, an intracellular complex that regulates maturation and release of interleukin (IL)-1β, is active in biopsies of metastatic melanoma. Here, we demonstrate that NLRP3 activation in melanoma cells drives tumor progression in mice. Subsequent to NLRP3 activation in melanoma cells, IL-1β induces melanoma-associated inflammation, resulting in immunosuppression. Oral administration of a single NLRP3 inhibitor (OLT1177) reduces melanoma growth and melanoma-associated myeloid-derived suppressor cell expansion. Inhibition of the NLRP3 signaling in combination with anti–PD-1 revealed augmented efficacy compared to monotherapy. These data propose that NLRP3 is a therapeutic target for human melanoma.

Interleukin-1β (IL-1β)–mediated inflammation suppresses antitumor immunity, leading to the generation of a tumor-permissive environment, tumor growth, and progression. Here, we demonstrate that nucleotide-binding domain, leucine-rich containing family, pyrin domain-containing-3 (NLRP3) inflammasome activation in melanoma is linked to IL-1β production, inflammation, and immunosuppression. Analysis of cancer genome datasets (TCGA and GTEx) revealed greater NLRP3 and IL-1β expression in cutaneous melanoma samples (n = 469) compared to normal skin (n = 324), with a highly significant correlation between NLRP3 and IL-1β (P < 0.0001). We show the formation of the NLRP3 inflammasome in biopsies of metastatic melanoma using fluorescent resonance energy transfer analysis for NLRP3 and apoptosis-associated speck-like protein containing a CARD. In vivo, tumor-associated NLRP3/IL-1 signaling induced expansion of myeloid-derived suppressor cells (MDSCs), leading to reduced natural killer and CD8⁺ T cell activity concomitant with an increased presence of regulatory T (Treg) cells in the primary tumors. Either genetic or pharmacological inhibition of tumor-derived NLRP3 by dapansutrile (OLT1177) was sufficient to reduce MDSCs expansion and to enhance antitumor immunity, resulting in reduced tumor growth. Additionally, we observed that the combination of NLRP3 inhibition and anti–PD-1 treatment significantly increased the antitumor efficacy of the monotherapy by limiting MDSC-mediated T cell suppression and tumor progression. These data show that NLRP3 activation in melanoma cells is a protumor mechanism, which induces MDSCs expansion and immune evasion. We conclude that inhibition of NLRP3 can augment the efficacy of anti–PD-1 therapy.

Precision Medicine and Melanoma: Multi-Omics Approaches to Monitoring the Immunotherapy Response

The treatment and management of patients with metastatic melanoma have evolved considerably in the "era" of personalized medicine. Melanoma was one of the first solid tumors to benefit from immunotherapy; life expectancy for patients in advanced stage of disease has improved. However, many progresses have yet to be made considering the (still) high number of patients who do not respond to therapies or who suffer adverse events. In this scenario, precision medicine appears fundamental to direct the most appropriate treatment to the single patient and to guide towards treatment decisions. The recent multi-omics analyses (genomics, transcriptomics, proteomics, metabolomics, radiomics, etc.) and the technological evolution of data interpretation have allowed to identify and understand several processes underlying the biology of cancer; therefore, improving the tumor clinical management. Specifically, these approaches have identified new pharmacological targets and potential biomarkers used to predict the response or adverse events to treatments. In this review, we will analyze and describe the most important omics approaches, by evaluating the methodological aspects and progress in melanoma precision medicine.

Investigating Mechanisms of Response or Resistance to Immune Checkpoint Inhibitors by Analyzing Cell-Cell Communications in Tumors Before and After Programmed Cell Death-1 (PD-1) Targeted Therapy: An Integrative Analysis Using Single-cell RNA and Bulk-RNA Sequencing Data

Currently, a significant proportion of cancer patients do not benefit from programmed cell death-1 (PD-1)-targeted therapy. Overcoming drug resistance remains a challenge. In this study, single-cell RNA sequencing and bulk RNA sequencing data from samples collected before and after anti-PD-1 therapy were analyzed. Cell-cell interaction analyses were performed to determine the differences between pretreatment responders and nonresponders and the relative differences in changes from pretreatment to posttreatment status between responders and nonresponders to ultimately investigate the specific mechanisms underlying response and resistance to anti-PD-1 therapy. Bulk-RNA sequencing data were used to validate our results. Furthermore, we analyzed the evolutionary trajectory of ligands/receptors in specific cell types in responders and nonresponders. Based on pretreatment data from responders and nonresponders, we identified several different cell-cell interactions, like WNT5A-PTPRK, EGFR-AREG, AXL-GAS6 and ACKR3-CXCL12. Furthermore, relative differences in the changes from pretreatment to posttreatment status between responders and nonresponders existed in SELE-PSGL-1, CXCR3-CCL19, CCL4-SLC7A1, CXCL12-CXCR3, EGFR-AREG, THBS1-a3b1 complex, TNF-TNFRSF1A, TNF-FAS and TNFSF10-TNFRSF10D interactions. In trajectory analyses of tumor-specific exhausted CD8 T cells using ligand/receptor genes, we identified a cluster of T cells that presented a distinct pattern of ligand/receptor expression. They highly expressed suppressive genes like HAVCR2 and KLRC1, cytotoxic genes like GZMB and FASLG and the tissue-residence-related gene CCL5. These cells had increased expression of survival-related and tissue-residence-related genes, like heat shock protein genes and the interleukin-7 receptor (IL-7R), CACYBP and IFITM3 genes, after anti-PD-1 therapy. These results reveal the mechanisms underlying anti-PD-1 therapy response and offer abundant clues for potential strategies to improve immunotherapy.

NF-κB in Cancer Immunity: Friend or Foe?

The emergence of immunotherapies has definitely proven the tight relationship between malignant and immune cells, its impact on cancer outcome and its therapeutic potential. In this context, it is undoubtedly critical to decipher the transcriptional regulation of these complex interactions. Following early observations demonstrating the roles of NF-κB in cancer initiation and progression, a series of studies converge to establish NF-κB as a master regulator of immune responses to cancer. Importantly, NF-κB is a family of transcriptional activators and repressors that can act at different stages of cancer immunity. In this review, we provide an overview of the selective cell-intrinsic contributions of NF-κB to the distinct cell types that compose the tumor immune environment. We also propose a new view of NF-κB targeting drugs as a new class of immunotherapies for cancer.