Background

According to the traditional tuberculosis (TB) treatment paradigm, the initial doses of treatment rapidly kill most Mycobacterium tuberculosis (Mtb) bacilli in sputum, yet many more months of daily treatment are required to eliminate a small, residual subpopulation of drug-tolerant bacilli. This paradigm has recently been challenged following the discovery that up to 90% of Mtb bacilli in sputum are culturable only with growth-factor supplementation. These "differentially culturable" bacilli are hypothesized to be more drug-tolerant than routinely culturable bacilli. This hypothesis implies an alternative paradigm in which TB treatment does not rapidly reduce the total Mtb population but only the small, routinely culturable subpopulation. To evaluate these competing paradigms, we developed a culture-independent method for quantifying the viable fraction of Mtb bacilli in sputum during treatment.

Methods

We used GeneXpert MTB/RIF to quantify Mtb DNA in sputa collected longitudinally from Ugandan adults taking standard 4-drug treatment for drug-susceptible pulmonary TB. We modeled GeneXpert cycle thresholds over time using nonlinear mixed-effects regression. We adjusted these models for clearance of DNA from killed-but-not-yet-degraded bacilli, assuming clearance half-lives ranging from 0 to 1.25 days. We used a convolution integral to quantify DNA from viable bacilli only, and converted cycle thresholds to Mtb genomic equivalents. We replicated our results in a South African cohort.

Results

We enrolled 41 TB patients in Uganda. Assuming a DNA-clearance half-life of 0 days, genomic equivalents of viable sputum bacilli decreased by 0.22 log/day until 8.8 days, then by 0.07 log/day afterwards. Assuming a DNA-clearance half-life of 1.25 days, genomic equivalents of viable bacilli decreased by 0.36 log/day until 5.0 days, then by 0.06 log/day afterwards. By day 7, viable Mtb had decreased by 97.2-98.8%. We found similar results for 19 TB patients in South Africa.

Discussion

Using a culture-independent method, we found that TB treatment rapidly eliminates most viable Mtb in sputum. These findings are incompatible with the hypothesis that differentially culturable bacilli are drug-tolerant.

Conclusions

A culture-independent method for measuring viable Mtb in sputum during treatment corroborates the traditional TB treatment paradigm in which a rapid bactericidal phase precedes slow, elimination of a small, residual bacillary subpopulation.

Droplet digital PCR is a particularly valuable tool for ratiometric assays because it provides simultaneous absolute quantification of two target sequences in a single assay. This manuscript addresses a challenge in establishing a new ratiometric droplet digital PCR assay for use in sputum, the rRNA synthesis ratio. In principle, the methods established to evaluate precision and determine the limit of quantification for a single measurand cannot be applied to a ratiometric assay. The precision of a ratio depends on precision in both the numerator and denominator. Here, we evaluated the MOVER approximated coefficient of variation as indicator of assay precision that does not require technical replicates. We estimated the MOVER approximated coefficient of variation in dilution series and routine assays and evaluated its agreement with the traditional coefficient of variation. We found that the MOVER approximated coefficient of variation was able to recapitulate the traditional coefficient of variation without the requirement for replicate assays. We also demonstrated that the MOVER approximated coefficient of variation threshold can be used to define the limit of quantification of the rRNA synthesis Ratio. In conclusion, the MOVER approximated coefficient of variation may be useful not only for the rRNA synthesis ratio but for other assays that measure ratios via droplet digital PCR.

Tuberculosis lung lesions are complex and harbor heterogeneous microenvironments that influence antibiotic effectiveness. Major strides have been made recently in understanding drug pharmacokinetics in pulmonary lesions, but the bacterial phenotypes that arise under these conditions and their contribution to drug tolerance are poorly understood. A pharmacodynamic marker called the RS ratio® quantifies ongoing rRNA synthesis based on the abundance of newly synthesized precursor rRNA relative to mature structural rRNA. Application of the RS ratio in the C3HeB/FeJ mouse model demonstrated that Mycobacterium tuberculosis populations residing in different tissue microenvironments are phenotypically distinct and respond differently to drug treatment with rifampin, isoniazid, or bedaquiline. This work provides a foundational basis required to address how anatomic and pathologic microenvironmental niches may contribute to long treatment duration and drug tolerance during the treatment of human tuberculosis.

Murine tuberculosis drug efficacy studies have historically monitored bacterial burden based on CFU of Mycobacterium tuberculosis in lung homogenate. In an alternative approach, a recently described molecular pharmacodynamic marker called the RS ratio quantifies drug effect on a fundamental cellular process, ongoing rRNA synthesis. Here, we evaluated the ability of different pharmacodynamic markers to distinguish between treatments in three BALB/c mouse experiments at two institutions. We confirmed that different pharmacodynamic markers measure distinct biological responses. We found that a combination of pharmacodynamic markers distinguishes between treatments better than any single marker. The combination of the RS ratio with CFU showed the greatest ability to recapitulate the rank order of regimen treatment-shortening activity, providing proof of concept that simultaneous assessment of pharmacodynamic markers measuring different properties will enhance insight gained from animal models and accelerate development of new combination regimens. These results suggest potential for a new era in which antimicrobial therapies are evaluated not only on culture-based measures of bacterial burden but also on molecular assays that indicate how drugs impact the physiological state of the pathogen.

