T cell receptor gene recombinations in human tumor specimen exome files: detection of T cell receptor-β VDJ recombinations associates with a favorable oncologic outcome for bladder cancer

IGL CDR3 Hydropathy and Antigen Chemical Complementarity Associated with Greater Disease-Free Survival in Lung Adenocarcinoma: Implications for Gender Disparities

With lung cancer remaining a challenging disease, new approaches to biomarker discovery and therapy development are needed. Recent immunogenomics, adaptive immune receptor approaches have indicated that it is very likely that B cells play an important role in mediating better overall outcomes. As such, we assessed physicochemical features of lung adenocarcinoma resident IGL complementarity determining region-3 (CDR3) amino acid (AA) sequences and determined that hydrophobic CDR3 AA sequences were associated with a better disease-free survival (DFS) probability. Further, using a recently developed chemical complementarity scoring algorithm particularly suitable for the evaluation of large patient datasets, we determined that IGL CDR3 chemical complementarity with certain cancer testis antigens was associated with better DFS. Chemical complementarity scores for IGL CDR3-MAGEC1 represented a gender bias, with an overrepresentation of males among the higher IGL-CDR3-CTA complementarity scores that were in turn associated with better DFS (logrank p < 0.065). Overall, this study pointed towards potential biomarkers for prognoses that, in some cases are likely gender-specific; and towards biomarkers for guiding therapy, e.g., IGL-based opportunities for antigen targeting in the lung cancer setting.

Delineation of a T-cell receptor CDR3-cancer mutanome aromaticity factor, assessable via blood samples, that facilitates the establishment of survival distinctions in bladder cancer

PURPOSE

A very large and still expanding collection of adaptive immune receptor (IR) recombination reads, representing many diseases, is becoming available for downstream analyses. Among the most productive approaches has been to establish risk stratification parameters via the chemical features of the IR complementarity determining region-3 (CDR3) amino acid (AA) sequences, particularly for large datasets where clinical information is available. Because the IR CDR3 AA sequences often play a large role in antigen binding, the chemistry of these AAs has the likelihood of representing a disease-related fingerprint as well as providing pre-screening information for candidate antigens. To approach this issue in a novel manner, we developed a bladder cancer, case evaluation approach based on CDR3 aromaticity.

METHODS

We developed and applied a simple and efficient algorithm for assessing aromatic, chemical complementarity between T-cell receptor (TCR) CDR3 AA sequences and the cancer specimen mutanome.

RESULTS

Results indicated a survival distinction for aromatic CDR3-aromatic mutanome complementary, versus non-complementary, bladder cancer case sets. This result applied to both tumor resident and blood TCR CDR3 AA sequences and was supported by CDR3 AA sequences represented by both exome and RNAseq files.

CONCLUSION

The described aromaticity factor algorithm has the potential of assisting in prognostic assessments and guiding immunotherapies for bladder cancer.

Tumor resident, TRA anti-viral CDR3 chemical sequence motifs are associated with a better breast cancer outcome

While for certain cancers, such as cervical cancer, the link to viral infections is very strong and very clear, other cancers represent a history of links to viral infections that are either co-morbidities or drive the cancer in ways that are not yet fully understood, for example the "hit and run" possibility. To further understand the connection of viral infections and the progress of breast cancer, we identified the chemical features of known anti-viral, T-cell receptor alpha chain (TRA) complementarity determining region-3 (CDR3) amino acid sequences among the CDR3s of breast cancer patient TRA recombinations and assessed the association of those features with patient outcomes. The application of this novel paradigm indicated consistent associations of tumor-derived, anti-CMV CDR3 chemical sequence motifs with better breast cancer patient outcomes but did not indicate an opportunity to establish risk stratifications for other cancer types. Interestingly, breast cancer samples with no detectable TRA recombinations represented a better outcome than samples with the non-anti-CMV CDR3s, further adding to a rapidly developing series of results allowing a distinction between positive and possibly harmful cancer immune responses.

Chemical complementarity between immune receptors and cancer mutants, independent of antigen presentation protein binding, is associated with increased survival rates

•

Establishment of an immunological distinction between endometrioid and serous uterine cancers.

•

High priority CDR3s, mutant amino acids (AA) for endometrioid cancer prognosis, therapy tools.

