Stochastic modelling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones.

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads.

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Exploring tumor-normal cross-talk with TranNet: Role of the environment in tumor progression.

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such dependencies. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights into the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations.

Noncoding RNA circuitry in melanoma onset, plasticity, and therapeutic response.

Melanoma, the cancer of the melanocyte, is the deadliest form of skin cancer with an aggressive nature, propensity to metastasize and tendency to resist therapeutic intervention. Studies have identified that the re-emergence of developmental pathways in melanoma contributes to melanoma onset, plasticity, and therapeutic response. Notably, it is well known that noncoding RNAs play a critical role in the development and stress response of tissues. In this review, we focus on the noncoding RNAs, including microRNAs, long non-coding RNAs, circular RNAs, and other small RNAs, for their functions in developmental mechanisms and plasticity, which drive onset, progression, therapeutic response and resistance in melanoma. Going forward, elucidation of noncoding RNA-mediated mechanisms may provide insights that accelerate development of novel melanoma therapies.

Melanoma clonal subline analysis uncovers heterogeneity-driven immunotherapy resistance mechanisms.

Intratumoral heterogeneity (ITH) can promote cancer progression and treatment failure, but the complexity of the regulatory programs and contextual factors involved complicates its study. To understand the specific contribution of ITH to immune checkpoint blockade (ICB) response, we generated single cell-derived clonal sublines from an ICB-sensitive and genetically and phenotypically heterogeneous mouse melanoma model, M4. Genomic and single cell transcriptomic analyses uncovered the diversity of the sublines and evidenced their plasticity. Moreover, a wide range of tumor growth kinetics were observed in vivo , in part associated with mutational profiles and dependent on T cell-response. Further inquiry into melanoma differentiation states and tumor microenvironment (TME) subtypes of untreated tumors from the clonal sublines demonstrated correlations between highly inflamed and differentiated phenotypes with the response to anti-CTLA-4 treatment. Our results demonstrate that M4 sublines generate intratumoral heterogeneity at both levels of intrinsic differentiation status and extrinsic TME profiles, thereby impacting tumor evolution during therapeutic treatment. These clonal sublines proved to be a valuable resource to study the complex determinants of response to ICB, and specifically the role of melanoma plasticity in immune evasion mechanisms.

Germline modifiers of the tumor immune microenvironment implicate drivers of cancer risk and immunotherapy response.

With the continued promise of immunotherapy for treating cancer, understanding how host genetics contributes to the tumor immune microenvironment (TIME) is essential to tailoring cancer screening and treatment strategies. Here, we study 1084 eQTLs affecting the TIME found through analysis of The Cancer Genome Atlas and literature curation. These TIME eQTLs are enriched in areas of active transcription, and associate with gene expression in specific immune cell subsets, such as macrophages and dendritic cells. Polygenic score models built with TIME eQTLs reproducibly stratify cancer risk, survival and immune checkpoint blockade (ICB) response across independent cohorts. To assess whether an eQTL-informed approach could reveal potential cancer immunotherapy targets, we inhibit CTSS, a gene implicated by cancer risk and ICB response-associated polygenic models; CTSS inhibition results in slowed tumor growth and extended survival in vivo. These results validate the potential of integrating germline variation and TIME characteristics for uncovering potential targets for immunotherapy.

Exploring tumor-normal cross-talk with TranNet: role of the environment in tumor progression.

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such regulation. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights for the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations. The TranNet method was implemented in python, source codes and the data sets used for and generated during this study are available at the Github site https://github.com/ncbi/TranNet .

Broad misappropriation of developmental splicing profile by cancer in multiple organs.

Oncogenesis mimics key aspects of embryonic development. However, the underlying mechanisms are incompletely understood. Here, we demonstrate that the splicing events specifically active during human organogenesis, are broadly reactivated in the organ-specific tumor. Such events are associated with key oncogenic processes and predict proliferation rates in cancer cell lines as well as patient survival. Such events preferentially target nitrosylation and transmembrane-region domains, whose coordinated splicing in multiple genes respectively affect intracellular transport and N-linked glycosylation. We infer critical splicing factors potentially regulating embryonic splicing events and show that such factors are potential oncogenic drivers and are upregulated specifically in malignant cells. Multiple complementary analyses point to MYC and FOXM1 as potential transcriptional regulators of critical splicing factors in brain and liver. Our study provides a comprehensive demonstration of a splicing-mediated link between development and cancer, and suggest anti-cancer targets including splicing events, and their upstream splicing and transcriptional regulators.

