Functional interaction of the ankylosing spondylitis-associated endoplasmic reticulum aminopeptidase 1 polymorphism and HLA-B27 in vivo

Ankylosing spondylitis: From pathogenesis to therapy

Ankylosing spondylitis (AS) is an autoimmune rheumatic disease that primarily affects the axial joints, with its etiology complex and still not fully understood. The unknown pathogenesis of AS limits the development of treatment strategies, so keeping up-to-date with the current research on AS can help in searching for potential therapeutic targets. In addition to the classic HLA-B27 genetic susceptibility and Th17-related inflammatory signals, increasing research is focusing on the influence of autoantigen-centered autoimmune responses and bone stromal cells on the onset of AS. Autoantigens derived from gut microbiota and preferential TCR both exacerbate the autoimmune response in patients with AS. Furthermore, dysregulated bone metabolism also promotes pathological new bone formation in AS. Current treatments approved for AS almost focus on the management of inflammation with inconsistent treatment results due to the heterogeneity of patients. In this review, we systematically summarized various pathogenesis and management of AS, meanwhile discussed the underlying risk factors and potential therapeutic targets.

ERAP1 Gene Variants and Haplotypes Associated With Psoriasis Vulgaris of Han Chinese in Inner Mongolia

BACKGROUND

This study aimed to investigate the association between genetic variants of ERAP1 (OMIM: 606832) and psoriasis vulgaris (PsV) susceptibility in Inner Mongolia Han nationality.

METHODS

For primary screening, the subjects included 142 PsV cases and 100 healthy controls without psoriasis. The 27 exons of ERAP1 gene were sequenced to screen significant genetic variants. For the validation study, the subjects included 1030 PsV cases and 965 healthy controls. A total of 18 mutations were detected for genetic variants of significance in primary screening and previously reported genetic variants.

RESULTS

In primary screening stage, 13 genetic variants of ERAP1 showed an association with psoriasis. A total of 18 genetic variants were typed for the validation, and 12 genetic variants were associated with PsV in Inner Mongolia Han population. Stratified analysis showed significant differences in the allele frequencies of 8 ERAP1 genetic variants in cases with positive family history, and significant differences in allele frequencies among 9 ERAP1 genetic variants in patients with negative family history. A risk haplotype (TCCCTCCAGACC) was significantly associated with PsV, and the most risk haplotype was E730/K528/R127/E56.

CONCLUSION

ERAP1 gene mutation may be associated with PsV and HLA-C*06:02 in Han nationality in Inner Mongolia. A risk haplotype of four-nonsynonymous mutation (E730/K528/R127/E56) is associated with PsV.

FASTMAP-a flexible and scalable immunopeptidomics pipeline for HLA- and antigen-specific T-cell epitope mapping based on artificial antigen-presenting cells

The study of peptide repertoires presented by major histocompatibility complex (MHC) molecules and the identification of potential T-cell epitopes contribute to a multitude of immunopeptidome-based treatment approaches. Epitope mapping is essential for the development of promising epitope-based approaches in vaccination as well as for innovative therapeutics for autoimmune diseases, infectious diseases, and cancer. It also plays a critical role in the immunogenicity assessment of protein therapeutics with regard to safety and efficacy concerns. The main challenge emerges from the highly polymorphic nature of the human leukocyte antigen (HLA) molecules leading to the requirement of a peptide mapping strategy for a single HLA allele. As many autoimmune diseases are linked to at least one specific antigen, we established FASTMAP, an innovative strategy to transiently co-transfect a single HLA allele combined with a disease-specific antigen into a human cell line. This approach allows the specific identification of HLA-bound peptides using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Using FASTMAP, we found a comparable spectrum of endogenous peptides presented by the most frequently expressed HLA alleles in the world's population compared to what has been described in literature. To ensure a reliable peptide mapping workflow, we combined the HLA alleles with well-known human model antigens like coagulation factor VIII, acetylcholine receptor subunit alpha, protein structures of the SARS-CoV-2 virus, and myelin basic protein. Using these model antigens, we have been able to identify a broad range of peptides that are in line with already published and in silico predicted T-cell epitopes of the specific HLA/model antigen combination. The transient co-expression of a single affinity-tagged MHC molecule combined with a disease-specific antigen in a human cell line in our FASTMAP pipeline provides the opportunity to identify potential T-cell epitopes/endogenously processed MHC-bound peptides in a very cost-effective, fast, and customizable system with high-throughput potential.

Non-mutational neoantigens in disease

The ability of mammals to mount adaptive immune responses culminating with the establishment of immunological memory is predicated on the ability of the mature T cell repertoire to recognize antigenic peptides presented by syngeneic MHC class I and II molecules. Although it is widely believed that mature T cells are highly skewed towards the recognition of antigenic peptides originating from genetically diverse (for example, foreign or mutated) protein-coding regions, preclinical and clinical data rather demonstrate that novel antigenic determinants efficiently recognized by mature T cells can emerge from a variety of non-mutational mechanisms. In this Review, we describe various mechanisms that underlie the formation of bona fide non-mutational neoantigens, such as epitope mimicry, upregulation of cryptic epitopes, usage of non-canonical initiation codons, alternative RNA splicing, and defective ribosomal RNA processing, as well as both enzymatic and non-enzymatic post-translational protein modifications. Moreover, we discuss the implications of the immune recognition of non-mutational neoantigens for human disease.

A large meta-analysis identifies genes associated with anterior uveitis

Anterior Uveitis (AU) is the inflammation of the anterior part of the eye, the iris and ciliary body and is strongly associated with HLA-B*27. We report AU exome sequencing results from eight independent cohorts consisting of 3,850 cases and 916,549 controls. We identify common genome-wide significant loci in HLA-B (OR = 3.37, p = 1.03e-196) and ERAP1 (OR = 0.86, p = 1.1e-08), and find IPMK (OR = 9.4, p = 4.42e-09) and IDO2 (OR = 3.61, p = 6.16e-08) as genome-wide significant genes based on the burden of rare coding variants. Dividing the cohort into HLA-B*27 positive and negative individuals, we find ERAP1 haplotype is strongly protective only for B*27-positive AU (OR = 0.73, p = 5.2e-10). Investigation of B*27-negative AU identifies a common signal near HLA-DPB1 (rs3117230, OR = 1.26, p = 2.7e-08), risk genes IPMK and IDO2, and several additional candidate risk genes, including ADGFR5, STXBP2, and ACHE. Taken together, we decipher the genetics underlying B*27-positive and -negative AU and identify rare and common genetic signals for both subtypes of disease.

The ER Aminopeptidases, ERAP1 and ERAP2, synergize to self-modulate their respective activities

Introduction

Critical steps in Major Histocompatibility Complex Class I (MHC-I) antigen presentation occur in the endoplasmic reticulum (ER). In general, peptides that enter the ER are longer than the optimal length for MHC-I binding. The final trimming of MHC-I epitopes is performed by two related aminopeptidases, ERAP1 and ERAP2 in humans that possess unique and complementary substrate trimming specificities. While ERAP1 efficiently trims peptides longer than 9 residues, ERAP2 preferentially trims peptides shorter than 9 residues.

Materials and Methods

Using a combination of biochemical and proteomic studies followed by biological verification.

Results

We demonstrate that the optimal ligands for either enzyme act as inhibitors of the other enzyme. Specifically, the presence of octamers reduced the trimming of long peptides by ERAP1, while peptides longer than nonomers inhibit ERAP2 activity.

Discussion

We propose a mechanism for how ERAP1 and ERAP2 synergize to modulate their respective activities and shape the MHC-I peptidome by generating optimal peptides for presentation.

Functional ERAP1 Variants Distinctively Associate with Ankylosing Spondylitis Susceptibility under the Influence of HLA-B27 in Taiwanese

