Abstract:Tuberculosis (TB), the disease caused by Mycobacterium tuberculosis (Mtb), remains a major health problem with 10.6 million cases of the disease and 1.6 million deaths in 2021. It is well understood that pulmonary TB is due to Mtb growth in the lung but quantitative estimates of rates of Mtb replication and death in lungs of patients or animals such as monkeys or rabbits remain largely unknown. We performed experiments with rabbits infected with a novel, virulent clinical Mtb isolate of the Beijing lineage, HN878, carrying an unstable plasmid pBP10. In our in vitro experiments we found that pBP10 is more stable in HN878 strain than in a more commonly used laboratory-adapted Mtb strain H37Rv (the segregation coefficient being $s=0.10$ in HN878 vs. $s=0.18$ in H37Rv). Interestingly, the kinetics of plasmid-bearing bacteria in lungs of Mtb-infected rabbits did not follow an expected monotonic decline; the percent of plasmid-bearing cells increased between 28 and 56 days post-infection and remained stable between 84 and 112 days post-infection despite a large increase in bacterial numbers in the lung at late time points. Mathematical modeling suggested that such a non-monotonic change in the percent of plasmid-bearing cells can be explained if the lung Mtb population consists of several (at least 2) sub-populations with different replication/death kinetics: one major population expanding early and being controlled/eliminated, while another, a smaller population expanding at later times causing a counterintuitive increase in the percent of plasmid-bearing cells. Importantly, a model with one kinetically homogeneous Mtb population could not explain the data including when the model was run stochastically. Given that in rabbits HN878 strain forms well circumscribed granulomas, our results suggest independent bacterial dynamics in subsets of such granulomas. Our model predictions can be tested in future experiments in which HN878-pBP10 dynamics in individual granulomas is followed over time. Taken together, our new data and mathematical modeling-based analyses illustrate differences in Mtb dynamics in mice and rabbits confirming a perhaps somewhat obvious observation that ``rabbits are not mice''.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the replication kinetics of Mycobacterium tuberculosis (Mtb) when it infects the lungs of rabbits. Specifically, the researchers focused on whether the replication rate and death rate of Mtb in different regions or lesions of the rabbit lungs are consistent, and how these kinetics affect the infection process. ### Background and Motivation Tuberculosis (TB) is a disease caused by Mycobacterium tuberculosis (Mtb), which causes millions of people to get sick and die every year. Although there is some understanding of the replication and mortality of Mtb in mouse models, it is not clear whether these results are applicable to other mammalian models (such as rabbits). In addition, the kinetics of Mtb in different hosts may be significantly different, which may affect the progress of the disease and treatment strategies. ### Research Methods The researchers used a clinical isolate HN878 carrying an unstable plasmid pBP10 to infect rabbits, and studied the kinetics of Mtb by tracking the total number of bacteria and the proportion of plasmid - carrying cells at different time points after infection. The plasmid pBP10 will be lost with a certain probability during cell division, so the replication and mortality of bacteria can be inferred by detecting the proportion of plasmid - carrying cells. ### Key Findings 1. **Plasmid Stability**: The researchers found that in vitro experiments, the stability of pBP10 in the HN878 strain was higher than that in the commonly used H37Rv strain in the laboratory (the segregation coefficient $ s $ was 0.10 vs. 0.18). 2. **Non - monotonic Change**: In the lungs of infected rabbits, the proportion of plasmid - carrying cells increased between the 28th and 56th days after infection, and then remained stable between the 84th and 112th days, although the total number of bacteria increased significantly during the same period. This phenomenon is inconsistent with the expected monotonic decreasing trend. 3. **Multi - subpopulation Model**: Mathematical modeling shows that this non - monotonic change can be explained by assuming that the Mtb population in the lungs consists of at least two subpopulations with different replication / death kinetics. One of the main subpopulations expands in the early stage but is then controlled or eliminated, while another smaller subpopulation expands in the later stage, resulting in an abnormal increase in the proportion of plasmid - carrying cells. ### Conclusions The research results show that the replication kinetics of Mtb in the rabbit lungs may be heterogeneous, that is, different lesion areas may have different replication and death rates. This finding emphasizes the importance of different mammalian models in studying Mtb kinetics, and provides a new direction for future research, especially by tracking the Mtb dynamics in a single lesion area to evaluate its heterogeneity. ### Formulas - The calculation formula of the segregation coefficient $ s $: \[ s=\frac{\text{Slope}_F}{\text{Slope}_N} \] where $\text{Slope}_F$ is the exponential decay rate of the proportion of plasmid - carrying cells, and $\text{Slope}_N$ is the exponential growth rate of the total number of bacteria. - The calculation formulas of Mtb replication rate $ r_i $ and mortality rate $ \delta_i $: \[ r_i = \frac{\text{Slope}_F}{s} \] \[ \delta_i=\frac{\text{Slope}_F}{s}-\text{Slope}_N \] These formulas are used to estimate the replication and mortality rates of Mtb from experimental data.

Mathematical modeling suggests heterogeneous replication of Mycobacterium tuberculosis in rabbits

The rabbit model for spinal tuberculosis: An overview

Establishment of mouse models of persistent tuberculosis and characteristics of that infection

Animal Models of Mycobacterium Tuberculosis Infection: an Update

Heterogeneity in immune cell composition is associated with Mycobacterium tuberculosis replication at the granuloma level

Local Inflammation, Dissemination and Coalescence of Lesions Are Key for the Progression toward Active Tuberculosis: The Bubble Model

Immunopathology of Pulmonary Mycobacterium tuberculosis Infection in a Humanized Mouse Model

Mathematical models of drug-resistant tuberculosis lack bacterial heterogeneity: A systematic review

Computational pharmacokinetics/pharmacodynamics of rifampin in a mouse tuberculosis infection model

Mycobacterium Tuberculosis from Chronic Murine Infections That Grows in Liquid but Not on Solid Medium

Establishment of a Rabbit Model of Spinal Tuberculosis Using Mycobacterium Tuberculosis Strain H37Rv.

Establishment of a new zealand rabbit model of spinal tuberculosis.

The Multistate Tuberculosis Pharmacometric Model: a Semi-Mechanistic Pharmacokinetic-Pharmacodynamic Model for Studying Drug Effects in an Acute Tuberculosis Mouse Model

The use of Mycobacterium tuberculosis H37Ra-infected immunocompetent mice as an in vivo model of persisters

An Examination of a Rabbit Model of Liquefaction Induced by Mycobacterium

One Size Fits All? Not in in Vivo Modeling of Tuberculosis Chemotherapeutics

Assessing the Impact of Vaccination Strategies and Re-infection Dynamics on Mycobacterium tuberculosis Transmission by Mathematical Modelling

Systems genetics uncover new loci containing functional gene candidates in Mycobacterium tuberculosis -infected Diversity Outbred mice

Experiment on Mycoplasma bovis Infection in Rabbits

Agent-based model predicts that layered structure and 3D movement work synergistically to reduce bacterial load in 3D in vitro models of tuberculosis granuloma

Efficient modelling of infectious diseases in wildlife: a case study of bovine tuberculosis in wild badgers