Mutations in mitochondrial (mt) DNA accumulate with age and can result in the generation of neopeptides. Immune surveillance of such neopeptides may allow suboptimal mitochondria to be eliminated, thereby avoiding mt-related diseases, but may also contribute to autoimmunity in susceptible individuals. To date, the direct recognition of neo-mtpeptides by the adaptive immune system has not been demonstrated. In this study we used bioinformatics approaches to predict MHC binding of neopeptides identified from known deletions in mtDNA. Six such peptides were confirmed experimentally to bind to HLA-A*02. Pre-existing human CD4+ and CD8+ T cells from healthy donors were shown to recognize and respond to these neopeptides. One remarkably promiscuous immunodominant peptide (P9) could be presented by diverse MHC molecules to CD4+ and/or CD8+ T cells from 75% of the healthy donors tested. The common soil microbe, Bacillus pumilus, encodes a 9-mer that differs by one amino acid from P9. Similarly, the ATP synthase F0 subunit 6 from normal human mitochondria encodes a 9-mer with a single amino acid difference from P9 with 89% homology to P9. T cells expanded from human PBMCs using the B. pumilus or self-mt peptide bound to P9/HLA-A2 tetramers, arguing for cross-reactivity between T cells with specificity for self and foreign homologs of the altered mt peptide. These findings provide proof of principal that the immune system can recognize peptides arising from spontaneous somatic mutations and that such responses might be primed by foreign peptides and/or be cross-reactive with self.

The hallmark of the immune system is its ability to distinguish self from nonself. The lack of response to self and the response to nonself is paramount for maintenance of human health, allowing the elimination of infection while minimizing harmful autoimmune reactions. The knowledge of mutation frequencies, coupled with the knowledge that mutated self-proteins can be maintained and presented, raises the question as to how the immune system manages its mutated self. Mitochondria are organelles with which this question may be answered. Mitochondria exist in many copies in a cell, they have multiple genomes, and these genomes can be (and are) somatically mutated to varying degrees and maintained within the cell (1, 2). Importantly, mitochondrial (mt) proteins can be presented to the immune system (3, 4).

The evolutionary origin of eukaryotic mitochondria can be traced to α-proteobacteria’s merger with the ancestors of present-day eukaryotic cells (5, 6). Mitochondrial peptides present a challenge to the immune system in regard to self/nonself discrimination because of this bacterial genetic heritage. Notably, peptides derived from native mt proteins can contribute to autoimmune diseases as evidenced in the disease primary biliary cirrhosis, where mt peptides are presented as autoantigens (7, 8).

mtDNA has a high mutation rate, ∼10- to 20-fold the rate of nuclear DNA (2, 911). This high mutation rate coupled with the multiple mt genomes per cell results in cells with mitochondria that have wild-type (WT) and mutated mtDNAs (point mutations, deletions, or insertions), a condition known as heteroplasmy. Each heteroplasmic cell can tolerate high levels of mutated mtDNA with the relative proportion of WT mtDNA to mutated mtDNA determining the fate of cell. Depending on the mutation, the cell can tolerate mutations in up to 60–90% of its mitochondria without any noticeable decline in respiratory efficiency. These somatic mutations accumulate with age and somatically mutated mitochondria have been associated with a myriad of metabolic, autoimmune, and neurodegenerative diseases (1214).

Research on the functional consequences of these mtDNA mutations has largely focused on loss-of-function, such as impairment of oxidative metabolism, whereas there has been minimal attention paid to their role in autoimmune diseases (15). mtDNA mutations lead to higher expression of MHC class I molecules in human cells, and β2-microglobulin–deficient mice have higher accumulation of mt mutations than do WT mice (16), providing indirect evidence of immune surveillance of mt mutations. The significantly higher mtDNA mutation load in rheumatoid arthritis patients compared with osteoarthritis patients is among the mounting evidence that neopeptides generated from mutated mtDNA have a role in autoimmune disease (17). The presence of large-scale deletions in mtDNA in multiple sclerosis patients also adds weight to this hypothesis (18). However, to date the direct recognition of such neopeptides by the immune system has not been reported.

In this study, to our knowledge we show for the first time that neopeptides from mtDNA harboring large-scale frameshift deletions bind to MHC molecules and can be recognized by T cells from normal individuals, providing strong evidence for immune surveillance of altered mitochondria, which in turn may function to limit the accumulation of damaged mitochondria with age. Alternatively, the close similarity of one such mutated mt peptide to both a self-peptide and an environmental microbe also leads us to hypothesize that this immune surveillance mechanism can result in autoimmunity in susceptible individuals.

Human mtDNA sequences were obtained from the National Center for Biotechnology Information (reference sequence ID NC_001807). Deletion mutations in the coding regions of CO1, ATP6, and ND4 were identified using http://www.mitomap.org (19).

The MATLAB Bioinformatics Toolbox was used for sequence manipulations, concatenations, and gene assembly.

SYFPEITHI is a database comprising >7000 known MHC class I and II binding peptide sequences compiled from published reports. SYFPEITHI was used to identify the 9-mers based on their potential to bind to HLA-A*0201 (20).

NetChop produces neural network predictions for cleavage sites of the human proteasome and was used to identify proteasomal cleavage sites in the polypeptides identified in the mtDNA-containing deletion mutations (21).

RankPep predicts peptides that bind to MHC classes I and II from protein sequence or sequence alignments using position-specific scoring matrices. RankPep predicts those MHC class I peptide ligands with a C terminus that is likely to be the result of proteasomal cleavage. It was used to identify potential peptides with C-terminal sequences appropriate for binding and presentation (22).

NetMHCpan predicts binding of peptides to any known MHC using artificial neural networks. The artificial neural networks of NetMHCpan were trained on >115,000 quantitative binding data covering >120 MHCs. Eighty-six percent of NetMHCpan predictions have been experimentally confirmed (23).

NetMHCII 2.2 predicts binding of peptides to MHC class II alleles using artificial neural networks. Predictions can be obtained for 14 HLA-DR alleles (covering the nine HLA-DR supertypes), 6 HLA-DQ, and 6 HLA-DP alleles (23).

Eighteen neopeptides derived from altered mtDNA sequences, the influenza matrix protein-derived peptide (GILGFVFTL), and HIV-gag p17 (SLYNTVATL) were synthesized by Proimmune (Sarasota, FL) at the recommended purity for T cell assays. Aliquots were stored at −80°C.

MHC stabilization assays were performed with T2 cells to test the predicted HLA-A2 binding of the neopeptides. T2 cells are TAP-deficient and thus have few class I MHC proteins on their surface (24, 25). Addition of exogenous peptides stabilizes MHC class I on the surface and can be used as evidence of peptide binding. TAP-deficient T2 cells were cultured in complete medium containing RPMI 1640 (Life Technologies) supplemented with 10% FBS, 0.1% 2-ME, 1× nonessential amino acids (Life Technologies), 1 mM l-glutamine, 1 mM sodium pyruvate, and antibiotics. Cells were incubated in 5% CO2 at 37°C conditions. T2 cells suspended in complete medium were distributed at 2 × 105 per well in a 96-well flat-bottom plate and incubated with each neopeptide overnight at varying concentrations. T2 cells were then washed once in FACS buffer (PBS with 1% [w/v] BSA and 0.1% [w/v] sodium azide) and stained with FITC-labeled anti–HLA-A2 Ab for 30 min at 4°C. The cells were washed twice with FACS buffer, and MHC stabilization was determined using a FACSCalibur. Influenza peptide-treated cells served as a positive control, while T2 cells without peptide were a negative control. The fluorescence indices (FI) were calculated as: FI = (MFI sample − MFI background)/MFI background, where MFI is mean FI, and background is the negative control value. FI is the average of at least two independent experiments.

The Human Participants Review Board and the Biosafety Committee of York University (Toronto, ON, Canada) approved these studies. Twelve healthy human donors between the ages of 22 and 66 y were recruited and provided informed consent. No donor had any recent illnesses (in the past 6 mo) and no donor had any history or symptoms of autoimmune diseases. Approximately 20 ml venous blood was drawn using the BD Vacutainer Safety-Lok blood collection set with a preattached holder (Source Medical Corporation). At the same time, ∼20 μl blood was obtained using BD Microtainer contact-activated lancets (Source Medical Corporation) for the QIAcard FTA spots (Qiagen, Toronto, ON, Canada) and, after drying, were stored in pouches at room temperature.

Immediately after blood draw, PBMCs were purified using Ficoll-Hypaque gradient centrifugation followed by cryopreservation in liquid nitrogen using 40% complete medium, 50% endotoxin-free FBS, and 10% DMSO.

Olerup (Westchester, PA) sequence-specific primer PCR HLA-A and HLA-DR low-resolution kits were used to determine HLA type. Punches from FTA discs were used directly. PCRs with primer sets specific for each allele were performed and the HLA type was determined by the presence and relative sizes of bands after electrophoresis in 0.8% agarose gels using lot-specific specificity tables.

