TimiGP-Response: the pan-cancer immune landscape associated with response to immunotherapy.

Accumulating evidence suggests that the tumor immune microenvironment (TIME) significantly influences the response to immunotherapy, yet this complex relationship remains elusive. To address this issue, we developed TimiGP-Response (TIME Illustration based on Gene Pairing designed for immunotherapy Response), a computational framework leveraging single-cell and bulk transcriptomic data, along with response information, to construct cell-cell interaction networks associated with responders and estimate the role of immune cells in treatment response. This framework was showcased in triple-negative breast cancer treated with immune checkpoint inhibitors targeting the PD-1:PD-L1 interaction, and orthogonally validated with imaging mass cytometry. As a result, we identified CD8+ GZMB+ T cells associated with responders and its interaction with regulatory T cells emerged as a potential feature for selecting patients who may benefit from these therapies. Subsequently, we analyzed 3,410 patients with seven cancer types (melanoma, non-small cell lung cancer, renal cell carcinoma, metastatic urothelial carcinoma, hepatocellular carcinoma, breast cancer, and esophageal cancer) treated with various immunotherapies and combination therapies, as well as several chemo- and targeted therapies as controls. Using TimiGP-Response, we depicted the pan-cancer immune landscape associated with immunotherapy response at different resolutions. At the TIME level, CD8 T cells and CD4 memory T cells were associated with responders, while anti-inflammatory (M2) macrophages and mast cells were linked to non-responders across most cancer types and datasets. Given that T cells are the primary targets of these immunotherapies and our TIME analysis highlights their importance in response to treatment, we portrayed the pan-caner landscape on 40 T cell subtypes. Notably, CD8+ and CD4+ GZMK+ effector memory T cells emerged as crucial across all cancer types and treatments, while IL-17-producing CD8+ T cells were top candidates associated with immunotherapy non-responders. In summary, this study provides a computational method to study the association between TIME and response across the pan-cancer immune landscape, offering resources and insights into immune cell interactions and their impact on treatment efficacy.

Integrative genomic and transcriptomic profiling of pulmonary sarcomatoid carcinoma identifies molecular subtypes associated with distinct immune features and clinical outcomes.

Background: Pulmonary sarcomatoid carcinoma (PSC) is a rare and aggressive subtype of non-small cell lung cancer (NSCLC), characterized by the presence of epithelial and sarcoma-like components. The molecular and immune landscape of PSC has not been well defined. Methods: Multiomics profiling of 21 pairs of PSCs with matched normal lung tissues was performed through targeted high-depth DNA panel, whole-exome, and RNA sequencing. We describe molecular and immune features that define subgroups of PSC with disparate genomic and immunogenic features as well as distinct clinical outcomes. Results: In total, 27 canonical cancer gene mutations were identified, with TP53 the most frequently mutated gene, followed by KRAS. Interestingly, most TP53 and KRAS mutations were earlier genomic events mapped to the trunks of the tumors, suggesting branching evolution in most PSC tumors. We identified two distinct molecular subtypes of PSC, driven primarily by immune infiltration and signaling. The Immune High (IM-H) subtype was associated with superior survival, highlighting the impact of immune infiltration on the biological and clinical features of localized PSCs. Conclusions: We provided detailed insight into the mutational landscape of PSC and identified two molecular subtypes associated with prognosis. IM-H tumors were associated with favorable recurrence-free survival and overall survival, highlighting the importance of tumor immune infiltration in the biological and clinical features of PSCs.

The combination of tumor mutational burden and T-cell receptor repertoire predicts the response to immunotherapy in patients with advanced non-small cell lung cancer.

