UVB-Induced Tumor Heterogeneity Diminishes Immune Response in Melanoma.

Although clonal neo-antigen burden is associated with improved response to immune therapy, the functional basis for this remains unclear. Here we study this question in a novel controlled mouse melanoma model that enables us to explore the effects of intra-tumor heterogeneity (ITH) on tumor aggressiveness and immunity independent of tumor mutational burden. Induction of UVB-derived mutations yields highly aggressive tumors with decreased anti-tumor activity. However, single-cell-derived tumors with reduced ITH are swiftly rejected. Their rejection is accompanied by increased T cell reactivity and a less suppressive microenvironment. Using phylogenetic analyses and mixing experiments of single-cell clones, we dissect two characteristics of ITH: the number of clones forming the tumor and their clonal diversity. Our analysis of melanoma patient tumor data recapitulates our results in terms of overall survival and response to immune checkpoint therapy. These findings highlight the importance of clonal mutations in robust immune surveillance and the need to quantify patient ITH to determine the response to checkpoint blockade.

Urea Cycle Dysregulation Generates Clinically Relevant Genomic and Biochemical Signatures.

The urea cycle (UC) is the main pathway by which mammals dispose of waste nitrogen. We find that specific alterations in the expression of most UC enzymes occur in many tumors, leading to a general metabolic hallmark termed "UC dysregulation" (UCD). UCD elicits nitrogen diversion toward carbamoyl-phosphate synthetase2, aspartate transcarbamylase, and dihydrooratase (CAD) activation and enhances pyrimidine synthesis, resulting in detectable changes in nitrogen metabolites in both patient tumors and their bio-fluids. The accompanying excess of pyrimidine versus purine nucleotides results in a genomic signature consisting of transversion mutations at the DNA, RNA, and protein levels. This mutational bias is associated with increased numbers of hydrophobic tumor antigens and a better response to immune checkpoint inhibitors independent of mutational load. Taken together, our findings demonstrate that UCD is a common feature of tumors that profoundly affects carcinogenesis, mutagenesis, and immunotherapy response.

Trade-off between Transcriptome Plasticity and Genome Evolution in Cephalopods.

RNA editing, a post-transcriptional process, allows the diversification of proteomes beyond the genomic blueprint; however it is infrequently used among animals for this purpose. Recent reports suggesting increased levels of RNA editing in squids thus raise the question of the nature and effects of these events. We here show that RNA editing is particularly common in behaviorally sophisticated coleoid cephalopods, with tens of thousands of evolutionarily conserved sites. Editing is enriched in the nervous system, affecting molecules pertinent for excitability and neuronal morphology. The genomic sequence flanking editing sites is highly conserved, suggesting that the process confers a selective advantage. Due to the large number of sites, the surrounding conservation greatly reduces the number of mutations and genomic polymorphisms in protein-coding regions. This trade-off between genome evolution and transcriptome plasticity highlights the importance of RNA recoding as a strategy for diversifying proteins, particularly those associated with neural function. PAPERCLIP.

KPC1-mediated ubiquitination and proteasomal processing of NF-κB1 p105 to p50 restricts tumor growth.

NF-κB is a key transcriptional regulator involved in inflammation and cell proliferation, survival, and transformation. Several key steps in its activation are mediated by the ubiquitin (Ub) system. One uncharacterized step is limited proteasomal processing of the NF-κB1 precursor p105 to the p50 active subunit. Here, we identify KPC1 as the Ub ligase (E3) that binds to the ankyrin repeats domain of p105, ubiquitinates it, and mediates its processing both under basal conditions and following signaling. Overexpression of KPC1 inhibits tumor growth likely mediated via excessive generation of p50. Also, overabundance of p50 downregulates p65, suggesting that a p50-p50 homodimer may modulate transcription in place of the tumorigenic p50-p65. Transcript analysis reveals increased expression of genes associated with tumor-suppressive signals. Overall, KPC1 regulation of NF-κB1 processing appears to constitute an important balancing step among the stimulatory and inhibitory activities of the transcription factor in cell growth control.

Identification of bacteria-derived HLA-bound peptides in melanoma.

