Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches.

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

Point-of-care C-reactive protein and Xpert MTB/RIF Ultra for tuberculosis screening and diagnosis in unselected antiretroviral therapy initiators: a prospective, cross-sectional, diagnostic accuracy study.

BACKGROUND: Tuberculosis, a major cause of death in people living with HIV, remains challenging to diagnose. Diagnostic accuracy data are scarce for promising triage and confirmatory tests such as C-reactive protein (CRP), sputum and urine Xpert MTB/RIF Ultra (Xpert Ultra), and urine Determine TB LAM Ag (a lateral flow lipoarabinomannan [LF-LAM] test), without symptom selection. We evaluated novel triage and confirmatory tests in ambulatory people with HIV initiating antiretroviral therapy (ART). METHODS: 897 ART-initiators were recruited irrespective of symptoms and sputum induction offered. For triage (n=800), we evaluated point-of-care blood-based CRP testing, compared with the WHO-recommended four-symptom screen (W4SS). For sputum-based confirmatory testing (n=787), we evaluated Xpert Ultra versus Xpert MTB/RIF (Xpert). For urine-based confirmatory testing (n=732), we evaluated Xpert Ultra and LF-LAM. We used a sputum culture reference standard. FINDINGS: 463 (52%) of 897 participants were female. The areas under the receiver operator characteristic curves for CRP was 0·78 (95% CI 0·73-0·83) and for number of W4SS symptoms was 0·70 (0·64-0·75). CRP (≥10 mg/L) had similar sensitivity to W4SS (77% [95% CI 68-85; 80/104] vs 77% [68-85; 80/104]; p>0·99] but higher specificity (64% [61-68; 445/696] vs 48% [45-52; 334/696]; p<0·0001]; reducing unnecessary confirmatory testing by 138 (95% CI 117-160) per 1000 people and number-needed-to-test from 6·91 (95% CI 6·25-7·81) to 4·87 (4·41-5·51). Sputum samples with Xpert Ultra, which required induction in 49 (31%) of 158 of people (95% CI 24-39), had higher sensitivity than Xpert (71% [95% CI 61-80; 74/104] vs 56% [46-66; 58/104]; p<0·0001). Of the people with one or more confirmatory sputum or urine test results that were positive, the proportion detected by Xpert Ultra increased from 45% (26-64) to 66% (46-82) with induction. Programmatically done haemoglobin, triage test combinations, and urine tests showed comparatively worse results. INTERPRETATION: CRP is a more specific triage test than W4SS in those initiating ART. Sputum induction improves diagnostic yield. Sputum samples with Xpert Ultra is a more accurate confirmatory test than with Xpert. FUNDING: South African Medical Research Council, EDCTP2, US National Institutes of Health-National Institute of Allergy and Infectious Diseases.

Latent tuberculosis infection is associated with an enrichment of short chain fatty acid producing bacteria in the stool of women living with HIV.

Background: Latent tuberculosis infection (LTBI) is common in people living with HIV (PLHIV) in high TB burden settings. Active TB is associated with specific stool taxa; however, little is known about the stool microbiota and LTBI, including in PLHIV. Method: Within a parent study that recruited adult females with HIV from Cape Town, South Africa into predefined age categories (18-25, 35-60 years), we characterised the stool microbiota of those with [interferon-γ release assay (IGRA)- and tuberculin skin test (TST)-positive] or without (IGRA- and TST- negative) LTBI (n=25 per group). 16S rRNA DNA sequences were analysed using QIIME2, Dirichlet Multinomial Mixtures, DESeq2 and PICRUSt2. Results: No α- or ß-diversity differences occurred by LTBI status; however, LTBI-positives were Faecalibacterium-, Blautia-, Gemmiger-, Bacteroides-enriched and Moryella-, Atopobium-, Corynebacterium-, Streptococcus-depleted. Inferred metagenome data showed LTBI-negative-enriched pathways included several involved in methylglyoxal degradation, L-arginine, putrescine, 4-aminobutanoate degradation and L-arginine and ornithine degradation. Stool from LTBI-positives demonstrated differential taxa abundance based on a quantitative response to antigen stimulation (Acidaminococcus-enrichment and Megamonas-, Alistipes-, and Paraprevotella-depletion associated with higher IGRA or TST responses, respectively). In LTBI-positives, older people had different ß-diversities than younger people whereas, in LTBI-negatives, no differences occurred across age groups. Conclusion: Amongst female PLHIV, those with LTBI had, vs. those without LTBI, Faecalibacterium, Blautia, Gemmiger, Bacteriodes-enriched, which are producers of short chain fatty acids. Taxonomic differences amongst people with LTBI occurred according to quantitative response to antigen stimulation and age. These data enhance our understanding of the microbiome's potential role in LTBI.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA.

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

