Integrated Modeling of Gene Regulatory and Metabolic Networks in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) is the causative bacterium of tuberculosis, a disease responsible for over a million deaths worldwide annually with a growing number of strains resistant to antibiotics. The development of better therapeutics would greatly benefit from improved understanding of the mechanisms associated with MTB responses to different genetic and environmental perturbations. Therefore, we expanded a genome-scale regulatory-metabolic model for MTB using the Probabilistic Regulation of Metabolism (PROM) framework. Our model, MTBPROM2.0, represents a substantial knowledge base update and extension of simulation capability. We incorporated a recent ChIP-seq based binding network of 2555 interactions linking to 104 transcription factors (TFs) (representing a 3.5-fold expansion of TF coverage). We integrated this expanded regulatory network with a refined genome-scale metabolic model that can correctly predict growth viability over 69 source metabolite conditions and predict metabolic gene essentiality more accurately than the original model. We used MTBPROM2.0 to simulate the metabolic consequences of knocking out and overexpressing each of the 104 TFs in the model. MTBPROM2.0 improves performance of knockout growth defect predictions compared to the original PROM MTB model, and it can successfully predict growth defects associated with TF overexpression. Moreover, condition-specific models of MTBPROM2.0 successfully predicted synergistic growth consequences of overexpressing the TF whiB4 in the presence of two standard anti-TB drugs. MTBPROM2.0 can screen in silico condition-specific transcription factor perturbations to generate putative targets of interest that can help prioritize future experiments for therapeutic development efforts.

Osmosensory signaling in Mycobacterium tuberculosis mediated by a eukaryotic-like Ser/Thr protein kinase.

Bacteria are able to adapt to dramatically different microenvironments, but in many organisms, the signaling pathways, transcriptional programs, and downstream physiological changes involved in adaptation are not well-understood. Here, we discovered that osmotic stress stimulates a signaling network in Mycobacterium tuberculosis regulated by the eukaryotic-like receptor Ser/Thr protein kinase PknD. Expression of the PknD substrate Rv0516c was highly induced by osmotic stress. Furthermore, Rv0516c disruption modified peptidoglycan thickness, enhanced antibiotic resistance, and activated genes in the regulon of the alternative σ-factor SigF. Phosphorylation of Rv0516c regulated the abundance of EspA, a virulence-associated substrate of the type VII ESX-1 secretion system. These findings identify an osmosensory pathway orchestrated by PknD, Rv0516c, and SigF that enables adaptation to osmotic stress through cell wall remodeling and virulence factor production. Given the widespread occurrence of eukaryotic-like Ser/Thr protein kinases in bacteria, these proteins may play a broad role in bacterial osmosensing.

Global analysis of mRNA stability in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) is a highly successful pathogen that infects over a billion people. As with most organisms, MTB adapts to stress by modifying its transcriptional profile. Remodeling of the transcriptome requires both altering the transcription rate and clearing away the existing mRNA through degradation, a process that can be directly regulated in response to stress. To understand better how MTB adapts to the harsh environs of the human host, we performed a global survey of the decay rates of MTB mRNA transcripts. Decay rates were measured for 2139 of the ~4000 MTB genes, which displayed an average half-life of 9.5 min. This is nearly twice the average mRNA half-life of other prokaryotic organisms where these measurements have been made. The transcriptome was further stabilized in response to lowered temperature and hypoxic stress. The generally stable transcriptome described here, and the additional stabilization in response to physiologically relevant stresses, has far-ranging implications for how this pathogen is able to adapt in its human host.

Bioinformatic and empirical analysis of novel hypoxia-inducible targets of the human antituberculosis T cell response.