A variety of species of bacteria are known to colonize human tumours1-11, proliferate within them and modulate immune function, which ultimately affects the survival of patients with cancer and their responses to treatment12-14. However, it is not known whether antigens derived from intracellular bacteria are presented by the human leukocyte antigen class I and II (HLA-I and HLA-II, respectively) molecules of tumour cells, or whether such antigens elicit a tumour-infiltrating T cell immune response. Here we used 16S rRNA gene sequencing and HLA peptidomics to identify a peptide repertoire derived from intracellular bacteria that was presented on HLA-I and HLA-II molecules in melanoma tumours. Our analysis of 17 melanoma metastases (derived from 9 patients) revealed 248 and 35 unique HLA-I and HLA-II peptides, respectively, that were derived from 41 species of bacteria. We identified recurrent bacterial peptides in tumours from different patients, as well as in different tumours from the same patient. Our study reveals that peptides derived from intracellular bacteria can be presented by tumour cells and elicit immune reactivity, and thus provides insight into a mechanism by which bacteria influence activation of the immune system and responses to therapy.

The biogenesis of the immunopeptidome.

The immunopeptidome is the repertoire of peptides bound and presented by the MHC class I, class II, and non-classical molecules. The peptides are produced by the degradation of most cellular proteins, and in some cases, peptides are produced from extracellular proteins taken up by the cells. This review attempts to first describe some of its known and well-accepted concepts, and next, raise some questions about a few of the established dogmas in this field: The production of novel peptides by splicing is questioned, suggesting here that spliced peptides are extremely rare, if existent at all. The degree of the contribution to the immunopeptidome by degradation of cellular protein by the proteasome is doubted, therefore this review attempts to explain why it is likely that this contribution to the immunopeptidome is possibly overstated. The contribution of defective ribosome products (DRiPs) and non-canonical peptides to the immunopeptidome is noted and methods are suggested to quantify them. In addition, the common misconception that the MHC class II peptidome is mostly derived from extracellular proteins is noted, and corrected. It is stressed that the confirmation of sequence assignments of non-canonical and spliced peptides should rely on targeted mass spectrometry using spiking-in of heavy isotope-labeled peptides. Finally, the new methodologies and modern instrumentation currently available for high throughput kinetics and quantitative immunopeptidomics are described. These advanced methods open up new possibilities for utilizing the big data generated and taking a fresh look at the established dogmas and reevaluating them critically.

Publisher Correction: Actively personalized vaccination trial for newly diagnosed glioblastoma.

The additional author support information was erroneously omitted from the Supplementary Information. This has been corrected online.

Actively personalized vaccination trial for newly diagnosed glioblastoma.

Patients with glioblastoma currently do not sufficiently benefit from recent breakthroughs in cancer treatment that use checkpoint inhibitors1,2. For treatments using checkpoint inhibitors to be successful, a high mutational load and responses to neoepitopes are thought to be essential3. There is limited intratumoural infiltration of immune cells4 in glioblastoma and these tumours contain only 30-50 non-synonymous mutations5. Exploitation of the full repertoire of tumour antigens-that is, both unmutated antigens and neoepitopes-may offer more effective immunotherapies, especially for tumours with a low mutational load. Here, in the phase I trial GAPVAC-101 of the Glioma Actively Personalized Vaccine Consortium (GAPVAC), we integrated highly individualized vaccinations with both types of tumour antigens into standard care to optimally exploit the limited target space for patients with newly diagnosed glioblastoma. Fifteen patients with glioblastomas positive for human leukocyte antigen (HLA)-A*02:01 or HLA-A*24:02 were treated with a vaccine (APVAC1) derived from a premanufactured library of unmutated antigens followed by treatment with APVAC2, which preferentially targeted neoepitopes. Personalization was based on mutations and analyses of the transcriptomes and immunopeptidomes of the individual tumours. The GAPVAC approach was feasible and vaccines that had poly-ICLC (polyriboinosinic-polyribocytidylic acid-poly-L-lysine carboxymethylcellulose) and granulocyte-macrophage colony-stimulating factor as adjuvants displayed favourable safety and strong immunogenicity. Unmutated APVAC1 antigens elicited sustained responses of central memory CD8+ T cells. APVAC2 induced predominantly CD4+ T cell responses of T helper 1 type against predicted neoepitopes.

Phosphorylation of PP2Ac by PKC is a key regulatory step in the PP2A-switch-dependent AKT dephosphorylation that leads to apoptosis.