Identification and molecular characterization of Mycobacterium bovis DNA in GeneXpert® MTB/RIF ultra-positive, culture-negative sputum from a rural community in South Africa.

This study investigated the presence of Mycobacterium bovis (M. bovis) DNA in archived human sputum samples previously collected from residents who reside adjacent to the M. bovis-endemic Hluhluwe-iMfolozi wildlife park, South Africa (SA). Sixty-eight sputum samples were GeneXpert MTB/RIF Ultra-positive for M. tuberculosis complex (MTBC) DNA but culture negative for M. tuberculosis. Amplification and Sanger sequencing of hsp65 and rpoB genes from DNA extracted from stored heat-inactivated sputum samples confirmed the presence of detectable amounts of MTBC from 20 out of the 68 sputum samples. Region of difference PCR, spoligotyping and gyrB long-read amplicon deep sequencing identified M. bovis (n = 10) and M. tuberculosis (n = 7). Notably, M. bovis spoligotypes SB0130 and SB1474 were identified in 4 samples, with SB0130 previously identified in local cattle and wildlife and SB1474 exclusively in African buffaloes in the adjacent park. M. bovis DNA in sputum, from people living near the park, underscores zoonotic transmission potential in SA. Identification of spoligotypes specifically associated with wildlife only and spoligotypes found in livestock as well as wildlife, highlights the complexity of TB epidemiology at wildlife-livestock-human interfaces. These findings support the need for integrated surveillance and control strategies to curb potential spillover and for the consideration of human M. bovis infection in SA patients with positive Ultra results.

Detection of Mycobacterium bovis in nasal swabs from communal goats (Capra hircus) in rural KwaZulu-Natal, South Africa.

Animal tuberculosis, caused by Mycobacterium bovis, presents a significant threat to both livestock industries and public health. Mycobacterium bovis tests rely on detecting antigen specific immune responses, which can be influenced by exposure to non-tuberculous mycobacteria, test technique, and duration and severity of infection. Despite advancements in direct M. bovis detection, mycobacterial culture remains the primary diagnostic standard. Recent efforts have explored culture-independent PCR-based methods for identifying mycobacterial DNA in respiratory samples. This study aimed to detect M. bovis in nasal swabs from goats (Capra hircus) cohabiting with M. bovis-infected cattle in KwaZulu-Natal, South Africa. Nasal swabs were collected from 137 communal goats exposed to M. bovis-positive cattle and 20 goats from a commercial dairy herd without M. bovis history. Swabs were divided into three aliquots for analysis. The first underwent GeneXpert® MTB/RIF Ultra assay (Ultra) screening. DNA from the second underwent mycobacterial genus-specific PCR and Sanger sequencing, while the third underwent mycobacterial culture followed by PCR and sequencing. Deep sequencing identified M. bovis DNA in selected Ultra-positive swabs, confirmed by region-of-difference (RD) PCR. Despite no other evidence of M. bovis infection, viable M. bovis was cultured from three communal goat swabs, confirmed by PCR and sequencing. Deep sequencing of DNA directly from swabs identified M. bovis in the same culture-positive swabs and eight additional communal goats. No M. bovis was found in commercial dairy goats, but various NTM species were detected. This highlights the risk of M. bovis exposure or infection in goats sharing pastures with infected cattle. Rapid Ultra screening shows promise for selecting goats for further M. bovis testing. These techniques may enhance M. bovis detection in paucibacillary samples and serve as valuable research tools.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing.

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

