Correlation of ERAP Gene Polymorphism with Antibody Response Induced by Sequential Immunization of Polio Vaccine
-
摘要:
目的 分析内质网氨肽酶(endoplasmic reticulum aminopeptidases,ERAP)基因多态性与脊髓灰质炎疫苗序贯免疫诱导的中和抗体应答的相关性。 方法 选取243名来自广西壮族自治区并完成2剂灭活脊髓灰质炎疫苗和1剂二价口服脊髓灰质炎减毒疫苗接种的壮族受试者,检测免前和基础免疫28 d血清中I、II、III型脊髓灰质炎中和抗体水平,采用TaqMan探针基因分型法对单核苷酸多态性(single nucleotide polymorphism,SNP)进行分型。共选择8个SNP位点,6个ERAP1基因(rs27037、rs27044、rs30187、rs26618、rs26653、rs3734016)和2个ERAP2基因(rs2549782、rs2248374),计算各SNPs的等位基因和基因型频率,分析各SNPs与各型抗体应答的相关性。 结果 在I型脊髓灰质炎抗体应答中,携带rs2549782-G和rs2248374-A等位基因个体抗体几何平均滴度(geometric mean titer,GMT)低于携带rs2549782-T和rs2248374-G等位基因个体[均为(11.590±1.979) vs (11.950±1.895),P = 0.031];rs2549782基因型GT和rs2248374基因型AG诱导的中和抗体低于rs2549782基因型TT和rs2248374基因型GG[均为(11.741±0.141) vs (12.378±0.157),P = 0.045]。 结论 ERAP2基因多态性可能影响脊髓灰质炎疫苗诱导的抗体水平。 Abstract:Objective To analyze the relationship between endoplasmic reticulum aminopeptidases (ERAP) gene polymorphisms and serum polio antibodies induced by sequential polio vaccine immunization. Methods A total of 243 Zhuang individuals were selected from the Guangxi Zhuang Autonomous Region who received two doses of inactivated polio vaccine and one dose of bivalent oral polio vaccine. Polio-neutralizing antibodies types I, II and III were tested from pre-immunization and 28 days after immunizations, and six ERAP1 and two ERAP2 SNPs were genotyped using TaqMan probe genotyping. The allele frequency and genotype frequency were calculated for each SNP, and the association between the SNPs and the polio antibody response was analyzed. Results It was found that individuals carrying the rs2549782-G andrs2248374-A allele of the ERAP1 gene had lower levels of type I polio-neutralizing antibodies compared to those carrying the T and G alleles, respectively (both of 11.590±1.979 vs 11.950±1.895, Padj = 0.031). In addition, it was observed that GT and AG genotypes of rs2549782 of rs2248374 exhibited lower GMT type I polio-neutralizing antibodies than those of TT and GG genotypes (both of 11.741±0.141 vs 12.378±0.157, Padj = 0.045). Conclusion Polio vaccine-induced antibody responses may be associated with ERAP2 gene polymorphism. -
女性宫颈癌(cervical cancer,CC)是较为常见的妇科恶性肿瘤之一,在世界范围内其发病率为第二,仅低于发病率第一的乳腺癌,在发展中国家和地区则排第一位[1]。宫颈癌的发生及发展是在多种因素作用下进行的,其中最重要因素为长时间存在感染人乳头瘤病毒(human papilloma virys,HPV)[2]。近年来研究发现,从正常组织到宫颈癌的变化是一个复杂过程,是多因素分阶段参与发生的,仅仅HPV感染不足以导致宫颈病变甚至宫颈癌的发生,多种其他基因功能和表达水平的改变也影响了宫颈癌的病理过程[3-4]。微小RNA (microRNA)是一种内源性非编码小 RNA分子,主要通过与其靶 mRNA3’ UTR结构完全或不完全配对,通过细胞内复杂的调控系统,对细胞代谢、生长、分化和凋亡等生物学过程发挥调节作用,近年来研究发现,microRNA多方面参与肿瘤细胞的生物化学过程,并且主要在肿瘤中发挥癌基因或抑癌基因的功能[5]。其中,microRNA-21 (miR-21)是较早发现、存在较广的人类 miRNA之一,在多种肿瘤中呈现高表达,如肺恶性肿瘤、肝恶性肿瘤、肝细胞胆管癌、乳腺癌、脑胶质瘤恶性肿瘤等,作为研究 miRNA在恶性肿瘤的发生及进展中作用的其中指标之一[6]。miR-21通过参与调控多种抑癌基因的表达,如PDCD4,BCL2和CCL20等,发挥生物学效应[7-8]。相关研究表明,miR-21通过反向调节多种肿瘤抑制相关基因的表达而在宫颈癌的发生中起到关键作用,在宫颈癌HeLa细胞株中降低miR-21表达可以使细胞的生长分化受到抑制[9-10]。有研究发现端粒酶逆转录酶(hTERT)在宫颈癌中表达明显上调,使用 hTERT的抑制剂 AZT和 BIBR1532或者通过 RNAi干扰技术敲除 hTERT的表达可以明显抑制宫颈癌细胞系 HeLa细胞的增殖和迁移,已有研究证实,在脑胶质瘤细胞中 miR-21表达水平的下调能够引起 hTERT表达水平的下降[11]。目前国内外暂未发现关于miR-21与hTERT的相关性及miR-21在宫颈癌中对hTERT所起的调控作用的相关报道。
本实验运用RT-PCR及West Blot法检测HeLa细胞里miR-21水平及hTERT蛋白的相对表达情况;利用脂质体包裹的技术方法进行细胞转染,使 miR-21模拟物及抑制剂存在HeLa细胞株中,并在上下调节细胞内miR-21表达水平的状态下,观察miR-21与 hTERT之间的联系,探讨miR-21是否也能靶向调节 hTERT而起调控作用,为宫颈癌的发病机制和新型治疗手段的研究提供创新性的研究思路。
1. 材料与方法
1.1 材料
长沙赢润生物技术有限公司买入人宫颈癌 HeLa细胞株,在 Invitrogen机构买入转染脂质体 liposome2000 (lipo2000),在 Hyclone机构买入 RPMI1640培养基,在Gibco机构买入 Opti-MEM培养基;由广州复能公司合成miR-21的模拟物、模拟对照、抑制剂和抑制剂对照;由MRC公司、Beyotime公司分别提供实验所需Trizol及BCA蛋白浓度测定试剂盒。
1.2 方法
1.2.1 细胞培养
用 RPMI1640培养液培养在37 ℃、5% CO2的恒温细胞培养箱中进行细胞培养,培养液中含1000 U/mL青霉素、10%胎牛血清、0.1 g/mL链霉素。利用胰酶消化细胞,从而进行传代培养,最佳时机为80%~90% 细胞密度时。
1.2.2 细胞转染
在转染24 h之前,在6孔培养板中接种一定数目的HeLa细胞,使细胞密度稳定维持在30%~50%之间,以便达到最佳转染条件。在lipo2000 使用说明书指导下进行转染细胞。用250 μL Opti-MEM稀释100 pmol miRNA模拟实验组、模拟对照实验组、抑制剂实验组及抑制剂实验对照组,混匀后在室温条件下孵育5 min;用250 μL Opti-MEM稀释5 μL lipo2000,轻轻混匀并室温孵育5 min;将两者混合,室温孵育20 min。将最终混合液加入准备的培养孔中,轻轻混匀;将培养板放置在CO2培养箱中进行培养(培养温度为37 ℃),24 h 后进行实验检测。
1.2.3 RT-PCR法测定miR-21、hTERT mRNA
在转染24 h后,在 Trizol试剂盒中进行数据提取,检测实验细胞RNA,测量出细胞的OD值,并将其进行电泳和逆转录。同样用 RT-PCR检测法检测脂质体组、模拟物组、模拟对照组、抑制剂组和抑制剂对照组中hTERT mRNA的表达。
1.2.4 测量各组实验细胞内 hTERT蛋白的相关表达
用 West Blot检测法进行细胞内蛋白含量的检测,采用BCA蛋白定量试剂盒将细胞中的蛋白含量进行检测,并进行电泳,膜转移,抗体温育和一抗的显影,孵育60 min后,用TBST洗至无脱脂奶粉状态。同时,将通用二抗在常温下温育60 min,孵育完成后,TBST 洗3 次,每次持续时间为5 min,即可进行显影。3 min后,准备保鲜薄膜,将完成的PVDF膜晾干后放置其上,滴加混匀的化学发光底物(取底物A和B分别500 µL混匀为1 mL液体,避光保存)。2 min后,把制作完成的PVDF膜覆盖好轻放于避光箱里,并固定。在观察到荧光效果后可停止压片,待片子在显影液中显出条带后,在蒸馏水中洗片,再次重复显影及冲洗过程一次,晾干片子,并将片子进行标记。测量得出细胞内 hTERT蛋白的相关表达。
1.3 统计学处理
利用SPSS 24.0 软件来完成对实验数据的分析,hTERT mRNA及蛋白在各组中的表达情况用(
$\bar{x}\pm s$ )表示,单因素多组比较采用单因素方差分析,两组间的比较采用两独立样本t检验,P < 0.05 为差异有统计学意义。2. 结果
2.1 miR-21及hTERTmRNA 在HeLa中的表达情况
miR-21及hTERTmRNA 在HeLa 细胞株中表达呈阳性,见图1。
2.2 检测hTERT mRNA 在各组别细胞中表达
转染后,对hTERT mRNA 在模拟物组、模拟物对照组、反义miR-21组别(抑制剂组)、阴性对照组(抑制剂对照组)及脂质体组中的 mRNA含量相对表达测定的结果进行各组间比较,发现模拟物组中的hTERT mRNA相对含量最低,反义miR-21组即抑制剂组hTERT mRNA 相对表达水平最高,明显高于hTERT mRNA 在阴性对照组及脂质体组的表达水平,差异有统计学意义(P < 0.05),见表1、图2。
表 1 hTERT mRNA 在各组别细胞中的表达情况($\bar{x}\pm s$ )Table 1. Expression of hTERT mRNA in different groups of cells ($\bar{x}\pm s$ )组别 hTERT mRNA t P 模拟物组 0.51 ± 0.33 4.478 0.011 模拟物对照组 0.69 ± 0.19 抑制剂组 1.33 ± 0.39 5.387 0.005 抑制剂对照组 0.83 ± 0.26 脂质体组 0.58 ± 0.16 注:单因素多组间比较(F = 5.403,P = 0.002)。 2.3 检测hTERT蛋白在各组别细胞中的含量
hTERT蛋白在模拟物组、模拟物对照组、反义miR-21组(抑制剂组)、阴性对照组(抑制剂对照组)及脂质体组中的蛋白相对表达量分别为(0.41±0.18)、(0.62±0.18)、(1.45±0.60)、(0.83±0.15) 和(0.82±0.30);将数据进行各组间比较,与其他四组相比,模拟物组中hTERT 蛋白的相对含量最低,表达显著下降,反义miR-21组即抑制剂组hTERT蛋白相对含量最高,明显高于hTERT蛋白在阴性对照组及脂质体组的表达水平,差异有统计学意义(P < 0.05),见表2、图3。
表 2 hTERT 蛋白在各组别细胞中的表达情况($\bar{x}\pm s$ )(1)Table 2. hTERT protein expression in each group of cells ($\bar{x}\pm s$ )(1)组别 实验结果 t P 模拟物组 0.41 ± 0.18 5.072 0.0071 模拟物对照组 0.62 ± 0.18 抑制剂组 1.45 ± 0.60 5.85 0.0043 抑制剂对照组 0.83 ± 0.15 脂质体组 0.82 ± 0.30 表 2 hTERT 蛋白在各组别细胞中的表达情况($\bar{x}\pm s$ )(2)Table 2. hTERT protein expression in each group of cells ($\bar{x}\pm s$ )(2)组别 hTERT蛋白 t P 模拟物组 0.41 ± 0.18 5.072 0.007 模拟物对照组 0.62 ± 0.18 抑制剂组 1.45 ± 0.60 5.85 0.004 抑制剂对照组 0.83 ± 0.15 脂质体组 0.82 ± 0.30 3. 讨论
宫颈癌是一种生殖系统的恶性肿瘤,因其高死亡率,极大地威胁着广大女性同胞的健康安全[12]。研究表明宫颈癌的发生与HPV的长时间感染直接相关,但并非唯一相关,而存在多方面因素多种条件的影响。近年来研究发现,miRNA非常规表达对宫颈病变产生影响,可致使宫颈癌的产生及影响宫颈癌的预后,例如 miR-21、miR-34a 等[13-14]。
miRNAs是一种小分子RNA,这种RNA是由内源性非蛋白质编码,它是由一段类似发卡环结构且长度为70~80 nt的单链RNA前体构成[15]。国内外相关研究发现,miR-21在肿瘤细胞中发生作用,参与细胞的生长、侵袭等生理生物学过程,从而导致肿瘤的发生及进展[16-17]。miR-21相关的靶基因包括程序性细胞死亡因子4(PDCD4)、PTEN、TPM1、STAT3 等[18-19]。苏俊玲等[20]的实验中,发现pre-miR-21对宫颈癌HeLa 细胞的增殖和生长有明显促进作用,而anti-miR-21则对宫颈癌HeLa 细胞的增长存在明显反向抑制作用。Xu L等[21]发现miR-21 通过抑制PTEN 的作用,致使宫颈癌HeLa 细胞的增长及生殖得到加快。miRNA与所调控基因并非单一的,一种miRNA可以对数个基因的表达进行调节,而同一基因调控也是由不同miRNA所完成。在本研究中,宫颈癌细胞株HeLa中检测出miR-21呈高表达状态,而Yao等[22]的实验也得出相同的结论。
端粒酶(Telomerase)是在细胞中能使端粒延长的一种端粒延长酶,是常见的一种的核蛋白逆转录酶。端粒酶主要有三个构成部分,即端粒酶RNA组分(hTERC)、假尿嘧啶合成酶和端粒酶逆转录酶(hTERT)[23]。在端粒酶中,hTERT是一个催化亚活性结构单位,是端粒酶的核心构成组分之一,hTERT 的表达水平直接影响端粒酶的活性,两者是正向平行相关的。限定hTERT 转录水平,端粒酶的活性将同时受到限定,因此认为hTERT在转录水平的活化致使端粒酶激活对肿瘤细胞内的生物化学效应产生影响,是肿瘤发生及进展的重要环节[24]。李春林等[25]的研究发现hTERT在宫颈癌细胞和80%以上宫颈癌组织中均有较高表达,在正常宫颈组织和细胞中为阴性,差异显著。本研究利用脂质体 liposome2000包裹宫颈癌细胞株 HeLa并予以转染,发现在模拟物组中hTERT mRNA 及hTERT 蛋白的相对含量最低,表达显著下降,反义miR-21即抑制剂组hTERT mRNA 及hTERT蛋白相对含量最高,明显高于阴性对照组、脂质体组的相对表达水平(P < 0.05),得出miR-21模拟物及miR-21抑制剂,分别促进和抑制宫颈癌细胞株HeLa内miR-21的表达水平,同时将细胞内miR-21表达水平进行上下调节,发现miR-21对hTERT目标基因有反向控制的作用,即miR- 21可通过下调hTERT mRNA和hTERT蛋白表达水平调控靶基因在翻译过程中的作用,从而产生负性调控的作用。Bachand 等[26] 研究证明,端粒酶 hTR 内有hTERT两个独立的结合位点,hTR 与 hTERT相互作用,二者均是端粒酶活性调节必不可少的因素。当我们以端粒酶模板区 RNA 反义核酸和端粒酶逆转录酶区反义核酸同时封闭 hTR和 hTERT,推测端粒酶模板区活性及逆转录酶区活性同时得以抑制,既能激活 RnaseH,降解 mRNA,又能减弱逆转录调控,表现为反义 hTR和反义 hTERT 对宫颈癌 HeLa细胞端粒酶活性下降的协同效应。
综上所述,可以发现miR-21与宫颈癌HeLa的生长及凋亡紧密关联。miR- 21可以调节hTERT mRNA 和hTERT蛋白表达水平,对hTERT靶基因有负向调节的作用,miR-21 调控hTERT 低表达使宫颈癌HeLa细胞生长受到抑制,影响肿瘤细胞的生长及肿瘤血管的生成,这一发现可指导宫颈癌基因治疗的研究。但是miR-21是否调控其他靶基因参与宫颈癌的发生发展及miR-21调控hTERT的具体机制需要更深入的实验与研究。
-
表 1 SNPs的位点信息
Table 1. SNPs information
基因 SNP位点 位置 突变类型 ERAP1 rs27037(G/T) chr5:96758990 intron rs27044(C/G) chr5:96783148 Gln730Glu(Q>E) rs30187(C/T) chr5:96788627 Lys528Arg(K>R) rs26618(T/C) chr5:96795133 Ile276Met(I>M) rs26653(C/G) chr5:96803547 Arg127Pro(R>P) rs3734016(C/T) chr5:96803761 Glu56Lys(E>K) ERAP2 rs2549782(T/G) chr5:96895296 Lys392Asn(K>N) rs2248374(G/A) chr5:96900192 splice region variant 表 2 研究对象基本特征(
$\bar x \pm s $ )Table 2. Basic characteristics of the subjects (
$\bar x \pm s $ )型别 性别(n) GMT t P I型 男(128) 11.628 ± 2.043 1.350 0.178 女(115) 11.964 ± 1.814 II型 男(128) 6.110 ± 1.806 0.945 0.345 女(115) 6.335 ± 1.910 III型 男(128) 11.409 ± 1.699 1.292 0.198 女(115) 11.675 ± 1.492 表 3 ERAP基因的SNPs位点等位基因与脊灰中和抗体(GMT)水平的相关性(
$\bar x \pm s $ )Table 3. Correlation between ERAP SNPs alleles and GMT of neutralizing antibody against poliovirus (
$\bar x \pm s $ )基因 SNP位点 脊灰抗体型 等位基因 n GMT P ERAP1 rs27037 I G 295 11.860 ± 1.873 0.351 T 191 11.670 ± 2.038 II G 295 6.205 ± 1.765 0.961 T 191 6.233 ± 1.986 III G 295 11.600 ± 1.463 0.838 T 191 11.430 ± 1.801 rs27044 I C 252 11.870 ± 1.913 0.282 G 234 11.700 ± 1.968 II C 252 6.272 ± 1.733 0.416 G 234 6.156 ± 1.937 III C 252 11.620 ± 1.450 0.760 G 234 11.450 ± 1.755 rs30187 I C 249 11.860 ± 1.915 0.338 T 237 11.710 ±1.966 II C 249 6.297 ± 1.784 0.245 T 237 6.131 ± 1.923 III C 249 11.610 ± 1.467 0.815 T 237 11.460 ± 1.738 rs26618 I C 157 11.740 ± 1.902 0.579 T 329 11.810 ± 1.960 II C 157 6.277 ± 1.804 0.497 T 329 6.187 ± 1.878 III C 157 11.650 ± 1.410 0.657 T 329 11.480 ± 1.690 rs26653 I C 247 11.760 ± 1.947 0.657 G 239 11.820 ± 1.936 II C 247 6.114 ± 1.899 0.155 G 239 6.322 ± 1.802 III C 247 11.470 ± 1.725 0.910 G 239 11.800 ± 1.472 rs3734016 I C 428 11.730 ± 1.954 0.077 T 58 12.190 ± 1.792 II C 428 6.170 ± 1.858 0.129 T 58 6.558 ± 1.790 III C 428 11.510 ± 1.634 0.582 T 58 11.720 ± 1.374 ERAP2 rs2549782 I G 221 11.590 ± 1.979 0.031* T 265 11.950 ± 1.895 II G 221 6.187 ± 1.836 0.992 T 265 6.241 ± 1.870 III G 221 11.560 ± 1.578 0.943 T 265 11.510 ± 1.630 rs2248374 I A 221 11.590 ± 1.979 0.031* G 265 11.950 ± 1.895 II A 221 6.187 ± 1.836 0.992 G 265 6.241 ± 1.870 III A 221 11.560 ± 1.578 0.943 G 265 11.510 ± 1.630 *P < 0.05。 表 4 ERAP基因SNPs的基因型与脊灰I型中和抗体GMT水平的相关性(
$\bar x \pm s $ )Table 4. Correlation between ERAP SNPs genotypes and GMT of neutralizing antibody against polio type1 (
$\bar x \pm s $ )ERAP SNP位点 抗体型别 基因型 n GMT P ERAP1 rs27037 I型 GG 86 12.099 ± 1.792 0.109 GT 123 11.535 ± 1.948 TT 34 11.910 ± 2.202 rs27044 I型 CC 58 12.010 ± 2.001 0.563 CG 136 11.750 ± 1.841 GG 49 11.630 ± 2.150 rs30187 I型 CC 56 11.990 ± 2.033 0.626 CT 137 11.750 ± 1.822 TT 50 11.650 ± 2.167 rs26618 I型 CC 20 11.930 ± 1.591 0.672 CT 117 11.670 ± 2.005 TT 106 11.890 ± 1.939 rs26653 I型 CC 57 11.780 ± 2.083 0.830 CG 133 11.730 ± 1.839 GG 53 11.930 ± 2.064 rs3734016 I型 CC 189 11.690 ± 1.968 0.235 CT 50 12.080 ± 1.854 TT 4 12.900 ± 1.286 ERAP2 rs2549782 I型 GG 48 11.600 ± 2.046 0.043* GT 125 11.590 ± 1.942 TT 70 12.280 ± 1.802 rs2248374 I型 AA 48 11.600 ± 2.046 0.043* AG 125 11.590 ± 1.942 GG 70 12.280 ± 1.802 *P < 0.05。 -
[1] Tan P L,Jacobson R M,Poland G A,et al. Twin studies of immunogenicity-determining the genetic contribution to vaccine failure[J]. Vaccine,2001,19(17-19):2434-2439. doi: 10.1016/S0264-410X(00)00468-0 [2] Newport M J,Goetghebuer T,Weiss H A,et al. Genetic regulation of immune responses to vaccines in early life[J]. Genes Immun,2004,5(2):122-129. doi: 10.1038/sj.gene.6364051 [3] Neefjes J,Jongsma M L,Paul P,et al. Towards a systems understanding of MHC class I and MHC class II antigen presentation[J]. Nat Rev Immunol,2011,11(12):823-836. doi: 10.1038/nri3084 [4] Vyas J M,Van Der Veen A G,Ploegh H L. The known unknowns of antigen processing and presentation[J]. Nat Rev Immunol,2008,8(8):607-618. doi: 10.1038/nri2368 [5] Yao Y,Liu N,Zhou Z,et al. Influence of ERAP1 and ERAP2 gene polymorphisms on disease susceptibility in different populations[J]. Hum Immunol,2019,80(5):325-334. doi: 10.1016/j.humimm.2019.02.011 [6] Hattori A,Tsujimoto M. Endoplasmic reticulum aminopeptidases:Biochemistry,physiology and pathology[J]. J Biochem,2013,154(3):219-228. doi: 10.1093/jb/mvt066 [7] Nguyen T T,Chang S C,Evnouchidou I,et al. Structural basis for antigenic peptide precursor processing by the endoplasmic reticulum aminopeptidase ERAP1[J]. Nat Struct Mol Biol,2011,18(5):604-613. doi: 10.1038/nsmb.2021 [8] Evnouchidou I,Birtley J,Seregin S,et al. A common single nucleotide polymorphism in endoplasmic reticulum aminopeptidase 2 induces a specificity switch that leads to altered antigen processing[J]. J Immunol,2012,189(5):2383-2392. doi: 10.4049/jimmunol.1200918 [9] 傅宇婷,黄腾,英志芳,等. 脊髓灰质炎免疫策略调整进程中不同免疫程序接种效果的比较[J]. 中国生物制品学杂志,2019,32(7):786-793. doi: 10.13200/j.cnki.cjb.002680 [10] 赵志梅,傅宇婷,施红媛,等. Sabin株脊髓灰质炎灭活疫苗与不同剂型Ⅰ型Ⅲ型脊髓灰质炎减毒活疫苗序贯接种的抗体滴度比较[J]. 中华疾病控制杂志,2019,23(4):402-406. doi: 10.16462/j.cnki.zhjbkz.2019.04.007 [11] Liao G,Li R,Li C,et al. Phase 3 Trial of a sabin strain-based inactivated poliovirus vaccine[J]. J Infect Dis,2016,214(11):1728-1734. doi: 10.1093/infdis/jiw433 [12] Liao G,Li R,Li C,et al. Safety and immunogenicity of inactivated poliovirus vaccine made from sabin strains:A phase II,randomized,positive-controlled trial[J]. J Infect Dis,2012,205(2):237-243. doi: 10.1093/infdis/jir723 [13] Chen H,Li L,Weimershaus M,et al. ERAP1-ERAP2 dimers trim MHC I-bound precursor peptides; implications for understanding peptide editing[J]. Sci Rep,2016,6(1):28902. doi: 10.1038/srep28902 [14] Stratikos E,Stamogiannos A,Zervoudi E,et al. A role for naturally occurring alleles of endoplasmic reticulum aminopeptidases in tumor immunity and cancer pre-disposition[J]. Front Oncol,2014,4(4):363. [15] Alvarez-Navarro C,López DE Castro J A. ERAP1 structure,function and pathogenetic role in ankylosing spondylitis and other MHC-associated diseases[J]. Mol Immunol,2014,57(1):12-21. [16] Stamogiannos A,Koumantou D,Papakyriakou A,et al. Effects of polymorphic variation on the mechanism of endoplasmic reticulum aminopeptidase 1[J]. Mol Immunol,2015,67(2 Pt B):426-435. [17] Coulombe-Huntington J,Lam K C,Dias C,et al. Fine-scale variation and genetic determinants of alternative splicing across individuals[J]. PLoS Genet,2009,5(12):e1000766. doi: 10.1371/journal.pgen.1000766 [18] Cagliani R,Riva S,Biasin M,et al. Genetic diversity at endoplasmic reticulum aminopeptidases is maintained by balancing selection and is associated with natural resistance to HIV-1 infection[J]. Hum Mol Genet,2010,19(23):4705-4714. doi: 10.1093/hmg/ddq401 [19] Biasin M,Sironi M,Saulle I,et al. Endoplasmic reticulum aminopeptidase 2 haplotypes play a role in modulating susceptibility to HIV infection[J]. AIDS,2013,27(11):1697-1706. doi: 10.1097/QAD.0b013e3283601cee [20] Liu S,Cao D,Shen Y,et al. The ERAP gene is associated with HCV chronic infection in a Chinese Han population[J]. Hum Immunol,2017,78(11-12):731-738. doi: 10.1016/j.humimm.2017.10.005 [21] Kuiper J J W,Setten J V,Devall M,et al. Functionally distinct ERAP1 and ERAP2 are a hallmark of HLA-A29-(Birdshot) uveitis[J]. Hum Mol Genet,2018,27(24):4333-4343. [22] Vanhille D L,Hill L D,Hilliard D D,et al. A novel ERAP2 haplotype structure in a chilean population: Implications for ERAP2 protein expression and preeclampsia risk[J]. Mol Genet Genomic Med,2013,1(2):98-107. doi: 10.1002/mgg3.13 [23] Castro-Santos P,Moro-García M A,Marcos-Fernández R,et al. ERAP1 and HLA-C interaction in inflammatory bowel disease in the Spanish population[J]. Innate Immun,2017,23(5):476-481. doi: 10.1177/1753425917716527 [24] Jiang R,Liu X,Sun X,et al. Immunogenicity and safety of the inactivated poliomyelitis vaccine made from sabin strains in a phase IV clinical trial for the vaccination of a large population[J]. Vaccine,2021,39(9):1463-1471. doi: 10.1016/j.vaccine.2021.01.027 [25] He H,Wang Y,Deng X,et al. Immunogenicity of three sequential schedules with sabin inactivated poliovirus vaccine and bivalent oral poliovirus vaccine in Zhejiang,China: an open-label,randomised,controlled trial[J]. Lancet Infect Dis,2020,20(9):1071-1079. doi: 10.1016/S1473-3099(19)30738-8 [26] Jacobson R M,Poland G A,Vierkant R A,et al. The association of class I HLA alleles and antibody levels after a single dose of measles vaccine[J]. Hum Immunol,2003,64(1):103-109. doi: 10.1016/S0198-8859(02)00741-3 [27] Ovsyannikova I G,Jacobson R M,Vierkant R A,et al. Replication of rubella vaccine population genetic studies:Validation of HLA genotype and humoral response associations[J]. Vaccine,2009,27(49):6926-6931. doi: 10.1016/j.vaccine.2009.08.109 [28] Mineta M,Tanimura M,Tana T,et al. Contribution of HLA class I and class II alleles to the regulation of antibody production to hepatitis B surface antigen in humans[J]. Int Immunol,1996,8(4):525-531. doi: 10.1093/intimm/8.4.525 [29] Hayney M S,Poland G A,Dimanlig P,et al. Polymorphisms of the TAP2 gene may influence antibody response to live measles vaccine virus[J]. Vaccine,1997,15(1):3-6. doi: 10.1016/S0264-410X(96)00133-8 [30] Yucesoy B,Talzhanov Y,Johnson V J,et al. Genetic variants within the MHC region are associated with immune responsiveness to childhood vaccinations[J]. Vaccine,2013,31(46):5381-5391. doi: 10.1016/j.vaccine.2013.09.026 [31] Lankat-Buttgereit B,Tampe R. The transporter associated with antigen processing:Function and implications in human diseases[J]. Physiol Rev,2002,82(1):187-204. doi: 10.1152/physrev.00025.2001 -