# Associations betweenERAP1 Gene Polymorphisms and Psoriasis Susceptibility: A Meta-Analysis of Case-Control Studies

Authors: Xiujuan Wu; Zongfeng Zhao
Journal: BioMed Research International
 (2021)

---

## Abstract

This study is to investigate the relationship of endoplasmic reticulum aminopeptidase 1 (ERAP1) gene polymorphisms with psoriasis. Five databases of PubMed, China National Knowledge Infrastructure (CNKI), Embase, Web of Science, and Cochrane Library were searched for potential studies until 25 December 2019. Newcastle-Ottawa Scale (NOS) was used to evaluate the quality of included studies. Meta-analysis was performed with PRISMA. A total of 9 case-control studies including 4858 psoriasis cases and 10,542 healthy controls were included. Three loci of ERAP1 gene polymorphisms (rs26653, rs30187, and rs27524) were evaluated in this meta-analysis. There was no significant association between rs26653 polymorphism and risk of psoriasis (C vs. G, OR=1.01, 95% CI: 0.80-1.28, P=0.93). However, there was a significant association between rs30187 polymorphisms and psoriasis susceptibility (T vs. C, OR=1.23, 95% CI: 1.15-1.32, P<0.00001) and a significant association between rs27524 polymorphisms and psoriasis susceptibility (A vs. G, OR=1.17, 95% CI: 1.09-1.25, P<0.00001). For there were insufficient data of rs27044, rs151823, rs2248374, and rs2910686, we only conducted a systematic review for them. The rs30187 (C/T) and rs27524 (G/A) polymorphisms of ERAP1 are associated with increased risk of psoriasis. However, no significant association is observed between rs26653 (G/C) polymorphism and risk of psoriasis.


---

## Body

## 1. Introduction

Psoriasis is a common chronic inflammatory disorder [1–3]. Psoriasis affects about 2%–3% of the population worldwide, and it is more prevalent in Western countries [4]. There are five types of psoriasis identified up to now, including plaque psoriasis, eruptive psoriasis, inverse psoriasis, pustular psoriasis, and erythrodermic psoriasis [1]. As a kind of skin disease, psoriasis brings physical and psychological burdens to patients. Therefore, it is urgent to improve the quality of life for patients.Many environmental factors have been reported to be associated with psoriasis, such as direct skin trauma, streptococcal throat infection, and human immunodeficiency virus infection [3]. In recent decades, powerful genome-wide association studies and many other candidate gene approaches have identified numerous genes associated with psoriasis risk [5]. Endoplasmic reticulum aminopeptidase 1 (ERAP1), 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 [6]. ERAP1 belongs to the oxytocinase subgroup of M1 zinc metallopeptidases, shares 49% sequence similarity, and can form heterodimers [7]. The human ERAP1 genes are encoded in the short arm chromosome 5q15 in a 167Kb region in the opposite direction, and probably, they have two shared regulatory elements [8]. Polymorphisms in the ERAP1 gene may affect susceptibility to psoriasis [9]. The rs26653, rs30187, and rs27524 polymorphisms are considered as the three most important loci of ERAP1 [10, 11]. One genome-wide association study confirmed that the polymorphism of ERAP1 rs26653 and rs27044 was associated with psoriasis in Chinese patients [12]. rs27524 polymorphism of ERAP1 is closely related to the susceptibility of psoriasis in Caucasians [13]. rs26653 and rs30187 are associated with early onset psoriasis in Caucasians [14]. Therefore, the polymorphism of ERAP1 rs27524, rs30187, rs26653, and rs27044 may play significant roles in patients with psoriasis. ERAP1 is upexpressed by TNF-α and IFN-γ stimulation, indicating its major function in antigen-presenting machinery [15]. ERAP1 plays a significant role in the migration and proliferation of endothelial cells as a vital factor for vessel regeneration [16]. The interrelation network of ERAP1 and its nearest associated functional protein partners are illustrated in Figure 1.Figure 1
Network ofERAP1 and its closest functional partners. These data were from the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://string-db.org/).Herein, considering the limited sample sizes and poor statistical power of each individual study, we conducted this meta-analysis. Our findings may provide a more comprehensive evaluation on the association ofERAP1 polymorphisms with psoriasis susceptibility.


## 2. Materials and Methods

### 2.1. Search Strategy

