A predicted operon map for Mycobacterium tuberculosis

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Global discovery the PstP interactions using Mtb proteome microarray and revealing novel connections with EthR

Mycobacterium tuberculosis (Mtb) serine/threonine protein phosphatase PstP plays an important role in regulating Mtb cell division and growth by reversible phosphorylation signaling. However, the substrates of Mtb with which the PstP interacts, and the underlying molecular mechanisms are still largely unknown. In this study, we performed an Mtb proteome microarray to globally identify the PstP bindings. In this way, we discovered 78 interactors between PstP and Mtb proteins, and found a novel connections with EthR. The interaction between PstP and EthR has been validated by Bio-Layer interferometry and Yeast-two-hybrid. And functional studies showed that PstP significantly enhances the binding between EthR and related DNA domain through its interaction with EthR. Phenotypically, overexpression of PstP promoted the resistance of Mycobacterium smegmatis with the antibiotic of ethionamide. Overall, we hopefully wish that the PstP interactors identified in this study will serve as a useful resource for further systematic studies of the roles that PstP plays in the regulation of Mtb dephosphorylation. SIGNIFICANCE: Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, which is responsible of ~1.5 million death per year. Understanding the knowledge about the basic biological regulation pathways in Mtb is an effective approach to discover the novel drug targets for cure TB. PstP is a serine/threonine protein phosphatase in Mtb, and plays important roles in regulating Mtb cell division and growth by reversible phosphorylation signaling. In this study, we identified 78 PstP interacting Mtb proteins using Mtb proteome microarray, which could preliminarily explain the roles of PstP played in Mtb. Moreover, functional analysis showed that a novel transcription factor EthR had been found regulated by PstP through binding, which could enhance the resistance to the antibiotic ETH. Overall, this network constructed with PstP-Mtb proteins could serve as a valuable resource for studying Mtb growth.

Mycobacterial phosphatase PstP regulates global serine threonine phosphorylation and cell division

Protein phosphatase PstP is conserved throughout the Actinobacteria in a genetic locus related to cell wall synthesis and cell division. In many Actinobacteria it is the sole annotated serine threonine protein phosphatase to counter the activity of multiple serine threonine protein kinases. We used transcriptional knockdown, electron microscopy and comparative phosphoproteomics to investigate the putative dual functions of PstP as a specific regulator of cell division and as a global regulator of protein phosphorylation. Comparative phosphoproteomics in the early stages of PstP depletion showed hyperphosphorylation of protein kinases and their substrates, confirming PstP as a negative regulator of kinase activity and global serine and threonine phosphorylation. Analysis of the 838 phosphorylation sites that changed significantly, suggested that PstP may regulate diverse phosphoproteins, preferentially at phosphothreonine near acidic residues, near the protein termini, and within membrane associated proteins. Increased phosphorylation of the activation loop of protein kinase B (PknB) and of the essential PknB substrate CwlM offer possible explanations for the requirement for pstP for growth and for cell wall defects when PstP was depleted.

Role of two-component regulatory systems in intracellular survival of Mycobacterium tuberculosis

The typical two-component regulatory systems (TCSs), consisting of response regulator and histidine kinase, play a central role in survival of pathogenic bacteria under stress conditions such as nutrient starvation, hypoxia, and nitrosative stress. A total of 11 complete paired two-component regulatory systems have been found in Mycobacterium tuberculosis, including a few isolated kinase and regulatory genes. Increasing evidence has shown that TCSs are closely associated with multiple physiological process like intracellular persistence, pathogenicity, and metabolism. This review gives the two-component signal transduction systems in M. tuberculosis and their signal transduction roles in adaption to the environment.

Rv1460, a SufR homologue, is a repressor of the suf operon in Mycobacterium tuberculosis

Iron–sulphur (Fe-S) clusters are ubiquitous co-factors which require multi-protein systems for their synthesis. In Mycobacterium tuberculosis, the Rv1460-Rv1461-Rv1462-Rv1463-csd-Rv1465-Rv1466 operon (suf operon) encodes the primary Fe-S cluster biogenesis system. The first gene in this operon, Rv1460, shares homology with the cyanobacterial SufR, which functions as a transcriptional repressor of the sufBCDS operon. Rv1460’s function in M. tuberculosis has however not been determined. In this study, we demonstrate that M. tuberculosis mutants lacking a functional Rv1460 protein are impaired for growth under standard culture conditions. Elevated expression of Rv1460 and Rv1461 was observed in the mutant, implicating Rv1460 in the regulation of the suf operon. Binding of an Fe-S cluster to purified recombinant Rv1460 was confirmed by UV-visible spectroscopy and circular dichroism. Furthermore, three conserved cysteine residues, C203, C216 and C244, proposed to provide ligands for the coordination of an Fe-S cluster, were shown to be required for the function of Rv1460 in M. tuberculosis. Rv1460 therefore seems to be functionally analogous to cyanobacterial SufR.

Computational operon prediction in whole-genomes and metagenomes

Microbial diversity in unique environmental settings enables abrupt responses catalysed by altering the gene regulation and formation of gene clusters called operons. Operons increases bacterial adaptability, which in turn increases their survival. This review article presents the emergence of computational operon prediction methods for whole microbial genomes and metagenomes, and discusses their strengths and limitations. Most of the whole-genome operon prediction methods struggle to generalize on unrelated genomes. The applicability of universal whole-genome operon prediction methods to metagenomic data is an interesting yet less investigated question. We have evaluated the potential of various operon prediction features for genomic and metagenomic data. Most of operon prediction methods with high accuracy have been compiled into databases. Despite of the high predictive performance, the data among many databases are not completely consistent for similar species. We performed a correlation analysis between the computationally predicted operon databases and experimentally validated data for Escherichia coli, Bacillus subtilis and Mycobacterium tuberculosis. Operon prediction for most of the less characterized microbes cannot be verified due to absence of experimentally validated operons. The generation of validated information for other microbes would test the authenticity of operon databases for other less annotated microbes as well. Advances in sequencing technologies and development of better analysis methods will help researchers to overcome the technological hurdles (such as long sequencing reads and improved contig size) and further improve operon predictions and better utilize operonic information.

Structural and genetic analysis of START superfamily protein MSMEG_0129 from Mycobacterium smegmatis

Mycobacterium tuberculosis is a notorious pathogen that continues to threaten human health. Rv0164, an antigen of both T- and B cells conserved across mycobacteria, and MSMEG_0129, its close homolog in Mycobacterium smegmatis, are predicted members of the START domain superfamily, but their molecular function is unknown. Here, gene knockout studies demonstrate MSMEG_0129 is essential for bacterial growth, suggesting Rv0164 may be a potential drug target. The MSMEG_0129 crystal structure determined at 1.95 Å reveals a fold similar to that in polyketide aromatase/cyclases ZhuI and TcmN from Streptomyces sp. Structural comparisons and docking simulations, however, infer that MSMEG_0129 and Rv0164 are unlikely to catalyze polyketide aromatization/cyclization, but probably play an irreplaceable role during mycobacterial growth, for example, in lipid transfer during cell envelope synthesis.

