The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

Identification of HLA-DPA1*01:03:01:57 and HLA-DPA1*02:01:01:29 from a case-control study of atopic dermatitis.

DPA1*01:03:01:57 and DPA1*02:01:01:29 differ by a single nucleotide from their closest references, DPA1*01:03:01:02 and DPA1*02:01:01:06.

Human abdominal aortic aneurysm (AAA): Evidence for an autoimmune antigen-driven disease.

Abdominal aortic aneurism (AAA) is a complex immunological disease with a strong genetic component, and one of the ten leading causes of death of individuals 55-74 years old worldwide. Strong evidence has been accumulated suggesting that AAA is an autoimmune specific antigen-driven disease. Mononuclear cells infiltrating AAA lesions comprised of T and B lymphocytes and other cells expressing early-, intermediate- and late-activation antigens, and the presence of antigen-presenting cells have been documented, demonstrating an ongoing immune response. The three components of the trimolecular complex, T-cell receptor (TCR)/peptide (antigen)/HLA have been identified in AAA, and specifically: (i) clonal expansions of T-cell clones in AAA lesions; (ii) the association of AAA with particular HLA Class I and Class II; and (iii) self or nonself putative AAA-associated antigens. IgG autoantibodies recognizing proteins present in normal aortic tissue have been reported in patients with AAA. Molecular mimicry, defined as the sharing of antigenic epitopes between microorganisms (bacteria, viruses) and self antigens, maybe is responsible for T-cell responses and antibody production in AAA. Also, the frequency and the suppressor activity of CD4+ CD25+ FOXP3+ Tregs and the expression of FOXP3 transcripts and protein have been reported to be significantly impaired in AAA patients vs normal donors.

Unrelated donor α/ß T cell- and B cell-depleted HSCT for the treatment of pediatric acute leukemia.

Unrelated donor (URD) hematopoietic stem cell transplant (HSCT) is associated with an increased risk of severe graft-versus-host disease (GVHD). TCRαß/CD19 depletion may reduce this risk, whereas maintaining graft-versus-leukemia. Outcome data with TCRαß/CD19 depletion generally describe haploidentical donors, with relatively few URDs. We hypothesized that TCRαß/CD19-depletion would attenuate the risks of GVHD and relapse for URD HSCT. Sixty pediatric and young adult (YA) patients with hematologic malignancies who lacked a matched-related donor were enrolled at 2 large pediatric transplantation centers between October 2014 and September 2019. All patients with acute leukemia had minimal residual disease testing, and DP typing was available for 77%. All patients received myeloablative total body irradiation- or busulfan-based conditioning with no posttransplant immune suppression. Engraftment occurred in 98%. Four-year overall survival was 69% (95% confidence interval [CI], 52%-81%), and leukemia-free survival was 64% (95% CI, 48%-76%), with no difference between lymphoid and myeloid malignancies (P = .6297 and P = .5441, respectively). One patient (1.7%) experienced primary graft failure. Relapse occurred in 11 patients (3-year cumulative incidence, 21%; 95% CI, 11-34), and 8 patients (cumulative incidence, 15%; 95% CI, 6.7-26) experienced nonrelapse mortality. Grade III to IV acute GVHD was seen in 8 patients (13%), and 14 patients (26%) developed chronic GVHD, of which 6 (11%) had extensive disease. Nonpermissive DP mismatch was associated with higher likelihood of acute GVHD (odds ratio, 16.50; 95% CI, 1.67-163.42; P = .0166) but not with the development of chronic GVHD. URD TCRαß/CD19-depleted peripheral HSCT is a safe and effective approach to transplantation for children/YAs with leukemia. This trial was registered at www.clinicaltrials.gov as #NCT02323867.

Association of KIR Genes and MHC Class I Ligands with Atopic Dermatitis.

Atopic dermatitis (AD) is a chronic illness that is associated with immune dysregulation. NK cell function has previously been associated with AD. NK cells directly interact with polymorphic HLA class I ligand variants using killer cell Ig-like receptors (KIRs). The purpose of this study was to identify potential associations between NK cell function and AD by evaluating variation in the presence of KIR genes as well as KIR gene interactions with the appropriate HLA class I KIR-specific ligands. Human DNA from the genetics of AD case-control study was used to genotype HLA class I KIR-specific ligands and the presence of KIR genes. In the full cohort, an increased risk of AD was noted for KIR2DL5 (1.51 [1.13, 2.01]), KIR2DS5 (1.72 [1.26, 2.34]), and KIR2DS1 (1.41 [1.04, 1.91]). Individuals with KIR2DS5 or KIR2DS1 and the HLA-C*C2 epitope were at an increased risk of AD (1.74 [1.21, 2.51] and 1.48 [1.04, 2.12], respectively). The HLA-B*-21T (TT) leader sequence increased the risk of AD across ethnicity. African Americans with KIR2DL2, KIR2DS1, KIR2DL5, and KIR2DS5 are more likely to have AD, and the risk increased for KIR2DS1 and KIR2DS5 in the presence of appropriate HLA-C C2 epitope. The risk of AD also increased for individuals with the HLA-B*-21T leader sequence. Future studies should focus on KIR gene allelic variation as well as consider cell-based measurements of KIR and the associated HLA class I epitopes.