The online databases of PubMed, China National Knowledge Infrastructure (CNKI), Embase, Web of Science, and Cochrane Library were searched for eligible studies until 25 December 2019. Medical Subject Headings (MeSH) together with free terms were as keywords for literature searching. The following searching strategies for PubMed were employed: (polymorphism OR single nucleotide polymorphism OR variation OR gene) AND (psoriasis OR “Psoriasis” [17]) AND (“ERAP1 protein, human” [Supplementary Concept] OR ERAP1 OR Adipocyte-derived leucine aminopeptidase, human OR Type 1 tumor necrosis factor receptor shedding aminopeptidase regulator, human OR endoplasmic-reticulum aminopeptidase-1, human).


### 2.2. Inclusion and Exclusion Criteria

Studies in accordance with the following criteria were included: case-control studies with subjects of human; studies with enough original data to calculate odds ratios (ORs); studies with psoriasis diagnosed according to clinical diagnosis criteria; studies on the relationship betweenERAP1 polymorphisms and psoriasis predisposition. Accordingly, studies that met the following criteria were excluded: reviews, conference abstracts, case reports, ascertained alleles, and overlapping data.


### 2.3. Quality Assessment

The methodological quality of each eligible study was evaluated by two investigators, respectively, in compliance with Newcastle-Ottawa Scale (NOS). NOS contains three aspects: selection (four items), comparability (two items), and the outcomes of case-control studies (three items). Each included study was scored as 0-9 points according to these items. A higher score indicated a better quality. Studies withscores≥5 points were considered to have a high quality for further analysis. When disagreements occurred between the two investigators, a third reviewer was consulted to make the final decision.


### 2.4. Data Extraction

The following information was extracted from each study by two independent investigators: first author, year of publication, ethnicity of study population, numbers of cases and controls, and the allele frequencies of theERAP1 polymorphisms, Hardy-Weinberg equilibrium (HWE) results.


### 2.5. Statistical Analysis

This meta-analysis was performed with the Preferred Reporting Items for Systematic Review and Meta-Analyses (PRISMA) [18]. The allele model of inheritance was used to analyze the data due to a lack of sufficient information. The association between ERAP1 gene polymorphisms and psoriasis was calculated by merging ORs with 95% CIs (confidence intervals). The Q-statistical test and I2 test were employed to evaluate the heterogeneity among all included studies [19]. The random-effects model was used to combine the data in the cases of heterogeneity (P<0.1 or I2>50%), or the fixed-effects model was used when it was out of heterogeneity (P>0.1 or I2<50%) [20, 21]. A leave-one-out sensitivity analysis was used to explore the stability of pooled results by omitting each investigated study at a time. A funnel plot of publication bias and forest plots of meta-analysis were plotted with Review Manager Version 5.3 software (Cochrane Collaboration, Software Update, Oxford, United Kingdom).


## 2.1. Search Strategy

The online databases of PubMed, China National Knowledge Infrastructure (CNKI), Embase, Web of Science, and Cochrane Library were searched for eligible studies until 25 December 2019. Medical Subject Headings (MeSH) together with free terms were as keywords for literature searching. The following searching strategies for PubMed were employed: (polymorphism OR single nucleotide polymorphism OR variation OR gene) AND (psoriasis OR “Psoriasis” [17]) AND (“ERAP1 protein, human” [Supplementary Concept] OR ERAP1 OR Adipocyte-derived leucine aminopeptidase, human OR Type 1 tumor necrosis factor receptor shedding aminopeptidase regulator, human OR endoplasmic-reticulum aminopeptidase-1, human).


## 2.2. Inclusion and Exclusion Criteria

Studies in accordance with the following criteria were included: case-control studies with subjects of human; studies with enough original data to calculate odds ratios (ORs); studies with psoriasis diagnosed according to clinical diagnosis criteria; studies on the relationship betweenERAP1 polymorphisms and psoriasis predisposition. Accordingly, studies that met the following criteria were excluded: reviews, conference abstracts, case reports, ascertained alleles, and overlapping data.


## 2.3. Quality Assessment

The methodological quality of each eligible study was evaluated by two investigators, respectively, in compliance with Newcastle-Ottawa Scale (NOS). NOS contains three aspects: selection (four items), comparability (two items), and the outcomes of case-control studies (three items). Each included study was scored as 0-9 points according to these items. A higher score indicated a better quality. Studies withscores≥5 points were considered to have a high quality for further analysis. When disagreements occurred between the two investigators, a third reviewer was consulted to make the final decision.


