利用抑制消减杂交技术研究结核分支杆菌强毒株H37Rv和弱毒株H37Ra的基因差异

为了解结核病的致病分子机理和筛选结核病致病菌的毒力基因,利用抑制消减杂交（SSH）技术分析了结核分支杆菌强毒株H37Rv和弱毒株H37Ra间的基因组DNA间差异。通过Southern杂交验证及序列分析得到仅在强毒株H37Rv基因组中有的DNA片段8个,其中一个编码已知的毒力因子mce蛋白,1个编码PE家族蛋白,1个编码purC合成酶,和4个潜在蛋白,另一个为非编码区片段。其中有2个基因经PCR方法已证实在强毒株H37Rv和临床分离的强毒株中存在,而在H37Ra和临床弱毒株中无;仅在弱毒株H37Ra基因组中的DNA片段3个,其中2个为新基因片段,已被GenBank收录。     

结核分枝杆菌拟核结合蛋白Lsr2阻遏调控rv1057基因转录的研究

探讨结核分枝杆菌拟核结合蛋白Lsr2是否是rv1057基因转录的阻遏蛋白。利用凝胶阻滞迁移试验分析Lsr2在体外与rv1057启动子片段的特异性结合能力;通过lacZ报告基因和荧光定量PCR检测rv1057启动子区域的不同突变形式对转录的影响;启动转录的能力;通过蛋白质免疫印迹法检测目标基因rv1057的表达。统计学处理采用t检验。Lsr2结合rv1057启动子;Lsr2与已知的阻遏蛋白TrcR都能结合rv1057启动子;缺失Lsr2结合位点的rv1057启动子能够在结核分枝杆菌生长周期中持续启动rv1057基因的转录;野生型菌株H37Rv生长早期(24 h)lsr2基因转录水平较低,而在H37Rv生长中期(48 h)、后期(120 h)lsr2基因的相对表达量变化倍数显著提高,差异具有统计学意义(t值分别为24.44、16.86,P0.05);H37Rv菌株和lsr2基因缺失菌株中trcR基因在生长早期(24 h)的相对表达量显著高于生长中期(48 h)、生长后期(120 h)的trcR表达量,差异具有统计学意义(t值分别为6.79、10.16,P0.05)。Lsr2是rv1057基因转录的阻遏蛋白,Lsr2和TrcR均可调控rv1057基因的表达。Lsr2在结核分枝杆菌生长周期的中后期大量表达并阻遏rv1057的转录,TrcR则在结核分枝杆菌生长周期的早期阻遏rv1057的转录。     

利用表型芯片系统高通量分析H37Ra与H37Rv差异表型

H37Rv是结核分枝杆菌标准有毒株,H37Ra是从H37Rv获得的稳定减毒株,但目前H37Ra毒力减弱原因尚不完全清楚。本研究利用表型芯片系统,高通量分析H37Ra生长表型,并与H37Rv表型比较,筛选两菌株表型差异,分析与H37Ra毒力减弱可能的相关表型及分子机制。结果发现,与H37Rv相比,H37Ra耐酸及耐渗透压能力显著下降,且不能利用丁二酸单甲酯和吐温40作为碳源。结核分枝杆菌耐酸能力直接影响其在吞噬体中的生存和代谢,耐高渗能力影响其必需营养物质的跨膜运输,代谢途径的改变影响其在宿主内的能量摄取,三者改变均可能与H37Ra毒力减弱相关。     

两株超级广泛耐药结核分枝杆菌全基因组比较分析

