耐药结核分枝杆菌基因组信息挖掘

<正>由于结核分枝杆菌的耐药性问题,使结核病这个古老的传染病死灰复燃,并成为全球性的严重公关问题。从1998年首个结核分枝杆菌(H37Rv)全基因组完成图的获得[1],到近年来采用新一代测序技术对多个菌株进行高通量的基因组测序与分析,对结核分枝杆菌进化及耐药机制的认识不断深入。关于结核分枝杆菌菌株的基因组比较分析有多篇报道,包括对中国12个省来源的161株结核分枝杆菌     

定量蛋白质组学分析链霉素耐药和敏感结核分枝杆菌临床分离株

[目的]发现结核分枝杆菌(Mycobacterium tuberculosis)链霉素耐药相关的潜在菌体蛋白.[方法]以结核分枝杆菌临床分离链霉素敏感株01105和结核分枝杆菌H37Rv为对照,采用iTRAQ技术和生物信息学鉴定并相对定量结核分枝杆菌临床分离链霉素耐药株01108菌体蛋白,并通过WEGO功能注释聚类分析01108菌株差异表达蛋白的细胞组分、分子功能和生物进程.[结果]01108菌株分别与01105菌株和H37Rv菌株比较差异表达蛋白为194个和146个,01108菌株与01105菌株和H37Rv比较均差异表达蛋白121个(共同差异表达蛋白).差异表达蛋白理论相对分子量和等电点分布广泛,其生物进程主要参与中间代谢、呼吸作用和脂质代谢,分子功能主要为催化活性功能和结合功能.共同差异表达蛋白:7个核糖体蛋白(Rv2785c,Rv0056,Rv0641,Rv0652,Rv0701,Rv1630和Rv2442c)在01108菌株中表达下调;7个蛋白在01108菌株中显著差异表达(上调大于1.20倍或下调小于0.55倍),分别为巯基过氧化物酶(Rv1932)、酰基载体蛋白脱氢酶(Rv0824c)、30S核糖体蛋白S15 (Rv2785c)、丙酮酸脱氢酶E2部分(Rv2215)、双组份转录调控蛋白(Rv3133c)以及假定未知蛋白(Rv2466c和Rv2626c).[结论]iTRAQ发现了链霉素耐药结核分枝杆菌相对于链霉素敏感结核分枝杆菌和H37Rv共同差异表达蛋白,为进一步探讨结核分枝杆菌链霉素耐药机制奠定了基础.     

结核分枝杆菌katG基因突变的研究

探讨编码过氧化氢-过氧化物酶的katG基因突变与结核分枝杆菌异烟肼(INH)耐药性的相关关系。根据结核分枝杆菌GenBank中的katG序列,自行设计特异性寡聚核苷酸引物,采用聚合酶链反应-单链构象多态性(PCR-SSCP)分析和直接测序法(DS)分析结核分枝杆菌中katG基因突变情况。以HR37Rv标准株为对照。所有23株敏感菌均未有SSCP结果异常;35株耐药菌中,有2株(5．7％)katG基因扩增阴性,且发生在高度耐药菌中。进一步分析发现,SSCP法突变检出23株(65．7％),测序法突变检出24株(68．6％),符合率为95．8％(23／24)。参照测序法对耐药菌突变序列的分析结果,PCR—SSCP敏感、特异,可快速检测结核分枝杆菌katG耐药基因突变,有利于耐药结核分枝杆菌耐药性的快速检测。     

结核分枝杆菌H37有毒株与无毒株的代谢相关基因pabB和lpdA启动子活性分析

【目的】通过分析结核分枝杆菌无毒株H37Ra的全基因组序列,并与H37Rv基因组序列比较,发现pabB和lpdA预测的启动子区发生了突变。我们利用报告基因,确认启动子突变与其基因转录水平的关系,探索结核分枝杆菌H37Ra毒力丧失的内在原因。【方法】利用生物信息学方法预测这两对基因的启动子区,采用PCR技术克隆这两对基因的启动子,与分枝杆菌启动子探针载体pMC210相连,DNA测序证实连接片段正确后,用电穿孔法将重组质粒转化至耻垢分枝杆菌mc2155。利用Quantitative Real-Time RT-PCR检测报告基因lacZ转录水平的差异,进一步验证这两对基因启动子的突变对相应基因转录水平的影响。【结果】Quantitative Real Time PCR检测结果显示H37RapabB启动子活性是H37Rv pabB启动子活性的6倍(p0.05),而H37Rv lpdA启动子的活性是H37RalpdA启动子的2倍(p0.05)。【结论】pabB,lpdA的启动子在H37Ra中的突变对其启动子的活性产生了影响,其中lpdA启动子的突变可能与结核分枝杆菌H37Ra的毒力丧失有关。     