## 2.4. Data Extraction

The following information was extracted from each study by two independent investigators: first author, year of publication, ethnicity of study population, numbers of cases and controls, and the allele frequencies of theERAP1 polymorphisms, Hardy-Weinberg equilibrium (HWE) results.


## 2.5. Statistical Analysis

This meta-analysis was performed with the Preferred Reporting Items for Systematic Review and Meta-Analyses (PRISMA) [18]. The allele model of inheritance was used to analyze the data due to a lack of sufficient information. The association between ERAP1 gene polymorphisms and psoriasis was calculated by merging ORs with 95% CIs (confidence intervals). The Q-statistical test and I2 test were employed to evaluate the heterogeneity among all included studies [19]. The random-effects model was used to combine the data in the cases of heterogeneity (P<0.1 or I2>50%), or the fixed-effects model was used when it was out of heterogeneity (P>0.1 or I2<50%) [20, 21]. A leave-one-out sensitivity analysis was used to explore the stability of pooled results by omitting each investigated study at a time. A funnel plot of publication bias and forest plots of meta-analysis were plotted with Review Manager Version 5.3 software (Cochrane Collaboration, Software Update, Oxford, United Kingdom).


## 3. Results

### 3.1. Literature Search

The detailed process of the literature screening is shown in Figure2. The primary searching in five databases yielded 186 records, consisting of 49 from PubMed, 41 from Embase, 0 from Cochrane Library, 65 from Web of Science, and 31 from CNKI. After removing 135 duplicate articles, the remaining 51 articles were subjected to screening by title and abstracts; after removing 7 articles irrelevant to ERAP1 and 19 reviews and letters, there were 25 articles that were eligible for full-text review; after excluding 5 article studies with insufficient data and 11 articles with irrelevant data, 9 original case-control studies consisting of 4858 psoriasis cases and 10542 healthy controls were finally included.Figure 2
Flow diagram of literature search and screening.


### 3.2. Study Characteristics

General characterizations of all included studies are demonstrated in Table1. All the included studies were published in English from 2011 to 2018. Among them, two of them were from China [22, 23], another two were from Poland [24, 25], the rest were from India [26], UK [27], Netherlands [28], Sweden [14], and Romania [29]. The quality of included studies was evaluated with NOS. The results showed that each included study had a score of more than 5 points (Table 2), indicating good quality of included studies.Table 1
Main characteristics of included studies.
Study
Year
Country
Ethnicity
Case/control
Case
Control
HWE
Mutant allele
Wild allele
Mutant allele
Wild allele
rs26653 (G/C)
Anamika
2017
India
Indian
911/1120
893
929
1030
1210
0.60
Andrzej
2018
Poland
Polish
461/454
252
670
256
652
0.02
Fu
2018
China
Han
143/149
141
145
122
176
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1155
2571
909
2587
NA
Marta
2018
Poland
Polish
127/101
44
210
70
132
NA
rs30187 (C/T)
Anamika
2017
India
Indian
911/1120
911
911
941
1299
0.59
Andrzej
2018
Poland
Polish
461/454
311
611
302
606
0.67
Anne
2011
UK
Caucasian
76/5200
58
94
3432
6968
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1416
2310
1189
2307
NA
Olivia
2016
Romania
Caucasian
88/139
69
107
89
189
>0.05
rs27524 (G/A)
Andrzej
2018
Poland
Polish
461/454
307
615
311
597
0.83
Bergboer
2012
Netherlands
Caucasian
215/450
172
258
295
605
>0.05
Fu
2018
China
Han
143/149
129
157
113
185
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1453
2273
1259
2237
NA
Yang
2013
China
Han
974/1181
885
1063
945
1417
0.58
Note: HWE: Hardy-Weinberg equilibrium.Table 2
Quality assessment of included studies according to the Newcastle-Ottawa Scale.
Item/study
Adequate definition of cases
Representativeness of cases
Selection of control subjects
Definition of control subjects
Control for important factor or additional factor
Exposure assessment
Same method of ascertainment for all subjects
Nonresponse rate
Total score
Anamika 2017
1
0
0
1
1
1
1
1
6
Andrzej 2018
1
0
1
1
1
1
1
1
7
Anne 2011
1
0
0
1
1
1
1
1
6
Bergboer 2012
1
0
0
1
1
1
1
1
6
Fu 2018
1
0
0
1
1
1
1
1
6
Josefin 2013
1
0
0
1
1
1
1
1
6
Marta 2018
1
0
0
1
1
1
1
1
6
Olivia 2016
1
0
0
1
1
1
1
1
6
Yang 2013
1
0
1
1
1
1
1
1
7


