疾病监测  2015, Vol. 30 Issue (7): 539-545

扩展功能

文章信息

仇艳, 蒋毅, 刘海灿, 李马超, 万康林
QIU Yan, JIANG Yi, LIU Hai-can, LI Ma-chao, WAN Kang-lin
结核分枝杆菌抗原PE 35和PPE 68基因多态性的初步研究
Polymorphism of antigens PE 35 and PPE 68 genes of Mycobacterium tuberculosis
疾病监测, 2015, 30(7): 539-545
Disease Surveillance, 2015, 30(7): 539-545
10.3784/j.issn.1003-9961.2015.07.005

文章历史

收稿日期:2015-02-11
结核分枝杆菌抗原PE 35和PPE 68基因多态性的初步研究
仇艳1, 2, 蒋毅1, 2, 刘海灿1, 2, 李马超1, 2, 万康林1, 2     
1. 中国疾病预防控制中心传染病预防控制所, 传染病预防控制国家重点实验室, 北京 102206;
2. 感染性疾病诊治协同创新中心, 杭州 310003
摘要目的 为研究结核分枝杆菌RD1区编码PE/PPE蛋白家族中的PE 35和PPE 68基因的多态性,并评估其基因多态性对免疫识别的影响。方法 笔者从本实验室保存的来自国内不同地区的2346株结核分枝杆菌复合群临床分离株中,选取代表不同Spoligotyping类型的161 株临床分离结核分枝杆菌菌株,基因扩增后进行序列比对,比较其在T细胞抗原表位的差异。结果 在161 株菌株中,23 株显示有PE 35基因多态性,8 株有PPE 68基因多态性。在PE 35基因中,除了菌株JL06018的突变外的其余突变都导致其T细胞抗原表位的改变。PPE 68基因中,62个T细胞表位中有58个(占93.5%)发生了改变。非北京家族PE 35基因的突变率明显高于北京家族 (P<0.05),而PPE 68基因突变虽是非北京家族高于北京家族,但差异无统计学意义。结论 PE 35和PPE 68的编码基因存在明显的多态性,并很有可能导致蛋白功能的明显改变。此外,这2个蛋白的人类T细胞表位是高度可变的,提示这2种蛋白可能通过参与分化选择来逃避宿主免疫。北京家族可能比非北京家族更容易被宿主T细胞识别,进而更容易传播。
关键词多态性    结核分枝杆菌    PE 35    PPE 68    
Polymorphism of antigens PE 35 and PPE 68 genes of Mycobacterium tuberculosis
QIU Yan1, 2, JIANG Yi1, 2, LIU Hai-can1, 2, LI Ma-chao1, 2, WAN Kang-lin     
1. State Key Laboratory for Communicable Disease Prevention and Control, Institute for Communicable Disease Prevention and Control, Chinese Center for Disease Control and Prevention, Beijing 102206, China;
2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou 310003, Zhejiang, China
Abstract:Objective To understand the polymorphism of antigens PE 35 and PPE 68 genes of Mycobacterium tuberculosis and evaluate the influence of the polymorphism on immune recognition. Methods A total of 161 clinical M. tuberculosis strains isolated in China were selected and comparative sequence analysis of PE 35 and PPE 68 genes was conducted after amplification. Results Among the 161 M. tuberculosis strains, 23 showed the polymorphism in the gene sequence of PE 35 and 8 showed the polymorphisms in the gene sequence of PPE 68. The mutations in all the PE 35 gene, except that in JL06018, influenced the T cell epitopes. In PPE 68 gene, 58 of 62 T cell epitopes changed, accounting for 93.5%. The mutation rate of PE 35 gene in non-Beijing family strains was significantly higher than that in Beijing family strains (P<0.05). Although the mutation rate of PPE 68 gene in non-Beijing family strains was higher than that in Beijing family strains, the difference had no statistical significance. Conclusion Polymorphisms were observed in PE 35 and PPE 68 genes, which might resulting in the drastic changes in the functions of these proteins and influence the virulence of the strains. Beijing family strains are more likely to be recognized than non-Beijing family strains by host T cell in PE 35 gene, which facilitates the spread of Beijing family strains.
Key words: Polymorphism    Mycobacterium tuberculosis    PE 35    PPE 68    