The transcriptome of Mycobacterium tuberculosis in a lipid-rich dormancy model through RNAseq analysis

Tuberculosis (TB) is currently the number one killer among infectious diseases worldwide. Lipids are abundant molecules during the infectious cycle of Mycobacterium tuberculosis (Mtb) and studies better mimicking its actual metabolic state during pathogenesis are needed. Though most studies have focused on the mycobacterial lipid metabolism under standard culture conditions, little is known about the transcriptome of Mtb in a lipid environment. Here we determined the transcriptome of Mtb H37Rv in a lipid-rich environment (cholesterol and fatty acid) under aerobic and hypoxic conditions, using RNAseq. Lipids significantly induced the expression of 368 genes. A main core lipid response was observed involving efflux systems, iron caption and sulfur reduction. In co-expression with ncRNAs and other genes discussed below, may act coordinately to prepare the machinery conferring drug tolerance and increasing a persistent population. Our findings could be useful to tag relevant pathways for the development of new drugs, vaccines and new strategies to control TB.

Probing Ligand Exchange in the P450 Enzyme CYP121 from Mycobacterium tuberculosis: Dynamic Equilibrium of the Distal Heme Ligand as a Function of pH and Temperature

CYP121 is a cytochrome P450 enzyme from Mycobacterium tuberculosis that catalyzes the formation of a C-C bond between the aromatic groups of its cyclodityrosine substrate (cYY). The crystal structure of CYP121 in complex with cYY reveals that the solvent-derived ligand remains bound to the ferric ion in the enzyme-substrate complex. Whereas in the generally accepted P450 mechanism, binding of the primary substrate in the active-site triggers the release of the solvent-derived ligand, priming the metal center for reduction and subsequent O2 binding. Here we employed sodium cyanide to probe the metal-ligand exchange of the enzyme and the enzyme-substrate complex. The cyano adducts were characterized by UV-vis, EPR, and ENDOR spectroscopies and X-ray crystallography. A 100-fold increase in the affinity of cyanide binding to the enzyme-substrate complex over the ligand-free enzyme was observed. The crystal structure of the [CYP121(cYY)CN] ternary complex showed a rearrangement of the substrate in the active-site, when compared to the structure of the binary [CYP121(cYY)] complex. Transient kinetic studies showed that cYY binding resulted in a lower second-order rate constant (kon (CN)) but a much more stable cyanide adduct with 3 orders of magnitude slower koff (CN) rate. A dynamic equilibrium between multiple high- and low-spin species for both the enzyme and enzyme-substrate complex was also observed, which is sensitive to changes in both pH and temperature. Our data reveal the chemical and physical properties of the solvent-derived ligand of the enzyme, which will help to understand the initial steps of the catalytic mechanism.

Construction of Mycobacterium tuberculosis cdd knockout and evaluation of invasion and growth in macrophages

Cytidine deaminase (MtCDA), encoded by cdd gene (Rv3315c), is the only enzyme identified in nucleotide biosynthesis pathway of Mycobacterium tuberculosis that is able to recycle cytidine and deoxycytidine. An M. tuberculosis knockout strain for cdd gene was obtained by allelic replacement. Evaluation of mRNA expression validated cdd deletion and showed the absence of polar effect. MudPIT LC-MS/MS data indicated thymidine phosphorylase expression was decreased in knockout and complemented strains. The cdd disruption does not affect M. tuberculosis growth both in Mid- dlebrook 7H9 and in RAW 264.7 cells, which indicates that cdd is not important for macrophage invasion and virulence.

A versatile vector for mycobacterial protein production with a functional minimized acetamidase regulon

Recombinant protein expression is a prerequisite for diverse investigations of proteins at the molecular level. For targets from Mycobacterium tuberculosis it is favorable to use M. smegmatis as an expression host, a species from the same genus. In the respective shuttle vectors, target gene expression is controlled by the complex tetra-cistronic acetamidase regulon. As a result, the size of those vectors is large, rendering them of limited use, especially when the target proteins are expressed from multi-cistronic operons. Therefore, in the current work we present a versatile new expression vector in which the acetamidase regulon has been minimized by deleting the two genes amiD and amiS. We assessed the functional properties of the resulting vector pMyCA and compared it with those of the existing vector pMyNT that contains the full-length acetamidase regulon. We analyzed the growth features and protein expression patterns of M. smegmatis cultures transformed with both vectors. In addition, we created mCherry expression constructs to spectroscopically monitor the expression properties of both vectors. Our experiments showed that the minimized vector exhibited several advantages over the pMyNT vector. First, the overall yield of expressed protein is higher due to the higher yield of bacterial mass. Second, the heterologous expression was regulated more tightly, offering an expression tool for diverse target proteins. Third, it is suitable for large multi-protein complexes that are expressed from multi-cistronic operons. Additionally, our results propose a new understanding of the regulation mechanism of the acetamidase regulon with the potential to construct more optimized vectors in the future.

Cross-linking of dicyclotyrosine by the cytochrome P450 enzyme CYP121 from Mycobacterium tuberculosis proceeds through a catalytic shunt pathway

CYP121, the cytochrome P450 enzyme in Mycobacterium tuberculosis that catalyzes a single intramolecular C-C cross-linking reaction in the biosynthesis of mycocyclosin, is crucial for the viability of this pathogen. This C-C coupling reaction represents an expansion of the activities carried out by P450 enzymes distinct from oxygen insertion. Although the traditional mechanism for P450 enzymes has been well studied, it is unclear whether CYP121 follows the general P450 mechanism or uses a different catalytic strategy for generating an iron-bound oxidant. To gain mechanistic insight into the CYP121-catalyzed reaction, we tested the peroxide shunt pathway by using rapid kinetic techniques to monitor the enzyme activity with its substrate dicyclotyrosine (cYY) and observed the formation of the cross-linked product mycocyclosin by LC-MS. In stopped-flow experiments, we observed that cYY binding to CYP121 proceeds in a two-step process, and EPR spectroscopy indicates that the binding induces active site reorganization and uniformity. Using rapid freeze-quenching EPR, we observed the formation of a high-spin intermediate upon the addition of peracetic acid to the enzyme-substrate complex. This intermediate exhibits a high-spin (S = 5/2) signal with g values of 2.00, 5.77, and 6.87. Likewise, iodosylbenzene could also produce mycocyclosin, implicating compound I as the initial oxidizing species. Moreover, we also demonstrated that CYP121 performs a standard peroxidase type of reaction by observing substrate-based radicals. On the basis of these results, we propose plausible free radical-based mechanisms for the C-C bond coupling reaction.

The alr-groEL1 operon in Mycobacterium tuberculosis: an interplay of multiple regulatory elements