ELISPOT assay kits (human IFN-γ ELISPOTPRO kits, Mabtech, Nacka Strand, Sweden) were used to determine the frequency of epitope-specific IFN-γ–secreting T cells. PBMCs (2 × 105/well/100 μl) were pulsed with peptides (10 μM) overnight in polyvinlylidene difluoride plates precoated with anti-human IFN-γ mAb at 37°C in 5% CO2. After overnight incubation, wells were washed with filtered PBS to remove the non adherent cells. Following washing, each well was treated with 100 μl one-step detection reagent (alkaline phosphatase–conjugated detection IFN-γ Ab) at 1:200 in filtered PBS containing 0.5% FBS for 2 h at 37°C. Plates were washed five more times with filtered PBS before chromogenic development with 100 μl/well ready-to-use 5-bromo-4-chloro-3-indolyl phosphate (BCIP)/NBT. BCIP/NBT-added plates were incubated for 15 min at room temperature, during which time the blue spots developed. Each spot-forming unit (SFU) corresponds to one IFN-γ–secreting T cell. The numbers of Ag-specific T cells are presented as SFU per total number of PBMCs added. The negative control was PBMCs in medium without peptide stimulation and was used to assess the spontaneous secretion of IFN-γ. To verify that the ELISPOT was working, we used the polyclonal activator anti-CD3 mAb as a positive control with each donor. Each peptide was incubated overnight in ELISPOT plates without PBMCs to detect nonspecific reactions. The negative control was PBMCs in medium without peptide stimulation and was used to assess the spontaneous secretion of IFN-γ. A well was considered positive when it contained at least 10 SFU and twice as many SFU as the negative control wells. Each datum point is the mean of duplicate wells.

T cells were expanded with peptide and costimulation, as previously described (26). Briefly, PBMCs were suspended at 7.5 × 106 cells/ml and plated at 0.1 ml/well in 96-well plates. Adherent monocytes were enriched at 37°C for 1 h and then infected with adenovirus containing CD80 (CD80-AdV) or empty adenovirus (CAdV) overnight. Nonadherent cells were removed and stored overnight at 37°C to be used as a source of T cells. CD80-AdV or CAdv was added to the adherent monocytes at a multiplicity of infection of 100, followed by centrifugation at 3000 rpm at 37°C for 1 h. The no Ag control cultures included monocytes treated with CD80-AdV alone. The neopeptides were added to the CD80-Adv and CAdV-modified monocytes at 10 μM. Following overnight incubation, the cultures were washed to remove excess AdV and peptide. T cells were purified from the nonadherent cell fraction by magnetic negative selection using Pan T cell isolation kit II from Miltenyi Biotec (MACS). T cells at 2 × 106 cells/ml were added to the activated monocytes and coincubated at 37°C for 7 d. For the Ab blocking experiments, adenovirus-modified monocytes were incubated with W6/32 (MHC class I blocking Ab) or L243 (HLA-DR blocking Ab) or HLA-DP blocking Ab (clone B7.21) at 10 μg/ml for 1 h prior to the addition of peptide. On day 8, expanded T cells were restimulated with 5 μM peptide for 6 h in the presence of GolgiPlug (BD Biosciences), a protein transport inhibitor. Restimulated cells were first stained for surface markers, fixed, and permeabilized using the Cytofix/Cytoperm kit (BD Biosciences), followed by intracellular cytokine staining. Cell fluorescence was measured using a FACSCalibur and data were analyzed with FlowJo software (Tree Star, Ashland, OR). The following mAbs were used: human (h)-specific anti-hCD3, -hCD8α, -hCD4, and –hIFN-γ (eBioscience). In all figures, fluorescence minus one (FMO) was used to determine background levels of staining. Approximately 1 × 105 to 4 × 105 total T cell events were collected per condition per donor.

MHC tetramers were provided by the National Institute of Allergy and Infectious Diseases MHC Tetramer Core Facility (Atlanta, GA). Tetramers were custom made for HLA-A*0201 with the P9 neopeptide (YVFTLLVPL), the influenza matrix protein-derived peptide (GILGFVFTL), and HIV-gag p17 (SLYNTVATL). T cells were cultured with peptides as described in the costimulation assay. On day 8, cells were first stained with the tetramer for 15 min at 37°C followed by staining for surface markers with anti-hCD3 and anti-hCD8Α (eBioscience) at 4°C for 20 min. In all figures, FMO was used to determine background levels of tetramer staining. FlowJo software (Tree Star) was used for data analysis. Approximately 1 × 105 to 4 × 105 total T cell events were collected per condition per donor.

PBMCs were thawed and rested for 2 h at a density of 2 × 106 cells/well in a round-bottom microtiter plate. Tetramers used were custom made for HLA-A*0201 with the P9 neopeptide (YVFTLLVPL). Cells were first stained with the tetramer for 15 min at 37°C followed by staining for surface markers (anti-hCD3, -hCD8α, -hCD45RA, and hCD27; eBioscience) at 4°C for 20 min. In all figures, FMO was used to determine background levels of tetramer staining. FlowJo software (Tree Star) was used for data analysis.

Unpaired Student t tests were performed using GraphPad Prism software (GraphPad Software, San Diego, CA) to compare T cell responses. A p value of <0.05 was considered significant.

To identify potential immunoreactive neopeptides, we used Mitomap as described in 2Materials and Methods. Of the 116 deletion mutations identified, 101 were selected based on the criteria that the deletions were within an open reading frame and that the fused sequences could be expressed and translated. The sequences were concatenated at the deletion points, and 300 nt flanking each side of the deletion breakpoint, including the mutation site, were computationally translated in the normally used reading frame using the mammalian mt genetic code in the MATLAB Bioinformatics Toolbox. Neopeptides were further screened for their absence in any of the 13 mitochondrially encoded native protein sequences as well as for their predicted binding to a class I MHC allele (using bioinformatics-based methods described in 2Materials and Methods and in Fig. 1). The 101 translated sequences were screened with SYFPEITHI to identify 9-mers that had the potential to bind to HLA-A*0201, the initially selected test allele. For each mutated mtDNA sequence, neopeptides with an SYFPEITHI score of at least 24 were selected as being potentially strong binding affinity neopeptides (27) and were subsequently assessed for proteosomal cleavage probability using NetChop and RankPep (detailed in 2Materials and Methods). Based on these criteria, 18 neopeptides were chosen for experimental testing. Characteristics of the neopeptides identified from the mutant mtDNA are shown in Supplemental Fig. 1. Full details on these mutations are listed in Supplemental Table I. The SYFPEITHI scores of each neopeptide are shown in Fig. 1A.

FIGURE 1.

Computationally predicted HLA-A*02 binding neopeptides from mutated mtDNA and experimental validation of binding to HLA-A*02 allele using MHC stabilization assays. (A) Details of 9-mer neopeptides generated from large-scale frameshift deletions in mtDNA. Information on mtDNA deletions was extracted from Mitomap (19). (B) Neopeptide stabilization of HLA-A2 molecule on TAP-deficient T2 cells was conducted as described in 2Materials and Methods. Data are reported as FI for the average of two independent experiments and each experiment was done in duplicates. FI = (MFI sample − MFI background)/MFI background, where background refers to no peptide stimulation (negative control). The HLA-A*02 binding peptide FluM1 (GILGFVFTL) is the positive control. Data are expressed as MFI ± SD. Representative histograms are shown for FluM1+ control peptide in the inset panel. The complex dashed/dotted, tinted, dashed, dotted, and thick black lines correspond to the MFI of unstained and 0, 5, 10, and 25 μM concentrations of FluM1 peptide, respectively.

FIGURE 1.

Computationally predicted HLA-A*02 binding neopeptides from mutated mtDNA and experimental validation of binding to HLA-A*02 allele using MHC stabilization assays. (A) Details of 9-mer neopeptides generated from large-scale frameshift deletions in mtDNA. Information on mtDNA deletions was extracted from Mitomap (19). (B) Neopeptide stabilization of HLA-A2 molecule on TAP-deficient T2 cells was conducted as described in 2Materials and Methods. Data are reported as FI for the average of two independent experiments and each experiment was done in duplicates. FI = (MFI sample − MFI background)/MFI background, where background refers to no peptide stimulation (negative control). The HLA-A*02 binding peptide FluM1 (GILGFVFTL) is the positive control. Data are expressed as MFI ± SD. Representative histograms are shown for FluM1+ control peptide in the inset panel. The complex dashed/dotted, tinted, dashed, dotted, and thick black lines correspond to the MFI of unstained and 0, 5, 10, and 25 μM concentrations of FluM1 peptide, respectively.

Close modal

To determine whether the 18 predicted neopeptides bind to HLA-A*02, their ability to stabilize the HLA-A*02 molecule on Ag-processing–deficient T2 cells was evaluated (described in 2Materials and Methods). Influenza matrix protein-derived peptide (58GILGFVFTL66; SYFPEITHI score of 30) was used as a positive control. Six of the 18 neopeptides (P4, P9, P10, P16, P17, and P18) demonstrated dose-dependent HLA-A*02 stabilization on T2 cells. The remaining 12 neopeptides displayed minimal effect on HLA-A2 levels and/or failed to show dose-dependent binding.