The MAGMA pipeline for comprehensive genomic analyses of clinical Mycobacterium tuberculosis samples.

BACKGROUND: Whole genome sequencing (WGS) holds great potential for the management and control of tuberculosis. Accurate analysis of samples with low mycobacterial burden, which are characterized by low (<20x) coverage and high (>40%) levels of contamination, is challenging. We created the MAGMA (Maximum Accessible Genome for Mtb Analysis) bioinformatics pipeline for analysis of clinical Mtb samples. METHODS AND RESULTS: High accuracy variant calling is achieved by using a long seedlength during read mapping to filter out contaminants, variant quality score recalibration with machine learning to identify genuine genomic variants, and joint variant calling for low Mtb coverage genomes. MAGMA automatically generates a standardized and comprehensive output of drug resistance information and resistance classification based on the WHO catalogue of Mtb mutations. MAGMA automatically generates phylogenetic trees with drug resistance annotations and trees that visualize the presence of clusters. Drug resistance and phylogeny outputs from sequencing data of 79 primary liquid cultures were compared between the MAGMA and MTBseq pipelines. The MTBseq pipeline reported only a proportion of the variants in candidate drug resistance genes that were reported by MAGMA. Notable differences were in structural variants, variants in highly conserved rrs and rrl genes, and variants in candidate resistance genes for bedaquiline, clofazmine, and delamanid. Phylogeny results were similar between pipelines but only MAGMA visualized clusters. CONCLUSION: The MAGMA pipeline could facilitate the integration of WGS into clinical care as it generates clinically relevant data on drug resistance and phylogeny in an automated, standardized, and reproducible manner.

"Spotting" Mycobacterium bovis infection in leopards (Panthera pardus) - novel application of diagnostic tools.

Background: Mycobacterium bovis (M. bovis) is the causative agent of animal tuberculosis (TB) which poses a threat to many of South Africa's most iconic wildlife species, including leopards (Panthera pardus). Due to limited tests for wildlife, the development of accurate ante-mortem tests for TB diagnosis in African big cat populations is urgently required. The aim of this study was to evaluate currently available immunological assays for their ability to detect M. bovis infection in leopards. Methods: Leopard whole blood (n=19) was stimulated using the QuantiFERON Gold Plus In-Tube System (QFT) to evaluate cytokine gene expression and protein production, along with serological assays. The GeneXpert® MTB/RIF Ultra (GXU®) qPCR assay, mycobacterial culture, and speciation by genomic regions of difference PCR, was used to confirm M. bovis infection in leopards. Results: Mycobacterium bovis infection was confirmed in six leopards and individuals that were tuberculin skin test (TST) negative were used for comparison. The GXU® assay was positive using all available tissue homogenates (n=5) from M. bovis culture positive animals. Mycobacterium bovis culture-confirmed leopards had greater antigen-specific responses, in the QFT interferon gamma release assay, CXCL9 and CXCL10 gene expression assays, compared to TST-negative individuals. One M. bovis culture-confirmed leopard had detectable antibodies using the DPP® Vet TB assay. Conclusion: Preliminary results demonstrated that immunoassays and TST may be potential tools to identify M. bovis-infected leopards. The GXU® assay provided rapid direct detection of infected leopards. Further studies should aim to improve TB diagnosis in wild felids, which will facilitate disease surveillance and screening.

The Detection of Mutations and Genotyping of Drug-Resistant Mycobacterium tuberculosis Strains Isolated from Patients in the Rural Eastern Cape Province.