Analysis of single nucleotide polymorphisms in chronic beryllium disease.

Sarcoidosis and chronic beryllium disease (CBD) are phenocopies, however the latter one has a clear trigger factor that is beryllium exposure. This study analyses single nucleotide polymorphisms (SNPs) in a large cohort for beryllium-exposed persons. SNPs were chosen for their relevance in sarcoidosis. Even though one of largest cohorts of beryllium-exposed persons was analysed, no statistically relevant association between any SNP and CBD could be verified. Notably, some SNPs exhibit inverse OR for beryllium sensitization and CBD with nominally statistical significance, which allows hypothesizing about pathophysiological role of genes for the disease triggering and development.

HLA Class I Polymorphisms Influencing Both Peptide Binding and KIR Interactions Are Associated with Remission among Children with Atopic Dermatitis: A Longitudinal Study.

Atopic dermatitis (AD) is a disease of immune dysregulation and skin barrier dysfunction with a relapsing, remitting course and has been associated with several different genetic risk variants. HLA represent a highly variable set of genes that code for cell surface protein molecules involved in the Ag-specific immune response, including the regulation or functioning of T cells, NK cells, and APCs. The purpose of this study was to evaluate associations between HLA class I polymorphisms and the progression of AD over time. We evaluated the associations of AD symptoms and HLA class I polymorphisms based on high-resolution two-field typing in a longitudinal cohort of children with AD (up to 10 y of follow-up). Seven hundred and ninety-two children were evaluated every 6 mo, resulting in 12,752 AD evaluations. Using generalized estimating equations and corrected p values, B*44:02 was found to be associated with AD remission (1.83 [1.35, 2.47]; p = 0.0015). The HLA-B residues at position 116 (d-aspartate) and 80 (T-threonine) were associated with remission (1.42 [1.13, 1.76], p = 0.003; corrected p = 0.028) and (1.45 [1.17, 1.80], p = 0.0008; corrected p = 0.0024), respectively. B80T is a killer-cell Ig-like receptor (KIR) site. Our findings reveal that two axes of immune response (T cell and NK cell) may influence disease progression. Identifying binding pocket changes in addition to other factors (e.g., allergens) that increase the risk or severity of AD can improve our understanding of the immunologic mechanisms associated with AD and may lead to personalized therapies for improving patient care.

Spatial constrains and information content of sub-genomic regions of the human genome.

Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes.

BACKGROUND: HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk. METHODS: One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed. RESULTS: A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1. CONCLUSIONS: The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

Genomic characterization of MICA gene using multiple next generation sequencing platforms: A validation study.

We have developed a protocol regarding the genomic characterization of the MICA gene by next generation sequencing (NGS). The amplicon includes the full length of the gene and is about 13 kb. A total of 156 samples were included in the study. Ninety-seven of these samples were previously characterized at MICA by legacy methods (Sanger or sequence specific oligonucleotide) and were used to evaluate the accuracy, precision, specificity, and sensitivity of the assay. An additional 59 DNA samples of unknown ethnicity volunteers from the United States were only genotyped by NGS. Samples were chosen to contain a diverse set of alleles. Our NGS approach included a first round of sequencing on the Illumina MiSeq platform and a second round of sequencing on the MinION platform by Oxford Nanopore Technology (ONT), on selected samples for the purpose of either characterizing new alleles or setting phase among multiple polymorphisms to resolve ambiguities or generate complete sequence for alleles that were only partially reported in the IMGT/HLA database. Complete consensus sequences were generated for every allele sequenced with ONT, extending from the 5' untranslated region (UTR) to the 3' UTR of the MICA gene. Thirty-two MICA sequences were submitted to the IMGT/HLA database including either new alleles or filling up the gaps (exonic, intronic and/or UTRs) of already reported alleles. Some of the challenges associated with the characterization of these samples are discussed.

Utilizing nanopore sequencing technology for the rapid and comprehensive characterization of eleven HLA loci; addressing the need for deceased donor expedited HLA typing.