Threonylcarbamoyladenosine is a universally conserved essential modification of tRNA that ensures translational fidelity in cellular milieu. TsaD, TsaB and TsaE are identified as tRNA-A₃₇-threonylcarbamoyl (t⁶A)-transferase enzymes that have been reconstituted in vitro, in few bacteria recently. However, transcriptional organization and regulation of these genes are not known in any of these organisms. This study describes the intricate architecture of a complex multicistronic alr-groEL1 operon, harboring essential genes, namely tsaD, tsaB, tsaE, groES, groEL1, and alr (required for cell wall synthesis), and rimI encoding an N-α- acetyltransferase in Mycobacterium tuberculosis. Using northern blotting, RT-PCR and in vivo fluorescence assays, genes alr to groEL1 were found to constitute an ~6.3 kb heptacistronic operon with multiple internal promoters and an I-shaped intrinsic hairpin-like cis-regulatory element. A strong promoter PtsaD within the coding sequence of rimI gene is identified in M. tuberculosis, in addition. The study further proposes an amendment in the known bicistronic groESL1 operon annotation by providing evidence that groESL1 is co-transcribed as sub-operon of alr-groEL1 operon. The architecture of alr-groEL1 operon, conservation of the genetic context and a mosaic transcriptional profile displayed under various stress conditions convincingly suggest the involvement of this operon in stress adaptation in M. tuberculosis.

REMap: Operon map of M. tuberculosis based on RNA sequence data

A map of the transcriptional organization of genes of an organism is a basic tool that is necessary to understand and facilitate a more accurate genetic manipulation of the organism. Operon maps are largely generated by computational prediction programs that rely on gene conservation and genome architecture and may not be physiologically relevant. With the widespread use of RNA sequencing (RNAseq), the prediction of operons based on actual transcriptome sequencing rather than computational genomics alone is much needed. Here, we report a validated operon map of Mycobacterium tuberculosis, developed using RNAseq data from both the exponential and stationary phases of growth. At least 58.4% of M. tuberculosis genes are organized into 749 operons. Our prediction algorithm, REMap (RNA Expression Mapping of operons), considers the many cases of transcription coverage of intergenic regions, and avoids dependencies on functional annotation and arbitrary assumptions about gene structure. As a result, we demonstrate that REMap is able to more accurately predict operons, especially those that contain long intergenic regions or functionally unrelated genes, than previous operon prediction programs. The REMap algorithm is publicly available as a user-friendly tool that can be readily modified to predict operons in other bacteria.

Integration of Metabolic Modeling with Gene Co-expression Reveals Transcriptionally Programmed Reactions Explaining Robustness in Mycobacterium tuberculosis

Robustness of metabolic networks is accomplished by gene regulation, modularity, re-routing of metabolites and plasticity. Here, we probed robustness against perturbations of biochemical reactions of M. tuberculosis in the form of predicting compensatory trends. In order to investigate the transcriptional programming of genes associated with correlated fluxes, we integrated with gene co-expression network. Knock down of the reactions NADH2r and ATPS responsible for producing the hub metabolites, and Central carbon metabolism had the highest proportion of their associated genes under transcriptional co-expression with genes of their flux correlated reactions. Reciprocal gene expression correlations were observed among compensatory routes, fresh activation of alternative routes and in the multi-copy genes of Cysteine synthase and of Phosphate transporter. Knock down of 46 reactions caused the activation of Isocitrate lyase or Malate synthase or both reactions, which are central to the persistent state of M. tuberculosis. A total of 30 new freshly activated routes including Cytochrome c oxidase, Lactate dehydrogenase, and Glycine cleavage system were predicted, which could be responsible for switching into dormant or persistent state. Thus, our integrated approach of exploring transcriptional programming of flux correlated reactions has the potential to unravel features of system architecture conferring robustness.

Mycobacterium tuberculosis folate metabolism and the mechanistic basis for para-aminosalicylic acid susceptibility and resistance

para-Aminosalicylic acid (PAS) entered clinical use in 1946 as the second exclusive drug for the treatment of tuberculosis (TB). While PAS was initially a first-line TB drug, the introduction of more potent antitubercular agents relegated PAS to the second-line tier of agents used for the treatment of drug-resistant Mycobacterium tuberculosis infections. Despite the long history of PAS usage, an understanding of the molecular and biochemical mechanisms governing the susceptibility and resistance of M. tuberculosis to this drug has lagged behind that of most other TB drugs. Herein, we discuss previous studies that demonstrate PAS-mediated disruption of iron acquisition, as well as recent genetic, biochemical, and metabolomic studies that have revealed that PAS is a prodrug that ultimately corrupts one-carbon metabolism through inhibition of the formation of reduced folate species. We also discuss findings from laboratory and clinical isolates that link alterations in folate metabolism to PAS resistance. These advancements in our understanding of the basis of the susceptibility and resistance of M. tuberculosis to PAS will enable the development of novel strategies to revitalize this and other antimicrobial agents for use in the global effort to eradicate TB.

The Psp system of Mycobacterium tuberculosis integrates envelope stress-sensing and envelope-preserving functions

The bacterial envelope integrates essential stress-sensing and adaptive functions; thus, envelope-preserving functions are important for survival. In Gram-negative bacteria, envelope integrity during stress is maintained by the multi-gene Psp response. Mycobacterium tuberculosis was thought to lack the Psp system since it encodes only pspA and no other psp ortholog. Intriguingly, pspA maps downstream from clgR, which encodes a transcription factor regulated by the MprAB-σ(E) envelope-stress-signaling system. clgR inactivation lowered ATP concentration during stress and protonophore treatment-induced clgR-pspA expression, suggesting that these genes express Psp-like functions. We identified a four-gene set - clgR, pspA (rv2744c), rv2743c, rv2742c - that is regulated by clgR and in turn regulates ClgR activity. Regulatory and protein-protein interactions within the set and a requirement of the four genes for functions associated with envelope integrity and surface-stress tolerance indicate that a Psp-like system has evolved in mycobacteria. Among Actinobacteria, the four-gene module occurred only in tuberculous mycobacteria and was required for intramacrophage growth, suggesting links between its function and mycobacterial virulence. Additionally, the four-gene module was required for MprAB-σ(E) stress-signaling activity. The positive feedback between envelope-stress-sensing and envelope-preserving functions allows sustained responses to multiple, envelope-perturbing signals during chronic infection, making the system uniquely suited to tuberculosis pathogenesis.

Differential expression of two members of Rv1922-LipD operon in Mycobacterium tuberculosis: Does rv1923 qualify for membership?

rv1922 and rv1923 (lipD) are members of Rv1922-LipD operon in the genome of Mycobacterium tuberculosis. rv1922 was expressed under aerobic and stress conditions, whereas rv1923 was not expressed in aerobically grown bacteria but expressed moderately under oxidative stress conditions. Different expression of both the operonic genes under normal and stress conditions posed apprehensions for the inclusion of rv1922 and rv1923 in the same operon. The results from this study indicated that although the genes were expressed in an operonic manner, there existed the possibility of differential regulation for rv1923 which has been supported by in silico analysis for the presence of putative internal regulatory sites in the operon.

Integration of heterogeneous molecular networks to unravel gene-regulation in Mycobacterium tuberculosis

BACKGROUND