【目的】结合现有数据,通过对两株临床超级广泛耐药的结核分枝杆菌全基因组的测序和分析,发现其型别相关的突变位点,解释发生广泛耐药的基因组突变机制。【方法】利用Solexa第二代测序技术对两株广泛耐药结核分枝杆菌(FJ05194和GuangZ0019)进行全基因组测序分析。以H37Rv为参考序列得到两株广泛耐药菌株的单核苷酸多态性(SNPs),构建系统发育树鉴定菌株型别,判断突变位点中型别相关和非型别相关的SNPs。定位SNPs所在的基因组区域,对型别相关的突变基因进行KEGG通路的富集分析,对非型别相关的突变基因和间隔区判断是否与耐药相关。【结果】两株广泛耐药菌株分别属于Lineage2和Lineage4型别,两菌株在碱基替换方面存在差异性,Lineage2型别相关的基因功能富集于ABC转运蛋白和核苷酸切除修复的通路。耐药方面,发现了已知的耐药相关基因的突变(rpoB、katG、rpsl、gyrA、gyrB、embB和ethA等),但卷曲霉素和卡那霉素相关的rrs、tlyA和eis启动子区域未发生突变,不足以解释其耐药性的产生。与最新报道的候选耐药基因比较,发现了卷曲霉素和卡那霉素相关的突变(Rv1393c、Rv0265c和narX等)和外排泵相关的pstB、Rv2333c和Rv2687c突变。【结论】结核分枝杆菌Lineage2型别相关的SNPs中含有影响结核分枝杆菌突变率和耐药性的突变。对于两株超级广泛耐药的结核菌,已知的激活药物或药靶相关的单耐药基因突变集合不能完全解释其广泛耐药性,还涉及新候选结核耐药基因、外排泵和补偿等其他潜在机制的相关基因突变。     

大肠杆菌yigP基因转录调控序列的鉴定

【目的】调查yigP基因启动子的活性,并对该转录调控序列进行分析。【方法】以lacZ为报告基因,克隆启动子片段至启动子探针质粒中,通过检测β-半乳糖苷酶活性判断启动子活性,并通过克隆一系列逐步缩短的启动子片段来确定启动子所在区域。利用定点突变技术,对启动子的重要序列进行定点突变,调查其对启动子活性的影响。【结果】确定了yigP基因启动子的区域,鉴定了启动子的-10区和-35区,并发现了启动子上游存在一个负调控序列,对该序列进行了初步的研究显示其中部分序列是这种负调控作用的核心序列。【结论】对yigP基因的转录调控序列进行了鉴定,丰富了我们对基因转录调控的认识。     

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

【目的】鉴定结核分枝杆菌基因组上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系统可能参与结核分枝杆菌在营养匮乏条件下的生长调控。     

结核分枝杆菌Rv1009结构域基因的克隆、表达及亲和层析纯化

目的:克隆结核分枝杆菌Rv1009结构域基因,经序列测定正确后进行融合表达和纯化。方法:采用PCR从结核分枝杆菌H37Rv基因组中扩增出Rv1009结构域基因,用限制性内切酶消化后插入pUC-19克隆载体中,经测序正确后亚克隆到融合表达载体pPro-EXHT中,转化大肠杆菌DH5α,目的基因经IPTG诱导,由T7启动子调控表达了N端带6个连续组氨酸残基的Rv1009结构域多肽,在变性条件下对目的蛋白进行纯化。结果:获得了结核分枝杆菌Rv1009结构域基因,得到融合6个组氨酸残基的Rv1009结构域多肽,纯化获得的蛋白纯度大于87%。结论:构建了结核分枝杆菌Rv1009结构域基因的重组表达载体,并获得了高纯度的融合表达蛋白,为后续深入研究奠定了基础。     

依赖PrfA转录调控的单核细胞增生李斯特菌毒力基因inlC启动子结构特点的初步研究

