文章信息
- 韩志英, 李静, 顾凯侃, 孙果梅, 江渊, 张阳奕, 徐飚.
- Han Zhiying, Li Jing, Gu Kaikan, Sun Guomei, Jiang Yuan, Zhang Yangyi, Xu Biao.
- 上海市静安区2010-2015年肺结核传播状况及影响因素分析
- Recent transmission of pulmonary tuberculosis and its influencing factors in Jing'an district, Shanghai, 2010-2015
- 中华流行病学杂志, 2018, 39(10): 1339-1345
- Chinese Journal of Epidemiology, 2018, 39(10): 1339-1345
- http://dx.doi.org/10.3760/cma.j.issn.0254-6450.2018.10.010
-
文章历史
收稿日期: 2018-03-13
2. 200336 上海市疾病预防控制中心结核病检测实验室;
3. 200032 上海, 复旦大学公共卫生学院教育部公共卫生安全重点实验室
2. Tuberculosis Laboratory, Shanghai Municipal Center for Disease Control and Prevention, Shanghai 200336, China;
3. School of Public Health, Key Laboratory of Public Health Safety of Ministry of Education, Fudan University, Shanghai 200032, China
结核病目前已成为全球受关注的重大公共卫生问题之一。据WHO全球结核病报告估算[1],2015年全球登记的新发结核病患者约1 040万,其中约48万为同时耐异烟肼和利福平的耐多药结核病(multiple drug-resistant tuberculosis,MDR-TB)患者。我国新发肺结核人数为91.8万(67/10万),MDR- TB患者约5.8万。2015年上海市报告肺结核病例6 699例,报告发病率为27.62/10万,总耐药率为22.90%,总耐多药率为6.50%。2010-2015年静安区结核病报告发病率居本区甲乙类传染病报告发病前3位。掌握结核病的传播规律,有助于控制结核病的传播和流行。
由于结核病潜伏期长,传统流行病学方法较难深入揭示其发病原因及传播状况,结合传统流行病学建立的分子流行病学能更好地解释结核病的传播特点[2-4]。目前,结核分枝杆菌可变数目串联重复序列(Mycobacterial interspersed repetitive unit-variable number tandem repeat,MIRU-VNTR)多态性分析方法已广泛用于研究鉴别内源性复发和近期外源性感染、耐药菌株的传播等。Supply等[5]2006年提出的15个或24个位点可变数目串联重复(VNTR)是应用最多的方案,但该方案对我国高度流行的北京家族菌株的分辨率较低[6]。本研究采用Luo等[6]2014年建立的一套位点数少且对北京家族菌株分辨率较高的12个位点MIRU-VNTR基因型分型法,结合现场流行病学调查从人群角度分析结核分枝杆菌成簇特征及其影响因素,了解静安区结核病的传播状况。
对象与方法1.样本来源:选择上海市CDC菌株库保存的2010-2015年静安区定点医院治疗管理的肺结核全部菌株92株,经菌型鉴定、剔除重复菌株和复苏失败菌株后获得有效菌株80株,进行药物敏感性试验和MIRU-VNTR基因型分型。
2.流行病学调查:采用电话访谈、上门调查与查阅病史资料、结核病登记管理资料相结合的方式收集患者基本情况、临床特征、治疗情况、副作用、转归等。调查对象首选患者本人,如果患者因某些原因无法完成调查,则选择患者的密切接触人员(如配偶、子女、父母等)。调查前均向患者或患者的密切接触人员告知并获得知情同意。结核病患者现场调查按《中国结核病防治规划实施工作指南》(2008版)相关规定实施。
对成簇病例开展现场流行病学调查,使用自制的调查问卷,收集成簇病例的密切接触者、人际关系、居住地点、常去地点、医院就诊情况等信息,以确定成簇病例之间的相互关联性。如果成簇病例之间为密切接触者或相互认识、或者居住在同一地点,则认为该组病例之间存在“确定的流行病学联系”;如果成簇病例之间居住地点较近、或有共同的常去地点,则认为该组病例之间存在“可能的流行病学联系”。
3.药物敏感性试验:按照中国防痨协会《结核病实验室检验规程》[7]进行,采用《结核病耐药监测指南》[8]推荐的比例法。采用L-J培养基,进行4种一线抗结核药物敏感性试验。抗结核药物的终浓度:异烟肼0.2 μg/ml、链霉素4.0 μg/ml、利福平40.0 μg/ml、乙胺丁醇2.0 μg/ml。耐药判定:根据菌落在培养基上生长情况,若含药培养基上生长的菌落数/对照培养基上生长的菌落数>1%,则认为受试菌对该抗结核药物耐药。
4. MIRU-VNTR基因型分型:
(1)结核分枝杆菌DNA制备:采用煮沸裂解法进行提取[9],从罗氏培养基上刮取1~2个接种环阳性培养物,重悬于200 μl TE(Tris-EDTA,三羟甲基氨基甲烷-乙二胺四乙酸)溶液中,80 ℃灭活30 min。灭活后的菌株DNA 100 ℃水浴10 min煮沸裂解,立即置冰上2 min,13 000×g离心15 min后,取上清置于另一无菌1.5 ml离心管中,-20 ℃保存备用。
(2)MIRU-VNTR基因型分型:采用12个位点MIRU-VNTR基因分型法[10]:包括9个位点:QUB11b、QUB18、Mtub21、Miru26、QUB26、Mtub04、Miru31、Miru40、VNTR2372和3个高变位点:VNTR3820、3232、4120。利用北京康为世纪生物科技有限公司生产的结核病鉴定试剂盒Ⅱ进行结核分枝杆菌VNTR多态性分析。PCR反应体系和条件参考试剂盒说明书。取1 μl PCR扩增产物用1%的琼脂糖凝胶进行电泳。电泳读数:根据各样本扩增条带与Marker相对位置,利用各VNTR位点重复单元读数表及VNTR各位点重复单元电泳图,计算出各菌株在该位点的重复单元数。以H37Rv标准株作为质控菌株。
5.定义:“簇(cluster)”菌株是指从不同病例分离得到的≥2个具有相同VNTR基因型的菌株,代表着该组病例很可能是由于“近期感染”而发病;而将从不同患者分离到的具有特异VNTR基因型的菌株定义为单一基因型菌株(unique),代表该病例很可能是由于潜伏感染的“内源性复燃”而发病[10-13]。成簇比(近期传播率):反映结核病近期传播情况,估算公式为C=n/N,其中C代表成簇比,n代表簇菌株数,N代表总菌株数(即簇菌株数+单一基因型菌株数)[10]。
初治患者是指从未用过抗结核药物或用抗结核药物不超过28 d的患者。复治患者是指接受抗结核治疗超过1个月的患者(包括复发、初治失败和返回的患者)。耐多药是指患者感染的结核分枝杆菌至少同时对异烟肼、利福平耐药。
6.统计学分析:采用EpiData 3.0软件建立数据库对资料进行双录入及一致性检验。应用SPSS 16.0软件进行数据整理及统计学分析。分类变量资料采用百分比进行描述,两样本率的比较采用χ2检验或Fisher确切概率法,多因素分析采用logistic回归方法,对80株菌株检测的12个MIRU-VNTR位点的重复次数采用Bio Numerics 5.0软件进行聚类分析。检验水准为α=0.05,双侧检验。
结果1.基本情况:共登记173例肺结核患者,其中结核分枝杆菌培养阳性92例,最终纳入获得有效菌株的患者80例。患者年龄17~89(47.15±20.02)岁,M=52岁;男性69例(86.25%),女性11例(13.75%);初治患者63例(78.75%),复治患者17例(21.25%)。职业以退休为主(25例,31.25%),待业(24例,30.00%)。
2.耐药情况:80株结核分枝杆菌菌株中,57株对全部4种一线抗结核药物(异烟肼、链霉素、利福平和乙胺丁醇)敏感,23株对至少1种抗结核药物耐药,总耐药率为28.75%(23/80),初治患者耐药率28.57%(18/63),复治患者耐药率为29.41%(5/17),差异无统计学意义(χ2=0.005,P=0.946)。耐多药结核分枝杆菌为13株,总耐多药率为16.25%(13/80),初治患者耐多药率为14.29%(9/63),复治患者耐多药率为23.53%(4/17),差异无统计学意义(χ2=0.841,P=0.359)。
3. MIRU-VNTR基因型分型结果及成簇情况(近期传播率):80株结核分枝杆菌菌株的12个位点MIRU-VNTR基因型分型结果和聚类分析树状结构,见图 1。VNTR基因型分析,共鉴定出65种基因型。58(72.50%)例患者的分离株为单一基因型,22例为成簇菌株,成簇百分比为27.50%(22/80)。22例成簇菌株共分7个簇,其中6个簇(85.71%)分别只包含2例菌株,1个簇包含10例菌株,均为MDR-TB患者。所有13例MDR-TB菌株共鉴定出4种VNTR基因型,其中10例菌株形成1个簇,耐多药菌株的成簇百分比为76.92%。耐多药菌株成簇百分比明显高于全敏感结核菌株成簇百分比(χ2=16.739,P<0.001)。
4.菌株成簇病例的分布特征:成簇病例中以初治(81.82%,18/22)、男性(86.36%,19/22)、<65岁(81.82%,18/22)的患者居多,职业主要为退休(31.82%,7/22)和待业(27.27%,6/22)。成簇病例中耐药比例为50.00%(11/22),MDR-TB为45.45%(10/22)。
5.近期传播的影响因素分析:
(1)成簇的单因素分析:有合并症患者成簇百分比明显高于无合并症患者(χ2=9.881,P=0.002)、耐药患者明显高于全敏感患者(χ2=6.689,P=0.010),MDR-TB患者明显高于非MDR-TB患者(χ2=16.739,P=0.000),痰菌涂阳患者明显高于涂阴患者(χ2=8.188,P=0.004),有既往住院史患者明显高于无既往住院史患者(χ2=6.180,P=0.013)。未发现年龄、性别、职业种类、文化程度、婚姻状况、化疗史、肺部有无空洞、是否遵医嘱服药、平常就诊定点医院以及患病前听说过肺结核与病例菌株成簇间有统计学联系(均P>0.05)。见表 1。
(2)多因素logistic回归分析:以患者的结核分枝杆菌分离株是否成簇作为因变量,成簇项为l,非成簇项为0。将可能的危险因素(合并症、耐多药、痰涂片、既往住院史、抗结核治疗史和年龄、性别)为自变量纳入logistic回归模型,使用ENTER法分析。结果表明耐多药(OR=35.799,95%CI:4.239~302.346)和合并症(OR=7.695,95%CI:1.421~41.658)是近期传播相关的危险因素,相比于敏感菌株,MDR-TB病例更容易发生近期传播。见表 1。
6.成簇患者的现场流行病学调查:对成簇患者进行现场流行病学调查结果发现,22例成簇患者中,11例患者存在“确定的流行病学联系”,9例患者存在“可能的流行病学联系”,2例患者暂未发现流行病学联系。最大的一簇包含10例MDR-TB患者,有9例患者互相认识(确定的流行病学联系),包括2对夫妻和1对父子,并且其中8例患者居住在一个弄堂里,第9例患者居住附近200 m以内;第10例患者与其他9例患者中的7例均有在附近同一棋牌室的高频率逗留史(可能的流行病学联系)。另有1个簇的2例患者生活在一套住房内,为房东和租客关系(确定的流行病学联系)。有4个簇(各含2例)患者每簇病例居住地间的直线距离在2 km以内,每簇间的病例有长期到附近同一菜场和超市购物的经历(可能的流行病学联系)。有1个簇的2例患者之间暂未发现直接或间接的流行病学联系。22例患者的住院时间和门诊就诊时间未发现存在交集。
讨论近期传播率是评价肺结核流行和控制效果的有效指标,结核分枝杆菌中的成簇百分比直接反映当地近期传播情况。本研究发现静安区结核分枝杆菌成簇百分比为27.5%,相当于近30%的患者是由结核分枝杆菌的近期传播所致。与杨崇广[10]研究上海松江区(30.7%)、普陀区[11](26.7%)的成簇百分比基本一致,稍高于宋红焕[12]的研究结果(22.8%),低于另外一项研究结果[13]。本研究中成簇患者以退休和待业为主,没有固定于某一特定人群,与国内的一些研究发现一致[10, 13-14],而与国外的一些研究存在很大区别,其肺结核发病和传播往往集中于特定的人群,如HIV感染者[15-16]、吸毒者[17]等。但由于国内活动性肺结核病的人群数量多,即使近30%的近期传播率,导致的二次病例数量也是非常可观[18]。
研究还发现静安区结核病患者中耐多药菌株的成簇百分比明显高于全敏感菌株。虽然国外有报道由于某些耐药基因突变使分枝杆菌的传播性和致病性降低,耐药菌株成簇的可能性小于敏感菌株[16],但也有专家研究提示耐药菌株的适应能力降低可能会通过代偿性突变来抵偿,一些未发生适应性降低的耐药菌株易于发生传播[19],还有研究报道[10, 20-22],MDR-TB更容易发生近期传播。MDR-TB患者由于诊断和治疗延迟,痰菌阴转慢、治疗时间长、难治愈等因素可能使患者具有较长的传染期,基本没有隔离限制措施等管理的一些问题[23-24],更增加了其发生近期传播的风险。本研究中的家庭和社区内MDR-TB二级病例的产生在国内很少报道,这一重要发现进一步证明了耐多药结核菌株可能具有相当强的传播能力。
对影响菌株成簇因素的logistic回归分析结果表明:耐多药和合并症是肺结核近期传播的主要危险因素。本研究中成簇患者中有近一半的患者为MDR-TB,且13例MDR-TB中10例出现成簇,且集中于同一个簇。这一发现与杨崇广[10]、胡屹等[25]研究结果类似,提示在我国,MDR-TB病例间的近期传播已趋向普遍,静安区MDR-TB近期传播的现象尤为令人关注。由于有合并症的患者身体状况不良,对药物副反应明显,导致治疗效果差,排菌时间长,更易传染。因此,其病例间发生近期传播的风险也较高。建议对有合并症的结核病患者要实行个体化治疗方案,加强管理,提高患者治疗的依从性和疗效。
有研究报道在低流行地区,同簇菌株患者间的联系难以被流行病学调查发现[26-27],本研究结果发现22例成簇患者中,11例患者存在“确定的流行病学联系”,9例存在“可能的流行病学联系”。有2簇患者有室内接触史,其中1个MDR簇的10例患者有9例患者互相认识,包括2对夫妻患者和1对父子患者。提示家庭成员间接触传播和机会接触导致的传播在本区都存在,提示结核病的传播并不是主要发生在传统意义上的长时间密切接触人群中。MDR-TB簇中的10例患者在不同时间先后发病,MIRU-VNTR基因型结果显示10例患者的分离株中6人在9位点上扩增目的片段大小一致,与其他3人相差1个拷贝数,还有1人与其他9人在QUB11b拷贝数相差1,10人在3个高变位点拷贝数均一致。出现这些微小差异很有可能是结核分枝杆菌在患者体内发生了微进化,进化后的结核菌作为传染源随之再传播,使得10例患者在VNTR个别位点存在差异,后续研究将进一步开展全基因组测序。MDR- TB簇中10例患者没有发现院内互相传染导致感染的机会,互相接触传染的可能性还是在小区及其附近的棋牌室等室内活动场所。某些公共场所存在通风差、空间相对狭小、人群比较集中等特点,如棋牌室、影剧院、快餐店、健身房、超市和图书馆等,结核病患者尤其是MDR-TB患者在传染期内到这类场所活动很可能引起传播,因此应采取措施限制传染期内结核病患者的活动范围,但是目前的门诊治疗方案在这方面没有可实际操作的条例规定,只有住院治疗,才可以真正限制患者的行动,因此建议结核病患者在传染期间接受住院治疗,直到痰菌转阴。
本研究存在不足,80株结核分枝杆菌分离株的样本量较小,可能影响研究结果的把握度,有待进一步扩大样本完善现有研究结论。
总之,2010-2015年上海市静安区结核病患者尤其是MDR-TB患者存在一定比例的近期传播,耐多药和合并症是结核病近期传播的重要危险因素。建议对MDR-TB患者活动频繁的场所的密切接触者进行主动筛查,尽早发现结核病患者,及时进行有效的治疗和管理,控制MDR-TB的流行。
利益冲突: 无
[1] |
WHO. Global Tuberculosis Control: WHO Report, 2016.
|
[2] |
WHO. Global tuberculosis control: surveillance, planning, financing[R]. WHO report. WHO/HTM/TB/2006, 362. Geneva: World Health Organization, 2006. http://apps.who.int/bookorders/anglais/detart1.jsp?codlan=1&codcol=15&codcch=3659
|
[3] |
乐军, 李瑶, 谢建平, 等. 寡核苷酸表达谱揭示异烟肼耐药对结核分枝杆菌致病性的影响[J]. 中华医院感染学杂志, 2004, 14(12): 1336-1340. Yue J, Li Y, Xie JP, et al. Effect of Isoniazid Resistance on Pathogenicity of Mycobacterium tuberculosis Revealed by Oligonucleotide Microarray Profiling[J]. Chin J Nosocomiol, 2004, 14(12): 1336-1340. DOI:10.3321/j.issn:1005-4529.2004.12.005 |
[4] |
匡红, 陈庆海, 府伟灵. 分子信标探针荧光芯片检测结核分枝杆菌的实验研究[J]. 中华医院感染学杂志, 2007, 17(5): 509-512. Kuang H, Chen QH, Fu WL. Detection of Mycobacterium tuberculosis by Fluorescence Chip of Molecular Beacon Probe[J]. Chin J Nosocomiol, 2007, 17(5): 509-512. DOI:10.3321/j.issn:1005-4529.2007.05.009 |
[5] |
Supply P, Allix C, Lesjean S, et al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis[J]. J Clin Microbiol, 2006, 44(12): 4498-4510. DOI:10.1128/JCM.01392-06 |
[6] |
Luo T, Yang C, Pang Y, et al. Development of a hierarchical variable-number tandem repeat typing scheme for Mycobacterium tuberculosis in China[J]. PLoS One, 2014, 9(2): e89726. DOI:10.1371/journal.pone.0089726 |
[7] |
中国防痨协会. 结核病实验室检验规程[M]. 北京: 人民卫生出版社, 2015: 59-65. Chinese Antituberculosis Association. Tuberculosis Laboratory Inspection Specification[M]. Beijing: People's Medical Publishing House, 2015: 59-65. |
[8] |
Bustreo F, Migliori GB, Nardini S, et al. Antituberculosis drug resistance:is it worth measuring? The WHO/IUATLD Working Group on Antituberculosis Drug Resistance Surveillance[J]. Monaldi Arch Chest Dis, 1996, 51(4): 299-302. |
[9] |
沈国妙, 查佳, 徐琳, 等. 结核分枝杆菌散在分布重复单位基因型分型法的应用研究[J]. 中华结核和呼吸杂志, 2005, 28(5): 292-296. Shen GM, Zha J, Xu L, et al. Evaluation of the mycobacterial interspersed repetitive units typing as a practical approach in molecular epidemiology of Mycobacterium tuberculosis[J]. Chin J Tubercul Respirat Dis, 2005, 28(5): 292-296. DOI:10.3760/j:issn:1001-0939.2005.05.002 |
[10] |
杨崇广.中国部分地区以人群为基础的结核病分子流行病学研究[D].上海: 复旦大学, 2013. DOI: 10.7666/d.Y2703749. Yang CG. Population-based Molecular Epidemiology of Tuberculosis in multi-regions of China[D]. Shanghai: Fudan University, 2013. DOI: 10.7666/d.Y2703749. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2703749 |
[11] |
郑亦慧, 陈俊, 张宇艳, 等. 上海市普陀区结核菌近期传播特征及影响因素研究[J]. 环境与职业医学, 2014, 31(5): 368-372. Zheng YH, Chen J, Zhang YY, et al. Pattern and Impact Factors of Recent Tuberculosis Transmission in Putuo District, Shanghai[J]. J Environ Occupa Med, 2014, 31(5): 368-372. DOI:10.13213/j.cnki.jeom.2014.0082 |
[12] |
宋红焕.江苏省耐多药结核分枝杆菌分子流行病学研究[D].东南大学硕士论文, 2016. Song HH. Study on Molecular Epidemiology of Mutidrug-resistant Mycobacterium tuberculosis in Jiangsu province[D]. Southeast University, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10286-1017116856.htm |
[13] |
罗丹.广西耐药结核病流行特征及结核分枝杆菌基因分型的研究[D].广西医科大学, 2013. Luo D. Epidemic Characteristics of drug-resistant tuberculosis and genotyping of Mycobacterium tuberculosis[D]. Guangxi Medical University, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10598-1013247448.htm |
[14] |
李霞.结核病及耐药结核病的传播[D].复旦大学, 2011. Li X. Transmission of Tuberculosis and Drug Resistant Tuberculosis[D]. Fudan University, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10246-1011183905.htm |
[15] |
Alland D, Kalkut GE, Moss AR, et al. Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and convention[J]. N Engl J Med, 1994, 330(24): 1710-1716. DOI:10.1056/NEJM199406163302403 |
[16] |
Godfrey-Faussett P, Sonnenberg P, Shearer SC, et al. Tuberculosis control and molecular epidemiology in a South African gold-mining community[J]. Lancet, 2000, 356(9235): 1066-1071. DOI:10.1016/S0140-6736(00)02730-6 |
[17] |
Cronin WA, Golub JE, Lathan MJ, et al. Molecular epidemiology of tuberculosis in a low-to moderate-incidence state:are contact investigations enough?[J]. Emerg Infect Dis, 2002, 8(11): 1271-1279. DOI:10.3201/eid0811.020261 |
[18] |
Collaboration CTC. The effect of tuberculosis control in China[J]. Lancet, 2004, 364(9432): 417-422. DOI:10.1016/S0140-6736(04)16764-0 |
[19] |
Hsu AH, Lin CB, Lee YS, et al. Molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in East Tai wan[J]. Int J Tuberc Lung Dis, 2010, 14(7): 924-926. |
[20] |
王伟炳.农村地区结核病传播机制及主动发现模式的流行病学研究[D].上海: 复旦大学, 2006. Wang WB. Pattern of Tuberculosis transmission and model of active case finding by conventional and molecular epidemiological method in rural areas of China[D]. Shanghai: Fudan University, 2006. http://cdmd.cnki.com.cn/Article/CDMD-10246-2007068737.htm |
[21] |
Li X, Zhang Y, Shen X, et al. Transmission of drug-resistant tuberculosis among treated patients in Shanghai, China[J]. J Infect Dis, 2007, 195(6): 864-869. DOI:10.1086/511985 |
[22] |
Moonan PK, Teeter LD, Salcedo K, et al. Transmission of multidrug-resistant tuberculosis in the USA:a cross-sectional study[J]. Lancet Infect Dis, 2013, 13(9): 777-784. DOI:10.1016/S1473-3099(13)70128-2 |
[23] |
He GX, Wang HY, Borgdorff MW, et al. Multidrug-resistant tuberculosis, People's Republic of China, 2007-2009[J]. Emerg Infect Dis, 2011, 17(10): 1831-1838. DOI:10.3201/eid1710.110546 |
[24] |
Zhao Y, Xu S, Wang L, et al. National survey of drug-resistant tuberculosis in China[J]. N Engl J Med, 2012, 366(23): 2161-2170. DOI:10.1056/NEJMoa1108789 |
[25] |
胡屹.中国华东农村地区耐药结核病流行和传播的分子流行病学研究[D].上海: 复旦大学, 2008. Hu Y. Molecular epidemiological study on the prevalence and transmission of drug resistant TB in rural areas of Eastern China[D]. Shanghai: Fudan University, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10246-2009018572.htm |
[26] |
梅建, 沈鑫, 查佳, 等. 上海市2000-2002年91株结核分枝杆菌分子流行病学分析[J]. 中华流行病学杂志, 2005, 26(9): 707-710. Mei J, Shen X, Zha J, et al. Study on the molecular epidemiology of Mycobacterium tuberculosis in Shanghai[J]. Chin J Epidemiol, 2005, 26(9): 707-710. DOI:10.3760/j.issn:0254-6450.2005.09.021 |
[27] |
Luo T, Yang C, Peng Y, et al. Whole-genome sequencing to detect recent transmission of Mycobacterium tuberculosis in settings with a high burden of tuberculosis[J]. Tuberculosis (Edinb), 2014, 94(4): 434-440. DOI:10.1016/j.tube.2014.04.005 |