Immortal genome assumption significantly underestimates replication and death rates of Mycobacterium tuberculosis in mice and monkeys

Immune correlates of protection against infection with Mycobacterium tuberculosis (Mtb) or against tuberculosis disease (TB) remain poorly defined. The ratio of colony forming units (CFUs) to chromosomal equivalents (CEQs), Z = CFUs/CEQs, recovered either from the whole lung (mice) or individual lesions (monkeys or rabbits) of Mtb-infected animals has been used as a metric for how effectively Mtb is killed in vivo. However, the contribution of bacterial killing to changes in the CFU/CEQ ratio during an infection has not been rigorously investigated. We developed alternative mathematical models to study the dynamics of CFUs, CEQs, and their ratio during an Mtb infection. We find that the ratio Z alone cannot be used to infer the death rate of bacteria, unless the dynamics of CEQs and CFUs are entirely uncoupled, which is biologically unreasonable and inconsistent with the view that CEQs reflect an accumulated burden of both viable and non-viable bacteria. Importantly, we estimate a decay rate of 3.6%/day (a half life of about 20 days) of Mtb H37Rv CEQs in B6 mice that is similar to 4%/day, previously found for Mtb Erdman in cynomolgus macaques. While the estimated Mtb DNA decay rate seems small, we found that estimated rates of Mtb replication and death/killing are still extremely sensitive even to slow decay of Mtb DNA, in part, because Mtb replication and death rates are also small especially during chronic phases of infection. By applying our models to data on Mtb dynamics during the first 3 weeks of infection in macaques, we provide evidence of substantial killing of Mtb bacteria, prior to arrival of adaptive immunity to the site of infection, challenging the previously established notion of non-dying bacteria in the absence of T cell immunity and granuloma formation. We also propose experiments that will allow more accurately to measure the rate of Mtb DNA loss, helping more rigorously quantify impact of immunity on within-host Mtb dynamics.