•

Further understanding of CDR3-mutant AA complementarity scoring factors, such as HLA binding.

Uterine cancer has been associated with a T-cell immune response that leads to increased survival. Therefore, we used several bioinformatics approaches to explore specific interactions between T-cell receptor (TCR) and tumor mutant peptide sequences. Using endometrioid uterine cancer exome files from the The Cancer Genome Atlas database, we obtained tumor resident V-J recombinations for the T-Cell Receptor alpha gene (TRA). The charged-based, chemical complementarity for each patient's LRP2 or TTN mutant amino acids (AAs) and the recovered, TRA complementarity determining region-3 (CDR3) sequences was calculated, allowing a division of patients into complementary and noncomplementary groups. Complementary groups with TTN mutants had increased disease-free survival and increased expression of complement genes. Furthermore, the survival distinction based on CDR3-mutant peptide complementarity was independent of programmatically assessed HLA class II binding and was not observable based on the CDR3 AA chemical features alone. The above approach provides a potential, highly efficient method for identifying TCR targets in uterine cancer and may aid in the development of novel prognostic tools.

Antiviral T Cell Receptor Complementarity Determining Region-3 Sequences Are Associated with a Worse Cancer Outcome: A Pancancer Analysis

Human papilloma virus has a clearly demonstrated role in cervical and head and neck cancers, but viral etiology for other solid tumors is less well understood. To expand this area of research, we obtained and analyzed the immune receptor recombinations available from both blood and tumor samples, through mining of exome files produced from those sources, for 32 cancer types represented by the cancer genome atlas (TCGA). Among TCGA data sets, the recovery frequency for antiviral complementarity determining region-3 sequences (CDR3s), for T cell receptor-alpha and T cell receptor-beta, ranged from 0% to 21% of the patients, for the different cancer types, with breast, lung, pancreatic, and thymus cancers representing the highest of that range, particularly for tumor tissue resident T cells. In several cases, recovery of the antiviral CDR3s associated with distinct survival rates, and in all of these cases, the recovery of an antiviral CDR3 associated with a worse survival rate.

A scoring system for the electrostatic complementarities of T-cell receptors and cancer-mutant amino acids: multi-cancer analyses of associated survival rates

The anti-tumor immune response is considered to be due to the T-cell receptor (TCR) binding to tumor antigens, which can be either wild-type, early stem cell proteins, presumably foreign to a developed immune system; or mutant peptides, foreign to the immune system because of a mutant amino acid (aa) or otherwise somatically altered aa sequence. Recently, very large numbers of TCR complementarity-determining region-3 (CDR3) aa sequences obtained from tumor specimens have become available. We developed a novel algorithm for assessing the complementarity of tumor mutant peptides and TCR CDR3s, based on the retrieval of TCR CDR3 aa sequences from both tumor specimen and patient blood exomes and by using an automated process of assessing CDR3 and mutant aa electrical charges. Results indicated many instances where high electrostatic complementarity was associated with a higher survival rate. In particular, our approach led to the identification of specific genes contributing significantly to the complementary, TCR CDR3-mutant aa. These results suggest a novel approach to tumor immunoscoring and may lead to the identification of high-priority neo-antigen, peptide vaccines; or to the identification of ex vivo stimulants of tumor-infiltrating lymphocytes.

TRAV gene segments further away from the TRAJ gene segment cluster appear more commonly in human tumor and blood samples

We considered the possibility that the greater the distance between an immune receptor V and J, the more likely the V usage. Such a hypothesis is supported by results from mouse experiments. And, such a hypothesis is consistent with the fundamental nature of recombination and genomic distance: the further the distance, the greater the chance of a DNA break. Thus, we exploited the vast dataset of V and J recombination reads available for the human TRA gene, particularly from cancer and blood specimens, to assess the frequency of TRAV usage with respect to distance from the TRAJ cluster. Results indicated that, indeed, over the entire TRAV cluster, there is a greater chance of V usage the further the distance from the J cluster. These results do not address causation, and are not consistent for certain individual V gene segments, but the results do indicate that overall, the larger the distance between the V and J gene segment cluster, the more likely the appearance of at least a subset of TRAV segments, particularly among tumor infiltrating lymphocytes. With a similar approach, the distal TRAV gene segments were also found to be more commonly associated with a subset of distal TRAJ segments. These results have implications for restrictions on the apparent TRA repertoire in disease settings.

