Proteogenomic analysis of Mycobacterium tuberculosis by high resolution mass spectrometry

The genome sequencing of H37Rv strain of Mycobacterium tuberculosis was completed in 1998 followed by the whole genome sequencing of a clinical isolate, CDC1551 in 2002. Since then, the genomic sequences of a number of other strains have become available making it one of the better studied pathogenic bacterial species at the genomic level. However, annotation of its genome remains challenging because of high GC content and dissimilarity to other model prokaryotes. To this end, we carried out an in-depth proteogenomic analysis of the M. tuberculosis H37Rv strain using Fourier transform mass spectrometry with high resolution at both MS and tandem MS levels. In all, we identified 3176 proteins from Mycobacterium tuberculosis representing ~80% of its total predicted gene count. In addition to protein database search, we carried out a genome database search, which led to identification of ~250 novel peptides. Based on these novel genome search-specific peptides, we discovered 41 novel protein coding genes in the H37Rv genome. Using peptide evidence and alternative gene prediction tools, we also corrected 79 gene models. Finally, mass spectrometric data from N terminus-derived peptides confirmed 727 existing annotations for translational start sites while correcting those for 33 proteins. We report creation of a high confidence set of protein coding regions in Mycobacterium tuberculosis genome obtained by high resolution tandem mass-spectrometry at both precursor and fragment detection steps for the first time. This proteogenomic approach should be generally applicable to other organisms whose genomes have already been sequenced for obtaining a more accurate catalogue of protein-coding genes.     

Independent large scale duplications in multiple M. tuberculosis lineages overlapping the same genomic region

Mycobacterium tuberculosis, the causative agent of most human tuberculosis, infects one third of the world's population and kills an estimated 1.7 million people a year. With the world-wide emergence of drug resistance, and the finding of more functional genetic diversity than previously expected, there is a renewed interest in understanding the forces driving genome evolution of this important pathogen. Genetic diversity in M. tuberculosis is dominated by single nucleotide polymorphisms and small scale gene deletion, with little or no evidence for large scale genome rearrangements seen in other bacteria. Recently, a single report described a large scale genome duplication that was suggested to be specific to the Beijing lineage. We report here multiple independent large-scale duplications of the same genomic region of M. tuberculosis detected through whole-genome sequencing. The duplications occur in strains belonging to both M. tuberculosis lineage 2 and 4, and are thus not limited to Beijing strains. The duplications occur in both drug-resistant and drug susceptible strains. The duplicated regions also have substantially different boundaries in different strains, indicating different originating duplication events. We further identify a smaller segmental duplication of a different genomic region of a lab strain of H37Rv. The presence of multiple independent duplications of the same genomic region suggests either instability in this region, a selective advantage conferred by the duplication, or both. The identified duplications suggest that large-scale gene duplication may be more common in M. tuberculosis than previously considered.     

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

Hemolytic phospholipase Rv0183 of Mycobacterium tuberculosis induces inflammatory response and apoptosis in alveolar macrophage RAW264.7 cells

The metabolic pathway of phospholipids is one of the most important physiologic pathways in Mycobacterium tuberculosis, a typical intracellular bacterium. The hemolytic phospholipase lip gene (Rv0183) is one of 24 phospholipase genes that have been demonstrated to play critical roles in the metabolism of phospholipids in M.?tuberculosis. Quantitative RT-PCR and flow cytometry were used to elucidate the immunological and pathogenic implications of the Rv0183 gene on the inflammatory response following persistent expression of Rv0183 in mouse alveolar macrophage RAW264.7 cells. Our results demonstrate that a time-course-dependent ectopic expression of Rv0183 significantly elevated the expression of IL-6, NF-κB, TLR-2, TLR-6, TNFα, and MyD88 in these alveolar macrophage cells. Furthermore, the persistent expression of Rv0183 induced RAW264.7 cell apoptosis in vitro. These findings demonstrate that the expression of Rv0183 induces an inflammatory response and cell apoptosis in the host cells, suggesting that Rv0183 may play an important role in the virulence and pathogenesis of M.?tuberculosis infection.     

A fragment of 21 ORFs around the direct repeat (DR) region of Mycobacterium tuberculosis is absent from the other sequenced mycobacterial genomes: implications for the evolution of the DR region

