Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping.

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

Binding peptide generation for MHC Class I proteins with deep reinforcement learning.

MOTIVATION: MHC Class I protein plays an important role in immunotherapy by presenting immunogenic peptides to anti-tumor immune cells. The repertoires of peptides for various MHC Class I proteins are distinct, which can be reflected by their diverse binding motifs. To characterize binding motifs for MHC Class I proteins, in vitro experiments have been conducted to screen peptides with high binding affinities to hundreds of given MHC Class I proteins. However, considering tens of thousands of known MHC Class I proteins, conducting in vitro experiments for extensive MHC proteins is infeasible, and thus a more efficient and scalable way to characterize binding motifs is needed. RESULTS: We presented a de novo generation framework, coined PepPPO, to characterize binding motif for any given MHC Class I proteins via generating repertoires of peptides presented by them. PepPPO leverages a reinforcement learning agent with a mutation policy to mutate random input peptides into positive presented ones. Using PepPPO, we characterized binding motifs for around 10 000 known human MHC Class I proteins with and without experimental data. These computed motifs demonstrated high similarities with those derived from experimental data. In addition, we found that the motifs could be used for the rapid screening of neoantigens at a much lower time cost than previous deep-learning methods. AVAILABILITY AND IMPLEMENTATION: The software can be found in https://github.com/minrq/pMHC. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of the T cell receptor repertoire and melanoma tumor microenvironment upon combined treatment with ipilimumab and hTERT vaccination.

BACKGROUND: This clinical trial evaluated a novel telomerase-targeting therapeutic cancer vaccine, UV1, in combination with ipilimumab, in patients with metastatic melanoma. Translational research was conducted on patient-derived blood and tissue samples with the goal of elucidating the effects of treatment on the T cell receptor repertoire and tumor microenvironment. METHODS: The trial was an open-label, single-center phase I/IIa study. Eligible patients had unresectable metastatic melanoma. Patients received up to 9 UV1 vaccinations and four ipilimumab infusions. Clinical responses were assessed according to RECIST 1.1. Patients were followed up for progression-free survival (PFS) and overall survival (OS). Whole-exome and RNA sequencing, and multiplex immunofluorescence were performed on the biopsies. T cell receptor (TCR) sequencing was performed on the peripheral blood and tumor tissues. RESULTS: Twelve patients were enrolled in the study. Vaccine-specific immune responses were detected in 91% of evaluable patients. Clinical responses were observed in four patients. The mPFS was 6.7 months, and the mOS was 66.3 months. There was no association between baseline tumor mutational burden, neoantigen load, IFN-γ gene signature, tumor-infiltrating lymphocytes, and response to therapy. Tumor telomerase expression was confirmed in all available biopsies. Vaccine-enriched TCR clones were detected in blood and biopsy, and an increase in the tumor IFN-γ gene signature was detected in clinically responding patients. CONCLUSION: Clinical responses were observed irrespective of established predictive biomarkers for checkpoint inhibitor efficacy, indicating an added benefit of the vaccine-induced T cells. The clinical and immunological read-out warrants further investigation of UV1 in combination with checkpoint inhibitors. Trial registration Clinicaltrials.gov identifier: NCT02275416. Registered October 27, 2014. https://clinicaltrials.gov/ct2/show/NCT02275416?term=uv1&draw=2&rank=6.

Characterization of the Small Intestinal Lesion in Celiac Disease by Label-Free Quantitative Mass Spectrometry.

Global characterization of tissue proteomes from small amounts of biopsy material has become feasible because of advances in mass spectrometry and bioinformatics tools. In celiac disease (CD), dietary gluten induces an immune response that is accompanied by pronounced remodeling of the small intestine. Removal of gluten from the diet abrogates the immune response, and the tissue architecture normalizes. In this study, differences in global protein expression of small intestinal biopsy specimens from CD patients were quantified by analyzing formalin-fixed, paraffin-embedded material using liquid chromatography-mass spectrometry and label-free protein quantitation. Protein expression was compared in biopsy specimens collected from the same patients before and after 1-year treatment with gluten-free diet (n = 10) or before and after 3-day gluten provocation (n = 4). Differential expression of proteins in particular from mature enterocytes, neutrophils, and plasma cells could distinguish untreated from treated CD mucosa, and Ig variable region IGHV5-51 expression was found to serve as a CD-specific marker of ongoing immune activation. In patients who had undergone gluten challenge, coordinated up-regulation of wound response proteins, including the CD autoantigen transglutaminase 2, was observed. Our study provides a global and unbiased assessment of antigen-driven changes in protein expression in the celiac intestinal mucosa.

