Integration and excision of the Mycobacterium tuberculosis prophage-like element,phiRv1

The genomes of Mycobacterium tuberculosis H37Rv and CDC1551 each contain two prophage-like elements, phiRv1 and phiRv2. The phiRv1 element is not only absent from Mycobacterium bovis BCG but is in different locations within the two sequenced M. tuberculosis genomes; in both cases phiRv1 is inserted into a REP13E12 repeated sequence, which presumably contains the bacterial attachment site, attB, for phiRv1. Although phiRv1 is probably too small to encode infectious phage particles, it may nevertheless have an active integration/excision system and be capable of moving from one chromosomal position to another. We show here that the M. tuberculosis H37Rv phiRv1 element does indeed encode an active site-specific recombination system in which an integrase of the serine recombinase family (Rv1586c) catalyses integration and excision and a small, basic phiRv1-encoded protein (Rv1584c) controls the directionality of re-combination. Integration-proficient plasmid vectors derived from phiRv1 efficiently transform BCG, can utilize four of the seven REP13E12 sites present in BCG as attachment sites, and can occupy more than one site simultaneously.     

Tracing the origin of functional and conserved domains in the human proteome: implications for protein evolution at the modular level

Background

The highly homologous PE_PGRS (Proline-glutamic acid_polymorphic GC-rich repetitive sequence) genes are members of the PE multigene family which is found only in mycobacteria. PE genes are particularly abundant within the genomes of pathogenic mycobacteria where they seem to have expanded as a result of gene duplication events. PE_PGRS genes are characterized by their high GC content and extensive repetitive sequences, making them prone to recombination events and genetic variability.

Results

Comparative sequence analysis of Mycobacterium tuberculosis genes PE_PGRS17 (Rv0978c) and PE_PGRS18 (Rv0980c) revealed a striking genetic variation associated with this typical tandem duplicate. In comparison to the M. tuberculosis reference strain H37Rv, the variation (named the 12/40 polymorphism) consists of an in-frame 12-bp insertion invariably accompanied by a set of 40 single nucleotide polymorphisms (SNPs) that occurs either in PE_PGRS17 or in both genes. Sequence analysis of the paralogous genes in a representative set of worldwide distributed tubercle bacilli isolates revealed data which supported previously proposed evolutionary scenarios for the M. tuberculosis complex (MTBC) and confirmed the very ancient origin of " M. canettii " and other smooth tubercle bacilli. Strikingly, the identified polymorphism appears to be coincident with the emergence of the post-bottleneck successful clone from which the MTBC expanded. Furthermore, the findings provide direct and clear evidence for the natural occurrence of gene conversion in mycobacteria, which appears to be restricted to modern M. tuberculosis strains.

Conclusion

This study provides a new perspective to explore the molecular events that accompanied the evolution, clonal expansion, and recent diversification of tubercle bacilli.     

Functional analysis of an intergenic non-coding sequence within <Emphasis Type="Italic">mce1</Emphasis> operon of <Emphasis Type="Italic">M.tuberculosis</Emphasis>

Background  

The mce operons play an important role in the entry of M. tuberculosis into macrophages and non-phagocytic cells. Their non-redundant function as well as complex regulation is implied by the phenotype of mce mutants. Recently, mce1 operon was found to extend over 13 genes, fadD5 (Rv0166) being the first gene of the operon. The presence of a non-coding sequence of 200 base pairs between Rv0166 and Rv0167 is peculiar to mce1 among the four mce operons of M.tuberculosis. We have examined the function of this region.     

Evidence for a rapid rate of molecular evolution at the hypervariable and immunogenic Mycobacterium tuberculosis PPE38 gene region

Background  

PPE38 (Rv 2352c) is a member of the large PPE gene family of Mycobacterium tuberculosis and related mycobacteria. The function of PPE proteins is unknown but evidence suggests that many are cell-surface associated and recognised by the host immune system. Previous studies targeting other PPE gene members suggest that some display high levels of polymorphism and it is thought that this might represent a means of providing antigenic variation. We have analysed the genetic variability of the PPE38 genomic region on a cohort of M. tuberculosis clinical isolates representing all of the major phylogenetic lineages, along with the ancestral M. tuberculosis complex (MTBC) member M. canettii, and supplemented this with analysis of publicly available whole genome sequences representing additional M. tuberculosis clinical isolates, other MTBC members and non tuberculous mycobacteria (NTM). Where possible we have extended this analysis to include the adjacent plcABC and PPE39/40 genomic regions.     