The direct repeat (DR) region is a singular locus of the Mycobacterium tuberculosis complex genome. This region consists of 36 bp repetitive sequences separated by non-repetitive unique spacer sequences. Around this region there are several genes coding for proteins of unknown function. To determine whether the M. smegmatis, M. avium, M. marinum and M. leprae genomes contain sequences and ORFs similar to those of the DR locus of the M. tuberculosis complex, we analysed the corresponding regions in these species. As a first step, some conserved genes that flank the DR genes [Rv2785c (rpsO), Rv2786c (ribF), Rv2790c (ltp1 ), Rv2793c (truB), Rv2800, Rv2825, Rv2828, Rv2831 (echA16 ), Rv2838 (rbfA) and Rv2845 (proS )] were used as markers to locate the corresponding orthologues in M. smegmatis, M. avium, M. marinum and M. leprae in silico. Most of these M. tuberculosis marker genes have highly similar orthologues located in the same order and orientation in the other mycobacteria. In contrast, no orthologues were found for ORFs Rv2801-Rv2824, suggesting that these genes are unique to M. tuberculosis within the genus Mycobacterium.We observed that in M. smegmatis and M. avium, Rv2800 and Rv2825 are adjacent. This observation was experimentally confirmed by PCR. In conclusion, as the DR locus and the ORFs around it are absent in M. smegmatis and M. avium and, as it is possible that these species are older than M. tuberculosis, we postulated that the DR locus was acquired by the M. tuberculosis complex species or by an ancestor bacterium.     

Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays

Whole-genome comparisons of the tubercle bacilli were undertaken using ordered bacterial artificial chromosome (BAC) libraries of Mycobacterium tuberculosis and the vaccine strain, Mycobacterium bovis BCG-Pasteur, together with the complete genome sequence of M. tuberculosis H37Rv. Restriction-digested BAC arrays of M. tuberculosis H37Rv were used in hybridization experiments with radiolabelled M. bovis BCG genomic DNA to reveal the presence of 10 deletions (RD1-RD10) relative to M. tuberculosis. Seven of these regions, RD4-RD10, were also found to be deleted from M. bovis, with the three M. bovis BCG-specific deletions being identical to the RD1-RD3 loci described previously. The distribution of RD4-RD10 in Mycobacterium africanum resembles that of M. tuberculosis more closely than that of M. bovis, whereas an intermediate arrangement was found in Mycobacterium microti, suggesting that the corresponding genes may affect host range and virulence of the various tubercle bacilli. Among the known products encoded by these loci are a copy of the proposed mycobacterial invasin Mce, three phospholipases, several PE, PPE and ESAT-6 proteins, epoxide hydrolase and an insertion sequence. In a complementary approach, direct comparison of BACs uncovered a third class of deletions consisting of two M. tuberculosis H37Rv loci, RvD1 and RvD2, deleted from the genome relative to M. bovis BCG and M. bovis. These deletions affect a further seven genes, including a fourth phospholipase, plcD. In summary, the insertions and deletions described here have important implications for our understanding of the evolution of the tubercle complex.     

差显技术分析结核杆菌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.     

A mycobacterial virulence gene cluster extending RD1 is required for cytolysis, bacterial spreading and ESAT-6 secretion

Initiation and maintenance of infection by mycobacteria in susceptible hosts are not well understood. A screen of Mycobacterium marinum transposon mutant library led to isolation of eight mutants that failed to cause haemolysis, all of which had transposon insertions in genes homologous to a region between Rv3866 and Rv3881c in Mycobacterium tuberculosis, which encompasses RD1 (Rv3871-Rv3879c), a known virulence gene cluster. The M. marinum mutants showed decreased virulence in vivo and failed to secrete ESAT-6, like M. tuberculosis RD1 mutants. M. marinum mutants in genes homologous to Rv3866-Rv3868 also failed to accumulate intracellular ESAT-6, suggesting a possible role for those genes in synthesis or stability of the protein. These transposon mutants and an ESAT-6/CFP-10 deletion mutant all showed reduced cytolysis and cytotoxicity to macrophages and significantly decreased intracellular growth at late stages of the infection only when the cells were infected at low multiplicity of infection, suggesting a defect in spreading. Direct evidence for cell-to-cell spread by wild-type M. marinum was obtained by microscopic detection in macrophage and epithelial monolayers, but the mutants all were defective in this assay. Expression of M. tuberculosis homologues complemented the corresponding M. marinum mutants, emphasizing the functional similarities between M. tuberculosis and M. marinum genes in this region that we designate extRD1 (extended RD1). We suggest that diminished membranolytic activity and defective spreading is a mechanism for the attenuation of the extRD1 mutants. These results extend recent findings on the genomic boundaries and functions of M. tuberculosis RD1 and establish a molecular cellular basis for the role that extRD1 plays in mycobacterial virulence. Disruption of the M. marinum homologue of Rv3881c, not previously implicated in virulence, led to a much more attenuated phenotype in macrophages and in vivo, suggesting that this gene plays additional roles in M. marinum survival in the host.     