Naive Donor NK Cell Repertoires Associated with Less Leukemia Relapse after Allogeneic Hematopoietic Stem Cell Transplantation.

Acute and latent human CMV cause profound changes in the NK cell repertoire, with expansion and differentiation of educated NK cells expressing self-specific inhibitory killer cell Ig-like receptors. In this study, we addressed whether such CMV-induced imprints on the donor NK cell repertoire influenced the outcome of allogeneic stem cell transplantation. Hierarchical clustering of high-resolution immunophenotyping data covering key NK cell parameters, including frequencies of CD56(bright), NKG2A(+), NKG2C(+), and CD57(+) NK cell subsets, as well as the size of the educated NK cell subset, was linked to clinical outcomes. Clusters defining naive (NKG2A(+)CD57(-)NKG2C(-)) NK cell repertoires in the donor were associated with decreased risk for relapse in recipients with acute myeloid leukemia and myelodysplastic syndrome (hazard ratio [HR], 0.09; 95% confidence interval [CI]: 0.03-0.27; p < 0.001). Furthermore, recipients with naive repertoires at 9-12 mo after hematopoietic stem cell transplantation had increased disease-free survival (HR, 7.2; 95% CI: 1.6-33; p = 0.01) and increased overall survival (HR, 9.3; 95% CI: 1.1-77, p = 0.04). Conversely, patients with a relative increase in differentiated NK cells at 9-12 mo displayed a higher rate of late relapses (HR, 8.41; 95% CI: 6.7-11; p = 0.02), reduced disease-free survival (HR, 0.12; 95% CI: 0.12-0.74; p = 0.02), and reduced overall survival (HR, 0.07; 95% CI: 0.01-0.69; p = 0.02). Thus, our data suggest that naive donor NK cell repertoires are associated with protection against leukemia relapse after allogeneic HSCT.

Profiling networks of distinct immune-cells in tumors.

BACKGROUND: It is now clearly evident that cancer outcome and response to therapy is guided by diverse immune-cell activity in tumors. Presently, a key challenge is to comprehensively identify networks of distinct immune-cell signatures present in complex tissue, at higher-resolution and at various stages of differentiation, activation or function. This is particularly so for closely related immune-cells with diminutive, yet critical, differences. RESULTS: To predict networks of infiltrated distinct immune-cell phenotypes at higher resolution, we explored an integrated knowledge-based approach to select immune-cell signature genes integrating not only expression enrichment across immune-cells, but also an automatic capture of relevant immune-cell signature genes from the literature. This knowledge-based approach was integrated with resources of immune-cell specific protein networks, to define signature genes of distinct immune-cell phenotypes. We demonstrate the utility of this approach by profiling signatures of distinct immune-cells, and networks of immune-cells, from metastatic melanoma patients who had undergone chemotherapy. The resultant bioinformatics strategy complements immunohistochemistry from these tumors, and predicts both tumor-killing and immunosuppressive networks of distinct immune-cells in responders and non-responders, respectively. The approach is also shown to capture differences in the immune-cell networks of BRAF versus NRAS mutated metastatic melanomas, and the dynamic changes in resistance to targeted kinase inhibitors in MAPK signalling. CONCLUSIONS: This integrative bioinformatics approach demonstrates that capturing the protein network signatures and ratios of distinct immune-cell in the tumor microenvironment maybe an important factor in predicting response to therapy. This may serve as a computational strategy to define network signatures of distinct immune-cells to guide immuno-pathological discovery.