BACKGROUND: Although GqPCR activation often leads to cell survival by activating the PI3K/AKT pathway, it was previously shown that in several cell types AKT activity is reduced and leads to JNK activation and apoptosis. The mechanism of AKT inactivation in these cells involves an IGBP1-coupled PP2Ac switch that induces the dephosphorylation and inactivation of both PI3K and AKT. However, the machinery involved in the initiation of PP2A switch is not known. METHODS: We used phospho-mass spectrometry to identify the phosphorylation site of PP2Ac, and raised specific antibodies to follow the regulation of this phosphorylation. Other phosphorylations were monitored by commercial antibodies. In addition, we used coimmunoprecipitation and proximity ligation assays to follow protein-protein interactions. Apoptosis was detected by a TUNEL assay as well as PARP1 cleavage using SDS-PAGE and Western blotting. RESULTS: We identified Ser24 as a phosphorylation site in PP2Ac. The phosphorylation is mediated mainly by classical PKCs (PKCα and PKCß) but not by novel PKCs (PKCδ and PKCε). By replacing the phosphorylated residue with either unphosphorylatable or phosphomimetic residues (S24A and S24E), we found that this phosphorylation event is necessary and sufficient to mediate the PP2A switch, which ultimately induces AKT inactivation, and a robust JNK-dependent apoptosis. CONCLUSION: Our results show that the PP2A switch is induced by PKC-mediated phosphorylation of Ser24-PP2Ac and that this phosphorylation leads to apoptosis upon GqPCR induction of various cells. We propose that this mechanism may provide an unexpected way to treat some cancer types or problems in the endocrine machinery.

Are There Indeed Spliced Peptides in the Immunopeptidome?

The claims that a large fraction of the immunopeptidome is composed of spliced major histocompatibility complex (MHC) peptides have stirred significant excitement and raised controversy. Here, I suggest that there are likely no spliced peptides in the immunopeptidome, and if they exist at all, they are extremely rare. I base this claim on both biochemical and bioinformatics considerations. First, as a reactant in normal proteolytic reactions, water will compete with transpeptidation, which has been suggested as the mechanism of peptide splicing. The high mobility and abundance of water in aqueous solutions renders transpeptidation very inefficient and therefore unlikely to occur. Second, new studies have refuted the bioinformatics assignments to spliced peptides of most of the immunopeptidome MS data, suggesting that the correct assignments are likely other canonical, noncanonical, and post-translationally modified peptides. Therefore, I call for rigorous experimental methodology using heavy stable isotope peptides spiking into the immunoaffinity-purified mixtures of natural MHC peptides and analysis by the highly reliable targeted MS, to claim that MHC peptides are indeed spliced.

The Effect of Interferons on Presentation of Defective Ribosomal Products as HLA Peptides.

A subset of class I major histocompatibility complex (MHC)-bound peptides is produced from immature proteins that are rapidly degraded after synthesis. These defective ribosomal products (DRiPs) have been implicated in early alert of the immune system about impending infections. Interferons are important cytokines, produced in response to viral infection, that modulate cellular metabolism and gene expression patterns, increase the presentation of MHC molecules, and induce rapid degradation of proteins and cell-surface presentation of their derived MHC peptides, thereby contributing to the battle against pathogen infections. This study evaluated the role of interferons in the induction of rapid degradation of DRiPs to modulate the repertoire of DRiP-derived MHC peptides. Cultured human breast cancer cells were treated with interferons, and the rates of synthesis and degradation of cellular protein and their degradation products were determined by LC-MS/MS analysis, following the rates of incorporation of heavy stable isotope-labeled amino acids (dynamic stable isotope labeling by amino acids in cell culture, dynamic SILAC) at several time points after the interferon application. Large numbers of MHC peptides that incorporated the heavy amino acids faster than their source proteins indicated that DRiP peptides were abundant in the MHC peptidome; interferon treatment increased by about twofold their relative proportions in the peptidome. Such typical DRiP-derived MHC peptides were from the surplus subunits of the proteasome and ribosome, which are degraded because of the transition to immunoproteasomes and a new composition of ribosomes incorporating protein subunits that are induced by the interferon. We conclude that degradation of surplus subunits induced by the interferon is a major source for DRiP-MHC peptides, a phenomenon relevant to coping with viral infections, where a rapid presentation of MHC peptides derived from excess viral proteins may help alert the immune system about the impending infection.

Differential Modulation of the Phosphoproteome by the MAP Kinases Isoforms p38α and p38ß.