Tumor mutational burden (TMB) and T-cell receptor (TCR) might predict the response to immunotherapy in patients with non-small cell lung cancer (NSCLC). However, the predictive value of the combination of TMB and TCR was not clear. Targeted DNA and TCR sequencing were performed on tumor biopsy specimens. We combined TMB and TCR diversity into a TMB-and-TCR (TMR) score using logistic regression. In total, 38 patients with advanced NSCLC were divided into a discovery set (n = 17) and validation set (n = 21). A higher TMR score was associated with better response and longer progression-free survival to immunotherapy in both the discovery set and validation set. The performance of TMR score was confirmed in the two external validation cohorts of 225 NSCLC patients and 306 NSCLC patients. Tumors with higher TMR scores were more likely to combine with LRP1B gene mutation (p = 0.027) and top 1% CDR3 sequences (p = 0.001). Furthermore, LRP1B allele frequency was negatively correlated with the top 1% CDR3 sequences (r = -0.55, p = 0.033) and positively correlated with tumor shrinkage (r = 0.68, p = 0.007). The TMR score could serve as a potential predictive biomarker for the response to immunotherapy in advanced NSCLC.

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion.

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

Maintenance Pembrolizumab Therapy in Patients with Metastatic HER2-negative Breast Cancer with Prior Response to Chemotherapy.

PURPOSE: Accumulating toxicities hinder indefinite chemotherapy for many patients with metastatic/recurrent HER2-negative breast cancer. We conducted a phase II trial of pembrolizumab monotherapy following induction chemotherapy to determine the efficacy of maintenance immunotherapy in patients with metastatic HER2-negative inflammatory breast cancer (IBC) and non-IBC triple-negative breast cancer (TNBC) and a biomarker study. PATIENTS AND METHODS: Patients with a complete response, partial response, or stable disease (SD) after at least three cycles of chemotherapy for HER2-negative breast cancer received pembrolizumab, regardless of programmed death-ligand 1 expression. Pembrolizumab (200 mg) was administered every 3 weeks until disease progression, intolerable toxicity, or 2 years of pembrolizumab exposure. The endpoints included the 4-month disease control rate (DCR), progression-free survival (PFS), overall survival, and response biomarkers in the blood. RESULTS: Of 43 treated patients, 11 had metastatic IBC and 32 non-IBC TNBC. The 4-month DCR was 58.1% [95% confidence interval (CI), 43.4-72.9]. For all patients, the median PFS was 4.8 months (95% CI, 3.0-7.1 months). The toxicity profile was similar to the previous pembrolizumab monotherapy study. Patients with high T-cell clonality at baseline had a longer PFS with pembrolizumab treatment than did those with low T-cell clonality (10.4 vs. 3.6 months, P = 0.04). Patients who achieved SD also demonstrated a significant increase in T-cell clonality during therapy compared with those who did not achieve SD (20% vs. 5.9% mean increase, respectively; P = 0.04). CONCLUSIONS: Pembrolizumab monotherapy achieved durable treatment responses. Patients with a high baseline T-cell clonality had prolonged disease control with pembrolizumab.

Assessment of human leukocyte antigen-based neoantigen presentation to determine pan-cancer response to immunotherapy.

Despite the central role of human leukocyte antigen class I (HLA-I) in tumor neoantigen presentation, quantitative determination of presentation capacity remains elusive. Based on a pooled pan-cancer genomic dataset of 885 patients treated with immune checkpoint inhibitors (ICIs), we developed a score integrating the binding affinity of neoantigens to HLA-I, as well as HLA-I allele divergence, termed the HLA tumor-Antigen Presentation Score (HAPS). Patients with a high HAPS were more likely to experience survival benefit following ICI treatment. Analysis of the tumor microenvironment indicated that the antigen presentation pathway was enriched in patients with a high HAPS. Finally, we built a neural network incorporating factors associated with neoantigen production, presentation, and recognition, which exhibited potential for differentiating cancer patients likely to benefit from ICIs. Our findings highlight the clinical utility of evaluating HLA-I tumor antigen presentation capacity and describe how ICI response may depend on HLA-mediated immunity.