新疆南疆MTB临床分离株异烟肼耐药情况及耐药相关基因的突变研究

目的:探索新疆南疆部分地区维吾尔族结核病人群中分离到的结核分枝杆菌基因组中的katG、inhA、kasA、ahpC基因突变与耐异烟肼的关联,揭示新疆肺结核病高患病率、高死亡率的可能原因.方法:收集新疆南疆部分地区维吾尔族结核病患者痰液标本,对MTB进行分离培养后,应用比例法检测其对异烟肼耐药性,运用PCR技术对所分离菌株的上述基因相关片段扩增并进行测序及序列分析.结果:筛选出实验菌株99株,其中敏感株63株,耐异烟肼菌株36株,耐药率为36.4%.耐药株中katG基因突变率为63.9%,315位点突变,突变类型AGC→ACC(Ser→Thr)、AGC→AAC(Ser→Asn);inhA基因突变率为47.2%,有缺失、同义突变和错义突变;kasA基因突变率为41.7%,突变类型为缺失和错义突变;ahpC基因突变率为8.3%,属于错义突变.结论:新疆南疆结核病高患病率、高死亡率的可能原因之一是结核分枝杆菌对INH产生了耐药,结核分枝杆菌对INH耐药机制可能与耐药菌株基因组内katG、inhA、kasA和ahpC基因发生了突变有关.     

结核分枝杆菌中毒素-抗毒素系统的鉴定

【目的】鉴定结核分枝杆菌基因组上MazF同源蛋白基因与其上游基因是否组成毒素-抗毒素系统,阐明毒素蛋白的作用机理,并初步探讨毒素-抗毒素系统在营养缺乏时的表达调控。【方法】在大肠杆菌和耻垢分枝杆菌中将MazF同源蛋白单独表达或与其对应的抗毒素蛋白共同表达,鉴定MazF同源蛋白对细菌生长的抑制作用以及其对应的抗毒素蛋白能否消除这种生长抑制;通过体外RNA切割实验,检测MazF同源蛋白是否具有RNA切割活性;检测正常生长条件下和饥饿条件下毒素-抗毒素系统的启动子活性,探讨其在应激条件下的表达调控。【结果】结核分枝杆菌MazF同源蛋白中,Rv0659c、Rv1495和Rv1942c不具有抑制细菌生长的毒素蛋白活性,Rv1991c、Rv2801c、Rv1102c和mtPemK能够抑制细菌生长,而且它们的抑制作用可以分别被其对应的抗毒素Rv1991a、Rv2801a、Rv1103c和mtPemI解除。Rv1991c、Rv2801c和Rv1102c具有RNA切割活性,mtPemK则不能切割RNA。Rv1991a-1991c和Rv2801a-2801c系统的启动子在饥饿条件下活性显著升高。【结论】结核分枝杆菌基因组上Rv1991a-1991c、Rv2801a-2801c、Rv1103c-1102c和mtPemI-mtPemK是毒素-抗毒素系统。毒素蛋白Rv1991c、Rv2801c和Rv1102c通过切割RNA发挥抑菌或杀菌活性,mtPemK具体作用机理目前还不清楚。Rv1991a-1991c和Rv2801a-2801c系统可能参与结核分枝杆菌在营养匮乏条件下的生长调控。     

结核分枝杆菌毒素-抗毒素系统 maz EF6缺失突变株的构建及其表型的初步探讨

为构建结核分枝杆菌毒素‐抗毒素系统 m azEF6缺失突变株,并对其表型进行初步探讨,首先用聚合酶链反应（PCR）分别从H37Rv标准株和PUC‐19K质粒扩增出 mazEF6基因的同源臂及卡那霉素抗性基因kan ;然后应用融合PCR技术将 mazEF6基因的同源臂与 kan基因进行杂交拼接,获得目的融合片段,将该融合片段克隆于pMD‐19T（simple）载体形成自杀质粒pMD‐19T‐ΔmazEF6‐kan ,并将自杀质粒转化至大肠埃希菌DH5α中;最后利用电穿孔技术将自杀质粒电转至H37Rv标准株中,在卡那霉素抗性改良罗氏培养基上筛选H37Rv ΔmazEF6缺失突变株单个菌落,提取阳性菌株全基因组DNA为模板,PCR扩增克隆片段并测序。将所获得的H37Rv ΔmazEF6缺失突变株进行遗传稳定性检测后,对其表型进行初步研究。结果显示,该缺失株在15代内未发生回复性突变;与野生株相比,缺失株生长速度缓慢且细菌形态短小。本研究证实,融合PCR技术便于快速获得结核分枝杆菌缺失突变株;结核分枝杆菌在缺失毒素‐抗毒素系统 m azEF6基因后生存能力下降,这为进一步研究毒素‐抗毒素系统的作用奠定了基础。     