比较基因组学已发现结核分枝杆菌中的有些基因区在牛分枝杆菌及卡介苗中缺失,也就是RDI区 和RD2区[1]。RD1区(Rv3871~Rv3879c)是在BCG菌株中缺失的区域[2]。因而该区编码的蛋白在研究结核分枝杆菌的毒力及致病方面有重要作用。RD1区编码的蛋白因其不与BCG有交叉反应可用于结核分枝杆菌的T细胞反应检测。在这些抗原中,ESAT6 (ESXA,Rv3875),CFP10 (ESXB,Rv3874),PPE 68 (Rv3873)及 PE 35 (Rv3872)是与机体免疫反应相关的[3, 4, 5]。以往对ESAT6和CFP10已有较为详细的研究,它们是结核分枝杆菌的主要毒力因子[6, 7]及结核病诊断的重要候选抗原[8]PE 35和PPE 68属于PE/PPE蛋白家族并与机体的免疫反应相关。PE/PPE蛋白由结核分枝杆菌细胞壁分泌或者与其相关,并调节宿主-病原相互作用[9, 10, 11],而 PE 35和 PPE 68与结核感染后的细胞免疫反应有关。

有研究表明结核分枝杆菌PE/PPE基因有高度序列变异[12, 13],进一步提示其可能参与抗原变异。为了更好地阐述RD1区PE-PPE基因家族PE 35和PPE 68基因的多态性,及基因突变对免疫识别的影响,我们从中国选取了161 株临床分离株;进行基因扩增,然后序列比对。分析其在T细胞抗原表位的多态性及有可能产生的功能改变。

1 材料与方法 1.1 菌株来源,培养及DNA提取

实验用的161 株菌株从中国疾病预防控制中心(CDC)传染病预防控制所结核室保存的2346株国内结核分枝杆菌复合群临床分离株中挑选,考虑到北京家族菌株在中国的优势性,研究挑选了北京家族菌株82 株和非北京家族菌株79 株,82 株北京家族菌株从1738 株的北京家族菌株中随机挑选,79 株非北京家族菌株从608 株非北京家族菌株中挑选,同时尽可能包括不同省份不同Spoligtyping 类型的菌株。这些菌株囊括了所有已经在我国发现的Spoligtyping 类型。161 株菌株的省份来源及Spoligtyping分型数据见表 12

表 1 不同省份的菌株数 Table 1 M. tuberculosis strains isolated from different provinces
省份菌株数
安徽省11
陕西省16
北京市11
福建省24
甘肃省12
广西壮族自治区29
四川省1
河南省12
湖南省7
西藏自治区4
新疆维吾尔自治区11
吉林省12
浙江省11
表 2 各型别的菌株数 Table 2 Type specific M. tuberculosis strains
Spoligotyping菌株数
Beijing82
T12
U27
MANU10
Haarlem4
EAI2
LAM2
S1
CAS3
new18

采用标准的罗氏培养基培养菌株3~4 周。用生理盐水从培养基的斜面上洗脱菌体,80 ℃孵育30 min灭活,离心收集菌体,用400 μl TE悬菌,于沸水中煮沸30 min,12 000 r/min 离心4 min,取上清即为DNA模板,-20 ℃保存备用。

1.2 目的基因引物及序列测定

引物是由DNAStar 软件根据H37Rv 标准菌株的基因序列设计获得:PE 35:5′-GTA ATC GAG TTC GGG CAA TG-3′,5′-AGG CTT CTC CCA GAG AGT T-3′; PPE 68:5′-GAC ATT GGC ACG CAA GTG AG-3′,5′-TAG CGG CAT CGG TCT TCA TC-3′。

PCR反应体系为:模板DNA(500 pg)1 μl,10×PCR buffer 10 μl,引物 100 nmol/L,dNTP 混合液 200 μmol/L,DNA Taq 酶(TaKaRa) 0.5 U,ddH2O 加至20 μl。PCR反应条件为:94 ℃预变性5 min,94 ℃变性45 s,62 ℃退火45 s,72 ℃延伸1 min,35个循环,72 ℃终延伸10 min。