Hypoxia-independent gene expression signature associated with radiosensitisation of prostate cancer cell lines by histone deacetylase inhibition.

BACKGROUND: Histone deacetylase inhibitors (HDACis) like vorinostat are promising radiosensitisers in prostate cancer, but their effect under hypoxia is not known. We investigated gene expression associated with radiosensitisation of normoxic and hypoxic prostate cancer cells by vorinostat. METHODS: Cells were exposed to vorinostat under normoxia or hypoxia and subjected to gene expression profiling before irradiation and clonogenic survival analysis. RESULTS: Pretreatment with vorinostat led to radiosensitisation of the intrinsically radioresistant DU 145 cells, but not the radiosensitive PC-3 and 22Rv1 cells, and was independent of hypoxia status. Knockdown experiments showed that the sensitisation was not caused by repression of hypoxia-inducible factor HIF1 or tumour protein TP53. Global deregulation of DNA repair and chromatin organisation genes was associated with radiosensitisation under both normoxia and hypoxia. A radiosensitisation signature with expression changes of 56 genes was generated and valid for both conditions. For eight signature genes, baseline expression also correlated with sensitisation, showing potential as pretreatment biomarker. The hypoxia independence of the signature was confirmed in a clinical data set. CONCLUSIONS: Pretreatment with HDACi may overcome radioresistance of hypoxic prostate tumours by similar mechanisms as under normoxia. We propose a gene signature to predict radiosensitising effects independent of hypoxia status.

Predicting physical interactions between protein complexes.

Protein complexes enact most biochemical functions in the cell. Dynamic interactions between protein complexes are frequent in many cellular processes. As they are often of a transient nature, they may be difficult to detect using current genome-wide screens. Here, we describe a method to computationally predict physical interactions between protein complexes, applied to both humans and yeast. We integrated manually curated protein complexes and physical protein interaction networks, and we designed a statistical method to identify pairs of protein complexes where the number of protein interactions between a complex pair is due to an actual physical interaction between the complexes. An evaluation against manually curated physical complex-complex interactions in yeast revealed that 50% of these interactions could be predicted in this manner. A community network analysis of the highest scoring pairs revealed a biologically sensible organization of physical complex-complex interactions in the cell. Such analyses of proteomes may serve as a guide to the discovery of novel functional cellular relationships.

The Genomic HyperBrowser: an analysis web server for genome-scale data.