Epistasis of ERAP1 single nucleotide variations (SNVs) and HLA-B27 has been linked to ankylosing spondylitis susceptibility (AS). The current study examined how prevalent ERAP1 allelic variants (SNV haplotypes) in Taiwan affect ERAP1 functions and AS susceptibility in the presence or absence of HLA-B27. Sanger sequencing was used to discover all ERAP1 coding SNVs and common allelic variants in Taiwanese full-length cDNAs from 45 human patients. For the genetic association investigation, TaqMan genotyping assays were utilized to establish the genotypes of ERAP1 SNVs in 863 AS patients and 1438 healthy controls. Ex vivo biological analysis of peripheral blood mononuclear cells from homozygous donors of two common-risk ERAP1 allelic variants was performed. Two common-risk ERAP1 allelic variants were also cloned and functionally studied. In Taiwanese, eleven frequent ERAP1 SNVs and six major ERAP1 allelic variants were discovered. We discovered that in Taiwanese, the most prevalent ERAP1-001 variant with 56E, 127R, 276I, 349M, 528K, 575D, 725R, and 730Q interacting with HLA-B27 significantly contributed to the development of AS. In HLA-B27 negative group, however, the second most prevalent ERAP1-002 variant with 56E, 127P, 276M, 349M, 528R, 575D, 725R, and 730E was substantially related with an increased risk of AS. Ex vivo and in vitro research demonstrated that ERAP1 allelic variants have a significant impact on ERAP1 functions, suggesting that ERAP1 plays a role in the development of AS. In an HLA-B27-dependent manner, common ERAP1 allelic variants are related with AS susceptibility.

T Cell Epitope Discovery in the Context of Distinct and Unique Indigenous HLA Profiles

CD8+ T cells are a pivotal part of the immune response to viruses, playing a key role in disease outcome and providing long-lasting immunity to conserved pathogen epitopes. Understanding CD8+ T cell immunity in humans is complex due to CD8+ T cell restriction by highly polymorphic Human Leukocyte Antigen (HLA) proteins, requiring T cell epitopes to be defined for different HLA allotypes across different ethnicities. Here we evaluate strategies that have been developed to facilitate epitope identification and study immunogenic T cell responses. We describe an immunopeptidomics approach to sequence HLA-bound peptides presented on virus-infected cells by liquid chromatography with tandem mass spectrometry (LC-MS/MS). Using antigen presenting cell lines that stably express the HLA alleles characteristic of Indigenous Australians, this approach has been successfully used to comprehensively identify influenza-specific CD8+ T cell epitopes restricted by HLA allotypes predominant in Indigenous Australians, including HLA-A*24:02 and HLA-A*11:01. This is an essential step in ensuring high vaccine coverage and efficacy in Indigenous populations globally, known to be at high risk from influenza disease and other respiratory infections.

A proposed HLA-B*27 screening method for ankylosing spondylitis detection based on tag-single nucleotide polymorphisms: a preliminary study

BACKGROUND

Ankylosing spondylitis (AS) is a type of inflammatory arthritis that affects primarily the spine. There is a strong association of the HLA-B*27 allele with AS pathogenesis, but recent studies have demonstrated the participation of ERAP1 gene in the genetic susceptibility. The aim of this study was to determine whether HLA-B tag-single nucleotide polymorphisms (SNPs) and ERAP1-related genetic variations associated with AS have equal or similarly performance in patients´ screening compared to HLA-B*27 standard genotyping in Mexican population.

METHODS AND RESULTS

Genomic DNA from patients with AS and population-based controls from Mexico City was analyzed for five single nucleotide polymorphisms (SNPs): rs4349859, rs13202464, rs116488202, tagging HLA-B*27; and rs30187 and rs27044 in ERAP1 gene. TaqMan genotype assay method was used for SNPs genotyping. We found a significant association between AS and the heterozygote genotypes and minor alleles of the HLA-B*27 tag-SNPs, as well as for their haplotypes. With respect to ERAP1 polymorphisms, no significant associations were observed (p > 0.05). The sensitivity and specificity analysis showed values of 0.96 and 1.00 for the rs4349859 SNP, and 0.96 and 0.94 for the rs116488202 SNP, respectively, in detecting HLA-B*27 compared to the B27 test as the gold standard.

CONCLUSIONS

HLA-B*27 tag-SNPs are associated with AS susceptibility; furthermore, the rs4349859 SNP by its own have an outstanding performance in detecting HLA-B*27 and therefore can be proposed as screening marker in the identification of HLA-B*27 in our population.

Genetic contribution of endoplasmic reticulum aminopeptidase 1 polymorphisms to liver fibrosis progression in patients with HCV infection

The endoplasmic reticulum aminopeptidase ERAP1 regulates innate and adaptive immune responses, trimming peptides and loading onto HLA class I molecules. Coding single nucleotide polymorphisms within ERAP1 are associated with autoimmune diseases, viral infections, and cancer development. Our purpose was to analyze the influence of ERAP1 variants on fibrogenesis in hepatitis C virus (HCV)-infected patients. A range of ERAP1 polymorphisms were genotyped in 722 unrelated Caucasian patients diagnosed with chronic HCV from two Spanish cohorts. Patients were classified according to their fibrosis stage. Paraffin-embedded tissue microarrays were constructed to assess ERAP1 expression (HCV = 38; alcoholic = 20) by immunohistochemistry. A statistical algorithm was applied to derive a fibrogenesis prediction model. The ERAP1 variants rs30187/T (K528, p_c < 0.001) and rs27044/G (Q730, p_c < 0.001) were related with severe fibrosis. These results were validated in the two independent cohorts. Furthermore, patients with the rs30187/T allele had stronger ERAP1 protein expression than those with the rs30187/C (p < 0.05). The statistical model showed that patients with rs30187 C/T and T/T genotypes took 15.58 years (median) to develop advanced fibrosis, but this value was 32.08 years in patients carrying C/C genotype (p < 0.005). ERAP1 variants may influence the clinical course of fibrogenesis in HCV-infected patients. These polymorphisms could be exploited as constitutive new markers of fibrosis evolution. The results highlight the possibility of using modulators of ERAP1 to generate a protective immune response against chronic HCV infection. KEY MESSAGES: What is known Several ERAP1 polymorphisms are associated with autoimmune diseases and cancer. ERAP1 trims peptides to HLA class I presentation. What is new here ERAP1 polymorphisms are associated with fibrogenesis. The ERAP1 polymorphisms genotype could help us in clinical management of patients. Potential translational impact The use of modulators of ERAP1 could generate a protective response depending on SNPs.

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.

Endoplasmic reticulum aminopeptidase 2 gene single nucleotide polymorphisms in association with susceptibility to ankylosing spondylitis in an Iranian population

BACKGROUND

Ankylosing spondylitis (AS) is a chronic autoimmune disease, in which genetic polymorphisms are critically important in establishing inflammatory state. Endoplasmic reticulum aminopeptidase (ERAP) 2 gene has been implied to be involved in AS etiopathogenesis. The current study evaluated the association of ERAP2 gene single nucleotide polymorphisms (SNPs) with susceptibility to AS in an Iranian population.

METHODS

Two hundred and forty AS patients and 240 healthy individuals were recruited. DNA extraction was performed from whole blood samples and RNA content was isolated from peripheral blood mononuclear cells (PBMCs). Real-time allelic discrimination approach was exerted to genotype all subjects for rs2910686, rs2248374, and rs2549782 SNPs. After cDNA synthesis, mRNA expression of cytokines was determined. Enzyme-linked immunosorbent assay (ELISA) was exerted to evaluate the cytokine levels in serum of participants.

RESULTS

