RACER-m leverages structural features for sparse T cell specificity prediction.

Reliable prediction of T cell specificity against antigenic signatures is a formidable task, complicated by the immense diversity of T cell receptor and antigen sequence space and the resulting limited availability of training sets for inferential models. Recent modeling efforts have demonstrated the advantage of incorporating structural information to overcome the need for extensive training sequence data, yet disentangling the heterogeneous TCR-antigen interface to accurately predict MHC-allele-restricted TCR-peptide interactions has remained challenging. Here, we present RACER-m, a coarse-grained structural model leveraging key biophysical information from the diversity of publicly available TCR-antigen crystal structures. Explicit inclusion of structural content substantially reduces the required number of training examples and maintains reliable predictions of TCR-recognition specificity and sensitivity across diverse biological contexts. Our model capably identifies biophysically meaningful point-mutant peptides that affect binding affinity, distinguishing its ability in predicting TCR specificity of point-mutants from alternative sequence-based methods. Its application is broadly applicable to studies involving both closely related and structurally diverse TCR-peptide pairs.

Protocol for inferring epithelial-to-mesenchymal transition trajectories from single-cell RNA sequencing data using R.

The epithelial-to-mesenchymal transition (EMT) provides crucial insights into the metastatic process and possesses prognostic value within the cancer context. Here, we present COMET, an R package for inferring EMT trajectories and inter-state transition rates from single-cell RNA sequencing data. We describe steps for finding the optimal number of EMT genes for a specific context, estimating EMT-related trajectories, optimal fitting of continuous-time Markov chain to inferred trajectories, and estimating inter-state transition rates.

Quantitative approaches for decoding the specificity of the human T cell repertoire.

T cell receptor (TCR)-peptide-major histocompatibility complex (pMHC) interactions play a vital role in initiating immune responses against pathogens, and the specificity of TCRpMHC interactions is crucial for developing optimized therapeutic strategies. The advent of high-throughput immunological and structural evaluation of TCR and pMHC has provided an abundance of data for computational approaches that aim to predict favorable TCR-pMHC interactions. Current models are constructed using information on protein sequence, structures, or a combination of both, and utilize a variety of statistical learning-based approaches for identifying the rules governing specificity. This review examines the current theoretical, computational, and deep learning approaches for identifying TCR-pMHC recognition pairs, placing emphasis on each method's mathematical approach, predictive performance, and limitations.

Optimal phenotypic adaptation in fluctuating environments.

Phenotypic adaptation is a universal feature of biological systems navigating highly variable environments. Recent empirical data support the role of memory-driven decision making in cellular systems navigating uncertain future nutrient landscapes, wherein a distinct growth phenotype emerges in fluctuating conditions. We develop a simple stochastic mathematical model to describe memory-driven cellular adaptation required for systems to optimally navigate such uncertainty. In this framework, adaptive populations traverse dynamic environments by inferring future variation from a memory of prior states, and memory capacity imposes a fundamental trade-off between the speed and accuracy of adaptation to new fluctuating environments. Our results suggest that the observed growth reductions that occur in fluctuating environments are a direct consequence of optimal decision making and result from bet hedging and occasional phenotypic-environmental mismatch. We anticipate that this modeling framework will be useful for studying the role of memory in phenotypic adaptation, including in the design of temporally varying therapies against adaptive systems.

Population dynamics of EMT elucidates the timing and distribution of phenotypic intra-tumoral heterogeneity.

The Epithelial-to-Mesenchymal Transition (EMT) is a hallmark of cancer metastasis and morbidity. EMT is a non-binary process, and cells can be stably arrested en route to EMT in an intermediate hybrid state associated with enhanced tumor aggressiveness and worse patient outcomes. Understanding EMT progression in detail will provide fundamental insights into the mechanisms underlying metastasis. Despite increasingly available single-cell RNA sequencing (scRNA-seq) data that enable in-depth analyses of EMT at the single-cell resolution, current inferential approaches are limited to bulk microarray data. There is thus a great need for computational frameworks to systematically infer and predict the timing and distribution of EMT-related states at single-cell resolution. Here, we develop a computational framework for reliable inference and prediction of EMT-related trajectories from scRNA-seq data. Our model can be utilized across a variety of applications to predict the timing and distribution of EMT from single-cell sequencing data.