Genomic Landscape of Lynch Syndrome Colorectal Neoplasia Identifies Shared Mutated Neoantigens for Immunoprevention.

BACKGROUND & AIMS: Lynch syndrome (LS) carriers develop mismatch repair-deficient neoplasia with high neoantigen (neoAg) rates. No detailed information on targetable neoAgs from LS precancers exists, which is crucial for vaccine development and immune-interception strategies. We report a focused somatic mutation and frameshift-neoAg landscape of microsatellite loci from colorectal polyps without malignant potential (PWOMP), precancers, and early-stage cancers in LS carriers. METHODS: We generated paired whole-exome and transcriptomic sequencing data from 8 colorectal PWOMP, 41 precancers, 8 advanced precancers, and 12 early-stage cancers of 43 LS carriers. A computational pipeline was developed to predict, rank, and prioritize the top 100 detected mutated neoAgs that were validated in vitro using ELISpot and tetramer assays. RESULTS: Mutation calling revealed >10 mut/Mb in 83% of cancers, 63% of advanced precancers, and 20% of precancers. Cancers displayed an average of 616 MHC-I neoAgs/sample, 294 in advanced precancers, and 107 in precancers. No neoAgs were detected in PWOMP. A total of 65% of our top 100 predicted neoAgs were immunogenic in vitro, and were present in 92% of cancers, 50% of advanced precancers, and 29% of precancers. We observed increased levels of naïve CD8+ and memory CD4+ T cells in mismatch repair-deficient cancers and precancers via transcriptomics analysis. CONCLUSIONS: Shared frameshift-neoAgs are generated within unstable microsatellite loci at initial stages of LS carcinogenesis and can induce T-cell responses, generating opportunities for vaccine development, targeting LS precancers and early-stage cancers.

pan-MHC and cross-Species Prediction of T Cell Receptor-Antigen Binding.

Profiling the binding of T cell receptors (TCRs) of T cells to antigenic peptides presented by MHC proteins is one of the most important unsolved problems in modern immunology. Experimental methods to probe TCR-antigen interactions are slow, labor-intensive, costly, and yield moderate throughput. To address this problem, we developed pMTnet-omni, an Artificial Intelligence (AI) system based on hybrid protein sequence and structure information, to predict the pairing of TCRs of αß T cells with peptide-MHC complexes (pMHCs). pMTnet-omni is capable of handling peptides presented by both class I and II pMHCs, and capable of handling both human and mouse TCR-pMHC pairs, through information sharing enabled this hybrid design. pMTnet-omni achieves a high overall Area Under the Curve of Receiver Operator Characteristics (AUROC) of 0.888, which surpasses competing tools by a large margin. We showed that pMTnet-omni can distinguish binding affinity of TCRs with similar sequences. Across a range of datasets from various biological contexts, pMTnet-omni characterized the longitudinal evolution and spatial heterogeneity of TCR-pMHC interactions and their functional impact. We successfully developed a biomarker based on pMTnet-omni for predicting immune-related adverse events of immune checkpoint inhibitor (ICI) treatment in a cohort of 57 ICI-treated patients. pMTnet-omni represents a major advance towards developing a clinically usable AI system for TCR-pMHC pairing prediction that can aid the design and implementation of TCR-based immunotherapeutics.

Spatial heterogeneity of T cell repertoire across NSCLC tumors, tumor edges, adjacent and distant lung tissues.