The p38 members of the mitogen-activated protein kinases (MAPKs) family mediate various cellular responses to stress conditions, inflammatory signals, and differentiation factors. They are constitutively active in chronic inflammatory diseases and some cancers. The differences between their transient effects in response to signals and the chronic effect in diseases are not known. The family is composed of four isoforms, of which p38α seems to be abnormally activated in diseases. p38α and p38ß are almost identical in sequence, structure, and biochemical and pharmacological properties, and the specific unique effects of each of them, if any, have not yet been revealed. This study aimed to reveal the specific effects induced by p38α and p38ß, both when transiently activated in response to stress and when chronically active. This was achieved via large-scale proteomics and phosphoproteomics analyses using stable isotope labeling of two experimental systems: one, mouse embryonic fibroblasts (MEFs) deficient in each of these p38 kinases and harboring either an empty vector or vectors expressing p38αWT, p38ßWT, or intrinsically active variants of these MAPKs; second, induction of transient stress by exposure of MEFs, p38α-/-, and p38ß-/- MEFs to anisomycin. Significant differences in the repertoire of the proteome and phosphoproteome between cells expressing active p38α and p38ß suggest distinct roles for each kinase. Interestingly, in both cases, the constitutive activation induced adaptations of the cells to the chronic activity so that known substrates of p38 were downregulated. Within the dramatic effect of p38s on the proteome and phosphoproteome, some interesting affected phosphorylation sites were those found in cancer-associated p53 and Hspb1 (HSP27) proteins and in cytoskeleton-associated proteins. Among these, was the stronger direct phosphorylation by p38α of p53-Ser309, which was validated on the Ser315 in human p53. In summary, this study sheds new light on the differences between chronic and transient p38α and p38ß signaling and on the specific targets of these two kinases.

Abundance, Betweenness Centrality, Hydrophobicity, and Isoelectric Points Are Relevant Factors in the Processing of Parental Proteins of the HLA Class II Ligandome.

Adaptive cellular and humoral immune responses to infectious agents require previous recognition of pathogenic peptides bound to human leukocyte antigen (HLA) class II molecules exposed on the surface of the professional antigen-presenting cells. Knowledge of how these peptide ligands are generated is essential to understand the basis for CD4+ T-cell-mediated immunity and tolerance. In this study, a high-throughput mass spectrometry analysis was used to identify more than 16,000 cell peptides bound to several HLA-DR and -DP class II molecules isolated from large amounts of uninfected and virus-infected human cells (ProteomeXchange accession: PXD028006). The analysis of the 1808 parental proteins containing HLA class II ligands revealed that these cell proteins were more acidic, abundant, and highly connected but less hydrophilic than non-parental proteomes. Therefore, the percentage of acidic residues was increased and hydroxyl and polar residues were decreased in the parental proteins for the HLA class II ligandomes versus the non-parental proteomes. This definition of the properties shared by parental proteins that constitute the source of the HLA class II ligandomes can serve as the basis for the development of bioinformatics tools to predict proteins that are most likely recognized by the immune system through the CD4+ helper T lymphocytes in both autoimmunity and infection.

Modulation of Natural HLA-B*27:05 Ligandome by Ankylosing Spondylitis-associated Endoplasmic Reticulum Aminopeptidase 2 (ERAP2).

The HLA-B*27:05 allele and the endoplasmic reticulum-resident aminopeptidases are strongly associated with AS, a chronic inflammatory spondyloarthropathy. This study examined the effect of ERAP2 in the generation of the natural HLA-B*27:05 ligandome in live cells. Complexes of HLA-B*27:05-bound peptide pools were isolated from human ERAP2-edited cell clones, and the peptides were identified using high-throughput mass spectrometry analyses. The relative abundance of a thousand ligands was established by quantitative tandem mass spectrometry and bioinformatics analysis. The residue frequencies at different peptide position, identified in the presence or absence of ERAP2, determined structural features of ligands and their interactions with specific pockets of the antigen-binding site of the HLA-B*27:05 molecule. Sequence alignment of ligands identified with species of bacteria associated with HLA-B*27-dependent reactive arthritis was performed. In the absence of ERAP2, peptides with N-terminal basic residues and minority canonical P2 residues are enriched in the natural ligandome. Further, alterations of residue frequencies and hydrophobicity profile at P3, P7, and PΩ positions were detected. In addition, several ERAP2-dependent cellular peptides were highly similar to protein sequences of arthritogenic bacteria, including one human HLA-B*27:05 ligand fully conserved in a protein from Campylobacter jejuni These findings highlight the pathogenic role of this aminopeptidase in the triggering of AS autoimmune disease.

Identification of Tumor Antigens in the HLA Peptidome of Patient-derived Xenograft Tumors in Mouse.

