利用抑制消减杂交技术研究结核分支杆菌强毒株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收录。     

结核分枝杆菌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的毒力丧失有关。     

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

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

Comparative genomic and proteomic analyses of PE/PPE multigene family of Mycobacterium tuberculosis H37Rv and H37Ra reveal novel and interesting differences with implications in virulence

Tuberculosis, caused by Mycobacterium tuberculosis, remains a leading infectious disease taking one human life every 15 s globally. The two well-characterized strains H₃₇Rv and H₃₇Ra, derived from the same parental strain M. tuberculosis H₃₇, show dramatically different pathogenic phenotypes. PE/PPE gene family, comprising of 176 open reading frames and present exclusively in genus Mycobacterium, accounts for ∼10% of the M. tuberculosis genome. Our comprehensive in silico analyses of PE/PPE family of H₃₇Ra and virulent H₃₇Rv strains revealed genetic differences between these strains in terms of several single nucleotide variations and InDels and these manifested in changes in physico-chemical properties, phosphorylation sites, and protein: protein interacting domains of the corresponding proteomes. Similar comparisons using the 13 sigma factor genes, 36 members of the mammalian cell entry family, 13 mycobacterial membrane protein large family members and 11 two-component signal transduction systems along with 5 orphaned response regulators and 2 orphaned sensor kinases failed to reveal very significant difference between H₃₇Rv and H₃₇Ra, reinforcing the importance of PE/PPE genes. Many of these changes between H₃₇Rv and H₃₇Ra can be correlated to differences in pathogenesis and virulence of the two strains.     

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.     

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

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

Genes of Mycobacterium tuberculosis H37Rv downregulated in the attenuated strain H37Ra are restricted to M. tuberculosis complex species

By comparing gene expression of virulent Mycobacterium tuberculosis H37Rv and attenuated strain H37Ra, we previously detected six genes that appear to be markedly downregulated in the attenuated strain compared with the virulent one. Three of these genes, i.e. Rv1345, Rv2770c, and Rv0288, code for proteins that can be predictively associated to immunological or pathogenetic aspects of M. tuberculosis infection; the other genes, i.e. Rv2336, Rv1320c, and Rv2819c, code for proteins with unknown functions (Rindi et al., 1999). In this paper we searched for the above mentioned genes in Pvu II-digested genomic DNA of a number of mycobacterial species by southern blot analysis employing PCR-generated probes in high-stringency conditions. Hybridization signals were only found in species belonging to the M. tuberculosis complex, i.e., M. tuberculosis, M. bovis, including the BCG strain, and M. microti, but not in other mycobacterial species, including M. avium, M. intracellulare, M. malmoense, M. xenopi, M. kansasii, M. simiae, M. marinum, M. scrofulaceum, M. gordonae, M. fortuitum, and M. smegmantis. These results indicate that genes Rv1345, Rv2770c, Rv0288, Rv2336, Rv1320c, and Rv2819c are associated with the most virulent mycobacteria and further support their potential role in M. tuberculosis virulence.     

Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains

Virulence and immunity are poorly understood in Mycobacterium tuberculosis. We sequenced the complete genome of the M. tuberculosis clinical strain CDC1551 and performed a whole-genome comparison with the laboratory strain H37Rv in order to identify polymorphic sequences with potential relevance to disease pathogenesis, immunity, and evolution. We found large-sequence and single-nucleotide polymorphisms in numerous genes. Polymorphic loci included a phospholipase C, a membrane lipoprotein, members of an adenylate cyclase gene family, and members of the PE/PPE gene family, some of which have been implicated in virulence or the host immune response. Several gene families, including the PE/PPE gene family, also had significantly higher synonymous and nonsynonymous substitution frequencies compared to the genome as a whole. We tested a large sample of M. tuberculosis clinical isolates for a subset of the large-sequence and single-nucleotide polymorphisms and found widespread genetic variability at many of these loci. We performed phylogenetic and epidemiological analysis to investigate the evolutionary relationships among isolates and the origins of specific polymorphic loci. A number of these polymorphisms appear to have occurred multiple times as independent events, suggesting that these changes may be under selective pressure. Together, these results demonstrate that polymorphisms among M. tuberculosis strains are more extensive than initially anticipated, and genetic variation may have an important role in disease pathogenesis and immunity.     

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.     

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.     

MTB-PCDB: Mycobacterium tuberculosis proteome comparison database

The Mycobacterium tuberculosis Proteome Comparison Database (MTB-PCDB) is an online database providing integrated access to proteome sequence comparison data for five strains of Mycobacterium tuberculosis (H37Rv, H37Ra, CDC 1551, F11 and KZN 1435) sequenced completely so far. MTB-PCDB currently hosts 40252 protein sequence comparison data obtained through inter-strain proteome comparison of five different strains of MTB. 2373 proteins were found to be identical in all 5 strains using MTB H(37)Rv as reference strain. To enable wide use of this data, MTB-PCDB provides a set of tools for searching, browsing, analyzing and downloading the data. By bringing together, M. tuberculosis proteome comparison among virulent & avirulent strains and also drug susceptible & drug resistance strains MTB-PCDB provides a unique discovery platform for comparative proteomics among these strains which may give insights into the discovery & development of TB drugs, vaccines and biomarkers. AVAILABILITY: The database is available for free at http://www.bicjbtdrc-mgims.in/MTB-PCDB/     

结核分枝杆菌蛋白Rv3425对细菌表型及毒力影响的研究