### 3.3. Meta-Analysis Results

Due to the limited number of studies on other loci, we only meta-analyzed three loci of rs26653 (G/C), rs30187 (C/T), and rs27524 (G/A) ofERAP1 polymorphisms.


### 3.4. Relationship between rs26653 (G/C) Polymorphism and Psoriasis Susceptibility

Meta-analysis of the relationship between rs26653 polymorphism and psoriasis susceptibility was conducted based on 5 studies consisting of 3505 cases and 3572 controls [14, 22, 24–26]. There was an enormous heterogeneity among the five included studies (P<0.00001, I2=87%). Therefore, the random-effects model was used. The results indicated that no significant association was observed between rs26653 polymorphism and risk of psoriasis (C vs. G, OR=1.01, 95% CI: 0.80-1.28, P=0.93, Figure 3).Figure 3
Forest plot of rs26653 in the ERAP1 gene and risk of psoriasis.


### 3.5. Relationship between rs30187 (C/T) Polymorphism and Psoriasis Susceptibility

Meta-analysis for the relationship between rs30187 (C/T) polymorphism and psoriasis susceptibility was on the basis of five studies [19–21, 23, 24], including 3399 cases and 8661 healthy controls. A marginal heterogeneity was observed among the five included studies (P=0.11, I2=47%), and therefore, the fixed-effects model was used. The results suggested that there was a significant association between rs30187 polymorphisms and psoriasis susceptibility (T vs. C, OR=1.23, 95% CI: 1.15-1.32, P<0.00001, Figure 4).Figure 4
Forest plot of rs30187 in the ERAP1 gene and risk of psoriasis.


### 3.6. Relationship between rs27524 (G/A) Polymorphism and Psoriasis Susceptibility

Analysis of the association of rs27524 polymorphism and psoriasis susceptibility was also based on five studies [14, 22, 23, 25, 28], including a total of 3656 psoriasis cases and 3982 healthy controls. There was a marginal heterogeneity among the five included studies, and thus, the fixed-effects model was used. The results displayed that there was a significant association between rs27524 polymorphisms and psoriasis susceptibility (A vs. G, OR=1.17, 95% CI: 1.09-1.25, P<0.00001, Figure 5).Figure 5
Forest plot of rs27524 in theERAP1 gene and risk of psoriasis.


### 3.7. Relationship between OtherERAP Gene Polymorphisms and Psoriasis Susceptibility

Since there were only a few studies that reported detailed information on rs27044, rs151823, rs2248374, and rs2910686, we only performed a narrative description for these loci. Three studies reported rs27044 polymorphisms in theERAP1 gene and psoriasis susceptibility. In Wiśniewski et al.’s study [25], rs27044 was not associated with psoriasis in the Polish population (G vs. C, OR=1.05, 95% CI: 0.86-1.29), which was consistent with the results by Popa et al. [29] (G vs. C, OR=1.40, 95% CI:0.92, 2.12). However, rs27044 polymorphism was found to be significantly associated with psoriasis risk in Das et al.’s study [26] (G vs. C, OR=1.24, 95% CI: 1.07–1.43).Two studies [30, 31] reported that rs151823 gene polymorphisms of ERAP1 may be also involved in the pathogenesis of psoriasis and treatment response of biologic therapies to psoriasis.The rs2248374 and rs2910686 are two loci located inERAP2, which plays the same roles as ERAP1 in biological process. In Popa et al.’s study [29], rs2248374 polymorphism was found to be significantly associated with psoriasis risk (A vs. G, OR=1.45, 95% CI: 1.01-2.10) while rs2910686 polymorphism was not associated with psoriasis risk (C vs. T, OR=1.41, 95% CI: 0.97-2.05).


### 3.8. Sensitivity Analysis and Publication Bias