阳性对照使用结核H37Rv菌株的500 pg DNA。阴性对照为不加模板DNA,只有PCR试剂。PCR产物的有无及大小用2%的琼脂糖凝胶电泳鉴定。所有的聚合酶链反应(PCR)实验都至少进行2次以确认实验的重复性。扩增产物由北京擎科生物技术有限公司进行双向测序和拼接。

1.3 PE 35和PPE 68基因中TB细胞表位区的确定

检索免疫表位数据库(Immune Epitope Database,IEDB),获得结核分枝杆菌中存在的T细胞抗原表位序列,再与H37Rv株的PE 35和PPE 68基因序列进行比对,最终获得在该基因中存在的人T细胞抗原表位信息,见表 3

表 3 PE 35和PPE 68人类T细胞表位氨基酸的改变(1) Table 3 Changes of amino acid in T cell epitopes of PE 35 and PPE 68 genes
表位ID表位肽序列Rv 位置基因氨基酸改变
注:(1)以H37Rv菌株的PE 35和PPE 68作为标准序列;阴影部分表示氨基酸改变的位置。
144881EGIQLLASNASAQRv3872PE 35GCT(A)-CCT(P); Frameshift
183AAGSSATGGAAPVGAGAMGQGAQSGRv3873PPE 68CCG(P)-TCG(S); Frameshift
191AAGWQTLSAALDAQAVELTARLNSLRv3873PPE 68Frameshift(5bp deletion)
265AALAMEVYQAETAVNTLFRv3873PPE 68ATG(M)-ATA(I);Frameshift
2434ALAMEVYQAETAVNTLFEKLEPMASRv3873PPE 68ATG(M)-ATA(I);Frameshift
2922ALTEMDYFIRMWNQAALAMEVYQAERv3873PPE 68ATG(M)-ATA(I);Frameshift
3098AMGQGAQSGGSTRPGLVARv3873PPE 68Frameshift
4186ARLMAGAGPAPMLAAAAGRv3873PPE 68No change
4187ARLMAGAGPAPMLAAAAGWQTLSAARv3873PPE 68No change
4727ASQSTTNPIFGMPSPGSSTPVGQLPRv3873PPE 68Frameshift
4969ATGGAAPVGAGAMGQGAQRv3873PPE 68CCG(P)-TCG(S);Frameshift
5063ATNFFGINTIPIALRv3873PPE 68Frameshift
11486EEAAQMGLLGTSPLSNHPRv3873PPE 68Frameshift
14339ETAVNTLFEKLEPMASILRv3873PPE 68Frameshift
15812FFGINTIPIARv3873PPE 68Frameshift
16010FGMPSPGSSTPVGQLPPARv3873PPE 68Frameshift
18685GAMGQGAQSGGSTRPGLVAPAPLAQRv3873PPE 68Frameshift
18776GASQSTTNPIFGMPSPGSRv3873PPE 68Frameshift
19860GGGSDKALAAATPMVVWLQTASTQARv3873PPE 68Frameshift
20016GGSGPSAGAGLLRAESLPRv3873PPE 68Frameshift
20048GGTGGGNPADEEAAQMGLRv3873PPE 68Frameshift
20997GLLGTSPLSNHPLAGGSGPSAGAGLRv3873PPE 68Frameshift
21179GLVAPAPLAQEREEDDEDDWDEEDDRv3873PPE 6818bp deletion; Frameshift
21707GPMQQLTQPLQQVTSLFSRv3873PPE 68GTG(V)-CTG(L); Frameshift
22351GSGPSAGAGLLRAESLPGAGGSLTRRv3873PPE 68Frameshift
22531GSSTPVGQLPPAATQTLGQLGEMSGRv3873PPE 68Frameshift
22657GTGGGNPADEEAAQMGLLGTSPLSNRv3873PPE 68Frameshift
29846KALAAATPMVVWLQTASTRv3873PPE 68Frameshift
35251LDPGASQSTTNPIFGRv3873PPE 68GAT(D)-AAT(N); Frameshift
35652LEPMASILDPGASQSTTNRv3873PPE 68GAT(D)-AAT(N); Frameshift
35819LFEKLEPMASILDPGASQSTTNPIFRv3873PPE 68GAT(D)-AAT(N); Frameshift
37727LLRAESLPGAGGSLTRTPRv3873PPE 68Frameshift
38448LPEIAANHITQAVLTATNRv3873PPE 68Frameshift
38492LPGAGGSLTRTPLMSQLIEKPVAPSRv3873PPE 68Frameshift
39817LTATNFFGINTIPIARv3873PPE 68Frameshift
41291MDYFIRMWNQAALAMEVYRv3873PPE 68ATG(M)-ATA(I); Frameshift
42095MLWHAMPPELNTARLMAGRv3873PPE 68No change
46130NTIPIALTEMDYFIRMWNRv3873PPE 68Frameshift
48567PMLAAAAGWQTLSAALDARv3873PPE 68No change
49875PVGQLPPAATQTLGQLGERv3873PPE 68Frameshift
50366QATAQAAAYTQAMATTPSLPEIAANRv3873PPE 68Frameshift
51367QLIEKPVAPSVMPAAAAGSSATGGARv3873PPE 68Frameshift
52167QQVTSLFSQVGGTGGGNPRv3873PPE 68GTG(V)-CTG(L); Frameshift
52556QTLGQLGEMSGPMQQLTQRv3873PPE 68Frameshift
60492SQLIEKPVAPSVMPAAAARv3873PPE 68Frameshift
62250SVMPAAAAGSSATGGAAPRv3873PPE 68CCG(P)-TCG(S); Frameshift
64822TLGQLGEMSGPMQQLTQPLQQVTSLRv3873PPE 68Frameshift
65054TLSAALDAQAVELTARLNRv3873PPE 68Frameshift (5bp deletion)
65767TPSLPEIAANHITQAVLTATNFFGIRv3873PPE 68Frameshift
65912TQPLQQVTSLFSQVGGTGGGNPADERv3873PPE 68GTG(V)-CTG(L); Frameshift
66074TRPGLVAPAPLAQEREEDRv3873PPE 68Frameshift
66364TSPLSNHPLAGGSGPSAGRv3873PPE 68Frameshift
68284VELTARLNSLGEAWTGGGRv3873PPE 68Frameshift(5 bp deletion)
68285VELTARLNSLGEAWTGGGSDKALAARv3873PPE 68Frameshift (5bp deletion)
69128VITMLWHAMPPELNTARLMAGAGPARv3873PPE 68Frameshift
69795VLTATNFFGINTIPIALTRv3873PPE 68Frameshift
69796VLTATNFFGINTIPIALTEMDYFIRRv3873PPE 68Frameshift
71944VWLQTASTQAKTRAMQATRv3873PPE 68Frameshift
71945VWLQTASTQAKTRAMQATAQAAAYTRv3873PPE 68Frameshift
73833YFIRMWNQAALAMEVRv3873PPE 68ATG(M)-ATA(I); Frameshift
121011TPMVVWLQTASTQAKTRRv3873PPE 68Frameshift
144907IALTEMDYFIRMWNQAALAMEVYRv3873PPE 68ATG(M)-ATA(I); Frameshift
144964TNFFGINTIPIALTRv3873PPE 68Frameshift
1.4 序列比对及数据分析

