The immune system is activated by MHC-peptide complexes, after which it eliminates infected and cancerous cells in various ways. The APM plays crucial roles in the initiation and development of various human diseases. As components of the APM, ERAP1 and ERAP2 are important determinants of the repertoire of peptides ultimately presented by HLA class I molecules (34-37). Our previous study showed that four polymorphisms in ERAP1 are associated with non-small cell lung cancer in a Chinese population (27). In the current study, we selected six SNPs in ERAP1 and five SNPs in ERAP2, and investigated their distribution in CIN patients, cervical cancer patients, and healthy individuals. Our results showed that rs26618, rs26653, and rs27044 in ERAP1 and rs2287988 in ERAP2 may be associated with cervical cancer risk.
The SNP, rs26618, in ERAP1 leads to an amino acid substitution (I276M) and the current study showed that the CC genotype of this SNP may be associated with an increased risk of cervical cancer (OR = 1.53; 95% CI: 1.14-2.05) compared with TT-CT genotypes. In 2016, Guasp et al. reported that I276M (rs26618) may affect the peptidome of ERAP1 by destroying peptides with p2 Ala, unless the p1 amino acid was resistant to ERAP1 trimming (38), which indicated that rs26618 may be associated with cervical cancer. However, in a Netherlands population, Mehta et al. reported no association between rs26618 and cervical carcinoma. This inconsistency may be due to the different sample sizes used. The sample size used by Mehta et al. was 251, while 2890 individuals were enrolled in the current study.
In 2007, Mehta et al. reported that the C allele of rs26653 in ERAP1 was associated with a higher cervical cancer risk in a Netherlands population (28). In the current study, the G allele (OR = 0.829; 95% CI: 0.738-0.930), compared to the C allele, and the 2GG+CG genotype, compared to the CC genotype (OR = 0.82; 95% CI: 0.73-0.93) of rs26653, were associated with lower cervical cancer risk. In 2014, Stratikos et al. and Alvarez-Navarro et al. reported that rs26653, which is a non-synonymous variant resulting in a P127R substitution, may be associated with ERAP expression (39, 40), and this substitution may also affect the enzymatic activity of ERAP1 in the editing of tumour antigen peptides. This finding may explain the association between rs26653 and cervical cancer risk; however, the mechanisms need to be determined in functional studies.
In 2007, Mehta et al. also found that rs27044 in ERAP1 was associated with cervical cancer risk. In the current study, rs27044 was found to be associated with cervical cancer risk (P = 0.003). The G allele of rs27044 (Q730) was found to be a risk factor for the progression from CIN to cervical cancer (OR = 1.193, 95% CI: 1.062-1.340), which was consistent with the results of Mehta’s study (28). The SNP, rs27044, another non-synonymous variant, leads to a Q730E substitution in the IV catalysis domain of ERAP1 (21) and may change the substrate length preferences of ERAP1 (40). Thus, rs27044 may play a role in cervical cancer by affecting ERAP1 function.
In the current study, we found an association between rs2287988 in ERAP2, which is responsible for a synonymous variant (Q563Q), and cervical cancer. The A allele may be associated with a lower risk of cervical cancer (P = 0.004; OR = 0.843, 95% CI: 0.751-0.946). Moreover, the GG-GA genotype was associated with an increased risk of cervical cancer (P = 0.002; OR = 1.33, 95% CI: 1.11-1.60). However, association studies of this SNP are rare. Previous studies have found that ERAP2 haplotypes containing rs2287988 affect ERAP2 splicing and expression (41, 42). Thus, additional association studies in different populations are necessary to investigate the role of this variant during the initiation and development of cervical cancer.
ERAPs are markedly polymorphic and ERAP haplotypes whose protein products differ at multiple amino acids may affect peptide editing by ERAPs (17, 18, 43, 44). In the current study, we also analysed haplotypes of ERAP SNPs in LD. The results showed that the ERAP1 haplotype, rs26618T-rs26653G-rs27044C-rs30187T-rs3734016C and the ERAP2 haplotypes, rs2248374G-rs2549782T-rs2287988A-rs2548538T-rs1056893T and rs2248374A-rs2549782G-rs2287988G-rs2548538A-rs1056893T may be associated with cervical cancer risk. These results indicated that SNPs in polymorphic genes may have combinatorial effects on disease susceptibility.