To address the ongoing global tuberculosis crisis, there is a need for shorter, more effective treatments. A major reason why tuberculosis requires prolonged treatment is that, following a short initial phase of rapid killing, the residual Mycobacterium tuberculosis withstands drug killing. Because existing methods lack sensitivity to quantify low-abundance mycobacterial RNA in drug-treated animals, cellular adaptations of drug-exposed bacterial phenotypes in vivo remain poorly understood. Here, we used a novel RNA-seq method called SEARCH-TB to elucidate the Mycobacterium tuberculosis transcriptome in mice treated for up to 28 days with standard doses of isoniazid, rifampin, pyrazinamide, and ethambutol. We compared murine results with in vitro SEARCH-TB results during exposure to the same regimen. Treatment suppressed genes associated with growth, transcription, translation, synthesis of rRNA proteins, and immunogenic secretory peptides. Bacteria that survived prolonged treatment appeared to transition from ATP-maximizing respiration toward lower-efficiency pathways and showed modification and recycling of cell wall components, large-scale regulatory reprogramming, and reconfiguration of efflux pump expression. Although the pre-treatment in vivo and in vitro transcriptomes differed profoundly, genes differentially expressed following treatment in vivo and in vitro were similar, with differences likely attributable to immunity and drug pharmacokinetics in mice. These results reveal cellular adaptations of Mycobacterium tuberculosis that withstand prolonged drug exposure in vivo, demonstrating proof of concept that SEARCH-TB is a highly granular pharmacodynamic readout. The surprising finding that differential expression is concordant in vivo and in vitro suggests that insights from transcriptional analyses in vitro may translate to the mouse. IMPORTANCE A major reason that curing tuberculosis requires prolonged treatment is that drug exposure changes bacterial phenotypes. The physiologic adaptations of Mycobacterium tuberculosis that survive drug exposure in vivo have been obscure due to low sensitivity of existing methods in drug-treated animals. Using the novel SEARCH-TB RNA-seq platform, we elucidated Mycobacterium tuberculosis phenotypes in mice treated for with the global standard 4-drug regimen and compared them with the effect of the same regimen in vitro. This first view of the transcriptome of the minority Mycobacterium tuberculosis population that withstands treatment in vivo reveals adaptation of a broad range of cellular processes, including a shift in metabolism and cell wall modification. Surprisingly, the change in gene expression induced by treatment in vivo and in vitro was largely similar. This apparent portability from in vitro to the mouse provides important new context for in vitro transcriptional analyses that may support early preclinical drug evaluation.

There is a critical need for improved pharmacodynamic markers for use in human tuberculosis (TB) drug trials. Pharmacodynamic monitoring in TB has conventionally used culture or molecular methods to enumerate the burden of Mycobacterium tuberculosis organisms in sputum. A recently proposed assay called the rRNA synthesis (RS) ratio measures a fundamentally novel property, how drugs impact ongoing bacterial rRNA synthesis. Here, we evaluated RS ratio as a potential pharmacodynamic monitoring tool by testing pretreatment sputa from 38 Ugandan adults with drug-susceptible pulmonary TB. We quantified the RS ratio in paired pretreatment sputa and evaluated the relationship between the RS ratio and microbiologic and molecular markers of M. tuberculosis burden. We found that the RS ratio was highly repeatable and reproducible in sputum samples. The RS ratio was independent of M. tuberculosis burden, confirming that it measures a distinct new property. In contrast, markers of M. tuberculosis burden were strongly associated with each other. These results indicate that the RS ratio is repeatable and reproducible and provides a distinct type of information from markers of M. tuberculosis burden. IMPORTANCE This study takes a major next step toward practical application of a novel pharmacodynamic marker that we believe will have transformative implications for tuberculosis. This article follows our recent report in Nature Communications that an assay called the rRNA synthesis (RS) ratio indicates the treatment-shortening of drugs and regimens. Distinct from traditional measures of bacterial burden, the RS ratio measures a fundamentally novel property, how drugs impact ongoing bacterial rRNA synthesis.

There is urgent need for new drug regimens that more rapidly cure tuberculosis (TB). Existing TB drugs and regimens vary in treatment-shortening activity, but the molecular basis of these differences is unclear, and no existing assay directly quantifies the ability of a drug or regimen to shorten treatment. Here, we show that drugs historically classified as sterilizing and non-sterilizing have distinct impacts on a fundamental aspect of Mycobacterium tuberculosis physiology: ribosomal RNA (rRNA) synthesis. In culture, in mice, and in human studies, measurement of precursor rRNA reveals that sterilizing drugs and highly effective drug regimens profoundly suppress M. tuberculosis rRNA synthesis, whereas non-sterilizing drugs and weaker regimens do not. The rRNA synthesis ratio provides a readout of drug effect that is orthogonal to traditional measures of bacterial burden. We propose that this metric of drug activity may accelerate the development of shorter TB regimens.