Background: A better understanding of T cells in lung cancer and their distribution across tumor-adjacent lungs and peripheral blood is needed to improve efficacy and minimize toxicity from immunotherapy to lung cancer patients. Methods: Here, we performed CDR3ß TCR sequencing of 136 samples from 20 patients with early-stage NSCLC including peripheral blood mononuclear cells, tumors, tumor edges (<1 cm from tumor), as well as adjacent lungs 1 cm, 2 cm, 5 cm, and 10 cm away from the tumor to gain insight into the spatial heterogeneity of T cells across the lungs in patients with NSCLC. PD-L1, CD4, and CD8 expression was assessed using immunohistochemical staining, and genomic features were derived by targeted sequencing of 1,021 cancer-related genes. Multiplex immunohistochemistry against PD-1, CTLA4, LAG3, and TIM3 was performed on four patients to assess T cell exhaustion. Results: Our study reveals a decreasing gradient in TIL Tumor Infiltrating Lymphocytes homology with tumor edge, adjacent lungs, and peripheral blood but no discernible distance-associated patterns of T cell trafficking within the adjacent lung itself. Furthermore, we show a decrease in pathogen-specific TCRs in regions with high T cell clonality and PD-L1 expression. Conclusions: Exclusion in T exhaustion cells at play across the lungs of patients with NSCLC may potentially be the mechanism for lung cancer occurrence.

Immune-active tumor-adjacent tissues are associated with favorable prognosis in stage I lung squamous cell carcinoma.

The immunogenomic features of tumor-adjacent lungs (TALs) in stage I lung squamous cell carcinoma (LUSC) are not clear. Multiomics analyses of tumor tissues and paired TALs from 59 stage I LUSC patients were performed. Compared to tumors, TALs exhibited a better-preserved immune contexture indicated by upregulation of immune pathways, increased immune infiltration, and higher expression of immune effector molecules. Notably, TALs had no mutations in PTEN and KEAP1, a lower incidence of human leukocyte antigen (HLA) loss and higher expression of HLA class I genes, major histocompatibility complex (MHC) I chaperones, and interferon (IFN)-γ-associated genes. Digital spatial profiling validated the generally higher immune infiltration in TALs and revealed a higher level of immune heterogeneity in LUSC tumors. Importantly, patients with higher immune infiltration in TALs had significantly longer survival, while high immune heterogeneity was associated with inferior patient survival. Our work can be considered in the selection of patients for adjuvant therapy, especially immunotherapy.

TimiGP: Inferring cell-cell interactions and prognostic associations in the tumor immune microenvironment through gene pairs.

Determining the prognostic association of different immune cell types in the tumor microenvironment is critical for understanding cancer biology and developing new therapeutic strategies. However, this is challenging in certain cancer types, where the abundance of different immune subsets is highly correlated. In this study, we develop a computational method named TimiGP to overcome this challenge. Based on bulk gene expression and survival data, TimiGP infers cell-cell interactions that reveal the association between immune cell relative abundance and prognosis. As demonstrated in metastatic melanoma, TimiGP prioritizes immune cells critical in prognosis based on the identified cell-cell interactions. Highly consistent results are obtained by TimiGP when applied to seven independent melanoma datasets and when different cell-type marker sets are used as inputs. Additionally, TimiGP can leverage single-cell RNA sequencing data to delineate the tumor immune microenvironment at high resolutions across a wide range of cancer types.

A novel statistical method for decontaminating T-cell receptor sequencing data.

The T-cell receptor (TCR) repertoire is highly diverse among the population and plays an essential role in initiating multiple immune processes. TCR sequencing (TCR-seq) has been developed to profile the T cell repertoire. Similar to other high-throughput experiments, contamination can happen during several steps of TCR-seq, including sample collection, preparation and sequencing. Such contamination creates artifacts in the data, leading to inaccurate or even biased results. Most existing methods assume 'clean' TCR-seq data as the starting point with no ability to handle data contamination. Here, we develop a novel statistical model to systematically detect and remove contamination in TCR-seq data. We summarize the observed contamination into two sources, pairwise and cross-cohort. For both sources, we provide visualizations and summary statistics to help users assess the severity of the contamination. Incorporating prior information from 14 existing TCR-seq datasets with minimum contamination, we develop a straightforward Bayesian model to statistically identify contaminated samples. We further provide strategies for removing the impacted sequences to allow for downstream analysis, thus avoiding any need to repeat experiments. Our proposed model shows robustness in contamination detection compared with a few off-the-shelf detection methods in simulation studies. We illustrate the use of our proposed method on two TCR-seq datasets generated locally.