A leave-one-out sensitivity analysis was used to explore the stability of pooled results. The overall effect size did not change after omitting one included study at a time. Publication bias detection was made by visual inspection of the funnel plots. The results did not identify substantial asymmetry (Figure6), indicating no publication bias.Figure 6
Analysis of publication bias. Funnel plot of included studies on rs26653 (a), rs30187 (b), and rs27524 (c) of theERAP1 gene.
(a)
(b)
(c)


## 3.1. Literature Search

The detailed process of the literature screening is shown in Figure2. The primary searching in five databases yielded 186 records, consisting of 49 from PubMed, 41 from Embase, 0 from Cochrane Library, 65 from Web of Science, and 31 from CNKI. After removing 135 duplicate articles, the remaining 51 articles were subjected to screening by title and abstracts; after removing 7 articles irrelevant to ERAP1 and 19 reviews and letters, there were 25 articles that were eligible for full-text review; after excluding 5 article studies with insufficient data and 11 articles with irrelevant data, 9 original case-control studies consisting of 4858 psoriasis cases and 10542 healthy controls were finally included.Figure 2
Flow diagram of literature search and screening.


## 3.2. Study Characteristics

General characterizations of all included studies are demonstrated in Table1. All the included studies were published in English from 2011 to 2018. Among them, two of them were from China [22, 23], another two were from Poland [24, 25], the rest were from India [26], UK [27], Netherlands [28], Sweden [14], and Romania [29]. The quality of included studies was evaluated with NOS. The results showed that each included study had a score of more than 5 points (Table 2), indicating good quality of included studies.Table 1
Main characteristics of included studies.
Study
Year
Country
Ethnicity
Case/control
Case
Control
HWE
Mutant allele
Wild allele
Mutant allele
Wild allele
rs26653 (G/C)
Anamika
2017
India
Indian
911/1120
893
929
1030
1210
0.60
Andrzej
2018
Poland
Polish
461/454
252
670
256
652
0.02
Fu
2018
China
Han
143/149
141
145
122
176
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1155
2571
909
2587
NA
Marta
2018
Poland
Polish
127/101
44
210
70
132
NA
rs30187 (C/T)
Anamika
2017
India
Indian
911/1120
911
911
941
1299
0.59
Andrzej
2018
Poland
Polish
461/454
311
611
302
606
0.67
Anne
2011
UK
Caucasian
76/5200
58
94
3432
6968
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1416
2310
1189
2307
NA
Olivia
2016
Romania
Caucasian
88/139
69
107
89
189
>0.05
rs27524 (G/A)
Andrzej
2018
Poland
Polish
461/454
307
615
311
597
0.83
Bergboer
2012
Netherlands
Caucasian
215/450
172
258
295
605
>0.05
Fu
2018
China
Han
143/149
129
157
113
185
>0.05
Josefin
2013
Sweden
Caucasian
1863/1748
1453
2273
1259
2237
NA
Yang
2013
China
Han
974/1181
885
1063
945
1417
0.58
Note: HWE: Hardy-Weinberg equilibrium.Table 2
Quality assessment of included studies according to the Newcastle-Ottawa Scale.
Item/study
Adequate definition of cases
Representativeness of cases
Selection of control subjects
Definition of control subjects
Control for important factor or additional factor
Exposure assessment
Same method of ascertainment for all subjects
Nonresponse rate
Total score
Anamika 2017
1
0
0
1
1
1
1
1
6
Andrzej 2018
1
0
1
1
1
1
1
1
7
Anne 2011
1
0
0
1
1
1
1
1
6
Bergboer 2012
1
0
0
1
1
1
1
1
6
Fu 2018
1
0
0
1
1
1
1
1
6
Josefin 2013
1
0
0
1
1
1
1
1
6
Marta 2018
1
0
0
1
1
1
1
1
6
Olivia 2016
1
0
0
1
1
1
1
1
6
Yang 2013
1
0
1
1
1
1
1
1
7


## 3.3. Meta-Analysis Results

Due to the limited number of studies on other loci, we only meta-analyzed three loci of rs26653 (G/C), rs30187 (C/T), and rs27524 (G/A) ofERAP1 polymorphisms.


## 3.4. Relationship between rs26653 (G/C) Polymorphism and Psoriasis Susceptibility