PCR产物的序列由ABI 3730xl DNA分析软件分析。序列连接用的是ClustalW软件[14]。序列比对和翻译使用Bioedit 7.2.3软件。SPSS 14.0 软件用于χ2分析,P<0.05为差异具有统计学意义。

2 结果 2.1 基因突变

经扩增PE 35和PPE 68基因序列及序列比对,所有的161 株菌株均获得相应的PCR产物。在这161 株菌株中,23 株显示有PE 35基因多态性,见图 1,8 株菌有PPE 68基因多态性。对于PE 35,21 株菌有A碱基缺失,另外2 株有非同义突变。而PPE 68,2 株有2种不同的缺失,6 株显示有5种非同义突变,见图 12

图 1 161株菌株PE 35抗原的序列比对结果 Fig. 1 Sequence alignment of PE 35 gene of 161 M. tuberculosis strains 注:23 株有基因水平的改变;另外138 株与H37Rv菌株相同。T细胞抗原表位区域在H37Rv中标出。阴影部分表示氨基酸改变的位置。
图 2 161株菌株PPE 68抗原的序列比对 Fig. 2 Sequence alignment of PPE 68 gene of 161 M. tuberculosis strains 注:8 株菌在基因水平有改变,另外153 株与标准株H37Rv相同,阴影部分表示氨基酸改变的位置。
2.2 蛋白水平的改变

