PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plospathPLoS PathogensSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Pathog. Jun 2012; 8(6): e1002769.
Published online Jun 21, 2012. doi:  10.1371/journal.ppat.1002769
PMCID: PMC3380936
Linking the Transcriptional Profiles and the Physiological States of Mycobacterium tuberculosis during an Extended Intracellular Infection
Kyle H. Rohde,1*¤ Diogo F. T. Veiga,2 Shannon Caldwell,1 Gábor Balázsi,2 and David G. Russell1
1Department of Microbiology and Immunology, Cornell University, Ithaca, New York, United States of America
2Department of Systems Biology-Unit 950, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Christopher M. Sassetti, Editor
University of Massachusetts, United States of America
* E-mail: Kyle.Rohde/at/ucf.edu
¤: Current address: Burnett School of Biomedical Sciences, College of Medicine, University of Central Florida, Orlando, Florida, United States of America.
Conceived and designed the experiments: KHR DFTV GB DGR. Performed the experiments: KHR SC. Analyzed the data: KHR DFTV SC. Contributed reagents/materials/analysis tools: KHR DFTV GB DGR. Wrote the paper: KHR DFTV. Revision and editing of final manuscript: KHR DFTV SC GB DGR.
Received October 9, 2011; Accepted May 7, 2012.
Intracellular pathogens such as Mycobacterium tuberculosis have evolved strategies for coping with the pressures encountered inside host cells. The ability to coordinate global gene expression in response to environmental and internal cues is one key to their success. Prolonged survival and replication within macrophages, a key virulence trait of M. tuberculosis, requires dynamic adaptation to diverse and changing conditions within its phagosomal niche. However, the physiological adaptations during the different phases of this infection process remain poorly understood. To address this knowledge gap, we have developed a multi-tiered approach to define the temporal patterns of gene expression in M. tuberculosis in a macrophage infection model that extends from infection, through intracellular adaptation, to the establishment of a productive infection. Using a clock plasmid to measure intracellular replication and death rates over a 14-day infection and electron microscopy to define bacterial integrity, we observed an initial period of rapid replication coupled with a high death rate. This was followed by period of slowed growth and enhanced intracellular survival, leading finally to an extended period of net growth. The transcriptional profiles of M. tuberculosis reflect these physiological transitions as the bacterium adapts to conditions within its host cell. Finally, analysis with a Transcriptional Regulatory Network model revealed linked genetic networks whereby M. tuberculosis coordinates global gene expression during intracellular survival. The integration of molecular and cellular biology together with transcriptional profiling and systems analysis offers unique insights into the host-driven responses of intracellular pathogens such as M. tuberculosis.
The impact of Mycobacterium tuberculosis on global health is undeniable, with ~2 million deaths and ~9 million new cases of tuberculosis each year. A key to the success of M. tuberculosis as a persistent, intracellular pathogen is its ability to survive for extended periods within professional phagocytes. Sustained growth within macrophage phagosomes requires avoiding or resisting antimicrobial mechanisms and adapting to replicate in a stressful, nutrient-restricted environment. Our understanding of the survival strategies, metabolism, and physiology of M. tuberculosis during intracellular growth remains incomplete. We employed multi-disciplinary approaches to gain new insights into adaptive responses that M. tuberculosis mobilizes to secure a productive infection. We simultaneously quantified replication and death rates, used electron microscopy to evaluate bacterial integrity, and determined the temporal changes in bacterial gene expression during a 14-day infection. By overlaying this temporal transcriptome dataset onto an extended Transcriptional Regulatory Network model, we identified regulatory pathways, stress responses, and metabolic adaptations activated during key physiological transitions over the 14 days of infection.
Mycobacterium tuberculosis (Mtb) has evolved successful strategies to survive in the human host and thwart the attempts of the immune system to eradicate it. The astounding morbidity and mortality caused by Mtb – with roughly 9 million new cases and 2 million deaths annually – is a testament to the success of this enduring pathogen. Its remarkable resistance to killing by both innate and acquired immune effectors and successful adaptation to constantly changing environments allows Mtb to maintain infections for decades. The ability of Mtb to survive and replicate within macrophages (M[var phi]) is central to the pathogenesis of tuberculosis [1].
A considerable body of literature has been devoted to understanding how Mtb survives its interaction with the M[var phi] [1][5], which serves as its niche for much of the infection cycle in the human host. The phagosome into which the bacterium is internalized fails to acquire the full complement of lysosomal characteristics due to its restricted maturation; nonetheless, this environment still represents a considerable challenge to bacterial survival. It has been suggested that Mtb survives better through escape into the cytosol [6], however recent work has indicated that accessing the cytoplasm leads rapidly to host cell death [7], which suggests that survival within the phagosome remains the more stable environment to which Mtb must adapt to succeed.
A characteristic essential for Mtb's survival is its ability to sense its environment and modulate genetic pathways controlling resistance mechanisms, metabolism, and other processes in a timely manner. The importance of gene regulation during intracellular survival and in vivo infection is illustrated by the attenuated phenotypes of numerous mutants in which key transcriptional regulators have been inactivated [8][11]. In previous microarray studies focusing on early stages of infection, we have characterized the initial transcriptional changes in Mtb immediately following M[var phi] invasion [12]. In the current study we focus on the long-term survival of Mtb within resting M[var phi] over the course of a 14-day infection. During this extended timeframe the bacteria appear to experience diverse and changing environments and adopt differing physiological states. Shortly after infection, the bacterial numbers decline indicating a period of acute stress as the macrophage attempts to eradicate the infection. This is followed by a second period during which the bacterial numbers remain relatively constant, presumably while the bacterium adapts to, and tries to modulate its environment. Finally, approximately 5–6 days post-infection (p.i.) the bacterial population increases, indicating that conditions to support survival and replication have been established. This model affords us unique access to unraveling Mtb's adaptation strategies to the changing stresses and environmental cues experienced during the establishment of a productive infection. We have exploited a range of technical approaches to define the growth states of the bacterium through the course of this 14-day infection, and to correlate these physiological states with transcriptional profiles of the surviving bacteria to identify genetic networks consistent with survival and the establishment of a productive infection.
The Dynamics of Bacterial Death and Replication during Extended Macrophage Infection
The establishment of a productive infection in the macrophage is not dependent merely on the short-term response following phagocytosis but the subsequent stages of adaptation that lead to a net increase in intracellular Mtb. To appreciate the adaptations required by Mtb to enter this replicative state one has to employ a temporal approach that accurately correlates bacterial numbers, including the assessment of death versus growth, with changing transcriptional profiles. We employed three different methods to probe the growth states of Mtb during a 14-day intracellular infection model.
First, resting primary M[var phi] isolated from C57BL/6 mice were infected with Mtb at a low multiplicity of infection (MOI) of 1[ratio]1. At two-day intervals, infected monolayers were lysed and dilutions plated to determine viable colony forming units (CFU). Supernatants were also diluted and plated (prior to lysis of M[var phi]) to monitor extracellular bacilli that could contribute to reinfection. An initial 0.5 log decrease in cell-associated viable CFU at day 2 was followed by a period of minimal change (days 2–6) and then a more steady increase in CFU (days 6–14) (Figure 1A). This profile, consistent with our previous observations [13], [14], suggested an early killing phase followed by delayed adaptation of surviving organisms to survive and replicate within M[var phi] phagosomes. At every time point, extracellular Mtb totaled less than 5% of the intracellular burden (data not shown) suggesting that re-infection is minimal.
Figure 1
Figure 1
Life and death dynamics during long-term intracellular survival of Mtb.
The intramacrophage growth of Mtb was also monitored by parallel transmission electron microscopy (TEM) analysis of fixed infected M[var phi] samples (Figure 1B). Enumeration of both morphologically normal and damaged Mtb in 100 randomly chosen macrophages revealed a rapid increase in the number of detectable bacteria per cell from 1.5 Mtb/cell at 2 hr post-infection (p.i.) to 10.3 Mtb/cell by day 4 (~1.7-fold increase/day). From day 4 to day 14, the average bacterial burden increased at a notably slower rate (~0.3-fold increase/day) (Figure 1B, Figure S1A). Interestingly, this early period of rapid Mtb replication following M[var phi] invasion coincided with the marked decrease in viable CFU (Figure 1A). This indicated that early in the infection bacterial replication is rapid, but is countered by effective bacterial killing by the M[var phi].
To further validate this conclusion, we calculated Mtb growth and death rates during long-term infection of resting M[var phi] based on the loss of an unstable replication clock plasmid [15]. The proportion of intracellular Mtb retaining the clock plasmid pBP10 was determined by plating bacteria on media with or without 25 µg/ml kanamycin. The profile of pBP10 plasmid loss (Figure 1C, red) indicated an initial period of rapid replication (day 0–day 2) followed by an extended phase of much slower cell division (day 2–day 9). The cumulative bacterial burden (CBB) – the total number of Mtb live, dead, or degraded within M[var phi] during the infection – was calculated based on the mathematical model developed by Gill et al. [15]. The predicted CBB (Figure 1C, black) at day 2 p.i., closely mirrored the quantitative TEM results enumerating all detectable bacteria regardless of viability, at ~12-fold higher than the number of viable CFU determined by plating. As shown in Figure 1D, this indicates substantial replication coinciding with even greater death rates during the early phase of infection. However, following the adoption of a slower growth rate from day 2 of infection, Mtb exhibited a steady increase in viable CFU.
These data indicate that, following invasion of M[var phi], Mtb encounters a bottleneck during which the rate of bacterial killing outpaces its relatively rapid rate of replication. A similar early phase of pronounced bacterial killing was noted in vivo by Gill et al. during the first 14 days in an Mtb-mouse model of infection [15]. Following day 2, a period of apparent adaptation ensues extending to approximately day 6, during which both the rate of bacterial replication and the rate of killing decrease. The subsequent overall increase in viable CFU reflects further enhanced survival and the successful establishment of a productive infection. The shift in the balance between replication and death over time (Figure 1D) suggests that the successful adaptation of some Mtb cells to avoid killing by M[var phi]-derived effectors is at least as important as mechanisms for sustained replication within phagosomes. It is tempting to speculate that the slower growth rate of Mtb at later time points may contribute to the shift in growth[ratio]death balance, perhaps by rendering Mtb more resistant to M[var phi]-derived pressures. Studies modeling the dynamics of Mtb-host interactions within human granulomas by Segovia-Juarez et al. support this idea, showing that slow intracellular growth rates are correlative with Mtb survival [16]. This novel insight into the life[ratio]death equilibrium of Mtb during a sustained model of infection provides a physiological context for the interpretation of the global gene expression profiles discussed hereafter.
EM Analysis Reveals Considerable Heterogeneity in the Bacterial Population during Adaptation
Descriptions of host-pathogen interactions based on population-level data, such as microarrays, are often interpreted without appreciation of the heterogeneity within that population. With this in mind, we performed detailed TEM image analysis of Mtb-infected M[var phi] at 2 day intervals over 14 days simultaneously with survival assays (Figure 1) and microarrays (see below).
Consistent with previous observations [17][20], at 2 hours p.i. most single bacteria resided in a vacuole surrounded by a phagosomal membrane tightly apposed to the Mtb cell wall (Figure 2A). While intracellular Mtb proved to be quite effective at resisting fusion with lysosomes, identified following endocytic uptake of colloidal gold (Figure 3A), ~25% of single Mtb did traffic to phagolysosomes (P-L) and colocalize with colloidal gold as early as 2 hr p.i. (Figure 3B). At 2 days p.i., image analysis of Mtb-infected M[var phi] revealed an increased frequency of morphologically-damaged Mtb in large granular lysosomes (Figure 2C–D), consistent with the early phase of killing of replicating Mtb indicated by the clock plasmid experiments. Whereas 99% of bacteria in the inoculum cultures were scored as intact, by day 2 only 75% of phagocytosed Mtb appeared morphologically normal. We also noted heavily damaged Mtb that appeared as hollowed out “ghosts”, but these were not scored as damaged. The discrepancy between the magnitude of killing, calculated based on clock plasmid data (Figure 1D), and the small proportion of visibly damaged bacteria as well as the constant ratio of normal to damaged bacilli (Figure S1A) suggests a relatively efficient clearance of nonviable organisms. At later time points, the bacterial population continued to increase (Figure 2E–F, Figure S1B) with up to ~150 Mtb/M[var phi]. In contrast to a recent report that a large subset of Mtb H37Rv translocated into the cytosol of M[var phi] within 48 hrs [6], phagosomal membranes were definitively detected surrounding ~90% (88% at day 4 p.i., 86.5% at day 8 p.i.) of intracellular Mtb CDC1551 throughout the duration of the infection (25 M[var phi] and >100 Mtb examined per time point).
Figure 2
Figure 2
Electron microscopy analysis of long-term Mtb-macrophage interactions.
Figure 3
Figure 3
Mtb occupies heterogeneous intracellular niches.
There was a surprising degree of heterogeneity of compartments in which visibly intact Mtb appeared to reside during long-term infections of M[var phi], including “replicating” Mtb in typical tight niches, small fused P-L, and putative double membrane-bound autophagosomes that colocalized with gold (Figure 3, Figure S2). The diversity of intra-host environments undoubtedly contributes to heterogeneity within the bacterial population.
Mtb Transcriptional Adaptation during Long-term Macrophage Infection
The transcriptional profile of Mtb during the course of the 14 day infection should mirror the physiological states and transitions defined in the previous section. Therefore, overlaying the temporal transcriptional changes with these physiological states will generate a correlative linkage between the two datasets. We conducted microarray analysis of RNA isolated from intracellular Mtb at 2 hr p.i. and at 2 day intervals over a 14 day period (GEO accession GSE35362). Fluorescent amplified RNA (aRNA) targets from each time-point were prepared from linearly amplified total Mtb RNA and hybridized against a 2 hr “no M[var phi]” control as previously described [12], [13]. This enabled us to determine dynamic expression ratios over time relative to a single common denominator mRNA sample, in this case an uninfected control time-matched with the earliest time-point (2 hr). Semi-quantitative RT-PCR of select target genes (aprAB, hspX, bfrB, icl, groEL2, katG, whiB7) was conducted to confirm temporal array profiles ([21] and data not shown). We identified 3626 genes with significant changes in expression during extended M[var phi] infection by combining both static statistical cutoffs (p<0.05 in at least one time-point) and EDGE (Extraction and analysis of Differential Gene Expression) methodology (Figure 4A) [22]. EDGE analysis identifies temporal changes in transcript levels that would not be deemed statistically significant at any single time-point by standard p-value measurements.
Figure 4
Figure 4
Dynamic Mtb transcriptome during long-term macrophage infection.
Temporal Gene Expression Profiles during Intracellular Adaptation
Our data indicate that Mtb encounters a bottleneck during the first two days p.i. during which a subpopulation is killed before surviving bacilli enter a sustained period of enhanced survival, albeit at a reduced replication rate. Temporal dynamics of the Mtb transcriptome during intracellular adaptation highlight genes responding, and thus likely active, during specific stages of infection. Given the short half-life of bacterial mRNA [23], [24] and steady clearance of nonviable Mtb, the transcriptional profiles observed are likely derived only from the surviving bacilli.
Up to day 2: Surviving the day 2 bottleneck
The transcriptome of bacilli sampled every 2 days largely reflected this transition, with gene expression changes of maximal magnitude and complexity occurring at day 2 (Figure 4B). We propose that most of these genes are regulated in immediate response to the stressful conditions encountered upon invasion of the macrophage. The continuity of the expression profiles over time beyond this early crisis phase argues that these are specific adaptive transcriptional trends. Transcript levels of a substantial portion of this geneset then gradually declined over the next 12 days (Figure 4B, Table S1) including genes involved in general stress responses (hsp, dnaK, groEL2), stress response regulation (sigB), anaerobic respiration (type II NADH dehydrogenase, ndh2), Fe-S cluster assembly (Rv1460, 1461, 1465), regulation of acetate and glyoxylate cycle metabolism (Rv0465c) [25], propionate metabolism (Rv1129c, prpC (Rv1131), prpD (Rv1130)), and the starvation response (lat). This temporal profile indicates that cues driving these responses are likely most acute as the bacteria try to maintain rapid early intracellular growth amidst an array of host-derived pressures.
Schnappinger et al. previously reported the intra M[var phi] gene expression of Mtb strain 1254 over a short-term infection out to 2 days p.i. [26]. The identification of 88 genes that were reliably detected and induced >1.5-fold in both studies at day 2 confirms key transcriptional responses of Mtb to the phagosomal environment, the majority of which remain elevated for the duration of the 14-day infection. These included genes known to be involved in general stress responses (hsp, dnaK, clpB), carbon metabolism (pckA, icl, prpC, prpD), oxidative stress (Rv1461, 1463), iron storage (bfrB), hypoxia and nitrosative stress (Rv1738, Rv2626c), and secreted antigens (mpt70, mpt83). This consistency with previously reported datasets from short-term infection models argues the general validity of the approach and the significance of the data from later time points in our extended infection model.
We also detected 137 additional genes upregulated at day 2 that were specific to this study including 17 members of the PhoPR regulon [27] (Figure S3A, Table S2) as well as numerous “MT genes”, comprised of 371 ORFs predicted in a later annotation of the Mtb CDC1551 genome whose intracellular expression has not been reported previously (Figure S3B, Table S3). The differential expression of MT genes during long-term M[var phi] infection suggests they may play important roles in adaptation to the phagosomal environment. This hypothesis is substantiated by our recent characterization of the MT2466-MT2467-Rv2396 locus (renamed aprABC), which showed rapid induction upon acid stress or M[var phi] invasion and is known to regulate carbon flux into cell wall lipids to support intracellular survival [21].
The transient induction of a geneset including whiB7 and adjacent ORF MT3290.2 at Day 2 followed by a return to baseline levels (Figure 5A, Table S4) reveals genes active specifically during this early stage of infection. WhiB7, a member of an Actinomycete-specific family of Fe-S cluster proteins, mediates enhanced antibiotic resistance by upregulation of a defined regulon in response to antibiotics and fatty acids [28]. The role of whiB7 during intracellular survival remains to be elucidated. It appears that as Mtb replication slows and the bacterium adapts to the phagosomal environment, many stressors and corresponding transcriptional responses return to control levels.
Figure 5
Figure 5
Profiles of Coordinated Temporal Gene Regulation.
Day 2 onwards: Survival and growth beyond the Bottleneck
The temporal changes in Mtb gene expression during the intracellular adaptation phase (following day 2) will indicate the metabolic pathways and survival mechanisms required for establishment of an environment supportive to bacterial growth. For example, Figure 5B (shown in red) shows 103 genes induced early (maximal by day 2 p.i.) that remain elevated throughout the 14-day timecourse (Table S5). This subset includes genes involved in fatty acid metabolism (pckA, icl, ald, Rv0448c-0449c), KstR-dependent cholesterol metabolism (igr), secreted antigens (mpt83-Rv2874-mpt70), and regulators (Rv2359). In a previous study, we characterized the functional impact of genetic diversity on Mtb-host interactions and gene expression by analyzing the intracellular transcriptome profiles of a globally diverse panel of M. tuberculosis complex (MTC) clinical isolates [13]. We discovered a set of genes conserved and similarly regulated within M[var phi] phagosomes across all MTC clinical isolates at 24 hr p.i. which we termed the “core transcriptome”. While the induction magnitudes and kinetics varied, the majority of the 215 genes belonging to the resting M[var phi] core transcriptome also exhibited early, sustained induction profiles (Figure S4, Table S6). The conserved nature of this response across bacterial isolates and its sustained expression throughout the duration of the infection implies strongly that these genes are important for maintaining infection within M[var phi].
In contrast, looking at genes with a reciprocal pattern of sustained down-regulation, we found ~30 genes (shown in black) with a profile of early, sustained repression including numerous ribosomal genes (rplF, rplI, rplJ, rplV) (Figure 5B, Table S7). This correlates well with the decreased replication rate revealed in the clock plasmid assay. Gene subsets were also identified that exhibited increasing (103 genes) or decreasing (78 genes) transcript levels throughout (Figure 5C, Table S8 and Figure 5D, Table S9 respectively), indicative of responses to stimuli that steadily intensified or diminished during the infection.
Although transcript levels for the vast majority of genes were altered by day 2, we did identify a set of genes that exhibited a delayed “bell-curve” pattern of expression (Figure 5E, Table S10). This profile suggests that Mtb may encounter novel host-derived pressures at later time points or that prolonged growth within the phagosomal environment requires additional bacterial adaptation. This genelist included several transcriptional regulators (Rv2989, Rv2516c) and fatty acid metabolism genes (fadE9, mmsA, Rv0762c, Rv0764c-0765c). In contrast to early responses likely triggered by immediate stress, Mtb may also encounter more insidious pressures linked to evolving metabolic requirements.
Temporal Expression of Signature Regulons during the Infection
In addition to the temporal alterations in the transcriptional profile, we also examined transcription profiles of genesets linked to some of the more well-characterized environmental responses to probe how they matched the physiological transitions we have proposed occur during the 14-day infection. We have focused on the DosR regulon, regulation of carbon flux, and the response to pH, as some of the better characterized physiological themes.
Expression of the DosR dormancy regulon
Figure 6A depicts the expression profiles of the so-called “dormancy” regulon, a set of ~50 genes induced by the DosRS/T two-component system in response to multiple signals including hypoxia, nitrosative stress, and carbon monoxide [29][32]. The DosR regulon appears to play a role during transitions into and out of stress-induced dormant states [33]. Although DosR induction is typically associated with hypoxic conditions within granulomas and reactive nitrogen intermediates (RNI) derived from activated M[var phi], we noted induction of ~60% of DosR-dependent genes in resting M[var phi] within 48 hr p.i., with upregulation of some transcripts detectable as early as 2 hr p.i. (Figure 6A). A subset of DosR-dependent genes triggered by hypoxia in vitro remained unchanged in resting M[var phi] phagosomes, perhaps reflecting distinct threshold levels of sensitivity to inducing signals. It is conceivable that all three known cues of DosRS/T signaling are acting in concert within the resting phagosome. Greiss assay analysis of supernatants from 14-day infected M[var phi] showed a spike in nitrate levels at day 2 p.i. that rapidly declined to baseline levels by day 4 (data not shown), implicating nitric oxide as a signal of early DosR regulon induction. Decreased intraphagosomal oxygen concentrations in macrophages [34] and/or carbon monoxide generated by infection-induced heme oxygenase-1 (HO-1) may influence the expression of DosR-dependent genes. Finally, the expression of DosR-dependent genes may be modulated by slow growth rates and nutrient limitation during intraphagosomal growth [35].
Figure 6
Figure 6
The transcriptome as a bioprobe of the Mtb intracellular environment.
The most striking aspect of dormancy gene expression during long-term M[var phi] infection was the sudden repression of DosR-dependent genes at ~day 8, despite steady transcript levels of dosR, dosS, and dosT (Figure 6A). The downregulation to sub-control levels of some DosR-dependent genes that were not significantly induced at earlier time points suggests active repression, as well as an intricate crosstalk between DosR and other transcriptional regulators to coordinate this stress response. Taking a “guilt-by association” approach, we utilized the temporal profile of hspX, a member of the DosRS-dependent dormancy regulon [29], to identify ~50 additional early-induced genes (Figure 6B, Table S11) that exhibited the striking down-regulation at ~Day 8 noted for the DosRS regulon (see below). The coordinated regulation of poorly characterized genes with known regulons or stimulons offers a starting point for elucidating their function. We suggest that this is a highly significant transcriptional pattern because it corresponds to the time point where Mtb transitions from a period of adaptation to net intracellular growth.
This shift in Mtb-M[var phi] interactions may reflect the sudden diminution of the cue(s) driving the response. Alternatively, this profile may reflect the staged adaptation of Mtb to phagosomal stresses such as hypoxia, NO, and/or CO. Balazsi et al. reported similar induction of the DosRS regulon followed by a downward shift at ~day 8 in an in vitro hypoxia model [36]. This profile is reminiscent of the enduring hypoxic response (EHR) described by Rustad et al. (Figure S5) in which DosRS is only active during the early phase of extended exposure to hypoxia [37]. Although unable to cope with rapid shifts to hypoxic conditions, Mtb can survive gradual oxygen depletion by coordinated reprogramming of gene expression and metabolism [38][40]. Early exposure to dosRS-inducing conditions within resting M[var phi] may allow Mtb to pre-adapt for harsher conditions within activated M[var phi] and hypoxic granulomas.
Adaptation of carbon metabolism during extended intracellular growth
Sequestration of Mtb within M[var phi] phagosomes undoubtedly imposes limitations on the nutrients available for growth. Our temporal gene expression analysis revealed that realignment of carbon metabolism pathways to allow utilization of host-derived fatty acids (FA) and cholesterol is a key component of Mtb intracellular adaptation. We noted early, sustained induction of pckA (phosphoenolpyruvate carboxykinase) and concomitant repression of pca (pyruvate carboxylase) (Figure 7A). pckA is induced by fatty acids, upregulated during mouse infection [41] and required for mycobacterial survival in both M[var phi] and mice [42], [43], highlighting the importance in vivo of gluconeogenesis of carbon derived from fatty acid β-oxidation. The rapid induction of icl, encoding the isocitrate lyase component of the glyoxylate shunt [44], is consistent with the availability of fatty acid carbon sources within resting phagosomes. The lack of induction of the second gene of the glyoxylate shunt, glcB (malate synthase) (Figure 7A), suggests an alternate role for icl (discussed below). Previous studies have demonstrated the ability of Mtb to utilize cholesterol as a source of energy and carbon [45][48]. The activation of the majority of the KstR regulon by day 2 p.i., including the mce4 [47], [49] and igr [46], [50] loci, indicates that Mtb likely has ready access to cholesterol within resting M[var phi] (Figure 7B–D, Figure S6). KstR is known to be induced by cholesterol [51] and multiple KstR-dependent genes have been implicated in cholesterol uptake and metabolism [45][49], [52], [53]. In particular, the mce4 locus encodes a dedicated transport system capable of internalizing cholesterol [47], [49]. Note that a subset of genes predicted to be directly controlled by KstR exhibit unique expression patterns (lack of induction) during extended M[var phi] infection (Figure S6D). This could be due to inaccurate identification of KstR-controlled genes or, more likely, additional layers of regulatory input acting within the vacuole.
Figure 7
Figure 7
Adaptation of Carbon Metabolism during Extended Intracellular Growth.
A potentially deleterious side-effect of cholesterol catabolism and β-oxidation of odd-chain length fatty acids is the generation of toxic metabolites including propionyl-CoA and downstream intermediates [46], [54], [55]. As shown in Figure 7A, Mtb is equipped with three pathways to detoxify and utilize propionyl-CoA as either a carbon source or building block for methyl-branched fatty acid components of the cell wall: methylcitrate cycle (MCC), methylmalonyl pathway (MM), and synthesis of methyl-branched fatty acids. To identify propionyl-CoA detoxification mechanisms activated during long-term intracellular survival, we overlayed temporal expression profiles of relevant metabolic genes (shown as heat maps) onto the pathway map (Figure 7A, adapted from Savvi et al. [54]). The concerted induction of icl1 and prpCD reflects enhanced activity of the MCC within M[var phi], serving to convert propionyl-CoA into pyruvate and succinate. In addition to its role in the glyoxylate shunt pathway as an isocitrate lyase (ICL), the bifunctional icl1 is also a 2-methylisocitrate lyase (MCL) that interacts with prpC (encoding methylcitrate synthase) and prpD (encoding methylcitrate dehydratase) in the MCC [56]. Mutants lacking icl1, prpC, and/or prpD were unable to grow on propionate as a sole carbon source. In contrast to the normal growth of ΔprpCD Mtb in vivo, a Δicl1 strain was markedly attenuated in mice [55], [56]. These paradoxical phenotypes are consistent with a model invoking accumulation of prpCD-generated toxic propionyl-CoA metabolites in the absence of MCL activity rather than the inability to utilize fatty acids due to interruption of the glyoxylate shunt. In a recent study, Griffin et al. confirmed the important role for the MCC enzymes prpCD regulated via Rv1129c for cholesterol utilization during intracellular survival [57].
Other than a small increase in mutAB (methylmalonyl-CoA mutase) transcript levels, we did not detect coordinated upregulation of the MM pathway during Mtb infection of resting M[var phi]. This may be due to inactivation of the MM pathway in our infection model due to the absence of vitamin B12, a cofactor required for MutAB activity [54]. Derivatives of propionyl-CoA can also be diverted into the biosynthesis of methyl-branched fatty acids such as phthiocerol dimycocerosate (PDIM), sulfolipid (SL), and di-and poly-acyl trehaloses (DAT, PAT). Whereas PDIM and mycolic acid synthesis appeared downregulated by intracellular Mtb (Figure S7), we noted marked induction of pks2-papA1 (SL) and pks3-pks3-papA3-mmpL10 (DAT, PAT) upon internalization of Mtb by M[var phi]. We have previously linked the rapid activation of these loci within minutes of phagocytosis to acidification of the vacuole [12]. Thus, enhanced synthesis of SL, DAT, and PAT may serve as a sink to prevent accumulation of toxic byproducts of intracellular fatty acid and cholesterol metabolism.
The Acid Response regulon
In a previous study, we demonstrated that the drop in pH experienced in transition from the extracellular milieu to ~pH 6.4 in the phagosome of a resting M[var phi] serves as a dominant signal to Mtb for early intracellular gene induction [12]. Genes of this acid regulon exhibited rapid and sustained induction throughout long-term M[var phi] infection (Figure 6C). Thus, acid pH serves as an early cue of M[var phi] invasion as well as a protracted stressor during the period of slowed replication. The specific roles of acid-induced genes in directly mediating pH homeostasis and/or intracellular survival are becoming clearer. Remodeling of the cell wall in response to phagosome acidification, indicated by the induction of pks2-papA1 (sulfolipid (SL) biosynthesis) and pks3-pks4-papA3 (di- and poly-acyltrehalose (DAT, PAT) biosynthesis) (Figure 6C), suggests one potential mechanism of coping with this intracellular stressor. The concomitant repression of genes involved in phthiocerol dimycocerosate (PDIM) during M[var phi] infection (Figure S7B) further validates the finding by Jain et al. [58] that SL and PDIM production are controlled by the availability of the common methyl malonyl CoA (MMCoA) precursor, and are linked to the bacterium's desire to regulate the propionyl CoA pool and limit its potential toxicity, as detailed in the previous section.
WhiB3, an Actinomycete-specific protein shown to modulate the incorporation of propionate into virulence lipids including SL and PDIM in response to redox signals [59], was also induced throughout the time-course (Figure 6C). The responsiveness of whiB3 to phagosomal acidification [12], [60] suggests a link between acid stress, internal redox balance, and virulence lipid biosynthesis. The acid responsive locus aprABC (MT2466-MT2467-Rv2396), which exhibits early sustained induction, mediates intracellular survival via unknown mechanisms by influencing gene expression and lipid metabolism [21]. Finally, Mtb may rely on phagosomal acidification as a cue of M[var phi] invasion to coordinate the timely expression of intracellular virulence mechanisms. For example, acid induction of espABC (Rv3614c–Rv3616c), a locus required for ESX-1 mediated secretion of virulence factors [61], [62], may facilitate the rapid secretion of preformed Type VII effectors into the phagosome (Figure 6C, Figure S8B). We also noted the induction of ESX-1 components encoded within the RD1 locus including EspR (Rv3849), a secreted regulator of espABC (Figure S8) [62].
Transcriptional Regulatory Network Response to Phagosome Cues
In addition to analyzing the transcriptional profiles according to temporal dynamics and known functional themes we also conducted a systems-level analysis to characterize the behavior of the Mtb Transcriptional Regulatory Network (TRN) underlying pathogen survival during the 14-day infection. This network-based approach incorporates extensive a priori information on Mtb gene regulation and network topology, combined with expression data, to assess network responses in surviving bacteria elicited by the M[var phi].
We started by expanding a large-scale Mtb TRN containing gene regulatory interactions extracted from both experimental and computational datasets [37]. The previous TRN comprised 738 genes (18% of the genome) and 937 regulatory links obtained from three sources: (i) literature mining; (ii) MtbRegList database [63]; (iii) inference from orthology with Escherichia coli [64]. To obtain the expanded network used in this study, we collected gene regulatory data from the following additional sources. (iv) We added orthology-based interactions inferred from the closely related Corynebacterium glutamicum available in the MycoRegNet database, which considerably expanded the TRN by adding 425 new interactions (Figure S9A). The extensive overlap with regulatory links supported by experimental data validated these interactions identified from orthologous Transcription Factor (TF)-Target Gene (TG) pairs in the two organisms (p = 10−5, Fisher's exact test, Figure S9A). (v) Further, we expanded the TRN by adding 114 protein-DNA interactions discovered by a new bacterial one-hybrid reporter system termed TB1Hybrid [65], being that 31 interactions were exclusively identified by this method. (vi) Finally, we performed operon-based network expansion, propagating a TF's regulatory effect to all members of the operon containing a given TG. The final expanded Mtb TRN contained 1133 genes (28% of the genome) and 1801 regulatory links, more than a half of which were experimentally determined (Figure S9B, C); the complete list of interactions is available in the Table S12.
The global TRN provides a static summary of all possible regulatory interactions that mycobacteria may use when facing a broad spectrum of environments, ranging from normal to stressful conditions inside the M[var phi]. However, previous work has suggested that only parts of the network are utilized in specific conditions [36]. Such parts (subnetworks) function as network modules regulated by a hierarchy of transcription factors in an environment-dependent fashion. To understand how the surviving subset of Mtb bacilli specifically utilizes the TRN modules during prolonged intracellular infection, we analyzed the temporal response of the TRN by overlaying the 14-day M[var phi] infection time-course array data on the extended Mtb regulatory network.
We improved the earlier method called NetReSFun (Network Response to Step Functions), which identifies responsive TF-regulated subnetworks from time course microarray data [36]. NetReSFun computes the Cov-score (Methods) to quantify the expression change within the module (subnetwork) that we define as the total genes directly regulated by a given TF. A significant Cov-score indicates subnetwork response, when the expression levels of the subnetwork's gene members are either down- or upregulated during consecutive time points (t, t+1). Alternatively, simultaneous change of a TF's direct target genes may also be a surrogate of the TF's activity. Under this assumption, NetReSFun may recognize TFs that are “turned on” through posttranslational modifications such as phosphorylation and metabolite binding, but may or may not show increased expression levels themselves.
The temporal map of network responses (Figure 8A) depicts specific TF-regulated subnetworks responsive during the time course, at a significance level of 0.05. The color scale indicates whether the overall trend of expression change within the subnetwork was positive (red) or negative (blue) at a given time interval. Intermediate colors denote subnetworks involving both up- and downregulated genes. The accompanying heatmap in Figure 8B indicates the source of regulatory links within the subnetwork (darker colors corresponding to higher fractions of links based on experimental evidence).
Figure 8
Figure 8
Temporal network response during macrophage infection.
Strikingly, the map reveals that the dynamic utilization of the TRN occurs in a defined pattern that can be mapped back to distinctive phases of intracellular growth. In the first 2 days p.i. - which corresponds to the stress phase in the growth curve of Figure 1A – we observed a high number of responsive TF-subnetworks (20 out of 83), mostly exhibiting increased expression of involved genes. Among these were DosR, HspR, KstR, members of the WhiB family (WhiB3, WhiB4), two-component response regulators (Rv0260c, Rv0818, RegX3), and alternative sigma factors (SigE, SigK, SigM). The sharp induction of a large number of subnetworks is reminiscent of the general Environmental Stress Response (ESR) in yeast [66] and in Bacillus subtilis [67].
In contrast, after ~6 days inside M[var phi], we observed a reciprocal scenario where the TRN reflects a significant downmodulation of target genes, many of which had been induced immediately after invasion. This pattern is especially evident for a number of stress-responsive subnetworks that shift into downregulation during the slow growth phase, including RegX3, HspR and DosR. This “repressive” transcription phase indicates that the surviving bacteria have either adapted to stress, or they reside in a less hostile niche. For example, the SigH-controlled subnetwork displays a strong negative response only late in the time course (~8 days and onwards). As a global regulator, SigH modulates the transcription of SigE and SigB, as well as its own promoter. Although SigH is not required for growth in M[var phi], mutants lacking sigH caused reduced immunopathology and lethality in mice [68], [69]. Alternatively, it is possible that these regulatory changes are associated with the surviving bacilli reprogramming their physiology to assume the slow growth phenotype observed from day 4 onwards.
Importantly, we observed responsive subnetworks throughout the entire time-course, which indicates their importance for establishing productive infection. The presence of sustained responders such as HspR and DosR can have two possible implications. First, the opposite trends in the early phase of the infection (primarily upregulation) and later phases of infection (primarily downregulation) may indicate that the stress to which these modules respond initially is ameliorated at later time points. Alternatively, sustained responders may be necessary both to counteract initial phagosomal stress during the early phase of infection as well as for driving the persistor phenotype encountered in later phases of infection. This is the case of DosR, which is crucial for maintaining redox balance and energy levels during transitions into and out of dormancy-like conditions that perturb aerobic respiration, electron transport, or menaquinone pools [70]. HspR, which activates a subset of the heat-shock general stress response upon M[var phi] invasion [71], is also necessary in the persistent phase since ΔhspR strains exhibited attenuated growth in the chronic infection [72].
Finally, TRN analyses revealed novel sustained responders that might be critical for Mtb adaptation within the intracellular compartment. For example, the Rv2034-controlled subnetwork (inferred from C. glutamicum orthology) contains multiple fadE homologs implicated in β-oxidation of fatty acids and redox homeostasis. Notably, Rv2034 was recently characterized as an activator of the phoP virulence regulator in mycobacteria [73], which makes this regulator an interesting candidate for follow-up studies.
By overlaying the genome-scale temporal expression data onto the TRN, we revealed the activation of additional regulons during macrophage survival not readily apparent in our supervised analyses. Thus, by leveraging the behavior of multiple TG as a readout of TF activity, TRN analysis of time-course microarray data further enhances the ability to detect adaptive changes in Mtb gene expression during productive infection of macrophages.
Conclusions
The process of infection is extremely dynamic as both host and pathogen seek to respond to the stimuli that they sense at their interface. In the current study we applied multiple analytical tools to establish a link between the transcriptional responses and physiological states through which Mtb transitions on its way to the establishment of a productive infection in its host macrophage. This analysis revealed several unexpected findings. Firstly, the initial phase of infection is marked by rapid bacterial replication coupled with effective bacterial killing by the macrophage. This is a period of marked stress for Mtb, which is illustrated by the greatest transcriptional response with respect to both up-regulated and down-regulated genes. Subsequently, the rate of replication slows and the bacterial number appears constant or at equilibrium, during which period the expression of many genes returns closer to control levels, whilst the divergent level of expression of others is sustained. Finally, the bacterial numbers start to increase indicating that the rate of replication exceeds that of death. At this time there is a marked down-regulation of many of the genes linked to general stress supporting the contention that Mtb has entered into a productive phase of infection characterized by enhanced intracellular survival. This is an important functional framework on which to hang the transcriptional profiles to determine which responses and/or metabolic themes impact which phase of infection.
Our examination of dynamic alternations of gene expression across the Mtb TRN also highlights stress responses and survival mechanisms deployed during distinct phases of host interaction. Among other things, our results indicate adaptive changes in lipid and energy metabolism akin to those observed in various dormancy models. Based on this, it is tempting to speculate that early infection of resting macrophages may serve to prepare Mtb for conditions encountered within granulomas after the onset of the adaptive immune response. Detailed understanding of the sensory and regulatory pathways required for Mtb virulence remains rudimentary at best. This point is illustrated by the large number of “genes of unknown function” which are actively regulated during intracellular survival.
In the discussion, we have presented the observed transcriptomic changes as the result of all Mtb cells sensing and responding in concert to phagosomal cues in the intracellular environment. Clearly, the microarray profiles in this study capture an average behavior over time of a population of intracellular bacilli that, as we have shown, exist in distinct vacuolar niches and presumably metabolic states. It could be argued that the changes in gene expression represent the minority of Mtb that fail to block P-L fusion mounting dramatic stress responses. However, our data appear to be inconsistent with this view. For example, we know that Mtb in P-L encounter Fe-deplete conditions whereas the average behavior of Fe-responsive signature genes reflects an Fe-replete environment (data not shown). Alternatively, the apparent up- and down-regulation may be due to the selective enrichment of pre-existing phenotypic variant cells with randomly upregulated stress response subnetworks, similar to the cell-to-cell heterogeneity reported by Aldridge et al. [74]. The rapid killing observed immediately after macrophage entry supports this scenario, suggesting that a subpopulation of bacilli capable of surviving the bottleneck may be present prior to infection, due to random phenotypic variation within the microbial population. However, the ability to block the induction of an acid regulon following M[var phi] invasion by chemical manipulation of phagosomal pH [12] would suggest that Mtb are altering gene expression in response to host-derived cues.
Thus, our data suggest that the dominant transcriptional profiles highlighted here represent the adaptive responses of the majority population. Analysis at the single-cell level will be required to explore the strategies employed by Mtb to survive across diverse host environments experienced by each individual bacterium. To address this, we have begun to exploit fluorescent reporter strains responsive to specific environmental cues to gain a high-resolution view of Mtb intracellular adaptation and cell biology. Interestingly, the expression of an acid-inducible locus required for normal intracellular survival, aprABC, is expressed at distinct levels by individual bacilli within the same host cell [21]. The impact of this type of heterogeneity on Mtb-host interactions and pathogenesis has yet to be determined. Previous work has shown that stochastic variation (intrinsic to individual cells) can aid survival in stressful environmental conditions [75]. However, a similar role of extrinsic, environmental variation is yet to be established.
The identification of mutants attenuated for intracellular survival is a popular and powerful tool for defining in vivo survival mechanisms, however, there are some limitations to this methodology that can be addressed by transcriptional profiling. Array-based mutant screens such as Transposon Site Hybridization (TraSH) [76], [77] are in essence end-point assays that are not readily amenable to quantitative kinetic analyses. In addition, there are several well-characterized examples of genes known to be required for M[var phi] survival, such as pckA [42], [43], icl ([44], [78], prpCD [56], katG [79], [80], or phoPR [9], [81], that have not been identified by TraSH screens indicating that the method is not comprehensive [82]. The TraSH methodology may preferentially identify mutants with severe survival defects, while being less effective at isolating mutants with less extreme phenotypes, such as ΔphoPR, that resist killing but fail to grow within phagosomes [9]. Finally, mutant screens are limited by their inability to query the in vivo role of genes that are essential in vitro or identify genes whose phenotype upon inactivation is masked by compensatory changes in gene expression. Future studies would benefit from the coordinated application of these two distinct but complementary approaches to identify genes contributing to the pathogenesis of Mtb.
We feel that the significance of this current study is that it transforms transcriptional profiling from a purely descriptive analysis to the generation of a predictive discovery tool that can be used to identify genes, and therefore metabolic pathways and physiological states, that are required to support distinct phases in the intracellular life cycle of M. tuberculosis.
Bacterial Strains and Cells
As previously described [83], the M. tuberculosis clinical isolate strain CDC1551 [84], [85] was cultured in Middlebrook 7H9-OADC medium in ventilated T-75 flasks without shaking. Bone marrow-derived macrophages (M[var phi]) were isolated from C57BL/6 mice and grown in DMEM supplemented with 20% L-cell conditioned medium, 10% fetal calf serum (FCS), and antibiotics (penicillin and streptomycin). Media lacking antibiotics was added to M[var phi] at least 24 hr prior to infection with M. tuberculosis.
Survival Assays
To monitor the survival/growth of Mtb in M[var phi], confluent M[var phi] monolayers in 24-well dishes were infected (MOI ~1[ratio]1) with Mtb CDC1551. After 2 hr (t = 0) and at two day intervals up to 14 days, intracellular Mtb released from monolayers by lysis with ddH2O+0.05% Tween-80 were serially diluted, and plated on 7H10+cycloheximide agar. CFU were enumerated after ~3 weeks incubation at 37°C. The overall integrity of M[var phi] monolayer was verified by microscopy.
Replication Clock Plasmid
In order to measure the dynamics of Mtb replication and death during M[var phi] infection, a replication clock plasmid, pBP10, was used (generous gift of Dr. David R. Sherman). As previously described [15], in the absence of kanamycin selection pBP10 is lost at a rate proportional to the growth rate. Macrophage survival assays using Mtb CDC1551 harboring pBP10 (introduced by electroporation) were conducted as described above except that aliquots were plated on 7H10+cycloheximide agar both with and without 25 µg/ml kanamycin. Based on the loss of kanamycin resistance over time, rates of replication and death were quantified using the mathematical model of Gill et al. [15].
Electron Microscopy
In parallel with M[var phi]-Mtb infections for survival assays and microarray analysis, confluent monolayers of M[var phi] were infected with Mtb CDC1551 at a low MOI (1[ratio]1) as described above. At select times p.i. (2 hr and day 2-day 14, alternating days), samples were fixed in buffered glutaraldehyde solution (2.5% glutaraldehyde in 0.1 M sodium cacodylate, 5 mM CaCl2, 5 mM MgCl2, 0.1 M sucrose, pH 7.2), rinsed with 0.1M sodium cacodylate buffer (see above), post-fixed with 1% osmium tetroxide (4% stock osmium diluted in 0.1 M sodium cacodylate, 5 mM CaCl2, 5 mM MgCl2, pH 7.2), rinsed with buffer, soaked in 1% aqueous uranyl acetate, and then rinsed with water. Samples were then dehydrated in a graded ethanol and propylene oxide series, followed by gradual infiltration of Spurr's resin. Blocks were then polymerized and ultrathin sections (~70 nm) were cut and contrasted with both lead citrate and uranyl acetate. In independent experiments to determine phagosome-lysosome fusion by the colocalization of Mtb with colloidal gold, at various time-points after infection (2 hr, day 2, 6, and 10) infected monolayers were pulsed with 15 nm colloidal gold (Aurion) for 2 hr, washed, and then chased for 45 min in infection media before fixation as above. In assessment of the morphology of intracellular Mtb, bacilli were considered intact if they maintained a rod shape (longitudinal sections) or circular shape (cross sections), if ultrastructural organization and electron opacity of the cytoplasm was preserved, and if no breaks in the cytoplasmic membrane or cell wall were detected. Otherwise, bacilli were counted as damaged.
Macrophage Infection, RNA Isolation and Linear Amplification
The methodology for transcriptional profiling of intramacrophage Mtb including RNA isolation, linear amplification, and hybridization has been described previously [12], [83]. Briefly, C57BL/6 bone marrow-derived M[var phi] were infected (MOI ~10[ratio]1) with M. tuberculosis CDC1551 from 2 hr up to 14 days (2 day intervals). In controls, aliquots of the same bacterial samples were incubated in flasks without a M[var phi] monolayer for 2 hr. Addition of guanidine thiocyanate-based lysis buffer selectively lysed M[var phi] and stabilized bacterial RNA while leaving mycobacteria intact. Pelleted bacilli were lysed in 65°C Trizol using a BeadBeater and 0.1 mm silicon beads. Total RNA was isolated from Trizol lysates by chloroform extraction and Qiagen RNeasy column purification. To generate array targets, 250 ng of total RNA was amplified using the MessageAmp-II Bacteria RNA Amplification system (Ambion). Amino-allyl UTP was incorporated into aRNA during transcription to allow labeling with Alexa dyes.
Target Preparation and Microarray Hybridization
Amino-allyl modified aRNA were labeled with Alexa Fluor 555 and Alexa Fluor 647 (Invitrogen) and purified using a MegaClear kit (Ambion). 10 µg of Alexa-labeled aRNA from paired samples was dried and resuspended in 75 µl of hybridization buffer (5× SSC, 25% formamide, 0.1% SDS, and 25 µg salmon sperm DNA). Slides were prehybridized for 1 hr in 25% formamide, 5× SSC, 0.1% SDS, 1% BSA and washed with H2O and isopropanol. Labeled targets were denatured at 95°C for 5 min, cooled to 60°C, and hybridized to microarrays at 45°C for 16–18 hr. Following hybridization, arrays were washed and processed for scanning as described previously [83]. The microarray platform used can be accessed via NCBI's Gene Expression Omnibus (GEO) database [86] under platform accession number GPL5754. This dataset has been deposited in GEO under series accession number GSE35362.
Microarray Data Analysis
Microarrays were scanned with a GenePix 4000B instrument (Axon Instruments, Inc.) with preliminary image analysis, spot intensity determination, background measurements, spot quality assessment and flagging conducted using Imagene software (version 6.0, Biodiscovery). Poor quality spots with signal intensities less than three standard deviations above background were excluded from further analysis. Genes that were not flagged as Present in at least 14 of 16 slides were omitted from further analyses. Subsequent normalization, statistical analysis, and visualization of array data were performed with Genespring 7.3 (Agilent). We utilized the EDGE (Extraction and analysis of Differential Gene Expression) methodology of Storey et al. to identify time-dependent transcriptional changes [22], allowing detection of genes whose expression exhibited significant temporal trends but failed to meet static statistical cutoffs at any single time point. Genes with significant changes in expression levels relative to controls were identified based on both static (p<0.05 for at least one time point) and EDGE analysis (cubic spline of 4, q<0.03).
Quantitative Real-Time RT-PCR (qRT-PCR)
RNA amplification and microarray methodology used in this study have been previously validated by qRT-PCR [12]. Additional qRT-PCR validation of temporal expression patterns of select genes during long-term M[var phi] infection was conducted by two-step real-time RT-PCR using iScript and iTaq SYBR Green reagents (Biorad). Each sample was analyzed in triplicate on an ABI 7500 starting with 100 ng of total RNA (amplified and unamplified). CT values were normalized to values obtained for sigA, a constitutively expressed Mtb gene, and relative changes in gene expression were calculated using the 2−ΔΔCT method [87].
Subnetwork Response Identification Using NetReSFun
We modified the tool NetReSFun (89) in order to identify TF-regulated subnetworks sequentially responsive during TB lifespan in the macrophage environment. Consider a set of t step functions each of which “jumps” from 0 to 1 at time point τ:
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e001.jpg
(1)
NetReSFun compares gene expression profiles to these pre-defined step functions using scaled covariance to detect whether a gene's expression changed in the interval [τ-1, τ]. The method starts by computing the scaled covariance An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e002.jpg between the expression (log10-ratio) profile An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e003.jpg of gene i and a step function An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e004.jpg that jumps at time point τ:
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e005.jpg
(2)
where brackets denote averaging over genes, horizontal bar is averaging over time, and σ is the standard deviation. Thus, An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e006.jpg is the response of gene i during interval [τ−1, τ]. The combined response of subnetwork I, or Cov-score, is the mean of the absolute covariances of subnetwork genes
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e007.jpg
(3)
To assess statistical significance, we compare the subnetwork response An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e008.jpg with a reference cumulative density function (c.d.f) An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e009.jpg, constructed with scores drawn of 1,000 random subnetworks (same size of subnetwork I, but assembled using nodes randomly chosen from the network). The density estimation of the c.d.f is done in a standard fashion by assigning probability mass 1/1000 for each observed random score; linear interpolation is applied to transform the discrete An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e010.jpg into a continuous function in the interval [0,1]. Finally, a p-value of observing An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e011.jpg by chance can be readily estimated as below, and a p-value≤0.05 was considered significant.
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e012.jpg
(4)
To determine the direction of transcriptional change displayed in the temporal map (Figure 8), we first compute the deviation index of the subnetwork as the positive/negative ratio of individual covariances
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e013.jpg
(5)
The z-score of the deviation index, which relates the change of a particular subnetwork I to all other subnetworks, is given by
A mathematical equation, expression, or formula.
 Object name is ppat.1002769.e014.jpg
(6)
where A denotes all TF-regulated subnetworks, brackets denote averaging over all subnetwork deviations, and σ is the standard deviation. Thus, the deviation z-score An external file that holds a picture, illustration, etc.
Object name is ppat.1002769.e015.jpg is the quantity plotted in the temporal map of Figure 8.
Figure S1
Quantitative electron microscopy analysis of long-term Mtb-macrophage interactions. At each timepoint, all visible intracellular Mtb in 100 macrophages were enumerated and categorized based on morphological criteria as described in Materials and Methods. (A) Morphology of intracellular Mtb. The proportion of intact, morphologically normal bacilli to damaged Mtb remained relatively constant over the 14-day infection. (B) Change in distribution of Mtb burden over time. The number of Mtb per macrophage increases steadily over time. Note, the large set of cells in which no Mtb were detected at early time points likely reflects the absence of detectable bacilli within the plane of section rather than uninfected cells.
(TIF)
Figure S2
Mtb occupies heterogeneous intracellular niches during macrophage infection. (A) Electron microscopy image shows dividing bacilli in “tight” vacuole, day 10 p.i. (B) Morphologically intact Mtb in spacious lysosome-like compartments surrounded by granular debris at day 6 p.i. This suggests Mtb can survive at least for a time in vacuoles that have fused with lysosomes. (C) Mtb surrounded by double-membrane (arrow) vacuole containing colloidal gold, consistent with an autophagosome (day 14 p.i.).
(TIF)
Figure S3
Temporal expression profiling reveals novel genes differentially regulated within macrophage phagosomes. (A) Gene tree (Euclidean distance measure) of 137 genes significantly upregulated at 48 hr p.i. novel to this study (compared to expression profiles at same time point reported in [26]). (B) Differential expression of “MT genes” during long-term macrophage infection. This geneset includes 292 predicted ORFs annotated in strain CDC1551 genome [88] not originally annotated in H37Rv genome [89] significantly regulated during long-term macrophage infection.
(TIF)
Figure S4
Long-term intracellular expression of the core transcriptome. Gene tree showing the distinct expression profiles of 215 genes comprising a core intracellular transcriptome previously defined based on their conserved presence and induction across a diverse panel of Mtb clinical isolates at 24 hr p.i. of resting macrophages [13].
(TIF)
Figure S5
Genes of the Enduring Hypoxic Response (EHR) exhibit distinct expression patterns during long-term adaptation within macrophages. Gene tree showing the transcriptional patterns of genes comprising the EHR, as identified by Rustad et al. [37], during long-term macrophage infection.
(TIF)
Figure S6
Selective up-regulation of members of the KstR-dependent cholesterol regulon during intracellular growth. (A) Gene tree showing the relative transcript levels of genes proposed to be directly controlled by KstR (Kendall et al. [51]) during long-term adaptation within macrophage phagosomes. Expanded views of distinct gene clusters are shown in (B–E). (B) Genes with early induction to moderate levels, most sustained throughout. (C) Genes with early induction to higher levels, with decrease at ~day12 p.i., (D) KstR-dependent genes minimally responsive to phagosomal cues. (E) Genes displaying sustained, high level induction.
(TIF)
Figure S7
Regulation of cell wall synthesis and composition in response to phagosomal cues. Gene trees showing the sustained downregulation of genes involved in mycolic acid synthesis (A) and pthiocerol dimycocerosate (PDIM) (B).
(TIF)
Figure S8
Regulation of ESX secretion systems during intracellular growth. (A,B) Early induction of the ESX-1 secretion system including the main ESX-1 locus (A) as well as the EspR regulator and accessory factors Rv3614c-3616c encoded outside the RD-1 locus (B). These gene products function coordinately to facilitate secretion of ESAT-6/CFP-10 complex. [90], [91].
(TIF)
Figure S9
The enlarged Mtb Transcriptional Regulatory Network (TRN). (A) Number of regulatory links from each data source and corresponding overlaps. Links with experimental evidence originate from literature, as well as from the MtbRegList database and the TB1H assay. (B) Distribution of interactions based on their inference method. (C) Overview of the TRN, depicting protein-DNA interactions as edges linking TFs (blue triangles) to TGs (yellow circles).
(TIF)
Table S1
Genes with peak induction at Day 2 followed by gradual decline toward control levels.
(XLSX)
Table S2
Genes significantly upregulated (>1.5X) at 48 hr novel to this study (compared to expression of >1.5X in both amplicon and oligo array data reported by Schnappinger et. al (ref 29).
(XLSX)
Table S3
Temporal expression profile of “MT genes”, predicted ORFs annotated in strain CDC1551 genome not originally annotated in the H37Rv genome. Note, the gene order does not match the accompanying gene tree in Fig. S3B.
(XLSX)
Table S4
Figure 5A - Genes with Early, transient induced temporal profile.
(XLSX)
Table S5
Figure 5B (red) - Genes with Early Induced, Sustained temporal profile.
(XLSX)
Table S6
Subsets of previously identified core intracellular transcriptome display distinct temporal expression profiles during 14-day macrophage infection.
(XLSX)
Table S7
Figure 5B (black) - Genes with Early Repressed, Sustained temporal profile.
(XLSX)
Table S8
Figure 5C -Genes with Steady, sustained induced temporal profile.
(XLSX)
Table S9
Figure 5D -Genes with Steady, sustained repression temporal profile.
(XLSX)
Table S10
Figure 5E - Genes with Delayed induced temporal profile.
(XLSX)
Table S11
Figure 6B - Genes with hspX-like temporal profiles.
(XLSX)
Table S12
M. tuberculosis Expanded Transcription Regulatory Network (TRN).
(XLSX)
Acknowledgments
We would like to thank Paul Debbie (Center for Gene Expression Profiling [CGEP], Boyce Thompson Institute, Ithaca, NY) and Dr. Mark Band (Functional Genomics Unit, Roy J. Carver Biotechnology Center, University of Illinois at Urbana-Champaign) for array printing and Dr. David Lin (Cornell Big Red Spots Core Facility, Cornell University) for array scanning. We are grateful to T.B.K. Reddy and Janos Demeter (Stanford University) for submission of microarray data to the TBDB and GEO databases. We also thank Dr. David Sherman (Seattle Biomedical Research Institute, Seattle, WA) for the generous gift of the pBP10 replication clock plasmid.
Footnotes
The authors have declared that no competing interests exist.
This work was supported by U.S. National Institutes of Health (NIH) grant AI067027 to D.G.R., European Community Seventh Framework Programme (FP7-223681), NIH Director's New Innovator Award Program (1DP2 OD006481-01), and NHLBI TB latency consortium (1-RO1-HL106788) to G.B.. D.F.T.V. was supported by a CAPES/Fulbright doctoral scholarship (Brazilian Ministry of Education). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. de Chastellier C. The many niches and strategies used by pathogenic mycobacteria for survival within host macrophages. Immunobiology. 2009;214:526–542. [PubMed]
2. Deretic V, Delgado M, Vergne I, Master S, De Haro S, et al. Autophagy in immunity against mycobacterium tuberculosis: a model system to dissect immunological roles of autophagy. Curr Top Microbiol Immunol. 2009;335:169–188. [PMC free article] [PubMed]
3. Roberts EA, Deretic V. The Mycobacterium tuberculosis phagosome. Methods Mol Biol. 2008;445:439–449. [PubMed]
4. Russell DG. Mycobacterium tuberculosis and the intimate discourse of a chronic infection. Immunol Rev. 2011;240:252–268. [PMC free article] [PubMed]
5. Russell DG, VanderVen BC, Lee W, Abramovitch RB, Kim MJ, et al. Mycobacterium tuberculosis wears what it eats. Cell Host Microbe. 2010;8:68–76. [PMC free article] [PubMed]
6. van der Wel N, Hava D, Houben D, Fluitsma D, van Zon M, et al. M. tuberculosis and M. leprae translocate from the phagolysosome to the cytosol in myeloid cells. Cell. 2007;129:1287–1298. [PubMed]
7. Simeone R, Bobard A, Lippmann J, Bitter W, Majlessi L, et al. Phagosomal Rupture by Mycobacterium tuberculosis Results in Toxicity and Host Cell Death. PLoS Pathog. 2012;8:e1002507. [PMC free article] [PubMed]
8. Converse PJ, Karakousis PC, Klinkenberg LG, Kesavan AK, Ly LH, et al. Role of the dosR-dosS two-component regulatory system in Mycobacterium tuberculosis virulence in three animal models. Infect Immun. 2009;77:1230–1237. [PMC free article] [PubMed]
9. Perez E, Samper S, Bordas Y, Guilhot C, Gicquel B, et al. An essential role for phoP in Mycobacterium tuberculosis virulence. Mol Microbiol. 2001;41:179–187. [PubMed]
10. Steyn AJ, Collins DM, Hondalus MK, Jacobs WR, Jr, Kawakami RP, et al. Mycobacterium tuberculosis WhiB3 interacts with RpoV to affect host survival but is dispensable for in vivo growth. Proc Natl Acad Sci U S A. 2002;99:3147–3152. [PubMed]
11. Zahrt TC, Deretic V. Mycobacterium tuberculosis signal transduction system required for persistent infections. Proc Natl Acad Sci U S A. 2001;98:12706–12711. [PubMed]
12. Rohde KH, Abramovitch RB, Russell DG. Mycobacterium tuberculosis invasion of macrophages: linking bacterial gene expression to environmental cues. Cell Host Microbe. 2007;2:352–364. [PubMed]
13. Homolka S, Niemann S, Russell DG, Rohde KH. Functional genetic diversity among Mycobacterium tuberculosis complex clinical isolates: delineation of conserved core and lineage-specific transcriptomes during intracellular survival. PLoS Pathog. 2010;6:e1000988. [PMC free article] [PubMed]
14. Pethe K, Swenson DL, Alonso S, Anderson J, Wang C, et al. Isolation of Mycobacterium tuberculosis mutants defective in the arrest of phagosome maturation. Proc Natl Acad Sci U S A. 2004;101:13642–13647. [PubMed]
15. Gill WP, Harik NS, Whiddon MR, Liao RP, Mittler JE, et al. A replication clock for Mycobacterium tuberculosis. Nat Med. 2009;15:211–214. [PMC free article] [PubMed]
16. Segovia-Juarez JL, Ganguli S, Kirschner D. Identifying control mechanisms of granuloma formation during M. tuberculosis infection using an agent-based model. J Theor Biol. 2004;231:357–376. [PubMed]
17. Armstrong JA, Hart PD. Response of cultured macrophages to Mycobacterium tuberculosis, with observations on fusion of lysosomes with phagosomes. J Exp Med. 1971;134:713–740. [PMC free article] [PubMed]
18. Armstrong JA, Hart PD. Phagosome-lysosome interactions in cultured macrophages infected with virulent tubercle bacilli. Reversal of the usual nonfusion pattern and observations on bacterial survival. J Exp Med. 1975;142:1–16. [PMC free article] [PubMed]
19. Xu S, Cooper A, Sturgill-Koszycki S, van Heyningen T, Chatterjee D, et al. Intracellular trafficking in Mycobacterium tuberculosis and Mycobacterium avium-infected macrophages. J Immunol. 1994;153:2568–2578. [PubMed]
20. de Chastellier C, Forquet F, Gordon A, Thilo L. Mycobacterium requires an all-around closely apposing phagosome membrane to maintain the maturation block and this apposition is re-established when it rescues itself from phagolysosomes. Cell Microbiol. 2009;11:1190–1207. [PubMed]
21. Abramovitch RB, Rohde KH, Hsu FF, Russell DG. aprABC: a Mycobacterium tuberculosis complex-specific locus that modulates pH-driven adaptation to the macrophage phagosome. Mol Microbiol. 2011;80:678–94. [PMC free article] [PubMed]
22. Storey JD, Xiao W, Leek JT, Tompkins RG, Davis RW. Significance analysis of time course microarray experiments. Proc Natl Acad Sci U S A. 2005;102:12837–12842. [PubMed]
23. Ingraham JL, Maaloe O, Neidhardt FC. The Growth of the Bacterial Cell. Sunderland MA: Sinauer Associates 1983
24. Selinger DW, Saxena RM, Cheung KJ, Church GM, Rosenow C. Global RNA half-life analysis in Escherichia coli reveals positional patterns of transcript degradation. Genome Res. 2003;13:216–223. [PubMed]
25. Micklinghoff JC, Breitinger KJ, Schmidt M, Geffers R, Eikmanns BJ, et al. Role of the transcriptional regulator RamB (Rv0465c) in the control of the glyoxylate cycle in Mycobacterium tuberculosis. J Bacteriol. 2009;191:7260–7269. [PMC free article] [PubMed]
26. Schnappinger D, Ehrt S, Voskuil MI, Liu Y, Mangan JA, et al. Transcriptional Adaptation of Mycobacterium tuberculosis within Macrophages: Insights into the Phagosomal Environment. J Exp Med. 2003;198:693–704. [PMC free article] [PubMed]
27. Walters SB, Dubnau E, Kolesnikova I, Laval F, Daffe M, et al. The Mycobacterium tuberculosis PhoPR two-component system regulates genes essential for virulence and complex lipid biosynthesis. Mol Microbiol. 2006;60:312–330. [PubMed]
28. Morris RP, Nguyen L, Gatfield J, Visconti K, Nguyen K, et al. Ancestral antibiotic resistance in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2005;102:12200–12205. [PubMed]
29. Park HD, Guinn KM, Harrell MI, Liao R, Voskuil MI, et al. Rv3133c/dosR is a transcription factor that mediates the hypoxic response of Mycobacterium tuberculosis. Mol Microbiol. 2003;48:833–843. [PMC free article] [PubMed]
30. Kendall SL, Movahedzadeh F, Rison SC, Wernisch L, Parish T, et al. The Mycobacterium tuberculosis dosRS two-component system is induced by multiple stresses. Tuberculosis (Edinb) 2004;84:247–255. [PubMed]
31. Roberts DM, Liao RP, Wisedchaisri G, Hol WG, Sherman DR. Two sensor kinases contribute to the hypoxic response of Mycobacterium tuberculosis. J Biol Chem. 2004;279:23082–23087. [PMC free article] [PubMed]
32. Kumar A, Deshane JS, Crossman DK, Bolisetty S, Yan BS, et al. Heme oxygenase-1-derived carbon monoxide induces the Mycobacterium tuberculosis dormancy regulon. J Biol Chem. 2008;283:18032–18039. [PubMed]
33. Leistikow RL, Morton RA, Bartek IL, Frimpong I, Wagner K, et al. The Mycobacterium tuberculosis DosR regulon assists in metabolic homeostasis and enables rapid recovery from nonrespiring dormancy. J Bacteriol. 2010;192:1662–1670. [PMC free article] [PubMed]
34. James PE, Grinberg OY, Michaels G, Swartz HM. Intraphagosomal oxygen in stimulated macrophages. J Cell Physiol. 1995;163:241–247. [PubMed]
35. Beste DJ, Laing E, Bonde B, Avignone-Rossa C, Bushell ME, et al. Transcriptomic analysis identifies growth rate modulation as a component of the adaptation of mycobacteria to survival inside the macrophage. J Bacteriol. 2007;189:3969–3976. [PMC free article] [PubMed]
36. Balazsi G, Heath AP, Shi L, Gennaro ML. The temporal response of the Mycobacterium tuberculosis gene regulatory network during growth arrest. Mol Syst Biol. 2008;4:225. [PMC free article] [PubMed]
37. Rustad TR, Harrell MI, Liao R, Sherman DR. The enduring hypoxic response of Mycobacterium tuberculosis. PLoS One. 2008;3:e1502. [PMC free article] [PubMed]
38. Rustad TR, Sherrid AM, Minch KJ, Sherman DR. Hypoxia: a window into Mycobacterium tuberculosis latency. Cell Microbiol. 2009;11:1151–1159. [PubMed]
39. Voskuil MI. Mycobacterium tuberculosis gene expression during environmental conditions associated with latency. Tuberculosis (Edinb) 2004;84:138–143. [PubMed]
40. Wayne LG, Hayes LG. An in vitro model for sequential study of shiftdown of Mycobacterium tuberculosis through two stages of nonreplicating persistence. Infect Immun. 1996;64:2062–2069. [PMC free article] [PubMed]
41. Timm J, Post FA, Bekker LG, Walther GB, Wainwright HC, et al. Differential expression of iron-, carbon-, and oxygen-responsive mycobacterial genes in the lungs of chronically infected mice and tuberculosis patients. Proc Natl Acad Sci U S A. 2003;100:14321–14326. [PubMed]
42. Liu K, Yu J, Russell DG. pckA-deficient Mycobacterium bovis BCG shows attenuated virulence in mice and in macrophages. Microbiology. 2003;149:1829–1835. [PubMed]
43. Marrero J, Rhee KY, Schnappinger D, Pethe K, Ehrt S. Gluconeogenic carbon flow of tricarboxylic acid cycle intermediates is critical for Mycobacterium tuberculosis to establish and maintain infection. Proc Natl Acad Sci U S A. 2010;107:9819–9824. [PubMed]
44. McKinney JD, Honer zu Bentrup K, Munoz-Elias EJ, Miczak A, Chen B, et al. Persistence of Mycobacterium tuberculosis in macrophages and mice requires the glyoxylate shunt enzyme isocitrate lyase. Nature. 2000;406:735–738. [PubMed]
45. Brzostek A, Pawelczyk J, Rumijowska-Galewicz A, Dziadek B, Dziadek J. Mycobacterium tuberculosis is able to accumulate and utilize cholesterol. J Bacteriol. 2009;191:6584–6591. [PMC free article] [PubMed]
46. Chang JC, Miner MD, Pandey AK, Gill WP, Harik NS, et al. igr Genes and Mycobacterium tuberculosis cholesterol metabolism. J Bacteriol. 2009;191:5232–5239. [PMC free article] [PubMed]
47. Pandey AK, Sassetti CM. Mycobacterial persistence requires the utilization of host cholesterol. Proc Natl Acad Sci U S A. 2008;105:4376–4380. [PubMed]
48. Van der Geize R, Yam K, Heuser T, Wilbrink MH, Hara H, et al. A gene cluster encoding cholesterol catabolism in a soil actinomycete provides insight into Mycobacterium tuberculosis survival in macrophages. Proc Natl Acad Sci U S A. 2007;104:1947–1952. [PubMed]
49. Mohn WW, van der Geize R, Stewart GR, Okamoto S, Liu J, et al. The actinobacterial mce4 locus encodes a steroid transporter. J Biol Chem. 2008;283:35368–35374. [PubMed]
50. Chang JC, Harik NS, Liao RP, Sherman DR. Identification of mycobacterial genes that alter growth and pathology in macrophages and in mice. J Infect Dis. 2007;196:788–795. [PubMed]
51. Kendall SL, Withers M, Soffair CN, Moreland NJ, Gurcha S, et al. A highly conserved transcriptional repressor controls a large regulon involved in lipid degradation in Mycobacterium smegmatis and Mycobacterium tuberculosis. Mol Microbiol. 2007;65:684–699. [PMC free article] [PubMed]
52. Yam KC, D'Angelo I, Kalscheuer R, Zhu H, Wang JX, et al. Studies of a ring-cleaving dioxygenase illuminate the role of cholesterol metabolism in the pathogenesis of Mycobacterium tuberculosis. PLoS Pathog. 2009;5:e1000344. [PMC free article] [PubMed]
53. Yang X, Gao J, Smith I, Dubnau E, Sampson NS. Cholesterol is not an essential source of nutrition for Mycobacterium tuberculosis during infection. J Bacteriol. 2011;193:1473–1476. [PMC free article] [PubMed]
54. Savvi S, Warner DF, Kana BD, McKinney JD, Mizrahi V, et al. Functional characterization of a vitamin B12-dependent methylmalonyl pathway in Mycobacterium tuberculosis: implications for propionate metabolism during growth on fatty acids. J Bacteriol. 2008;190:3886–3895. [PMC free article] [PubMed]
55. Upton AM, McKinney JD. Role of the methylcitrate cycle in propionate metabolism and detoxification in Mycobacterium smegmatis. Microbiology. 2007;153:3973–3982. [PubMed]
56. Munoz-Elias EJ, Upton AM, Cherian J, McKinney JD. Role of the methylcitrate cycle in Mycobacterium tuberculosis metabolism, intracellular growth, and virulence. Mol Microbiol. 2006;60:1109–1122. [PubMed]
57. Griffin JE, Pandey AK, Gilmore SA, Mizrahi V, McKinney JD, et al. Cholesterol catabolism by Mycobacterium tuberculosis requires transcriptional and metabolic adaptations. Chem Biol. 2012;19:218–227. [PMC free article] [PubMed]
58. Jain M, Petzold CJ, Schelle MW, Leavell MD, Mougous JD, et al. Lipidomics reveals control of Mycobacterium tuberculosis virulence lipids via metabolic coupling. Proc Natl Acad Sci U S A. 2007;104:5133–5138. [PubMed]
59. Singh A, Crossman DK, Mai D, Guidry L, Voskuil MI, et al. Mycobacterium tuberculosis WhiB3 maintains redox homeostasis by regulating virulence lipid anabolism to modulate macrophage response. PLoS Pathog. 2009;5:e1000545. [PMC free article] [PubMed]
60. Geiman DE, Raghunand TR, Agarwal N, Bishai WR. Differential gene expression in response to exposure to antimycobacterial agents and other stress conditions among seven Mycobacterium tuberculosis whiB-like genes. Antimicrob Agents Chemother. 2006;50:2836–2841. [PMC free article] [PubMed]
61. Fortune SM, Jaeger A, Sarracino DA, Chase MR, Sassetti CM, et al. Mutually dependent secretion of proteins required for mycobacterial virulence. Proc Natl Acad Sci U S A. 2005;102:10676–10681. [PubMed]
62. Raghavan S, Manzanillo P, Chan K, Dovey C, Cox JS. Secreted transcription factor controls Mycobacterium tuberculosis virulence. Nature. 2008;454:717–721. [PMC free article] [PubMed]
63. Jacques PE, Gervais AL, Cantin M, Lucier JF, Dallaire G, et al. MtbRegList, a database dedicated to the analysis of transcriptional regulation in Mycobacterium tuberculosis. Bioinformatics. 2005;21:2563–2565. [PubMed]
64. Madan Babu M, Teichmann SA, Aravind L. Evolutionary dynamics of prokaryotic transcriptional regulatory networks. J Mol Biol. 2006;358:614–633. [PubMed]
65. Guo M, Feng H, Zhang J, Wang W, Wang Y, et al. Dissecting transcription regulatory pathways through a new bacterial one-hybrid reporter system. Genome Res. 2009;19:1301–1308. [PubMed]
66. Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, et al. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000;11:4241–4257. [PMC free article] [PubMed]
67. Hecker M, Volker U. General stress response of Bacillus subtilis and other bacteria. Adv Microb Physiol. 2001;44:35–91. [PubMed]
68. Kaushal D, Schroeder BG, Tyagi S, Yoshimatsu T, Scott C, et al. Reduced immunopathology and mortality despite tissue persistence in a Mycobacterium tuberculosis mutant lacking alternative sigma factor, SigH. Proc Natl Acad Sci U S A. 2002;99:8330–8335. [PubMed]
69. Manganelli R, Voskuil MI, Schoolnik GK, Dubnau E, Gomez M, et al. Role of the extracytoplasmic-function sigma factor sigma(H) in Mycobacterium tuberculosis global gene expression. Mol Microbiol. 2002;45:365–374. [PubMed]
70. Honaker RW, Dhiman RK, Narayanasamy P, Crick DC, Voskuil MI. DosS responds to a reduced electron transport system to induce the Mycobacterium tuberculosis DosR regulon. J Bacteriol. 2010;192:6447–6455. [PMC free article] [PubMed]
71. Stewart GR, Wernisch L, Stabler R, Mangan JA, Hinds J, et al. Dissection of the heat-shock response in Mycobacterium tuberculosis using mutants and microarrays. Microbiology. 2002;148:3129–3138. [PubMed]
72. Stewart GR, Snewin VA, Walzl G, Hussell T, Tormay P, et al. Overexpression of heat-shock proteins reduces survival of Mycobacterium tuberculosis in the chronic phase of infection. Nat Med. 2001;7:732–737. [PubMed]
73. Gao CH, Yang M, He ZG. An ArsR-like transcriptional factor recognizes a conserved sequence motif and positively regulates the expression of phoP in mycobacteria. Biochem Biophys Res Commun. 2011;411:726–731. [PubMed]
74. Aldridge BB, Fernandez-Suarez M, Heller D, Ambravaneswaran V, Irimia D, et al. Asymmetry and aging of mycobacterial cells lead to variable growth and antibiotic susceptibility. Science. 2012;335:100–104. [PMC free article] [PubMed]
75. Blake WJ, Balazsi G, Kohanski MA, Isaacs FJ, Murphy KF, et al. Phenotypic consequences of promoter-mediated transcriptional noise. Mol Cell. 2006;24:853–865. [PubMed]
76. Sassetti CM, Boyd DH, Rubin EJ. Genes required for mycobacterial growth defined by high density mutagenesis. Mol Microbiol. 2003;48:77–84. [PubMed]
77. Sassetti CM, Rubin EJ. Genetic requirements for mycobacterial survival during infection. Proc Natl Acad Sci U S A. 2003;100:12989–12994. [PubMed]
78. Munoz-Elias EJ, McKinney JD. Mycobacterium tuberculosis isocitrate lyases 1 and 2 are jointly required for in vivo growth and virulence. Nat Med. 2005;11:638–644. [PMC free article] [PubMed]
79. Heym B, Stavropoulos E, Honore N, Domenech P, Saint-Joanis B, et al. Effects of overexpression of the alkyl hydroperoxide reductase AhpC on the virulence and isoniazid resistance of Mycobacterium tuberculosis. Infect Immun. 1997;65:1395–1401. [PMC free article] [PubMed]
80. Ng VH, Cox JS, Sousa AO, MacMicking JD, McKinney JD. Role of KatG catalase-peroxidase in mycobacterial pathogenesis: countering the phagocyte oxidative burst. Mol Microbiol. 2004;52:1291–1302. [PubMed]
81. Ferrer NL, Gomez AB, Neyrolles O, Gicquel B, Martin C. Interactions of attenuated Mycobacterium tuberculosis phoP mutant with human macrophages. PLoS One. 2010;5:e12978. [PMC free article] [PubMed]
82. Rengarajan J, Bloom BR, Rubin EJ. Genome-wide requirements for Mycobacterium tuberculosis adaptation and survival in macrophages. Proc Natl Acad Sci U S A. 2005;102:8327–8332. [PubMed]
83. Homolka S, Niemann S, Russell DG, Rohde K. Functional Genetic Diversity among M. tuberculosis Clinical Isolates: Delineation of Conserved Core and Lineage-Specific Transcriptomes during Intracellular Survival. PLoS Pathog. 2010;6:e1000988. [PMC free article] [PubMed]
84. Bishai WR, Dannenberg AM, Jr, Parrish N, Ruiz R, Chen P, et al. Virulence of Mycobacterium tuberculosis CDC1551 and H37Rv in rabbits evaluated by Lurie's pulmonary tubercle count method. Infect Immun. 1999;67:4931–4934. [PMC free article] [PubMed]
85. Manca C, Tsenova L, Barry CE, 3rd, Bergtold A, Freeman S, et al. Mycobacterium tuberculosis CDC1551 induces a more vigorous host response in vivo and in vitro, but is not more virulent than other clinical isolates. J Immunol. 1999;162:6740–6746. [PubMed]
86. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–210. [PMC free article] [PubMed]
87. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25:402–408. [PubMed]
88. Fleischmann RD, Alland D, Eisen JA, Carpenter L, White O, et al. Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains. J Bacteriol. 2002;184:5479–5490. [PMC free article] [PubMed]
89. Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998;393:537–544. [PubMed]
90. Abdallah AM, Gey van Pittius NC, Champion PA, Cox J, Luirink J, et al. Type VII secretion–mycobacteria show the way. Nat Rev Microbiol. 2007;5:883–891. [PubMed]
91. McLaughlin B, Chon JS, MacGurn JA, Carlsson F, Cheng TL, et al. A mycobacterium ESX-1-secreted virulence factor with unique requirements for export. PLoS Pathog. 2007;3:e105. [PMC free article] [PubMed]
Articles from PLoS Pathogens are provided here courtesy of
Public Library of Science