Optimal cancer evasion in a dynamic immune microenvironment generates diverse post-escape tumor antigenicity profiles.

The failure of cancer treatments, including immunotherapy, continues to be a major obstacle in preventing durable remission. This failure often results from tumor evolution, both genotypic and phenotypic, away from sensitive cell states. Here, we propose a mathematical framework for studying the dynamics of adaptive immune evasion that tracks the number of tumor-associated antigens available for immune targeting. We solve for the unique optimal cancer evasion strategy using stochastic dynamic programming and demonstrate that this policy results in increased cancer evasion rates compared to a passive, fixed strategy. Our foundational model relates the likelihood and temporal dynamics of cancer evasion to features of the immune microenvironment, where tumor immunogenicity reflects a balance between cancer adaptation and host recognition. In contrast with a passive strategy, optimally adaptive evaders navigating varying selective environments result in substantially heterogeneous post-escape tumor antigenicity, giving rise to immunogenically hot and cold tumors.

Epigenetic memory acquired during long-term EMT induction governs the recovery to the epithelial state.

Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFß need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFß withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.

Contact map dependence of a T-cell receptor binding repertoire.

The T-cell arm of the adaptive immune system provides the host protection against unknown pathogens by discriminating between host and foreign material. This discriminatory capability is achieved by the creation of a repertoire of cells each carrying a T-cell receptor (TCR) specific to non-self-antigens displayed as peptides bound to the major histocompatibility complex (pMHC). The understanding of the dynamics of the adaptive immune system at a repertoire level is complex, due to both the nuanced interaction of a TCR-pMHC pair and to the number of different possible TCR-pMHC pairings, making computationally exact solutions currently unfeasible. To gain some insight into this problem, we study an affinity-based model for TCR-pMHC binding in which a crystal structure is used to generate a distance-based contact map that weights the pairwise amino acid interactions. We find that the TCR-pMHC binding energy distribution strongly depends both on the number of contacts and the repeat structure allowed by the topology of the contact map of choice; this in turn influences T-cell recognition probability during negative selection, with higher variances leading to higher survival probabilities. In addition, we quantify the degree to which neoantigens with mutations in sites with higher contacts are recognized at a higher rate.

Changes in Triple-Negative Breast Cancer Molecular Subtypes in Patients Without Pathologic Complete Response After Neoadjuvant Systemic Chemotherapy.

PURPOSE: Lehmann et al have identified four molecular subtypes of triple-negative breast cancer (TNBC)-basal-like (BL) 1, BL2, mesenchymal (M), and luminal androgen receptor-and an immunomodulatory (IM) gene expression signature modifier. Our group previously showed that the response of TNBC to neoadjuvant systemic chemotherapy (NST) differs by molecular subtype, but whether NST affects the subtype was unknown. Here, we tested the hypothesis that in patients without pathologic complete response, TNBC subtypes can change after NST. Moreover, in cases with the changed subtype, we determined whether epithelial-to-mesenchymal transition (EMT) had occurred. MATERIALS AND METHODS: From the Pan-Pacific TNBC Consortium data set containing TNBC patient samples from four countries, we examined 64 formalin-fixed, paraffin-embedded pairs of matched pre- and post-NST tumor samples. The TNBC subtype was determined using the TNBCtype-IM assay. We analyzed a partial EMT gene expression scoring metric using mRNA data. RESULTS: Of the 64 matched pairs, 36 (56%) showed a change in the TNBC subtype after NST. The most frequent change was from BL1 to M subtypes (38%). No tumors changed from M to BL1. The IM signature was positive in 14 (22%) patients before NST and eight (12.5%) patients after NST. The EMT score increased after NST in 28 (78%) of the 36 patients with the changed subtype (v 39% of the 28 patients without change; P = .002254). CONCLUSION: We report, to our knowledge, for the first time that the TNBC molecular subtype and IM signature frequently change after NST. Our results also suggest that EMT is promoted by NST. Our findings may lead to innovative adjuvant therapy strategies in TNBC cases with residual tumor after NST.