The comprehensive characterization of human leukocyte antigen (HLA) genomic sequences remains a challenging problem. Despite the significant advantages of next-generation sequencing (NGS) in the field of Immunogenetics, there has yet to be a single solution for unambiguous, accurate, simple, cost-effective, and timely genotyping necessary for all clinical applications. This report demonstrates the benefits of nanopore sequencing introduced by Oxford Nanopore Technologies (ONT) for HLA genotyping. Samples (n = 120) previously characterized at high-resolution three-field (HR-3F) for 11 loci were assessed using ONT sequencing paired to a single-plex PCR protocol (Holotype) and to two multiplex protocols OmniType (Omixon) and NGSgo®-MX6-1 (GenDx). The results demonstrate the potential of nanopore sequencing for delivering accurate HR-3F typing with a simple, rapid, and cost-effective protocol. The protocol is applicable to time-sensitive applications, such as deceased donor typings, enabling better assessments of compatibility and epitope analysis. The technology also allows significantly shorter turnaround time for multiple samples at a lower cost. Overall, the nanopore technology appears to offer a significant advancement over current next-generation sequencing platforms as a single solution for all HLA genotyping needs.

Lifetime and Stability of Silicon Nitride Nanopores and Nanopore Arrays for Ionic Measurements.

Nanopores are promising for many applications including DNA sequencing and molecular filtration. Solid-state nanopores are preferable over their biological counterparts for applications requiring durability and operation under a wider range of external parameters, yet few studies have focused on optimizing their robustness. We report the lifetime and durability of pores and porous arrays in 10 to 100 nm-thick, low-stress silicon nitride (SiNx) membranes. Pores are fabricated using a transmission electron microscope (TEM) and/or electron beam lithography (EBL) and reactive ion etching (RIE), with diameters from 2 to 80 nm. We store them in various electrolyte solutions (KCl, LiCl, MgCl2) and record open pore conductance over months to quantify pore stability. Pore diameters increase with time, and diameter etch rate increases with electrolyte concentration from Δd/Δt ∼ 0.2 to ∼ 3 nm/day for 0.01 to 3 M KCl, respectively. TEM confirms the range of diameter etch rates from ionic measurements. Using electron energy loss spectroscopy (EELS), we observe a N-deficient region around the edges of TEM-drilled pores. Pore expansion is caused by etching of the Si/SiO2 pore walls, which resembles the dissolution of silicon found in minerals such as silica (SiO2) in salty ocean water. The etching process occurs where the membrane was exposed to the electron beam and can result in pore formation. However, coating pores with a conformal 1 nm-thick hafnium oxide layer prevents expansion in 1 M KCl, in stark contrast to bare SiNx pores (∼ 1.7 nm/day). EELS data reveal the atomic composition of bare and HfO2-coated pores.

Partially CD3+-Depleted Unrelated and Haploidentical Donor Peripheral Stem Cell Transplantation Has Favorable Graft-versus-Host Disease and Survival Rates in Pediatric Hematologic Malignancy.

Most children who may benefit from stem cell transplantation lack a matched related donor. Alternative donor transplantations with an unrelated donor (URD) or a partially matched related donor (PMRD) carry an increased risk of graft-versus-host-disease (GVHD) and mortality compared with matched related donor transplantations. We hypothesized that a strategy of partial CD3+/CD19+ depletion for URD or PMRD peripheral stem cell transplantation (PSCT) would attenuate the risks of GVHD and mortality. We enrolled 84 pediatric patients with hematologic malignancies at the Children's Hospital of Philadelphia and the Children's Hospital of Wisconsin between April 2005 and February 2015. Two patients (2.4%) experienced primary graft failure. Relapse occurred in 23 patients (27.4%; cumulative incidence 26.3%), and 17 patients (20.2%) experienced nonrelapse mortality (NRM). Grade III-IV acute GVHD was observed in 18 patients (21.4%), and chronic GVHD was observed and graded as limited in 24 patients (35.3%) and extensive in 8 (11.7%). Three-year overall survival (OS) was 61.8% (95% confidence interval [CI], 50.2% to 71.4%) and event-free survival (EFS) was 52.0% (95% CI, 40.3% to 62.4%). Age ≥15 years was associated with decreased OS (P= .05) and EFS (P= .05). Relapse was more common in children in second complete remission (P = .03). Partially CD3+-depleted alternative donor PSCT NRM, OS, and EFS compare favorably with previously published studies of T cell-replete PSCT. Historically, T cell-replete PSCT has been associated with a higher incidence of extensive chronic GVHD compared with limited chronic GVHD, which may explain the comparatively low relapse and NRM rates in our study cohort despite similar overall rates of chronic GVHD. Partial T cell depletion may expand donor options for children with malignant transplantation indications lacking a matched related donor by mitigating, but not eliminating, chronic GVHD.

Clonally expanded alpha-chain T-cell receptor (TCR) transcripts are present in aneurysmal lesions of patients with Abdominal Aortic Aneurysm (AAA).