T cell responses in PBMCs of 12 healthy individuals were tested by overnight stimulation with synthesized neopeptides using the IFN-γ ELISPOTPRO assay (Supplemental Table II shows the details of donors). Because naive T cells cannot be activated to produce IFN-γ within the 18–24 h stimulation time frame, this overnight incubation time ensures that the immune response to these neopeptides is from pre-existing memory T cells in the donors. Fig. 2 shows the response to peptides, divided into A2+ and A2 donors as well as into older and younger donors. All four of the older (>60 y) A2+ donors showed a response of >10 SFU/200,000 cells to at least one of the peptides, with several peptides inducing responses of >40 SFU/200,000 cells in at least one of the older donors (Fig. 2A). Donor OA1 (66 y of age) responded to 6 of 18 neopeptides (33%) tested. Two neopeptides, P9 and P13, stimulated IFN-γ responses in PBMCs from other older donors. Alternatively, little or no response to the neopeptides was observed in the younger HLA-A*02+ donors (mean age of 23.75 y) with the exception of neopeptide P9, which induced a weak response in two younger donors (Fig. 2A). Analysis of age versus frequency of IFN-γ responses to P9 showed a positive correlation with age, although not statistically significant (Supplemental Fig. 2). Although a larger sample size would be required to confirm an age-related response, the results are consistent with T cell responses to neopeptides accumulating with age.

FIGURE 2.

Response of pre-existing memory T cells to neopeptides in IFN-γ ELISPOT assay. (A) HLA-A*02+ older donors and HLA-A*02+ younger donors. (B) HLA-A*02 donors. T cell memory responses were determined in the IFN-γ ELISPOT assay by overnight incubation of donors’ PBMCs with peptides in anti-IFN-γ–coated microtiter plates followed by treatment with anti–IFN-γ detection Ab conjugated to alkaline phosphatase followed with substrate BCIP/NBT treatment. Each blue spot on microtiter plates corresponds to one IFN-γ–secreting T cell (SFU). The positive control polyclonal activator anti-CD3 mAb was used in the functionality of the immunoassay. A well was considered positive when it contained at least 10 SFU and twice as many as SFU as the negative control wells. ELISPOT assays were set up in duplicate. All data are expressed as means ± SD.

FIGURE 2.

Response of pre-existing memory T cells to neopeptides in IFN-γ ELISPOT assay. (A) HLA-A*02+ older donors and HLA-A*02+ younger donors. (B) HLA-A*02 donors. T cell memory responses were determined in the IFN-γ ELISPOT assay by overnight incubation of donors’ PBMCs with peptides in anti-IFN-γ–coated microtiter plates followed by treatment with anti–IFN-γ detection Ab conjugated to alkaline phosphatase followed with substrate BCIP/NBT treatment. Each blue spot on microtiter plates corresponds to one IFN-γ–secreting T cell (SFU). The positive control polyclonal activator anti-CD3 mAb was used in the functionality of the immunoassay. A well was considered positive when it contained at least 10 SFU and twice as many as SFU as the negative control wells. ELISPOT assays were set up in duplicate. All data are expressed as means ± SD.

Close modal

HLA-A*02 donors (mean age of 42 y) also responded to some of the identified neopeptides (Fig. 2B). T cells from donor NA1 (34 y of age) responded to 6 of 18 neopeptides tested, albeit at lower magnitude when compared with the HLA-A*02+ older donor OA1. Other HLA-A*02 donors had T cells with detectable responses to P9 and P13.

Because T cells from HLA-A*02 donors showed responses to the HLA-A*02 binding neopeptides, bioinformatics prediction and HLA typing was used to evaluate the possible restriction of these reactive neopeptides to HLA alleles other than A*02. The analysis revealed that the predicted MHC restriction of neopeptides P9 and P13 matched the HLA-type of A*02 donors NA1 and NA2, likely explaining their response to these peptides (Table I). When neopeptide P9 was screened for its potential binding to HLA-B and HLA-C alleles using the NetMHCpan server, it was predicted to strongly bind to a number of HLA-B and HLA-C alleles (Table II). Additionally, using the NetMHCII 2.2 server, neopeptide P9 was predicted to be presented by class II alleles HLA-DRB1*0101, as well as HLA-DPA1*03:01-DPB1*04:02, HLA-DPA1*01:03-HLA-DPB1*03:01, and HLA-DPA1*01:03-HLA-DPB1*04:01 alleles, albeit with modest predicted binding affinity. These results suggest that some neopeptides could bind promiscuously to several MHC molecules, including restriction of one of the neopeptides (P9) to both an MHC class I and a class II allele.

Table I.
Comparison of HLA-A allele of donor with the HLA-A allele predicted for peptides P9 and P13 using the NetMHCpan server
DonorMHC-I (HLA-A by Typing)MHC Class I Allele Predicted for Peptide P9 (Binding Affinity in nM)MHC Class I Allele Predicted for Peptide P13 (Binding Affinity in nM)
NA1 HLA-A*32 HLA-A*32:01 (42.4) HLA-A*32:05 (74.6) 
HLA-A*32:05 (10.4) 
HLA-A*32:08 (39.8) 
HLA-A*32:12 (42.4) 
NA2 HLA-A*01 HLA-A*31:07 (167.2)  
HLA-A*31 
DonorMHC-I (HLA-A by Typing)MHC Class I Allele Predicted for Peptide P9 (Binding Affinity in nM)MHC Class I Allele Predicted for Peptide P13 (Binding Affinity in nM)
NA1 HLA-A*32 HLA-A*32:01 (42.4) HLA-A*32:05 (74.6) 
HLA-A*32:05 (10.4) 
HLA-A*32:08 (39.8) 
HLA-A*32:12 (42.4) 
NA2 HLA-A*01 HLA-A*31:07 (167.2)  
HLA-A*31 

Binding affinity thresholds of peptides/MHC alleles: strong, ≤50.0 nM; weak, >50.0 to ≤500.0 nM.

Table II.
Predicted MHC class I B and C alleles of peptide P9 using the NetMHCpan server
Peptide NameSequenceMHC Class I B and C Alleles (Predicted Binding Affinity in nM)
P9 YVFTLLVPL HLA-B*15:03 (27.8) 
HLA-B*39:12 (24.8) 
HLA-B*39:14 (22.8) 
HLA-B*39:15 (28.4) 
HLA-B*39:22 (29.6) 
HLA-B*39:27 (17.9) 
HLA-B*39:29 (18.2) 
 
HLA-C*01:04 (22.6) 
HLA-C*02:27 (16.4) 
HLA-C*03:02 (23.3) 
HLA-C*03:06 (22.6) 
HLA-C*03:14 (25.9) 
Peptide NameSequenceMHC Class I B and C Alleles (Predicted Binding Affinity in nM)
P9 YVFTLLVPL HLA-B*15:03 (27.8) 
HLA-B*39:12 (24.8) 
HLA-B*39:14 (22.8) 
HLA-B*39:15 (28.4) 
HLA-B*39:22 (29.6) 
HLA-B*39:27 (17.9) 
HLA-B*39:29 (18.2) 
 
HLA-C*01:04 (22.6) 
HLA-C*02:27 (16.4) 
HLA-C*03:02 (23.3) 
HLA-C*03:06 (22.6) 
HLA-C*03:14 (25.9) 

Binding affinity thresholds of peptides-MHC alleles: strong, ≤ 50.0 nM; weak, >50.0 to ≤500.0 nM.

Because peptide P9 was immunogenic in most of the donors (with 75% of donors showing a response in the ELISPOT assay), we further evaluated this response. T cells from HLA-A*02+ and HLA*A-02 P9-responsive donors were expanded for 7 d in the presence of neopeptide P9 and costimulation. On day 8, T cells were briefly restimulated with peptide and the IFN-γ–secreting cells were analyzed by flow cytometry. The assays demonstrated that both CD3+CD8+ and CD3+CD4+ T cells from donors OA1 and NA1 responded to neopeptide P9 (Fig. 3A and 3B, respectively).

FIGURE 3.

Characterization of T cells. IFN-γ secretion by CD3+CD4+ and CD3+CD8+ T cells responding to neopeptide P9 by HLA-A*02+ older donor, OA1 (A) or (B) HLA-A*02 donor, NA1 in the presence or absence of the MHC class I blocking Ab, W6/32, or class II HLA-DR blocking Ab, L243. Representative flow cytometry plots from two independent experiments are shown for each donor. Numbers in the gated areas correspond to percentages. Each experiment was set up in duplicate. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) and relative to culture treated with class I HLA blocking Ab (W6/32; 10 μg/ml) or class II HLA-DR blocking Ab (L243; 10 μg/ml) prior to the addition of neopeptide but including costimulation (P9 + CD80 + W6/32 or L243). *p < 0.05, **p < 0.01, ***p < 0.001 by unpaired t test. T cells were cultured on peptide-pulsed autologous monocytes modified with replication defective adenovirus expressing CD80 or with control adenovirus (CAdv). Additional controls included monocytes modified with adenovirus expressing CD80 without Ag (NoAg + CD80). Cells were harvested on day 8 and restimulated for 6 h with P9 followed by surface staining, fixation and permeabilization, and intracellular cytokine staining (details in 2Materials and Methods). Flow cytometric analyses were performed on a lymphocyte-gated population. Gating was based on FMO controls for each stain.

