Abstract
Mixed infections comprising multiple Mycobacterium tuberculosis Complex (MTBC) strains are observed in populations with high incidence rates of tuberculosis (TB), yet the difficulty to detect these via conventional diagnostic approaches has resulted in their contribution to TB epidemiology and treatment outcomes being vastly underrecognised. In endemic regions, detection of all component strains is crucial for accurate reconstruction of TB transmission dynamics. Currently available tools for detecting mixed infections from whole genome sequencing (WGS) data have insufficient sensitivity to detect low-frequency mixtures with less than 10% minor strain fraction, leading to a systematic underestimation of the frequency of mixed infection. Our R package, TBtypeR, identifies mixed infections from whole genome sequencing by comparing sample data to an expansive phylogenetic SNP panel of over 10,000 sites and 164 MTBC strains. A statistical likelihood is derived for putative strain mixtures based on the observed reference and alternative allele counts at each site under the binomial distribution. This provides robust and high-resolution sublineage classification for both single- and mixed-infections with as low as 1% minor strain frequency. Benchmarking with simulated in silico and in vitro mixture data demonstrates the superior performance of TBtypeR over existing tools, particularly in detecting low frequency mixtures. We apply TBtypeR to 5,000 MTBC WGS from a published dataset and find a 6-fold higher rate of mixed infection than existing methods. The TBtypeR R package and accompanying end-to-end Nextflow pipeline are available at github.com/bahlolab/TBtypeR.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was supported by a Walter and Eliza Hall Institute Innovation Fund grant to A.K.C and M.B. M.B. was supported by an NHMRC Investigator grant (GNT1195236). A.K.C was partially supported by an NHMRC Ideas grant (GNT2020750). Additional funding was provided by the NHMRC Independent Research Institute Infrastructure Support Scheme and the Victorian State Government Operational Infrastructure Program.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
All data and source code produced are available online in the TBtypeR GitHub repository