We analyzed whole genome-based transcriptional profiles of Mycobacterium tuberculosis subjected to prolonged hypoxia to guide the discovery of novel potential Ags, by a combined bioinformatic and empirical approach. We analyzed the fold induction of the 100 most highly induced genes at 7 d of hypoxia, as well as transcript abundance, peptide-binding prediction (ProPred) adjusted for population-specific MHC class II allele frequency, and by literature search. Twenty-six candidate genes were selected by this bioinformatic approach and evaluated empirically using IFN-γ and IL-2 ELISPOT using immunodominant Ags (Acr-1, CFP-10, ESAT-6) as references. Twenty-three of twenty-six proteins induced an IFN-γ response in PBMCs of persons with active or latent tuberculosis. Five novel immunodominant proteins-Rv1957, Rv1954c, Rv1955, Rv2022c, and Rv1471-were identified that induced responses similar to CFP-10 and ESAT-6 in both magnitude and frequency. IL-2 responses were of lower magnitude than were those of IFN-γ. Only moderate evidence of infection stage-specific recognition of Ags was observed. Reconciliation of bioinformatic and empirical hierarchies of immunodominance revealed that Ags could be predicted, providing transcriptomic data were combined with peptide-binding prediction adjusted by population-specific MHC class II allele frequency.

Hypoxia induces an immunodominant target of tuberculosis specific T cells absent from common BCG vaccines.

M. tuberculosis (MTB) species-specific antigenic determinants of the human T cell response are important for immunodiagnosis and vaccination. As hypoxia is a stimulus in chronic tuberculosis infection, we analyzed transcriptional profiles of MTB subject to 168 hours of hypoxia to test the hypothesis that upregulation by hypoxia might result in gene products being recognized as antigens. We identified upregulation of two region of difference (RD) 11 (Rv2658C and Rv2659c), and one RD2 (Rv1986) absent from commonly used BCG strains. In MTB infected persons, the IL-2 ELISpot response to Rv1986 peptides was several times greater than the corresponding IFN-γ response to the reference immunodominant ESAT-6 or CFP-10 antigens. The IL-2 response was confined to two epitopic regions containing residues 61-80 and 161-180. The biggest population of IL-2 secreting T cells was single cytokine positive central memory T cells. The IL-2 response to live MTB bacilli lacking Rv1986 was significantly lower than the response to wild type or mutant complemented with Rv1986. In addition, the IL-2 response to Rv1986 was significantly lower in HIV-TB co-infected persons than in HIV uninfected persons, and significantly increased during antiretroviral therapy. These findings demonstrate that Rv1986 is an immunodominant target of memory T cells and is therefore of relevance when considering the partial efficacy of currently used BCG vaccines and provide evidence for a clinical trial comparing BCG strains.

Transcriptional regulator-induced phenotype screen reveals drug potentiators in Mycobacterium tuberculosis.

Transposon-based strategies provide a powerful and unbiased way to study the bacterial stress response1-8, but these approaches cannot fully capture the complexities of network-based behaviour. Here, we present a network-based genetic screening approach: the transcriptional regulator-induced phenotype (TRIP) screen, which we used to identify previously uncharacterized network adaptations of Mycobacterium tuberculosis to the first-line anti-tuberculosis drug isoniazid (INH). We found regulators that alter INH susceptibility when induced, several of which could not be identified by standard gene disruption approaches. We then focused on a specific regulator, mce3R, which potentiated INH activity when induced. We compared mce3R-regulated genes with baseline INH transcriptional responses and implicated the gene ctpD (Rv1469) as a putative INH effector. Evaluating a ctpD disruption mutant demonstrated a previously unknown role for this gene in INH susceptibility. Integrating TRIP screening with network information can uncover sophisticated molecular response programs.

Ultra-low Dose Aerosol Infection of Mice with Mycobacterium tuberculosis More Closely Models Human Tuberculosis.

Tuberculosis (TB) is a heterogeneous disease manifesting in a subset of individuals infected with aerosolized Mycobacterium tuberculosis (Mtb). Unlike human TB, murine infection results in uniformly high lung bacterial burdens and poorly organized granulomas. To develop a TB model that more closely resembles human disease, we infected mice with an ultra-low dose (ULD) of between 1-3 founding bacteria, reflecting a physiologic inoculum. ULD-infected mice exhibited highly heterogeneous bacterial burdens, well-circumscribed granulomas that shared features with human granulomas, and prolonged Mtb containment with unilateral pulmonary infection in some mice. We identified blood RNA signatures in mice infected with an ULD or a conventional Mtb dose (50-100 CFU) that correlated with lung bacterial burdens and predicted Mtb infection outcomes across species, including risk of progression to active TB in humans. Overall, these findings highlight the potential of the murine TB model and show that ULD infection recapitulates key features of human TB.