None of the SNPs were associated with AS risk in the whole population. However, allele and heterozygote genotype of rs2910686 SNP were associated significantly with higher risk of AS in Human leukocyte antigen (HLA)-B27 positive group. mRNA expression and serum concentrations of interleukin (IL)-17A, IL-23, interferon (IFN)-γ, and tumor necrosis factor (TNF)-α was increased in AS patients compared with controls. Nonetheless, mRNA expression and serum levels of cytokines was not significantly different among HLA-B27 positive AS patients with different three genotypes for rs2910686 SNP.

CONCLUSIONS

AlthoughERAP2 gene rs2910686 polymorphism was significantly associated with increased risk of AS susceptibility, it might not be involved in regulation of the inflammatory cytokines during AS pathogenesis.

The roles of ERAP1 and ERAP2 in autoimmunity and cancer immunity: New insights and perspective

Autoimmunity and cancer affect millions worldwide and both, in principal, result from dysregulated immune responses. There are many well-known molecules involved in immunological process playing as a double-edged sword, by which associating autoimmune diseases and cancer. In this regard, Endoplasmic reticulum aminopeptidases (ERAP) 1, which belongs to the M1 family of aminopeptidases, plays a central role as a "molecular ruler", proteolyzing of N-terminal of the antigenic peptides before their loading onto HLA-I molecules for antigen presentation in the Endoplasmic Reticulum (ER). Several genome-wide association studies (GWAS) highlighted the significance of ERAP1 and ERAP2 in autoimmune diseases, including Ankylosing spondylitis, Psoriasis, Bechet's disease, and Birdshot chorioretinopathy, as well as in cancers. The expression of ERAP1/2 is mostly altered in different cancers compared to normal cells, but how this affects anti-cancer immune responses and cancer growth has been little explored. Recent studies on the immunological outcomes and the catalytic functions of ERAP1 and ERAP2 have provided a better understanding of their potential pathogenetic role in autoimmunity and cancer. In this review, we summarize the role of ERAP1 and ERAP2 in the autoimmune diseases and cancer immunity based on the recent advances in GWAS studies.

Variant in ERAP1 promoter region is associated with low expression in a patient with a Behçet-like MHC-I-opathy

Behçet disease (BD) is an immune-mediated disease. The cause of BD remains unknown, but the existence of multiple pathological pathways is suspected, including different genetic factors. Polymorphisms in ERAP1 gene have been associated with an increased risk of BD. However, while current BD-associated ERAP1 variants are suggested to contribute to disease by altering the activity of the encoded protein, there is no knowledge of variants that alter the expression level of ERAP1, despite previous associations between ERAP1 expression and BD. Here, we used whole-exome sequencing of a patient with a Behçet-like MHC-I-opathy to identify that the patient, unlike its healthy parents, was homozygous for a rare 1-bp deletion, rs140416843, in the promoter region of ERAP1. rs140416843 has not previously been associated with disease, but is linked to ERAP1 haplotype Hap10 which is associated with BD. The expression of ERAP1 by both RT-qPCR and RNA sequencing showed that ERAP1 mRNA expression correlated with the zygosity for the identified deletion and was decreased in comparison to a healthy cohort. In conclusion, we diagnosed the patient as having BD, and hypothesize that rs140416843-mediated changes in ERAP1 expression play a causative role in BD and that this risk factor is contributing to the association between Hap10 and BD. This is the first report to identify a variant that may cause BD by altering the expression of ERAP1, and our findings suggest that downregulation of ERAP1 expression can serve as a diagnostic marker for BD.

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.

Advances in nanomedicine for the treatment of ankylosing spondylitis

Ankylosing spondylitis (AS) is a complex disease characterized by inflammation and ankylosis primarily at the cartilage-bone interface. The disease is more common in young males and risk factors include both genetic and environmental. While the pathogenesis of AS is not completely understood, it is thought to be an immune-mediated disease involving inflammatory cellular infiltrates, and human leukocyte antigen-B27. Currently, there is no specific diagnostic technique available for this disease; therefore conventional diagnostic approaches such as clinical symptoms, laboratory tests and imaging techniques are used. There are various review papers that have been published on conventional treatment approaches, and in this review work, we focus on the more promising nanomedicine-based treatment modalities to move this field forward.

Inhibitors of ER Aminopeptidase 1 and 2: From Design to Clinical Application

Endoplasmic Reticulum aminopeptidase 1 and 2 are two homologous enzymes that help generate peptide ligands for presentation by Major Histocompatibility Class I molecules. Their enzymatic activity influences the antigenic peptide repertoire and indirectly controls adaptive immune responses. Accumulating evidence suggests that these two enzymes are tractable targets for the regulation of immune responses with possible applications ranging from cancer immunotherapy to treating inflammatory autoimmune diseases. Here, we review the state-of-the-art in the development of inhibitors of ERAP1 and ERAP2 as well as their potential and limitations for clinical applications.

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.

ERAP1 allotypes shape the epitope repertoire of virus-specific CD8+ T cell responses in acute hepatitis C virus infection

BACKGROUND & AIMS

Endoplasmic reticulum aminopeptidase 1 (ERAP1) polymorphisms are linked with human leukocyte antigen (HLA) class I-associated autoinflammatory disorders, including ankylosing spondylitis and Behçet's disease. Disease-associated ERAP1 allotypes exhibit distinct functional properties, but it remains unclear how differential peptide trimming in vivo affects the repertoire of epitopes presented to CD8⁺ T cells. The aim of this study was to determine the impact of ERAP1 allotypes on the virus-specific CD8⁺ T cell epitope repertoire in an HLA-B*27:05⁺ individual with acute hepatitis C virus (HCV) infection.

METHODS

We performed genetic and functional analyses of ERAP1 allotypes and characterized the HCV-specific CD8⁺ T cell repertoire at the level of fine epitope specificity and HLA class I restriction, in a patient who had acquired an HCV genotype 1a infection through a needle-stick injury.

RESULTS

Two hypoactive allotypic variants of ERAP1 were identified in an individual with acute HCV infection. The associated repertoire of virus-derived epitopes recognized by CD8⁺ T cells was uncommon in a couple of respects. Firstly, reactivity was directed away from classically immunodominant epitopes, preferentially targeting either novel or subdominant epitopes. Secondly, reactivity was biased towards longer epitopes (10-11-mers). Despite the patient exhibiting favorable prognostic indicators, these atypical immune responses failed to clear the virus and the patient developed persistent low-level infection with HCV.

CONCLUSIONS

ERAP1 allotypes modify the virus-specific CD8⁺ T cell epitope repertoire in vivo, leading to altered immunodominance patterns that may contribute to the failure of antiviral immunity after infection with HCV.

LAY SUMMARY

Endoplasmic reticulum aminopeptidase 1 (ERAP1) plays a key role in antigen presentation. Genetic variants of ERAP1 (leading to distinct allotypes) are linked with specific autoinflammatory disorders, such as ankylosing spondylitis and Behçet's disease. We found that ERAP1 allotypes modified the repertoire of virus-specific CD8⁺ T cell epitopes in a patient with hepatitis C virus, leading to an altered pattern of immunodominance that may have contributed to the failure of antiviral immunity in this patient.

Association between IL23R and ERAP1 polymorphisms and sacroiliac or spinal MRI inflammation in spondyloarthritis: DESIR cohort data

Background

To investigate the association between 12 single nucleotide polymorphisms (SNPs) located on ERAP1 and IL23R with the presence of inflammation on the sacroiliac joint (SIJ) or spinal magnetic resonance imagery (MRI) in an early onset spondyloarthritis (SpA) cohort.

Methods

All the patients included in the DESIR cohort with an axial SpA and available DNA at baseline were enrolled in this study (n = 645 patients) and underwent a clinical examination, CRP assay, SIJ and spinal MRI scans. Six SNPs located on ERAP1 (rs30187, rs27044, rs27434, rs17482078, rs10050860, rs2287987) and six SNPs located on IL23R (rs1004819, rs10489629, rs1343151, rs2201841, rs10889677, rs11209032) were genotyped. Univariable analyses were performed to test the association between the genotypes and SIJ and spinal MRI inflammation, as well as disease activity based on Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Ankylosing Spondylitis Disease Activity Score-C-Reactive Protein (ASDAS-CRP) and CRP.