为了研究单核细胞增生李斯特菌毒力基因启动子的结构特点与转录调控因子PrfA蛋白之间的关系,应用PCR定点突变和重组PCR技术缺失了该菌毒力基因inlC启动子上可能与PrfA蛋白结合以及诱发转录起始相关的碱基序列,构建了一系列突变启动子与lacZ报告基因融合表达质粒,使lacZ基因的表达置于inlC突变启动子下,并分别电转化单核细胞增生李斯特菌野生株P14、PrfA蛋白高表达突变株P14a和prfA基因等位缺失突变株A42中,检测相应的β-半乳糖苷酶活性。结果表明:位于inlC启动子转录起始点下游22bp处的一段17bp的类似PrfA蛋白结合序列TTAACAGCGTTTGTTAA并没有增强和抑制PrfA转录调控活性的功能;甚至将其改造成“完美的”PrfA蛋白结合序列TTAACATTTGTTAA后,也不影响inlC依赖于PrfA的转录活性地表达;但是,如果缺失inlC启动子上原始的PrfA蛋白结合序列,则使inlC依赖于PrfA的转录活性完全丧失;另外,单核细胞增生李斯特菌毒力基因inlC和plcA依赖于PrfA的转录活性的表达也与启动子上PrfA蛋白结合区(PrfA-box)距离-10区的碱基个数有关:最适为22或23bp,长于23bp或短于22bp的突变启动子的依赖PrfA的转录活性大大降低,甚至没有活性。说明除PrfA蛋白结合序列外,受PrfA调控的毒力基因启动子上还可能存在其它尚未阐明的结构和序列影响PrfA蛋白的结合以及启动转录表达。     

依赖PrfA转录调控的单核细胞增生李斯特菌毒力基因inlC启动子结构特点的初步研究

为了研究单核细胞增生李斯特菌毒力基因启动子的结构特点与转录调控因子PrfA蛋白之间的关系,应用PCR定点突变和重组PCR技术缺失了该菌毒力基因inlC启动子上可能与PrfA蛋白结合以及诱发转录起始相关的碱基序列,构建了一系列突变启动子与lacZ报告基因融合表达质粒, 使lacZ基因的表达置于inlC突变启动子下,并分别电转化单核细胞增生李斯特菌野生株P14、PrfA蛋白高表达突变株P14a 和prfA基因等位缺失突变株A42中,检测相应的β-半乳糖苷酶活性。结果表明：位于inlC启动子转录起始点下游22bp 处的一段17bp的类似PrfA蛋白结合序列TTAACAGCGTTTGTTAA并没有增强和抑制PrfA转录调控活性的功能;甚至将其改造成“完美的” PrfA蛋白结合序列TTAACATTTGTTAA后,也不影响inlC依赖于PrfA的转录活性地表达;但是,如果缺失inlC启动子上原始的PrfA蛋白结合序列,则使inlC依赖于PrfA的转录活性完全丧失;另外,单核细胞增生李斯特菌毒力基因inlC和plcA 依赖于PrfA的转录活性的表达也与启动子上PrfA蛋白结合区(PrfA-box)距离-10区的碱基个数有关:最适为22或23bp,长于23bp或短于22bp的突变启动子的依赖PrfA的转录活性大大降低,甚至没有活性。说明除PrfA蛋白结合序列外,受PrfA调控的毒力基因启动子上还可能存在其它尚未阐明的结构和序列影响PrfA蛋白的结合以及启动转录表达。     

鸟分枝杆菌PhoP功能分析及其基因突变株的构建

【目的】对鸟分枝杆菌PhoP的功能进行分析及构建PhoP基因突变株,为深入研究PhoP的调控机制打下基础。【方法】利用PCR扩增出鸟分枝杆菌PhoP DNA结合区(PhoPC)编码序列,与表达载体p GEX-4T-3连接后,转化入大肠杆菌BL21(DE3)中表达GST-PhoPC融合蛋白。用凝血酶去除GST标签,制备PhoPC蛋白;利用PCR扩增出鸟分枝杆菌PhoP基因及其下游基因MAV0127、PhoU和Amt的启动子片段,采用凝胶迁移率移动试验(EMSA)分别检测PhoPC与PhoP、MAV0127、PhoU和Amt的启动子结合的情况。通过PCR扩增PhoP基因上、下游片段,构建PhoP基因缺失性同源核苷酸片段,与自杀质粒p GMB151连接后,通过电转化导入鸟分枝杆菌进行同源交换,利用PCR筛选出PhoP基因缺失突变株。【结果】EMSA结果显示,鸟分枝杆菌PhoP能与PhoP、MAV0127及Amt基因启动子结合,不能与PhoU结合。通过PCR和序列分析证实基因突变株的PhoP基因缺失了309个碱基。【结论】PhoP不仅可调控其下游基因MAV0127和Amt的转录水平,还可调控其自身基因的转录,但不参与调节PhoU二元调控系统。构建了PhoP基因缺失突变株,为进一步研究其在鸟分枝杆菌的调控功能奠定了基础。     