MCT4-dependent lactate secretion suppresses antitumor immunity in LKB1-deficient lung adenocarcinoma.

Inactivating STK11/LKB1 mutations are genomic drivers of primary resistance to immunotherapy in KRAS-mutated lung adenocarcinoma (LUAD), although the underlying mechanisms remain unelucidated. We find that LKB1 loss results in enhanced lactate production and secretion via the MCT4 transporter. Single-cell RNA profiling of murine models indicates that LKB1-deficient tumors have increased M2 macrophage polarization and hypofunctional T cells, effects that could be recapitulated by the addition of exogenous lactate and abrogated by MCT4 knockdown or therapeutic blockade of the lactate receptor GPR81 expressed on immune cells. Furthermore, MCT4 knockout reverses the resistance to PD-1 blockade induced by LKB1 loss in syngeneic murine models. Finally, tumors from STK11/LKB1 mutant LUAD patients demonstrate a similar phenotype of enhanced M2-macrophages polarization and hypofunctional T cells. These data provide evidence that lactate suppresses antitumor immunity and therapeutic targeting of this pathway is a promising strategy to reversing immunotherapy resistance in STK11/LKB1 mutant LUAD.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance.

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.

Addition of Metastasis-Directed Therapy to Intermittent Hormone Therapy for Oligometastatic Prostate Cancer: The EXTEND Phase 2 Randomized Clinical Trial.

Importance: Despite evidence demonstrating an overall survival benefit with up-front hormone therapy in addition to established synergy between hormone therapy and radiation, the addition of metastasis-directed therapy (MDT) to hormone therapy for oligometastatic prostate cancer, to date, has not been evaluated in a randomized clinical trial. Objective: To determine in men with oligometastatic prostate cancer whether the addition of MDT to intermittent hormone therapy improves oncologic outcomes and preserves time with eugonadal testosterone compared with intermittent hormone therapy alone. Design, Setting, Participants: The External Beam Radiation to Eliminate Nominal Metastatic Disease (EXTEND) trial is a phase 2, basket randomized clinical trial for multiple solid tumors testing the addition of MDT to standard-of-care systemic therapy. Men aged 18 years or older with oligometastatic prostate cancer who had 5 or fewer metastases and were treated with hormone therapy for 2 or more months were enrolled to the prostate intermittent hormone therapy basket at multicenter tertiary cancer centers from September 2018 to November 2020. The cutoff date for the primary analysis was January 7, 2022. Interventions: Patients were randomized 1:1 to MDT, consisting of definitive radiation therapy to all sites of disease and intermittent hormone therapy (combined therapy arm; n = 43) or to hormone therapy only (n = 44). A planned break in hormone therapy occurred 6 months after enrollment, after which hormone therapy was withheld until progression. Main Outcomes and Measures: The primary end point was disease progression, defined as death or radiographic, clinical, or biochemical progression. A key predefined secondary end point was eugonadal progression-free survival (PFS), defined as the time from achieving a eugonadal testosterone level (≥150 ng/dL; to convert to nanomoles per liter, multiply by 0.0347) until progression. Exploratory measures included quality of life and systemic immune evaluation using flow cytometry and T-cell receptor sequencing. Results: The study included 87 men (median age, 67 years [IQR, 63-72 years]). Median follow-up was 22.0 months (range, 11.6-39.2 months). Progression-free survival was improved in the combined therapy arm (median not reached) compared with the hormone therapy only arm (median, 15.8 months; 95% CI, 13.6-21.2 months) (hazard ratio, 0.25; 95% CI, 0.12-0.55; P < .001). Eugonadal PFS was also improved with MDT (median not reached) compared with the hormone therapy only (6.1 months; 95% CI, 3.7 months to not estimable) (hazard ratio, 0.32; 95% CI, 0.11-0.91; P = .03). Flow cytometry and T-cell receptor sequencing demonstrated increased markers of T-cell activation, proliferation, and clonal expansion limited to the combined therapy arm. Conclusions and Relevance: In this randomized clinical trial, PFS and eugonadal PFS were significantly improved with combination treatment compared with hormone treatment only in men with oligometastatic prostate cancer. Combination of MDT with intermittent hormone therapy may allow for excellent disease control while facilitating prolonged eugonadal testosterone intervals. Trial Registration: ClinicalTrials.gov Identifier: NCT03599765.