Different methods have been developed to infer regulatory networks from heterogeneous omics datasets and to construct co-expression networks. Each algorithm produces different networks and efforts have been devoted to automatically integrate them into consensus sets. However each separate set has an intrinsic value that is diluted and partly lost when building a consensus network. Here we present a methodology to generate co-expression networks and, instead of a consensus network, we propose an integration framework where the different networks are kept and analysed with additional tools to efficiently combine the information extracted from each network.

RESULTS

We developed a workflow to efficiently analyse information generated by different inference and prediction methods. Our methodology relies on providing the user the means to simultaneously visualise and analyse the coexisting networks generated by different algorithms, heterogeneous datasets, and a suite of analysis tools. As a show case, we have analysed the gene co-expression networks of Mycobacterium tuberculosis generated using over 600 expression experiments. Regarding DNA damage repair, we identified SigC as a key control element, 12 new targets for LexA, an updated LexA binding motif, and a potential mismatch repair system. We expanded the DevR regulon with 27 genes while identifying 9 targets wrongly assigned to this regulon. We discovered 10 new genes linked to zinc uptake and a new regulatory mechanism for ZuR. The use of co-expression networks to perform system level analysis allows the development of custom made methodologies. As show cases we implemented a pipeline to integrate ChIP-seq data and another method to uncover multiple regulatory layers.

CONCLUSIONS

Our workflow is based on representing the multiple types of information as network representations and presenting these networks in a synchronous framework that allows their simultaneous visualization while keeping specific associations from the different networks. By simultaneously exploring these networks and metadata, we gained insights into regulatory mechanisms in M. tuberculosis that could not be obtained through the separate analysis of each data type.

Specialized transduction designed for precise high-throughput unmarked deletions in Mycobacterium tuberculosis

Specialized transduction has proven to be useful for generating deletion mutants in most mycobacteria, including virulent Mycobacterium tuberculosis. We have improved this system by developing (i) a single-step strategy for the construction of allelic exchange substrates (AES), (ii) a temperature-sensitive shuttle phasmid with a greater cloning capacity than phAE87, and (iii) bacteriophage-mediated transient expression of site-specific recombinase to precisely excise antibiotic markers. The methods ameliorate rate-limiting steps in strain construction in these difficult-to-manipulate bacteria. The new methods for strain construction were demonstrated to generalize to all classes of genes and chromosomal loci by generating more than 100 targeted single- or multiple-deletion substitutions. These improved methods pave the way for the generation of a complete ordered library of M. tuberculosis null strains, where each strain is deleted for a single defined open reading frame in M. tuberculosis.

This work reports major advances in the methods of genetics applicable to all mycobacteria, including but not limited to virulent M. tuberculosis, which would facilitate comparative genomics to identify drug targets, genetic validation of proposed pathways, and development of an effective vaccine. This study presents all the new methods developed and the improvements to existing methods in an integrated way. The work presented in this study could increase the pace of mycobacterial genetics significantly and will immediately be of wide use. These new methods are transformative and allow for the undertaking of construction of what has been one of the most fruitful resources in model systems: a comprehensive, ordered library set of the strains, each of which is deleted for a single defined open reading frame.

The PhoP-dependent ncRNA Mcr7 modulates the TAT secretion system in Mycobacterium tuberculosis

The PhoPR two-component system is essential for virulence in Mycobacterium tuberculosis where it controls expression of approximately 2% of the genes, including those for the ESX-1 secretion apparatus, a major virulence determinant. Mutations in phoP lead to compromised production of pathogen-specific cell wall components and attenuation both ex vivo and in vivo. Using antibodies against the native protein in ChIP-seq experiments (chromatin immunoprecipitation followed by high-throughput sequencing) we demonstrated that PhoP binds to at least 35 loci on the M. tuberculosis genome. The PhoP regulon comprises several transcriptional regulators as well as genes for polyketide synthases and PE/PPE proteins. Integration of ChIP-seq results with high-resolution transcriptomic analysis (RNA-seq) revealed that PhoP controls 30 genes directly, whilst regulatory cascades are responsible for signal amplification and downstream effects through proteins like EspR, which controls Esx1 function, via regulation of the espACD operon. The most prominent site of PhoP regulation was located in the intergenic region between rv2395 and PE_PGRS41, where the mcr7 gene codes for a small non-coding RNA (ncRNA). Northern blot experiments confirmed the absence of Mcr7 in an M. tuberculosis phoP mutant as well as low-level expression of the ncRNA in M. tuberculosis complex members other than M. tuberculosis. By means of genetic and proteomic analyses we demonstrated that Mcr7 modulates translation of the tatC mRNA thereby impacting the activity of the Twin Arginine Translocation (Tat) protein secretion apparatus. As a result, secretion of the immunodominant Ag85 complex and the beta-lactamase BlaC is affected, among others. Mcr7, the first ncRNA of M. tuberculosis whose function has been established, therefore represents a missing link between the PhoPR two-component system and the downstream functions necessary for successful infection of the host.

One of the best characterized two-component systems in Mycobacterium tuberculosis is represented by the PhoPR pair, with PhoR being the transmembrane sensor kinase and PhoP playing an essential part in controlling expression of virulence-associated genes, such as those encoding the ESX-1 secretion apparatus. Previous studies showed that mutations in phoP resulted in attenuation in the mouse model of infection, thus providing the basis for the development of a novel live attenuated Mycobacterium tuberculosis vaccine carrying a deletion in phoP which is today in clinical trials. To thoroughly investigate the role of PhoP in M. tuberculosis, we undertook a systems biology approach comprising ChIP-seq and RNA-seq technologies. We demonstrated binding of PhoP to at least 35 targets on the M. tuberculosis chromosome and direct impact on expression of 30 genes, while further amplification of the signal is provided by regulators acting downstream. The strongest binding site was located between rv2395 and PE_PGRS41, where transcription of the non-coding RNA (ncRNA) Mcr7 was demonstrated. Expression of Mcr7 was found to be restricted to M. tuberculosis species and totally silenced in a phoP mutant. Genetics and proteomics approaches proved that Mcr7 controls activity of the Twin Arginine (Tat) secretion system, thus modulating secretion of the immunodominant antigen Ag85 complex and the BlaC beta-lactamase.

WXG100 protein superfamily consists of three subfamilies and exhibits an α-helical C-terminal conserved residue pattern

Members of the WXG100 protein superfamily form homo- or heterodimeric complexes. The most studied proteins among them are the secreted T-cell antigens CFP-10 (10 kDa culture filtrate protein, EsxB) and ESAT-6 (6 kDa early secreted antigen target, EsxA) from Mycobacterium tuberculosis. They are encoded on an operon within a gene cluster, named as ESX-1, that encodes for the Type VII secretion system (T7SS). WXG100 proteins are secreted in a full-length form and it is known that they adopt a four-helix bundle structure. In the current work we discuss the evolutionary relationship between the homo- and heterodimeric WXG100 proteins, the basis of the oligomeric state and the key structural features of the conserved sequence pattern of WXG100 proteins. We performed an iterative bioinformatics analysis of the WXG100 protein superfamily and correlated this with the atomic structures of the representative WXG100 proteins. We find, firstly, that the WXG100 protein superfamily consists of three subfamilies: CFP-10-, ESAT-6- and sagEsxA-like proteins (EsxA proteins similar to that of Streptococcus agalactiae). Secondly, that the heterodimeric complexes probably evolved from a homodimeric precursor. Thirdly, that the genes of hetero-dimeric WXG100 proteins are always encoded in bi-cistronic operons and finally, by combining the sequence alignments with the X-ray data we identify a conserved C-terminal sequence pattern. The side chains of these conserved residues decorate the same side of the C-terminal α-helix and therefore form a distinct surface. Our results lead to a putatively extended T7SS secretion signal which combines two reported T7SS recognition characteristics: Firstly that the T7SS secretion signal is localized at the C-terminus of T7SS substrates and secondly that the conserved residues YxxxD/E are essential for T7SS activity. Furthermore, we propose that the specific α-helical surface formed by the conserved sequence pattern including YxxxD/E motif is a key component of T7SS-substrate recognition.