FIGURE 3.

Characterization of T cells. IFN-γ secretion by CD3+CD4+ and CD3+CD8+ T cells responding to neopeptide P9 by HLA-A*02+ older donor, OA1 (A) or (B) HLA-A*02 donor, NA1 in the presence or absence of the MHC class I blocking Ab, W6/32, or class II HLA-DR blocking Ab, L243. Representative flow cytometry plots from two independent experiments are shown for each donor. Numbers in the gated areas correspond to percentages. Each experiment was set up in duplicate. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) and relative to culture treated with class I HLA blocking Ab (W6/32; 10 μg/ml) or class II HLA-DR blocking Ab (L243; 10 μg/ml) prior to the addition of neopeptide but including costimulation (P9 + CD80 + W6/32 or L243). *p < 0.05, **p < 0.01, ***p < 0.001 by unpaired t test. T cells were cultured on peptide-pulsed autologous monocytes modified with replication defective adenovirus expressing CD80 or with control adenovirus (CAdv). Additional controls included monocytes modified with adenovirus expressing CD80 without Ag (NoAg + CD80). Cells were harvested on day 8 and restimulated for 6 h with P9 followed by surface staining, fixation and permeabilization, and intracellular cytokine staining (details in 2Materials and Methods). Flow cytometric analyses were performed on a lymphocyte-gated population. Gating was based on FMO controls for each stain.

Close modal

To confirm that the response was MHC-specific, we included the pan-MHC class I–reactive Ab W6/32 in the cultures to block MHC class I recognition. Indeed, W6/32 addition greatly reduced the expansion of CD3+CD8+ T cells in response to P9, whereas, as expected, the CD4+ T cell response to this peptide was not blocked (Fig. 3). Alternatively, the anti-DR Ab L243 did not significantly reduce the CD4+ T cell response, indicating restriction to alleles other than HLA-DR (Fig. 3). Based on our bioinformatics analysis described above, it is likely that the CD4+ T cell response to these peptides is DP-restricted. Treatment of APCs with class II HLA-DP blocking Ab (clone B7.21), prior to the addition of neopeptide P9, resulted in significant reduction in the expansion of CD3+CD4+ T cells in response to P9 (Fig. 4). This observation confirms the potential class II HLA restriction of CD3+CD4+ T cells recognizing P9. Taken together, these data confirm the MHC class I and II–restricted recognition of P9 in both A*02+ and A*02 donors and provide evidence that P9 can be promiscuously presented to both CD4+ and CD8+ T cells.

FIGURE 4.

The CD3+CD4+ T cells response to neopeptide P9 is HLA-DP restricted. IFN-γ production in response to neopeptide P9 by HLA-A*02+ older donor, OA1, is shown. Representative flow cytometry data as well as summary bar graphs from two experiments are shown. Numbers in the gated areas correspond to percentages. Each experiment was set up in duplicate. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) and relative to culture treated with class II HLA-DP blocking Ab (HLA-DP (clone B7.21; 10 μg/ml) prior to the addition of neopeptide but including costimulation (P9 + CD80 + HLA-DP). *p < 0.05 by unpaired t test. T cells from donor were cultured on peptide-pulsed autologous monocytes modified with either adenovirus expressing CD80 or control Adv (CAdv) as described in Fig. 3 and 2Materials and Methods.

FIGURE 4.

The CD3+CD4+ T cells response to neopeptide P9 is HLA-DP restricted. IFN-γ production in response to neopeptide P9 by HLA-A*02+ older donor, OA1, is shown. Representative flow cytometry data as well as summary bar graphs from two experiments are shown. Numbers in the gated areas correspond to percentages. Each experiment was set up in duplicate. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) and relative to culture treated with class II HLA-DP blocking Ab (HLA-DP (clone B7.21; 10 μg/ml) prior to the addition of neopeptide but including costimulation (P9 + CD80 + HLA-DP). *p < 0.05 by unpaired t test. T cells from donor were cultured on peptide-pulsed autologous monocytes modified with either adenovirus expressing CD80 or control Adv (CAdv) as described in Fig. 3 and 2Materials and Methods.

Close modal

To further verify that the T cells in the cultures were responding specifically to the neopeptide P9, we analyzed the binding of expanded T cells to HLA-A*02 tetramers containing the P9 peptide. The tetramer specificity was verified using PBMCs from an HLA-A*02 donor (NA1) (Fig. 5). CD3+CD8+ T cells expanded using P9 peptide reacted with the P9/A*02 tetramer (1.59%) whereas there was minimal tetramer binding to T cells expanded from a young healthy donor (YA3) who was nonresponsive to P9 in the IFN-γ ELISPOTPRO (0.422%), and there was minimal tetramer binding to T cells from an A2 donor (NA1, 0.282%). Moreover, P9-expanded T cells did not bind an HIV-gag p17 (SLYNTVATL)/A*02 tetramer. Thus, healthy donor OA1 has pre-existing memory CD8+ T cells that can be expanded by restimulation with peptide P9, and these cells bind specifically to the A*02/P9 tetramer.

FIGURE 5.

CD8+ T cells expanded in response to P9 neopeptide show specific binding to the HLA-A*0201/YVFTLLVPL (P9) tetramer. (A) Representative flow cytometry plot of three independent experiments. Each experiment was set up in duplicate. Numbers in the gated area of the dot plots correspond to percentages. (B) Summary of CD3+CD8+ tetramer+ T cells in response to peptide P9. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) for OA1 donor. CD3+CD8+ P9 tetramer responses were compared between OA1 and YA3 donors and OA1 and NA1 donors. T cells were cultured as described in costimulation assay (see Materials and Methods) and Fig. 3. On day 8, cells were harvested and first stained with the tetramer followed by staining for surface markers. Flu-M1 peptide (GILGFVFTL) tetramer was used as a positive control for HLA-A*02 donors. HIV-gag p17 (SLYNTVATL) tetramer was used as a negative control. Flow cytometric analyses were performed on a lymphocyte-gated population. In all figures, FMO control was used to determine background levels of tetramer staining. *p < 0.05, **p < 0.01 by unpaired t test.

FIGURE 5.

CD8+ T cells expanded in response to P9 neopeptide show specific binding to the HLA-A*0201/YVFTLLVPL (P9) tetramer. (A) Representative flow cytometry plot of three independent experiments. Each experiment was set up in duplicate. Numbers in the gated area of the dot plots correspond to percentages. (B) Summary of CD3+CD8+ tetramer+ T cells in response to peptide P9. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) for OA1 donor. CD3+CD8+ P9 tetramer responses were compared between OA1 and YA3 donors and OA1 and NA1 donors. T cells were cultured as described in costimulation assay (see Materials and Methods) and Fig. 3. On day 8, cells were harvested and first stained with the tetramer followed by staining for surface markers. Flu-M1 peptide (GILGFVFTL) tetramer was used as a positive control for HLA-A*02 donors. HIV-gag p17 (SLYNTVATL) tetramer was used as a negative control. Flow cytometric analyses were performed on a lymphocyte-gated population. In all figures, FMO control was used to determine background levels of tetramer staining. *p < 0.05, **p < 0.01 by unpaired t test.

Close modal

Having established the tetramer specificity of HLA-A*02 tetramer of P9, we next sought to characterize the surface phenotype of P9 T cells immediately ex vivo to determine whether the tetramer reactivity was due to naive or Ag-experienced cells. The expression of the leukocyte common Ag isoform, CD45RA, as well as the costimulatory molecule CD27, has been used to characterize CD8 T cell subsets. The CD45RA+CD27+ subset of CD8 T cells is thought to represent the naive compartment whereas the memory compartment is represented by CD45RACD27+ (central memory cells), CD45RACD27- (effector memory T cells), and CD45RA+CD27 (highly differentiated effector memory population) cells, respectively (2830).

P9/A*02 tetramer+CD8+ T cells largely consisted of CD45RA memory cells, with most exhibiting a CD45RACD27+ memory phenotype as shown in Fig. 6. These results support the conclusion that P9-reactive cells in donor OA1 consist of Ag-experienced T cells with a memory phenotype.

FIGURE 6.

Characterization of native phenotype of T cells recognizing neopeptide P9. (A) Representative flow cytometry plot of two independent experiments. Each experimental condition was done in duplicate in each of the two experiments. Numbers in the gated area of the dot plots correspond to percentages. (B) FMO control was used to determine background levels of tetramer staining. PBMCs were thawed and rested for 2 h at a density of 2 × 106 cells/well in a round-bottom microtiter plate. Cells were first stained with tetramer followed by staining for surface markers as described in 2Materials and Methods.

FIGURE 6.