差显技术分析结核杆菌H37Rv与H37Ra差异表达的基因

采用差异显示技术比较了结核分支杆菌强毒株H37Rv和弱毒株H37Ra在体外培养条件下的基因表达差异。通过20种引物组合进行mRNA差异显示,克隆到了两菌株间的20余个差异表达基因,经序列分析及杂交鉴定发现其中2个基因仅在H37Rv中表达。其中Rv1894c基因编码的可能是H37Rv中的一个新蛋白。而在H37Ra的基因组中含有这2个基因的编码序列,但均未检测到基因的表达。     

Comparative proteome analysis of culture supernatant proteins of Mycobacterium tuberculosis H37Rv and H37Ra

To examine the virulence factors of Mycobacterium tuberculosis H37Rv, the proteome was used to characterize the differences in protein expression between virulent M. tuberculosis H37Rv and attenuated M. tuberculosis H37Ra. Two-dimensional gel electrophoresis was performed to separate culture supernatant proteins extracted from M. tuberculosis H37Rv and M. tuberculosis H37Ra. The protein spots of interest were identified by mass spectrometry, and then the genes encoding the identified proteins were cloned and sequenced. Comparison of silver-stained gels showed that three well-resolved protein spots were present in M. tuberculosis H37Rv but absent from M. tuberculosis H37Ra. Protein spot no. 1 was identified as Rv2346c. Protein spot no. 2 was identified as Rv2347c, Rv1197, Rv1038c, and Rv3620c, which shared significant homology and had the same peptide fingerprinting using tryptic digestion. No M. tuberculosis protein matched protein spot no. 3. Rv2346c, Rv2347c, Rv1038c, and Rv3620c of M. tuberculosis H37Rv were located on the M. tuberculosis H37Ra chromosome, and multiple mutations were observed in the corresponding areas of M. tuberculosis H37Ra. Codon 59 (CAG, Gln) of Rv2347c and Rv3620c was replaced by termination codon (TAG) in M. tuberculosis H37Ra, which probably terminated the polypeptide elongation. These results demonstrate the importance of studying the gene products of M. tuberculosis and show that subtle differences in isogenic mutant strains might play an important role in identifying the attenuating mutations.     

A point mutation in the two-component regulator PhoP-PhoR accounts for the absence of polyketide-derived acyltrehaloses but not that of phthiocerol dimycocerosates in Mycobacterium tuberculosis H37Ra

Similarities between Mycobacterium tuberculosis phoP-phoR mutants and the attenuated laboratory strain M. tuberculosis H37Ra in terms of morphological and cytochemical properties, lipid content, gene expression and virulence attenuation prompted us to analyze the functionality of this two-component regulator in the latter strain. Sequence analysis revealed a base substitution resulting in a one-amino-acid change in the likely DNA-binding region of PhoP in H37Ra relative to H37Rv. Using gel-shift assays, we show that this mutation abrogates the ability of the H37Ra PhoP protein to bind to a 40-bp segment of its own promoter. Consistent with this result, the phoP gene from H37Rv but not that from H37Ra was able to restore the synthesis of sulfolipids, diacyltrehaloses and polyacyltrehaloses in an isogenic phoP-phoR knock-out mutant of M. tuberculosis Moreover, complementation of H37Ra with phoP from H37Rv fully restored sulfolipid, diacyltrehalose and polyacyltrehalose synthesis, clearly indicating that the lack of production of these lipids in H37Ra is solely due to the point mutation in phoP. Using a pks2-3/4 knock-out mutant of M. tuberculosis H37Rv, evidence is further provided that the above-mentioned polyketide-derived acyltrehaloses do not significantly contribute to the virulence of the tubercle bacillus in a mouse model of infection. Reasons for the attenuation of H37Ra thus most likely stand in other virulence factors, many of which are expected to belong to the PhoP regulon and another of which, unrelated to PhoP, appears to be the lack of production of phthiocerol dimycocerosates in this strain.     