IpsA, a novel LacI-type regulator, is required for inositol-derived lipid formation in Corynebacteria and Mycobacteria

Background

The development of new drugs against tuberculosis and diphtheria is focused on disrupting the biogenesis of the cell wall, the unique architecture of which confers resistance against current therapies. The enzymatic pathways involved in the synthesis of the cell wall by these pathogens are well understood, but the underlying regulatory mechanisms are largely unknown.

Results

Here, we characterize IpsA, a LacI-type transcriptional regulator conserved among Mycobacteria and Corynebacteria that plays a role in the regulation of cell wall biogenesis. IpsA triggers myo-inositol formation by activating ino1, which encodes inositol phosphate synthase. An ipsA deletion mutant of Corynebacterium glutamicum cultured on glucose displayed significantly impaired growth and presented an elongated cell morphology. Further studies revealed the absence of inositol-derived lipids in the cell wall and a complete loss of mycothiol biosynthesis. The phenotype of the C. glutamicum ipsA deletion mutant was complemented to different extend by homologs from Corynebacterium diphtheriae (dip1969) and Mycobacterium tuberculosis (rv3575), indicating the conserved function of IpsA in the pathogenic species. Additional targets of IpsA with putative functions in cell wall biogenesis were identified and IpsA was shown to bind to a conserved palindromic motif within the corresponding promoter regions. Myo-inositol was identified as an effector of IpsA, causing the dissociation of the IpsA-DNA complex in vitro.

Conclusions

This characterization of IpsA function and of its regulon sheds light on the complex transcriptional control of cell wall biogenesis in the mycolata taxon and generates novel targets for drug development.

Revealing of Mycobacterium marinum transcriptome by RNA-seq

Transcriptome analysis has played an essential role for revealing gene expression and the complexity of regulations at transcriptional level. RNA-seq is a powerful tool for transcriptome profiling, which uses deep-sequencing technologies to directly determine the cDNA sequence. Here, we utilized RNA-seq to explore the transcriptome of Mycobacteriummarinum (M. marinum), which is a useful model to study the pathogenesis of Mycobacterium tuberculosis (Mtb). Two profiles of exponential and early stationary phase cultures were generated after a physical ribosome RNA removal step. We systematically described the transcriptome and analyzed the functions for the differentiated expressed genes between the two phases. Furthermore, we predicted 360 operons throughout the whole genome, and 13 out of 17 randomly selected operons were validated by qRT-PCR. In general, our study has primarily uncovered M. marinum transcriptome, which could help to gain a better understanding of the regulation system in Mtb that underlines disease pathogenesis.

Operon prediction using chaos embedded particle swarm optimization

Operons contain valuable information for drug design and determining protein functions. Genes within an operon are co-transcribed to a single-strand mRNA and must be coregulated. The identification of operons is, thus, critical for a detailed understanding of the gene regulations. However, currently used experimental methods for operon detection are generally difficult to implement and time consuming. In this paper, we propose a chaotic binary particle swarm optimization (CBPSO) to predict operons in bacterial genomes. The intergenic distance, participation in the same metabolic pathway and the cluster of orthologous groups (COG) properties of the Escherichia coli genome are used to design a fitness function. Furthermore, the Bacillus subtilis, Pseudomonas aeruginosa PA01, Staphylococcus aureus and Mycobacterium tuberculosis genomes are tested and evaluated for accuracy, sensitivity, and specificity. The computational results indicate that the proposed method works effectively in terms of enhancing the performance of the operon prediction. The proposed method also achieved a good balance between sensitivity and specificity when compared to methods from the literature.

Substrate and reaction specificity of Mycobacterium tuberculosis cytochrome P450 CYP121: insights from biochemical studies and crystal structures

Cytochrome P450 CYP121 is essential for the viability of Mycobacterium tuberculosis. Studies in vitro show that it can use the cyclodipeptide cyclo(l-Tyr-l-Tyr) (cYY) as a substrate. We report an investigation of the substrate and reaction specificities of CYP121 involving analysis of the interaction between CYP121 and 14 cYY analogues with various modifications of the side chains or the diketopiperazine (DKP) ring. Spectral titration experiments show that CYP121 significantly bound only cyclodipeptides with a conserved DKP ring carrying two aryl side chains in l-configuration. CYP121 did not efficiently or selectively transform any of the cYY analogues tested, indicating a high specificity for cYY. The molecular determinants of this specificity were inferred from both crystal structures of CYP121-analog complexes solved at high resolution and solution NMR spectroscopy of the analogues. Bound cYY or its analogues all displayed a similar set of contacts with CYP121 residues Asn(85), Phe(168), and Trp(182). The propensity of the cYY tyrosyl to point toward Arg(386) was dependent on the presence of the DKP ring that limits the conformational freedom of the ligand. The correct positioning of the hydroxyl of this tyrosyl was essential for conversion of cYY. Thus, the specificity of CYP121 results from both a restricted binding specificity and a fine-tuned P450 substrate relationship. These results document the catalytic mechanism of CYP121 and improve our understanding of its function in vivo. This work contributes to progress toward the design of inhibitors of this essential protein of M. tuberculosis that could be used for antituberculosis therapy.

Evidence classification of high-throughput protocols and confidence integration in RegulonDB

RegulonDB provides curated information on the transcriptional regulatory network of Escherichia coli and contains both experimental data and computationally predicted objects. To account for the heterogeneity of these data, we introduced in version 6.0, a two-tier rating system for the strength of evidence, classifying evidence as either ‘weak’ or ‘strong’ (Gama-Castro,S., Jimenez-Jacinto,V., Peralta-Gil,M. et al. RegulonDB (Version 6.0): gene regulation model of Escherichia Coli K-12 beyond transcription, active (experimental) annotated promoters and textpresso navigation. Nucleic Acids Res., 2008;36:D120–D124.). We now add to our classification scheme the classification of high-throughput evidence, including chromatin immunoprecipitation (ChIP) and RNA-seq technologies. To integrate these data into RegulonDB, we present two strategies for the evaluation of confidence, statistical validation and independent cross-validation. Statistical validation involves verification of ChIP data for transcription factor-binding sites, using tools for motif discovery and quality assessment of the discovered matrices. Independent cross-validation combines independent evidence with the intention to mutually exclude false positives. Both statistical validation and cross-validation allow to upgrade subsets of data that are supported by weak evidence to a higher confidence level. Likewise, cross-validation of strong confidence data extends our two-tier rating system to a three-tier system by introducing a third confidence score ‘confirmed’.