表 2表 4显示了PE 35和PPE 68 AA碱基的改变及位置。所有改变均导致氨基酸的改变。PE 35基因中有21 株菌株A碱基的缺失引起移码突变,进而导致蛋白翻译的提前终止,必然会影响蛋白功能。PPE 68基因中,菌株HeN06041 5 bp的缺失也导致了蛋白翻译的提前,可能对蛋白的功能有明显的影响。菌株HuN06004 18 bp的缺失引起PPE 68蛋白6个氨基酸的缺失。

2.3 突变菌株的Spoligotyping分型

对于PE 35基因,23 株突变菌株包括4 株北京家族,13 株U家族,3 株T家族,1 株MANU家族和2 株新型家族。2 株非同义突变的菌株均是北京家族。对于PPE 68,8 株突变菌株包括2 株北京家族菌株,2 株T家族菌株,2 株EAI家族菌株,1株U家族菌株和1 株 MANU家族菌株。2株EAI家族菌株,FJ06051和 FJ05406在PPE 68基因229(V-L)氨基酸位置都有相同的突变,这很有可能代表EAI家族一种特殊的突变。菌株HuN06004在PE 35和PPE 68基因均有呈现出基因多态性。

非北京家族PE 35基因的突变率明显高于北京家族(表 4P<0.05),而PPE 68基因突变虽是非北京家族高于北京家族,但差异无统计学意义。

表 4 北京家族与非北京家族菌株突变的比较 Table 4 Comparison of mutation between Beijing family strains and non-Beijing family strains
突变菌株北京家族非北京家族
注:北京家族与非北京家族的比较: (1)P=0.001;(2)P=0.253。
PE 35(1)419
PPE 68(2)26
全部8279
2.4 T细胞表位的改变

根据免疫细胞表位数据库(IEDB)[15]PE 35和PPE 68共有63个人类T细胞表位。表 4显示了这两种抗原的T细胞表位的改变。在PE 35基因中,除了菌株JL06018的突变外的其余突变都导致其T细胞抗原表位的改变。在PPE 68基因,62个T细胞抗原表位覆盖了所有的基因序列,没有非T细胞抗原表位区,这也提示PPE 68对于T细胞感染免疫反应的重要性。PPE 68基因中,62个T细胞表位中有58个(占93.5%)发生了改变。

菌株HeN06041 5 bp的缺失引起了PPE 68蛋白编码位置的移码突变,这导致相应的T细胞表位的改变,包括IEDB_ID 191,65054,68284,68285及其下游的表位。菌株HuN06004 18 bp的缺失导致6个氨基酸的缺失,引起表位IEDB_ID21179的改变。

3 讨论

本研究根据地区和基因型分布等,随机从国内不同地区来源,不同Spoligtyping基因型的2346 株临床分离的结核分枝杆菌复合群菌株中,选取了161 株,采用PCR和DNA测序的方法,对PE 35和PPE 68基因及其T细胞抗原表位进行了分析。菌株代表性好,得出的结果有较高的可信度。