Characterization of native phenotype of T cells recognizing neopeptide P9. (A) Representative flow cytometry plot of two independent experiments. Each experimental condition was done in duplicate in each of the two experiments. Numbers in the gated area of the dot plots correspond to percentages. (B) FMO control was used to determine background levels of tetramer staining. PBMCs were thawed and rested for 2 h at a density of 2 × 106 cells/well in a round-bottom microtiter plate. Cells were first stained with tetramer followed by staining for surface markers as described in 2Materials and Methods.

Close modal

Most donors tested (75%, 9 of 12; including HLA-A*02 and HLA–non-A*02) showed reactivity to neopeptide P9 in the ELISPOT assay and a substantial proportion of this reactivity was due to CD27 memory T cells. How could such reactivity arise? Reactivity could be due to an initial recognition of mutant mt peptides and the ensuing maintenance of memory. The reactivity could also be the result of priming of T cells with peptides derived from environmental organisms. We searched for homology to P9 within microbes (using BLASTP, National Center for Biotechnology Information). Remarkably, neopeptide P9 has 89% similarity with Bacillus pumilus fructose PTS peptide, YLFTLLVPL (P9: YVFTLLVPL), differing only at position 2, with an SYEPEITHI score of 30 for HLA-A*02. B. pumilus is a soil-dwelling, Gram+, spore-forming bacterium and has been found as a bloodstream pathogen during infancy (31, 32). A homology search with human self-protein also revealed that neopeptide P9 shares 89% similarity with mt self-protein, ATP synthase F0 subunit 6 (YVFTLLVSL; P9: YVFTLLVPL), again differing at only one amino acid and with an SYEPEITHI score of 26.

This analysis raises the possibility that P9 peptide–specific memory T cells could have been initially primed via cross-reactive responses to bacterial peptides. The findings also raise the possibility that T cells activated by a neopeptide could also cross-react with the corresponding host self-proteins.

The PBMC T cell responses of the HLA-A*02 OA1 healthy individual (66 y of age) were evaluated by overnight stimulation with the synthesized self-peptide and the bacterial homolog peptide of neopeptide (P9) using the IFN-γ ELISPOT assay. This short incubation time ensures that measured responses are recall responses from memory T cells. IFN-γ T cell response was measured as number of SFU, where each spot corresponds to a T cell. PBMCs stimulated with self-peptide and bacterial homolog peptide of P9 resulted in positive IFN-γ T cell responses (47.5 SFU/200,000 PBMCs, SEM 3.227 and 15.25 SFU/200,000 PBMCs, SEM 1.25). The bacterial homolog peptide was able to generate a T cell response, which was perhaps not unexpected given the fact that it is a foreign peptide. More surprising, however, was the finding that the self-peptide generated T cell responses. This result demonstrates that there were pre-existing T cells in the peripheral blood that recognized self-peptide from mtDNA. OA1 donor also responded to P9 and was shown to have P9-specific T cells based on tetramer binding. This observation raises a testable hypothesis of cross-reactivity among P9 and its self and the bacterial homolog.

To determine whether an immunological cross-reactivity exists among P9 and its homolog peptides, additional stimulation assays were carried out. T cells from donor OA1 were expanded in the presence of self and bacterial homolog peptides, respectively, and subsequently stained with the P9-specific tetramer. T cells expanded with the self and bacterial homolog peptides showed staining with P9-specific HLA-A*02 tetramer (Fig. 7). This observation led to the conclusion that T cells recognizing self and bacterial homologs of P9 can cross-react with neopeptide P9-expressing APCs.

FIGURE 7.

CD8+ T cells expanded by self and bacterial homologs of neopeptide P9 peptide bind A*02/P9 tetramers. (A) Representative flow cytometry plot of two independent experiments. Each experimental condition was done in duplicate. Numbers in the gated area of the dot plots correspond to percentages. (B) Summary of CD3+CD8+ A*02/P9 tetramer+ T cells. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) for each of the peptides (P9, self, and bacterial homologs of P9). T cells from donor OA1 were expanded in the presence of self and bacterial homolog peptides, respectively, and subsequently stained with the P9-specific tetramer, as described in costimulation assay (see 2Materials and Methods) and Fig. 3. Flow cytometric analyses were performed on a lymphocyte-gated population. In all figures, FMO control was used to determine background levels of tetramer staining. *p < 0.05, **p < 0.01 by unpaired t test.

FIGURE 7.

CD8+ T cells expanded by self and bacterial homologs of neopeptide P9 peptide bind A*02/P9 tetramers. (A) Representative flow cytometry plot of two independent experiments. Each experimental condition was done in duplicate. Numbers in the gated area of the dot plots correspond to percentages. (B) Summary of CD3+CD8+ A*02/P9 tetramer+ T cells. All data are expressed as means ± SEM. Statistical significance was analyzed relative to no Ag control, but including costimulation (NoAg + CD80) for each of the peptides (P9, self, and bacterial homologs of P9). T cells from donor OA1 were expanded in the presence of self and bacterial homolog peptides, respectively, and subsequently stained with the P9-specific tetramer, as described in costimulation assay (see 2Materials and Methods) and Fig. 3. Flow cytometric analyses were performed on a lymphocyte-gated population. In all figures, FMO control was used to determine background levels of tetramer staining. *p < 0.05, **p < 0.01 by unpaired t test.

Close modal

These data provide proof of principal that mitochondria neopeptides can be recognized by the human immune system and raise the possibility that this recognition could be initially due to molecular mimicry and could ultimately lead to recognition of the homologous self-peptide.

Human health, in large part, is the result of a functional and balanced immune system, an immune system that surveys presented epitopes and eliminates pathogens while minimizing harmful autoimmune activation. The combined bioinformatics/experimental investigation presented in this study provides strong evidence that neopeptides generated from large-scale deletion mutations in mtDNA can be presented in an MHC-restricted but promiscuous fashion, and can elicit recall T cell responses from healthy individuals. Collectively, this study provides a new paradigm of immune recognition of mt neopeptides, namely, that altered mitochondria are subject to immune surveillance. The close similarity of one of the peptides to an epitope from nonmutated mt proteins raises the possibility that mt neopeptides could contribute to autoimmune diseases.

One of the tested peptides, neopeptide P9, was recognized by 75% of the healthy donors tested, with a trend toward more robust responses in older donors. P9 was promiscuous in its restriction, as it was presented to CD8+ T cells by HLA-A*02 and other MHC I alleles, as well as stimulating CD4+ T cell responses. The promiscuity of P9 notwithstanding, the question remains as to how pre-existing immunity to P9 (and other neopeptides) arose in the healthy donors. Bioinformatics analyses and experimental investigations identified possible sources. The common soil microbe, B. pumilus, has an encoded 9-mer that differs by one amino acid from P9. Similarly, the ATP synthase F0 subunit 6 from normal human mitochondria encodes a 9-mer with a single amino acid difference from P9 with 89% homology to P9. Both peptides have SYEPEITHI scores predictive of strong HLA-A*0201 binding. These observations raise a possible scenario whereby T cells respond to B. pumilus and give rise to immune responses that then cross-react with mt neopeptides that arise from mutation. These T cells in turn may then respond to the peptides from nonmutated mt proteins. Alternatively, the mt neopeptides could be the initiating immunogens. In support of cross-reactivity being an instrumental mechanism in priming self-recognition, PBMCs from OA1 generated IFN-γ T cell response when stimulated with the self and the bacterial homolog peptides of P9. Furthermore, T cells expanded on the self and bacterial homolog peptides were able to bind HLA-A*02 tetramers loaded with P9 peptide, raising the potential for T cell cross-reactivity.

Although the precise mechanism of processing and MHC loading of the peptides was not investigated, the detection of mitochondrially derived peptides on mammalian cells in the context of MHC suggests that such peptides released into cytosol enter conventional MHC Ag presentation pathways, leading to recognition by T cells (3, 4). The finding that ABC transporters have been shown to aid in the efflux of peptides from mitochondria supports this concept. Peptides from mitochondria can be released either by passive diffusion through porins or by the translocase of the outer mt membrane (33, 34). The neopeptides would be delivered to the cytosol and to the endoplasmic reticulum by the TAP transporters in the same way as virally derived peptides. Indeed, Gu et al. (16) demonstrated increased expression of MHC class I in the cells harboring mtDNA deletions. In the case of the class II pathway, once the mt peptides enter the cytosol, they could be recruited into the endosome pathway by autophagy, a well-established mechanism by which viral Ags enter the MHC class II pathway in infected cells (3537). Autophagy is the major mechanism by which thymic epithelial cells present Ags, including soluble mt Ags in the class II pathway for negative selection (3840). It would be expected that the immune system would be tolerant to normal mt peptides delivered to the class II molecule by autophagy. Thus, we surmise that exposure to either the altered mt peptides or the bacterial homolog may allow recognition of these peptides by T cells and potentially break tolerance to the highly similar self–mt peptides.

Direct repeats flanking deleted mtDNA (4143), the presence of deletion breakpoints within or near the predicted non-B DNA structures (44), and the formation of long-distance stable imperfect duplexes between long-distance homologous mtDNA segments (11) have been implicated in generating mtDNA deletions. Additional mechanisms include models implicating replication errors in mtDNA (4547) and deletions during repair of damaged mtDNA (48).

