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.

The presence of intratumoral Porphyromonas gingivalis correlates with a previously defined pancreatic adenocarcinoma, immune cell expression phenotype and with tumor resident, adaptive immune receptor features.

The association between pancreatic adenocarcinoma (PAAD) and the pancreatic microbiome is not fully understood, although bacteria may decrease the effectiveness of chemotherapy and lead to anti-apoptotic, pro-inflammatory microenvironments. To better understand the relationship between the PAAD microbiome and the microenvironment, we identified Porphyromonas gingivalis-positive PAAD samples and found a strong association between intratumoral P. gingivalis and: (i) an immune cell gene expression phenotype previously defined by others as gene program 7; and (ii) recovery of immunoglobulin recombination, sequencing reads. We applied a novel chemical complementarity scoring algorithm, suitable for a big data setting, and determined that the previously established P. gingivalis antigen, rpgB had a reduced chemical complementarity with T-cell receptor (TCR) complementarity-determining region-3 amino acid sequences recovered from PAAD samples with P. gingivalis in comparison to TCR-rpgB chemical complementarity represented by the PAAD samples that lacked P. gingivalis. This finding strengthens the existing body of evidence correlating P. gingivalis with PAAD, which may have implications for the treatment and prognosis of patients. Furthermore, demonstrating the correlation of P. gingivalis and gene program 7 raises the question of whether P. gingivalis infection is responsible for the gene program 7 subdivision of PAAD?

A comparison of immune receptor recombination databases sourced from tumour exome or RNAseq files: Verifications of immunological distinctions between primary and metastatic melanoma.

It became apparent several years ago that RNAseq and exome files prepared from tissue could be mined for adaptive immune receptor (IR) recombinations, which has given extra value to datasets originally intended for gene expression or mutation studies. For example, recovery of IR recombination reads from tumour specimen genomics files can correlate with survival rates. In particular, many benchmarking processes have been applied to the two sets of the IR recombination reads obtained from the cancer genome atlas files, but these two sets have never been directly compared. Here we show that both sets largely agree regarding several parameters. For example, recovery of TRB recombination reads from both WXS and RNAseq files representing metastatic melanoma was associated with a better outcome (p < .0004 in both cases); and T-cell receptor recombination read recovery, for both genomics file types, associated very strongly with T-cell gene expression markers. However, the use of CDR3 chemical features for survival distinctions was not consistent. This topic, and the surprising result that both datasets indicated that primary melanoma with recovery of IR recombination reads, in stark contrast to metastatic melanoma, represents a worse outcome, are discussed.

Electrostatic Complementarity of T-Cell Receptor-Alpha CDR3 Domains and Mutant Amino Acids Is Associated with Better Survival Rates for Sarcomas.

While sarcoma immunology has advanced with regard to basic, and even some applied topics, this disease has not been subject to more recent immunogenomics approaches. Thus, we assessed the immune receptor recombinations available from the cancer genome atlas (TCGA) sarcoma database via tumor sample exome and RNASeq files. Results indicated that recovery of T-cell receptor-alpha recombination reads (TRA) correlated with a better survival rate, with the expression of T-cell biomarkers, and with tumor sample apoptosis signatures consistent with the longer patient survival times. Furthermore, samples representing TRA complementarity determining region-3 (CDR3) net charge per residue (NCPR) based complementarity with the corresponding sarcoma mutanome had a better survival rate, and more granzyme expression, than samples lacking such complementarity. By specifically using RNASeq-recovered TRA CDR3s and related NCPR assessments, three genes, TP53, ATRX, and RB1, were identified as being key components of the mutanome-based complementarity. Thus, these genes may represent key immune system targets for soft tissue sarcomas. Also, several key results from above were reproduced with a pediatric osteosarcoma dataset, work that led to identification of MUC6 mutations as potentially linked to a strong immune response. In sum, TRA CDR3s are likely to be important prognostic indicators, and possibly a beginning tool for immunotherapy development strategies, for adult and pediatric sarcomas.

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.

High Level Immunoglobulin Gene Expression and Reduced Intra-tumoral Bacteria Associated With the Transition from Initial to Progressive Diffuse Intrinsic Pontine Glioma.

