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

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

Computational Prediction and Experimental Assessment of Secreted/Surface Proteins from Mycobacterium tuberculosis H37Rv

The mycobacterial cell envelope has been implicated in the pathogenicity of tuberculosis and therefore has been a prime target for the identification and characterization of surface proteins with potential application in drug and vaccine development. In this study, the genome of Mycobacterium tuberculosis H37Rv was screened using Machine Learning tools that included feature-based predictors, general localizers and transmembrane topology predictors to identify proteins that are potentially secreted to the surface of M. tuberculosis, or to the extracellular milieu through different secretory pathways. The subcellular localization of a set of 8 hypothetically secreted/surface candidate proteins was experimentally assessed by cellular fractionation and immunoelectron microscopy (IEM) to determine the reliability of the computational methodology proposed here, using 4 secreted/surface proteins with experimental confirmation as positive controls and 2 cytoplasmic proteins as negative controls. Subcellular fractionation and IEM studies provided evidence that the candidate proteins Rv0403c, Rv3630, Rv1022, Rv0835, Rv0361 and Rv0178 are secreted either to the mycobacterial surface or to the extracellular milieu. Surface localization was also confirmed for the positive controls, whereas negative controls were located on the cytoplasm. Based on statistical learning methods, we obtained computational subcellular localization predictions that were experimentally assessed and allowed us to construct a computational protocol with experimental support that allowed us to identify a new set of secreted/surface proteins as potential vaccine candidates.     

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.     

Revisiting Host Preference in the Mycobacterium tuberculosis Complex: Experimental Infection Shows M. tuberculosis H37Rv to Be Avirulent in Cattle

Experiments in the late 19th century sought to define the host specificity of the causative agents of tuberculosis in mammals. Mycobacterium tuberculosis, the human tubercle bacillus, was independently shown by Smith, Koch, and von Behring to be avirulent in cattle. This finding was erroneously used by Koch to argue the converse, namely that Mycobacterium bovis, the agent of bovine tuberculosis, was avirulent for man, a view that was subsequently discredited. However, reports in the literature of M. tuberculosis isolation from cattle with tuberculoid lesions suggests that the virulence of M. tuberculosis for cattle needs to be readdressed. We used an experimental bovine infection model to test the virulence of well-characterized strains of M. tuberculosis and M. bovis in cattle, choosing the genome-sequenced strains M. tuberculosis H37Rv and M. bovis 2122/97. Cattle were infected with approximately 10⁶ CFU of M. tuberculosis H37Rv or M. bovis 2122/97, and sacrificed 17 weeks post-infection. IFN-γ and tuberculin skin tests indicated that both M. bovis 2122 and M. tuberculosis H37Rv were equally infective and triggered strong cell-mediated immune responses, albeit with some indication of differential antigen-specific responses. Postmortem examination revealed that while M. bovis 2122/97–infected animals all showed clear pathology indicative of bovine tuberculosis, the M. tuberculosis–infected animals showed no pathology. Culturing of infected tissues revealed that M. tuberculosis was able to persist in the majority of animals, albeit at relatively low bacillary loads. In revisiting the early work on host preference across the M. tuberculosis complex, we have shown M. tuberculosis H37Rv is avirulent for cattle, and propose that the immune status of the animal, or genotype of the infecting bacillus, may have significant bearing on the virulence of a strain for cattle. This work will serve as a baseline for future studies into the genetic basis of host preference, and in particular the molecular basis of virulence in M. bovis.     

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.     

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

Rv1915 and Rv1916 from Mycobacterium tuberculosis H37Rv form in vitro protein-protein complex