Analysis of Immune Intratumor Heterogeneity Highlights Immunoregulatory and Coinhibitory Lymphocytes as Hallmarks of Recurrence in Stage I Non-Small Cell Lung Cancer.

Our understanding of the molecular mechanisms underlying postsurgical recurrence of non-small cell lung cancer (NSCLC) is rudimentary. Molecular and T cell repertoire intratumor heterogeneity (ITH) have been reported to be associated with postsurgical relapse; however, how ITH at the cellular level impacts survival is largely unknown. Here we report the analysis of 2880 multispectral images representing 14.2% to 27% of tumor areas from 33 patients with stage I NSCLC, including 17 cases (relapsed within 3 years after surgery) and 16 controls (without recurrence ≥5 years after surgery) using multiplex immunofluorescence. Spatial analysis was conducted to quantify the minimum distance between different cell types and immune cell infiltration around malignant cells. Immune ITH was defined as the variance of immune cells from 3 intratumor regions. We found that tumors from patients having relapsed display different immune biology compared with nonrecurrent tumors, with a higher percentage of tumor cells and macrophages expressing PD-L1 (P =.031 and P =.024, respectively), along with an increase in regulatory T cells (Treg) (P =.018), antigen-experienced T cells (P =.025), and effector-memory T cells (P =.041). Spatial analysis revealed that a higher level of infiltration of PD-L1+ macrophages (CD68+PD-L1+) or antigen-experienced cytotoxic T cells (CD3+CD8+PD-1+) in the tumor was associated with poor overall survival (P =.021 and P =.006, respectively). A higher degree of Treg ITH was associated with inferior recurrence-free survival regardless of tumor mutational burden (P =.022), neoantigen burden (P =.021), genomic ITH (P =.012) and T cell repertoire ITH (P =.001). Using multiregion multiplex immunofluorescence, we characterized ITH at the immune cell level along with whole exome and T cell repertoire sequencing from the same tumor regions. This approach highlights the role of immunoregulatory and coinhibitory signals as well as their spatial distribution and ITH that define the hallmarks of tumor relapse of stage I NSCLC.

Analysis of genomic and immune intratumor heterogeneity in linitis plastica via multiregional exome and T-cell receptor sequencing.

The molecular landscape and the intratumor heterogeneity (ITH) architecture of gastric linitis plastica (LP) are poorly understood. We performed whole-exome sequencing (WES) and T-cell receptor (TCR) sequencing on 40 tumor regions from four LP patients. The landscape and ITH at the genomic and immunological levels in LP tumors were compared with multiple cancers that have previously been reported. The lymphocyte infiltration was further assessed by immunohistochemistry (IHC) in LP tumors. In total, we identified 6339 non-silent mutations from multi-samples, with a median tumor mutation burden (TMB) of 3.30 mutations per Mb, comparable to gastric adenocarcinoma from the Cancer Genome Atlas (TCGA) cohort (P = 0.53). An extremely high level of genomic ITH was observed, with only 12.42%, 5.37%, 5.35%, and 30.67% of mutations detectable across 10 regions within the same tumors of each patient, respectively. TCR sequencing revealed that TCR clonality was substantially lower in LP than in multi-cancers. IHC using antibodies against CD4, CD8, and PD-L1 demonstrated scant T-cell infiltration in the four LP tumors. Furthermore, profound TCR ITH was observed in all LP tumors, with no T-cell clones shared across tumor regions in any of the patients, while over 94% of T-cell clones were restricted to individual tumor regions. The Morisita overlap index (MOI) ranged from 0.21 to 0.66 among multi-regions within the same tumors, significantly lower than that of lung cancer (P = 0.002). Our results show that LP harbored extremely high genomic and TCR ITH and suppressed T-cell infiltration, suggesting a potential contribution to the frequent recurrence and poor therapeutic response of this adenocarcinoma.