Results

One SNP located on ERAP1 (rs27434) and haplotype CCT of ERAP1 were associated with SIJ inflammation detected by MRI, but these associations were below the Bonferroni corrected threshold of significance. However, one SNP (rs1004819) located on IL23R was associated with SIJ MRI inflammation (rs1004819: TT 42.3%, CT 40.5%, CC 26.5%, p = 0.0005). This locus was also significantly associated with Spondyloarthritis Research Consortium of Canada scores while no association with another inflammatory parameter such as BASDAI, ASDAS-CRP, CRP or Berlin MRI score was identified in this population.

Conclusion

One locus of the IL23R gene was associated with SIJ MRI inflammation and might be a marker of more active disease in recent onset SpA.

Trial registration

clinicaltrials.gov, NCTO 164 8907

Association analysis of ERAP1 gene single nucleotide polymorphism in susceptibility to ankylosing spondylitis in Iranian population

Background Ankylosing spondylitis (AS) is a debilitating spondyloarthropathy that has been associated with variation in several genes. Human leukocyte antigen (HLA)-B27 constructs an impaired structure, culminating in recognition and activation of immune system. Impaired function of Endoplasmic reticulum aminopeptidase (ERAP) 1, which primes peptides to be loaded in HLA molecules, has strongly been associated with AS proneness. Here, we intended to investigate the possible association of ERAP1 gene single nucleotide polymorphisms (SNPs) with AS susceptibility in Iranian patients. Methods Two-hundred and twenty AS patients and 220 healthy controls were enrolled in this study. DNA was extracted from blood samples and then was genotyped for rs27044, rs17482078, and rs10050860 polymorphism by SSP-PCR approach. Results It was seen that G allele and GG genotype of rs27044 SNP significantly increased the risk of AS that was even stronger in HLA-B27 positive patients. Moreover, the T allele and TT genotype of rs10050860 polymorphism were associated with increased risk of the disease in both all and HLA-B27 positive AS group. Two haplotypes were associated with the risk of AS and there was linkage disequilibrium between SNPs. Two SNPs were associated with clinicopathological manifestations of AS subjects. Conclusions This association study replicated the role ofERAP1 gene polymorphisms with the risk of AS in an Iranian population.

Associations between ERAP1 polymorphisms and susceptibility to ankylosing spondylitis: A meta-analysis of East Asian Population

BACKGROUND

The genetic factor is importantly enrolled in the pathogenesis of ankylosing spondylitis (AS) and haplotype leukocyte antigen (HLA)-B27 is the most well-known. However, only 1% to 5% of B27-positive individuals will develop AS, and it confers only 20% to 30% of the overall genetic risks, indicating more genes other than HLA-B27 may play important roles in AS pathologies. The present study aims to investigate whether the polymorphisms of endoplasmic reticulum aminopeptidase 1 (ERAP1) is associated with increased risk of AS susceptibility.

METHODS

The Cochrane library, Pubmed, and Embase databases were carefully searched for potential researches published before May 30, 2018. The title, abstract, and full text were assessed to determine whether the paper was suitable for inclusion. The pooled odds ratios (ORs) with 95% confidence intervals (CIs) were presented to assess the associations between ERAP1 polymorphisms and AS susceptibility.

RESULTS

The study finally enrolled 10 papers, 4 matched single nucleotide polymorphisms (SNPs) of ERAP1 (rs27044, rs27434, rs30187, and rs27037), and a total of 30552 patients (12492 with AS and 18060 for control). No significant difference was found between the AS susceptibility and polymorphisms of rs27044 and rs27434. However, there was a significant association between ERAP1 polymorphisms rs30187 and rs27037 (T vs C, OR, 1.322, 95% CI = 1.240-10410, P <.05; T vs G, OR, 1.247, 95% CI = 1.149-1.353; P <.05; respectively) and AS susceptibility.

CONCLUSION

There was a significant association between ERAP1 polymorphisms (rs30187 and rs27037) and increased risk of AS susceptibility.

How ERAP1 and ERAP2 Shape the Peptidomes of Disease-Associated MHC-I Proteins

Four inflammatory diseases are strongly associated with Major Histocompatibility Complex class I (MHC-I) molecules: birdshot chorioretinopathy (HLA-A^*29:02), ankylosing spondylitis (HLA-B^*27), Behçet's disease (HLA-B^*51), and psoriasis (HLA-C^*06:02). The endoplasmic reticulum aminopeptidases (ERAP) 1 and 2 are also risk factors for these diseases. Since both enzymes are involved in the final processing steps of MHC-I ligands it is reasonable to assume that MHC-I-bound peptides play a significant pathogenetic role. This review will mainly focus on recent studies concerning the effects of ERAP1 and ERAP2 polymorphism and expression on shaping the peptidome of disease-associated MHC-I molecules in live cells. These studies will be discussed in the context of the distinct mechanisms and substrate preferences of both enzymes, their different patterns of genetic association with various diseases, the role of polymorphisms determining changes in enzymatic activity or expression levels, and the distinct peptidomes of disease-associated MHC-I allotypes. ERAP1 and ERAP2 polymorphism and expression induce significant changes in multiple MHC-I-bound peptidomes. These changes are MHC allotype-specific and, without excluding a degree of functional inter-dependence between both enzymes, reflect largely separate roles in their processing of MHC-I ligands. The studies reviewed here provide a molecular basis for the distinct patterns of genetic association of ERAP1 and ERAP2 with disease and for the pathogenetic role of peptides. The allotype-dependent alterations induced on distinct peptidomes may explain that the joint association of both enzymes and unrelated MHC-I alleles influence different pathological outcomes.

Translating GWAS in rheumatic disease: approaches to establishing mechanism and function for genetic associations with ankylosing spondylitis

Ankylosing spondylitis (AS) is a highly heritable chronic inflammatory arthritis characterized by osteoproliferation, fusion of affected joints and systemic manifestations. Many disease associations for AS have been reported through genome-wide association studies; however, identifying modulated genes and functional mechanism remains challenging. This review summarizes current genetic associations involving AS and describes strategic approaches for functional follow-up of disease-associated variants. Fine mapping using methods leveraging Bayesian approaches are outlined. Evidence highlighting the importance of context specificity for regulatory variants is reviewed, noting current evidence in AS for the relevant cell and tissue type to conduct such analyses. Technological advances for understanding the regulatory landscape within which functional variants may act are discussed using exemplars. Approaches include defining regulatory elements based on chromatin accessibility, effects of variants on genes at a distance through evidence of physical interactions (chromatin conformation capture), expression quantitative trait loci mapping and single-cell methodologies. Opportunities for mechanistic studies to investigate the function of specific variants, regulatory elements and genes enabled by genome editing using clustered regularly interspaced short palindromic repeats/Cas9 are also described. Further progress in our understanding of the genetics of AS through functional genomic and epigenomic approaches offers new opportunities to understand mechanism and develop innovative treatments.

Allele-specific Alterations in the Peptidome Underlie the Joint Association of HLA-A*29:02 and Endoplasmic Reticulum Aminopeptidase 2 (ERAP2) with Birdshot Chorioretinopathy

