High-throughput quantitative assessments of the chemical complementarity of celiac disease related IGH CDR3s and a gliadin epitope.

The long-term value of efficient antigen discovery includes gaining insights into the variety of potential cancer neoantigens, effective vaccines lacking adverse effects, and adaptive immune receptor (IR) targets for blocking adaptive IR-antigen interactions in autoimmunity. While the preceding goals have been partially addressed via big data approaches to HLA-epitope binding, there has been little such progress in the big data setting for adaptive IR-epitope binding. This delay in progress for the latter is likely due to, among other things, the much more complicated adaptive IR repertoire in an individual compared to individual HLA alleles. Thus, results described here represent the application of an algorithm for efficient assessment of IGH CDR3-gliadin epitope interactions, with a focus on epitopes known to be associated with an immune response in celiac disease. The hydrophobic, chemical complementarity between celiac case IGH CDR3s and known celiac epitopes was found to be greater in comparison to the hydrophobic, chemical complementarity between the same celiac case IGH CDR3s and a series of control epitopes. Thus, the approaches indicated here likely offer guidance for the development of conveniently applied algorithms for antigen verification and discovery.

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.

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.

Exploiting big data survival information to unify risk-stratification related, adaptive immune receptor parameters for multiple myeloma.

With the improvement of treatment options, multiple myeloma related life expectancy has been prolonged, but the disease remains largely incurable. Immunotherapy is a growing field that shows promise in advancements for treatment, and recent work has demonstrated an opportunity to use immune receptor, complementarity determining region-3 (CDR3)-candidate antigen chemical complementarity scores to identify survival distinctions among subgroups of patients. Here, we have applied the complementarity scoring algorithm to identify multiple myeloma related, CDR3-cancer testis antigen (CTA) relationships associated with survival distinctions. Furthermore, we have overlapped these immune receptor features with a previous study that showed a dramatic survival distinction based on T-cell receptor, V- and J-gene segment usage, HLA allele combinations, whereby 100% of the patients in certain combination groups had no mortality related to multiple myeloma, during the study period. This overlap evaluation was consistent with the idea that there are likely considerable constraints on productive TRB-antigen-HLA combinations but more flexibility, and unpredictability, for the TRA-antigen-HLA combinations. Also, the approaches in this reported indicated the potential importance of the CTA, IGSF11, as a multiple myeloma antigen, an antigen previously, independently considered as a vaccine candidate in other settings.

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?

Adaptive immune receptor features related to breast cancer tissue in Kenyan patients: high immunoglobulin gene expression and high levels of gamma-delta T-cells.

BACKGROUND: Characterization of the breast cancer (BC) immune response may provide information for a point of intervention, such as application of immunotherapeutic treatments. In this study, we sought to recover and characterize the adaptive immune receptor (IR) recombination reads from genomics files representing Kenyan patients, to better understand the immune response specifically related to those patients. METHODS: We used a previously applied algorithm and software to obtain productive IR recombination reads from cancer and adjacent normal tissue samples representing 22 Kenyan BC patients. RESULTS: From both the RNAseq and exome files, there were significantly more T-cell receptor (TCR) recombination reads recovered from tumor samples compared to marginal tissue samples. Also, the immunoglobulin (IG) genes were expressed at a much higher level than the TCR genes (p-value = 0.0183) in the tumor samples. And, the tumor IG CDR3s consistently represented more positively charged amino acid R-groups, in comparison to the marginal tissue, IG CDR3s. CONCLUSION: For Kenyan patients, a high level of IG expression, representing specific CDR3 chemistries, was associated with BC. These results lay the foundation for studies that could support specific immunotherapeutic interventions for Kenyan BC patients.

TRB CDR3 chemical complementarity with HBV epitopes correlates with increased hepatocellular carcinoma, disease-free survival.