Sex Biases in Cancer and Autoimmune Disease Incidence Are Strongly Positively Correlated with Mitochondrial Gene Expression across Human Tissues.

Cancer occurs more frequently in men while autoimmune diseases (AIDs) occur more frequently in women. To explore whether these sex biases have a common basis, we collected 167 AID incidence studies from many countries for tissues that have both a cancer type and an AID that arise from that tissue. Analyzing a total of 182 country-specific, tissue-matched cancer-AID incidence rate sex bias data pairs, we find that, indeed, the sex biases observed in the incidence of AIDs and cancers that occur in the same tissue are positively correlated across human tissues. The common key factor whose levels across human tissues are most strongly associated with these incidence rate sex biases is the sex bias in the expression of the 37 genes encoded in the mitochondrial genome.

The Hitchhiker's Guide across a Waddington's landscape of melanoma.

Melanoma evolution may recapitulate the embryonic development of its progenitor tissue, neural crest (NC), but the exact process is unclear. In a recent issue of Nature, Karras et al. (2022) demonstrate that melanoma expansion mirrors the hierarchic process of NC differentiation, generating cell subpopulations, each with distinct function, including growth and metastasis.

Cross-species identification of cancer resistance-associated genes that may mediate human cancer risk.

Cancer is a predominant disease across animals. We applied a comparative genomics approach to systematically characterize genes whose conservation levels correlate positively (PC) or negatively (NC) with cancer resistance estimates across 193 vertebrates. Pathway analysis reveals that NC genes are enriched for metabolic functions and PC genes in cell cycle regulation, DNA repair, and immune response, pointing to their corresponding roles in mediating cancer risk. We find that PC genes are less tolerant to loss-of-function (LoF) mutations, are enriched in cancer driver genes, and are associated with germline mutations that increase human cancer risk. Their relevance to cancer risk is further supported via the analysis of mouse functional genomics and cancer mortality of zoo mammals' data. In sum, our study describes a cross-species genomic analysis pointing to candidate genes that may mediate human cancer risk.

Making a mouse out of a molehill: how precision modeling repurposes drugs for congenital giant nevi.

Patients with congenital giant nevi (CGN), which can compromise quality of life and progress to melanoma, have limited treatment options. Choi et al. have demonstrated that topical application of a proinflammatory hapten for alopecia treatment [squaric acid dibutylester (SADBE)] caused nevus regression and prevented melanoma in an Nras mouse CGN model. Their results demonstrate the promise of repurposing drugs through precision modeling.

A T cell resilience model associated with response to immunotherapy in multiple tumor types.

Despite breakthroughs in cancer immunotherapy, most tumor-reactive T cells cannot persist in solid tumors due to an immunosuppressive environment. We developed Tres (tumor-resilient T cell), a computational model utilizing single-cell transcriptomic data to identify signatures of T cells that are resilient to immunosuppressive signals, such as transforming growth factor-ß1, tumor necrosis factor-related apoptosis-inducing ligand and prostaglandin E2. Tres reliably predicts clinical responses to immunotherapy in melanoma, lung cancer, triple-negative breast cancer and B cell malignancies using bulk T cell transcriptomic data from pre-treatment tumors from patients who received immune-checkpoint inhibitors (n = 38), infusion products for chimeric antigen receptor T cell therapies (n = 34) and pre-manufacture samples for chimeric antigen receptor T cell or tumor-infiltrating lymphocyte therapies (n = 84). Further, Tres identified FIBP, whose functions are largely unknown, as the top negative marker of tumor-resilient T cells across many solid tumor types. FIBP knockouts in murine and human donor CD8+ T cells significantly enhanced T cell-mediated cancer killing in in vitro co-cultures. Further, Fibp knockout in murine T cells potentiated the in vivo efficacy of adoptive cell transfer in the B16 tumor model. Fibp knockout T cells exhibit reduced cholesterol metabolism, which inhibits effector T cell function. These results demonstrate the utility of Tres in identifying biomarkers of T cell effectiveness and potential therapeutic targets for immunotherapies in solid tumors.

A novel PD-L1-targeted shark VNAR single-domain-based CAR-T cell strategy for treating breast cancer and liver cancer.