Gene expression profiles of inflammatory breast cancer reveal high heterogeneity across the epithelial-hybrid-mesenchymal spectrum.

Inflammatory breast cancer (IBC) is a highly aggressive breast cancer that metastasizes largely via tumor emboli, and has a 5-year survival rate of less than 30%. No unique genomic signature has yet been identified for IBC nor has any specific molecular therapeutic been developed to manage the disease. Thus, identifying gene expression signatures specific to IBC remains crucial. Here, we compare various gene lists that have been proposed as molecular footprints of IBC using different clinical samples as training and validation sets and using independent training algorithms, and determine their accuracy in identifying IBC samples in three independent datasets. We show that these gene lists have little to no mutual overlap, and have limited predictive accuracy in identifying IBC samples. Despite this inconsistency, single-sample gene set enrichment analysis (ssGSEA) of IBC samples correlate with their position on the epithelial-hybrid-mesenchymal spectrum. This positioning, together with ssGSEA scores, improves the accuracy of IBC identification across the three independent datasets. Finally, we observed that IBC samples robustly displayed a higher coefficient of variation in terms of EMT scores, as compared to non-IBC samples. Pending verification that this patient-to-patient variability extends to intratumor heterogeneity within a single patient, these results suggest that higher heterogeneity along the epithelial-hybrid-mesenchymal spectrum can be regarded to be a hallmark of IBC and a possibly useful biomarker.

Implications of Tumor-Immune Coevolution on Cancer Evasion and Optimized Immunotherapy.

Cancer represents a diverse collection of diseases characterized by heterogeneous cell populations that dynamically evolve in their environment. As painfully evident in cases of treatment failure and recurrence, this general feature makes identifying long-term successful therapies difficult. It is now well-established that the adaptive immune system recognizes and eliminates cancer cells, and various immunotherapeutic strategies have emerged to augment this effect. These therapies, while promising, often fail as a result of immune-specific cancer evasion. Increasingly available empirical evidence details both cancer and immune system populations pre- and post-treatment, providing rich opportunity for mathematical models of the tumor-immune interaction and subsequent co-evolution. Integrated mathematical and experimental efforts bear immediate relevance for optimized therapies and will undoubtedly accelerate our understanding of this emergent field.

Transcriptomic-Based Quantification of the Epithelial-Hybrid-Mesenchymal Spectrum across Biological Contexts.

Epithelial-mesenchymal plasticity (EMP) underlies embryonic development, wound healing, and cancer metastasis and fibrosis. Cancer cells exhibiting EMP often have more aggressive behavior, characterized by drug resistance, and tumor-initiating and immuno-evasive traits. Thus, the EMP status of cancer cells can be a critical indicator of patient prognosis. Here, we compare three distinct transcriptomic-based metrics-each derived using a different gene list and algorithm-that quantify the EMP spectrum. Our results for over 80 cancer-related RNA-seq datasets reveal a high degree of concordance among these metrics in quantifying the extent of EMP. Moreover, each metric, despite being trained on cancer expression profiles, recapitulates the expected changes in EMP scores for non-cancer contexts such as lung fibrosis and cellular reprogramming into induced pluripotent stem cells. Thus, we offer a scoring platform to quantify the extent of EMP in vitro and in vivo for diverse biological applications including cancer.