BackgroundMycobacterium tuberculosis (Mtb) isocitrate lyase (ICL) is an established drug target that facilitates Mtb persistence. Unlike other mycobacterial strains, where ICL2 is a single gene product, H37Rv has a split event, resulting in two tandemly coded icls - rv1915 and rv1916. Our recent report on functionality of individual Rv1915 and Rv1916, led to postulate the cooperative role of these proteins in pathogen's survival under nutrient-limiting conditions. This study investigates the possibility of Rv1915 and Rv1916 interacting and forming a complex.MethodsPull down assay, activity assay, mass spectrometry and site directed mutagenesis was employed to investigate and validate Rv1915-Rv1916 complex formation.ResultsRv1915 and Rv1916 form a stable complex in vitro, with enhanced ICL/MICL activities as opposed to individual proteins. Further, activities monitored in the presence of acetyl-CoA show significant increase for Rv1916 and the complex but not of Rv0467 and Rv1915Δ90CT. Both full length and truncated Rv1915Δ90CT can form complex, implying the absence of its C-terminal disordered region in complex formation. Further, in silico analysis and site-directed mutagenesis studies reveal Y64 and Y65 to be crucial residues for Rv1915-Rv1916 complex formation.ConclusionsThis study uncovers the association between Rv1915 and Rv1916 and supports the role of acetyl-CoA in escalating the ICL/MICL activities of Rv1916 and Rv1915Δ90CT-Rv1916 complex.General significancePartitioning of ICL2 into Rv1915 and Rv1916 that associates to form a complex in Mtb H37Rv, suggests its importance in signaling and regulation of metabolic pathway particularly in carbon assimilation.     

Pathways of carbohydrate metabolism in mycobacterium tuberculosis H37Rv1.

Radiorespirometric studies using glucose labelled at 1, 2, 3-4, and 6 positions and enzymatic studies were conducted to determine the primary pathways of glucose dissimilation in Mycobacterium tuberculosis H37Rv. The pattern of 14CO2 recovery was C3-4 greater than C1 greater than C6 = C2. The Embden-Meyerhof pathway was found to be the predominant pathway for glucose oxidation, operative to the extent of 94%. The pentose phosphate pathway accounted for the remaining 6%. Maximum incorporation of 14C into cellular components was from C2 and C6 labelled glucose.     

Up regulated virulence genes in M. tuberculosis H37Rv gleaned from genome wide expression profiles

Identification of up regulated virulence genes in M. tuberculosis H37Rv using genome wide expression profiles is of interest in drug discovery for the disease. Hence, we report 17 up-regulated PPIN (Protein-Protein Interaction Network) enriched potential virulence linked genes using expression data available at the Gene Expression Omnibus (GEO) database for further consideration.     

Mining Genomic Patterns in Mycobacterium tuberculosis H37Rv Using a Web Server Tuber-Gene

Mycobacterium tuberculosis (MTB), causative agent of tuberculosis, is one of the most dreaded diseases of the century. It has long been studied by researchers throughout the world using various wet-lab and dry-lab techniques. In this study, we focus on mining useful patterns at genomic level that can be applied for in silico functional characterization of genes from the MTB complex. The model developed on the basis of the patterns found in this study can correctly identify 99.77% of the input genes from the genome of MTB strain H37Rv. The model was tested against four other MTB strains and the homologue M. bovis to further evaluate its generalization capability. The mean prediction accuracy was 85.76%. It was also observed that the GC content remained fairly constant throughout the genome, implicating the absence of any pathogenicity island transferred from other organisms. This study reveals that dinucleotide composition is an efficient functional class discriminator for MTB complex. To facilitate the application of this model, a web server Tuber-Gene has been developed, which can be freely ac- cessed at http://www.bifmanit.org/tb2/.     

Fatty acid regulation of adenylyl cyclase Rv2212 from Mycobacterium tuberculosis H37Rv

Adenylyl cyclase Rv2212 from Mycobacterium tuberculosis has a domain composition identical to the pH-sensing isoform Rv1264, an N-terminal regulatory domain and a C-terminal catalytic domain. The maximal velocity of Rv2212 was the highest of all 10 mycobacterial cyclases investigated to date (3.9 micromol cAMP.mg(-1).min(-1)), whereas ATP substrate affinity was low (SC(50) = 2.1 mm ATP). Guanylyl cyclase side activity was absent. The activities and kinetics of the holoenzyme and of the catalytic domain alone were similar, i.e. in distinct contrast to the Rv1264 adenylyl cyclase, in which the N-terminal domain is autoinhibitory. Unsaturated fatty acids strongly stimulated Rv2212 activity by increasing substrate affinity. In addition, fatty acids greatly enhanced the pH sensitivity of the holoenzyme, thus converting Rv2212 to a pH sensor adenylyl cyclase. Fatty acid binding to Rv2212 was modelled by homology to a recent structure of the N-terminal domain of Rv1264, in which a fatty acid-binding pocket is defined. Rv2212 appears to integrate three cellular parameters: ATP concentration, presence of unsaturated fatty acids, and pH. These regulatory properties open the possibility that novel modes of cAMP-mediated signal transduction exist in the pathogen.     