Virtually all patients of the rare inflammatory eye disease birdshot chorioretinopathy (BSCR) carry the HLA-A*29:02 allele. BSCR is also associated with endoplasmic reticulum aminopeptidase 2 (ERAP2), an enzyme involved in processing HLA class I ligands, thus implicating the A*29:02 peptidome in this disease. To investigate the relationship between both risk factors we employed label-free quantitative mass spectrometry to characterize the effects of ERAP2 on the A*29:02-bound peptidome. An ERAP2-negative cell line was transduced with lentiviral constructs containing GFP-ERAP2 or GFP alone, and the A*29:02 peptidomes from both transduced cells were compared. A similar analysis was performed with two additional A*29:02-positive, ERAP1-concordant, cell lines expressing or not ERAP2. In both comparisons the presence of ERAP2 affected the following features of the A*29:02 peptidome: 1) Length, with increased amounts of peptides >9-mers, and 2) N-terminal residues, with less ERAP2-susceptible and more hydrophobic ones. The paradoxical effects on peptide length suggest that unproductive binding to ERAP2 might protect some peptides from ERAP1 over-trimming. The influence on N-terminal residues can be explained by a direct effect of ERAP2 on trimming, without ruling out and improved processing in concert with ERAP1. The alterations in the A*29:02 peptidome suggest that the association of ERAP2 with BSCR is through its effects on peptide processing. These differ from those on the ankylosing spondylitis-associated HLA-B*27. Thus, ERAP2 alters the peptidome of distinct HLA molecules as a function of their specific binding preferences, influencing different pathological outcomes in an allele-dependent way.

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.

Genetic Variants in ERAP1 and ERAP2 Associated With Immune-Mediated Diseases Influence Protein Expression and the Isoform Profile

OBJECTIVE

Endoplasmic reticulum aminopeptidase 1 (ERAP-1) and ERAP-2, encoded on chromosome 5q15, trim endogenous peptides for HLA-mediated presentation to the immune system. Polymorphisms in ERAP1 and/or ERAP2 are strongly associated with several immune-mediated diseases with specific HLA backgrounds, implicating altered peptide handling and presentation as prerequisites for autoreactivity against an arthritogenic peptide. Given the thorough characterization of disease risk-associated polymorphisms that alter ERAP activity, this study aimed instead to interrogate the expression effect of chromosome 5q15 polymorphisms to determine their effect on ERAP isoform and protein expression.

METHODS

RNA sequencing and genotyping across chromosome 5q15 were performed to detect genetic variants in ERAP1 and ERAP2 associated with altered total gene and isoform-specific expression. The functional implication of a putative messenger RNA splice-altering variant on ERAP-1 protein levels was validated using mass spectrometry.

RESULTS

Polymorphisms associated with ankylosing spondylitis (AS) significantly influenced the transcript and protein expression of ERAP-1 and ERAP-2. Disease risk-associated polymorphisms in and around both genes were also associated with increased gene expression. Furthermore, key risk-associated ERAP1 variants were associated with altered transcript splicing, leading to allele-dependent alternate expression of 2 distinct isoforms and significant differences in the type of ERAP-1 protein produced.

CONCLUSION

In accordance with studies demonstrating that polymorphisms that increase aminopeptidase activity predispose to immune disease, the increased risk also attributed to increased expression of ERAP1 and ERAP2 supports the notion of using aminopeptidase inhibition to treat AS and other ERAP-associated conditions.

Associations of ERAP1 coding variants and domain specific interaction with HLA-C∗06 in the early onset psoriasis patients of India

Interferon-γ-induced aminopeptidase ERAP1 trims peptides within the endoplasmic reticulum so that they can be loaded onto MHC class I and presented to the CD8⁺ T-cells. ERAP1 association and its interaction with HLA-C∗06 is controversial across different populations. We have investigated the association and possible functional role of non-synonymous SNPs at different exons of ERAP1 (rs26653: Arg127Pro, rs30187: Lys528Arg and rs27044: Gln730Glu) and their interactions with HLA-C∗06 in psoriasis. Significant associations of HLA-C∗06 (OR=5.47, P<2.2×10^-16), rs30187 (OR 1.35, P=7.4×10^-4) and rs27044 (OR=1.24, P=5.8×10^-3) were observed. All three ERAP1 SNPs showed significant association only for HLA-C∗06 positive patients, while rs30187 and rs27044 showed significant association only for early onset patients (rs30187: OR=1.47, P=9.6×10^-5; rs27044: OR=1.36, P=3.3×10^-4). No differential expression of ERAP1 was observed either between paired uninvolved and involved skin tissues of psoriasis patients or between non-risk and risk variants in the involved skin. Significant epistatic interaction was observed between HLA-C∗06 and the SNP (rs27044) located at the peptide-binding cavity of ERAP1. Evolutionary conservation analysis among mammals showed confinement of Lys528 and Gln730 within highly conserved regions of ERAP1 and suggested the possible detrimental effect of this allele in ERAP1 regulation.

The interplay between HLA-B27 and ERAP1/ERAP2 aminopeptidases: from anti-viral protection to spondyloarthritis

The human leukocyte antigen class I gene HLA-B27 is the strongest risk factor for ankylosing spondylitis (AS), a chronic inflammatory arthritic disorder. More recently, the Endoplasmic Reticulum Aminopeptidase (ERAP) 1 and 2 genes have been identified by genome wide association studies (GWAS) as additional susceptibility factors. In the ER, these aminopeptidases trim the peptides to a length suitable to fit into the groove of the major histocompatibility complex (MHC) class I molecules. It is noteworthy that an epistatic interaction between HLA-B27 and ERAP1, but not between HLA-B27 and ERAP2, has been highlighted. However, these observations suggest a paramount centrality for the HLA-B27 peptide repertoire that determines the natural B27 immunological function, i.e. the T cell antigen presentation and, as a by-product, elicits HLA-B27 aberrant behaviours: (i) the misfolding leading to ER stress responses and autophagy and (ii) the surface expression of homodimers acting as ligands for innate immune receptors. In this context, it has been observed that the HLA-B27 carriers, besides being prone to autoimmunity, display a far better surveillance to some viral infections. This review focuses on the ambivalent role of HLA-B27 in autoimmunity and viral protection correlating its functions to the quantitative and qualitative effects of ERAP1 and ERAP2 polymorphisms on their enzymatic activity.

Functional Interaction of the Ankylosing Spondylitis-Associated Endoplasmic Reticulum Aminopeptidase 2 With the HLA-B*27 Peptidome in Human Cells

OBJECTIVE

To determine the influence of endoplasmic reticulum aminopeptidase 2 (ERAP-2) expression on the HLA-B*27 peptidome in live cells.

METHODS

Using immunoaffinity chromatography and acid extraction, HLA-B*27:05-bound peptides were isolated from 2 ERAP-2-negative lymphoblastoid cell lines and 1 ERAP-2-positive lymphoblastoid cell line expressing functionally indistinguishable ERAP-1 variants. More than 2,000-4,000 B*27:05 ligands were identified from each cell line, and their relative abundance was established by quantitative tandem mass spectrometry and MaxQuant-based peptide analyses. Pairwise comparisons were used to determine the structural features of peptides whose relative abundance was dependent on the presence of ERAP-2. Synthetic peptide digestions were performed with recombinant ERAP-1 and ERAP-2. Peptide affinity was estimated with standard algorithms.

RESULTS

The B*27:05 peptidome from ERAP-2-positive cells showed 3-4% fewer peptides with N-terminal basic residues than did the peptidome from ERAP-2-negative cells. Among the shared peptides, those most abundant in the presence of ERAP-2 included more nonamers, fewer decamers, and fewer N-terminal basic residues than the peptides predominant in ERAP-2-negative cells. These ERAP-2-dependent changes did not alter the global affinity of the B*27:05 peptidome.

CONCLUSION

ERAP-2 significantly influences the B*27:05-bound peptidome by destroying some ligands and decreasing the abundance of many more ligands with N-terminal basic residues, while increasing the abundance of nonamers. The former effects are best explained by direct ERAP-2 trimming. The effects on peptide length might be attributed to ERAP-2-induced activation of ERAP-1 trimming. These data support the notion of a peptide-mediated mechanism as the basis for the association of ERAP-2 with ankylosing spondylitis. Analogous effects on other major histocompatibility complex class I peptidomes might explain the involvement of ERAP-2 in HLA-B27-negative spondyloarthritis.