为探索蛋白Rv3425在结核分枝杆菌(Mycobacterium tuberculosis,M. tuberculosis)中的功能,本研究以耻垢分枝杆菌(Mycobacterium smegmatis,M. smegmatis)为模式菌株,构建重组了耻垢分枝杆菌Ms-Rv3425。分别将构建菌株(Ms-Rv3425)、野生株(Ms)及空载对照(Ms-Pact)接种于7H9-OADC培养基中37 ℃培养,观察Ms-Rv3425与Ms及Ms-Pact之间在生长速率、菌落形态、生物膜以及聚集度方面的差异。分别用低pH值以及含有十二烷基磺酸钠(sodium dodecyl sulfate,SDS)、氨苄西林、异烟肼及利福平的培养基进行培养,计算存活率以分析抗逆和抗药能力;用上述压力条件培养结核分枝杆菌标准株H37Ra,分析Rv3425内源表达量的变化;进行THP-1细胞感染和BALB/c小鼠攻毒实验分析菌株的毒性变化。结果显示,与Ms及Ms-Pact相比,Ms-Rv3425的菌落形态更为粗糙且隆起,成膜及聚集能力增强;在压力条件下,Ms-Rv3425表现出更高的抗逆和抗药能力,H37Ra中Rv3425的表达量也显著上调;胞内存活率及小鼠致死率更高,各脏器病理损伤更为严重。综上所述,过表达Rv3425能够改变耻垢分枝杆菌的表型,提高抗逆性、抗药性和毒力。深入探讨PPE家族蛋白Rv3425的功能,将为结核病的防治带来新的视角。     

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.     

Requirement of gene fadD33 for the growth of Mycobacterium tuberculosis in a hepatocyte cell line

Gene fadD33 of Mycobacterium tuberculosis, one of the 36 homologues of gene fadD of Escherichia coli identified in the M. tuberculosis genome, predictively encodes an acyl-CoA synthase, an enzyme involved in fatty acids metabolism. The gene is underexpressed in the attenuated strain M. tuberculosis H37Ra relative to virulent H37Rv and plays a role in M. tuberculosis virulence in BALB/c mice by supporting mycobacterial replication in the liver. In the present paper, we investigated the role of fadD33 expression in bacterial growth within the hepatocyte cell line HepG2, as well as in human monocyte-derived THP-1 cells and peripheral blood mononuclear cells. M. tuberculosis H37Rv proved able to grow within HepG2 cells, while the intracellular replication of M. tuberculosis H37Ra was markedly impaired; complementation of strain H37Ra with gene fadD33 restored its replication to the levels of H37Rv. Moreover, disruption of gene fadD33 by allelic exchange mutagenesis reduced the intracellular growth of M. tuberculosis H37Rv, and complementation of the fadD33-disrupted mutant with gene fadD33 restored bacterial replication. Conversely, fadD33 expression proved unable to influence M. tuberculosis growth in human phagocytes, as fadD33-disrupted M. tuberculosis H37Rv mutant, as well as fadD33-complemented M. tuberculosis H37Ra, grew within THP-1 cells and peripheral monocytes basically at the same rates as parent H37Rv and H37Ra strains. The results of these experiments indicate that gene fadD33 expression confers growth advantage to M. tuberculosis in immortalized hepatocytes, but not in macrophages, thus emphasizing the importance of fadD33 in liver-specific replication of M. tuberculosis.     

Comparative profiles of intramacrophage behavior of Mycobacterium tuberculosis and Mycobacterium avium complex with different levels of virulence

Mycobacterium tuberculosis (MTB) and M. avium complex (MAC) strains with different levels of virulence in mice were examined for profiles of interaction with murine peritoneal macrophages (Mphis). Their growth rates in Mphis were in these orders: H37Ra strain (attenuated) > H37Rv strain (virulent) for MTB, and N-260 strain (moderate virulence) > MAC N-444 strain (low virulence) for MAC. MTB but not MAC caused the necrotic death of host Mphis in terms of increased release of lactate dehydrogenase from infected Mphis. The MTB H37Ra strain induced a greater production of reactive nitrogen intermediates (RNI) by Mphis than the MTB H37Rv strain did. However, this phenomenon was not observed with MAC, implying less important roles of RNI in the expression of Mphi antimicrobial activity against MAC organisms.     

Mycobacterium tuberculosis CDC1551 induces a more vigorous host response in vivo and in vitro, but is not more virulent than other clinical isolates.

Mycobacterium tuberculosis CDC1551, a clinical isolate reported to be hypervirulent and to grow faster than other isolates, was compared with two other clinical isolates (HN60 and HN878) and two laboratory strains (H37Rv and Erdman). The initial (1-14 days) growth of CDC1551, HN60, HN878, and H37Rv was similar in the lungs of aerosol-infected mice, but growth of Erdman was slower. Thereafter, the growth rate of CDC1551 decreased relative to the other strains which continued to grow at comparable rates up to day 21. In the lungs of CDC1551-infected mice, small well-organized granulomas with high levels of TNF-alpha, IL-6, IL-10, IL-12, and IFN-gamma mRNA were apparent sooner than in lungs of mice infected with the other strains. CDC1551-infected mice survived significantly longer. These findings were confirmed in vitro. The growth rates of H37Rv and CDC1551 in human monocytes were the same, but higher levels of TNF-alpha, IL-10, IL-6, and IL-12 were induced in monocytes after infection with CDC1551 or by exposure of monocytes to lipid fractions from CDC1551. CD14 expression on the surface of the monocytes was up-regulated to a greater extent by exposure to the lipids of CDC1551. Thus, CDC1551 is not more virulent than other M. tuberculosis isolates in terms of growth in vivo and in vitro, but it induces a more rapid and robust host response.