The liver is a site of immune privilege, compared with the bladder and skin, for example. To study this attenuation of the immune response in the cancer setting, we compared quantities and features of adaptive immune receptor (IR) recombination reads obtained from hepatocellular carcinoma (HCC) and six other cancers. Of these cancers, HCC had the lowest numbers of IR recombination reads and was the only cancer with a greater number immunoglobulin rather than T-cell receptor recombination reads. To better understand the role of adaptive IRs obtained from the tumor microenvironment in shaping the outcome of HCC cases, we quantified the chemical complementarity between HCC tumor TRB and IGH complementarity determining region-3 (CDR3) amino acid (AA) sequences, and known hepatitis B virus (HBV) epitopes. High chemical complementarity between HCC-resident CDR3s and three HBV epitopes correlated with increased survival probabilities, for two sources of CDR3s representing different CDR3 recovery algorithms. These results suggest the potential of CDR3 AA sequences as biomarkers for HCC patient stratification and as guides for future development of therapeutics.

TCR gene segment usage and HLA alleles that are associated with cancer survival rates also represent racial disparities.

Understanding racial disparities in cancer outcomes continues to be a challenge, with likely many factors at play, including socioeconomic factors and genetic polymorphisms impacting basic cellular and molecular functions. Additionally, it is possible that specific combinations of environment and genetics have specific impacts. T-cell receptor (TCR) gene segment usage, HLA allele combinations have been associated with autoimmune and infectious disease courses, and more recently, TCR gene segment usage, HLA allele combinations have been associated with distinct survival outcomes in cancer as well. We examined several such, previously reported cancer-related TCR gene segment usage, HLA allele combinations for evidence of racial disparities, with regard to the prevalence of the combination in different racial groups. Results indicated that TCR gene segment usage, potentially reflecting environmental factors related to previous pathogen exposure, in combination with certain HLA alleles or independently, may represent a novel explanation for racial disparities in cancer outcomes. Overall, at this point, a genetic connection to racial disparities in cancer outcomes is detectable but remains modest, suggesting that other factors, such as socioeconomic factors, remain as important considerations.

Chemical features of blood-borne TRG CDR3s associated with an increased overall survival in breast cancer.

PURPOSE: Immunogenomics and earlier, pioneering studies, particularly by Whiteside and colleagues, have indicated a positive role for B-cells in breast cancer, as well as a positive role for gamma-delta T-cells. However, these studies have been completely limited to assessing breast cancer tumor tissue. METHODS AND RESULTS: Our analyses here has shown that blood-borne T-cell receptor gamma (TRG) chain sequences were associated with greater overall survival, of particular note due to the comparative longevity of primary breast cancer patients, whereby assessments of disease-free, but rarely overall survival parameters are possible. Additional immunogenomics approaches narrowed the overall survival correlations to specific, TRG complementarity determining region-3, amino acid (AA) sequence chemical features, independently of many common, confounding variables in the breast cancer setting, such as estrogen or progesterone receptor status. CONCLUSIONS: These results are discussed in the context of patient age and with regard to potential antigenic targets, based on the chemistry of the TRG CDR3 AA sequences associated with the higher survival rates.

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.

Potential MMP2-mediated availability of HLA binding, mutant ECM peptides reflects better melanoma survival rates and greater T-cell infiltrates.

Proteases in the cancer microenvironment have been studied for some time, with a general conclusion that such proteases facilitate the spread of cancer, although there is some controversy regarding that conclusion in later-stage cancer development. More recently, a very large collection of data regarding mutant amino acids in the potential substrates of cancer microenvironment proteases have become available. To better understand the potential impact of these mutant amino acids on protease function and cancer progression, we established a bioinformatics approach to assessing the impact of melanoma mutants, among a previously defined set of extracellular matrix (ECM) structural proteins, on the sensitivity of matrix metalloproteinase-2 (MMP2), extensively associated with melanoma. The results indicated that tumor samples with mutant amino acids adjacent to the ECM structural protein, MMP2 sites also represented a better survival rate and a larger proportion of mutant peptides with high HLA class I-binding affinities, particularly in comparison with melanoma samples with a reduced or absent T-cell infiltrate. Furthermore, even better HLA class I binders were identified among the samples representing the ECM structural protein mutants resistant to MMP2. Samples representing only MMP2-resistant mutants also represented a worse overall survival. Overall, this analysis suggested that MMP2 has the capacity of freeing mutant peptides that could facilitate an anti-tumor response and a better survival rate, and this analysis has the potential of resolving some of the controversy surrounding the role of cancer proteases in cancer progression.

