pncA、rpsA和panD基因突變對(duì)結(jié)核分枝桿菌吡嗪酰胺耐藥性的預(yù)測(cè)效能分析

趙東陽，石 潔，蘇茹月，鄭丹薇，朱巖昆，馬曉光，王少華，李 輝，孫國清，孫定勇

河南省疾病預(yù)防控制中心 鄭州 450016

吡嗪酰胺(pyrazinamide，PZA)作為一種前沿抗結(jié)核藥物在結(jié)核病一線和二線治療方案中均發(fā)揮重要作用[1]。該藥物能有效對(duì)抗體內(nèi)巨噬細(xì)胞中滯留的不能被其他的抗結(jié)核藥物所殺死的結(jié)核分枝桿菌。PZA需要在吡嗪酰胺酶(PZase)的作用下轉(zhuǎn)化為有活性的吡嗪酸發(fā)揮作用。PZase是由561個(gè)核苷酸編碼的蛋白。pncA基因突變導(dǎo)致PZase失活，是引起結(jié)核分枝桿菌對(duì)PZA耐藥的主要機(jī)制。pncA突變具有高度可變性，且散在分布于開放閱讀框及上游基因調(diào)控區(qū)域[2]。除了pncA基因突變外，編碼30S核糖體蛋白S1的基因rpsA突變后，會(huì)改變吡嗪酸的結(jié)合位點(diǎn)，也可引起PZA耐藥[3]。有研究[4]表明panD和結(jié)核分枝桿菌PZA的耐藥相關(guān)。到目前為止，rpsA和panD突變是否引起結(jié)核分枝桿菌PZA耐藥仍有爭議性，因此需要更多實(shí)驗(yàn)數(shù)據(jù)來闡明兩者與PZA耐藥的關(guān)系。

由于PZA只有在低pH值條件下才具有抗菌活性，因此在做PZA藥敏測(cè)定時(shí)需在酸性環(huán)境下進(jìn)行。且PZA藥敏測(cè)定操作較復(fù)雜，實(shí)驗(yàn)失敗可能性較大，同時(shí)結(jié)核分枝桿菌生長周期長，常規(guī)PZA藥敏檢測(cè)需要2～3周，因此近年來分子診斷成為檢測(cè)PZA耐藥的一種重要工具。本研究對(duì)152株耐多藥(MDR)結(jié)核分枝桿菌菌株進(jìn)行了PZA藥敏測(cè)定，同時(shí)對(duì)pncA、rpsA和panD突變預(yù)測(cè)PZA耐藥性的效能進(jìn)行了評(píng)價(jià)。

1 材料與方法

1.1樣本來源收集2018年1月至12月于新密市、嵩縣、扶溝縣、焦作市疾病預(yù)防控制中心，中牟縣衛(wèi)生防疫站，鄧州市、南陽市、開封市、安陽市結(jié)核病防治所，鶴壁市傳染病醫(yī)院就診的152例MDR肺結(jié)核病患者的痰標(biāo)本，進(jìn)行痰培養(yǎng)。同時(shí)對(duì)患者進(jìn)行問卷調(diào)查，將菌株標(biāo)本和調(diào)查問卷集中送至河南省疾病預(yù)防控制中心進(jìn)行后續(xù)研究。該項(xiàng)目由河南省疾病預(yù)防控制中心醫(yī)學(xué)倫理委員會(huì)批準(zhǔn)。從病歷記錄中獲取信息。就診前接受過不規(guī)則化療超過1個(gè)月的結(jié)核病患者被定義為復(fù)治患者。

1.2pncA、rpsA和panD基因突變檢測(cè)從培養(yǎng)陽性的羅氏培養(yǎng)基上刮取新鮮的培養(yǎng)菌至500 μL TE緩存液中，85 ℃滅活30 min，煮沸5 min后離心待用[5]。取2 μL DNA加入含有引物的PCR預(yù)混液(購自北京康為世紀(jì)生物科技有限公司)，使用粗提的DNA為模板進(jìn)行擴(kuò)增，pncA上游引物5’-AA CAGTTCATCCCGGTTC-3’，下游引物5’-GCGTCAT GGACCCTATATC-3’；rpsA上游引物5’-CGGAG CAACCCAACAATA-3’，下游引物5’-GTGGACAG CAACGACTTC-3’；panD上游引物5’-TCAACGGT TCCGGTCGGCTGCT-3’，下游引物5’-TATCCGC CACTGCTGCACGACCTT-3’。 擴(kuò)增體系20 μL，PCR擴(kuò)增條件為95 ℃預(yù)變性5 min，95 ℃ 30 s，60 ℃ 30 s，72 ℃ 1 min，30個(gè)循環(huán)；最后72 ℃10 min。PCR產(chǎn)物送至生工生物工程(上海)股份有限公司進(jìn)行測(cè)序分析，并與野生型菌株H37Rv序列比較。