Personalized cancer immunotherapy targeting patient-specific cancer/testis antigens (CTA) and neoantigens may benefit from large-scale tumor human leukocyte antigen (HLA) peptidome (immunopeptidome) analysis, which aims to accurately identify antigens presented by tumor cells. Although significant efforts have been invested in analyzing the HLA peptidomes of fresh tumors, it is often impossible to obtain sufficient volumes of tumor tissues for comprehensive HLA peptidome characterization. This work attempted to overcome some of these obstacles by using patient-derived xenograft tumors (PDX) in mice as the tissue sources for HLA peptidome analysis. PDX tumors provide a proxy for the expansion of the patient tumor by re-grafting them through several passages to immune-compromised mice. The HLA peptidomes of human biopsies were compared with those derived from PDX tumors. Larger HLA peptidomes were obtained from the significantly larger PDX tumors as compared with the patient biopsies. The HLA peptidomes of different PDX tumors derived from the same source tumor biopsy were very reproducible, even following subsequent passages to new naïve mice. Many CTA-derived HLA peptides were discovered, as well as several potential neoantigens/variant sequences. Taken together, the use of PDX tumors for HLA peptidome analysis serves as a highly expandable and stable source of reproducible and authentic peptidomes, opening up new opportunities for defining large HLA peptidomes when only small tumor biopsies are available. This approach provides a large source for tumor antigens identification, potentially useful for personalized immunotherapy.

Substantial Influence of ERAP2 on the HLA-B*40:02 Peptidome: Implications for HLA-B*27-Negative Ankylosing Spondylitis.

HLA-B*40:02 is one of a few major histocompatibility complex class I (MHC-I) molecules associated with ankylosing spondylitis (AS) independently of HLA-B*27. The endoplasmic reticulum aminopeptidase 2 (ERAP2), an enzyme that process MHC-I ligands and preferentially trims N-terminal basic residues, is also a risk factor for this disease. Like HLA-B*27 and other AS-associated MHC-I molecules, HLA-B*40:02 binds a relatively high percentage of peptides with ERAP2-susceptible residues. In this study, the effects of ERAP2 depletion on the HLA-B*40:02 peptidome were analyzed. ERAP2 protein expression was knocked out by CRISPR in the transfectant cell line C1R-B*40:02, and the differences between the peptidomes from the wild-type and ERAP2-KO cells were determined by label-free quantitative comparisons. The qualitative changes dependent on ERAP2 affected about 5% of the peptidome, but quantitative changes in peptide amounts were much more substantial, reflecting a significant influence of this enzyme on the generation/destruction balance of HLA-B*40:02 ligands. As in HLA-B*27, a major effect was on the frequencies of N-terminal residues. In this position, basic and small residues were increased, and aliphatic/aromatic ones decreased in the ERAP2 knockout. Other peptide positions were also affected. Because most of the non-B*27 MHC-I molecules associated with AS risk bind a relatively high percentage of peptides with N-terminal basic residues, we hypothesize that the non-epistatic association of ERAP2 with AS might be related to the processing of peptides with these residues, thus affecting the peptidomes of AS-associated MHC-I molecules.

Redundancy and Complementarity between ERAP1 and ERAP2 Revealed by their Effects on the Behcet's Disease-associated HLA-B*51 Peptidome.

The endoplasmic reticulum aminopeptidases ERAP1 and ERAP2 trim peptides to be loaded onto HLA molecules, including the main risk factor for Behçet's disease HLA-B*51. ERAP1 is also a risk factor among HLA-B*51-positive individuals, whereas no association is known with ERAP2. This study addressed the mutual relationships between both enzymes in the processing of an HLA-bound peptidome, interrogating their differential association with Behçet's disease. CRISPR/Cas9 was used to generate knock outs of ERAP1, ERAP2 or both from transfectant 721.221-HLA-B*51:01 cells. The surface expression of HLA-B*51 was reduced in all cases. The effects of depleting each or both enzymes on the B*51:01 peptidome were analyzed by quantitative label-free mass spectrometry. Substantial quantitative alterations of peptide length, subpeptidome balance, N-terminal residue usage, affinity and presentation of noncanonical ligands were observed. These effects were often different in the presence or absence of the other enzyme, revealing their mutual dependence. In the absence of ERAP1, ERAP2 showed similar and significant processing of B*51:01 ligands, indicating functional redundancy. The high overlap between the peptidomes of wildtype and double KO cells indicates that a large majority of B*51:01 ligands are present in the ER even in the absence of ERAP1/ERAP2. These results indicate that both enzymes have distinct, but complementary and partially redundant effects on the B*51:01 peptidome, leading to its optimization and maximal surface expression. The distinct effects of both enzymes on the HLA-B*51 peptidome provide a basis for their differential association with Behçet's disease and suggest a pathogenetic role of the B*51:01 peptidome.