Hypoxia: a window into Mycobacterium tuberculosis latency.

Tuberculosis is a massive public health problem on a global scale and the success of Mycobacterium tuberculosis is linked to its ability to persist within humans for long periods without causing any overt disease symptoms. Hypoxia is predicted to be a key host-induced stress limiting growth of the pathogen in vivo. However, multiple studies in vitro and in vivo indicate that M. tuberculosis adapts to oxygen limitation by entering into a metabolically altered state, while awaiting the opportunity to reactivate. Molecular signatures of bacteria adapted to hypoxia in vitro are accumulating, although correlations to human disease are only now being established. Similarly, defining the mechanisms that control this adaptation is an active area of research. In this review we discuss the historical precedents linking hypoxia and latency, and the gathering knowledge of M. tuberculosis hypoxic responses. We also examine the role of these responses in tuberculosis latency, and identify promising avenues for future studies.

A comprehensive map of genome-wide gene regulation in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) is a pathogenic bacterium responsible for 12 million active cases of tuberculosis (TB) worldwide. The complexity and critical regulatory components of MTB pathogenicity are still poorly understood despite extensive research efforts. In this study, we constructed the first systems-scale map of transcription factor (TF) binding sites and their regulatory target proteins in MTB. We constructed FLAG-tagged overexpression constructs for 206 TFs in MTB, used ChIP-seq to identify genome-wide binding events and surveyed global transcriptomic changes for each overexpressed TF. Here we present data for the most comprehensive map of MTB gene regulation to date. We also define elaborate quality control measures, extensive filtering steps, and the gene-level overlap between ChIP-seq and microarray datasets. Further, we describe the use of TF overexpression datasets to validate a global gene regulatory network model of MTB and describe an online source to explore the datasets.

The DNA-binding network of Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) infects 30% of all humans and kills someone every 20-30 s. Here we report genome-wide binding for ~80% of all predicted MTB transcription factors (TFs), and assayed global expression following induction of each TF. The MTB DNA-binding network consists of ~16,000 binding events from 154 TFs. We identify >50 TF-DNA consensus motifs and >1,150 promoter-binding events directly associated with proximal gene regulation. An additional ~4,200 binding events are in promoter windows and represent strong candidates for direct transcriptional regulation under appropriate environmental conditions. However, we also identify >10,000 'dormant' DNA-binding events that cannot be linked directly with proximal transcriptional control, suggesting that widespread DNA binding may be a common feature that should be considered when developing global models of coordinated gene expression.

Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network.

BACKGROUND: Mycobacterium tuberculosis senses and responds to the shifting and hostile landscape of the host. To characterize the underlying intertwined gene regulatory network governed by approximately 200 transcription factors of M. tuberculosis, we have assayed the global transcriptional consequences of overexpressing each transcription factor from an inducible promoter. RESULTS: We cloned and overexpressed 206 transcription factors in M. tuberculosis to identify the regulatory signature of each. We identified 9,335 regulatory consequences of overexpressing each of 183 transcription factors, providing evidence of regulation for 70% of the M. tuberculosis genome. These transcriptional signatures agree well with previously described M. tuberculosis regulons. The number of genes differentially regulated by transcription factor overexpression varied from hundreds of genes to none, with the majority of expression changes repressing basal transcription. Exploring the global transcriptional maps of transcription factor overexpressing (TFOE) strains, we predicted and validated the phenotype of a regulator that reduces susceptibility to a first line anti-tubercular drug, isoniazid. We also combined the TFOE data with an existing model of M. tuberculosis metabolism to predict the growth rates of individual TFOE strains with high fidelity. CONCLUSION: This work has led to a systems-level framework describing the transcriptome of a devastating bacterial pathogen, characterized the transcriptional influence of nearly all individual transcription factors in M. tuberculosis, and demonstrated the utility of this resource. These results will stimulate additional systems-level and hypothesis-driven efforts to understand M. tuberculosis adaptations that promote disease.