Identification and characterization of Rv3281 as a novel subunit of a biotin-dependent acyl-CoA Carboxylase in Mycobacterium tuberculosis H37Rv

Mycobacterium tuberculosis produces a large number of structurally diverse lipids generated from the carboxylation products of acetyl-CoA and propionyl-CoA. A biotin-dependent acyl-CoA carboxylase was purified from M. tuberculosis H37Rv by avidin affinity chromatography, and the three major protein components were determined by N-terminal sequencing to be the 63-kDa alpha3-subunit (AccA3, Rv3285), the 59-kDa beta5-subunit (AccD5, Rv3280), and the 56-kDa beta4-subunit (AccD4, Rv3799). A minor protein of about 24 kDa that co-purified with the above subunits was identified by matrix-assisted laser desorption/ionization-time of flight mass spectrometry to be the product of Rv3281 that is located immediately downstream of the open reading frame encoding the beta5-subunit. This protein displays identity over a short stretch of amino acids with the recently discovered epsilon-subunits of Streptomyces coelicolor, suggesting that it might be an epsilon-subunit of the mycobacterial acyl-CoA carboxylase. To test this hypothesis, the carboxylase subunits were expressed in Escherichia coli and purified. Acyl-CoA carboxylase activity was successfully reconstituted for the first time from purified subunits of the acyl-CoA carboxylase of M. tuberculosis. The reconstituted alpha3-beta5 showed higher activity with propionyl-CoA than with acetyl-CoA, and the addition of the epsilon-subunit stimulated the carboxylation by 3.2- and 6.3-fold, respectively. The alpha3-beta4 showed very low activity with the above substrates but carboxylated long chain acyl-CoA. This epsilon-subunit contains five sets of tandem repeats at the N terminus that are required for maximal enhancement of carboxylase activity. The Rv3281 open reading frame is co-transcribed with Rv3280 in the mycobacterial cell, and the level of epsilon-protein was highest during the log phase and decreased during the stationary phase.     

Characterization of an acid inducible lipase Rv3203 from Mycobacterium tuberculosis H37Rv

The Rv3203 (LipV) of Mycobacterium tuberculosis (Mtb) H37Rv, is annotated as a member of Lip family based on the presence of characteristic consensus esterase motif ‘GXSXG’. In vitro culture studies of Mtb H37Ra indicated that expression of Rv3203 gene was up-regulated during acidic stress as compared to normal whereas no expression was observed under nutrient and oxidative stress conditions. Therefore, detailed characterization of Rv3203 was done by gene cloning and its further expression and purification as his-tagged protein in microbial expression system. The enzyme was purified to homogeneity by affinity chromatography. It demonstrated broad substrate specificity and preferentially hydrolyzed p-nitrophenyl myristate. The purified enzyme demonstrated an optimum activity at pH 8.0 and temperature 50 °C. The specific activity, K _m and V _max of enzyme was determined to be 21.29 U mg^?1 protein, 714.28 μM and 62.5 μmol ml^?1 min^?1, respectively. The pH stability assay and circular dichroism spectroscopic analysis revealed that Rv3203 protein is more stable in acidic condition. Tetrahydrolipstatin, a specific lipase inhibitor and RHC80267, a diacylglycerol lipase inhibitor abolished the activity of this enzyme. The catalytic triad residues were determined to be Ser⁵⁰, Asp¹⁸⁰ and His²⁰³ residues by site-directed mutagenesis.