Identification of Tumor Antigens Among the HLA Peptidomes of Glioblastoma Tumors and Plasma.

Glioblastoma multiforme (GBM) is the most aggressive brain tumor with poor prognosis to most patients. Immunotherapy of GBM is a potentially beneficial treatment option, whose optimal implementation may depend on familiarity with tumor specific antigens, presented as HLA peptides by the GBM cells. Further, early detection of GBM, such as by a routine blood test, may improve survival, even with the current treatment modalities. This study includes large-scale analyses of the HLA peptidome (immunopeptidome) of the plasma-soluble HLA molecules (sHLA) of 142 plasma samples, and the membranal HLA of GBM tumors of 10 of these patients' tumor samples. Tumor samples were fresh-frozen immediately after surgery and the plasma samples were collected before, and at multiple visits after surgery. In total, this HLA peptidome analysis involved 52 different HLA allotypes and resulted in the identification of more than 35,000 different HLA peptides. Strong correlations were observed in the signal intensities and in the repertoires of identified peptides between the tumors and plasma-soluble HLA peptidomes of the individual patients, whereas low correlations were observed between these HLA peptidomes and the tumors' proteomes. HLA peptides derived from Cancer/Testis Antigens (CTAs) were selected based on their presence among the HLA peptidomes of the patients and absence of expression of their source genes from any healthy and essential human tissues, except from immune-privileged sites. Additionally, peptides were selected as potential biomarkers if their levels in the plasma-sHLA peptidome were significantly reduced after the removal of tumor mass. The CTAs identified among the analyzed HLA peptidomes provide new opportunities for personalized immunotherapy and for early diagnosis of GBM.

The direction of protein entry into the proteasome determines the variety of products and depends on the force needed to unfold its two termini.

Poorly structured domains in proteins enhance their susceptibility to proteasomal degradation. To learn whether the presence of such a domain near either end of a protein determines its direction of entry into the proteasome, directional translocation was enforced on several proteasome substrates. Using archaeal PAN-20S complexes, mammalian 26S proteasomes, and cultured cells, we identified proteins that are degraded exclusively from either the C or N terminus and some showing no directional preference. This property results from interactions of the substrate's termini with the regulatory ATPase and could be predicted based on the calculated relative stabilities of the N and C termini. Surprisingly, the direction of entry into the proteasome affected markedly the spectrum of peptides released and consequently influenced the efficiency of MHC class I presentation. Thus, easily unfolded termini are translocated first, and the direction of translocation influences the peptides generated and presented to the immune system.

Ranking the Contribution of Ankylosing Spondylitis-associated Endoplasmic Reticulum Aminopeptidase 1 (ERAP1) Polymorphisms to Shaping the HLA-B*27 Peptidome.

The Endoplasmic reticulum aminopeptidase I (ERAP1) trims peptides to their optimal size for binding to Major Histocompatibility Complex class I proteins. The natural polymorphism of this enzyme is associated with ankylosing spondylitis (AS) in epistasis with the major risk factor for this disease, HLA-B*27, suggesting a direct relationship between AS and HLA-B*27-bound peptides. Three polymorphisms that affect peptide trimming protect from AS: K528R, D575N/R725Q, and Q730E. We characterized and ranked the effects of each mutation, and their various combinations, by quantitative comparisons of the HLA-B*27 peptidomes from cells expressing distinct ERAP1 variants. Five features were examined: peptide length, N-terminal flanking residues, N-terminal residues of the natural ligands, internal sequences and affinity for B*27:05. Polymorphism at residue 528 showed the largest influence, affecting all five features regardless of peptide length. D575N/R725Q showed a much smaller effect. Yet, when co-occurring with K528R, it further decreased ERAP1 activity. Polymorphism at residue 730 showed a significant influence on peptide length, because of distinct effects on trimming of nonamers compared with longer peptides. Accordingly, multiple features were affected by the Q730E mutation in a length-dependent way. The alterations induced in the B*27:05 peptidome by natural ERAP1 variants with different K528R/Q730E combinations reflected separate and additive effects of both mutations. Thus, the influence of ERAP1 on HLA-B*27 is very diverse at the population level, because of the multiplicity and complexity of ERAP1 variants, and to the distinct effects of their co-occurring polymorphisms, leading to significant modulation of disease risk among HLA-B*27-positive individuals.