Rapid Assessment of T-Cell Receptor Specificity of the Immune Repertoire.

Accurate assessment of TCR-antigen specificity at the whole immune repertoire level lies at the heart of improved cancer immunotherapy, but predictive models capable of high-throughput assessment of TCR-peptide pairs are lacking. Recent advances in deep sequencing and crystallography have enriched the data available for studying TCR-p-MHC systems. Here, we introduce a pairwise energy model, RACER, for rapid assessment of TCR-peptide affinity at the immune repertoire level. RACER applies supervised machine learning to efficiently and accurately resolve strong TCR-peptide binding pairs from weak ones. The trained parameters further enable a physical interpretation of interacting patterns encoded in each specific TCR-p-MHC system. When applied to simulate thymic selection of an MHC-restricted T-cell repertoire, RACER accurately estimates recognition rates for tumor-associated neoantigens and foreign peptides, thus demonstrating its utility in helping address the large computational challenge of reliably identifying the properties of tumor antigen-specific T-cells at the level of an individual patient's immune repertoire.

Targeting the Id1-Kif11 Axis in Triple-Negative Breast Cancer Using Combination Therapy.

The basic helix-loop-helix (bHLH) transcription factors inhibitor of differentiation 1 (Id1) and inhibitor of differentiation 3 (Id3) (referred to as Id) have an important role in maintaining the cancer stem cell (CSC) phenotype in the triple-negative breast cancer (TNBC) subtype. In this study, we aimed to understand the molecular mechanism underlying Id control of CSC phenotype and exploit it for therapeutic purposes. We used two different TNBC tumor models marked by either Id depletion or Id1 expression in order to identify Id targets using a combinatorial analysis of RNA sequencing and microarray data. Phenotypically, Id protein depletion leads to cell cycle arrest in the G0/G1 phase, which we demonstrate is reversible. In order to understand the molecular underpinning of Id proteins on the cell cycle phenotype, we carried out a large-scale small interfering RNA (siRNA) screen of 61 putative targets identified by using genomic analysis of two Id TNBC tumor models. Kinesin Family Member 11 (Kif11) and Aurora Kinase A (Aurka), which are critical cell cycle regulators, were further validated as Id targets. Interestingly, unlike in Id depletion conditions, Kif11 and Aurka knockdown leads to a G2/M arrest, suggesting a novel Id cell cycle mechanism, which we will explore in further studies. Therapeutic targeting of Kif11 to block the Id1-Kif11 axis was carried out using small molecular inhibitor ispinesib. We finally leveraged our findings to target the Id/Kif11 pathway using the small molecule inhibitor ispinesib in the Id+ CSC results combined with chemotherapy for better response in TNBC subtypes. This work opens up exciting new possibilities of targeting Id targets such as Kif11 in the TNBC subtype, which is currently refractory to chemotherapy. Targeting the Id1-Kif11 molecular pathway in the Id1+ CSCs in combination with chemotherapy and small molecular inhibitor results in more effective debulking of TNBC.

Comparative Study of Transcriptomics-Based Scoring Metrics for the Epithelial-Hybrid-Mesenchymal Spectrum.