Drug-resistant tuberculosis (DR-TB) is still a major public health concern in South Africa. Mutations in M. tuberculosis can cause varying levels of phenotypic resistance to anti-TB medications. There have been no prior studies on gene mutations and the genotyping of DR-TB in the rural Eastern Cape Province; hence, we aimed to identify DR-TB mutations, genetic diversity, and allocated lineages among patients in this area. Using Xpert® MTB/RIF, we assessed the rifampin resistance of sputum samples collected from 1157 patients suspected of having tuberculosis. GenoType MTBDR plus VER 2.0 was used for the detection of mutations causing resistance to anti-TB medications. The next step was to spoligotype 441 isolates. The most prevalent rifampin resistance-conferring mutations were in rpoB codon S531L in INH-resistant strains; the katG gene at codon S315TB and the inhA gene at codon C-15TB had the most mutations; 54.5% and 24.7%, respectively. In addition, 24.6% of strains showed mutations in both the rpoB and inhA genes, while 69.9% of strains showed mutations in both the katG and rpoB genes. Heteroresistance was seen in 17.9% of all cases in the study. According to spoligotyping analysis, Beijing families predominated. Investigation of the evolutionary lineages of M. tuberculosis isolates can be carried out using the information provided by the study's diversity of mutations. In locations wherein these mutations have been discovered, decision-making regarding the standardization of treatment regimens or individualized treatment may be aided by the detection frequency of rpoB, katG, and inhA mutations in various study areas.

Site-directed mutagenesis of Mycobacterium tuberculosis and functional validation to investigate potential bedaquiline resistance-causing mutations.

Molecular detection of bedaquiline resistant tuberculosis is challenging as only a small proportion of mutations in candidate bedaquiline resistance genes have been statistically associated with phenotypic resistance. We introduced two mutations, atpE Ile66Val and Rv0678 Thr33Ala, in the Mycobacterium tuberculosis H37Rv reference strain using homologous recombineering or recombination to investigate the phenotypic effect of these mutations. The genotype of the resulting strains was confirmed by Sanger- and whole genome sequencing, and bedaquiline susceptibility was assessed by minimal inhibitory concentration (MIC) assays. The impact of the mutations on protein stability and interactions was predicted using mutation Cutoff Scanning Matrix (mCSM) tools. The atpE Ile66Val mutation did not elevate the MIC above the critical concentration (MIC 0.25-0.5 µg/ml), while the MIC of the Rv0678 Thr33Ala mutant strains (> 1.0 µg/ml) classifies the strain as resistant, confirming clinical findings. In silico analyses confirmed that the atpE Ile66Val mutation minimally disrupts the bedaquiline-ATP synthase interaction, while the Rv0678 Thr33Ala mutation substantially affects the DNA binding affinity of the MmpR transcriptional repressor. Based on a combination of wet-lab and computational methods, our results suggest that the Rv0678 Thr33Ala mutation confers resistance to BDQ, while the atpE Ile66Val mutation does not, but definite proof can only be provided by complementation studies given the presence of secondary mutations.

The impact of genotype on the phenotype of Mycobacterium tuberculosis ΔsufR mutants.

Iron-sulphur (FeS) cluster biogenesis is a tightly regulated process in vivo. In Mycobacterium tuberculosis (Mtb), SufR functions as a transcriptional repressor of the operon encoding the primary FeS cluster biogenesis system. Previously, three independently isolated mutants (ΔRv1460stop_1.19, ΔRv1460stop _5.19 and ΔRv1460stop _5.20) harbouring the same deletion in sufR, displayed different growth kinetics in OADC supplemented 7H9 media. To investigate this discrepancy, we performed whole genome sequencing of the 3 mutants and the wild-type progenitor. Single nucleotide polymorphisms (SNPs) were identified in 3 genes in the ΔRv1460stop_1.19 mutant and one gene in the ΔRv1460stop_5.20 mutant. Phenotyping of the ΔRv1460stop_5.19 mutant, which had no additional SNPs, revealed increased susceptibility to clofazimine, DMNQ and menadione, while uptake and survival in THP-1 cells were not significantly different from the wild-type strain. Given that these results differ from those reported for other sufR deletion mutants (ΔSufRMTB and MtbΔSufR), they suggest that the position of the sufR deletion and the genotype of the progenitor strain impact the resulting phenotype.

Effect of compensatory evolution in the emergence and transmission of rifampicin-resistant Mycobacterium tuberculosis in Cape Town, South Africa: a genomic epidemiology study.

BACKGROUND: Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. METHODS: We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains. FINDINGS: Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors. INTERPRETATION: Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented. FUNDING: Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030_188888, CRSII5_177163, and IZLSZ3_170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council.