Meta-analysis of the relationship between rs26653 polymorphism and psoriasis susceptibility was conducted based on 5 studies consisting of 3505 cases and 3572 controls [14, 22, 24–26]. There was an enormous heterogeneity among the five included studies (P<0.00001, I2=87%). Therefore, the random-effects model was used. The results indicated that no significant association was observed between rs26653 polymorphism and risk of psoriasis (C vs. G, OR=1.01, 95% CI: 0.80-1.28, P=0.93, Figure 3).Figure 3
Forest plot of rs26653 in the ERAP1 gene and risk of psoriasis.


## 3.5. Relationship between rs30187 (C/T) Polymorphism and Psoriasis Susceptibility

Meta-analysis for the relationship between rs30187 (C/T) polymorphism and psoriasis susceptibility was on the basis of five studies [19–21, 23, 24], including 3399 cases and 8661 healthy controls. A marginal heterogeneity was observed among the five included studies (P=0.11, I2=47%), and therefore, the fixed-effects model was used. The results suggested that there was a significant association between rs30187 polymorphisms and psoriasis susceptibility (T vs. C, OR=1.23, 95% CI: 1.15-1.32, P<0.00001, Figure 4).Figure 4
Forest plot of rs30187 in the ERAP1 gene and risk of psoriasis.


## 3.6. Relationship between rs27524 (G/A) Polymorphism and Psoriasis Susceptibility

Analysis of the association of rs27524 polymorphism and psoriasis susceptibility was also based on five studies [14, 22, 23, 25, 28], including a total of 3656 psoriasis cases and 3982 healthy controls. There was a marginal heterogeneity among the five included studies, and thus, the fixed-effects model was used. The results displayed that there was a significant association between rs27524 polymorphisms and psoriasis susceptibility (A vs. G, OR=1.17, 95% CI: 1.09-1.25, P<0.00001, Figure 5).Figure 5
Forest plot of rs27524 in theERAP1 gene and risk of psoriasis.


## 3.7. Relationship between OtherERAP Gene Polymorphisms and Psoriasis Susceptibility

Since there were only a few studies that reported detailed information on rs27044, rs151823, rs2248374, and rs2910686, we only performed a narrative description for these loci. Three studies reported rs27044 polymorphisms in theERAP1 gene and psoriasis susceptibility. In Wiśniewski et al.’s study [25], rs27044 was not associated with psoriasis in the Polish population (G vs. C, OR=1.05, 95% CI: 0.86-1.29), which was consistent with the results by Popa et al. [29] (G vs. C, OR=1.40, 95% CI:0.92, 2.12). However, rs27044 polymorphism was found to be significantly associated with psoriasis risk in Das et al.’s study [26] (G vs. C, OR=1.24, 95% CI: 1.07–1.43).Two studies [30, 31] reported that rs151823 gene polymorphisms of ERAP1 may be also involved in the pathogenesis of psoriasis and treatment response of biologic therapies to psoriasis.The rs2248374 and rs2910686 are two loci located inERAP2, which plays the same roles as ERAP1 in biological process. In Popa et al.’s study [29], rs2248374 polymorphism was found to be significantly associated with psoriasis risk (A vs. G, OR=1.45, 95% CI: 1.01-2.10) while rs2910686 polymorphism was not associated with psoriasis risk (C vs. T, OR=1.41, 95% CI: 0.97-2.05).


## 3.8. Sensitivity Analysis and Publication Bias

A leave-one-out sensitivity analysis was used to explore the stability of pooled results. The overall effect size did not change after omitting one included study at a time. Publication bias detection was made by visual inspection of the funnel plots. The results did not identify substantial asymmetry (Figure6), indicating no publication bias.Figure 6
Analysis of publication bias. Funnel plot of included studies on rs26653 (a), rs30187 (b), and rs27524 (c) of theERAP1 gene.
(a)
(b)
(c)


## 4. Discussion