基因组学及生物信息学分析表明RD 1区编码的蛋白是参与分泌ESAT 6和CFP 10分泌系统的一部分,即EAST-6分泌系统(ESX-1)。PE 35,CFP 10及ESAT 6是刺激人类外周血单核淋巴细胞参与保护性Th 1细胞反应检测的重要抗原,也就是刺激增殖并产生γ-干扰素的抗原[4]PPE 68是主要的细胞壁相关抗原[16],并与RD 1区的蛋白Rv 3866,Rv 3868,CFP10,ESAT-6形成免疫复合物[17, 18] 。最近研究表明PE 35及PPE 68在RD 1-相关的发病机制中发挥重要作用,可能促成及维持M.tb的感染[19]。在161 株菌种,有14.3%的菌株在PE 35基因中有一A碱基的缺失导致蛋白编码的提前终止,从全长99氨基酸缩短为16个氨基酸。此改变可能会完全改变蛋白结构,同时可引起PE 35功能的缺失。菌株HeN 06041 PPE 68蛋白5 bp的缺失同样引起蛋白编码提前终止,这一改变也有可能对蛋白功能有明显的影响。考虑到PE 35和PPE 68与结核分枝杆菌的毒力和免疫相关,PE 35和PPE 68蛋白的突变菌株很有可能会引起菌株毒力的改变。此外,因PE 35对EsxA/B分泌系统及RD 1的毒力是必需的,PE 35的无效突变(即蛋白编码提前终止的菌株)有可能会影响ESX-1抗原的分泌(CFP 10,ESAT 6)[20]。而PPE 68对 ESX-1抗原的分泌具有调节作用[20],因此我们推测PPE 68的突变可能会影响ESAT 6,CFP 10的分泌。下一步工作需进一步研究PPE 68和PE 35蛋白的突变菌株与野生菌株抗原毒力的差异。

PE/PPE基因具是多变的,并编码细胞表面蛋白,该基因可能参与了抗原变异并被认为是逃避宿主免疫的主要因子[21]。2010年,Comas等[22]研究报道结核分枝杆菌人类T细胞抗原表位的进化是高度保守的,因此推测其缺少抗原变异及免疫逃逸。但因由于一些技术原因,该研究没有包括PE/PPE基因。在本研究中,根据免疫表位数据库(IEDB),在PE 35和PPE 68中共有63个人类T细胞表位[15]。在所有的菌株中,63个T细胞抗原表位中有59个(占 93.7%)发生了改变,这些T细胞抗原表位大量氨基酸的改变可能反映出正在进行的免疫逃逸。本研究的实验数据说明一些PE/PPE家族蛋白基因有高度的序列变异,这有可能引起相关抗原的改变,进而导致免疫逃逸。

PE 35的基因突变在非北京家族中明显高于北京家族,这意味着北京家族菌株比非北京家族菌株更为保守,非北京家族更有可能应对宿主的免疫压力而产生抗原变异。这与之前的研究结果一致,即北京家族菌株相对于其他的结核分枝杆菌菌株基因更为保守,比其他型别的结核分枝杆菌有选择优势[23]。有证据显示T细胞的识别可能会直接促成人与人之间结核分枝杆菌复合群的传播[24]。本研究结果提示北京家族可能比非北京家族更容易被宿主T细胞识别,进而更容易传播。

据报道PE 35有可能作为血清学诊断的抗原[25]。本研究发现PE 35具有相当多的氨基酸改变,因此在选用PE 35作为血清学诊断候选抗原时应充分考虑到菌株突变的影响。