Ankylosing spondylitis: beyond genome-wide association studies

PURPOSE OF REVIEW

This article discusses genomic investigations in ankylosing spondylitis (AS) beyond genome-wide association (GWA) studies, but prior to this, genetic variants achieving genome-wide significance will be summarized highlighting key pathways contributing to disease pathogenesis.

RECENT FINDINGS

Evidence suggests that disease pathogenesis is attributed to a complex interplay of genetic, environmental and immunological factors. GWA studies have greatly enhanced our understanding of AS pathogenesis by illuminating distinct immunomodulatory pathways affecting innate and acquired immunity, most notably the interleukin-23/interleukin-17 pathway. However, despite the wealth of new information gleaned from such studies, a fraction of the heritability (24.4%) has been explained. This review will focus on investigations beyond GWA studies including copy number variants, gene expression profiling, including microRNA (miRNA), epigenetics, rare variants and gene-gene interactions.

SUMMARY

To address the 'missing heritability' and advance beyond GWA studies, a concerted effort involving rethinking of study design and implementation of newer technologies will be required. The coming of age of next-generation sequencing and advancements in epigenetic and miRNA technologies, combined with familial-focused investigations using well-characterized cohorts, is likely to reveal some of the hidden genomic mysteries associated with AS.

Pathogenesis of ankylosing spondylitis - recent advances and future directions

Over the past 5 years, advances in high-throughput techniques and studies involving large cohorts of patients have led to considerable advances in the identification of novel genetic associations and immune pathways involved in ankylosing spondylitis (AS). These discoveries include genes encoding cytokine receptors, transcription factors, signalling molecules and transport proteins. Although progress has been made in understanding the functions and potential pathogenic roles of some of these molecules, much work remains to be done to comprehend their complex interactions and therapeutic potential in AS. In this Review, we outline the current knowledge of AS pathogenesis, including genetic risk associations, HLA-B27-mediated pathology, perturbations in antigen-presentation pathways and the contribution of the type 3 immune response.

The Behçet's disease-associated variant of the aminopeptidase ERAP1 shapes a low-affinity HLA-B*51 peptidome by differential subpeptidome processing

A low-activity variant of endoplasmic reticulum aminopeptidase 1 (ERAP1), Hap10, is associated with the autoinflammatory disorder Behçet's disease (BD) in epistasis with HLA-B*51, which is the main risk factor for this disorder. The role of Hap10 in BD pathogenesis is unknown. We sought to define the effects of Hap10 on the HLA-B*51 peptidome and to distinguish these effects from those due to HLA-B*51 polymorphisms unrelated to disease. The peptidome of the BD-associated HLA-B*51:08 subtype expressed in a Hap10-positive cell line was isolated, characterized by mass spectrometry, and compared with the HLA-B*51:01 peptidome from cells expressing more active ERAP1 allotypes. We additionally performed synthetic peptide digestions with recombinant ERAP1 variants and estimated peptide-binding affinity with standard algorithms. In the BD-associated ERAP1 context of B*51:08, longer peptides were generated; of the two major HLA-B*51 subpeptidomes with Pro-2 and Ala-2, the former one was significantly reduced, and the latter was increased and showed more ERAP1-susceptible N-terminal residues. These effects were readily explained by the low activity of Hap10 and the differential susceptibility of X-Pro and X-Ala bonds to ERAP1 trimming and together resulted in a significantly altered peptidome with lower affinity. The differences due to ERAP1 were clearly distinguished from those due to HLA-B*51 subtype polymorphism, which affected residue frequencies at internal positions of the peptide ligands. The alterations in the nature and affinity of HLA-B*51·peptide complexes probably affect T-cell and natural killer cell recognition, providing a sound basis for the joint association of ERAP1 and HLA-B*51 with BD.

The Human Leukocyte Antigen (HLA)-B27 Peptidome in Vivo, in Spondyloarthritis-susceptible HLA-B27 Transgenic Rats and the Effect of Erap1 Deletion

HLA-B27 is a class I major histocompatibility (MHC-I) allele that confers susceptibility to the rheumatic disease ankylosing spondylitis (AS) by an unknown mechanism. ERAP1 is an aminopeptidase that trims peptides in the endoplasmic reticulum for binding to MHC-I molecules. ERAP1 shows genetic epistasis with HLA-B27 in conferring susceptibility to AS. Male HLA-B27 transgenic rats develop arthritis and serve as an animal model of AS, whereas female B27 transgenic rats remain healthy. We used large scale quantitative mass spectrometry to identify over 15,000 unique HLA-B27 peptide ligands, isolated after immunoaffinity purification of the B27 molecules from the spleens of HLA-B27 transgenic rats. Heterozygous deletion of Erap1, which reduced the Erap1 level to less than half, had no qualitative or quantitative effects on the B27 peptidome. Homozygous deletion of Erap1 affected approximately one-third of the B27 peptidome but left most of the B27 peptidome unchanged, suggesting the possibility that some of the HLA-B27 immunopeptidome is not processed in the presence of Erap1. Deletion of Erap1 was permissive for the AS-like phenotype, increased mean peptide length and increased the frequency of C-terminal hydrophobic residues and of N-terminal Ala, Ser, or Lys. The presence of Erap1 increased the frequency of C-terminal Lys and Arg, of Glu and Asp at intermediate residues, and of N-terminal Gly. Several peptides of potential interest in AS pathogenesis, previously identified in human cell lines, were isolated. However, rats susceptible to arthritis had B27 peptidomes similar to those of non-susceptible rats, and no peptides were found to be uniquely associated with arthritis. Whether specific B27-bound peptides are required for AS pathogenesis remains to be determined. Data are available via ProteomeXchange with identifier PXD005502.

Separate effects of the ankylosing spondylitis associated ERAP1 and ERAP2 aminopeptidases determine the influence of their combined phenotype on the HLA-B*27 peptidome

Ankylosing spondylitis (AS) is an inflammatory disease strongly associated with the Major Histocompatibility Complex class I (MHC-I) allotype HLA-B*27. The endoplasmic reticulum aminopeptidases (ERAP)1 and 2, which trim peptides to their optimal length for MHC-I binding, are also susceptibility factors for this disease. Both highly active ERAP1 variants and ERAP2 expression favor AS, whereas loss-of-function ERAP1 and loss-of-expression ERAP2 variants are protective. Yet, only ERAP1 is in epistasis with HLA-B*27. We addressed two issues concerning the functional interaction of ERAP1 and ERAP2 with the HLA-B*27 peptidome in human cells: 1) distinguishing the effects of ERAP1 from those of ERAP2, and 2) determining the influence of ERAP2 in distinct ERAP1 contexts. Quantitative comparisons of the HLA-B*27:05 peptidomes from cells with various ERAP1/ERAP2 phenotypes were carried out. When cells expressing ERAP2 and either high or low activity ERAP1 variants were compared, increased amounts of nonamers, relative to longer ligands, and decreased amounts of peptides with Ala1, were observed in the more active ERAP1 context. When cells expressing ERAP2 in a low activity ERAP1 context or lacking ERAP2 but expressing a highly active ERAP1 variant were compared, the same effects on peptide length and Ala1, but also significantly lower amounts of peptides with N-terminal basic residues and lower affinity of the peptidome, were observed in the ERAP2-positive context. Thus, ERAP1 and ERAP2 have significant and distinct effects on the HLA-B*27 peptidome, suggesting that both enzymes largely act as separate entities in vivo. This may explain their different patterns of association with AS.

A single endoplasmic reticulum aminopeptidase-1 protein allotype is a strong risk factor for Behçet's disease in HLA-B*51 carriers

INTRODUCTION