ERAP1 was identified by GWAS to be preferentially associated with early onset psoriasis among Chinese [12]. ERAP1 polymorphism was first reported to be associated with the disease of ankylosing spondylitis in 2007 [32]. After that, the polymorphisms in ERAP1 have also been found to be associated with Behcet’s disease by genome-wide association studies [13, 33, 34], and since then, replication studies were performed by researchers all over the world. Many polymorphisms in the ERAP1 gene have been shown to be associated with genetic susceptibility to psoriasis [14]. ERAP1 had genetic heterogeneity in a Swedish cohort with early onset ankylosing spondylitis [35]. ERAP1 SNPs are identified as intronic variants in most studies. ERAP1 haplotype is predisposing to psoriasis in East Asians [8]. ERAP2 was established upon controlling for the contribution of ERAP1 in psoriasis [36]. The effect of ERAP2 on disease susceptibility is often masked by ERAP1 because of ERAP1 and ERAP2 within the same haplotype [37].HighERAP1 expression is a risk factor for psoriasis [38]. To the best of our knowledge, the present study is the first and most comprehensive meta-analysis on the association of ERAP1 gene polymorphisms with psoriasis risk. Of the 9 included studies, only Asian and Caucasian population was eligible. Since there was not enough information on the exact number for all genotypes of every individual involved in the original study, we only analyzed the data using an allele model. The meta-analysis results indicated that carriers with the minor allele of rs30187 and rs27524 had an increased risk of psoriasis, while carriers with the minor allele of rs26653 seemed to have neither increased nor decreased risk of psoriasis. Subgroup analyses were not conducted due to a lack of enough original studies.The role of ERAP1 variations is underlined by the association of rs27524 in HLA-Cw6 positive individuals [13, 39]. Studies [40, 41] have shown that rs26653 is associated with a risk of psoriasis independent of hla-c06 expression. rs26653 SNP may be involved in psoriasis by targeting important components such as IL-17/IL-23 [14]. rs26653 (R127 P) in the junction domains could indirectly affect either enzymatic activity or specificity by switching between open and close conformations [14]. rs27044 SNP is the most common SNP in almost all populations [33]. It confers a protective effect due to a significant reduction in aminopeptidase activity [42]. Five studies reporting rs26653 polymorphism showed no significant association with psoriasis susceptibility [14, 22, 24–26]. In addition, five studies identified significant association of rs30187 (C/T) polymorphism with psoriasis [14, 25–27, 29]. The five studies reporting rs27524 polymorphism demonstrated a significant susceptibility to psoriasis [14, 22, 23, 25, 28]. Other studies have shown that rs27044 [25, 26, 29] and rs151823 [30, 31] are associated with susceptibility to psoriasis. The rs2248374 and rs2910686 are two loci located in ERAP2, which plays the same roles as ERAP1 [29].ERAP1 together with ERAP2 belongs to the oxytocinase subfamily of M1 zinc metallopeptidases [43] and contains 20 exons and 19 introns [44]. ERAP1 participates in the final processing of major histocompatibility complex (MHC) class I ligands in the endoplasmatic reticulum by trimming the peptides to a required length for presentation [45, 46]. Apart from the antigen processing pathway of MHC-I molecules, ERAP1 is also involved in the indirect presentation of cytoplasmic MHC-II ligands [47, 48]. The downregulated expression level of ERAP1 will upregulate the antigen-presenting properties and deregulate the immunological effect produced by MHC-I molecules in the host defense against pathogens [49–51]. ERAP1 also participates in the process of oncogenesis [52]. The effect of ERAP1 on angiotensin II modulates angiogenesis and therefore leads to the occurrence of endometrial carcinoma [52]. ERAP1 regulates the immune defense against tumorigenesis by affecting the balance of activating and inhibitory signals through MHC-I proteins [53]. ERAP1 also has other biological functions, such as regulating blood pressure and angiogenesis [41, 43].Several limitations of the present meta-analysis should not be overlooked. First, the study number and the sample size of included studies were relatively small, which might provide a limited conclusion. Second, the heterogeneity between studies was relatively considerable, and thus, the conclusion of this meta-analysis should be quoted with caution. Third, only studies on Caucasian and Asian populations were included. Whether the conclusions are consistent with those on other populations needs further investigation. Fourth, due to the limited number of included studies, meta-regression analysis was not performed. Last, we only searched English and Chinese databases to perform this work. Potentially relevant articles published in other languages may be ignored.


## 5. Conclusions

In summary,ERAP1 polymorphism of rs30187 (C/T) and rs27524 (G/A) was associated with an increased risk of psoriasis. However, no significant association between rs26653 (G/C) polymorphism and psoriasis risk was observed. The results of this study may provide a more comprehensive evaluation of the relationship between ERAP1 polymorphism and psoriasis susceptibility.



---

*Source: 5515868-2021-08-03.xml*