Characterization of a Clp protease gene regulator and the reaeration response in Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB) enters a non-replicating state when exposed to low oxygen tension, a condition the bacillus encounters in granulomas during infection. Determining how mycobacteria enter and maintain this state is a major focus of research. However, from a public health standpoint the importance of latent TB is its ability to reactivate. The mechanism by which mycobacteria return to a replicating state upon re-exposure to favorable conditions is not understood. In this study, we utilized reaeration from a defined hypoxia model to characterize the adaptive response of MTB following a return to favorable growth conditions. Global transcriptional analysis identified the approximately 100 gene Reaeration Response, induced relative to both log-phase and hypoxic MTB. This response includes chaperones and proteases, as well as the transcription factor Rv2745c, which we characterize as a Clp protease gene regulator (ClgR) orthologue. During reaeration, genes repressed during hypoxia are also upregulated in a wave of transcription that includes genes crucial to transcription, translation and oxidative phosphorylation and culminates in bacterial replication. In sum, this study defines a new transcriptional response of MTB with potential relevance to disease, and implicates ClgR as a regulator involved in resumption of replication following hypoxia.

Isolation of mycobacterial RNA.

This chapter describes two protocols for isolating total RNA from mycobacteria: one for extraction from in vitro cultures and one for extraction from in vivo. In these protocols, RNA is liberated from mycobacteria by disruption with small glass beads in the presence of Trizol to stabilize the RNA. The RNA is further purified with DNAse treatment and RNeasy columns. This protocol leads to microgram quantities of RNA from log-phase cultures.

The enduring hypoxic response of Mycobacterium tuberculosis.

BACKGROUND: A significant body of evidence accumulated over the last century suggests a link between hypoxic microenvironments within the infected host and the latent phase of tuberculosis. Studies to test this correlation have identified the M. tuberculosis initial hypoxic response, controlled by the two-component response regulator DosR. The initial hypoxic response is completely blocked in a dosR deletion mutant. METHODOLOGY/PRINCIPAL FINDINGS: We show here that a dosR deletion mutant enters bacteriostasis in response to in vitro hypoxia with only a relatively mild decrease in viability. In the murine infection model, the phenotype of the mutant was indistinguishable from that of the parent strain. These results suggested that additional genes may be essential for entry into and maintenance of bacteriostasis. Detailed microarray analysis of oxygen starved cultures revealed that DosR regulon induction is transient, with induction of nearly half the genes returning to baseline within 24 hours. In addition, a larger, sustained wave of gene expression follows the DosR-mediated initial hypoxic response. This Enduring Hypoxic Response (EHR) consists of 230 genes significantly induced at four and seven days of hypoxia but not at initial time points. These genes include a surprising number of transcriptional regulators that could control the program of bacteriostasis. We found that the EHR is independent of the DosR-mediated initial hypoxic response, as EHR expression is virtually unaltered in the dosR mutant. CONCLUSIONS/SIGNIFICANCE: Our results suggest a reassessment of the role of DosR and the initial hypoxic response in MTB physiology. Instead of a primary role in survival of hypoxia induced bacteriostasis, DosR may regulate a response that is largely optional in vitro and in mouse infections. Analysis of the EHR should help elucidate the key regulatory factors and enzymatic machinery exploited by M. tuberculosis for long-term bacteriostasis in the face of oxygen deprivation.

The Candida albicans mating type like locus [MTL] is not involved in chlamydospore formation.