The Epithelial-mesenchymal transition (EMT) is a cellular process implicated in embryonic development, wound healing, and pathological conditions such as cancer metastasis and fibrosis. Cancer cells undergoing EMT exhibit enhanced aggressive behavior characterized by drug resistance, tumor-initiation potential, and the ability to evade the immune system. Recent in silico, in vitro, and in vivo evidence indicates that EMT is not an all-or-none process; instead, cells can stably acquire one or more hybrid epithelial/mesenchymal (E/M) phenotypes which often can be more aggressive than purely E or M cell populations. Thus, the EMT status of cancer cells can prove to be a critical estimate of patient prognosis. Recent attempts have employed different transcriptomics signatures to quantify EMT status in cell lines and patient tumors. However, a comprehensive comparison of these methods, including their accuracy in identifying cells in the hybrid E/M phenotype(s), is lacking. Here, we compare three distinct metrics that score EMT on a continuum, based on the transcriptomics signature of individual samples. Our results demonstrate that these methods exhibit good concordance among themselves in quantifying the extent of EMT in a given sample. Moreover, scoring EMT using any of the three methods discerned that cells can undergo varying extents of EMT across tumor types. Separately, our analysis also identified tumor types with maximum variability in terms of EMT and associated an enrichment of hybrid E/M signatures in these samples. Moreover, we also found that the multinomial logistic regression (MLR)-based metric was capable of distinguishing between "pure" individual hybrid E/M vs. mixtures of E and M cells. Our results, thus, suggest that while any of the three methods can indicate a generic trend in the EMT status of a given cell, the MLR method has two additional advantages: (a) it uses a small number of predictors to calculate the EMT score and (b) it can predict from the transcriptomic signature of a population whether it is comprised of "pure" hybrid E/M cells at the single-cell level or is instead an ensemble of E and M cell subpopulations.

Histone deacetylases, Mbd3/NuRD, and Tet2 hydroxylase are crucial regulators of epithelial-mesenchymal plasticity and tumor metastasis.

An epithelial-mesenchymal transition (EMT) represents a basic morphogenetic process of high cell plasticity underlying embryogenesis, wound healing, cancer metastasis and drug resistance. It involves a profound transcriptional and epigenetic reprogramming of cells. A critical role of epigenetic modifiers and their specific chromatin modifications has been demonstrated during EMT. However, it has remained elusive whether epigenetic control differs between the dynamic cell state transitions of reversible EMT and the fixed differentiation status of irreversible EMT. We have employed varying EMT models of murine breast cancer cells to identify the key players establishing epithelial-mesenchymal cell plasticity during reversible and irreversible EMT. We demonstrate that the Mbd3/NuRD complex and the activities of histone deacetylases (HDACs), and Tet2 hydroxylase play a critical role in keeping cancer cells in a highly metastatic mesenchymal state. Combinatorial interference with their functions leads to mesenchymal-epithelial transition (MET) and efficiently represses metastasis formation by invasive murine and human breast cancer cells in vivo.

Sustained Coevolution in a Stochastic Model of Cancer-Immune Interaction.

The dynamic interactions between an evolving malignancy and the adaptive immune system generate diverse evolutionary trajectories that ultimately result in tumor clearance or immune escape. Here, we create a simple mathematical model coupling T-cell recognition with an evolving cancer population that may randomly produce evasive subclones, imparting transient protection against the effector T cells. T-cell turnover declines and evasion rates together explained differences in early incidence data across almost all cancer types. Fitting the model to TRACERx evolutionary data argued in favor of substantial and sustained immune pressure exerted upon a developing tumor, suggesting that clinically observed incidence is a small proportion of all cancer initiation events. This dynamical model promises to increase our quantitative understanding of many immune escape contexts, including cancer progression and intracellular pathogenic infections. SIGNIFICANCE: The early cancer-immune interaction sculpts intratumor heterogeneity through the selection of immune-evasive clones. This study provides a mathematical framework for investigating the coevolution between an immune-evasive cancer population and the adaptive immune system.

Differential Contributions of Pre- and Post-EMT Tumor Cells in Breast Cancer Metastasis.