Chimeric antigen receptor (CAR)-T cell therapy shows excellent potency against hematological malignancies, but it remains challenging to treat solid tumors, mainly because of a lack of appropriate antigenic targets and an immunosuppressive tumor microenvironment (TME). The checkpoint molecule programmed death-ligand 1 (PD-L1) is widely overexpressed in multiple tumor types, and the programmed death-ligand 1 (PD-1)/PD-L1 interaction is a crucial mediator of immunosuppression in the TME. Here we constructed a semi-synthetic shark VNAR phage library and isolated anti-PD-L1 single-domain antibodies. Among these VNARs, B2 showed cross-reactivity to human, mouse, and canine PD-L1, and it partially blocked the interaction of human PD-1 with PD-L1. CAR (B2) T cells specifically lysed human breast cancer and liver cancer cells by targeting constitutive and inducible expression of PD-L1 and hindered tumor metastasis. Combination of PD-L1 CAR (B2) T cells with CAR T cells targeted by GPC3 (a liver cancer-specific antigen) regresses liver tumors in mice. We concluded that PD-L1-targeted shark VNAR single-domain-based CAR-T therapy is a novel strategy to treat breast and liver cancer. This study provides a rationale for potential use of PD-L1 CAR-T cells as a monotherapy or in combination with a tumor-specific therapy in clinical studies.

Fast intratumor heterogeneity inference from single-cell sequencing data.

We introduce HUNTRESS, a computational method for mutational intratumor heterogeneity inference from noisy genotype matrices derived from single-cell sequencing data, the running time of which is linear with the number of cells and quadratic with the number of mutations. We prove that, under reasonable conditions, HUNTRESS computes the true progression history of a tumor with high probability. On simulated and real tumor sequencing data, HUNTRESS is demonstrated to be faster than available alternatives with comparable or better accuracy. Additionally, the progression histories of tumors inferred by HUNTRESS on real single-cell sequencing datasets agree with the best known evolution scenarios for the associated tumors.

Dietary fiber and probiotics influence the gut microbiome and melanoma immunotherapy response.

Gut bacteria modulate the response to immune checkpoint blockade (ICB) treatment in cancer, but the effect of diet and supplements on this interaction is not well studied. We assessed fecal microbiota profiles, dietary habits, and commercially available probiotic supplement use in melanoma patients and performed parallel preclinical studies. Higher dietary fiber was associated with significantly improved progression-free survival in 128 patients on ICB, with the most pronounced benefit observed in patients with sufficient dietary fiber intake and no probiotic use. Findings were recapitulated in preclinical models, which demonstrated impaired treatment response to anti­programmed cell death 1 (anti­PD-1)­based therapy in mice receiving a low-fiber diet or probiotics, with a lower frequency of interferon-γ­positive cytotoxic T cells in the tumor microenvironment. Together, these data have clinical implications for patients receiving ICB for cancer.

Myeloid-Derived Suppressive Cell Expansion Promotes Melanoma Growth and Autoimmunity by Inhibiting CD40/IL27 Regulation in Macrophages.

The relationship between cancer and autoimmunity is complex. However, the incidence of solid tumors such as melanoma has increased significantly among patients with previous or newly diagnosed systemic autoimmune disease (AID). At the same time, immune checkpoint blockade (ICB) therapy of cancer induces de novo autoinflammation and exacerbates underlying AID, even without evident antitumor responses. Recently, systemic lupus erythematosus (SLE) activity was found to drive myeloid-derived suppressor cell (MDSC) formation in patients, a known barrier to healthy immune surveillance and successful cancer immunotherapy. Cross-talk between MDSCs and macrophages generally drives immune suppressive activity in the tumor microenvironment. However, it remains unclear how peripheral pregenerated MDSC under chronic inflammatory conditions modulates global macrophage immune functions and the impact it could have on existing tumors and underlying lupus nephritis. Here we show that pathogenic expansion of SLE-generated MDSCs by melanoma drives global macrophage polarization and simultaneously impacts the severity of lupus nephritis and tumor progression in SLE-prone mice. Molecular and functional data showed that MDSCs interact with autoimmune macrophages and inhibit cell surface expression of CD40 and the production of IL27. Moreover, low CD40/IL27 signaling in tumors correlated with high tumor-associated macrophage infiltration and ICB therapy resistance both in murine and human melanoma exhibiting active IFNγ signatures. These results suggest that preventing global macrophage reprogramming induced by MDSC-mediated inhibition of CD40/IL27 signaling provides a precision melanoma immunotherapy strategy, supporting an original and advantageous approach to treat solid tumors within established autoimmune landscapes. SIGNIFICANCE: Myeloid-derived suppressor cells induce macrophage reprogramming by suppressing CD40/IL27 signaling to drive melanoma progression, simultaneously affecting underlying autoimmune disease and facilitating resistance to immunotherapy within preexisting autoimmune landscapes.