Candida albicans produces chlamydospores, which can be used as a diagnostic tool for species identification. It has been suggested that these chlamydospores are degenerate spores. If so, then their production might be linked to the mating loci, and clinical strains that are homozygous for the C. albicans mating locus MTL may be altered in chlamydospore formation, which could cause problems in diagnostics and species identification. In Saccharomyces cerevisiae diploid cells, the heterodimeric transcriptional repressor formed by the products of the mating genes MATa1 and MATalpha2 is an important regulator of sporulation. It was therefore of interest to determine if the disruptions of the MATa1 and MATalpha2 homologs in C. albicans, MTLa1 and MTLalpha2, result in inhibition of chlamydospore formation. Laboratory strains containing disruptions of either the entire MTL locus or specific genes within the locus were assayed for their ability to form chlamydospores. Clinical strains that are homozygous for one of the two MTL loci were also assayed. No change in chlamydospore formation was seen in these strains compared to the standard laboratory strain.

Homozygosity at the Candida albicans MTL locus associated with azole resistance.

Antifungal drug resistance in the pathogenic fungus Candida albicans is a serious threat to the growing population of immunocompromised patients. This study describes a significant correlation between loss of heterozygosity at the C. albicans mating-type-like (MTL) locus and resistance to azole antifungals. A pool of 96 clinical isolates consisting of 50 azole-resistant or susceptible dose-dependent isolates and 46 azole-susceptible isolates was screened by PCR for the presence of MTLa1 and MTLalpha1. These genes were used as markers for the MTLa and MTLalpha loci. Both loci were present in 84 of the isolates. Six isolates failed to amplify MTLa1 and six failed to amplify MTLalpha1. Further PCR analysis demonstrated that loss of the MTLa1 and MTLalpha1 genes corresponded to loss of all of the loci-specific genes, resulting in homozygosity at the MTL locus. Southern analysis and single nucleotide polymorphism (SNP) analysis were used to determine that this loss of heterogeneity was due to replacement of one of the MTL loci with a duplicate of the other locus resulting in two homozygous copies of the MTL locus. Of the 12 homozygous isolates, one isolate was sensitive to azole drugs. Statistical analysis of the data demonstrates a strong correlation between homozygosity at the MTL locus and azole resistance (P<0 small middle dot003). In a set of serial isolates, an increase in azole resistance correlated with the loss of heterozygosity at the MTL locus, lending further strength to the correlation. Gene disruptions of the MTL loci were found to have no effect on azole susceptibility.

Genes required for intrinsic multidrug resistance in Mycobacterium avium.

Genes required for intrinsic multidrug resistance by Mycobacterium avium were identified by screening a library of transposon insertion mutants for the inability to grow in the presence of ciprofloxacin, clarithromycin, and penicillin at subinhibitory concentrations. Two genes, pks12 and Maa2520, were disrupted in multiple drug-susceptible mutants. The pks12 gene (Maa1979), which may be cotranscribed with a downstream gene (Maa1980), is widely conserved in the actinomycetes. Its ortholog in Mycobacterium tuberculosis is a polyketide synthase required for the synthesis of dimycocerosyl phthiocerol, a major cell wall lipid. Mutants of M. avium with insertions into pks12 exhibited altered colony morphology and were drug susceptible, but they grew as well as the wild type did in vitro and intracellularly within THP-1 cells. A pks12 mutant of M. tuberculosis was moderately more susceptible to clarithromycin than was its parent strain; however, susceptibility to ciprofloxacin and penicillin was not altered. M. avium complex (MAC) and M. tuberculosis appear to have different genetic mechanisms for resisting the effects of these antibiotics, with pks12 playing a relatively more significant role in MAC. The second genetic locus identified in this study, Maa2520, is a conserved hypothetical gene with orthologs in M. tuberculosis and Mycobacterium leprae. It is immediately upstream of Maa2521, which may code for an exported protein. Mutants with insertions at this locus were susceptible to multiple antibiotics and slow growing in vitro and were unable to survive intracellularly within THP-1 cells. Like pks12 mutants, they exhibited increased Congo red binding, an indirect indication of cell wall modifications. Maa2520 and pks12 are the first genes to be linked by mutation to intrinsic drug resistance in MAC.