High Sequence Variability of the ppE18 Gene of Clinical Mycobacterium tuberculosis Complex Strains Potentially Impacts Effectivity of Vaccine Candidate M72/AS01E

The development of an effective vaccine is urgently needed to fight tuberculosis (TB) which is still the leading cause of death from a single infectious agent worldwide. One of the promising vaccine candidates M72/AS01E consists of two proteins subunits PepA and PPE18 coded by Rv0125 and Rv1196. However, preliminary data indicate a high level of sequence variability among clinical Mycobacterium tuberculosis complex (MTBC) strains that might have an impact on the vaccine efficacy. To further investigate this finding, we determined ppE18 sequence variability in a well-characterized reference collection of 71 MTBC strains from 23 phylogenetic lineages representing the global MTBC diversity. In total, 100 sequence variations consisting of 96 single nucleotide polymorphisms (SNPs), three insertions and one deletion were detected resulting in 141 variable positions distributed over the entire gene. The majority of SNPs detected were non-synonymous (n = 68 vs. n = 28 synonymous). Strains from animal adapted lineages, e.g., M. bovis, showed a significant higher diversity than the human pathogens such as M. tuberculosis Haarlem. SNP patterns specific for different lineages as well as for deeper branches in the phylogeny could be identified. The results of our study demonstrate a high variability of the ppE18 gene even in the N-terminal domains that is normally highly conserved in ppe genes. As the N-terminal region interacts with TLR2 receptor inducing a protective anti-inflammatory immune response, genetic heterogeneity has a potential impact on the vaccine efficiency, however, this has to be investigated in future studies.     

The Mycobacterium DosR regulon structure and diversity revealed by comparative genomic analysis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), which claims approximately two million people annually, remains a global health concern. The non‐replicating or dormancy like state of this pathogen which is impervious to anti‐tuberculosis drugs is widely recognized as the culprit for this scenario. The dormancy survival regulator (DosR) regulon, composed of 48 co‐regulated genes, is held as essential for Mtb persistence. The DosR regulon is regulated by a two‐component regulatory system consisting of two sensor kinases—DosS (Rv3132c) and DosT (Rv2027c), and a response regulator DosR (Rv3133c). The underlying regulatory mechanism of DosR regulon expression is very complex. Many factors are involved, particularly the oxygen tension. The DosR regulon enables the pathogen to persist during lengthy hypoxia. Comparative genomic analysis demonstrated that the DosR regulon is widely distributed among the mycobacterial genomes, ranging from the pathogenic strains to the environmental strains. In‐depth studies on the DosR response should provide insights into its role in TB latency in vivo and shape new measures to combat this exceeding recalcitrant pathogen. J. Cell. Biochem. 114: 1–6, 2012. © 2012 Wiley Periodicals, Inc.     

Characterization of a chromosomal toxin-antitoxin, Rv1102c-Rv1103c system in Mycobacterium tuberculosis

Toxin-antitoxin systems, ubiquitous in prokaryotic genomes, have been proposed to play an important role in several stress responses. While Mycobacterium tuberculosis contains more than 80 putative TA loci, the roles they play in this pathogen are yet to be studied. Here, we characterize a chromosomal Rv1102c-Rv1103c TA system in M. tuberculosis. We found that the Rv1102c toxin interacts with the Rv1103c antitoxin in a pull-down assay and the yeast two-hybrid system. Rv1102c cleaved the era mRNA in Escherichia coli, and cleavage was inhibited by co-expression of Rv1103c. Heterologous expression of Rv1102c led to growth arrest in E. coli, which was fully recovered only when Rv1103c was co-expressed in cis with Rv1102c, suggesting that the production and assembly of Rv1102c and Rv1103c are tightly linked. Our additional results indicate that translational coupling of the Rv1102c and Rv1103c genes is important for Rv1102c-Rv1103c binding. Finally, we discovered that the expression of Rv1102c induced growth arrest and increased the level of persister cells in Mycobacterium smegmatis. These results suggest that the Rv1102c-Rv1103c TA system could play a role in M. tuberculosis pathogenesis via generating bacilli that survive in the face of multidrug therapy.     