Point-of-care C-reactive protein and Xpert MTB/RIF Ultra for tuberculosis screening and diagnosis in unselected antiretroviral therapy initiators: a prospective diagnostic accuracy study.

Background: Tuberculosis (TB), a major cause of death in people living with HIV (PLHIV), remains challenging to diagnose. Diagnostic accuracy data are lacking for promising triage tests, such as C-reactive protein (CRP), and confirmatory tests, such as sputum and urine Xpert MTB/RIF Ultra (Ultra), and urine LAM, without prior symptom selection. Methods: 897 PLHIV initiating antiretroviral therapy were consecutively recruited in settings with high TB incidence, irrespective of symptoms. Participants were offered sputum induction, with a liquid culture reference standard. First, we evaluated point-of-care CRP testing on blood, compared to the World Health Organization (WHO)-recommended four-symptom screen (W4SS) for triage (n=800). Second, we evaluated Xpert MTB/RIF Ultra (Ultra) versus Xpert MTB/RIF (Xpert) for sputum-based confirmatory testing (n=787), with or without sputum induction. Third, we evaluated Ultra and Determine LF-LAM for urine-based confirmatory testing (n=732). Findings: CRP and number of W4SS symptoms had areas under the receiver operator characteristic curve of 0.78 (95% confidence interval 0.73, 0.83) and 0.70 (0.64, 0.75), respectively. For triage, CRP (≥10 mg/l) has similar sensitivity to W4SS [77% (68, 85) vs. 77% (68, 85); p>0.999] but higher specificity [64% (61, 68) vs. 48% (45, 52); p<0.001]; reducing unnecessary confirmatory testing by 138 per 1000 people and the number-needed-to-test from 6.91 (6.25, 7.81) to 4.87 (4.41, 5.51). Using sputum, which required induction in 31% (24, 39) of people, Ultra had higher sensitivity than Xpert [71% (61, 80) vs. 56% (46, 66); p<0.001] but lower specificity [98% (96, 100) vs. 99% (98, 100); p<0.001]. The proportion of people with ≥1 positive confirmatory result detected by Ultra increased from 45% (26, 64) to 66% (46, 82) when induction was done. Programmatically-done haemoglobin, triage test combinations, and urine tests showed comparatively worse performance. Interpretation: Among ART-initiators in a high burden setting, CRP is a more specific triage test than W4SS. Sputum induction improves yield. Sputum Ultra is a more accurate confirmatory test than Xpert. Funding: SAMRC (MRC-RFA-IFSP-01-2013), EDCTP2 (SF1401, OPTIMAL DIAGNOSIS), NIH/NIAD (U01AI152087). Research in context: Evidence before this study: Novel triage and confirmatory tests are urgently needed for TB, especially in key risk groups like PLHIV. Many TB cases do not meet World Health Organization (WHO)-recommended four-symptom screen (W4SS) criteria despite accounting for significant transmission and morbidity. W4SS also lacks specificity, which makes onward referral of triage-positive people for expensive confirmatory testing inefficient and hampers diagnostic scale-up. Alternative triage approaches like CRP have promise, but have comparatively little data in ART-initiators, especially when done without syndromic preselection and using point-of-care (POC) tools. After triage, confirmatory testing can be challenging due to sputum scarcity and paucibacillary early-stage disease. Next generation WHO-endorsed rapid molecular tests (including Xpert MTB/RIF Ultra; Ultra) are a standard-of-care for confirmatory testing. However, there are no supporting data in ART-initiators, among whom Ultra may offer large sensitivity gains over predecessors like Xpert MTB/RIF (Xpert). The added value of sputum induction to augment diagnostic sampling for confirmatory testing is also unclear. Lastly, the performance of urine tests (Ultra, Determine LF-LAM) in this population requires more data.Added value of this study: We evaluated repurposed and new tests for triage and confirmatory testing using a rigorous microbiological reference standard in a highly vulnerable high-priority patient population (ART-initiators) regardless of symptoms and ability to naturally expectorate sputum. We showed POC CRP triage is feasible, performs better than W4SS, and that combinations of different triage approaches offer no advantages over CRP alone. Sputum Ultra has superior sensitivity to Xpert; often detecting W4SS-negative TB. Furthermore, without induction, confirmatory sputum-based testing would not be possible in a third of people. Urine tests had poor performance. This study contributed unpublished data to systematic reviews and meta-analyses used by the WHO to inform global policy supporting use of CRP triage and Ultra in PLHIV.Implication of all the available evidence: POC CRP triage testing is feasible and superior to W4SS and, together with sputum induction in people who triage CRP-positive should, after appropriate cost and implementation research, be considered for roll-out in ART-initiators in high burden settings. Such people should be offered Ultra, which outperforms Xpert.