Abdominal aortic aneurysm (AAA) is a life-threatening immunological disease responsible for 1 to 2% of all deaths in 65 year old or older individuals. Although mononuclear cell infiltrates have been demonstrated in AAA lesions and autoimmunity may be responsible for the initiation and account for the propagation of the disease, the information available about the pathogenesis of AAA is limited. To examine whether AAA lesions from patients with AAA contain clonally expanded α-chain TCR transcripts, we amplified by the non-palindromic adaptor-PCR (NPA-PCR)/Vα-specific PCR and/or the Vα-specific PCR these α-chain TCR transcripts. The amplified transcripts were cloned and sequenced. Substantial proportions of identical α-chain TCR transcripts were identified in AAA lesions of 4 of 5 patients, demonstrating that clonally expanded T cells are present in these AAA lesions. These results were statistically significant by the bimodal distribution. Three of 5 of these patients were typed by DNA-based HLA-typing and all three expressed DRB1 alleles containing the DRßGln70 amino acid residue that has been demonstrated to be associated with AAA. All three patients exhibited clonally expanded T cells in AAA lesions. Four of the 5 patients with AAA who exhibited clonal expansions of α-chain TCR transcripts, also exhibited clonal expansions of ß-chain TCR transcripts in AAA lesions, as we have demonstrated previously (J Immunol 192:4897, 2014). αß TCR-expressing T cells infiltrating AAA lesions contain T-cell clones which have undergone proliferation and clonal expansion in vivo in response to as yet unidentified specific antigens that may be self or nonself. These results provide additional evidence supporting the hypothesis that AAA is a specific antigen-driven T-cell autoimmune disease.

Resolving MiSeq-Generated Ambiguities in HLA-DPB1 Typing by Using the Oxford Nanopore Technology.

The technical limitations of current next-generation sequencing technologies, combined with an ever-increasing number of human leukocyte antigen (HLA) alleles, form the basis for the additional ambiguities encountered at an increasing rate in clinical practice. HLA-DPB1 characterization, particularly, generates a significant percentage of ambiguities (25.5%), posing a challenge for accurate and unambiguous HLA-DPB1 genotyping. Phasing of exonic heterozygous positions between exon 2 and all other downstream exons has been the major cause of ambiguities. In this study, the Oxford Nanopore MinION, a third-generation sequencing technology, was used to resolve the phasing. The accurate MiSeq sequencing data, combined with the long reads obtained from the MinION platform, allow for the resolution of the tested ambiguities.

Personalized cancer vaccine effectively mobilizes antitumor T cell immunity in ovarian cancer.

We conducted a pilot clinical trial testing a personalized vaccine generated by autologous dendritic cells (DCs) pulsed with oxidized autologous whole-tumor cell lysate (OCDC), which was injected intranodally in platinum-treated, immunotherapy-naïve, recurrent ovarian cancer patients. OCDC was administered alone (cohort 1, n = 5), in combination with bevacizumab (cohort 2, n = 10), or bevacizumab plus low-dose intravenous cyclophosphamide (cohort 3, n = 10) until disease progression or vaccine exhaustion. A total of 392 vaccine doses were administered without serious adverse events. Vaccination induced T cell responses to autologous tumor antigen, which were associated with significantly prolonged survival. Vaccination also amplified T cell responses against mutated neoepitopes derived from nonsynonymous somatic tumor mutations, and this included priming of T cells against previously unrecognized neoepitopes, as well as novel T cell clones of markedly higher avidity against previously recognized neoepitopes. We conclude that the use of oxidized whole-tumor lysate DC vaccine is safe and effective in eliciting a broad antitumor immunity, including private neoantigens, and warrants further clinical testing.

Sensitive and frequent identification of high avidity neo-epitope specific CD8 + T cells in immunotherapy-naive ovarian cancer.

Immunotherapy directed against private tumor neo-antigens derived from non-synonymous somatic mutations is a promising strategy of personalized cancer immunotherapy. However, feasibility in low mutational load tumor types remains unknown. Comprehensive and deep analysis of circulating and tumor-infiltrating lymphocytes (TILs) for neo-epitope specific CD8+ T cells has allowed prompt identification of oligoclonal and polyfunctional such cells from most immunotherapy-naive patients with advanced epithelial ovarian cancer studied. Neo-epitope recognition is discordant between circulating T cells and TILs, and is more likely to be found among TILs, which display higher functional avidity and unique TCRs with higher predicted affinity than their blood counterparts. Our results imply that identification of neo-epitope specific CD8+ T cells is achievable even in tumors with relatively low number of somatic mutations, and neo-epitope validation in TILs extends opportunities for mutanome-based personalized immunotherapies to such tumors.