Endoplasmic reticulum aminopeptidase-1 (ERAP1) protein is highly polymorphic with numerous missense amino acid variants. We sought to determine the naturally occurring ERAP1 protein allotypes and their contribution to Behçet's disease.

METHODS

Genotypes of all reported missense ERAP1 gene variants with 1000 Genomes Project EUR superpopulation frequency >1% were determined in 1900 Behçet's disease cases and 1779 controls from Turkey. ERAP1 protein allotypes and their contributions to Behçet's disease risk were determined by haplotype identification and disease association analyses.

RESULTS

One ERAP1 protein allotype with five non-ancestral amino acids was recessively associated with disease (p=3.13×10-6, OR 2.55, 95% CI 1.70 to 3.82). The ERAP1 association was absent in individuals who lacked HLA-B*51. Individuals who carry HLA-B*51 and who are also homozygous for the haplotype had an increased disease odds compared with those with neither risk factor (p=4.80×10-20, OR 10.96, 95% CI 5.91 to 20.32).

DISCUSSION

The Behçet's disease-associated ERAP1 protein allotype was previously shown to have poor peptide trimming activity. Combined with its requirement for HLA-B*51, these data suggest that a hypoactive ERAP1 allotype contributes to Behçet's disease risk by altering the peptides available for binding to HLA-B*51.

A Comparative Analysis of the Peptide Repertoires of HLA-DR Molecules Differentially Associated With Rheumatoid Arthritis

OBJECTIVE

To evaluate similarity of the peptide repertoires bound to HLA-DR molecules that are differentially associated with rheumatoid arthritis (RA), and to define structural features of the shared peptides.

METHODS

Peptide pools bound to HLA-DRB1*01:01, HLA-DRB1*04:01, and HLA-DRB1*10:01 (RA associated) and those bound to HLA-DRB1*15:01 (non-RA-associated) were purified and analyzed by liquid chromatography (LC) matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MS) and LC-ion-trap MS. Peptide pools from each allotype were compared in terms of size, protein origin, composition, and affinity (both theoretical and experimental with some peptides). Finally, 1 peptide sequenced from DR1, DR4, and DR10, but not from DR15, was modeled in complex with all 4 HLA-DRB1 molecules and HLA-DRB5*01:01.

RESULTS

A total of 6,309 masses and 962 unique peptide sequences were compared. DR10 shared 29 peptides with DR1, 9 with DR4, and 1 with DR15; DR1 shared 6 peptides with DR4 and 9 with DR15; and DR4 and DR15 shared 4 peptides. The direct identification of peptide ligands indicated that DR1 and DR10 were the most similar molecules regarding the peptides that they could share. The peptides common to these molecules contained a high proportion of Leu at P4 and basic residues at P8 binding core positions.

CONCLUSION

The degree of overlap between peptide repertoires associated with different HLA-DR molecules is low. The repertoires associated with DR1 and DR10 have the highest similarity among the molecules analyzed (∼10% overlap). Among the peptides shared between DR1 and DR10, a high proportion contained Leu(4) and basic residues at the P8 position of the binding core.

New insights toward the pathogenesis of ankylosing spondylitis; genetic variations and epigenetic modifications

Ankylosing spondylitis (AS) is a chronic inflammatory autoimmune disease, characterized by typically an axial arthritis. AS is the prototype of a group of disorders called spondyloarthropathies, which is believed to have common clinical manifestations and genetic predisposition. To date, the exact etiology of AS remains unclear. Over the past few years, however, the role of genetic susceptibility and epigenetic modifications caused through environmental factors have been extensively surveyed with respect to the pathogenesis of AS, resulted in important advances. This review article focuses on the recent advances in the field of AS research, including HLA and non-HLA susceptibility genes identified in genome-wide association studies (GWAS), and aberrant epigenetic modifications of gene loci associated with AS. HLA genes most significantly linked with AS susceptibility include HLA-B27 and its subtypes. Numerous non-HLA genes such as those in ubiquitination, aminopeptidases and MHC class I presentation molecules like ERAP-1 were also reported. Moreover, epigenetic modifications occurred in AS has been summarized. Taken together, the findings presented in this review attempt to explain the circumstance by which both genetic variations and epigenetic modifications are involved in triggering and development of AS. Nonetheless, several unanswered dark sides continue to clog our exhaustive understanding of AS. Future researches in the field of epigenetics should be carried out to extend our vision of AS etiopathogenesis.

ERAAP Shapes the Peptidome Associated with Classical and Nonclassical MHC Class I Molecules

The peptide repertoire presented by classical as well as nonclassical MHC class I (MHC I) molecules is altered in the absence of the endoplasmic reticulum aminopeptidase associated with Ag processing (ERAAP). To characterize the extent of these changes, peptides from cells lacking ERAAP were eluted from the cell surface and analyzed by high-throughput mass spectrometry. We found that most peptides found in wild-type (WT) cells were retained in the absence of ERAAP. In contrast, a subset of "ERAAP-edited" peptides was lost in WT cells, and ERAAP-deficient cells presented a unique "unedited" repertoire. A substantial fraction of MHC-associated peptides from ERAAP-deficient cells contained N-terminal extensions and had a different molecular composition than did those from WT cells. We found that the number and immunogenicity of peptides associated with nonclassical MHC I was increased in the absence of ERAAP. Conversely, only peptides presented by classical MHC I were immunogenic in ERAAP-sufficient cells. Finally, MHC I peptides were also derived from different intracellular sources in ERAAP-deficient cells.

The Peptidome of Behçet's Disease-Associated HLA-B*51:01 Includes Two Subpeptidomes Differentially Shaped by Endoplasmic Reticulum Aminopeptidase 1

OBJECTIVE

To characterize the peptidome of the Behçet's disease-associated HLA-B*51:01 allotype as well as the differential features of major peptide subsets and their distinct endoplasmic reticulum aminopeptidase 1 (ERAP-1)-mediated processing.

METHODS

The endogenous B*51:01-bound peptidome was characterized from 721.221 transfectant cells, after affinity chromatography and acid extraction, by tandem mass spectrometry. Recombinant ERAP-1 variants were used to digest synthetic B*51:01 ligands. HLA and transporter associated with antigen processing (TAP) binding affinities of peptide ligands were calculated with well-established algorithms. ERAP-1 and ERAP-2 from 721.221 cells were characterized by genomic sequencing and Western blotting.

RESULTS

The B*51:01 peptidome consisted of 29.5% octamers, 61.7% nonamers, 4.8% decamers, and 4.0% longer peptides. The major peptide motif consisted of Pro and Ala at position 2, aliphatic/aromatic position 3 residues, and Val and Ile at the C-terminal position. The ligands with Pro or Ala at position 2 constituted 2 distinct subpeptidomes. Peptides with Pro at position 2 showed higher affinity for B*51:01 and lower affinity for TAP than those with Ala at position 2. Most important, both peptide subsets differed drastically in the susceptibility of their position 1 residues to ERAP-1, revealing a distinct influence of this enzyme on both subpeptidomes, which may alter their balance, affecting the global affinity of B*51:01-peptide complexes.

CONCLUSION

ERAP-1 has a significant influence on the B*51:01 peptidome and its affinity. This influence is based on very distinct effects on the 2 subpeptidomes, whereby only peptides in the subpeptidome with Ala at position 2 are extensively destroyed, except when their position 1 residues are ERAP-1 resistant. This pattern provides a mechanism for the epistatic association of ERAP-1 and B*51:01 in Behçet's disease.

Silencing or inhibition of endoplasmic reticulum aminopeptidase 1 (ERAP1) suppresses free heavy chain expression and Th17 responses in ankylosing spondylitis

OBJECTIVE

Human leucocyte antigen (HLA)-B27 and endoplasmic reticulum aminopeptidase 1 (ERAP1) are strongly associated with ankylosing spondylitis (AS). ERAP1 is a key aminopeptidase in HLA class I presentation and can potentially alter surface expression of HLA-B27 free heavy chains (FHCs). We studied the effects of ERAP1 silencing/inhibition/variations on HLA-B27 FHC expression and Th17 responses in AS.