Mycobacterium tuberculosis Rv3034c regulates mTORC1 and PPAR‐γ dependant pexophagy mechanism to control redox levels in macrophages

Mycobacterium tuberculosis survives inside the macrophages by employing several host immune evasion strategies. Here, we reported a novel mechanism in which M. tuberculosis acetyltransferase, encoded by Rv3034c, induces peroxisome homeostasis to regulate host oxidative stress levels to facilitate intracellular mycobacterial infection. Presence of M. tuberculosis Rv3034c induces the expression of peroxisome biogenesis and proliferation factors such as Pex3, Pex5, Pex19, Pex11b, Fis‐1 and DLP‐1; while depletion of Rv3034c decreased the expression of these molecules, thereby selective degradation of peroxisomes via pexophagy. Further studies revealed that M. tuberculosis Rv3034c inhibit induction of pexophagy mechanism by down‐regulating the expression of pexophagy associated proteins (p‐AMPKα, p‐ULK‐1, Atg5, Atg7, Beclin‐1, LC3‐II, TFEB and Keap‐1) and adaptor molecules (NBR1 and p62). Inhibition was found to be dependent on the phosphorylation of mTORC1 and activation of peroxisome proliferator activated receptor‐γ. In order to maintain intracellular homeostasis during oxidative stress, M. tuberculosis Rv3034c was found to induce degradation of dysfunctional and damaged peroxisomes through activation of Pex14 in infected macrophages. In conclusion, this is the first report which demonstrated that M. tuberculosis acetyltransferase regulate peroxisome homeostasis in response to intracellular redox levels to favour mycobacterial infection in macrophage.     

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

