- Research article
- Open Access
- Open Peer Review
Profiling persistent tubercule bacilli from patient sputa during therapy predicts early drug efficacy
BMC Medicine volume 14, Article number: 68 (2016)
New treatment options are needed to maintain and improve therapy for tuberculosis, which caused the death of 1.5 million people in 2013 despite potential for an 86 % treatment success rate. A greater understanding of Mycobacterium tuberculosis (M.tb) bacilli that persist through drug therapy will aid drug development programs. Predictive biomarkers for treatment efficacy are also a research priority.
Methods and Results
Genome-wide transcriptional profiling was used to map the mRNA signatures of M.tb from the sputa of 15 patients before and 3, 7 and 14 days after the start of standard regimen drug treatment. The mRNA profiles of bacilli through the first 2 weeks of therapy reflected drug activity at 3 days with transcriptional signatures at days 7 and 14 consistent with reduced M.tb metabolic activity similar to the profile of pre-chemotherapy bacilli. These results suggest that a pre-existing drug-tolerant M.tb population dominates sputum before and after early drug treatment, and that the mRNA signature at day 3 marks the killing of a drug-sensitive sub-population of bacilli. Modelling patient indices of disease severity with bacterial gene expression patterns demonstrated that both microbiological and clinical parameters were reflected in the divergent M.tb responses and provided evidence that factors such as bacterial load and disease pathology influence the host-pathogen interplay and the phenotypic state of bacilli. Transcriptional signatures were also defined that predicted measures of early treatment success (rate of decline in bacterial load over 3 days, TB test positivity at 2 months, and bacterial load at 2 months).
This study defines the transcriptional signature of M.tb bacilli that have been expectorated in sputum after two weeks of drug therapy, characterizing the phenotypic state of bacilli that persist through treatment. We demonstrate that variability in clinical manifestations of disease are detectable in bacterial sputa signatures, and that the changing M.tb mRNA profiles 0–2 weeks into chemotherapy predict the efficacy of treatment 6 weeks later. These observations advocate assaying dynamic bacterial phenotypes through drug therapy as biomarkers for treatment success.
Tuberculosis (TB) caused the death of 1.5 million people in 2013 despite potential for an 86 % treatment success rate . New drug regimens are needed to maintain and improve therapy for tuberculosis, shortening treatment duration and targeting drug-resistant bacilli which complicate 3.5 % of new and 20.5 % of previously-treated TB cases . The standard chemotherapy regimen for drug-sensitive TB uses combinations of 4 drugs over 6 months. The recommended treatment for multidrug-resistant TB lasts 18–24 months or more, with increasingly toxic combinations of second-line drugs. Extended periods of chemotherapy are required to remove sub-populations of Mycobacterium tuberculosis (M.tb) bacilli that persist through the early phase of antimicrobial drug treatment . A recent treatment shortening trial failed to show non-inferiority despite evidence that the experimental regimens were more bactericidal in the first 4 months . Evidence for the existence of persister populations recalcitrant to treatment is accumulating; quantitation of mycobacteria in serial sputum samples during treatment with regimens containing isoniazid display a characteristic biphasic pattern of killing, indicative of the presence of multiple populations of M.tb . The exposure of sputum-derived bacilli to resuscitation-promoting factors unmasks a previously-unculturable drug-tolerant population of M.tb in sputum , and sub-populations of bacilli that only grow in liquid culture and not on solid media may represent 90 % of bacilli in sputum . Similar M.tb populations have also been identified in chronic murine tuberculosis models , and generated in vitro . Microfluidic systems have revealed heterogeneous mycobacterial responses to antibiotic exposure in genetically-homogenous populations, providing a basis for the generation of such sub-populations in vivo . It is unlikely that the duration of tuberculosis chemotherapy will be reduced until drug regimens are identified that can kill these persister sub-populations.
Success in tuberculosis chemotherapy is measured by the proportion of patients who fail therapy or who relapse after treatment is completed; therapy is monitored by counting M.tb in sputum and assaying markers of drug toxicity. Clinical or microbiological parameters, such as number of lung cavities or extent of lung cavitation, M.tb culture positivity at 2 months  or bacterial load in sputum before the start of treatment , are associated with early treatment success but are poor predictors of treatment outcome. For example, use of a combination of biomarkers, culture negativity at 2 months and low extent of cavitation by X-ray, failed to be predictive of individuals where treatment could be shortened from 6 to 4 months . However, molecular profiling assays, such as those used to identify tuberculosis disease from patient blood transcriptional signatures [13, 14], have not been applied to bacteria during patient drug therapy to assess the predictive power of dynamic bacterial responses to antimicrobial drug exposure.
The transcriptional signature of M.tb reflects the bacterial physiological state and offers insight into the mechanisms required to survive [15–17]. An increased understanding of which M.tb bacterial phenotypes survive chemotherapy during natural infection will aid the design of novel intervention strategies. To this end, M.tb bacilli have been profiled from in vitro models that mimic hypothesized features of M.tb persister populations [8, 18, 19]. Investigation of the mRNA signature of M.tb derived from human sputa revealed a slow/non-growing gene expression pattern alongside an accumulation of lipid bodies, a ‘fat and lazy’ phenotype . In addition, by mapping the differential expression of M.tb respiratory pathways, the microenvironment surrounding bacilli was predicted to be, at least in part, hypoxic . This observation was confirmed by 3-dimensional PET-CT imaging of human lungs that highlighted the hypoxic and dynamic nature of lesions within an individual . Transcriptional profiling the response of bacilli to drug exposure in vitro has helped to define antimicrobial drug class of action and mechanisms that may influence drug efficacy [22–24]. These drug-inducible signatures, used as a bioprobe, may also identify drug-tolerant M.tb populations by classifying divergent responses to drug exposure [25, 26].
This study profiles sputum M.tb transcriptional responses during the first 14 days of standard anti-tuberculous therapy, testing the hypothesis that persister-type bacilli are the dominant population in human sputum. Drug-induced changes in M.tb gene expression were observed 3 days after the start of chemotherapy that were not evident at 7 or 14 days. Furthermore, the profile of bacilli derived from sputum one or two weeks after drug therapy resembled pre-treatment sputum bacilli. This suggests that bacilli with a phenotype able to survive drug therapy were already present prior to the commencement of treatment. Bacteria with a drug-responsive phenotype during days 0 to 3 were no-longer present by day 7 and were presumably killed. Importantly, we demonstrate for the first time that the diverse pathology of human disease influences the phenotype of M.tb bacilli in sputum. We show that microbiological (bacterial load) and clinical (number of cavities) measures of disease were predicted from the changing M.tb sputum signature over time, as were parameters with prognostic value (positive TTP or MBL TB test, and bacterial load at 2 months).
Patient sample collection, clinical and molecular parameters
Subjects with active pulmonary tuberculosis (HIV negative, sputum smear-positive pulmonary TB) were recruited in primary healthcare tuberculosis clinics in the Western Cape Province, South Africa following local ethical approval (Stellenbosch University Health Research Ethics Committee, Study no 99/039). Patients consented to be involved in the study. Details of the patients in this cohort have been reported previously . Clinical parameters measuring severity of disease, such as average chest radiograph (CXR) score and number of observable cavities, were recorded. Briefly, full-sized CXRs (postero-anterior and lateral) were obtained and read by a pulmonologist blinded to patient clinical history using a standardized protocol . Expectorated sputum was immediately collected into 4 M GTC solution at the clinic as previously described  and frozen at -80 °C. For each patient, samples were collected before the start of chemotherapy, and 3, 7 and 14 days after initiation of treatment. Standard 5 day/week clinic-based directly observed treatment (DOT) was given by routine clinic nurses using fixed-dose combinations, with dose adjustment based on patient body weight. Treatment consisted of a 2-month intensive phase of rifampicin, isoniazid, pyrazinamide, and ethambutol, followed by a 4-month continuation phase of rifampicin and isoniazid. Treatment was monitored using microbiological (BACTEC 460) time-to-positivity scores (TTP) at diagnosis (day 0), day 7 and day 14, and by the molecular bacterial load (MBL) assay  at diagnosis (day 0), day 3, 7, 14 and 56. Participants’ molecular, microbiological and clinical parameters are detailed in Additional file 1: Table S1.
M.tb RNA extraction and amplification
Mycobacterial RNA was extracted from tuberculous sputa as previously described using the GTC/Trizol method . Briefly, sputum was thawed and bacterial pellets recovered from GTC by centrifugation at 1800 g for 30 minutes. Bacterial pellets were resuspended in Trizol (Life Technologies), disrupted using a ribolyzer (MP Biomedicals) and the nucleic acid recovered in the aqueous phase after addition of chloroform. The RNA preparations were purified and DNase-treated using RNeasy columns (Qiagen). Mycobacterial RNA yield and quality were assayed using the Nano-Drop ND-1000 Spectrophotometer (NanoDrop Technologies) and Agilent 2100 Bioanalyser (Agilent Technologies). RNA samples were amplified from 100 ng total RNA using the MessageAmp II Bacteria system (Life Technologies) [16, 28]. All sputum samples were extracted and amplified together to minimize technical variation.
Transcriptional profiling M.tb from sputa
Amplified M.tb RNA derived from 15 subjects at multiple time intervals before and during chemotherapy (totalling 52 samples) was profiled alongside M.tb H37Rv RNA extracted from in vitro log phase bacilli (two biological replicates hybridized in duplicate) as a standardized comparator [16, 20]. Amplified mycobacterial RNA (2 μg) was directly labelled with Cy3 fluorophore using the Universal Linkage System (ULS, Kreatech Diagnostics). Microarray hybridizations were conducted as previously described [16, 20] using an M.tb complex pan-genome microarray generated by the Bacterial Microarray Group at St. George’s (ArrayExpress accession number A-BUGS-41). Significantly differentially expressed genes were identified using moderated t-tests (p-value <0.05 with Benjamini and Hochberg multiple testing correction), and fold change >2 (for aerobic comparisons) and >1.5 (for sputum temporal responses) in GeneSpring 12.6 (Agilent Technologies). Samples were hierarchically clustered using Cluster and Treeview . Hypergeometric probability was used to identify significantly enriched transcriptional signatures from functional classifications or from published datasets. Short time-series expression miner (STEM) was used to determine significantly represented temporal gene expression profiles (p <0.05 after Bonferroni multiple testing correction), identifying time-dependent transcriptional patterns in bacilli extracted from human sputa . Significant genes identified in each comparison are detailed in Additional file 2: Table S2 and Additional file 3: Table S3.
Quantitative RT-PCR to verify transcriptional signatures
M.tb cDNA (20 μL) was prepared for each sample using 1 μg RNA per reaction (Maxima 1st strand cDNA synthesis kit for RT-qPCR, Thermo Scientific) and amplified according to the manufacturer’s instructions for high GC templates. For pooled patient analyses, 10 μL first strand synthesis products were combined for each timepoint. The two M.tb log phase samples were amplified in triplicate and combined. Quantitative 5-plex PCR was performed on a Rotor-Gene Q platform using QuantiTect Multiplex PCR NoROX Kit (Qiagen). A total of 5 μL pooled cDNA was added per 25 μL reaction with 2× QuantiTect reaction mix and 0.2 μM primer and probes (Additional file 4: Table S4). Multiplex qPCR was run for 40 cycles according to the manufacturer’s instructions. In each 5-plex set sigA was used to normalize input cDNA concentration. Reactions were run in duplicate and accepted if replicate Cq values were within 1 cycle. Cq values were used to determine changes in gene expression either relative to day 0 (before treatment) or to log phase bacilli using the 2-(delta delta Cq) method .
Modelling associations between RNA signatures and microbiological/clinical parameters
Gene expression data were gathered into a matrix X q ∈ ℝ52 × 4456 for quantile normalized and X d ∈ ℝ37 × 4456 for day 0 normalized data, with samples (patient, timepoint) in rows and genes in columns. Both matrices were standardized so each column had a mean of 0 and standard deviation of 1. Unsupervised dimensionality reduction technique principal component analysis (PCA)  was applied to the gene expression matrix X q to visualize temporal movement of the samples in principal component space. To investigate whether the gene expression trajectories over the first two weeks of therapy followed classifiable patterns, the patients (k) were divided into two classes based on the variance of the PC data points along the principal component axis: y k = − 1, if var(PC1) > var(PC2) and y k = 1, otherwise. That is, the class corresponds to the principal component axis along which the PC data points of a patient has larger variance. These target classes were modelled using support vector machine (SVM)  with gene expression data or clinical/microbiological variables as the predictors to define gene signatures or clinical variables, after applying stability selection or forward feature selection, respectively, that explained the patient trajectories.
The relationships between the M.tb transcriptional profile and patient clinical/microbiological variables were defined using machine learning methods. Time-to-positivity was modelled with X q and other clinical variables with X d . Firstly, stability selection , implemented according to the SCoRS framework, was used for feature selection . A total of 500 sub-sampling iterations were performed with a sub-sample that consisted of 500 features (columns in X) and 2/3 of the samples (rows in X). Stability selection yielded a set of stable, most predictive features of a particular clinical variable, which was used to train a regression or classification model. To optimize the number of features and to determine the predictive performance of the models, nested leave-one-out cross-validation was carried out, in which the inner loop determined the optimal number of variables and the outer loop the error of the model (Additional file 5: Figure S1). The reported errors are the average prediction errors of the test samples in the outer loop. Continuous clinical variables time-to-positivity and average chest X-ray score were modelled with L1-norm regression (Lasso) , and the performance evaluated with root mean squared error (RMSE) and Pearson correlation coefficient between the fitted and true values. Binary clinical variables, TB test positivity at week 8 (MBL or TTP positive test vs. negative tests at week 8), bacterial load at day 0 (MBL >= 106 bacilli defined as high, MBL<106 classed as low), and rate of decline in bacterial load from day 0 to 3 (ratio of MBL at day 3 relative to day 0, MBL ratio >= 10 defined as high, MBL ratio <10 classed as low) were modelled with SVM. Molecular bacterial load at week 8 was modelled with binary SVM by dividing the response values into two bins based on the distribution of the values (MBL week 8 <300 defined as low, MBL week 8 >=300 classed as high). Performance of the SVM models was measured with average classification error, the fraction of misclassified test samples in cross-validation. All analyses, except the SCoRS algorithm, were implemented with Matlab and R using their built-in functions.
The transcriptional states of M.tb bacilli derived from human sputa were mapped through the early stages of chemotherapy. M.tb bacilli were isolated from the sputa of 15 patients before chemotherapy was begun and then at 3, 7 and 14 days of therapy with the standard regimen (isoniazid, ethambutol, rifampicin and pyrazinamide). All patients had non recurrent active tuberculosis and were treated successfully for drug-sensitive tuberculosis and all patients were culture negative at 3 months.
The ‘fat and lazy’ transcriptome of pre-chemotherapy sputum-derived bacilli
The transcriptional profile of M.tb bacilli from sputa before the start of chemotherapy (day 0) was defined relative to bacilli grown in vitro in axenic log phase culture. A total of 1083 genes were significantly differentially expressed in sputum-derived bacilli (day 0) compared to log phase growth, which revealed adaptations to respiratory and metabolic pathways (Fig. 1a/b, Additional file 2: Table S2, Additional file 5: Figure S2 and Figure S3). The sputum transcriptional signature was dominated by the repression of genes involved in intermediate metabolism (II.A) and ribosome synthesis (II.A.1). These functional categories (defined by Cole et al. ) were significantly over-represented amongst down-regulated genes (hypergeometric probabilities of 9.4x10-15 and 3.1x10-30, respectively). This drop in markers of cellular activity was accompanied with the repression of central metabolism and lipid biosynthesis, including the citric acid cycle (gltA2, kgd, mdh, korA/B, sucC, rv0247c/48c, fumC and mqo), FAS-1 (fas), FAS-II (fabD, acpM, kasA/B, fabG1 and inhA), and mycolic acid synthesis and modification (mmaA2/3/4, cmaA2, pcaA, fadD32 and pks13) pathways [38, 39]. Conversely, genes involved in the glyoxylate shunt pathway and methylcitrate cycle (icl, prpC and rv1129c), catabolism of cholesterol and fatty acids , and the predicted triacylglycerol synthases (rv1425, rv1760, rv3087 and rv3371)  were induced in sputum-derived bacilli. Thus, the M.tb sputum transcriptome suggests that central carbon metabolism and general metabolic activity is reduced in bacilli expectorated from the lungs of tuberculosis patients, with an increased emphasis on lipolytic pathways. The transcriptomic evidence also predicted that the respiratory state of bacilli was altered in sputum with the down-regulation of NADH dehydrogenase I (nuoD/E/F/G/K), cytochrome C reductase (qcrA/B/C) and aa3 cytochrome C reductase (ctaC/D/E) pathways that are utilized in aerobic and microaerophilic conditions . Correspondingly, genes involved in nitrate reduction (narK2/3) were induced, indicating the potential for a switch to alternative electron acceptors and anaerobic respiration. These adaptations are likely to be mediated in part by the transcriptional regulator DosR (DevR) that is induced by hypoxia and nitric oxide , with several DosR-regulated genes significantly up-regulated in sputa (nrdZ, narK2, rv1738, pfkB, hspX, hrp1, rv3126c and rv3128c). Energy metabolism was also affected as ATP synthase genes (atpA-G) were repressed in sputum-derived bacilli; unsurprisingly the functional categories for aerobic respiration (1.B.6.a) and ATP-proton motive force (1.B.8)  were significantly enriched in those genes down-regulated in sputum (hg p-values 6.5x10-5 and 3.2x10-8, respectively). The differential expression of these pathways and gene families are summarized in Fig. 1a and Additional file 5: Figure S3. In addition, mRNA profiles of 10 metabolic and respiratory indicator genes (icl, hspX, sigG, tgs1, prpC, atpE, kasA, nuoA, qcrC and ctaD) were confirmed by quantitative RT-PCR (Fig. 1b), the high concordance of expression ratios between contrasting assays validating our approach.
The M.tb mRNA signature 7 and 14 days into chemotherapy resembles untreated bacilli
We were able to successfully extract, amplify and profile mycobacterial RNA from sputa, after 3, 7 and 14 days of standard chemotherapy despite falling bacterial viable counts. Hierarchical clustering revealed that the M.tb signatures from sputum-derived bacilli 14 days after the start of drug therapy were most similar to pre-chemotherapy (day 0) bacilli and day 7 bacilli, and that profiles from day 3 clustered away from the other time intervals (Fig. 2a, Additional file 5: Figure S4). These observations from unsupervised hierarchical clustering were supported by the results of pairwise significance testing. A total of 109 genes were identified as significantly differentially expressed at day 3 compared to day 0; in comparison to 37 or 42 genes that were significantly different between day 0 and days 7 or 14, respectively (Fig. 2b/c, Additional file 3: Table S3). No genes were significantly divergently expressed between day 7 and day 14 time intervals. The variation between sample profiles at each time point were similar (Spearman’s rank correlations of 0.871 to 0.901, Fig. 2c), suggesting that the number of statistically significant genes identified in each comparison could be used reliably as a measure of similarity. A pattern emerged, analogous to the unsupervised hierarchical clustering, demonstrating that M.tb responses at 7 and 14 days during chemotherapy were most similar to that of bacilli before drug therapy had begun.
The transcriptional signature from M.tb bacilli that have persisted in patients through 14 days of standard drug therapy was similar to the pre-chemotherapy profile; for example, 501 of the 528 genes induced in sputum at day 0 relative to log phase bacilli were also up-regulated in sputum at day 14 (Fig. 1a, Additional file 2: Table S2). Thus, the gross changes to metabolic and respiratory pathways are comparable between sputum-derived bacilli before and during chemotherapy. However, the magnitude of transcriptional adaptations was elevated at day 14 compared to day 0; using aerobic log phase bacilli as a comparator, 594 genes were significantly induced (605 genes repressed) at 14 days in contrast to 528 up-regulated genes (555 down-regulated genes) at day 0 (Additional file 5: Figure S2, Additional file 2: Table S2). A direct comparison of day 7 or day 14 transcriptional profiles to day 0 revealed that M.tb gene expression was predominately repressed with time (Fig. 2b/c); 37 genes were significantly down-regulated at day 7 compared to day 0, and 37 genes were repressed at day 14 (29 overlapping with day 7 signature), with five genes induced (Fig. 2b/c, Additional file 3: Table S3). Of the genes that were down-regulated over time, 88 % were also repressed in sputum relative to log phase bacilli, suggesting that the non/slow-growing phenotype of sputum-derived bacilli was further enhanced over time and through chemotherapy. These genes encode ribosomal proteins (rpmH and rpsJ), ribonucleotide reductase subunits (nrdH and nrdI) involved in the generation of precursors for DNA synthesis, and type VII secretion system elements (esxA, esxB, esxH, esxN, esxO, esxV, espD and espG3). In contrast, prpC, a methylcitrate synthase up-regulated in day 0 bacilli relative to log phase bacilli was down-regulated with drug therapy, indicating that specific modifications to metabolic pathways may occur over time.
Evidence of anti-mycobacterial drug action 3 days into treatment
An inflection point at day 3 of chemotherapy was identified by mapping significantly represented temporal gene expression profiles (Fig. 3a). Of the eight significantly represented gene curves describing the M.tb response to early drug therapy, six were modified at day 3 in comparison to day 0, 7 or 14 time intervals. The inference from this observation, that the mRNA signature 3 days after the start of chemotherapy was distinct from the other sputa time points, was reinforced by hierarchical clustering (Fig. 2a, Additional file 5: Figure S4) and significance testing (Fig. 2b/c). A total of 83 genes were significantly up-regulated 3 days after the start of drug treatment compared to pre-chemotherapy bacilli. The majority of these genes were also induced at day 3 relative to both 7 and 14 days of drug therapy (Additional file 5: Figure S5) suggesting that the day 3 transcriptional pattern represented a short-lived response to the start of drug treatment. This mRNA signature consisted of a diverse subset of genes involved in intermediary metabolism (argF, bkdC, galU, hisA, hisF, icl1, moaD1, pfkB, phoT and pstC1), cell wall metabolism (alr, fbpD and murC) and response to oxidative stress (alkB, cyp136 and rv0547c). The alternative sigma factor sigG, induced as part of the RecA-independent DNA damage response and implicated in the regulation of detoxification systems , was also up-regulated. Perhaps importantly, three genes involved in the export of antimicrobial drugs were induced at day 3: 1) rv1218c, encoding a putative ATP-dependent efflux pump regulated by RaaS (rv1219c) that is functionally significant in the response to antimicrobial drugs in non-permissive growth conditions ; 2) rv2688c, a predicted ABC fluoroquinolone efflux pump ; and 3) rv3066, a mycobacterial transcriptional regulator of the multidrug efflux pump Mmr (Rv3065) that has been shown to influence mycobacterial resistance to multiple toxic compounds . Quantitative RT-PCR confirmed the induction of bkdC, involved in branched amino acid metabolism, ndhA, encoding a nonproton-pumping type II NADH dehydrogenase, and nadA, a quinolinate synthetase involved in NAD biosynthesis, at day 3 relative to day 0 (Fig. 3b). The maximal induction of sigG, hspX and tgs1 at day 3 was also verified by qRT-PCR (Fig. 1b).
To further explore the relationship between sputum mRNA signatures and drug action in vivo, the transcriptional adaptations to a range of antimicrobial drugs in vitro  were mapped to the M.tb sputa profiles (Fig. 3c). Ten gene clusters (responsive to drug exposure) were identified to overlap with sputum transcriptional signatures exclusively after drug treatment had started. Transcriptional patterns reflecting exposure of bacilli to pyrazinamide and rifampicin were observed providing molecular evidence of in vivo drug action. Of the ten enriched drug responsive gene clusters, six were only significant at day 3. Quantitative RT-PCR of genes that are highly up-regulated after M.tb exposure to cell-wall inhibitors [22, 23], exemplified by iniB and efpA, were specifically induced at day 3 compared to drug-free day 0 bacilli (Fig. 3b). The expression of these benchmark genes for the activity of isoniazid and ethambutol decreased at 7 and 14 days relative to day 3. Thus, the transcriptional signature of bacilli 3 days into chemotherapy may reflect in vivo drug-induced changes; however, many of these responses to drug action were short-lived and were not evident 7 or 14 days after the start of drug therapy. We hypothesize that these signatures represent the killing of a drug-sensitive M.tb population after 3 days chemotherapy, and reveal the presence of a pre-existing drug-tolerant M.tb population that persists through early drug treatment.
M.tb transcriptional signatures reflect patient disease and predict treatment progress
All patients in this cohort were treated successfully for tuberculosis: nine of 15 patients were culture negative at 2 months, and all patients were culture negative at 3 months. Microbiological time-to-positivity scores (TTP) and molecular bacterial load (MBL) estimates were recorded alongside clinical parameters measuring severity of disease (chest radiograph score and the number of observable cavities) (Additional file 1: Table S1). Unsupervised principal component analysis (PCA) plotting the first and second principle components of the M.tb transcriptional signatures from each patient over the first 2 weeks of treatment time showed that M.tb responses followed distinct trajectories (Fig. 4a). This suggested that M.tb responses to drug therapy may differ between patients and that these bacterial mRNA signatures may reflect patient disease severity or predict treatment progress. Interestingly, two directions of travel emerged at 0–3 days with South-North and East-West trajectories most common (Fig. 4a and Additional file 5: Figure S6). Therefore, the PC trajectories were grouped into two classes according to the direction of movement in principal component space (South-North or East-West). Unsurprisingly, since the PCA patterns are a representation of the transcriptional responses, support vector machine classification (SVM) confirmed that PC trajectory class could be predicted from the gene signatures (classification error of 8 %, linear kernel). Furthermore, SVM was performed to test whether clinical or microbiological parameters collected for each patient would affect M.tb profiles such that the patient’s membership in the South-North or East-West class could be forecasted. The number of observed cavities in the lung (at diagnosis) and the bacterial load at 8 weeks (MBL measurement) predicted PC trajectory class with a classification error of 13 %, a success rate of 87 %. These observations demonstrated that the expression profile of M.tb in sputum was associated with specific measurable patient parameters, that variability in clinical and microbiological manifestations of disease could be detected in the bacterial sputa signatures, which may also be of prognostic value.
These findings were expanded using L1-norm regression to model the impact of continuous clinical and microbiological variables on M.tb sputa transcriptional profiles. Gene signatures were characterized that optimally designated time-to-positivity (day 0, day 7, day 14) with a RMSE (root mean squared error) of 3.6 hours. The fitted compared to true data are plotted in Fig. 4b/c, alongside the predictor gene set. The Pearson correlation coefficient between true and predicted values was 0.59 with p-value 4.4x10-5. Transcriptional patterns were also identified that correctly reflected patient chest X-ray score (at day 0) with a RMSE of 13.9 CXR score, Pearson correlation coefficient between true and predicted values of 0.73 with p-value 3.7x10-7 (Fig. 4d/e). M.tb expression profiles were also defined that were able to discriminate between high (>=106 bacilli measured by MBL) and low (<106 MBL bacilli) bacterial load at day 0 using a linear SVM (with classification error of 5 %). Thus, these clinical and microbiological parameters representing severity of disease and number of bacteria correlated with measurable changes in the M.tb transcriptional profile in sputa over the first 14 days of drug therapy. Notably, M.tb sputa signatures also predicted measures of early treatment response. The rate of decline in bacterial load from day 0 to day 3 (measured by MBL) was forecast with a Gaussian SVM model with a classification error of 3 %. TB test positivity at week 8 (positive wk 8 TTP or MBL assay) was correctly classified using a Gaussian SVM with a classification error of 11 % (successfully calling treatment progress 89 % of the time). Similarly, a high or low molecular bacterial load at week 8 could be determined from the sputa transcriptional dataset using SVM with a classification error of 0 %, 100 % test accuracy. Thus, the changing transcriptional profile of M.tb derived from human sputa reflects features relevant to patient disease and may be used to predict early treatment success.
M. tuberculosis bacilli derived from human sputa 2 weeks after the start of treatment should be portrayed as ‘fat, lazy and indifferent’ to drug therapy. In this study we were able to use genome-wide transcriptional profiling to map the mRNA signatures of M.tb from the sputa of clinically well-defined patients through the first 2 weeks of treatment, offering insight into the phenotypic state of bacilli that persist through chemotherapy. This contributes to our understanding of the in vivo efficacy of combination drug regimens and may aid novel anti-mycobacterial drug development programs targeting drug-persistent bacilli. The pre-chemotherapy M.tb signature defined here extends and validates the findings of a previous microarray study , confirming the M.tb sputa profile in a larger cohort of clinically-defined patients collected independently in a different country and assayed using a complementary profiling methodology. In addition, single gene inferences were verified by qRT-PCR using a panel of benchmark indicator genes.
The M.tb mRNA signatures 7 and 14 days after the start of drug therapy were most similar to that of bacilli before drug therapy had begun, as evidenced by hierarchical clustering and differential gene expression analysis. Transcriptional changes at 7 and 14 days also indicated that the M.tb day 0 sputum phenotype was enhanced with chemotherapy, with transcriptional markers of an active metabolic state further repressed over time. This mirrors findings from a recent study of M.tb in sputa using multiplex qRT-PCR that showed transition to a slow-growing low-metabolic activity phenotype after the start of drug treatment . The broad changes in M.tb sputum physiology reflected in these transcriptional patterns were conserved between studies, where day 14 signatures were more similar to day 7 than day 2  or day 3 (this study). However, we argue here that by using a well-defined log phase M.tb population as a comparator (as applied in many studies) enabled the pre-chemotherapy (day 0) M.tb population to be characterized , allowing us to propose that a pre-existing, slow/non-replicating and likely drug-tolerant M.tb population dominates sputum before, and after, early drug treatment. This conclusion is supported by evidence from other studies identifying large non-culturable on solid media or rpf-dependent M.tb sub-populations in sputum prior to treatment [5, 6].
An inflection point in the mRNA profile 3 days after the start of drug treatment, of stress-responsive genes including mediators of drug efflux, likely represents the effects of anti-TB drug action, using bacterial mRNA responses as a bioprobe for cidal drug action during natural disease. This transcriptional signature was only observed 3 days into therapy and was absent after 7 or 14 days, where the mRNA patterns were more similar to sputa bacilli before treatment. We hypothesize that the mRNA signature at day 3 marks the killing of a drug-sensitive sub-population of bacilli, since these responses are not detected after drug exposure in phenotypically drug-tolerant bacilli . By mapping previously-characterized in vitro M.tb transcriptional responses to antimicrobial drugs, we were able to define the in vivo action of standard regimen chemotherapy in a clinical setting. The expression of isoniazid-inducible benchmark genes (iniB, efpA) at day 3 alone suggest that the sputum bacterial population profiled here 1 week after the start of chemotherapy no longer respond to the antimicrobial actions of isoniazid. This finding supports the Walter et al. study that detected an isoniazid response in sputa bacilli 2 days after the start of drug therapy that disappeared 7 days into treatment , and mirrors evidence describing the bimodal early bactericidal activity of isoniazid in combination regimens [2, 4].
This study should be seen as capturing changes in M.tb mRNA abundance rather than differential transcriptional regulation since the structure of the underlying M.tb population is unknown. Therefore, the transcriptional read-out described here represents the average changes in gene expression from shifting mycobacterial sub-populations found in sputum. As such, lack of a response to drug therapy might reflect limited exposure of bacilli to antimicrobial drugs  and/or the presence of drug-tolerant bacterial populations that are able to survive drug treatment. The development of mycobacterial phenotypic drug tolerance has been described in many in vitro models for tuberculosis [8, 18, 19] and is often associated with a slow/non-replicating bacterial state, such as that inferred from the sputum M.tb mRNA signature. Moreover, recent observations using single cell reporter technology described the development of non-growing but metabolically-active mycobacterial sub-populations during murine infection . Our study further emphasizes the significance of drug-tolerant bacilli in human tuberculosis, identifying a transcriptionally-active M.tb population that persists through 2 weeks of standard chemotherapy.
Modelling M.tb gene expression in the sputa of patients over time revealed that the transcriptional pattern of bacilli varied from patient to patient with drug treatment. Correlation of patient indices of disease severity with bacterial mRNA signatures showed that basic microbiological (time-to-positivity, bacterial load at day 0) and clinical (chest radiograph score) parameters were reflected in the divergent M.tb responses. These associations suggest that factors such as bacterial load and disease pathology influence the host-pathogen interplay and, thus, the phenotypic state of bacilli, which in turn might affect the natural history of patient disease. Importantly, this demonstrates for the first time that the diverse pathology of tuberculosis affects measurable changes on the phenotype of M.tb bacilli in sputa. Notably, transcriptional signatures were also identified that predicted measures of early treatment success (rate of decline in bacterial load over 3 days, MBL or TTP positivity at 2 months, bacterial load at 2 months). Although this study was not significantly powered to test these signatures, these observations highlight the potential use of assaying dynamic bacterial phenotypes through drug therapy as biomarkers for treatment efficacy. These data are supported by recent evidence suggesting that higher percentages of lipid-body-positive acid-fast bacilli 3–4 weeks after the start of treatment, rather than initial baseline counts, correlated with treatment failure or relapse . We demonstrate here a novel concept and proof-of-principle that proposes using the changing transcriptional state of infecting bacilli to monitor treatment success.
This study defines the transcriptional signature of bacilli in sputum after 2 weeks of drug treatment, mapping the molecular phenotype of persister-type bacilli. We demonstrate for the first time that variability in clinical manifestations of disease are detectable in bacterial sputa signatures, and that the changing M.tb mRNA profiles 0–2 weeks into chemotherapy predict the efficacy of treatment 6 weeks later. These findings advocate a novel biomarker discovery strategy, profiling the phenotype of infecting bacteria, to find predictive markers of treatment success that are desperately needed in clinical trials and to stratify at-risk patients in the clinic.
Availability of data and materials
Fully annotated microarray data have been deposited in ArrayExpress, accession number: E-MTAB-3872.
- hg p-values:
- M.tb :
Molecular bacterial load
Principal component analysis
Root mean squared error
Support vector machine
World Health Organisation. Global tuberculosis report. Switzerland: WHO Press; 2014.
Mitchison D, Davies G. The chemotherapy of tuberculosis: past, present and future. Int J Tuberc Lung Dis. 2012;16(6):724–32.
Gillespie SH, Crook AM, McHugh TD, Mendel CM, Meredith SK, Murray SR, et al. Four-month moxifloxacin-based regimens for drug-sensitive tuberculosis. N Engl J Med. 2014;371(17):1577–87.
Jindani A, Dore CJ, Mitchison DA. Bactericidal and sterilizing activities of antituberculosis drugs during the first 14 days. Am J Respir Crit Care Med. 2003;167(10):1348–54.
Mukamolova GV, Turapov O, Malkin J, Woltmann G, Barer MR. Resuscitation-promoting factors reveal an occult population of tubercle bacilli in sputum. Am J Respir Crit Care Med. 2010;181(2):174–80.
Dhillon J, Fourie PB, Mitchison DA. Persister populations of Mycobacterium tuberculosis in sputum that grow in liquid but not on solid culture media. J Antimicrob Chemother. 2014;69(2):437–40.
Dhillon J, Lowrie DB, Mitchison DA. Mycobacterium tuberculosis from chronic murine infections that grows in liquid but not on solid medium. BMC Infect Dis. 2004;4:51.
Salina EG, Waddell SJ, Hoffmann N, Rosenkrands I, Butcher PD, Kaprelyants AS. Potassium availability triggers Mycobacterium tuberculosis transition to, and resuscitation from, non-culturable (dormant) states. Open Biol. 2014;4:140106. http://dx.doi.org/10.1098/rsob.140106
Wakamoto Y, Dhar N, Chait R, Schneider K, Signorino-Gelo F, Leibler S, et al. Dynamic persistence of antibiotic-stressed mycobacteria. Science. 2013;339(6115):91–5.
Benator D, Bhattacharya M, Bozeman L, Burman W, Cantazaro A, Chaisson R, et al. Rifapentine and isoniazid once a week versus rifampicin and isoniazid twice a week for treatment of drug-susceptible pulmonary tuberculosis in HIV-negative patients: a randomised clinical trial. Lancet. 2002;360(9332):528–34.
Hesseling AC, Walzl G, Enarson DA, Carroll NM, Duncan K, Lukey PT, et al. Baseline sputum time to detection predicts month two culture conversion and relapse in non-HIV-infected patients. Int J Tuberc Lung Dis. 2010;14(5):560–70.
Johnson JL, Hadad DJ, Dietze R, Maciel EL, Sewali B, Gitta P, et al. Shortening treatment in adults with noncavitary tuberculosis and 2-month culture conversion. Am J Respir Crit Care Med. 2009;180(6):558–63.
Berry MP, Graham CM, McNab FW, Xu Z, Bloch SA, Oni T, et al. An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature. 2010;466(7309):973–7.
Anderson ST, Kaforou M, Brent AJ, Wright VJ, Banwell CM, Chagaluka G, et al. Diagnosis of childhood tuberculosis and host RNA expression in Africa. N Engl J Med. 2014;370(18):1712–23.
Schnappinger D, Ehrt S, Voskuil MI, Liu Y, Mangan JA, Monahan IM, et al. Transcriptional adaptation of Mycobacterium tuberculosis within macrophages: insights into the phagosomal environment. J Exp Med. 2003;198(5):693–704.
Tailleux L, Waddell SJ, Pelizzola M, Mortellaro A, Withers M, Tanne A, et al. Probing host pathogen cross-talk by transcriptional profiling of both Mycobacterium tuberculosis and infected human dendritic cells and macrophages. PLoS One. 2008;3(1):e1403.
Rohde KH, Veiga DF, Caldwell S, Balazsi G, Russell DG. Linking the transcriptional profiles and the physiological states of Mycobacterium tuberculosis during an extended intracellular infection. PLoS Pathog. 2012;8(6):e1002769.
Betts JC, Lukey PT, Robb LC, McAdam RA, Duncan K. Evaluation of a nutrient starvation model of Mycobacterium tuberculosis persistence by gene and protein expression profiling. Mol Microbiol. 2002;43(3):717–31.
Deb C, Lee CM, Dubey VS, Daniel J, Abomoelak B, Sirakova TD, et al. A novel in vitro multiple-stress dormancy model for Mycobacterium tuberculosis generates a lipid-loaded, drug-tolerant, dormant pathogen. PLoS One. 2009;4(6):e6077.
Garton NJ, Waddell SJ, Sherratt AL, Lee SM, Smith RJ, Senner C, et al. Cytological and transcript analyses reveal fat and lazy persister-like bacilli in tuberculous sputum. PLoS Med. 2008;5(4):e75.
Barry III CE, Boshoff HI, Dartois V, Dick T, Ehrt S, Flynn J, et al. The spectrum of latent tuberculosis: rethinking the biology and intervention strategies. Nat Rev Microbiol. 2009;7(12):845–55.
Wilson M, DeRisi J, Kristensen HH, Imboden P, Rane S, Brown PO, et al. Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization. Proc Natl Acad Sci U S A. 1999;96(22):12833–8.
Boshoff HI, Myers TG, Copp BR, McNeil MR, Wilson MA, Barry III CE. The transcriptional responses of Mycobacterium tuberculosis to inhibitors of metabolism: novel insights into drug mechanisms of action. J Biol Chem. 2004;279(38):40174–84.
Makarov V, Manina G, Mikusova K, Mollmann U, Ryabova O, Saint-Joanis B, et al. Benzothiazinones kill Mycobacterium tuberculosis by blocking arabinan synthesis. Science. 2009;324(5928):801–4.
Karakousis PC, Williams EP, Bishai WR. Altered expression of isoniazid-regulated genes in drug-treated dormant Mycobacterium tuberculosis. J Antimicrob Chemother. 2008;61(2):323–31.
Tudo G, Laing K, Mitchison DA, Butcher PD, Waddell SJ. Examining the basis of isoniazid tolerance in nonreplicating Mycobacterium tuberculosis using transcriptional profiling. Future Med Chem. 2010;2(8):1371–83.
Honeyborne I, McHugh TD, Phillips PP, Bannoo S, Bateson A, Carroll N, et al. Molecular bacterial load assay, a culture-free biomarker for rapid and accurate quantification of sputum Mycobacterium tuberculosis bacillary load during treatment. J Clin Microbiol. 2011;49(11):3905–11.
Waddell SJ, Laing K, Senner C, Butcher PD. Microarray analysis of defined Mycobacterium tuberculosis populations using RNA amplification strategies. BMC Genomics. 2008;9(1):94.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95(25):14863–8.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006;7:191.
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(4):402–8.
Jolliffe IT. Principal Component Analysis. 2nd ed. New York: Springer-Verlag; 2002.
Cortes C, Vapnik V. Support-vector networks. Mach Learn. 1995;20(3):273–97.
Meinshausen N, Bühlmann P. Stability selection. J R Stat Soc Ser B (Stat Methodol). 2010;72:417–73.
Rondina J, Hahn T, de Oliveira L, Marquand A, Dresler T, Leitner T, Fallgatter A, Shawe-Taylor J, Mourao-Miranda J. SCoRS - a method based on stability for feature selection and mapping in neuroimaging. IEEE Trans Med Imaging. 2014;33(1):85–98.
Tibshirani R. Regression shrinkage and selection via the Lasso. J R Stat Soc Ser B (Stat Methodol). 1996;58(1):267–88.
Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998;393(6685):537–44.
Shi L, Sohaskey CD, Pfeiffer C, Datta P, Parks M, McFadden J, et al. Carbon flux rerouting during Mycobacterium tuberculosis growth arrest. Mol Microbiol. 2010;78(5):1199–215.
Rhee KY, de Carvalho LP, Bryk R, Ehrt S, Marrero J, Park SW, et al. Central carbon metabolism in Mycobacterium tuberculosis: an unexpected frontier. Trends Microbiol. 2011;19(7):307–14.
Van der Geize R, Yam K, Heuser T, Wilbrink MH, Hara H, Anderton MC, 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(6):1947–52.
Daniel J, Deb C, Dubey VS, Sirakova TD, Abomoelak B, Morbidoni HR, et al. Induction of a novel class of diacylglycerol acyltransferases and triacylglycerol accumulation in Mycobacterium tuberculosis as it goes into a dormancy-like state in culture. J Bacteriol. 2004;186(15):5017–30.
Shi L, Sohaskey CD, Kana BD, Dawes S, North RJ, Mizrahi V, et al. Changes in energy metabolism of Mycobacterium tuberculosis in mouse lung and under in vitro conditions affecting aerobic respiration. Proc Natl Acad Sci U S A. 2005;102(43):15629–34.
Kendall SL, Movahedzadeh F, Rison SC, Wernisch L, Parish T, Duncan K, et al. The Mycobacterium tuberculosis dosRS two-component system is induced by multiple stresses. Tuberculosis (Edinb). 2004;84(3-4):247–55.
Gaudion A, Dawson L, Davis E, Smollett K. Characterisation of the Mycobacterium tuberculosis alternative sigma factor SigG: its operon and regulon. Tuberculosis (Edinb). 2013;93(5):482–91.
Turapov O, Waddell SJ, Burke B, Glenn S, Sarybaeva AA, Tudo G, et al. Antimicrobial treatment improves mycobacterial survival in nonpermissive growth conditions. Antimicrob Agents Chemother. 2014;58(5):2798–806.
Pasca MR, Guglierame P, Arcesi F, Bellinzoni M, De Rossi E, Riccardi G. Rv2686c-Rv2687c-Rv2688c, an ABC fluoroquinolone efflux pump in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2004;48(8):3175–8.
Bolla JR, Do SV, Long F, Dai L, Su CC, Lei HT, et al. Structural and functional analysis of the transcriptional regulator Rv3066 of Mycobacterium tuberculosis. Nucleic Acids Res. 2012;40(18):9340–55.
Walter ND, Dolganov GM, Garcia BJ, Worodria W, Andama A, Musisi E, Ayakaka I, Van TT, Voskuil MI, de Jong BC. Transcriptional adaptation of drug-tolerant Mycobacterium tuberculosis during treatment of human tuberculosis. J Infect Dis. 2015;212(6):990–8.
Prideaux B, Via LE, Zimmerman MD, Eum S, Sarathy J, O'Brien P, Chen C, Kaya F, Weiner DM, Chen PY. The association between sterilizing activity and drug distribution into tuberculosis lesions. Nat Med. 2015;21(10):1223–7.
Manina G, Dhar N, McKinney JD. Stress and host immunity amplify Mycobacterium tuberculosis phenotypic heterogeneity and induce nongrowing metabolically active forms. Cell Host Microbe. 2015;17(1):32–46.
Sloan DJ, Mwandumba HC, Garton NJ, Khoo SH, Butterworth AE, Allain TJ, et al. Pharmacodynamic modeling of bacillary elimination rates and detection of bacterial lipid bodies in sputum to predict and understand outcomes in treatment of pulmonary tuberculosis. Clin Infect Dis. 2015;61(1):1–8.
SJW, PDB, SGH and TMcH are part of the PreDiCT-TB consortium (http://www.predict-tb.eu) which is funded from the Innovative Medicines Initiative Joint Undertaking (http://www.imi.europa.eu) under grant agreement No 115337, resources of which are composed of financial contribution from the European Union’s Seventh Framework Programme (FP7/2007-2013) and EFPIA companies’ in kind contribution. We would like to thank Kate Gould in the Bacterial Microarray Group at St. George’s for assistance with the microarray work. PDB acknowledges funding from the Wellcome Trust for the Bacterial Microarray Group at St. George’s (grant numbers 062511, 080039, and 086547) for access to M.tb microarrays. IH, TMcH and SGH acknowledge funding from the Medical Research Council (G0601466) and the European Metrology Research Programme (EMRP) INFECT-MET (HLT-08); the EMRP is jointly funded by the EMRP participating countries within EURAMET and the European Union. The samples originate from a sub-study of the GlaxoSmithKline Action TB Initiative, cofounded by the South African Technology for Human Resources and Industry Programme (THRIP). We thank the SA-MRC and research team at the Desmond Tutu TB Centre at Stellenbosch University, clinic staff and patients for their participation, and Cape Town City Health for permission to conduct this study.
The authors declare that they have no competing interests.
IH, TMcH, DFR, JR, PDB and SJW conceived and designed the experiments; KR, PDvH and GW enrolled patients; IH, IK, AC, DE and SJW performed the experiments; IH, TMcH, IK, AC, DFR, JR, PDB and SJW analysed the data; IH, TMcH, IK, AC, DE, KR, PDvH, SHG, DFR, GW, JR, PDB and SJW contributed to the writing of the manuscript. All authors read and approved the final manuscript.
Participants’ molecular, microbiological and clinical parameters. (XLSX 37 kb)
M.tuberculosis genes significantly differentially expressed in bacilli derived from sputum before the start of chemotherapy (day 0) and 3, 7 and 14 days during treatment compared to log phase aerobic growth. (XLSX 352 kb)
M.tuberculosis genes significantly differentially expressed in sputa over time. (XLSX 48 kb)
Quantitative RT-PCR multiplex primer and probe sequences. (XLSX 11 kb)
Figure S1. An illustration of the computation process to test associations between transcriptional signatures and patient variables. Figure S2. Venn diagrams describing the overlapping transcriptional signatures of bacilli in sputum relative to aerobic log phase growth. Figure S3. Box and whisker plots mapping the differential expression of gene families. Figure S4. Hierarchical clustering of the mean M.tb transcriptional profiles derived from sputa. Figure S5. Venn diagram highlighting genes significantly induced 3 days after the start of drug therapy. Figure S6. Contrasting patient trajectories as defined by principle component analysis plotting day 0 and day 3 timepoints only. (PDF 526 kb)