Database URL:http://regulondb.ccg.unam.mx/

Reannotation of translational start sites in the genome of Mycobacterium tuberculosis

Identification and correction of incorrect ORF start sites is important for a variety of experimental and analytical purposes, ranging from cloning to inference of operon structure. The genome of the H37Rv reference strain of Mycobacterium tuberculosis (Mtb) was originally annotated when it was first sequenced nearly 15 years ago. While this annotation has served the TB research community well as a standard of reference for over a decade, it has been demonstrated experimentally that the actual start sites for an estimated 5-10% of open reading frames differ from the annotation. In this paper, we present a comprehensive bioinformatic analysis of all 3989 ORFs (open reading frames) in the M. tuberculosis H37Rv genome. Our method combines information from comparative analysis (alignment to start sites of orthologs in other Actinobacteria), sequence conservation, "protein likeness", putative ribosome binding sites, and other data to identify translational start sites. The features are combined in a linear model that is trained on dataset of known start sites verified by mass spectrometry, with a cross-validated accuracy of 94%. The method can be viewed as an augmentation of Hidden Markov Model-based tools such as Glimmer and GeneMark by incorporating more information than just the raw genomic sequence to decide which position is the legitimate translational start site for each ORF. Using this analysis, we identify 269 genes that most likely need to be re-annotated, and identify the best alterative translational start site for each. These revised ORF definitions could be used in the reannotation of the H37Rv genome, as well as to prioritize genes for experimental start-site validation.

High-resolution transcriptome and genome-wide dynamics of RNA polymerase and NusA in Mycobacterium tuberculosis

To construct a regulatory map of the genome of the human pathogen, Mycobacterium tuberculosis, we applied two complementary high-resolution approaches: strand-specific RNA-seq, to survey the global transcriptome, and ChIP-seq, to monitor the genome-wide dynamics of RNA polymerase (RNAP) and the anti-terminator NusA. Although NusA does not bind directly to DNA, but rather to RNAP and/or to the nascent transcript, we demonstrate that NusA interacts with RNAP ubiquitously throughout the chromosome, and that its profile mirrors RNAP distribution in both the exponential and stationary phases of growth. Generally, promoter-proximal peaks for RNAP and NusA were observed, followed by a decrease in signal strength reflecting transcriptional polarity. Differential binding of RNAP and NusA in the two growth conditions correlated with transcriptional activity as reflected by RNA abundance. Indeed, a significant association between expression levels and the presence of NusA throughout the gene body was detected, confirming the peculiar transcription-promoting role of NusA. Integration of the data sets pinpointed transcriptional units, mapped promoters and uncovered new anti-sense and non-coding transcripts. Highly expressed transcriptional units were situated mainly on the leading strand, despite the relatively unbiased distribution of genes throughout the genome, thus helping the replicative and transcriptional complexes to align.

Topological effects of data incompleteness of gene regulatory networks

BACKGROUND

The topological analysis of biological networks has been a prolific topic in network science during the last decade. A persistent problem with this approach is the inherent uncertainty and noisy nature of the data. One of the cases in which this situation is more marked is that of transcriptional regulatory networks (TRNs) in bacteria. The datasets are incomplete because regulatory pathways associated to a relevant fraction of bacterial genes remain unknown. Furthermore, direction, strengths and signs of the links are sometimes unknown or simply overlooked. Finally, the experimental approaches to infer the regulations are highly heterogeneous, in a way that induces the appearance of systematic experimental-topological correlations. And yet, the quality of the available data increases constantly.

RESULTS

In this work we capitalize on these advances to point out the influence of data (in)completeness and quality on some classical results on topological analysis of TRNs, specially regarding modularity at different levels.

CONCLUSIONS

In doing so, we identify the most relevant factors affecting the validity of previous findings, highlighting important caveats to future prokaryotic TRNs topological analysis.

Long-range transcriptional control of an operon necessary for virulence-critical ESX-1 secretion in Mycobacterium tuberculosis

The ESX-1 secretion system of Mycobacterium tuberculosis has to be precisely regulated since the secreted proteins, although required for a successful virulent infection, are highly antigenic and their continued secretion would alert the immune system to the infection. The transcription of a five-gene operon containing espACD-Rv3613c-Rv3612c, which is required for ESX-1 secretion and is essential for virulence, was shown to be positively regulated by the EspR transcription factor. Thus, transcription from the start site, found to be located 67 bp upstream of espA, was dependent upon EspR enhancer-like sequences far upstream (between 884 and 1,004 bp), which we term the espA activating region (EAR). The EAR contains one of the known binding sites for EspR, providing the first in vivo evidence that transcriptional activation at the espA promoter occurs by EspR binding to the EAR and looping out DNA between this site and the promoter. Regulation of transcription of this operon thus takes place over long regions of the chromosome. This regulation may differ in some members of the M. tuberculosis complex, including Mycobacterium bovis, since deletions of the intergenic region have removed the upstream sequence containing the EAR, resulting in lowered espA expression. Consequent differences in expression of ESX-1 in these bacteria may contribute to their various pathologies and host ranges. The virulence-critical nature of this operon means that transcription factors controlling its expression are possible drug targets.

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 accessed at http://www.bifmanit.org/tb2/.

Mycoplasma hyopneumoniae transcription unit organization: genome survey and prediction

Mycoplasma hyopneumoniae is associated with swine respiratory diseases. Although gene organization and regulation are well known in many prokaryotic organisms, knowledge on mycoplasma is limited. This study performed a comparative analysis of three strains of M. hyopneumoniae (7448, J and 232), with a focus on genome organization and gene comparison for open read frame (ORF) cluster (OC) identification. An in silico analysis of gene organization demonstrated 117 OCs and 34 single ORFs in M. hyopneumoniae 7448 and J, while 116 OCs and 36 single ORFs were identified in M. hyopneumoniae 232. Genomic comparison revealed high synteny and conservation of gene order between the OCs defined for 7448 and J strains as well as for 7448 and 232 strains. Twenty-one OCs were chosen and experimentally confirmed by reverse transcription–PCR from M. hyopneumoniae 7448 genome, validating our prediction. A subset of the ORFs within an OC could be independently transcribed due to the presence of internal promoters. Our results suggest that transcription occurs in ‘run-on’ from an upstream promoter in M. hyopneumoniae, thus forming large ORF clusters (from 2 to 29 ORFs in the same orientation) and indicating a complex transcriptional organization.

The transcriptional regulatory network of Mycobacterium tuberculosis