参考文献
[1] Liu XQ, Dosanjh D, Varia H, et al. Evaluation of T-cell responses to novel RD1-and RD2-encoded Mycobacterium tuberculosis gene products for specific detection of human tuberculosis infection[J]. Infect Immun,2004,72(5):2574-2581.
[2] Mahairas GG, Sabo PJ, Hickey MJ, et al. Molecular analysis of genetic differences between Mycobacterium bovis BCG and virulent M. bovis[J]. J Bacteriol,1996,178(5):1274-1282.
[3] Hanif SN, El-Shammy AM, Al-Attiyah R, et al. Whole blood assays to identify Th1 cell antigens and peptides encoded by Mycobacterium tuberculosis-specific RD1 genes[J]. Med Princ Pract,2008,17(3):244-249.
[4] Mustafa AS, Al-Attiyah R, Hanif SN, et al. Efficient testing of large pools of Mycobacterium tuberculosis RD1 peptides and identification of major antigens and immunodominant peptides recognized by human Th1 cells[J]. Clin Vaccine Immunol,2008,15(6):916-924.
[5] Hanif SN, Al-Attiyah R, Mustafa AS. Species-specific antigenic Mycobacterium tuberculosis proteins tested by delayed-type hypersensitivity response[J]. Int J Tuberc Lung Dis,2010,14(4):489-494.
[6] Pym AS, Brodin P, Brosch R, et al. Loss of RD1 contributed to the attenuation of the live tuberculosis vaccines Mycobacterium bovis BCG and Mycobacterium microti[J]. Mol Microbiol,2002,46(3):709-717.
[7] Lewis KN, Liao R, Guinn KM, et al. Deletion of RD1 from Mycobacterium tuberculosis mimics bacille Calmette-Guerin attenuation[J]. J Infect Dis,2003,187(1):117-123.
[8] Malaghini M, Thomaz-Soccol V, Probst CM, et al. Recombinant antigen production for assays of intradermoreaction for diagnosis and surveillance of tuberculosis[J]. J Biotechnol,2011,156(1):56-58.
[9] Mukhopadhyay S, Balaji KN. The PE and PPE proteins of Mycobacterium tuberculosis[J]. Tuberculosis,2011,91(5):441-447.
[10] Sampson SL. Mycobacterial PE/PPE proteins at the host-pathogen interface[J]. Clin Dev Immunol,2011,2011:497203.
[11] Akhter Y, Ehebauer MT, Mukhopadhyay S, et al. The PE/PPE multigene family codes for virulence factors and is a possible source of mycobacterial antigenic variation: perhaps more?[J]. Biochimie,2012,94(1):110-116.
[12] McEvoy CRE, Cloete R, Muller B, et al. Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints[J]. PLoS One,2012,7(4):e30593.
[13] Copin R, Coscollá M, Seiffert SN, et al. Sequence diversity in the pe_pgrs genes of Mycobacterium tuberculosis is independent of human T cell recognition[J]. mBio,2014,5(1):e00960-13.
[14] Larkin MA, Blackshields G, Brown NP, et al. Clustal W and Clustal X version 2.0[J]. Bioinformatics,2007,23(21):2947-2948.
[15] Kim Y, Ponomarenko J, Zhu Z, et al. Immune epitope database analysis resource[J]. Nucleic Acids Res,2012,40(Web Server issue):W525-W530.
[16] Okkels LM, Brock I, Follmann F, et al. PPE protein (Rv3873) from DNA segment RD1 of Mycobacterium tuberculosis: strong recognition of both specific T-cell epitopes and epitopes conserved within the PPE family[J]. Infect Immun,2003,71(11):6116-6123.
[17] Okkels LM, Andersen P. Protein-protein interactions of proteins from the ESAT-6 family of Mycobacterium tuberculosis[J]. J Bacteriol,2004,186(8):2487-2491.
[18] Teutschbein J, Schumann G, Mollmann U, et al. A protein linkage map of the ESAT-6 secretion system 1 (ESX-1) of Mycobacterium tuberculosis[J]. Microbiol Res,2009,164(3):253-259.
[19] Tiwari B, Soory A, Raghunand TR. An immunomodulatory role for the Mycobacterium tuberculosis region of difference 1 locus proteins PE 35 (Rv3872) and PPE 68 (Rv3873)[J]. Febs J,2014,281(6):1556-1570.
[20] Brodin P, Majlessi L, Marsollier L, et al. Dissection of ESAT-6 system 1 of Mycobacterium tuberculosis and impact on immunogenicity and virulence[J]. Infect Immun,2006,74(1):88-98.
[21] Vordermeier HM, Hewinson RG, Wilkinson RJ, et al. Conserved immune recognition hierarchy of mycobacterial PE/PPE proteins during infection in natural hosts[J]. PLoS One,2012,7(8):e40890.
[22] Comas I, Chakravartti J, Small PM, et al. Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved[J]. Nat Genet,2010,42(6):498-503.
[23] Parwati DI, van Crevel R, van Soolingen D. Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains[J]. Lancet Infect Dis,2010,10(2):103-111.
[24] Kwan CK, Ernst JD. HIV and tuberculosis: a deadly human syndemic[J]. Clin Microbiol Rev,2011,24(2):351-376.
[25] Mukherjee P, Dutta M, Datta P, et al. The RD1-encoded antigen Rv3872 of Mycobacterium tuberculosis as a potential candidate for serodiagnosis of tuberculosis[J]. Clin Microbiol Infect,2007,13(2):146-152.