The Mycobacterium smegmatis HesB Protein, MSMEG_4272, Is Required for In Vitro Growth and Iron Homeostasis.

A-type carrier (ATC) proteins are proposed to function in the biogenesis of Fe-S clusters, although their exact role remains controversial. The genome of Mycobacterium smegmatis encodes a single ATC protein, MSMEG_4272, which belongs to the HesB/YadR/YfhF family of proteins. Attempts to generate an MSMEG_4272 deletion mutant by two-step allelic exchange were unsuccessful, suggesting that the gene is essential for in vitro growth. CRISPRi-mediated transcriptional knock-down of MSMEG_4272 resulted in a growth defect under standard culture conditions, which was exacerbated in mineral-defined media. The knockdown strain displayed reduced intracellular iron levels under iron-replete conditions and increased susceptibility to clofazimine, 2,3-dimethoxy-1,4-naphthoquinone (DMNQ), and isoniazid, while the activity of the Fe-S containing enzymes, succinate dehydrogenase, and aconitase were not affected. This study suggests that MSMEG_4272 plays a role in the regulation of intracellular iron levels and is required for in vitro growth of M. smegmatis, particularly during exponential growth.

Microfluidic Capture of Mycobacterium tuberculosis from Clinical Samples for Culture-Free Whole-Genome Sequencing.

Mycobacterium tuberculosis whole-genome sequencing (WGS) is a powerful tool as it can provide data on population diversity, drug resistance, disease transmission, and mixed infections. Successful WGS is still reliant on high concentrations of DNA obtained through M. tuberculosis culture. Microfluidics technology plays a valuable role in single-cell research but has not yet been assessed as a bacterial enrichment strategy for culture-free WGS of M. tuberculosis. In a proof-of-principle study, we evaluated the use of Capture-XT, a microfluidic lab-on-chip cleanup and pathogen concentration platform to enrich M. tuberculosis bacilli from clinical sputum specimens for downstream DNA extraction and WGS. Three of the four (75%) samples processed by the microfluidics application passed the library preparation quality control, compared to only one of the four (25%) samples not enriched by the microfluidics M. tuberculosis capture application. WGS data were of sufficient quality, with mapping depth of ≥25× and 9 to 27% of reads mapping to the reference genome. These results suggest that microfluidics-based M. tuberculosis cell capture might be a promising method for M. tuberculosis enrichment in clinical sputum samples, which could facilitate culture-free M. tuberculosis WGS. IMPORTANCE Diagnosis of tuberculosis is effective using molecular methods; however, a comprehensive characterization of the resistance profile of Mycobacterium tuberculosis often requires culturing and phenotypic drug susceptibility testing or culturing followed by whole-genome sequencing (WGS). The phenotypic route can take anywhere from 1 to >3 months to result, by which point the patient may have acquired additional drug resistance. The WGS route is a very attractive option; however, culturing is the rate-limiting step. In this original article, we provide proof-of-principle evidence that microfluidics-based cell capture can be used on high-bacillary-load clinical samples for culture-free WGS.

Variation in missed doses and reasons for discontinuation of anti-tuberculosis drugs during hospital treatment for drug-resistant tuberculosis in South Africa.