Immune receptor recombinations from breast cancer exome files, independently and in combination with specific HLA alleles, correlate with better survival rates.

PURPOSE: Immune characterizations of cancers, including breast cancer, have led to information useful for prognoses and are considered to be important in the future of refining the use of immunotherapies, including immune checkpoint inhibitor therapies. In this study, we sought to extend these characterizations with genomics approaches, particularly with cost-effective employment of exome files. METHODS: By recovery of immune receptor recombination reads from the cancer genome atlas (TCGA) breast cancer dataset, we observed associations of these recombinations with T-cell and B-cell biomarkers and with distinct survival rates. RESULTS: Recovery of TRD or IGH recombination reads was associated with an improved disease-free survival (p = 0.047 and 0.045, respectively). Determination of the HLA types using the exome files allowed matching of T-cell receptor V- and J-gene segment usage with specific HLA alleles, in turn allowing a refinement of the association of immune receptor recombination read recoveries with survival. For example, the TRBV7, HLA-C*07:01 combination represented a significantly worse, disease-free outcome (p = 0.014) compared to all other breast cancer samples. By direct comparisons of distinct TRB gene segment usage, HLA allele combinations revealed breast cancer subgroups, within the entire TCGA breast cancer dataset with even more dramatic survival distinctions. CONCLUSIONS: In sum, the use of exome files for recovery of adaptive immune receptor recombination reads, and the simultaneous determination of HLA types, has the potential of advancing the use of immunogenomics for immune characterization of breast tumor samples.

An age-based, RNA expression paradigm for survival biomarker identification for pediatric neuroblastoma and acute lymphoblastic leukemia.

BACKGROUND: Pediatric cancer survival rates overall have been improving, but neuroblastoma (NBL) and acute lymphoblastic leukemia (ALL), two of the more prevalent pediatric cancers, remain particularly challenging. One issue not yet fully addressed is distinctions attributable to age of diagnosis. METHODS: In this report, we verified a survival difference based on diagnostic age for both pediatric NBL and pediatric ALL datasets, with younger patients surviving longer for both diseases. We identified several gene expression markers that correlated with age, along a continuum, and then used a series of age-independent survival metrics to filter these initial correlations. RESULTS: For pediatric NBL, we identified 2 genes that are expressed at a higher level in lower surviving patients with an older diagnostic age; and 4 genes that are expressed at a higher level in longer surviving patients with a younger diagnostic age. For pediatric ALL, we identified 3 genes expressed at a higher level in lower surviving patients with an older diagnostic age; and 17 genes expressed at a higher level in longer surviving patients with a younger diagnostic age. CONCLUSIONS: This process implicated pan-chromosome effects for chromosomes 11 and 17 in NBL; and for the X chromosome in ALL.

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.

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.

Cytoskeleton and ECM tumor mutant peptides: Increased protease sensitivities and potential consequences for the HLA class I mutant epitope reservoir.

Cytoskeleton and extracellular matrix-related proteins (CECMPs) represent the most common class of cancer mutants, owing to the large size of their coding regions and the randomness of mutagenesis. We used a bioinformatics approach to assess the impact of amino acid (AA) substitutions on the sensitivity of CECMPs to proteases relevant to melanoma and on the binding affinities for HLA class I. CECMP peptides with AA substitutions overwhelmingly reflect increased sensitivity to proteases implicated in melanoma development (MME, CTSS, MMP2, CTSD, CTSL) in comparison to the wild-type peptide sequences. Furthermore, peptides with AA substitutions representing increased peptide protease sensitivity also represented relatively high binding affinities for HLA class I allelic variants. These analyses raise the question of whether increased protease sensitivity, of mutant cancer peptides represents a significant increase in the availability of cancer mutant, HLA class I epitopes and a hitherto unappreciated aspect of cancer cell immunogenicity, particularly in the case of melanoma?

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.