T-Cell Receptor Repertoire Sequencing in the Era of Cancer Immunotherapy.

T cells are integral components of the adaptive immune system, and their responses are mediated by unique T-cell receptors (TCR) that recognize specific antigens from a variety of biological contexts. As a result, analyzing the T-cell repertoire offers a better understanding of immune responses and of diseases like cancer. Next-generation sequencing technologies have greatly enabled the high-throughput analysis of the TCR repertoire. On the basis of our extensive experience in the field from the past decade, we provide an overview of TCR sequencing, from the initial library preparation steps to sequencing and analysis methods and finally to functional validation techniques. With regards to data analysis, we detail important TCR repertoire metrics and present several computational tools for predicting antigen specificity. Finally, we highlight important applications of TCR sequencing and repertoire analysis to understanding tumor biology and developing cancer immunotherapies.

Bayesian hierarchical quantile regression with application to characterizing the immune architecture of lung cancer.

The successful development and implementation of precision immuno-oncology therapies requires a deeper understanding of the immune architecture at a patient level. T-cell receptor (TCR) repertoire sequencing is a relatively new technology that enables monitoring of T-cells, a subset of immune cells that play a central role in modulating immune response. These immunologic relationships are complex and are governed by various distributional aspects of an individual patient's tumor profile. We propose Bayesian QUANTIle regression for hierarchical COvariates (QUANTICO) that allows simultaneous modeling of hierarchical relationships between multilevel covariates, conducts explicit variable selection, estimates quantile and patient-specific coefficient effects, to induce individualized inference. We show QUANTICO outperforms existing approaches in multiple simulation scenarios. We demonstrate the utility of QUANTICO to investigate the effect of TCR variables on immune response in a cohort of lung cancer patients. At population level, our analyses reveal the mechanistic role of T-cell proportion on the immune cell abundance, with tumor mutation burden as an important factor modulating this relationship. At a patient level, we find several outlier patients based on their quantile-specific coefficient functions, who have higher mutational rates and different smoking history.

Netie: inferring the evolution of neoantigen-T cell interactions in tumors.

Neoantigens are the key targets of antitumor immune responses from cytotoxic T cells and play a critical role in affecting tumor progressions and immunotherapy treatment responses. However, little is known about how the interaction between neoantigens and T cells ultimately affects the evolution of cancerous masses. Here, we develop a hierarchical Bayesian model, named neoantigen-T cell interaction estimation (netie) to infer the history of neoantigen-CD8+ T cell interactions in tumors. Netie was systematically validated and applied to examine the molecular patterns of 3,219 tumors, compiled from a panel of 18 cancer types. We showed that tumors with an increase in immune selection pressure over time are associated with T cells that have an activation-related expression signature. We also identified a subset of exhausted cytotoxic T cells postimmunotherapy associated with tumor clones that newly arise after treatment. These analyses demonstrate how netie enables the interrogation of the relationship between individual neoantigen repertoires and the tumor molecular profiles. We found that a T cell inflammation gene expression profile (TIGEP) is more predictive of patient outcomes in the tumors with an increase in immune pressure over time, which reveals a curious synergy between T cells and neoantigen distributions. Overall, we provide a new tool that is capable of revealing the imprints left by neoantigens during each tumor's developmental process and of predicting how tumors will progress under further pressure of the host's immune system.