1.3Spoligotyping和NTF分析按照文獻(xiàn)[5]進(jìn)行Spoligotyping分析。應(yīng)用DR區(qū)上游引物5’-CCGAGAGGGGACGGAAAC-3’和下游引物5’ -GGTTTTGGGTCTGACGAC-3’對(duì)抽提的基因組進(jìn)行PCR擴(kuò)增，然后將PCR產(chǎn)物與預(yù)先包被43條間隔寡核苷酸的膜雜交，根據(jù)43個(gè)不同間隔區(qū)存在情況判讀結(jié)果。根據(jù)結(jié)核分枝桿菌NTF區(qū)設(shè)計(jì)上游引物5’-CCAGATATCGGGTGTGTCGAC-3’，下游引物5’-TGCCGTTCTCGAAATCTAAACAA-3’[6]；若NTF區(qū)域存在IS6110片段插入為現(xiàn)代型，否則為古典型。

1.4PZA藥敏及最小抑菌濃度(MIC)測(cè)定使用購自美國BD公司的結(jié)核分枝桿菌液體培養(yǎng)基(Bactec MIGT 960)培養(yǎng)細(xì)菌，加入終濃度為100 mg/L的PZA，7～10 d后讀取結(jié)果，報(bào)告陽性的為耐藥菌，陰性的為敏感菌[7]。使用微孔板法進(jìn)行MIC檢測(cè)。在96孔板第1列加入含有PZA濃度400 mg/L的pH5.5的7H9液體培養(yǎng)基，倍比稀釋，直至倒數(shù)第2列PZA濃度為3.13 mg/L。將1個(gè)麥?zhǔn)蠞舛鹊慕Y(jié)核分枝桿菌稀釋100倍后，每孔加入100 μL。37 ℃培養(yǎng)7 d，每孔加入3 μL 體積分?jǐn)?shù)為1%的刃天青和12 μL 體積分?jǐn)?shù)為5% 的Tween-80，繼續(xù)培養(yǎng)24 h，觀察顏色變化。以藍(lán)色孔為無菌生長，粉色孔為有菌生長，藍(lán)色孔對(duì)應(yīng)的最低藥物濃度為MIC。

1.5數(shù)據(jù)分析使用SPSS 19.0進(jìn)行比值比和95%CI的分析來判斷不同組別間PZA耐藥性是否存在差異。根據(jù)標(biāo)準(zhǔn)差計(jì)算公式使用MedCalc在線計(jì)算工具(https://www.MedCalc.net/)計(jì)算靈敏度和特異度的95%CI。

2 結(jié)果

2.1不同特征的MDR結(jié)核患者PZA耐藥分析見表1。152株MDR結(jié)核分枝桿菌菌株中有105株P(guān)ZA耐藥，耐藥率為69.1%。比較不同特征MDR 結(jié)核患者PZA的耐藥性發(fā)現(xiàn)北京家族更易產(chǎn)生PZA耐藥，且現(xiàn)代型菌株對(duì)PZA具有更高的耐藥頻率。

表1 不同特征的MDR結(jié)核患者PZA的耐藥性分析

2.2pncA、rpsA和panD基因突變情況152株MDR結(jié)核分歧桿菌菌株中有102株發(fā)生pncA基因突變，突變率為67.1%(102/152)，其中3個(gè)為同義突變；82株發(fā)生堿基突變(包含3個(gè)同義突變)，20株菌發(fā)生移碼突變，有9種突變?yōu)榘l(fā)生在多個(gè)菌株的共享突變。有133株P(guān)ZA藥敏結(jié)果和pncA測(cè)序結(jié)果一致；有8個(gè)PZA敏感菌株的pncA基因發(fā)生了堿基突變，其中3個(gè)是同義突變；11株P(guān)ZA耐藥菌株的pncA基因未發(fā)生突變。對(duì)發(fā)生了pncA基因突變的8株P(guān)ZA敏感菌菌株進(jìn)行MIC測(cè)定，結(jié)果見表2。

表2 pncA突變的8株P(guān)ZA敏感菌的MIC