Prima facie large deletions on the order of several kilobases appear to render the mitochondrion translationally incompetent. However, a single mitochondrion contains multiple coexisting genomes, and thus deleted mt genomes can be maintained by transcomplementation (4951). Pioneering studies on mt deletions have demonstrated that these deleted regions continue to be transcribed (49). Because the major method of degrading nuclear transcripts, nonsense-mediated decay, does not seem to apply to mt transcripts, it is likely that the mutant transcripts persist (5255). The question remains then as to whether these regions are translated. In support of the presence of the altered polypeptides, Tiranti et al. (56) observed that the mutant transcripts of mtDNA COIII gene formed by a novel frameshift mutation generated prematurely truncated polypeptide (110 aa). Thus, mutant mtDNAs can be retained in the mitochondria, and they will generate translational products that can elicit an immune response.

The high frequency of responses in this study suggests that mt mutations and their recognition by the human immune system are frequent occurrences. However, we were unable to detect mtDNA with these deletion mutations in the blood of healthy individuals by PCR amplification of mtDNA sequences (including OA1; data not shown). Although other tissues might harbor higher levels of mutated mitochondria, the lack of detection of these mutations in healthy individuals is consistent with their rapid elimination by the immune system. Because T cells are exquisitely sensitive to low numbers of MHC/peptide complexes (5759), it may take very few mutant mitochondria to initiate an immune response.

Based on IFN-γ T cell responses detected after overnight stimulation with neopeptides, there appeared to be pre-existing T cell responses to neopeptide P9. Consistently, CD8+tetramer+ T cells for P9 analyzed immediately ex vivo were predominantly of a CD45RA memory phenotype. Mitochondrial mutations accumulate with age (60). It is plausible that exposure to neopeptide P9 arising from mt mutations in the older donor (OA1, 66 y of age) led to the induction of the P9-reactive T cells. A small fraction of the P9 tetramer-bright cells were of the CD45RA+CD27+ “naive” phenotype. Although it is possible that this represents the naive repertoire, it has also been reported that CD45RA+CD27+ T cells from older subjects include a non-naive compartment as shown by the presence of CD11a+CD127 cells in this subset of T cells (61). Regardless, it is clear that >80% of the P9-reactive T cells has a phenotype consistent with prior activation.

The frequency and magnitude of T cell memory responses toward mt neopeptides showed a weak, but statistically nonsignificant correlation with age. Such a correlation would be consistent with other studies demonstrating the accumulation of mt mutations with age; however, further studies with a larger number of donors would be required to substantiate this conclusion.

The fact that T cell memory responses to mt neopeptides are from healthy individuals suggests a scenario whereby a robust immune surveillance mechanism is functioning to detect neopeptides arising from mutant mtDNA. Surveillance would play a significant role in preventing the accumulation of cells with inefficient mitochondria, which have the potential to lead to mt-related diseases. The robustness of the immune system in the recognition of mtDNA neopeptides is exemplified by the demonstration of MHC promiscuity and CD8+, CD4+ T cell responses toward these peptides. A single frameshift deletion can lead to many potential immunogenic neopeptides. A continuous immune response to these mtDNA neopeptides could result in and/or perpetuate an ongoing chronic inflammation state, thus leading to considerable damage to affected and unaffected tissues. The inflammation could in turn cause damage to mtDNA, resulting in a vicious cycle of an increased production of reactive oxygen species and further damage to mtDNA and the maintenance of immunogenic neopeptides, leading to potential immune-mediated damage.

A number of studies have linked mt deletions with cell death and cell death with inflammation, a characteristic of some autoimmune diseases (62, 63). Mirabella et al. (64) found apoptotic changes in the muscle fibers of 36 assayed muscle biopsies from patients with mtDNA deletions and tRNA point mutations. Similar findings in vitro were reported by Asoh et al. (65). Because apoptotic cell death is an active, energy-dependent process, defects in mt function render cells more susceptible to apoptosis. Moreover, cell death has been shown to cause the release of mt Ags that can induce inflammation. A link between inflammatory mt Ags and autoimmunity has been uncovered by a number of researchers (62, 63). Campbell et al. (18) found abundant mt deletions in multiple sclerosis patients. They observed that neurons with high heteroplasmy and mtDNA deletions were evident in both young and aged multiple sclerosis patients. Neurodegeneration would be the result of respiratory deficiency in the neurons and/or chronic inflammation caused by the neopeptides encoded by mutated DNAs (mtDNA, pathogens’ DNA). An increased mutational burden in autoimmune rheumatoid arthritis patients when compared with osteoarthritis patients further adds support to the concept that neopeptides from mutated mtDNA are implicated in the pathogenesis of inflammatory autoimmunity (17, 66, 67).

Most studies to date implicating mitochondria in inflammation and autoimmunity have focused on the role of innate responses of monocytes/macrophages and neutrophils to the mt components released during cell death (62, 63, 67). The data presented in the present study provide evidence that the adaptive immune system, namely CD8+ and CD4+ T cells, respond to mtDNA-derived neopeptides. Mitochondria, unlike microbial pathogens, are present within the host cells, and the mutations in mtDNA can result in the persistence of immunogenic peptides and thus contribute to the perpetuation of inflammation.

These data provide evidence of a robust immune surveillance mechanism in humans that can detect changes in mtDNA through the recognition of mtDNA-derived neopeptides by effector T cells. This immune surveillance could be preventing the pathological accumulation of cells with defective mitochondria but may also contribute to the onset or exacerbation of autoimmunity in susceptible individuals due to immunological cross-reactivity between mt neopeptides and self-homologs.

This work was supported by Arthritis Society Grant 560102, as well as by Canadian Institute of Health Research grants to G.E.W., Canadian Institute of Health Research Grant MOP-74492 (to T.H.W.), a Natural Sciences and Engineering Research Council of Canada grant to J.G., a York University research grant to V.J., a Natural Sciences and Engineering Research Council of Canada grant and Mathematics of Information Technology and Complex Systems/Mprime grant to V.R.D. and J.W., and a Canadian Research Chairs grant to J.W. B.D. was supported by the President Susan Mann Dissertation Scholarship, the Ontario Graduate Scholarship, and The Arthritis Society/Canadian Arthritis Network Post Doctoral Fellowship. T.H.W. holds the Sanofi Pasteur Chair in Human Immunology at the University of Toronto. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

The online version of this article contains supplemental material.

Abbreviations used in this article:

BCIP

5-bromo-4-chloro-3-indolyl phosphate

CAdV

empty adenovirus

CD80-AdV

adenovirus containing CD80

FI

fluorescence index

FMO

fluorescence minus one

h

human

MFI

mean FI

mt

mitochondrial

SFU

spot-forming unit

WT

wild-type.