T-cell receptor-α CDR3 domain chemical features correlate with survival rates in bladder cancer

PURPOSE

In certain cancer settings, a T-cell response to cancer represents a relatively favorable outcome. Thus, the near-future challenges include a better understanding of exactly which T-cell features contribute to a response to which cancer antigen-groups, to maximize the opportunities for tumor-infiltrating lymphocyte (TIL)-based therapies and other immunotherapies.

METHODS

The immune receptor complementarity determining region-3 (CDR3) is considered to be important for antigen binding, hence, in this report, we evaluated the chemical features of the CDR3 of 846 T-cell receptor-α (TCR-α) coding regions associated with bladder tumor tissue, using bioinformatics databases.

RESULTS

Results indicated that statistically significantly distinct, low value, CDR3 region isoelectric points associate with a better outcome (log rank p < 0.027, overall survival). Moreover, in samples representing the more favorable isoelectric points, known driver mutations, for example, PIK3CA (E → K) with chemically complementary features overlap the better-outcome, low isoelectric point samples. Further work extended these results, i.e., survival rate-CDR3 associations, to other CDR3 chemical features and other cancers, consistent with the initial isoelectric point-related, bladder cancer findings.

CONCLUSIONS

A bioinformatics assessment of cancer-associated TCR biochemical features may improve the accuracy of the predictions of which TILs will be best for ex-vivo amplification and which patients will benefit from other immunotherapies.

TcR-α recombinations in renal cell carcinoma exome files correlate with an intermediate level of T-cell exhaustion biomarkers

Renal cell carcinoma exome-derived, V(D)J recombination reads had an elevated presence and variability, for both TcR-α and -β, when compared to marginal tissue, reflecting an opportunity to assess tumor immunogenicity by comparison with marginal tissue T cells. PD-1, PD-L2, CTLA4 and FOXP3, all of which are implicated in the evasion of an anti-tumor immune response, had a significantly higher expression for samples representing co-detection of productive TcR-α and -β recombination reads. Samples representing tumors with productive TcR-α recombination reads but no detectable, productive TcR-β recombination reads, reflected a 20% survival advantage, and RNASeq data indicated an intermediate level of immune checkpoint gene expression for those samples. These results raise the question of whether relatively high levels of detection of productive TcR-α recombination reads, in comparison with detection of reads representing the TcR-β gene, identify a microenvironment that has not yet entered a T-cell exhaustion phase and may thereby represent conditions for immune enhancements that do not require anti-immune checkpoint therapies.

T cell receptor-β J usage, in combination with particular HLA class II alleles, correlates with better cancer survival rates

T cell receptor (TCR) β V and J usage correlates with either the HLA class I or HLA class II major histocompatibility subtypes, and in both infectious diseases and autoimmune settings, the use of particular TCR-β V and J's, in persons with specific HLA alleles, represents either better outcomes or certain clinical features. However, the relationship of TCR V and J usage, HLA alleles, and clinical parameters in the cancer setting has been less well studied. Here, we have evaluated the relationship of what is likely dominant TCR-β V and J usage among tissue-resident lymphocytes for lung, head and neck, kidney, stomach, ovarian, and endometrial cancers, with patient HLA class II alleles. The most striking indication is that TCR-β J subgroup usage, in combination with particular patient HLA class II alleles, correlated with either better or worse outcomes for lung cancer. One combination, TCR-β J2 segment usage and the HLA-DRB1*1501 allele, correlated with a better survival rate for both lung and head and neck cancers. These results fill a gap in knowledge regarding the relevance of HLA typing to cancer and indicate that HLA typing, along with an indication of dominant TCR-β J usage among tissue-resident lymphocytes, can be useful for prognosis.

TRB-J1 usage, in combination with the HLA-A*01:01 allele, represents an apparent survival advantage for uterine corpus endometrial carcinoma: Comparisons with microscopic assessments of lymphocyte infiltrates