【目的】探讨异烟肼(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耐药有关。     

Mammalian heterotrimeric G-protein-like proteins in mycobacteria: implications for cell signalling and survival in eukaryotic host cells

Mammalian heterotrimeric GTP-binding proteins (G proteins) are involved in transmembrane signalling that couples a number of receptors to effectors mediating various physiological processes in mammalian cells. We demonstrate that bacterial proteins such as a Ras-like protein from Pseudomonas aeruginosa or a 65 kDa protein from Mycobacterium smegmatis can form complexes with human or yeast nucleoside diphosphate kinase (Ndk) to modulate their nucleoside triphosphate synthesizing specificity to GTP or UTP. In addition, we demonstrate that bacteria such as M. smegmatis or Mycobacterium tuberculosis harbour proteins that cross react with antibodies against the α-, β- or the γ-subunits of heterotrimeric G proteins. Such antibodies also alter the GTP synthesizing ability of specific membrane fractions isolated from glycerol gradients of such cells, suggesting that a membrane-associated Ndk–G-protein homologue complex is responsible for part of GTP synthesis in these bacteria. Indeed, purified Ndk from human erythrocytes and M. tuberculosis showed extensive complex formation with the purified mammalian α and β G-protein subunits and allowed specific GTP synthesis, suggesting that such complexes may participate in transmembrane signalling in the eukaryotic host. We have purified the α-, β- and γ-subunit homologues from M. tuberculosis and we present their internal amino acid sequences as well as their putative homologies with mammalian subunits and the localization of their genes on the M. tuberculosis genome. Using oligonucleotide probes from the conserved regions of the α- and γ-subunit of M. tuberculosis G-protein homologue, we demonstrate hybridization of these probes with the genomic digest of M. tuberculosis H37Rv but not with that of M. smegmatis, suggesting that M. smegmatis might lack the genes present in M. tuberculosis H37Rv. Interestingly, the avirulent strain H37Ra showed weak hybridization with these two probes, suggesting that these genes might have been deleted in the avirulent strain or are present in limited copy numbers as opposed to those in the virulent strain H37Rv.     

Comparative Whole-Genome Analysis of Clinical Isolates Reveals Characteristic Architecture of Mycobacterium tuberculosis Pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.     

Characterization of a bifunctional β-lactamase/ribonuclease and its interaction with a chaperone-like protein in the pathogen <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> H37Rv

Most mycobacteria appear to be naturally resistant to β-lactam antibiotics such as penicillin. However, very few β-lactamases and their regulation have been clearly characterized in Mycobacterium tuberculosis H37Rv. In this study, a unique bifunctional protein, Rv2752c, from M. tuberculosis showed both β-lactamase and RNase activities. Two residues, D184 and H397, appear to be involved in Zn²⁺-binding and are essential for the dual functions. Both activities are lost upon deletion of the C-terminal 100 a.a. long Rv2752c tail, which contains an additional loop when compared with the RNase J of Bacillus subtilis. A chaperone-like protein, Rv2373c, physically interacted with Rv2752c and inhibited both activities. This is the first report of characterization of a bifunctional β-lactamase and its regulation in mycobacteria. These data offered important clues for further investigation of the structure and function of microbial β-lactamases. Increased understanding of this protein will provide further insights into the mechanism of microbial drug resistance.     

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.     

Association of the Rv0679c protein with lipids and carbohydrates in <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>/<Emphasis Type="Italic">Mycobacterium bovis</Emphasis> BCG

The Rv0679c gene in Mycobacterium tuberculosis H37Rv encodes a protein with a predicted molecular mass of 16,586 Da consisting of 165 amino acids which contains a putative N-terminal signal sequence and a consensus lipoprotein-processing motif. Globomycin treatment, Triton X-114 separation and mass spectrometry analyses clarified a property of the Rv0679c protein as a lipoprotein. In addition, trifluoromethanesulphonic acid treatment of the lysate revealed an association of the recombinant Rv0679c protein with carbohydrates. The Rv0679c protein homolog of Mycobacterium bovis BCG was also expressed as the protein associated with lipids and carbohydrates. In Western blot analysis, each of the protein homolog and Lipoarabinomannan (LAM) was detected as a similar pattern by anti-Rv0679c and anti-LAM antibodies, respectively. Interestingly, the Rv0679c protein was detected in commercially available LAM purified from M. tuberculosis. Inhibition assay of LAM synthesis in M. bovis BCG by ethambutol showed an altered migration pattern of the Rv0679c protein to low molecular mass similar to that of LAM. The results suggest that the Rv0679c protein exists as a tight complex with LAM in M. tuberculosis/M. bovis BCG.     

One‐plasmid tunable coexpression for mycobacterial protein–protein interaction studies

A single plasmid that allows controlled coexpression has been developed for use in mycobacteria. The tetracycline inducible promoter, PtetO, was used to provide tetracycline‐dependent induction of one gene, while the Psmyc, Pimyc, or Phsp promoters were used to provide three different levels of constitutive expression of a second gene. The functions of these four individual promoters were established using green fluorescent protein (GFP) and a newly identified red fluorescence inducible protein from Geobacillus sterothermophilus strain G1.13 (RFIP) as reporters. The tandem use of GFP and RFIP as reporter genes allowed optimization of the tunable coexpression in Mycobacterium smegmatis; either time at a fixed inducer concentration or changes in inducer concentration could be used to control the protein:protein ratio. This single vector system was used to coexpress the two‐protein Mycobacterium tuberculosis stearoyl‐CoA Δ⁹ desaturase complex (integral membrane desaturase Rv3229c and NADPH oxidoreductase Rv3230c) in M. smegmatis. The catalytic activity was found to increase in a manner corresponding to increasing the level of Rv3230c relative to a fixed level of Rv3229c. This system, which can yield finely tuned coexpression of the fatty acid desaturase complex in mycobacteria, may be useful for study of other multicomponent complexes. Furthermore, the tunable coexpression strategy used herein should also be applicable in other species with minor modifications.     

15N, 13C and 1H resonance assignments and secondary structure determination of the Mycobacterium tuberculosis Rv0287–Rv0288 protein complex

Comprehensive sequence specific ¹H, ¹⁵N, and ¹³C resonance assignments are reported for the Mycobacterium tuberculosis Rv0287–Rv0288 protein complex. Analysis of the chemical shift data obtained indicates that each protein in the complex contains two relatively long helical regions joined by an irregular loop.