METHODS

Flow cytometry was used to measure surface expression of HLA class I in peripheral blood mononuclear cells (PBMCs) from patients with AS carrying different ERAP1 genotypes (rs2287987, rs30187 and rs27044) and in ERAP1-silenced/inhibited/mutated HLA-B27-expressing antigen presenting cells (APCs). ERAP1-silenced/inhibited APCs were cocultured with KIR3DL2CD3ε-reporter cells or AS CD4+ T cells. Th17 responses of AS CD4+ T cells were measured by interleukin (IL)-17A ELISA and Th17 intracellular cytokine staining. FHC cell surface expression and Th17 responses were also measured in AS PBMCs following ERAP1 inhibition.

RESULTS

The AS-protective ERAP1 variants, K528R and Q730E, were associated with reduced surface FHC expression by monocytes from patients with AS and HLA-B27-expressing APCs. ERAP1 silencing or inhibition in APCs downregulated HLA-B27 FHC surface expression, reduced IL-2 production by KIR3DL2CD3ε-reporter cells and suppressed the Th17 expansion and IL-17A secretion by AS CD4+ T cells. ERAP1 inhibition of AS PBMCs reduced HLA class I FHC surface expression by monocytes and B cells, and suppressed Th17 expansion.

CONCLUSIONS

ERAP1 activity determines surface expression of HLA-B27 FHCs and potentially promotes Th17 responses in AS through binding of HLA-B27 FHCs to KIR3DL2. Our data suggest that ERAP1 inhibition has potential for AS treatment.

Endoplasmic reticulum-associated amino-peptidase 1 and rheumatic disease: genetics

PURPOSE OF REVIEW

This article will review the genetic evidence implicating ERAP1, which encodes the endoplasmic reticulum-associated amino-peptidase 1, in susceptibility to rheumatic disease.

RECENT FINDINGS

Genetic variants and haplotypes of ERAP1 are associated with AS, psoriasis, and Behçet's disease in people of varying ancestries. In each of these diseases, disease-associated variants of ERAP1 have been shown to interact with disease-associated class I human leukocyte antigen alleles to influence disease risk. Functionally, disease-associated missense variants of ERAP1 concertedly alter ERAP1 enzymatic function, both quantitatively and qualitatively, whereas other disease-associated variants influence ERAP1 expression. Therefore, ERAP1 haplotypes (or allotypes) should be examined as functional units. Biologically, this amounts to an examination of the gene regulation and function of the protein encoded by each allotype. Genetically, the relationship between disease risk and ERAP1 allotypes should be examined to determine whether allotypes or individual variants produce the most parsimonious risk models.

SUMMARY

Future investigations of ERAP1 should focus on comprehensively characterizing naturally occurring ERAP1 allotypes, examining the enzymatic function and gene expression of each allotype, and identifying specific allotypes that influence disease susceptibility.

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

INTRODUCTION

Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS.

METHODS

Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR.

RESULTS

There was no significant difference in HLA-class I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls.

DISCUSSION

Large differences were not seen in HLA-B27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

Genetics of ankylosing spondylitis--insights into pathogenesis

Ankylosing spondylitis (AS), an immune-mediated arthritis, is the prototypic member of a group of conditions known as spondyloarthropathies that also includes reactive arthritis, psoriatic arthritis and enteropathic arthritis. Patients with these conditions share a clinical predisposition for spinal and pelvic joint dysfunction, as well as genetic associations, notably with HLA-B(*)27. Spondyloarthropathies are characterized by histopathological inflammation in entheses (regions of high mechanical stress where tendons and ligaments insert into bone) and in the subchondral bone marrow, and by abnormal osteoproliferation at involved sites. The association of AS with HLA-B(*)27, first described >40 years ago, led to hope that the cause of the disease would be rapidly established. However, even though many theories have been advanced to explain how HLA-B(*)27 is involved in AS, no consensus about the answers to this question has been reached, and no successful treatments have yet been developed that target HLA-B27 or its functional pathways. Over the past decade, rapid progress has been made in discovering further genetic associations with AS that have shed new light on the aetiopathogenesis of the disease. Some of these discoveries have driven translational ideas, such as the repurposing of therapeutics targeting the cytokines IL-12 and IL-23 and other factors downstream of this pathway. AS provides an excellent example of how hypothesis-free research can lead to major advances in understanding pathogenesis and to the development of innovative therapeutic strategies.

Structural Basis for Antigenic Peptide Recognition and Processing by Endoplasmic Reticulum (ER) Aminopeptidase 2

Endoplasmic reticulum (ER) aminopeptidases process antigenic peptide precursors to generate epitopes for presentation by MHC class I molecules and help shape the antigenic peptide repertoire and cytotoxic T-cell responses. To perform this function, ER aminopeptidases have to recognize and process a vast variety of peptide sequences. To understand how these enzymes recognize substrates, we determined crystal structures of ER aminopeptidase 2 (ERAP2) in complex with a substrate analogue and a peptidic product to 2.5 and 2.7 Å, respectively, and compared them to the apo-form structure determined to 3.0 Å. The peptides were found within the internal cavity of the enzyme with no direct access to the outside solvent. The substrate analogue extends away from the catalytic center toward the distal end of the internal cavity, making interactions with several shallow pockets along the path. A similar configuration was evident for the peptidic product, although decreasing electron density toward its C terminus indicated progressive disorder. Enzymatic analysis confirmed that visualized interactions can either positively or negatively impact in vitro trimming rates. Opportunistic side-chain interactions and lack of deep specificity pockets support a limited-selectivity model for antigenic peptide processing by ERAP2. In contrast to proposed models for the homologous ERAP1, no specific recognition of the peptide C terminus by ERAP2 was evident, consistent with functional differences in length selection and self-activation between these two enzymes. Our results suggest that ERAP2 selects substrates by sequestering them in its internal cavity and allowing opportunistic interactions to determine trimming rates, thus combining substrate permissiveness with sequence bias.

Association study of ankylosing spondylitis and polymorphisms in ERAP1 gene in Zhejiang Han Chinese population

The susceptibility loci of ERAP1 polymorphisms have been found to be strongly associated with ankylosing spondylitis (AS). The researches in multiple ethnic cohorts suggested that the population attributable risk in ERAP1 polymorphisms is at a high significance level. This study was undertaken to estimate the prevalence and incidence of subsets of AS and investigate the specific variants of ERAP1 polymorphisms in AS susceptibility, in the Han ethnic Chinese population in Zhejiang Province. AS patients were selected, diagnosed, and confirmed by a qualified rheumatologist. The basal clinical and demographic characteristics were compared with all subjects. Genotypes for eight selected single nucleotide polymorphisms (SNPs) in ERAP1 gene (rs27038, rs27037, rs27434, rs27980, rs7711564, rs30187, rs10050860, and rs17482078) were determined by using the Sequenom MassARRAY iPLEX platform in Zhejiang Han Chinese population. Association analyses were performed on the whole genotyped data set in 707 unrelated ankylosing spondylitis cases and 837 ethnically matched controls. We observed the strongest association between AS and HLA-B27, which confers over 90 % of ankylosing spondylitis cases. Moreover, we found three loci of ERAP1 polymorphisms were at a high significance level (rs27037 P = 0.00451; rs27434 P = 0.00012; rs27980 P = 0.00682) with AS in Zhejiang population. We also confirmed polymorphism locus of ERAP1 previously reported association with AS (rs27434; P = 5.3 × 10(-12)). Our results indicated a difference in the mechanism of susceptibility loci in subsets of Zhejiang Han Chinese population and provided further evidence that rs27434 is the key polymorphism associated with AS in ERAP1 gene.