Comparison of membrane proteins of Mycobacterium tuberculosis H37Rv and H37Ra strains

Background  

The potential causes for variation in virulence between distinct M. tuberculosis strains are still not fully known. However, differences in protein expression are probably an important factor. In this study we used a label-free quantitative proteomic approach to estimate differences in protein abundance between two closely related M. tuberculosis strains; the virulent H37Rv strain and its attenuated counterpart H37Ra.     

Variation among Genome Sequences of H37Rv Strains of Mycobacterium tuberculosis from Multiple Laboratories

The publication of the complete genome sequence for Mycobacterium tuberculosis H37Rv in 1998 has had a great impact on the research community. Nonetheless, it is suspected that genetic differences have arisen in stocks of H37Rv that are maintained in different laboratories. In order to assess the consistency of the genome sequences among H37Rv strains in use and the extent to which they have diverged from the original strain sequenced, we carried out whole-genome sequencing on six strains of H37Rv from different laboratories. Polymorphisms at 73 sites were observed, which were shared among the lab strains, though 72 of these were also shared with H37Ra and are likely to be due to sequencing errors in the original H37Rv reference sequence. An updated H37Rv genome sequence should be valuable to the tuberculosis research community as well as the broader microbial research community. In addition, several polymorphisms unique to individual strains and several shared polymorphisms were identified and shown to be consistent with the known provenance of these strains. Aside from nucleotide substitutions and insertion/deletions, multiple IS6110 transposition events were observed, supporting the theory that they play a significant role in plasticity of the M. tuberculosis genome. This genome-wide catalog of genetic differences can help explain any phenotypic differences that might be found, including a frameshift mutation in the mycocerosic acid synthase gene which causes two of the strains to be deficient in biosynthesis of the surface glycolipid phthiocerol dimycocerosate (PDIM). The resequencing of these six lab strains represents a fortuitous “in vitro evolution” experiment that demonstrates how the M. tuberculosis genome continues to evolve even in a controlled environment.Publication of the whole genome sequence of the H37Rv strain of Mycobacterium tuberculosis by Stewart Cole and colleagues in 1998 provided a breakthrough in tuberculosis (TB) research (8), leading to insights into the biology, metabolism, and evolution of this infectious pathogen. Large protein families related to fatty acid and polyketide biosynthesis, regulation (e.g., sigma factors and two-component sensor systems), drug efflux pumps and transporters, and the PE_PGRS proteins (a large duplicated family unique to the M. tuberculosis group of mycobacteria) were identified. In addition, transposons, prophage-like elements, and other repetitive and/or mobile genetic elements were identified (18). This genomic information has played an essential role in interpreting gene expression studies, modeling persistence, and identifying essential proteins as putative targets for drug discovery. However, to date the functions of only half of the genes (1,756/4,066) have been determined or predicted, and the rest remain annotated as “hypothetical proteins” (6).The H37Rv strain was initially selected for sequencing because it is a widely used laboratory strain that has retained its virulence. H37Rv was initially derived from a clinical isolate, H37, obtained from a patient with pulmonary tuberculosis in 1905. H37Rv falls in the T clade (5) and single-nucleotide polymorphism (SNP) cluster group SCG-6b (12). The virulence of H37Rv can be demonstrated in a number of animal models. For example, SCID mice infected with H37Rv typically have a mean time to death of 30 to 35 days, depending on the dose and route of inoculation (13).An avirulent strain, H37Ra, was also derived from H37 by culturing on solid egg medium and selecting for resistance to lysis (42). The strain was found not to cause disease in guinea pigs (43) or in mice (27). It has a colony morphology (smooth) different from that of H37Rv (rough) and several other phenotypic differences (14, 29). The H37Rv (ATCC 25618) and H37Ra (ATCC 25177) strains are maintained at the Trudeau Institute in New York (3), although unfortunately, the original H37 clinical isolate has been lost. Strain ATCC 27294 (TMC 102) is also frequently used as a representative of H37Rv in studies and treated equivalently in the literature. ATCC 25618 and ATCC 27294 were both isolated from the same patient in different years, and both are fully drug susceptible.The complete genome of H37Ra has been sequenced by Zheng et al. (48), who found 272 polymorphisms compared to the genome sequence determined by Cole et al. (8) for H37Rv. However, a subset of the polymorphic sites were found to match CDC1551, and upon resequencing of 85 such sites in H37Rv, 79 were determined to be errors in the H37Rv reference sequence. In addition, H37Ra has insertions of IS6110 at two novel sites and a loss of one, compared to the 16 sites in H37Rv. The 130 genuine H37Ra-specific polymorphisms found were divided into those in coding regions, those in upstream regulatory regions, and those in noncoding, nonregulatory intergenic regions in order to assess potential relevance to virulence. Polymorphisms in the promoter regions of sigC, nrdH (glutaredoxin-like electron transporter), and pabB (para-amino benzoate synthase), as well as nonsynonymous substitutions in mazG (regulator of stringent response), phoP (two-component sensor regulating biosynthesis of cell surface lipid antigens), pks12 (polyketide synthase involved in biosynthesis of mycoketides), and nrp (nonribosomal peptide synthetase potentially involved in phthiocerol dimycocerosate [PDIM] biosynthesis), were highlighted as possible causes of the loss of virulence. H37Ra does not synthesize a number of cell surface antigens, including sulfolipid-1, trehalose mycolates, and PDIM (7). The roles of mutations in phoP and sigM, both of which regulate expression of genes involved in biosynthesis of cell surface antigens, have been subsequently investigated, though neither seems to be singularly responsible for the avirulence of H37Ra (17, 35). Multiple mutations in PPE and PE_PGRS genes are also observed in H37Ra, and there has been speculation about the role of these genes in virulence (39). However, the RvD2 region (an 8-kb region present in H37Ra but deleted in H37Rv, including an IS6110 insertion element, mmpL14, and several hypothetical genes) is known not to be responsible for differences in virulence (25).Because of its importance as a model strain used in laboratory studies, it is essential to determine how consistent different stocks of H37Rv in different laboratories are with the reference genome sequence and with each other. Different stocks could accumulate independent polymorphisms over time, and such inconsistencies could potentially make results of studies obtained with H37Rv cultures from different labs difficult to compare, particularly if they affect virulence, drug tolerance, metabolism, cell wall constitution, etc. Furthermore, sequencing errors in the original genome sequence are possible. In order to evaluate differences among currently used variants of H37Rv, we resequenced the complete genomes of six extant H37Rv strains (two samples of ATCC 25618 and four of ATCC 27294) using Illumina sequencing technology. We compared differences among them and differences from the reference sequences for H37Rv and H37Ra available from GenBank. The results of this study identify a common set of 73 polymorphisms shared among all six sequenced strains relative to the H37Rv reference strain. Most (72) of these are shared with H37Ra and likely correspond to sequencing errors in the original H37Rv genome sequence. However, there are several sites where additional polymorphisms are shared among a subset of strains, and several strains have a small number of unique polymorphisms. Furthermore, examination of insertion sites of the IS6110 transposable element reveals several changes that have occurred among these strains. These results illustrate the ongoing evolution of this strain and divergence from the sequenced reference strain of H37Rv and highlight the importance of understanding the genetic differences unique to the stock used in each laboratory.