1
Anderson
S.
,
Bankier
A. T.
,
Barrell
B. G.
,
de Bruijn
M. H.
,
Coulson
A. R.
,
Drouin
J.
,
Eperon
I. C.
,
Nierlich
D. P.
,
Roe
B. A.
,
Sanger
F.
, et al
.
1981
.
Sequence and organization of the human mitochondrial genome.
Nature
290
:
457
465
.
2
Legros
F.
,
Malka
F.
,
Frachon
P.
,
Lombès
A.
,
Rojo
M.
.
2004
.
Organization and dynamics of human mitochondrial DNA.
J. Cell Sci.
117
:
2653
2662
.
3
Loveland
B.
,
Wang
C. R.
,
Yonekawa
H.
,
Hermel
E.
,
Lindahl
K. F.
.
1990
.
Maternally transmitted histocompatibility antigen of mice: a hydrophobic peptide of a mitochondrially encoded protein.
Cell
60
:
971
980
.
4
Fischer Lindahl
K.
,
Hermel
E.
,
Loveland
B. E.
,
Wang
C. R.
.
1991
.
Maternally transmitted antigen of mice: a model transplantation antigen.
Annu. Rev. Immunol.
9
:
351
372
.
5
Gray
M. W.
,
Burger
G.
,
Lang
B. F.
.
1999
.
Mitochondrial evolution.
Science
283
:
1476
1481
.
6
Timmis
J. N.
,
Ayliffe
M. A.
,
Huang
C. Y.
,
Martin
W.
.
2004
.
Endosymbiotic gene transfer: organelle genomes forge eukaryotic chromosomes.
Nat. Rev. Genet.
5
:
123
135
.
7
Lindenborn-Fotinos
J.
,
Sayers
T. J.
,
Berg
P. A.
.
1982
.
Mitochondrial antibodies in primary biliary cirrhosis. VI. Association of the complement fixing antigen with a component of the mitochondrial F1-ATPase complex.
Clin. Exp. Immunol.
50
:
267
274
.
8
Lindenborn-Fotinos
J.
,
Baum
H.
,
Berg
P. A.
.
1985
.
Mitochondrial antibodies in primary biliary cirrhosis: species and nonspecies specific determinants of M2 antigen.
Hepatology
5
:
763
769
.
9
Yakes
F. M.
,
Van Houten
B.
.
1997
.
Mitochondrial DNA damage is more extensive and persists longer than nuclear DNA damage in human cells following oxidative stress.
Proc. Natl. Acad. Sci. USA
94
:
514
519
.
10
Zastawny
T. H.
,
Dabrowska
M.
,
Jaskolski
T.
,
Klimarczyk
M.
,
Kulinski
L.
,
Koszela
A.
,
Szczesniewicz
M.
,
Sliwinska
M.
,
Witkowski
P.
,
Olinski
R.
.
1998
.
Comparison of oxidative base damage in mitochondrial and nuclear DNA.
Free Radic. Biol. Med.
24
:
722
725
.
11
Guo
X.
,
Popadin
K. Y.
,
Markuzon
N.
,
Orlov
Y. L.
,
Kraytsberg
Y.
,
Krishnan
K. J.
,
Zsurka
G.
,
Turnbull
D. M.
,
Kunz
W. S.
,
Khrapko
K.
.
2010
.
Repeats, longevity and the sources of mtDNA deletions: evidence from “deletional spectra”.
Trends Genet.
26
:
340
343
.
12
Cortopassi
G. A.
,
Shibata
D.
,
Soong
N. W.
,
Arnheim
N.
.
1992
.
A pattern of accumulation of a somatic deletion of mitochondrial DNA in aging human tissues.
Proc. Natl. Acad. Sci. USA
89
:
7370
7374
.
13
Reeve
A. K.
,
Krishnan
K. J.
,
Turnbull
D.
.
2008
.
Mitochondrial DNA mutations in disease, aging, and neurodegeneration.
Ann. N. Y. Acad. Sci.
1147
:
21
29
.
14
Larsson
N. G.
2010
.
Somatic mitochondrial DNA mutations in mammalian aging.
Annu. Rev. Biochem.
79
:
683
706
.
15
Craig
I.
1975
.
Letter: Mitochondria, autoimmunity, and oncogenesis.
Lancet
305
:
218
219
.
16
Gu
Y.
,
Wang
C.
,
Roifman
C. M.
,
Cohen
A.
.
2003
.
Role of MHC class I in immune surveillance of mitochondrial DNA integrity.
J. Immunol.
170
:
3603
3607
.
17
Da Sylva
T. R.
,
Connor
A.
,
Mburu
Y.
,
Keystone
E.
,
Wu
G. E.
.
2005
.
Somatic mutations in the mitochondria of rheumatoid arthritis synoviocytes.
Arthritis Res. Ther.
7
:
R844
R851
.
18
Campbell
G. R.
,
Ziabreva
I.
,
Reeve
A. K.
,
Krishnan
K. J.
,
Reynolds
R.
,
Howell
O.
,
Lassmann
H.
,
Turnbull
D. M.
,
Mahad
D. J.
.
2011
.
Mitochondrial DNA deletions and neurodegeneration in multiple sclerosis.
Ann. Neurol.
69
:
481
492
.
19
Ruiz-Pesini
E.
,
Lott
M. T.
,
Procaccio
V.
,
Poole
J. C.
,
Brandon
M. C.
,
Mishmar
D.
,
Yi
C.
,
Kreuziger
J.
,
Baldi
P.
,
Wallace
D. C.
.
2007
.
An enhanced MITOMAP with a global mtDNA mutational phylogeny.
Nucleic Acids Res.
35
(
Database issue
):
D823
D828
.
20
Rammensee
H.
,
Bachmann
J.
,
Emmerich
N. P.
,
Bachor
O. A.
,
Stevanović
S.
.
1999
.
SYFPEITHI: database for MHC ligands and peptide motifs.
Immunogenetics
50
:
213
219
.
21
Nielsen
M.
,
Lundegaard
C.
,
Lund
O.
,
Keşmir
C.
.
2005
.
The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage.
Immunogenetics
57
:
33
41
.
22
Reche
P. A.
,
Glutting
J. P.
,
Reinherz
E. L.
.
2002
.
Prediction of MHC class I binding peptides using profile motifs.
Hum. Immunol.
63
:
701
709
.
23
Nielsen
M.
,
Lundegaard
C.
,
Blicher
T.
,
Lamberth
K.
,
Harndahl
M.
,
Justesen
S.
,
Røder
G.
,
Peters
B.
,
Sette
A.
,
Lund
O.
,
Buus
S.
.
2007
.
NetMHCpan, a method for quantitative predictions of peptide binding to any HLA-A and -B locus protein of known sequence.
PLoS ONE
2
:
e796
.
24
Henderson
R. A.
,
Michel
H.
,
Sakaguchi
K.
,
Shabanowitz
J.
,
Appella
E.
,
Hunt
D. F.
,
Engelhard
V. H.
.
1992
.
HLA-A2.1-associated peptides from a mutant cell line: a second pathway of antigen presentation.
Science
255
:
1264
1266
.
25
Luft
T.
,
Rizkalla
M.
,
Tai
T. Y.
,
Chen
Q.
,
MacFarlan
R. I.
,
Davis
I. D.
,
Maraskovsky
E.
,
Cebon
J.
.
2001
.
Exogenous peptides presented by transporter associated with antigen processing (TAP)-deficient and TAP-competent cells: intracellular loading and kinetics of presentation.
J. Immunol.
167
:
2529
2537
.
26
Bukczynski
J.
,
Wen
T.
,
Ellefsen
K.
,
Gauldie
J.
,
Watts
T. H.
.
2004
.
Costimulatory ligand 4-1BBL (CD137L) as an efficient adjuvant for human antiviral cytotoxic T cell responses.
Proc. Natl. Acad. Sci. USA
101
:
1291
1296
.
27
Williams
D. S.
,
Bird
M. J.
,
Jorissen
R. N.
,
Yu
Y. L.
,
Walker
F.
,
Zhang
H. H.
,
Nice
E. C.
,
Burgess
A. W.
.
2010
.
Nonsense mediated decay resistant mutations are a source of expressed mutant proteins in colon cancer cell lines with microsatellite instability.
PLoS ONE
5
:
e16012
.
28
Hamann
D.
,
Baars
P. A.
,
Rep
M. H.
,
Hooibrink
B.
,
Kerkhof-Garde
S. R.
,
Klein
M. R.
,
van Lier
R. A.
.
1997
.
Phenotypic and functional separation of memory and effector human CD8+ T cells.
J. Exp. Med.
186
:
1407
1418
.
29
Sallusto
F.
,
Geginat
J.
,
Lanzavecchia
A.
.
2004
.
Central memory and effector memory T cell subsets: function, generation, and maintenance.
Annu. Rev. Immunol.
22
:
745
763
.
30
Mallard
E.
,
Vernel-Pauillac
F.
,
Velu
T.
,
Lehmann
F.
,
Abastado
J. P.
,
Salcedo
M.
,
Bercovici
N.
.
2004
.
IL-2 production by virus- and tumor-specific human CD8 T cells is determined by their fine specificity.
J. Immunol.
172
:
3963
3970
.
31
Priest
F. G.
1993
.
Systematics and ecology of Bacillus
. In
Bacillus subtilis and Other Gram-Positive Bacteria: Biochemistry, Physiology, and Molecular Genetics.
Sonenshein
A. L.
,
Hoch
J. A.
,
Losick
R.
, eds.
ASM Press
,
Washington, D.C.
, p.
3
16
.
32
Kimouli
M.
,
Vrioni
G.
,
Papadopoulou
M.
,
Koumaki
V.
,
Petropoulou
D.
,
Gounaris
A.
,
Friedrich
A. W.
,
Tsakris
A.
.
2012
.
Two cases of severe sepsis caused by Bacillus pumilus in neonatal infants.
J. Med. Microbiol.
61
:
596
599
.
33
Mannella
C. A.
1992
.
The “ins” and “outs” of mitochondrial membrane channels.
Trends Biochem. Sci.
17
:
315
320
.
34
Hill
K.
,
Model
K.
,
Ryan
M. T.
,
Dietmeier
K.
,
Martin
F.
,
Wagner
R.
,
Pfanner
N.
.
1998
.
Tom40 forms the hydrophilic channel of the mitochondrial import pore for preproteins.
Nature
395
:
516
521
.
35
Dengjel
J.
,
Schoor
O.
,
Fischer
R.
,
Reich
M.
,
Kraus
M.
,
Müller
M.
,
Kreymborg
K.
,
Altenberend
F.
,
Brandenburg
J.
,
Kalbacher
H.
, et al
.
2005
.
Autophagy promotes MHC class II presentation of peptides from intracellular source proteins.
Proc. Natl. Acad. Sci. USA
102
:
7922
7927
.
36
Levine
B.
,
Deretic
V.
.
2007
.
Unveiling the roles of autophagy in innate and adaptive immunity.
Nat. Rev. Immunol.
7
:
767
777
.
37
Crotzer
V. L.
,
Blum
J. S.
.
2009
.
Autophagy and its role in MHC-mediated antigen presentation.
J. Immunol.
182
:
3335
3341
.
38
Nedjic
J.
,
Aichinger
M.
,
Emmerich
J.
,
Mizushima
N.
,
Klein
L.
.
2008
.
Autophagy in thymic epithelium shapes the T-cell repertoire and is essential for tolerance.
Nature
455
:
396
400
.
39
Levine
B.
,
Mizushima
N.
,
Virgin
H. W.
.
2011
.
Autophagy in immunity and inflammation.
Nature
469
:
323
335
.
40
Aichinger
M.
,
Wu
C.
,
Nedjic
J.
,
Klein
L.
.
2013
.
Macroautophagy substrates are loaded onto MHC class II of medullary thymic epithelial cells for central tolerance.
J. Exp. Med.
210
:
287
300
.
41
Schon
E. A.
,
Rizzuto
R.
,
Moraes
C. T.
,
Nakase
H.
,
Zeviani
M.
,
DiMauro
S.
.
1989
.
A direct repeat is a hotspot for large-scale deletion of human mitochondrial DNA.
Science
244
:
346
349
.
42
Samuels
D. C.
,
Schon
E. A.
,
Chinnery
P. F.
.
2004
.
Two direct repeats cause most human mtDNA deletions.
Trends Genet.
20
:
393
398
.
43
Khaidakov
M.
,
Siegel
E. R.
,
Shmookler Reis
R. J.
.
2006
.
Direct repeats in mitochondrial DNA and mammalian lifespan.
Mech. Ageing Dev.
127
:
808
812
.
44
Almasan
A.
,
Mishra
N. C.
.
1991
.
Recombination by sequence repeats with formation of suppressive or residual mitochondrial DNA in Neurospora.
Proc. Natl. Acad. Sci. USA
88
:
7684
7688
.
45
Robberson
D. L.
,
Clayton
D. A.
.
1972
.
Replication of mitochondrial DNA in mouse L cells and their thymidine kinase derivatives: displacement replication on a covalently-closed circular template.
Proc. Natl. Acad. Sci. USA
69
:
3810
3814
.
46
Yasukawa
T.
,
Reyes
A.
,
Cluett
T. J.
,
Yang
M. Y.
,
Bowmaker
M.
,
Jacobs
H. T.
,
Holt
I. J.
.
2006
.
Replication of vertebrate mitochondrial DNA entails transient ribonucleotide incorporation throughout the lagging strand.
EMBO J.
25
:
5358
5371
.
47
Holt
I. J.
,
Lorimer
H. E.
,
Jacobs
H. T.
.
2000
.
Coupled leading- and lagging-strand synthesis of mammalian mitochondrial DNA.
Cell
100
:
515
524
.
48
Krishnan
K. J.
,
Reeve
A. K.
,
Samuels
D. C.
,
Chinnery
P. F.
,
Blackwood
J. K.
,
Taylor
R. W.
,
Wanrooij
S.
,
Spelbrink
J. N.
,
Lightowlers
R. N.
,
Turnbull
D. M.
.
2008
.
What causes mitochondrial DNA deletions in human cells?
Nat. Genet.
40
:
275
279
.
49
Nakase
H.
,
Moraes
C. T.
,
Rizzuto
R.
,
Lombes
A.
,
DiMauro
S.
,
Schon
E. A.
.
1990
.
Transcription and translation of deleted mitochondrial genomes in Kearns-Sayre syndrome: implications for pathogenesis.
Am. J. Hum. Genet.
46
:
418
427
.
50
Nakada
K.
,
Inoue
K.
,
Ono
T.
,
Isobe
K.
,
Ogura
A.
,
Goto
Y. I.
,
Nonaka
I.
,
Hayashi
J. I.
.
2001
.
Inter-mitochondrial complementation: mitochondria-specific system preventing mice from expression of disease phenotypes by mutant mtDNA.
Nat. Med.
7
:
934
940
.
51
Gilkerson
R. W.
,
Schon
E. A.
,
Hernandez
E.
,
Davidson
M. M.
.
2008
.
Mitochondrial nucleoids maintain genetic autonomy but allow for functional complementation.
J. Cell Biol.
181
:
1117
1128
.
52
Chang
Y. F.
,
Imam
J. S.
,
Wilkinson
M. F.
.
2007
.
The nonsense-mediated decay RNA surveillance pathway.
Annu. Rev. Biochem.
76
:
51
74
.
53
Zhang
J.
,
Sun
X.
,
Qian
Y.
,
LaDuca
J. P.
,
Maquat
L. E.
.
1998
.
At least one intron is required for the nonsense-mediated decay of triosephosphate isomerase mRNA: a possible link between nuclear splicing and cytoplasmic translation.
Mol. Cell. Biol.
18
:
5272
5283
.
54
Thermann
R.
,
Neu-Yilik
G.
,
Deters
A.
,
Frede
U.
,
Wehr
K.
,
Hagemeier
C.
,
Hentze
M. W.
,
Kulozik
A. E.
.
1998
.
Binary specification of nonsense codons by splicing and cytoplasmic translation.
EMBO J.
17
:
3484
3494
.
55
Nagy
E.
,
Maquat
L. E.
.
1998
.
A rule for termination-codon position within intron-containing genes: when nonsense affects RNA abundance.
Trends Biochem. Sci.
23
:
198
199
.
56
Tiranti
V.
,
Corona
P.
,
Greco
M.
,
Taanman
J. W.
,
Carrara
F.
,
Lamantea
E.
,
Nijtmans
L.
,
Uziel
G.
,
Zeviani
M.
.
2000
.
A novel frameshift mutation of the mtDNA COIII gene leads to impaired assembly of cytochrome c oxidase in a patient affected by Leigh-like syndrome.
Hum. Mol. Genet.
9
:
2733
2742
.
57
Brower
R. C.
,
England
R.
,
Takeshita
T.
,
Kozlowski
S.
,
Margulies
D. H.
,
Berzofsky
J. A.
,
Delisi
C.
.
1994
.
Minimal requirements for peptide mediated activation of CD8+ CTL.
Mol. Immunol.
31
:
1285
1293
.
58
Sykulev
Y.
,
Cohen
R. J.
,
Eisen
H. N.
.
1995
.
The law of mass action governs antigen-stimulated cytolytic activity of CD8+ cytotoxic T lymphocytes.
Proc. Natl. Acad. Sci. USA
92
:
11990
11992
.
59
Sykulev
Y.
,
Joo
M.
,
Vturina
I.
,
Tsomides
T. J.
,
Eisen
H. N.
.
1996
.
Evidence that a single peptide-MHC complex on a target cell can elicit a cytolytic T cell response.
Immunity
4
:
565
571
.
60
Wallace
D. C.
2010
.
Mitochondrial DNA mutations in disease and aging.
Environ. Mol. Mutagen.
51
:
440
450
.
61
Dolfi
D. V.
,
Mansfield
K. D.
,
Polley
A. M.
,
Doyle
S. A.
,
Freeman
G. J.
,
Pircher
H.
,
Schmader
K. E.
,
Wherry
E. J.
.
2013
.
Increased T-bet is associated with senescence of influenza virus-specific CD8 T cells in aged humans.
J. Leukoc. Biol.
93
:
825
836
.
62
Zhang
Q.
,
Raoof
M.
,
Chen
Y.
,
Sumi
Y.
,
Sursal
T.
,
Junger
W.
,
Brohi
K.
,
Itagaki
K.
,
Hauser
C. J.
.
2010
.
Circulating mitochondrial DAMPs cause inflammatory responses to injury.
Nature
464
:
104
107
.
63
Krysko
D. V.
,
Agostinis
P.
,
Krysko
O.
,
Garg
A. D.
,
Bachert
C.
,
Lambrecht
B. N.
,
Vandenabeele
P.
.
2011
.
Emerging role of damage-associated molecular patterns derived from mitochondria in inflammation.
Trends Immunol.
32
:
157
164
.
64
Mirabella
M.
,
Di Giovanni
S.
,
Silvestri
G.
,
Tonali
P.
,
Servidei
S.
.
2000
.
Apoptosis in mitochondrial encephalomyopathies with mitochondrial DNA mutations: a potential pathogenic mechanism.
Brain
123
:
93
104
.
65
Asoh
S.
,
Mori
T.
,
Hayashi
J.
,
Ohta
S.
.
1996
.
Expression of the apoptosis-mediator Fas is enhanced by dysfunctional mitochondria.
J. Biochem.
120
:
600
607
.
66
Collins
L. V.
,
Hajizadeh
S.
,
Holme
E.
,
Jonsson
I. M.
,
Tarkowski
A.
.
2004
.
Endogenously oxidized mitochondrial DNA induces in vivo and in vitro inflammatory responses.
J. Leukoc. Biol.
75
:
995
1000
.
67
Hajizadeh
S.
,
DeGroot
J.
,
TeKoppele
J. M.
,
Tarkowski
A.
,
Collins
L. V.
.
2003
.
Extracellular mitochondrial DNA and oxidatively damaged DNA in synovial fluid of patients with rheumatoid arthritis.
Arthritis Res. Ther.
5
:
R234
R240
.

The authors have no financial conflicts of interest.

Supplementary data