BACKGROUND: Updated World Health Organization (WHO) treatment guidelines prioritize all-oral drug-resistant tuberculosis (DR-TB) regimens. Several poorly tolerated drugs, such as amikacin and para-aminosalicylic acid (PAS), remain treatment options for DR-TB in WHO-recommended longer regimens as Group C drugs. Incomplete treatment with anti-TB drugs increases the risk of treatment failure, relapse, and death. We determined whether missed doses of individual anti-TB drugs, and reasons for their discontinuation, varied in closely monitored hospital settings prior to the 2020 WHO DR-TB treatment guideline updates. METHODS: We collected retrospective data on adult patients with microbiologically confirmed DR-TB between 2008 and 2015 who were selected for a study of acquired drug resistance in the Western Cape Province of South Africa. Medical records through mid-2017 were reviewed. Patients received directly observed treatment during hospitalization at specialized DR-TB hospitals. Incomplete treatment with individual anti-TB drugs, defined as the failure to take medication as prescribed, regardless of reason, was determined by comparing percent missed doses, stratified by HIV status and DR-TB regimen. We applied a generalized mixed effects model. RESULTS: Among 242 patients, 131 (54%) were male, 97 (40%) were living with HIV, 175 (72%) received second-line treatment prior to first hospitalization, and 191 (79%) died during the study period. At initial hospitalization, 134 (55%) patients had Mycobacterium tuberculosis with resistance to rifampicin and isoniazid (multidrug-resistant TB [MDR-TB]) without resistance to ofloxacin or amikacin, and 102 (42%) had resistance to ofloxacin and/or amikacin. Most patients (129 [53%]) had multiple hospitalizations and DST changes occurred in 146 (60%) by the end of their last hospital discharge. Incomplete treatment was significantly higher for amikacin (18%), capreomycin (18%), PAS (17%) and kanamycin (16%) than other DR-TB drugs (P<0.001), including ethionamide (8%), moxifloxacin (7%), terizidone (7%), ethambutol (7%), and pyrazinamide (6%). Among the most frequently prescribed drugs, second-line injectables had the highest rates of discontinuation for adverse events (range 0.56-1.02 events per year follow-up), while amikacin, PAS and ethionamide had the highest rates of discontinuation for patient refusal (range 0.51-0.68 events per year follow-up). Missed doses did not differ according to HIV status or anti-TB drug combinations. CONCLUSION: We found that incomplete treatment for second-line injectables and PAS during hospitalization was higher than for other anti-TB drugs. To maximize treatment success, interventions to improve person-centered care and mitigate adverse events may be necessary in cases when PAS or amikacin (2020 WHO recommended Group C drugs) are needed.

A machine-learning based model for automated recommendation of individualized treatment of rifampicin-resistant tuberculosis.

Background: Rifampicin resistant tuberculosis remains a global health problem with almost half a million new cases annually. In high-income countries patients empirically start a standardized treatment regimen, followed by an individualized regimen guided by drug susceptibility test (DST) results. In most settings, DST information is not available or is limited to isoniazid and fluoroquinolones. Whole genome sequencing could more accurately guide individualized treatment as the full drug resistance profile is obtained with a single test. Whole genome sequencing has not reached its full potential for patient care, in part due to the complexity of translating a resistance profile into the most effective individualized regimen. Methods: We developed a treatment recommender clinical decision support system (CDSS) and an accompanying web application for user-friendly recommendation of the optimal individualized treatment regimen to a clinician. Results: Following expert stakeholder meetings and literature review, nine drug features and 14 treatment regimen features were identified and quantified. Using machine learning, a model was developed to predict the optimal treatment regimen based on a training set of 3895 treatment regimen-expert feedback pairs. The acceptability of the treatment recommender CDSS was assessed as part of a clinical trial and in a routine care setting. Within the clinical trial setting, all patients received the CDSS recommended treatment. In 8 of 20 cases, the initial recommendation was recomputed because of stock out, clinical contra-indication or toxicity. In routine care setting, physicians rejected the treatment recommendation in 7 out of 15 cases because it deviated from the national TB treatment guidelines. A survey indicated that the treatment recommender CDSS is easy to use and useful in clinical practice but requires digital infrastructure support and training. Conclusions: Our findings suggest that global implementation of the novel treatment recommender CDSS holds the potential to improve treatment outcomes of rifampicin resistant tuberculosis.