BACKGROUND/AIM: In the past decade, diffuse intrinsic pontine glioma (DIPG), the most common childhood brainstem glioma, has benefitted from an increase in tissue-based research because of improved biopsy collection techniques. However, the adaptive immune receptor (IR) features represented by tumor material and tumor infiltrating lymphocytes have remained poorly understood. MATERIALS AND METHODS: Herein, we characterized the adaptive immune parameters of DIPG through the recovery of IR recombination reads from RNAseq files representing initial and progressive DIPG samples. RESULTS: An elevated level of immunoglobulin gene expression in the progressive DIPG sample files and a reduced number of bacterial sequencing read recoveries in comparison to RNAseq files representing the initial form of DIPG, was found. Furthermore, the RNAseq files representing both initial and progressive DIPG samples had significant numbers of reads representing Cutibacterium acnes, a bacterium previously linked to prostate cancer development. Results also indicated an opportunity to distinguish overall survival probabilities based on IGL complementarity determining region-3 amino acid sequence physicochemical parameters. CONCLUSION: Genomics analyses allow for a better understanding of adaptive IR features and bacterial infections in the DIPG setting.

Neuroblastoma and Glioblastoma Cases With Amplified Oncogenes Have Reduced Numbers of Tumor-Resident Adaptive Immune Receptor Recombinations.

PURPOSE: In certain cancers, oncogene amplification is correlated with an immunologically cold or noninflamed, tumor immune microenvironment (TIME) and a worse prognosis, for example, in the case of MYCN-amplified neuroblastoma (NBL). However, for other cancer types, the relationship between oncogene amplification and immune response is more complicated or unresolved. One such cancer is glioblastoma multiforme (GBM), in which the epidermal growth factor receptor (EGFR) oncogene is commonly amplified. Unlike MYCN-amplified NBL, EGFR-amplified GBM has not been shown to correlate with a distinct survival probability. METHODS: Given this contrasting state for NBL and GBM, we sought to apply a genomics approach to evaluating the immune response for cases with gene amplification. RESULTS: Our results confirmed and added further specificity to the cold TIME of MYCN-amplified NBL. Moreover, we demonstrated a novel state of immunologically cold EGFR-amplified GBM tumors. CONCLUSION: This approach to using copy number variation and immune receptor recombination read recovery levels to assess gene amplification and TIME, respectively, may be particularly efficient for the rapid evaluation of many other cancer types.

Blood-based T cell receptor anti-viral CDR3s are associated with worse overall survival for neuroblastoma.

With the advent of large collections of adaptive immune receptor recombination reads representing cancer, there is the opportunity to further investigate the adaptive immune response to viruses in the cancer setting. This is a particularly important goal due to longstanding but still not well-resolved questions about viral etiologies in cancer and viral infections as comorbidities. In this report, we assessed the T cell receptor complementarity determining region-3 (CDR3) amino acid (AA) sequences, for blood-sourced TCRs from neuroblastoma (NBL) cases, for exact AA sequence matches to previously identified anti-viral TCR CDR3 AA sequences. Results indicated the presence of anti-viral TCR CDR3 AA sequences in the NBL blood samples highly significantly correlated with worse overall survival. Furthermore, the TCR CDR3 AA sequences demonstrating chemical complementarity to many cytomegalovirus antigens represented cases with a worse outcome, including cases where such CDR3s were obtained from tumor samples. Overall, these results indicate a significant need for, and provide a novel strategy for assessing viral infection complications in NBL patients.

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.

A Bioinformatics Assessment Indicating Better Outcomes With Breast Cancer Resident, Immunoglobulin CDR3-MMP2 Binding.

BACKGROUND/AIM: The recombination of V, D, and J immunoglobulin (IG) gene segments leads to many variations in the amino acids (AAs) encoded at that site, the complementarity determining region-3 (CDR3). Thus, cancer patients may have varying degrees of CDR3 AA binding specificity for cancer proteases, for example, matrix metalloproteinase 2 (MMP2). MMP2 in breast cancer has been found to contribute to metastasis and is used as a marker for tumor staging. Thus, this report evaluated the tumor resident, patient specific IG CDR3 binding affinities to cancer proteases to test the hypothesis that greater binding affinities would be associated with a better outcome. MATERIALS AND METHODS: Using two independent bioinformatics tools, we evaluated the IG CDR3-MMP2 binding affinities throughout the cancer genome atlas breast cancer (TCGA-BRCA) dataset. RESULTS: Results indicated that the better the CDR3-MMP2 binding, the better the survival probability. An analogous evaluation for four other proteases, including calpain-1 and thermolysin, displayed no such associations with survival probabilities. CONCLUSION: This study is consistent with the possibility that patient IG-cancer protease interactions could impact outcomes and raises the question of whether therapeutic antibody targeting of MMP2 would reduce breast cancer mediated tissue destruction and breast cancer mortality rates.