3-Ketosteroid 9α-hydroxylase is an essential factor in the pathogenesis of Mycobacterium tuberculosis

Mycobacterium tuberculosis H37Rv contains the kshA ( Rv3526 ) and kshB ( Rv3571 ) genes, encoding 3-ketosteroid 9α-hydroxylase (KSH). Consistent with their predicted roles, the Δ kshA and Δ kshB deletion mutants of M. tuberculosis H37Rv were unable to use cholesterol and 4-androstene-3,17-dione as primary carbon and energy sources. Interestingly, Δ kshA and Δ kshB mutants were also unable to metabolize the steroid substrate 5α-androstane-3,17-dione, whereas wild-type M. tuberculosis H37Rv could. The deletion of either of these genes lead to rapid death of the microorganism in murine infection models and in macrophages, showing that kshA and kshB are essential factors for M. tuberculosis pathogenesis. Penta-acylated trehalose (PAT) biosynthesis was altered in the Δ kshB mutant, but not the Δ kshA mutant. The Δ kshB mutant synthesizes all other types of lipids. The Δ kshB mutant had a thickened outer layer in its cell wall. KshB thus appears to be involved in multiple processes, probably as a reductase of different oxygenases. We conclude that an impaired 3-ketosteroid 9α-hydroxylase activity is the cause of the highly attenuated phenotype of our M. tuberculosis H37Rv mutants.     

A new high-throughput AFLP approach for identification of new genetic polymorphism in the genome of the clonal microorganism Mycobacterium tuberculosis

We have here applied high-throughput amplified fragment length polymorphism (htAFLP) analysis to strains belonging to the five classical species of the Mycobacterium tuberculosis complex. Using 20 strains, three enzyme combinations and eight selective amplification primer pairs, 24 AFLP reactions were performed per strain. Overall, this resulted in 480 DNA fingerprints and more than 1200 htAFLP-amplified PCR fragments were visualised per strain. The cumulative dendrogram correctly clustered strains from the various species, albeit within a distance of 6.5% for most of them. The single isolate of Mycobacterium canettii presented separately at 19% distance. All over, 169 fragments (14%) appeared to be polymorphic. Sixty-eight were specific for M. canetti and forty-five for Mycobacterium bovis. For the 10 different M. tuberculosis strains included in the present analysis, 56 polymorphic markers were identified. Upon sequencing 20 of these marker regions and comparisons with the H37Rv genome sequence, 25% appeared to share homology to members of the antigenically variable PE/PPE surface protein encoding gene family confirming previous findings on the genetic heterogeneity within these genes. In addition, homologues for phage genes and insertion element-encoded genes were detected. Forty-five percent of the sequences derived from ORFs with a currently unknown function, which was corroborated by genome sequence comparison for the clinical M. tuberculosis CD 1551 isolate. Sequence variation in M. tuberculosis was assessed in more detail for a subset of these loci by newly designed PCR restriction fragment length polymorphism (RFLP) tests and direct sequencing. Fourteen novel PCR RFLP tests were developed and twelve novel single nucleotide polymorphisms (SNPs) were identified, all suited for epidemiological analysis of M. tuberculosis. The tests allowed for identification of the major Mycobacterium species and M. tuberculosis variants and clones.     

tRNA genes in mycobacteria: organization and molecular cloning.

DNAs from nine mycobacteria cleaved with restriction endonucleases were hybridized with cDNA probes synthesized to tRNAs from Mycobacterium tuberculosis and Mycobacterium smegmatis. The tRNA genes are conserved, but their gross genomic organization has diverged in six of the nine species examined. Organisms of the M. tuberculosis H37Ra and H37Rv-M. bovis BCG complex appeared to have identical tRNA gene organization and were indistinguishable from each other. M. tuberculosis and M. smegmatis tRNA-derived cDNA probes hybridized differentially to tRNA-coding DNA segments in five of the species examined, suggesting the existence of qualitatively different tRNA pools in these slow- and fast-growing species. Mycobacterial DNAs hybridized with cDNA synthesized to 23S plus 16S rRNAs from Escherichia coli, and the data suggested that the tRNA genes map close to the rRNA genes. A gene bank of M. tuberculosis H37Rv DNA was constructed, and a recombinant plasmid, pSB2, coding for tRNA(s) and rRNA(s) was partially characterized. Plasmid pSB2 recognized a SalI restriction fragment length polymorphism (RFLP) in M. tuberculosis H37Rv and H37Ra; however, the RFLP is not linked to the tRNA-coding region. To the best of our knowledge, this is the first report of an RFLP which distinguishes the pathogenic strain M. tuberculosis H37Rv from its avirulent derivative H37Ra.     