相对定量分析异烟肼、链霉素单耐药与多耐药结核分枝杆菌临床分离株的蛋白质组差异

【目的】探讨异烟肼(isoniazid,INH)、链霉素(streptomycin,SM)单耐药结核分枝杆菌(Mycobacterium tuberculosis,MTB)与INH/SM多耐药MTB蛋白质组差异。【方法】应用i TRAQ结合Nano LC-MS/MS定量蛋白质组学技术,分析临床分离INH、SM或INH/SM耐药MTB与H37Rv标准株间均表达差异蛋白;并以INH/SM耐药MTB与H37Rv比值为对照,相对定量分析单耐药与多耐药MTB蛋白表达差异倍数;运用DAVID 6.7分析差异蛋白生物功能;STITCH 5.0分析差异蛋白与INH和SM相互作用。【结果】与H37Rv标准株比较,58个蛋白在INH、SM耐药与INH/SM耐药MTB间均有表达差异,共同差异蛋白生物功能主要为氧化还原酶活性和转移酶活性;主要参与丙酸代谢信号通路。共同差异蛋白中,与INH/SM耐药MTB比较,Rv2986c和Rv1908c在INH、SM耐药MTB均表达上调1.25倍;Rv3133c和Rv0577则均表达下调0.7倍;生物信息学预测发现以上4种蛋白可直接或间接与INH、SM进行相互作用。【结论】INH、SM单耐药和INH/SM多耐药MTB蛋白表达谱有较大差异,蛋白Rv2986c、Rv1908c、Rv3133c和Rv0577表达水平及相互作用可能与INH和SM耐药有关。     

通辽地区犊牛腹泻大肠杆菌耐药性检测及一株多重耐药菌全基因组测序分析

【背景】大肠杆菌(Escherichia coli)是引起犊牛腹泻的最主要病原菌,其耐药性菌株的不断出现引起广泛关注。【目的】了解内蒙古自治区通辽市犊牛腹泻大肠杆菌耐药性及耐药基因流行情况。【方法】从通辽市多个旗县采集犊牛腹泻样品40份,经细菌分离纯化及16S rRNA基因测序,最终鉴定出20株大肠杆菌。采用药敏试验和PCR方法对分离菌进行耐药性及耐药基因检测分析,并对其中1株多重耐药菌株进行全基因组测序。【结果】20株分离菌均具有多重耐药性,对链霉素、环丙沙星、恩诺沙星和复方新诺明的耐药率达80%以上。所检耐药基因中,aphA1、strB、TEM-1和qnrS检出率达100%。通过对代表性菌株TL-13全基因组测序发现,其基因组大小为4897185bp,GC含量为50.68%,同时携带2个质粒,大小分别为108288bp(pTL13-1)和64018bp(pTL13-2)。质粒中共携带18个可移动耐药基因。【结论】通辽地区犊牛腹泻大肠杆菌多重耐药性普遍存在,4种常见耐药基因普遍流行。     

雄烯二酮高产菌新金分枝杆菌MN4的全基因组测序及序列分析北大核心CSCD

【目的】新金分枝杆菌(Mycobacterium neoaurum)MN4是1株经诱变育种获得的高产雄烯二酮,并且能够耐受高浓度底物植物甾醇的突变菌株。为深入研究MN4菌株耐底物的机制及雄烯二酮的生物合成途径,有必要解析MN4菌株的全基因组序列信息。【方法】本研究采用高通量测序技术对MN4进行全基因组测序,然后使用相关软件对测序数据进行基因组装、基因预测与功能注释、COG聚类分析以及次级代谢产物合成基因簇预测等。【结果】新金分枝杆菌MN4基因组装获得33个Contigs,整个基因组大小为5.39 Mb,GC含量为66.9%,编码4920个蛋白基因,序列提交至Gen Bank数据库,登录号为JXYZ00000000。【结论】本研究首次报道了1株高产雄烯二酮菌株MN4的全基因组序列,分析了基因组的基本特征,初步解析了该菌株降解植物甾醇生产雄烯二酮的关键基因,将为MN4的功能基因组学研究及相关次级代谢产物的生物合成途径与异源表达研究提供基础。     