The immense increase in availability of genomic scale datasets, such as those provided by the ENCODE and Roadmap Epigenomics projects, presents unprecedented opportunities for individual researchers to pose novel falsifiable biological questions. With this opportunity, however, researchers are faced with the challenge of how to best analyze and interpret their genome-scale datasets. A powerful way of representing genome-scale data is as feature-specific coordinates relative to reference genome assemblies, i.e. as genomic tracks. The Genomic HyperBrowser (http://hyperbrowser.uio.no) is an open-ended web server for the analysis of genomic track data. Through the provision of several highly customizable components for processing and statistical analysis of genomic tracks, the HyperBrowser opens for a range of genomic investigations, related to, e.g., gene regulation, disease association or epigenetic modifications of the genome.

From proteomes to complexomes in the era of systems biology.

Protein complexes carry out almost the entire signaling and functional processes in the cell. The protein complex complement of a cell, and its network of complex-complex interactions, is referred to here as the complexome. Computational methods to predict protein complexes from proteomics data, resulting in network representations of complexomes, have recently being developed. In addition, key advances have been made toward understanding the network and structural organization of complexomes. We review these bioinformatics advances, and their discovery-potential, as well as the merits of integrating proteomics data with emerging methods in systems biology to study protein complex signaling. It is envisioned that improved integration of proteomics and systems biology, incorporating the dynamics of protein complexes in space and time, may lead to more predictive models of cell signaling networks for effective modulation.

Identifying pathogenic processes by integrating microarray data with prior knowledge.

BACKGROUND: It is of great importance to identify molecular processes and pathways that are involved in disease etiology. Although there has been an extensive use of various high-throughput methods for this task, pathogenic pathways are still not completely understood. Often the set of genes or proteins identified as altered in genome-wide screens show a poor overlap with canonical disease pathways. These findings are difficult to interpret, yet crucial in order to improve the understanding of the molecular processes underlying the disease progression. We present a novel method for identifying groups of connected molecules from a set of differentially expressed genes. These groups represent functional modules sharing common cellular function and involve signaling and regulatory events. Specifically, our method makes use of Bayesian statistics to identify groups of co-regulated genes based on the microarray data, where external information about molecular interactions and connections are used as priors in the group assignments. Markov chain Monte Carlo sampling is used to search for the most reliable grouping. RESULTS: Simulation results showed that the method improved the ability of identifying correct groups compared to traditional clustering, especially for small sample sizes. Applied to a microarray heart failure dataset the method found one large cluster with several genes important for the structure of the extracellular matrix and a smaller group with many genes involved in carbohydrate metabolism. The method was also applied to a microarray dataset on melanoma cancer patients with or without metastasis, where the main cluster was dominated by genes related to keratinocyte differentiation. CONCLUSION: Our method found clusters overlapping with known pathogenic processes, but also pointed to new connections extending beyond the classical pathways.

Identification of eight candidate target genes of the recurrent 3p12-p14 loss in cervical cancer by integrative genomic profiling.

The pathogenetic role, including its target genes, of the recurrent 3p12-p14 loss in cervical cancer has remained unclear. To determine the onset of the event during carcinogenesis, we used microarray techniques and found that the loss was the most frequent 3p event, occurring in 61% of 92 invasive carcinomas, in only 2% of 43 high-grade intraepithelial lesions (CIN2/3), and in 33% of 6 CIN3 lesions adjacent to invasive carcinomas, suggesting a role in acquisition of invasiveness or early during the invasive phase. We performed an integrative DNA copy number and expression analysis of 77 invasive carcinomas, where all genes within the recurrent region were included. We selected eight genes, THOC7, PSMD6, SLC25A26, TMF1, RYBP, SHQ1, EBLN2, and GBE1, which were highly down-regulated in cases with loss, as confirmed at the protein level for RYBP and TMF1 by immunohistochemistry. The eight genes were subjected to network analysis based on the expression profiles, revealing interaction partners of proteins encoded by the genes that were coordinately regulated in tumours with loss. Several partners were shared among the eight genes, indicating crosstalk in their signalling. Gene ontology analysis showed enrichment of biological processes such as apoptosis, proliferation, and stress response in the network and suggested a relationship between down-regulation of the eight genes and activation of tumourigenic pathways. Survival analysis showed prognostic impact of the eight-gene signature that was confirmed in a validation cohort of 74 patients and was independent of clinical parameters. These results support the role of the eight candidate genes as targets of the 3p12-p14 loss in cervical cancer and suggest that the strong selection advantage of the loss during carcinogenesis might be caused by a synergetic effect of several tumourigenic processes controlled by these targets.

Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response.

BACKGROUND: The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer. METHODS: RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles. RESULTS: The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls. CONCLUSIONS: The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.

LTX-315 and adoptive cell therapy using tumor-infiltrating lymphocytes generate tumor specific T cells in patients with metastatic soft tissue sarcoma.

LTX-315 is an oncolytic peptide that elicits both local and systemic immune responses upon intratumoral injection. In the present pilot trial, we treated patients with metastatic soft tissue sarcoma with the combination of LTX-315 and adoptive T-cell therapy using in vitro expanded tumor-infiltrating lymphocytes. Six heavily pretreated patients were included in the trial and treated with LTX-315 of which four patients proceeded to adoptive T-cell therapy. Overall, the treatment was considered safe with only expected and manageable toxicity. The best overall clinical response was stable disease for 208 days, and in this patient, we detected tumor-reactive T cells in the blood that lasted until disease progression. In three patients T-cell reactivity against in silico predicted neoantigens was demonstrated. Additionally, de novo T-cell clones were generated and expanded in the blood following LTX-315 injections. In conclusion, this pilot study provides proof that it is feasible to combine LTX-315 and adoptive T-cell therapy, and that this treatment can induce systemic immune responses that resulted in stabilization of the disease in sarcoma patients with otherwise progressive disease. Further optimization of the treatment protocol is warranted to increase clinical activity. ClinicalTrials.gov Identifier: NCT03725605.

Corrigendum: Experimental validation of immunogenic SARS-CoV-2 T cell epitopes identified by artificial intelligence.

[This corrects the article DOI: 10.3389/fimmu.2023.1265044.].

Candidate tumor-specific CD8+ T cell subsets identified in the malignant pleural effusion of advanced lung cancer patients by single-cell analysis.

Isolation of tumor-specific T cells and their antigen receptors (TCRs) from malignant pleural effusions (MPE) may facilitate the development of TCR-transduced adoptive cellular immunotherapy products for advanced lung cancer patients. However, the characteristics and markers of tumor-specific T-cells in MPE are largely undefined. To this end, to establish the phenotypes and antigen specificities of CD8+ T cells, we performed single-cell RNA and TCR sequencing of samples from three advanced lung cancer patients. Dimensionality reduction on a total of 4,983 CD8+ T cells revealed 10 clusters including naïve, memory, and exhausted phenotypes. We focused particularly on exhausted T cell clusters and tested their TCR reactivity against neoantigens predicted from autologous cancer cell lines. Four different TCRs specific for the same neoantigen and one orphan TCR specific for the autologous cell line were identified from one of the patients. Differential gene expression analysis in tumor-specific T cells relative to the other T cells identified CXCL13, as a candidate gene expressed by tumor-specific T cells. In addition to expressing CXCL13, tumor-specific T cells were present in a higher proportion of T cells co-expressing PDCD1(PD-1)/TNFRSF9(4-1BB). Furthermore, flow cytometric analyses in advanced lung cancer patients with MPE documented that those with high PD-1/4-1BB expression have a better prognosis in the subset of 57 adenocarcinoma patients (p = .039). These data suggest that PD-1/4-1BB co-expression might identify tumor-specific CD8+ T cells in MPE, which are associated with patients' prognosis. (233 words).

Clinical Activity of Combined Telomerase Vaccination and Pembrolizumab in Advanced Melanoma: Results from a Phase I Trial.

PURPOSE: Cancer vaccines represent a novel treatment modality with a complementary mode of action addressing a crucial bottleneck for checkpoint inhibitor (CPI) efficacy. CPIs are expected to release brakes in T-cell responses elicited by vaccination, leading to more robust immune responses. Increased antitumor T-cell responses may confer increased antitumor activity in patients with less immunogenic tumors, a subgroup expected to achieve reduced benefit from CPIs alone. In this trial, a telomerase-based vaccine was combined with pembrolizumab to assess the safety and clinical activity in patients with melanoma. PATIENTS AND METHODS: Thirty treatment-naïve patients with advanced melanoma were enrolled. Patients received intradermal injections of UV1 with adjuvant GM-CSF at two dose levels, and pembrolizumab according to the label. Blood samples were assessed for vaccine-induced T-cell responses, and tumor tissues were collected for translational analyses. The primary endpoint was safety, with secondary objectives including progression-free survival (PFS), overall survival (OS), and objective response rate (ORR). RESULTS: The combination was considered safe and well-tolerated. Grade 3 adverse events were observed in 20% of patients, with no grade 4 or 5 adverse events reported. Vaccination-related adverse events were mostly mild injection site reactions. The median PFS was 18.9 months, and the 1- and 2-year OS rates were 86.7% and 73.3%, respectively. The ORR was 56.7%, with 33.3% achieving complete responses. Vaccine-induced immune responses were observed in evaluable patients, and inflammatory changes were detected in posttreatment biopsies. CONCLUSIONS: Encouraging safety and preliminary efficacy were observed. Randomized phase II trials are currently ongoing.

Experimental validation of immunogenic SARS-CoV-2 T cell epitopes identified by artificial intelligence.

During the COVID-19 pandemic we utilized an AI-driven T cell epitope prediction tool, the NEC Immune Profiler (NIP) to scrutinize and predict regions of T cell immunogenicity (hotspots) from the entire SARS-CoV-2 viral proteome. These immunogenic regions offer potential for the development of universally protective T cell vaccine candidates. Here, we validated and characterized T cell responses to a set of minimal epitopes from these AI-identified universal hotspots. Utilizing a flow cytometry-based T cell activation-induced marker (AIM) assay, we identified 59 validated screening hits, of which 56% (33 peptides) have not been previously reported. Notably, we found that most of these novel epitopes were derived from the non-spike regions of SARS-CoV-2 (Orf1ab, Orf3a, and E). In addition, ex vivo stimulation with NIP-predicted peptides from the spike protein elicited CD8+ T cell response in PBMC isolated from most vaccinated donors. Our data confirm the predictive accuracy of AI platforms modelling bona fide immunogenicity and provide a novel framework for the evaluation of vaccine-induced T cell responses.

The interplay between neoantigens and immune cells in sarcomas treated with checkpoint inhibition.

Introduction: Sarcomas are comprised of diverse bone and connective tissue tumors with few effective therapeutic options for locally advanced unresectable and/or metastatic disease. Recent advances in immunotherapy, in particular immune checkpoint inhibition (ICI), have shown promising outcomes in several cancer indications. Unfortunately, ICI therapy has provided only modest clinical responses and seems moderately effective in a subset of the diverse subtypes. Methods: To explore the immune parameters governing ICI therapy resistance or immune escape, we performed whole exome sequencing (WES) on tumors and their matched normal blood, in addition to RNA-seq from tumors of 31 sarcoma patients treated with pembrolizumab. We used advanced computational methods to investigate key immune properties, such as neoantigens and immune cell composition in the tumor microenvironment (TME). Results: A multifactorial analysis suggested that expression of high quality neoantigens in the context of specific immune cells in the TME are key prognostic markers of progression-free survival (PFS). The presence of several types of immune cells, including T cells, B cells and macrophages, in the TME were associated with improved PFS. Importantly, we also found the presence of both CD8+ T cells and neoantigens together was associated with improved survival compared to the presence of CD8+ T cells or neoantigens alone. Interestingly, this trend was not identified with the combined presence of CD8+ T cells and TMB; suggesting that a combined CD8+ T cell and neoantigen effect on PFS was important. Discussion: The outcome of this study may inform future trials that may lead to improved outcomes for sarcoma patients treated with ICI.

Robust spike-specific CD4+ and CD8+ T cell responses in SARS-CoV-2 vaccinated hematopoietic cell transplantation recipients: a prospective, cohort study.

Poor overall survival of hematopoietic stem cell transplantation (HSCT) recipients who developed COVID-19 underlies the importance of SARS-CoV-2 vaccination. Previous studies of vaccine efficacy have reported weak humoral responses but conflicting results on T cell immunity. Here, we have examined the relationship between humoral and T cell response in 48 HSCT recipients who received two doses of Moderna's mRNA-1273 or Pfizer/BioNTech's BNT162b2 vaccines. Nearly all HSCT patients had robust T cell immunity regardless of protective humoral responses, with 18/48 (37%, IQR 8.679-5601 BAU/mL) displaying protective IgG anti-receptor binding domain (RBD) levels (>2000 BAU/mL). Flow cytometry analysis of activation induced markers (AIMs) revealed that 90% and 74% of HSCT patients showed reactivity towards immunodominant spike peptides in CD8+ and CD4+ T cells, respectively. The response rate increased to 90% for CD4+ T cells as well when we challenged the cells with a complete set of overlapping peptides spanning the entire spike protein. T cell response was detectable as early as 3 months after transplant, but only CD4+ T cell reactivity correlated with IgG anti-RBD level and time after transplantation. Boosting increased seroconversion rate, while only one patient developed COVID-19 requiring hospitalization. Our data suggest that HSCT recipients with poor serological responses were protected from severe COVID-19 by vaccine-induced T cell responses.