The opportunity for the highly efficient recovery of immune receptor recombination data from cancer specimens, including the ready assessment of immune receptor V and J usage, raises the issue of establishing precise values of assessing the immune receptor status as opposed to obtaining basic information regarding lymphocyte infiltration, in the cancer setting. In this report, we obtained the lymphocyte infiltration percentages from the cancer digital slide archive representing uterine corpus endometrial carcinoma (UCEC) and correlated these data with recovery of the immune receptor recombination reads from corresponding UCEC exome files. Results indicated a basic correlation of the recovery of productive T-cell receptor beta (TRB) recombination reads with lymphocyte infiltration percentages. However, the recovery of specific immune receptor recombination reads did not indicate the same survival outcomes as microscope detection of lymphocyte infiltrate percentages. To further exploit the value of recovery of the TRB recombination reads from the UCEC exome files, we determined the survival outcomes for combinations of TRB gene segment usage and HLA class I alleles, with the most important result being that the combination of HLA-A*01:01 and TRB-J1 segment usage reflected a strikingly high survival rate. Overall, this report emphasized the increased value of the knowledge of the immune receptor recombinations, in comparison with basic lymphocyte infiltration percentages, in assessing cancer survival rates.

Mutant cytoskeletal and ECM peptides sensitive to the ST14 protease are associated with a worse outcome for glioblastoma multiforme

We previously identified a set of the most frequently mutated cytoskeleton- and extracellular matrix-related proteins (CECMPs) in numerous cancer datasets. In this report, we used a bioinformatics approach to assess the impact of amino acid (AA) substitutions on the sensitivity of CECMPs to the ST14 protease (matriptase I), a transmembrane serine protease previously implicated in cancer development. Results indicated that AA substitutions in glioblastoma multiforme (GBM) CECMPs are skewed toward increased resistance to the ST14 protease, in comparison to the wild-type peptide sequence. Furthermore, the protease resistant AA substitutions represent relatively high binding affinities to HLA class I proteins, when assessing the binding specificities using HLA class I alleles matched to the source of the mutant AA. Moreover, samples representing AA substitutions that increased protease sensitivity also represented reduced overall and disease-free survival periods for patients with glioblastoma. To assess tumor specimen immunogenicity, we identified T-cell receptor (TCR) V(D)J recombinations in GBM exome files. The overlap between ST14 protease sensitive mutant barcodes and the TCR V(D)J recombination read positive barcodes represented significantly reduced survival.

T-cell receptor-β V and J usage, in combination with particular HLA class I and class II alleles, correlates with cancer survival patterns

Class I and class II HLA proteins, respectively, have been associated with subsets of V(D)J usage resulting from recombination of the T-cell receptor (TCR) genes. Additionally, particular HLA alleles, in combination with dominant TCR V(D)J recombinations, have been associated with several autoimmune diseases. The recovery of TCR recombination reads from tumor specimen exome files has allowed rapid and extensive assessments of V(D)J usage, likely for cancer resident T-cells, across relatively large cancer datasets. The results from this approach, in this report, have permitted an extensive alignment of TCR-β VDJ usage and HLA class I and II alleles. Results indicate the correlation of both better and worse cancer survival rates with particular TCR-β, V and J usage-HLA allele combinations, with differences in median survival times ranging from 7 to 130 months, depending on the cancer and the specific TCR-β V and J usage/HLA class allele combination.

Recovery of Immunoglobulin VJ Recombinations from Pancreatic Cancer Exome Files Strongly Correlates with Reduced Survival

We assessed pancreatic cancer, lymphocyte infiltrates with a computational genomics approach. We took advantage of tumor-specimen exome files available from the cancer genome atlas to mine T- and B-cell immune receptor recombinations, using highly efficient, scripted algorithms established in several previous reports. Surprisingly, the results indicated that pancreatic cancer exomes represent one of the highest level yields for immune receptor recombinations, significantly higher than two comparison cancers used in this study, head and neck and bladder cancer. In particular, pancreatic cancer exomes have very large numbers of immunoglobulin light chain recombinations, both with regard to number of samples characterized by recovery of such recombinations and with regard to numbers of recombination reads per sample. These results were consistent with B-cell biomarkers, which emphasized the Th2 nature of the pancreatic lymphocyte infiltrate. The tumor specimen exomes with B-cell immune receptor recombination reads represented a dramatically poor outcome, a result not detected with either the head and neck or bladder cancer datasets. The results presented here support the potential value of immunotherapies designed to engineer a Th2 to Th1 shift in treating certain forms of pancreatic cancer.