Efflux Pump Gene Expression in Multidrug-Resistant Mycobacterium tuberculosis Clinical Isolates

Isoniazid (INH) and rifampicin (RIF) are the two most effective drugs in tuberculosis therapy. Understanding the molecular mechanisms of resistance to these two drugs is essential to quickly diagnose multidrug-resistant (MDR) tuberculosis and extensive drug-resistant tuberculosis. Nine clinical Mycobacterium tuberculosis isolates resistant to only INH and RIF and 10 clinical pan-sensitive isolates were included to evaluate the expression of 20 putative drug efflux pump genes and sequence mutations in rpoB (RIF), katG (INH), the inhA promoter (INH), and oxyR-ahpC (INH). Nine and three MDR isolates were induced to overexpress efflux pump genes by INH and RIF, respectively. Eight and two efflux pump genes were induced to overexpress by INH and RIF in MDR isolates, respectively. drrA, drrB, efpA, jefA (Rv2459), mmr, Rv0849, Rv1634, and Rv1250 were overexpressed under INH or RIF stress. Most efflux pump genes were overexpressed under INH stress in a MDR isolates that carried the wild-type katG, inhA, and oxyR-ahpC associated with INH resistance than in those that carried mutations. The expression levels of 11 genes (efpA, Rv0849, Rv1250, P55 (Rv1410c), Rv1634, Rv2994, stp, Rv2459, pstB, drrA, and drrB) without drug inducement were significantly higher (P < 0.05) in nine MDR isolates than in 10 pan-sensitive isolates. In conclusion, efflux pumps may play an important role in INH acquired resistance in MDR M. tuberculosis, especially in those strains having no mutations in genes associated with INH resistance; basal expression levels of some efflux pump genes are higher in MDR isolates than in pan-sensitive isolates and the basal expressional differences may be helpful to diagnose and treat resistant tuberculosis.     

Efflux pump inhibition by 11H-pyrido[2,1-b]quinazolin-11-one analogues in mycobacteria

Mycobacterium tuberculosis infection causes 1.8 million deaths worldwide, of which half a million has been diagnosed with resistant tuberculosis (TB). Emergence of multi drug resistant and extensive drug resistant strains has made all the existing anti-TB therapy futile. The major involvement of efflux pump in drug resistance has made it a direct approach for therapeutic exploration against resistant M. tuberculosis. This study demarcates the role of 11H-pyrido[2,1-b]quinazolin-11-one (quinazolinone) analogues as efflux pump inhibitor in Mycobacterium smegmatis. Sixteen quinazolinone analogues were synthesized by treating 2-aminopyridine and 2-fluorobenzonitrile with K^tOBu. Analogues were tested, and 3a, 3b, 3c, 3g, 3j, 3l, 3m, and 3p were found to modulate EtBr MIC by >4 whereas 3a, 3g, 3i and 3o showed >4 modulation on norfloxacin MIC. 3l and 3o in addition to their very low toxicity they showed high EtBr and norfloxacin accumulation respectively. Time kill curve showed effective log reduction in colony forming unit in presence of these analogues, thus confirming their role as efflux pump inhibitor. Through docking and alignment studies, we have also shown that the LfrA amino acid residues that the analogues are interacting with are present in Rv2333c and Rv2846c of M. tuberculosis. This study have shown for the first time the possibility of developing the 11H-pyrido[2,1-b]quinazolin-11-one analogues as efflux pump inhibitors for M. smegmatis and hence unbolts the scope to advance this study against resistant M. tuberculosis as well.     

Contribution of efflux to the emergence of isoniazid and multidrug resistance in Mycobacterium tuberculosis

Multidrug resistant (MDR) tuberculosis is caused by Mycobacterium tuberculosis resistant to isoniazid and rifampicin, the two most effective drugs used in tuberculosis therapy. Here, we investigated the mechanism by which resistance towards isoniazid develops and how overexpression of efflux pumps favors accumulation of mutations in isoniazid targets, thus establishing a MDR phenotype. The study was based on the in vitro induction of an isoniazid resistant phenotype by prolonged serial exposure of M. tuberculosis strains to the critical concentration of isoniazid employed for determination of drug susceptibility testing in clinical isolates. Results show that susceptible and rifampicin monoresistant strains exposed to this concentration become resistant to isoniazid after three weeks; and that resistance observed for the majority of these strains could be reduced by means of efflux pumps inhibitors. RT-qPCR assessment of efflux pump genes expression showed overexpression of all tested genes. Enhanced real-time efflux of ethidium bromide, a common efflux pump substrate, was also observed, showing a clear relation between overexpression of the genes and increased efflux pump function. Further exposure to isoniazid resulted in the selection and stabilization of spontaneous mutations and deletions in the katG gene along with sustained increased efflux activity. Together, results demonstrate the relevance of efflux pumps as one of the factors of isoniazid resistance in M. tuberculosis. These results support the hypothesis that activity of efflux pumps allows the maintenance of an isoniazid resistant population in a sub-optimally treated patient from which isoniazid genetically resistant mutants emerge. Therefore, the use of inhibitors of efflux should be considered in the development of new therapeutic strategies for preventing the emergence of MDR-TB during treatment.     