Exploiting adaptive immune receptor recombination read recoveries from exome files to identify subsets of ALL and to establish TCR features that correlate with better outcomes.

INTRODUCTION: The recovery of adaptive immune receptor (IR) recombination reads from tumour-derived genomics files has advanced the understanding of the immune system's interaction with cancer. This approach has been largely limited to solid tumours, where genomics file preparation allows for the recovery of adaptive IR reads corresponding to the T-cells and B-cells found in the solid tumour microenvironment. In this study, we sought to determine whether IR recombination reads from liquid tumour genomics files could also be informative. METHODS: We recovered the adaptive IR recombination reads from acute lymphoblastic leukaemia (ALL) normal and pathological genomics files of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET)-ALL, phase 2 project. RESULTS: In the bone marrow setting, results indicated that there was little or no B-cell response to ALL. However, results did show survival distinctions for B-cell ALL, in cases with specific T-cell complementarity determining region-3 chemical features, potentially reflecting specificity of the adaptive T-cell response against ALL. Furthermore, we found that the B-cell form of ALL, as well as what is likely TRD clonotypic, T-cell ALL, could likely be diagnosed via the recovery of B-cell receptor and TRD recombination reads, respectively, from pathological bone marrow exome files. CONCLUSIONS: Recovery of IR recombination reads from ALL exomes could aid in sub-type diagnoses and prognoses.

Immunogenomics Parameters for Patient Stratification in Alzheimer's Disease.

BACKGROUND: Despite the fact that only modest adaptive immune system related approaches to treating Alzheimer's disease (AD) are available, an immunogenomics approach to the study of AD has not yet substantially advanced. OBJECTIVE: Thus, we sought to better understand adaptive immune receptor chemical features in the AD setting. METHODS: We characterized T-cell receptor alpha (TRA) complementarity determining region-3 (CDR3) physicochemical features and identified TRA CDR3 homology groups, represented by TRA recombination reads extracted from 2,665 AD-related, blood- and brain-derived exome files. RESULTS: We found that a higher isoelectric value for the brain TRA CDR3s was associated with a higher (clinically worse) Braak stage and that a number of TRA CDR3 chemical homology groups, in particular representing bloodborne TRA CDR3s, were associated with higher or lower Braak stages. Lastly, greater chemical complementarity of both blood- and brain-derived TRA CDR3s and tau, based on a recently described CDR3-candidate antigen chemical complementarity scoring process (https://adaptivematch.com), was associated with higher Braak stages. CONCLUSION: Overall, the data reported here raise the questions of (a) whether progression of AD is facilitated by the adaptive immune response to tau; and (b) whether assessment of such an anti-tau immune response could potentially serve as a basis for adaptive immune receptor related, AD risk stratification?

Immune receptor CDR3 chemical features that preserve sequence information are highly efficient in reflecting survival distinctions: A pan-cancer analysis.

Adaptive immune receptor (IR) chemical features have been used as signatures of an immune response for numerous medical conditions, raising the question of whether certain approaches to assessing the IR chemical features are more robust than others? In the cancer setting, a very large dataset of IR complementarity determining region-3 (CDR3) amino acid (AA) sequences has become available via the mining of cancer specimen and blood genomics files for IR recombination reads. The IR CDR3 AA sequences have been evaluated for chemical features, and survival rates have been correlated with distinct chemical features. Two common approaches have been i) to assign a single value to the CDR3, representing a chemical attribute, such as aromaticity; or ii) to reduce the actual CDR3 AA sequence to a chemical sequence motif, which merges similar CDR3 chemistries represented by distinct AA sequences but preserves potential functional aspects of the order of the AAs in the sequence. While a controlled comparison of the two approaches is not possible, the application of the two approaches to the same clinical datasets offers the opportunity to appreciate a trend with regard to the overall potential in distinguishing survival probabilities. We demonstrate that application of the chemical sequence motif approach is more likely to identify survival distinctions within cancer datasets, for both tumor specimen and blood sourced, adaptive IR CDR3 AA sequences.

Tumor Resident, B-Cell Receptor Chemical Characteristics Associated with Better Overall Survival for Neuroblastoma.