In silico analysis of DosR regulon proteins of Mycobacterium tuberculosis

One of the challenges faced by Mycobacterium tuberculosis (M. tuberculosis) in dormancy is hypoxia. DosR/DevR of M. tuberculosis is a two component dormancy survival response regulator which induces the expression of 48 genes. In this study, we have used DosR regulon proteins of M. tuberculosis H37Rv as the query set and performed a comprehensive homology search against the non-redundant database. Homologs were found in environmental mycobacteria, environmental bacteria and archaebacteria. Analysis of genomic context of DosR regulon revealed that they are distributed as nine blocks in the genome of M. tuberculosis with many transposases and integrases in their vicinity. Further, we classified DosR regulon proteins into eight functional categories. One of the hypothetical proteins Rv1998c could probably be a methylisocitrate lyase or a phosphonomutase. Another hypothetical protein, Rv0572 was found only in mycobacteria. Insights gained in this study can potentially aid in the development of novel therapeutic interventions.     

The crystal structure of Rv0793, a hypothetical monooxygenase from <Emphasis Type="BoldItalic">M.␣tuberculosis</Emphasis>

Mycobacterium tuberculosis infects millions worldwide. The Structural Genomics Consortium for M. tuberculosis has targeted all genes from this bacterium in hopes of discovering and developing new therapeutic agents. Open reading frame Rv0793 from M. tuberculosis was annotated with an unknown function. The 3-dimensional structure of Rv0793 has been solved to 1.6 A resolution. Its structure is very similar to that of Streptomyces coelicolor ActVA-Orf6, a monooxygenase that participates in tailoring of polyketide antibiotics in the absence of a cofactor. It is also similar to the recently solved structure of YgiN, a quinol monooxygenase from Escherichia coli. In addition, the structure of Rv0793 is similar to several structures of other proteins with unknown function. These latter structures have been determined recently as a result of structural genomic projects for various bacterial species. In M. tuberculosis, Rv0793 and its homologs may represent a class of monooygenases acting as reactive oxygen species scavengers that are essential for evading host defenses. Since the most prevalent mode of attack by the host defense on M. tuberculosis is by reactive oxygen species and reactive nitrogen species, Rv0793 may provide a novel target to combat infection by M. tuberculosis.     

Comparative proteome analysis of Mycobacterium tuberculosis and Mycobacterium bovis BCG strains: towards functional genomics of microbial pathogens

In 1993, the WHO declared tuberculosis a global emergency on the basis that there are 8 million new cases per year. The complete genome of the strain H37Rv of the causative microorganism, Mycobacterium tuberculosis, comprising 3924 genes has been sequenced. We compared the proteomes of two non-virulent vaccine strains of M. bovis BCG (Chicago and Copenhagen) with two virulent strains of M. tuberculosis (H37Rv and Erdman) to identify protein candidates of value for the development of vaccines, diagnostics and therapeutics. The mycobacterial strains were analysed by two-dimensional electrophoresis (2-DE) combining non-equilibrium pH gradient electrophoresis (NEPHGE) with SDS-PAGE. Distinct and characteristic proteins were identified by mass spectrometry and introduced into a dynamic 2-DE database (http://www.mpiib-berlin.mpg.de/2D-PAGE). Silver-stained 2-DE patterns of mycobacterial cell proteins or culture supernatants contained 1800 or 800 spots, respectively, from which 263 were identified. Of these, 54 belong to the culture supernatant. Sixteen and 25 proteins differing in intensity or position between M. tuberculosis H37Rv and Erdman, and H37Rv and M. bovis BCG Chicago, respectively, were identified and categorized into protein classes. It is to be hoped that the availability of the mycobacterial proteome will facilitate the design of novel measures for prevention and therapy of one of the great health threats, tuberculosis.