QTL‐seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L.)

Rust and late leaf spot (LLS) are the two major foliar fungal diseases in groundnut, and their co‐occurrence leads to significant yield loss in addition to the deterioration of fodder quality. To identify candidate genomic regions controlling resistance to rust and LLS, whole‐genome resequencing (WGRS)‐based approach referred as ‘QTL‐seq’ was deployed. A total of 231.67 Gb raw and 192.10 Gb of clean sequence data were generated through WGRS of resistant parent and the resistant and susceptible bulks for rust and LLS. Sequence analysis of bulks for rust and LLS with reference‐guided resistant parent assembly identified 3136 single‐nucleotide polymorphisms (SNPs) for rust and 66 SNPs for LLS with the read depth of ≥7 in the identified genomic region on pseudomolecule A03. Detailed analysis identified 30 nonsynonymous SNPs affecting 25 candidate genes for rust resistance, while 14 intronic and three synonymous SNPs affecting nine candidate genes for LLS resistance. Subsequently, allele‐specific diagnostic markers were identified for three SNPs for rust resistance and one SNP for LLS resistance. Genotyping of one RIL population (TAG 24 × GPBD 4) with these four diagnostic markers revealed higher phenotypic variation for these two diseases. These results suggest usefulness of QTL‐seq approach in precise and rapid identification of candidate genomic regions and development of diagnostic markers for breeding applications.     

The role of ABC efflux pump,Rv1456c-Rv1457c-Rv1458c,from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> clinical isolates in China

Recently the ATP-binding cassette (ABC) efflux pumps have been proved to be a major component of drug resistance in Mycobacterium tuberculosis. The objective of this study was to investigate the expression profiles of Rv1456c-Rv1457c-Rv1458c efflux system in clinical isolates of M. tuberculosis and its involvement in drug-resistance mechanisms. Significantly increased mRNA expression of Rv1456c, Rv1457c, and Rv1458c appeared among the clinical isolates (P < 0.05), which are resistant to at least one of the four first-line drugs including rifampin, isoniazid, streptomycin, and ethambutol. In addition, overexpression of this efflux system was more frequently found in multidrug-resistant and extensively drug-resistant M. tuberculosis strains. Therefore, Rv1456c-Rv1457c-Rv1458c efflux pumps may play an important role in drug resistance of treatment of M. tuberculosis. Further investigation of this gene may lead to the development of countermeasures against M. tuberculosis drug resistance.     

Molecular cloning, overexpression and biochemical characterization of hypothetical beta-lactamases of Mycobacterium tuberculosis H37Rv

Aim:  Molecular cloning, overexpression and biochemical characterization of the genes from the Mycobacterium tuberculosis H37Rv genome having hypothetical β-lactamases activity.
Methods and Results:  Analysis of the M. tuberculosis H37Rv genome revealed that Rv 2068c , Rv 0406 c and Rv 3677 c gene products were predicted to exhibit β-lactamases activity. All the three genes were cloned in pET28a vector and overexpressed in C41 (DE3) Escherichia coli cells. The His-tagged recombinant proteins were confirmed by immunoblotting and were shown to have β-lactamase activity by the hydrolysis of nitrocefin and other β-lactams. Catalytic parameters for all the recombinant proteins were derived followed by the enzyme inhibition studies. Antibiotic susceptibility studies using the recombinant strains showed an increased resistance against different classes of β-lactam antibiotics.
Conclusion:  The study revealed the possibility of more than one gene in M. tuberculosis , encoding proteins having β-lactamase or β-lactamase-like activity, giving wide spectrum of resistance against β-lactams.
Significance and Impact of the Study:  Systematic study of hypothetical β-lactamases of M. tuberculosis and related species and their correlation with β-lactam and inhibitor susceptibility profile might be useful in developing new antibiotic regime for the treatment of tuberculosis caused by multiple drug resistant (MDR) strains.