rpsA和panD基因突變情況見表3。rpsA和panD的啟動(dòng)子區(qū)均未發(fā)現(xiàn)突變位點(diǎn)，并且大多數(shù)rpsA和panD發(fā)生突變的PZA耐藥菌株均同時(shí)伴隨著pncA的突變。152株MDR結(jié)核分枝桿菌菌株中有76株發(fā)生rpsA基因突變，突變率為72.4%。rpsA的突變散在分布，但主要集中在N末端。在rpsA基因357位和532位核苷酸突變的PZA耐藥菌株中未檢測(cè)到pncA突變；rpsA基因636位核苷酸的同義突變?yōu)榉翘禺愋酝蛔儯?8株P(guān)ZA耐藥菌株和35株P(guān)ZA敏感菌株中均發(fā)生了該突變。在47株P(guān)ZA敏感菌株中未發(fā)現(xiàn)panD基因突變。在105株P(guān)ZA耐藥菌株中有2株菌發(fā)生panD基因突變，其中1株未檢測(cè)到pncA的突變。

表3 rpsA和panD基因突變情況

以PZA的藥敏表型為參照，使用pncA突變檢測(cè)PZA耐藥的靈敏度為89.52%(95%CI為81.64%～94.39%)，特異度為89.36% (95%CI為76.10%～96.01%)，見表4。rpsA/panD突變與PZA的耐藥表型的一致性較差，但rpsA和panD突變與pncA聯(lián)合后，可將僅用pncA檢測(cè)PZA耐藥的靈敏度提高。

表4 MDR結(jié)核分枝桿菌菌株P(guān)ZA耐藥性的效能評(píng)價(jià)

3 討論

PZA在MDR結(jié)核患者治療中發(fā)揮了重要作用[8]。本研究調(diào)查了河南MDR結(jié)核患者的PZA耐藥情況，結(jié)果顯示69.1%的MDR結(jié)核分枝桿菌菌株對(duì)PZA耐藥,與重慶地區(qū)的耐藥率(62.4%)相近[9]，略高于之前的研究結(jié)果[10-13]。

有研究[14]認(rèn)為結(jié)核桿菌的基因譜系與PZA的耐藥性不存在相關(guān)性。本研究結(jié)果表明北京型結(jié)核分枝桿菌的PZA耐藥比例明顯高于非北京型，且現(xiàn)代型高于古典型。但由于樣本量的欠缺，尚不能認(rèn)為PZA的耐藥性和結(jié)核桿菌基因型相關(guān)。

不同地區(qū)PZA耐藥菌株中pncA基因發(fā)生突變的比率有所不同，巴西為45.7%[15]，伊朗為70.6%[16]，臺(tái)灣為75.0%[17]，重慶為88%[9]，瑞典為94.1%[18]。本研究pncA基因突變發(fā)生率為67.1%(102/152)。雖然之前文獻(xiàn)[2,19-20]報(bào)道pncA基因突變與菌株成簇情況有一定關(guān)聯(lián)，由于本研究中納入菌株較少且主要為北京家族基因型，不能判斷這些菌株是否具有成簇性。

與之前的研究[7，12]比較，本研究中的MDR結(jié)核分期桿菌菌株rpsA突變頻率較高，有72.4%(76/105)PZA耐藥的菌株發(fā)生rpsA的突變。636位核苷酸在PZA耐藥和敏感菌株中均有較高比例的同義突變。但在PZA的敏感株中未發(fā)現(xiàn)rpsA的錯(cuò)義或移碼突變。本研究發(fā)現(xiàn)rpsA基因的突變主要集中于C末端，與之前研究[7,10]相一致。此外rpsA的突變位點(diǎn)中有5個(gè)為首次報(bào)道，這些rpsA的突變位點(diǎn)對(duì)PZA耐藥性的影響還需進(jìn)一步驗(yàn)證。本研究結(jié)果表明PZA耐藥菌株中有2株發(fā)生panD基因突變，有1株為pncA野生型耐藥菌株。

本研究PZA的耐藥表型和pncA突變序列相一致的有89.52%。8株菌的pncA基因發(fā)生突變但其藥敏結(jié)果顯示對(duì)PZA敏感(其中3株為同義突變)，有8株pncA野生型分離株卻表現(xiàn)出PZA耐藥。作者對(duì)pncA發(fā)生突變的PZA敏感株進(jìn)行了MIC測(cè)定，除了3株發(fā)生同義突變的菌株，其余5株pncA發(fā)生突變菌株均表現(xiàn)為PZA低水平耐藥(12.5 mg/L

本研究聯(lián)合使用3種基因測(cè)序結(jié)果評(píng)估PZA耐藥性的靈敏度和特異度分別為92.38%和89.36%，與其他研究相一致[23]。

納入的樣本量較少是本研究主要限制，并且所有的菌株均來自同一地區(qū)也降低了結(jié)果的代表性，日后將收集更多樣本進(jìn)行研究。

綜上所述，通過對(duì)河南省MDR結(jié)核患者的pncA、rpsA和panD基因突變情況進(jìn)行檢測(cè)和分析，可快速、準(zhǔn)確地預(yù)測(cè)PZA耐藥性。