Under the perspectives of network science and systems biology, the characterization of transcriptional regulatory (TR) networks beyond the context of model organisms offers a versatile tool whose potential remains yet mainly unexplored. In this work, we present an updated version of the TR network of Mycobacterium tuberculosis (M.tb), which incorporates newly characterized transcriptional regulations coming from 31 recent, different experimental works available in the literature. As a result of the incorporation of these data, the new network doubles the size of previous data collections, incorporating more than a third of the entire genome of the bacterium. We also present an exhaustive topological analysis of the new assembled network, focusing on the statistical characterization of motifs significances and the comparison with other model organisms. The expanded M.tb transcriptional regulatory network, considering its volume and completeness, constitutes an important resource for diverse tasks such as dynamic modeling of gene expression and signaling processes, computational reliability determination or protein function prediction, being the latter of particular relevance, given that the function of only a small percent of the proteins of M.tb is known.

Comprehensive insights into Mycobacterium tuberculosis DevR (DosR) regulon activation switch

DevR regulon function is believed to be crucial for the survival of Mycobacterium tuberculosis during dormancy. In this study, we undertook a comprehensive analysis of the DevR regulon. All the regulon promoters were assigned to four classes based on the number of DevR binding sites (Dev boxes). A minimum of two boxes are essential for complete interaction and their tandem arrangement is an architectural hallmark at all promoters. Initial interaction of DevR with the conserved box is essential for its cooperative binding to adjacent sites bearing low to very poor sequence conservation and is the universal mechanism underlying DevR-mediated transcriptional induction. The functional importance of tandem arrangement was established by analyzing promoter variants harboring Dev boxes with altered spacing. Conserved sequence logos were generated from 47 binding sequences which included 24 newly discovered Dev boxes. In each half site of an 18-bp binding motif, G₅ and C₇ are essential for DevR binding. Finally, we show that DevR regulon induction occurs in a temporal manner and genes that are induced early are also usually powerfully induced. The information theory-based approach along with binding and temporal expression studies provide us with comprehensive insights into the complex pattern of DevR regulon activation.

Characterization of a novel heat shock protein (Hsp22.5) involved in the pathogenesis of Mycobacterium tuberculosis

Tuberculosis is a worldwide health problem, given that one-third of the world's population is currently infected with Mycobacterium tuberculosis. Understanding the regulation of virulence on the molecular level will provide a better understanding of how M. tuberculosis can establish chronic infection. Using in vivo microarray analysis (IVMA), we previously identified a group of genes that are activated in BALB/c mouse lungs compared to in vitro cultures, including the rv0990c gene. Our analysis indicated that this gene is a member of the heat shock regulon and was activated under other stress conditions, including survival in macrophages or during the late phase of chronic tuberculosis in the murine lungs. Deletion of rv0990c from the genome of M. tuberculosis strain H37Rv affected the transcriptional profiles of many genes (n = 382) and operons involved in mycobacterial survival, including the dormancy regulon, ATP synthesis, respiration, protein synthesis, and lipid metabolism. Comparison of the proteomes of the mutant to those of the wild-type strain further confirmed the differential expression of 15 proteins, especially those involved in the heat shock response (e.g., DnaK and GrpE). Finally, the rv0990c mutant strain showed survival equivalent to that of the isogenic wild-type strain during active tuberculosis in guinea pigs, despite showing significant attenuation in BALB/c mice during the chronic phase of the disease. Overall, we suggest that rv0990c encodes a heat shock protein that plays an important role in mycobacterial virulence. Hence, we renamed rv0990c heat shock protein 22.5 (hsp22.5), reflecting its molecular mass.

ClgR regulation of chaperone and protease systems is essential for Mycobacterium tuberculosis parasitism of the macrophage

Chaperone and protease systems play essential roles in cellular homeostasis and have vital functions in controlling the abundance of specific cellular proteins involved in processes such as transcription, replication, metabolism and virulence. Bacteria have evolved accurate regulatory systems to control the expression and function of chaperones and potentially destructive proteases. Here, we have used a combination of transcriptomics, proteomics and targeted mutagenesis to reveal that the clp gene regulator (ClgR) of Mycobacterium tuberculosis activates the transcription of at least ten genes, including four that encode protease systems (ClpP1/C, ClpP2/C, PtrB and HtrA-like protease Rv1043c) and three that encode chaperones (Acr2, ClpB and the chaperonin Rv3269). Thus, M. tuberculosis ClgR controls a larger network of protein homeostatic and regulatory systems than ClgR in any other bacterium studied to date. We demonstrate that ClgR-regulated transcriptional activation of these systems is essential for M. tuberculosis to replicate in macrophages. Furthermore, we observe that this defect is manifest early in infection, as M. tuberculosis lacking ClgR is deficient in the ability to control phagosome pH 1 h post-phagocytosis.

High content phenotypic cell-based visual screen identifies Mycobacterium tuberculosis acyltrehalose-containing glycolipids involved in phagosome remodeling

The ability of the tubercle bacillus to arrest phagosome maturation is considered one major mechanism that allows its survival within host macrophages. To identify mycobacterial genes involved in this process, we developed a high throughput phenotypic cell-based assay enabling individual sub-cellular analysis of over 11,000 Mycobacterium tuberculosis mutants. This very stringent assay makes use of fluorescent staining for intracellular acidic compartments, and automated confocal microscopy to quantitatively determine the intracellular localization of M. tuberculosis. We characterised the ten mutants that traffic most frequently into acidified compartments early after phagocytosis, suggesting that they had lost their ability to arrest phagosomal maturation. Molecular analysis of these mutants revealed mainly disruptions in genes involved in cell envelope biogenesis (fadD28), the ESX-1 secretion system (espL/Rv3880), molybdopterin biosynthesis (moaC1 and moaD1), as well as in genes from a novel locus, Rv1503c-Rv1506c. Most interestingly, the mutants in Rv1503c and Rv1506c were perturbed in the biosynthesis of acyltrehalose-containing glycolipids. Our results suggest that such glycolipids indeed play a critical role in the early intracellular fate of the tubercle bacillus. The unbiased approach developed here can be easily adapted for functional genomics study of intracellular pathogens, together with focused discovery of new anti-microbials.

Mycobacterium tuberculosis MycP1 protease plays a dual role in regulation of ESX-1 secretion and virulence

Mycobacterium tuberculosis uses the ESX-1 secretion system to deliver virulence proteins during infection of host cells. Here we report a mechanism of posttranscriptional control of ESX-1 mediated by MycP1, a M. tuberculosis serine protease. We show that MycP1 is required for ESX-1 secretion but that, unexpectedly, genetic inactivation of MycP1 protease activity increases secretion of ESX-1 substrates. We demonstrate that EspB, an ESX-1 substrate required for secretion, is a target of MycP1 in vitro and in vivo. During macrophage infection, an inactive MycP1 protease mutant causes hyperactivation of ESX-1-stimulated innate signaling pathways. MycP1 is required for growth in mice during acute infection, while loss of its protease activity leads to attenuated virulence during chronic infection. As the key ESX-1 substrates ESAT-6 and CFP-10 are highly immunogenic, fine-tuning of their secretion by MycP1 may balance virulence and immune detection and be essential for successful maintenance of long-term M. tuberculosis infection.

Functional analysis of an intergenic non-coding sequence within mce1 operon of M.tuberculosis

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.

Results