Immunogenomics: A Negative Prostate Cancer Outcome Associated with TcR-γ/δ Recombinations

We developed a scripted algorithm, based on previous, earlier editions of the algorithm, to mine prostate cancer exome files for T-cell receptor (TcR) recombination reads: Reads representing TcR gene recombinations were identified in 497 prostate cancer exome files from the cancer genome atlas (TCGA). As has been reported for melanoma, co-detection of productive TcR-α and TcR-β recombination reads correlated with an RNA expression signature representing T-cell exhaustion, particularly with high RNA levels for PD-1 and PD-L1, in comparison to several different control sets of samples. Co-detection of TcR-α and TcR-β recombination reads also correlated with high level expression of genes representing antigen presenting functions, further supporting the conclusion that co-detection of TcR-α and TcR-β recombination reads represents an immunologically relevant microenvironment. Finally, detection of unproductive TcR-δ recombinations, and unproductive and productive TcR-γ recombinations, strongly correlated with, and may represent a convenient biomarker for a poor clinical outcome. These results underscore the value of the genomics-based assessment of unproductive TcR recombinations and raise questions about the impact of tumor microenvironment lymphocytes in the absence of antigenicity.

Lung tumor exome files with T-cell receptor recombinations: a mouse model of T-cell infiltrates reflecting mutation burdens

Tumor exomes and RNASeq data were originally intended for obtaining tumor mutations and gene expression profiles, respectively. However, recent work has determined that tumor exome and RNAseq read files contain reads representing T-cell and B-cell receptor (TcR and BcR) recombinations, presumably due to infiltrating lymphocytes. Furthermore, the recovery of immune receptor recombination reads has demonstrated correlations with specific, previously appreciated aspects of tumor immunology. To further understand the usefulness of recovering TcR and BcR recombinations from tumor exome files, we developed a scripted algorithm for recovery of reads representing these recombinations from a previously described mouse model of lung tumorigenesis. Results indicated that exomes representing lung adenomas reveal significantly more TcR recombinations than do exomes from lung adenocarcinomas; and that exome files representing high mutation adenomas, arising from chemical mutagens, have more TcR recombinations than do exome files from low mutation adenomas arising from an activating Kras mutation. The latter results were also consistent with a similar analysis performed on human lung adenocarcinoma exomes. The mouse and human results for obtaining TcR recombination reads from tumor specimen exomes are consistent with human tumor biology results indicating that adenomas and high mutation cancers are sites of high immune activity. The results indicate hitherto unappreciated opportunities for the use of tumor specimen exome files, particularly from experimental animal models, to study the connection between the adenoma stage of tumorigenesis, or high cancer mutation rates, and high level lymphocyte infiltrates.

Protected cytoskeletal-related proteins: Towards a resolution of contradictions regarding the role of the cytoskeleton in cancer

Initial reports of the role of the cytoskeleton in cancer indicated that tumor cells with a more disorganized cytoskeleton were more tumorigenic. These reports were based on stains for the F-actin cytoskeleton, for example, using phalloidin or anti-F-actin antibody reagents, and gave a basic impression of F-actin-based cytoskeletal integrity. Later developments emphasized the significance of the cytoskeletal elements in cell migration, presumably associated with either basement membrane invasion or metastasis, or both, with several specific proteins implicated in the formation of cell invadopodia. With the advent of genomics approaches, it has become clear that cytoskeletal related proteins are indeed common targets of mutagenesis in cancer and commonly rank among the most mutated proteins in cancers, presumably due to large coding region sizes and the significant stochastic component to human mutagenesis. This cytoskeletal genomics result is consistent with the loss of cytoskeleton integrity as a hallmark of tumor development, but raises the question of whether such mutational sensitivity relates to the migration and invadopodia aspects of tumor progression. In the present study, the authors report that it is possible to identify a set of cytoskeletal related proteins protected from mutation, in comparison to the commonly mutated cytoskeleton related proteins in certain, but not all cancer, datasets.