Whole genome sequencing improves the discrimination between Mycobacterium bovis strains on the southern border of Kruger National Park, South Africa.

Background: Mycobacterium bovis forms part of the Mycobacterium tuberculosis complex and has an extensive host range and zoonotic potential. Various genotyping methods (e.g., spoligotyping) have been used to describe the molecular epidemiology of M. bovis. Advances in whole genome sequencing (WGS) have increased resolution to enable detection of genomic variants to the level of single nucleotide polymorphisms. This is especially relevant to One Health research on tuberculosis which benefits by being able to use WGS to identify epidemiologically linked cases, especially recent transmission. The use of WGS in molecular epidemiology has been extensively used in humans and cattle but is limited in wildlife. This approach appears to overcome the limitations of conventional genotyping methods due to lack of genetic diversity in M. bovis. Methods: This pilot study investigated the spoligotype and WGS of M. bovis isolates (n = 7) from wildlife in Marloth Park (MP) and compared these with WGS data from other South African M. bovis isolates. In addition, the greater resolution of WGS was used to explore the phylogenetic relatedness of M. bovis isolates in neighbouring wildlife populations. Results: The phylogenetic analyses showed the closest relatives to the seven isolates from MP were isolates from wildlife in Kruger National Park (KNP), which shares a border with MP. However, WGS data indicated that the KNP and MP isolates formed two distinct clades, even though they had similar spoligotypes and identical in silico genetic regions of difference profiles. Conclusions: Mycobacterium bovis isolates from MP were hypothesized to be directly linked to KNP wildlife, based on spoligotyping. However, WGS indicated more complex epidemiology. The presence of two distinct clades which were genetically distinct (SNP distance of 19-47) and suggested multiple transmission events. Therefore, WGS provided new insight into the molecular epidemiology of the M. bovis isolates from MP and their relationship to isolates from KNP. This approach will facilitate greater understanding of M. bovis transmission at wildlife-livestock-human interfaces and advances One Health research on tuberculosis, especially across different host species.

Methods for handling missing data in serially sampled sputum specimens for mycobacterial culture conversion calculation.

BACKGROUND: The occurrence and timing of mycobacterial culture conversion is used as a proxy for tuberculosis treatment response. When researchers serially sample sputum during tuberculosis studies, contamination or missed visits leads to missing data points. Traditionally, this is managed by ignoring missing data or simple carry-forward techniques. Statistically advanced multiple imputation methods potentially decrease bias and retain sample size and statistical power. METHODS: We analyzed data from 261 participants who provided weekly sputa for the first 12 weeks of tuberculosis treatment. We compared methods for handling missing data points in a longitudinal study with a time-to-event outcome. Our primary outcome was time to culture conversion, defined as two consecutive weeks with no Mycobacterium tuberculosis growth. Methods used to address missing data included: 1) available case analysis, 2) last observation carried forward, and 3) multiple imputation by fully conditional specification. For each method, we calculated the proportion culture converted and used survival analysis to estimate Kaplan-Meier curves, hazard ratios, and restricted mean survival times. We compared methods based on point estimates, confidence intervals, and conclusions to specific research questions. RESULTS: The three missing data methods lead to differences in the number of participants achieving conversion; 78 (32.8%) participants converted with available case analysis, 154 (64.7%) converted with last observation carried forward, and 184 (77.1%) converted with multiple imputation. Multiple imputation resulted in smaller point estimates than simple approaches with narrower confidence intervals. The adjusted hazard ratio for smear negative participants was 3.4 (95% CI 2.3, 5.1) using multiple imputation compared to 5.2 (95% CI 3.1, 8.7) using last observation carried forward and 5.0 (95% CI 2.4, 10.6) using available case analysis. CONCLUSION: We showed that accounting for missing sputum data through multiple imputation, a statistically valid approach under certain conditions, can lead to different conclusions than naïve methods. Careful consideration for how to handle missing data must be taken and be pre-specified prior to analysis. We used data from a TB study to demonstrate these concepts, however, the methods we described are broadly applicable to longitudinal missing data. We provide valuable statistical guidance and code for researchers to appropriately handle missing data in longitudinal studies.