Metastases are responsible for the majority of breast cancer-associated deaths. The contribution of epithelial-to-mesenchymal transition (EMT) in the establishment of metastases is still controversial. To obtain in vivo evidence of EMT in metastasis, we established an EMT lineage tracing (Tri-PyMT) model, in which tumor cells undergoing EMT would irreversibly switch their fluorescent marker from RFP+ to GFP+ due to mesenchymal-specific Cre expression. Surprisingly, we found that lung metastases were predominantly derived from the epithelial compartment of breast tumors. However, concerns were raised on the fidelity and sensitivity of RFP-to-GFP switch of this model in reporting EMT of metastatic tumor cells. Here, we evaluated Tri-PyMT cells at the single-cell level using single-cell RNA-sequencing and found that the Tri-PyMT cells exhibited a spectrum of EMT phenotypes, with EMT-related genes concomitantly expressed with the activation of GFP. The fluorescent color switch in these cells precisely marked an unequivocal change in EMT status, defining the pre-EMT and post-EMT compartments within the tumor. Consistently, the pre-EMT cells played dominant roles in metastasis, while the post-EMT cells were supportive in promoting tumor invasion and angiogenesis. Importantly, the post-EMT (GFP+) cells in the Tri-PyMT model were not permanently committed to the mesenchymal phenotype; they were still capable of reverting to the epithelial phenotype and giving rise to secondary tumors, suggesting their persistent EMT plasticity. Our study addressed major concerns with the Tri-PyMT EMT lineage tracing model, which provides us with a powerful tool to investigate the dynamic EMT process in tumor biology. SIGNIFICANCE: These findings confirm the fidelity and sensitivity of the EMT lineage tracing (Tri-PyMT) model and highlight the differential contributions of pre- and post-EMT tumor cells in breast cancer metastasis.See related commentary by Bunz, p. 153.

A CTC-Cluster-Specific Signature Derived from OMICS Analysis of Patient-Derived Xenograft Tumors Predicts Outcomes in Basal-Like Breast Cancer.

Circulating tumor cell clusters (CTCcl) have a higher metastatic potential compared to single CTCs and predict long-term outcomes in breast cancer (BC) patients. Because of the rarity of CTCcls, molecular characterization of primary tumors that give rise to CTCcl hold significant promise for better diagnosis and target discovery to combat metastatic BC. In our study, we utilized the reverse-phase protein array (RPPA) and transcriptomic (RNA-Seq) data of 10 triple-negative BC patient-derived xenograft (TNBC PDX) transplantable models with CTCs and evaluated expression of upregulated candidate protein Bcl2 (B-cell lymphoma 2) by immunohistochemistry (IHC). The sample-set consisted of six CTCcl-negative (CTCcl-) and four CTCcl-positive (CTCcl+) models. We analyzed the RPPA and transcriptomic profiles of CTCcl- and CTCcl+ TNBC PDX models. In addition, we derived a CTCcl-specific gene signature for testing if it predicted outcomes using a publicly available dataset from 360 patients with basal-like BC. The RPPA analysis of CTCcl+ vs. CTCcl- TNBC PDX tumors revealed elevated expression of Bcl2 (false discovery rate (FDR) < 0.0001, fold change (FC) = 3.5) and reduced acetyl coenzyme A carboxylase-1 (ACC1) (FDR = 0.0005, FC = 0.3) in CTCcl+ compared to CTCcl- tumors. Genome-wide transcriptomic analysis of CTCcl+ vs. CTCcl- tumors revealed 549 differentially expressed genes associated with the presence of CTCcls. Apoptosis was one of the significantly downregulated pathways (normalized enrichment score (NES) = -1.69; FDR < 0.05) in TNBC PDX tumors associated with CTCcl positivity. Two out of four CTCcl+ TNBC PDX primary tumors had high Bcl2 expression by IHC (H-score > 34); whereas, only one of six CTCcl- TNBC PDX primary tumors met this criterion. Evaluation of epithelial-mesenchymal transition (EMT)-specific signature did not show significant differences between CTCcl+ and CTCcl- tumors. However, a gene signature associated with the presence of CTCcls in TNBC PDX models was associated with worse relapse-free survival in the publicly available dataset from 360 patients with basal-like BC. In summary, we identified the multigene signature of primary PDX tumors associated with the presence of CTCcls. Evaluation of additional TNBC PDX models and patients can further illuminate cellular and molecular pathways facilitating CTCcl formation.