We predicted putative promoter activity of the 200 base pairs of non-coding, intergenic region between Rv0166 and Rv0167 in silico using MEME software and designate it as intergenic promoter, IGPr. We demonstrate both promoter activity and a putative negative regulatory function of this fragment by reporter assays carried out in the surrogate host M.smegmatis. We find that the repressive elements not only control the native promoter but also repress a heterologous promoter of M.smegmatis. The higher activity of the intergenic promoter in a clinical isolate in comparison with the wild type sequence from M.tuberculosis H37Rv could be correlated with a point mutation within the negative element. We have mapped two transcription start sites for mce1 operon both of which are utilized in M.tuberculosis H37Rv as well as the clinical isolate VPCI591. Our studies show that the promoter activity in the non-coding region is relevant not only in reporter gene expression but also in the expression of mce1 operon in M. tuberculosis cells grown in synthetic medium.

Conclusion

The mce operon of M.tuberculosis H37Rv potentially can be transcribed from two promoters P1 and P2, former mapping upstream of Rv0166 and the latter in the non-coding intergenic region between Rv0166 and Rv0167. The transcription initiation from P1 results in a transcript with Rv0166 while that from P2 will be without it. The sequences between the translation start site of Rv0167 and the promoter P2 have a negative regulatory role, as point mutation within the sequence leads to enhanced activity of P2 as well as a heterologous promoter from M.smegmatis. The mutation detected in the clinical isolate VPCI591 therefore behaves like a gain-of-function mutation.

The transcriptome of Mycobacterium tuberculosis

In this review, we summarize the present understanding of the Mycobacterium tuberculosis transcriptome and catalog both experimentally verified findings and computationally derived predictions. We also provide a new analysis for a range of discoveries by comparing the results of previously independent research papers. Bringing these data together and improving their accessibility should help catalyze further discoveries. This minireview also provides some general insights that may be valuable to those working to characterize the transcriptome of less-studied prokaryotes.

Binary particle swarm optimization for operon prediction

An operon is a fundamental unit of transcription and contains specific functional genes for the construction and regulation of networks at the entire genome level. The correct prediction of operons is vital for understanding gene regulations and functions in newly sequenced genomes. As experimental methods for operon detection tend to be nontrivial and time consuming, various methods for operon prediction have been proposed in the literature. In this study, a binary particle swarm optimization is used for operon prediction in bacterial genomes. The intergenic distance, participation in the same metabolic pathway, the cluster of orthologous groups, the gene length ratio and the operon length are used to design a fitness function. We trained the proper values on the Escherichia coli genome, and used the above five properties to implement feature selection. Finally, our study used the intergenic distance, metabolic pathway and the gene length ratio property to predict operons. Experimental results show that the prediction accuracy of this method reached 92.1%, 93.3% and 95.9% on the Bacillus subtilis genome, the Pseudomonas aeruginosa PA01 genome and the Staphylococcus aureus genome, respectively. This method has enabled us to predict operons with high accuracy for these three genomes, for which only limited data on the properties of the operon structure exists.

An anaerobic-type alpha-ketoglutarate ferredoxin oxidoreductase completes the oxidative tricarboxylic acid cycle of Mycobacterium tuberculosis

Aerobic organisms have a tricarboxylic acid (TCA) cycle that is functionally distinct from those found in anaerobic organisms. Previous reports indicate that the aerobic pathogen Mycobacterium tuberculosis lacks detectable α-ketoglutarate (KG) dehydrogenase activity and drives a variant TCA cycle in which succinyl-CoA is replaced by succinic semialdehyde. Here, we show that M. tuberculosis expresses a CoA-dependent KG dehydrogenase activity, albeit one that is typically found in anaerobic bacteria. Unlike most enzymes of this family, the M. tuberculosis KG: ferredoxin oxidoreductase (KOR) is extremely stable under aerobic conditions. This activity is absent in a mutant strain deleted for genes encoding a previously uncharacterized oxidoreductase, and this strain is impaired for aerobic growth in the absence of sufficient amounts of CO₂. Interestingly, inhibition of the glyoxylate shunt or exclusion of exogenous fatty acids alleviates this growth defect, indicating the presence of an alternate pathway that operates in the absence of β-oxidation. Simultaneous disruption of KOR and the first enzyme of the succinic semialdehyde pathway (KG decarboxylase; KGD) results in strict dependence upon the glyoxylate shunt for growth, demonstrating that KG decarboxylase is also functional in M. tuberculosis intermediary metabolism. These observations demonstrate that unlike most organisms M. tuberculosis utilizes two distinct TCA pathways from KG, one that functions concurrently with β-oxidation (KOR-dependent), and one that functions in the absence of β-oxidation (KGD-dependent). As these pathways are regulated by metabolic cues, we predict that their differential utilization provides an advantage for growth in different environments within the host.

Knowledge of the basic biology of Mycobacterium tuberculosis is essential to identifying novel ways to combat the emerging threat of drug-resistant tuberculosis. Since the tricarboxylic acid (TCA) cycle is a cornerstone of metabolism and M. tuberculosis does not possess a “typical” TCA cycle enzyme set, much effort has been focused on elucidating the components of this pathway. Previous reports indicate that M. tuberculosis possesses a variant TCA cycle in which succinic semialdehyde replaces succinyl-CoA. Since this pathway does not conserve as much metabolic energy as the canonical pathway, we considered an alternative hypothesis: that M. tuberculosis might possess an anaerobic type α-ketoglutarate dehydrogenase. In this manuscript, we investigate this previously unknown activity for mycobacteria using a combination of genetic and biochemical approaches, and demonstrate that M. tuberculosis is capable of driving a conventional TCA cycle in an unconventional way. We also validate the existence of the previously described variant pathway and provide evidence that these two pathways are differentially utilized in response to a metabolic signal, fatty acid catabolism.

Unusual diheme conformation of the heme-degrading protein from Mycobacterium tuberculosis

Heme degradation plays a pivotal role in the availability of the essential nutrient, iron, in pathogenic bacteria. A previously unannotated protein from Mycobacterium tuberculosis, Rv3592, which shares homology to heme-degrading enzymes, has been identified. Biochemical analyses confirm that Rv3592, which we have termed MhuD (mycobacterial heme utilization, degrader), is able to bind and degrade heme. Interestingly, contrary to previously reported stoichiometry for the Staphylococcus aureus heme degraders, iron-regulated surface determinant (Isd)G and IsdI, MhuD has the ability to bind heme in a 1:2 protein-to-heme ratio, although the MhuD-diheme complex is inactive. Furthermore, the 1.75-A crystal structure of the MhuD-diheme complex reveals two stacked hemes forming extensive contacts with residues in the active site. In particular, the solvent-exposed heme is axially liganded by His75 and is stacked planar upon the solvent-protected heme. The solvent-protected heme is coordinated by a chloride ion, which is, in turn, stabilized by Asn7. Structural comparison between MhuD-diheme and inactive IsdG and IsdI bound to only one highly distorted metalloporphyrin ring reveals that several residues located in alpha-helix 2 and the subsequent loop appear to be responsible for heme stoichiometric differences and suggest open and closed conformations for substrate entry and product exit.