Pediatric neuroblastoma (NBL) is one of the most common pediatric cancers, and it can often be aggressive. Genetic and demographic factors can correlate with the severity of NBL, but the variations in the B-cell receptors (BCRs) or immunoglobulin proteins present in the NBL tumors, and their relationships to survival, are not well understood. BCRs contain variations in their complementary determining region-3 (CDR3s) amino acid sequences, due to variable recombinations of the V- and J-gene segments. Accordingly, these variations in CDR3s may represent different antigen interactions and thereby different survival probabilities. Thus, we mined the TARGET project, NBL tumor RNAseq files for BCR recombination reads. Evaluations of the physicochemical properties of IGK, IGL, and IGH CDR3s from these tumors pointed to properties of IGK and IGL in particular as associated with survival distinctions, based on several independent bioinformatics approaches, including a novel homology grouping approach facilitated by a recently developed web tool, adaptivematch.com. In conclusion, tumor resident BCR chemical features are likely useful for better risk stratification and for guiding therapy, and the availability of a user-friendly web tool will likely facilitate using BCR chemical features to meet those goals.

Systemic Adaptive Immune Parameters Associated with Neuroblastoma Outcomes: the Significance of Gamma-Delta T Cells.

We mined a set of neuroblastoma (NBL) exomes for immune receptor recombinations representing all seven, human adaptive immune receptor genes, using a very high standard for identifications of V- and J-gene segments in sequencing reads. Results indicated an unusually large number of TRD recombination reads in the NBL samples, possibly related to the younger age of the patients. In general, recovery of immune receptor (IR) recombination reads representing any of the immune receptors, from either blood or tumor samples, was associated with a lower overall survival rate, consistent with an emerging literature indicating that systemic immunology parameters can be informative for cancer evaluations. Despite the overall negative association of IR recombination frequencies and outcomes, survival rate distinctions could still be established as associated with certain chemical features of IR complementarity determining region-3 (CDR3) amino acid sequences, thereby likely revealing a distinction between the negative impacts of a general adaptive immune response versus the positive aspects of specific CDR3 chemical interaction potentials. These data underscore the relevance of gamma-delta T cells in the development of cancer in younger patients. And for the first time, these data allow a distinction within an NBL cohort with active disease, between two contrasting systemic immune states: (i) general and likely harmful adaptive immunity development versus (ii) a likely positive, adaptive immune response with particular antigenic specificities.

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

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.

Specific HLA alleles, paired with TCR V- and J-gene segment usage, link to distinct multiple myeloma survival rates.

Multiple myeloma (MM) immunogenomics studies related to T-cell characterizations and involving large patient sets have been lacking, particularly in comparison to solid tumor types. Thus, we evaluated (i) HLA alleles, and (ii) T-Cell Receptor (TCR) V- and J-gene segment, HLA allele combinations, based on TCR recombinations in blood samples, for their potential associations with overall survival distinctions among an MM cohort. Two HLA alleles, and seven TCR V- or J-gene segment, HLA allele combinations were found to be associated with distinct overall survival rates. For examples, HLA-C*08:02, and the TRAV19, HLA-C*07:01 combination, were found to be associated with negative outcomes. In addition, anti-cytomegalovirus immune receptor sequences, from blood samples, were found to be associated with a positive outcome (p = 0.012, n = 278). These data, and other related immunogenomics data, indicate a potential opportunity to use personal immunogenetics parameters as guides to prognosis and therapies.

TRBV and TRBJ usage, when paired with specific HLA alleles, associates with distinct head and neck cancer survival rates.

Common or dominant, T-cell receptor (TCR), V and J usage, in combination with particular human leukocyte antigen (HLA) alleles, has been associated with differing outcomes in viral infections, autoimmunity, and more recently, in cancer. Cervical cancer in particular represents the most dramatic series of distinctions of outcomes associated with differing combinations of dominant V or J usage and HLA alleles, possibly because of the strong association of cervical cancer with human papilloma virus (HPV), in turn leading to a likely molecular consistency in the mechanism of HPV antigen presentation. Thus, we considered assessing TRB V and J usage, HLA allele combinations, for their associations with survival rates and related data, in the cancer genome atlas head and neck cancer dataset. We obtained the TRB VDJ recombination reads from both the blood and tumor exome files and determined the V and J identities. We then established case ID (patient) subsets of V or J usage, HLA alleles, and determined, for example, that the TRBJ2-7, HLA-B*40:01 combination was associated with a better disease free survival rate than were either the TRBJ1-3, HLA-DPB1*03:01 or the TRBJ2-1, HLA-DPB1*02:01 combinations. Furthermore, these analyses led to the conclusion that TRBJ1-5 usage, and the HLA-C*08:02 and HLA-DRB1*03:01 alleles, had independent associations with distinct overall survival rates. In sum, the results suggest that dominant V or J usage, HLA allele combinations, and in certain cases, dominant V or J usage independently of HLA, could be useful in prognosis and in guiding immunotherapies.