Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease.

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Bridging the gap between molecular and genomic epidemiology in tuberculosis: inferring MIRU-VNTR patterns from genomic data.

The transition from MIRU-VNTR-based epidemiology studies in tuberculosis (TB) to genomic epidemiology has transformed how we track transmission. However, short-read sequencing is poor at analyzing repetitive regions such as the MIRU-VNTR loci. This causes a gap between the new genomic data and the large amount of information stored in historical databases. Long-read sequencing could bridge this knowledge gap by allowing analysis of repetitive regions. However, the feasibility of extracting MIRU-VNTRs from long reads and linking them to historical data has not been evaluated. In our study, an in silico arm, consisting of inference of MIRU patterns from long-read sequences (using MIRUReader program), was compared with an experimental arm, involving standard amplification and fragment sizing. We analyzed overall performance on 39 isolates from South Africa and confirmed reproducibility in a sample enriched with 62 clustered cases from Spain. Finally, we ran 25 consecutive incident cases, demonstrating the feasibility of correctly assigning new clustered/orphan cases by linking data inferred from genomic analysis to MIRU-VNTR databases. Of the 3,024 loci analyzed, only 11 discrepancies (0.36%) were found between the two arms: three attributed to experimental error and eight to misassigned alleles from long-read sequencing. A second round of analysis of these discrepancies resulted in agreement between the experimental and in silico arms in all but one locus. Adjusting the MIRUReader program code allowed us to flag potential in silico misassignments due to suboptimal coverage or unfixed double alleles. Our study indicates that long-read sequencing could help address potential chronological and geographical gaps arising from the transition from molecular to genomic epidemiology of tuberculosis. IMPORTANCE: The transition from molecular epidemiology in tuberculosis (TB), based on the analysis of repetitive regions (VNTR-based genotyping), to genomic epidemiology transforms in the precision with which we track transmission. However, short-read sequencing, the most common method for performing genomic analysis, is poor at analyzing repetitive regions. This means that we face a gap between the new genomic data and the large amount of information stored in historical databases, which is also an obstacle to cross-national surveillance involving settings where only molecular data are available. Long-read sequencing could help bridge this knowledge gap by allowing analysis of repetitive regions. Our study demonstrates that MIRU-VNTR patterns can be successfully inferred from long-read sequences, allowing the correct assignment of new cases as clustered/orphan by linking new data extracted from genomic analysis to historical MIRU-VNTR databases. Our data may provide a starting point for bridging the knowledge gap between the molecular and genomic eras in tuberculosis epidemiology.

The impact of alcohol and illicit substance use on the pharmacokinetics of first-line TB drugs.

BACKGROUND: In South Africa, an estimated 11% of the population have high alcohol use, a major risk factor for TB. Alcohol and other substance use are also associated with poor treatment response, with a potential mechanism being altered TB drug pharmacokinetics. OBJECTIVES: To investigate the impact of alcohol and illicit substance use on the pharmacokinetics of first-line TB drugs in participants with pulmonary TB. METHODS: We prospectively enrolled participants ≥15 years old, without HIV, and initiating drug-susceptible TB treatment in Worcester, South Africa. Alcohol use was measured via self-report and blood biomarkers. Other illicit substances were captured through a urine drug test. Plasma samples were drawn 1 month into treatment pre-dose, and 1.5, 3, 5 and 8 h post-dose. Non-linear mixed-effects modelling was used to describe the pharmacokinetics of rifampicin, isoniazid, pyrazinamide and ethambutol. Alcohol and drug use were tested as covariates. RESULTS: The study included 104 participants, of whom 70% were male, with a median age of 37 years (IQR 27-48). Alcohol use was high, with 42% and 28% of participants having moderate and high alcohol use, respectively. Rifampicin and isoniazid had slightly lower pharmacokinetics compared with previous reports, whereas pyrazinamide and ethambutol were consistent. No significant alcohol use effect was detected, other than 13% higher ethambutol clearance in participants with high alcohol use. Methaqualone use reduced rifampicin bioavailability by 19%. CONCLUSION: No clinically relevant effect of alcohol use was observed on the pharmacokinetics of first-line TB drugs, suggesting that poor treatment outcome is unlikely due to pharmacokinetic alterations. That methaqualone reduced rifampicin means dose adjustment may be beneficial.

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.

PneumoniaCheck, a novel aerosol collection device, permits capture of airborne Mycobacterium tuberculosis and characterisation of the cough aeromicrobiome in people with tuberculosis.

BACKGROUND: Tuberculosis (TB), a major cause of disease and antimicrobial resistance, is spread via aerosols. Aerosols have diagnostic potential and airborne-microbes other than Mycobacterium tuberculosis complex (MTBC) may influence transmission. We evaluated whether PneumoniaCheck (PMC), a commercial aerosol collection device, captures MTBC and the aeromicrobiome of people with TB. METHODS: PMC was done in sputum culture-positive people (≥ 30 forced coughs each, n = 16) pre-treatment and PMC air reservoir (bag, corresponding to upper airways) and filter (lower airways) washes underwent Xpert MTB/RIF Ultra (Ultra) and 16S rRNA gene sequencing (sequencing also done on sputum). In a subset (n = 6), PMC microbiota (bag, filter) was compared to oral washes and bronchoalveolar lavage fluid (BALF). FINDINGS: 54% (7/13) bags and 46% (6/14) filters were Ultra-positive. Sequencing read counts and microbial diversity did not differ across bags, filters, and sputum. However, microbial composition in bags (Sphingobium-, Corynebacterium-, Novosphingobium-enriched) and filters (Mycobacterium-, Sphingobium-, Corynebacterium-enriched) each differed vs. sputum. Furthermore, sequencing only detected Mycobacterium in bags and filters but not sputum. In the subset, bag and filter microbial diversity did not differ vs. oral washes or BALF but microbial composition differed. Bags vs. BALF were Sphingobium-enriched and Mycobacterium-, Streptococcus-, and Anaerosinus-depleted (Anaerosinus also depleted in filters vs. BALF). Compared to BALF, none of the aerosol-enriched taxa were enriched in oral washes or sputum. INTERPRETATION: PMC captures aerosols with Ultra-detectable MTBC and MTBC is more detectable in aerosols than sputum by sequencing. The aeromicrobiome is distinct from sputum, oral washes and BALF and contains differentially-enriched lower respiratory tract microbes.

More than Mycobacterium tuberculosis: site-of-disease microbial communities, and their functional and clinical profiles in tuberculous lymphadenitis.

BACKGROUND: Lymphadenitis is the most common extrapulmonary tuberculosis (EPTB) manifestation. The microbiome is important to human health but uninvestigated in EPTB. We profiled the site-of-disease lymph node microbiome in tuberculosis lymphadenitis (TBL). METHODS: Fine-needle aspiration biopsies were collected from 158 pretreatment presumptive TBL patients in Cape Town, South Africa. 16S Illumina MiSeq rRNA gene sequencing was done. RESULTS: We analysed 89 definite TBLs (dTBLs) and 61 non-TBLs (nTBLs), which had similar α- but different ß-diversities (p=0.001). Clustering identified five lymphotypes prior to TB status stratification: Mycobacterium-dominant, Prevotella-dominant and Streptococcus-dominant lymphotypes were more frequent in dTBLs whereas a Corynebacterium-dominant lymphotype and a fifth lymphotype (no dominant taxon) were more frequent in nTBLs. When restricted to dTBLs, clustering identified a Mycobacterium-dominant lymphotype with low α-diversity and non-Mycobacterium-dominated lymphotypes (termed Prevotella-Corynebacterium, Prevotella-Streptococcus). The Mycobacterium dTBL lymphotype was associated with HIV-positivity and features characteristic of severe lymphadenitis (eg, larger nodes). dTBL microbial communities were enriched with potentially proinflammatory microbial short-chain fatty acid metabolic pathways (propanoate, butanoate) vs nTBLs. 11% (7/61) of nTBLs had Mycobacterium reads BLAST-confirmed as Mycobacterium tuberculosis complex. CONCLUSIONS: TBL at the site-of-disease is not microbially homogeneous. Distinct microbial community clusters exist that, in our setting, are associated with different clinical characteristics, and immunomodulatory potentials. Non-Mycobacterium-dominated dTBL lymphotypes, which contain taxa potentially targeted by TB treatment, were associated with milder, potentially earlier stage disease. These investigations lay foundations for studying the microbiome's role in lymphatic TB. The long-term clinical significance of these lymphotypes requires prospective validation.

An All-Oral 6-Month Regimen for Multidrug-Resistant Tuberculosis: A Multicenter, Randomized Controlled Clinical Trial (the NExT Study).

Rationale: Improving treatment outcomes while reducing drug toxicity and shortening the treatment duration to ∼6 months remains an aspirational goal for the treatment of multidrug-resistant/rifampicin-resistant tuberculosis (MDR/RR-TB). Objectives: To conduct a multicenter randomized controlled trial in adults with MDR/RR-TB (i.e., without resistance to fluoroquinolones or aminoglycosides). Methods: Participants were randomly assigned (1:1 ratio) to a ∼6-month all-oral regimen that included levofloxacin, bedaquiline, and linezolid, or the standard-of-care (SOC) ⩾9-month World Health Organization (WHO)-approved injectable-based regimen. The primary endpoint was a favorable WHO-defined treatment outcome (which mandates that prespecified drug substitution is counted as an unfavorable outcome) 24 months after treatment initiation. The trial was stopped prematurely when bedaquiline-based therapy became the standard of care in South Africa. Measurements and Main Results: In total, 93 of 111 randomized participants (44 in the comparator arm and 49 in the interventional arm) were included in the modified intention-to-treat analysis; 51 (55%) were HIV coinfected (median CD4 count, 158 cells/ml). Participants in the intervention arm were 2.2 times more likely to experience a favorable 24-month outcome than participants in the SOC arm (51% [25 of 49] vs. 22.7% [10 of 44]; risk ratio, 2.2 [1.2-4.1]; P = 0.006). Toxicity-related drug substitution occurred more frequently in the SOC arm (65.9% [29 of 44] vs. 34.7% [17 of 49]; P = 0.001)], 82.8% (24 of 29) owing to kanamycin (mainly hearing loss; replaced by bedaquiline) in the SOC arm, and 64.7% (11 of 17) owing to linezolid (mainly anemia) in the interventional arm. Adverse event-related treatment discontinuation in the safety population was more common in the SOC arm (56.4% [31 of 55] vs. 32.1% [17 of 56]; P = 0.007). However, grade 3 adverse events were more common in the interventional arm (55.4% [31 of 56] vs. 32.7 [18 of 55]; P = 0.022). Culture conversion was significantly better in the intervention arm (hazard ratio, 2.6 [1.4-4.9]; P = 0.003) after censoring those with bedaquiline replacement in the SOC arm (and this pattern remained consistent after censoring for drug replacement in both arms; P = 0.01). Conclusions: Compared with traditional injectable-containing regimens, an all-oral 6-month levofloxacin, bedaquiline, and linezolid-containing MDR/RR-TB regimen was associated with a significantly improved 24-month WHO-defined treatment outcome (predominantly owing to toxicity-related drug substitution). However, drug toxicity occurred frequently in both arms. These findings inform strategies to develop future regimens for MDR/RR-TB.Clinical trial registered with www.clinicaltrials.gov (NCT02454205).

Variants in Bedaquiline-Candidate-Resistance Genes: Prevalence in Bedaquiline-Naive Patients, Effect on MIC, and Association with Mycobacterium tuberculosis Lineage.

Studies have shown that variants in bedaquiline-resistance genes can occur in isolates from bedaquiline-naive patients. We assessed the prevalence of variants in all bedaquiline-candidate-resistance genes in bedaquiline-naive patients, investigated the association between these variants and lineage, and the effect on phenotype. We used whole-genome sequencing to identify variants in bedaquiline-resistance genes in isolates from 509 bedaquiline treatment naive South African tuberculosis patients. A phylogenetic tree was constructed to investigate the association with the isolate lineage background. Bedaquiline MIC was determined using the UKMYC6 microtiter assay. Variants were identified in 502 of 509 isolates (98.6%), with the highest (85%) prevalence of variants in the Rv0676c (mmpL5) gene. We identified 36 unique variants, including 19 variants not reported previously. Only four isolates had a bedaquiline MIC equal to or above the epidemiological cut-off value of 0.25 µg/mL. Phylogenetic analysis showed that 14 of the 15 variants observed more than once occurred monophyletically in one Mycobacterium tuberculosis (sub)lineage. The bedaquiline MIC differed between isolates belonging to lineage 2 and 4 (Fisher's exact test, P = 0.0004). The prevalence of variants in bedaquiline-resistance genes in isolates from bedaquiline-naive patients is high, but very few (<2%) isolates were phenotypically resistant. We found an association between variants in bedaquiline resistance genes and Mycobacterium tuberculosis (sub)lineage, resulting in a lineage-dependent difference in bedaquiline phenotype. Future studies should investigate the impact of the presence of variants on bedaquiline-resistance acquisition and treatment outcome.

Nanopore Sequencing for Mycobacterium tuberculosis: a Critical Review of the Literature, New Developments, and Future Opportunities.

The next-generation, short-read sequencing technologies that generate comprehensive, whole-genome data with single nucleotide resolution have already advanced tuberculosis diagnosis, treatment, surveillance, and source investigation. Their high costs, tedious and lengthy processes, and large equipment remain major hurdles for research use in high tuberculosis burden countries and implementation into routine care. The portable next-generation sequencing devices developed by Oxford Nanopore Technologies (ONT) are attractive alternatives due to their long-read sequence capability, compact low-cost hardware, and continued improvements in accuracy and throughput. A systematic review of the published literature demonstrated limited uptake of ONT sequencing in tuberculosis research and clinical care. Of the 12 eligible articles presenting ONT sequencing data on at least one Mycobacterium tuberculosis sample, four addressed software development for long-read ONT sequencing data with potential applications for M. tuberculosis. Only eight studies presented results of ONT sequencing of M. tuberculosis, of which five performed whole-genome and three did targeted sequencing. Based on these findings, we summarize the standard processes, reflect on the current limitations of ONT sequencing technology, and the research needed to overcome the main hurdles. The low capital cost, portable nature and continued improvement in the performance of ONT sequencing make it an attractive option for sequencing for research and clinical care, but limited data are available on its application in the tuberculosis field. Important research investment is needed to unleash the full potential of ONT sequencing for tuberculosis research and care.

Comparative Performance of Genomic Methods for the Detection of Pyrazinamide Resistance and Heteroresistance in Mycobacterium tuberculosis.

Pyrazinamide is an important component of both drug-susceptible and drug-resistant tuberculosis treatment regimens. Although approximately 50% of rifampin-resistant isolates are also resistant to pyrazinamide, pyrazinamide susceptibility testing is not routinely performed due to the challenging nature of the assay. We investigated the diagnostic accuracy of genotypic and phenotypic methods and explored the occurrence of pyrazinamide heteroresistance. We assessed pyrazinamide susceptibility among 358 individuals enrolled in the South African EXIT-RIF cohort using Sanger and targeted deep sequencing (TDS) of the pncA gene, whole-genome sequencing (WGS), and phenotypic drug susceptibility testing. We calculated the diagnostic accuracy of the different methods and investigated the prevalence and clinical impact of pncA heteroresistance. True pyrazinamide susceptibility status was assigned to each isolate using the Köser classification and expert rules. We observed 100% agreement across genotypic methods for detection of pncA fixed mutations; only TDS confidently identified three isolates (0.8%) with minor variants. For the 355 (99.2%) isolates that could be assigned true pyrazinamide status with confidence, phenotypic DST had a sensitivity of 96.5% (95% confidence interval [CI], 93.8 to 99.3%) and specificity of 100% (95% CI, 100 to 100%), both Sanger sequencing and WGS had a sensitivity of 97.1% (95% CI, 94.6 to 99.6%) and specificity of 97.8% (95% CI, 95.7 to 99.9%), and TDS had sensitivity of 98.8% (95% CI, 97.2 to 100%) and specificity of 97.8% (95% CI, 95.7 to 99.9%). We demonstrate high sensitivity and specificity for pyrazinamide susceptibility testing among all assessed genotypic methods. The prevalence of pyrazinamide heteroresistance in Mycobacterium tuberculosis isolates was lower than that identified for other first-line drugs.

Whole-Genome Sequencing Has the Potential To Improve Treatment for Rifampicin-Resistant Tuberculosis in High-Burden Settings: a Retrospective Cohort Study.

Treatment of multidrug-resistant or rifampicin-resistant tuberculosis (MDR/RR-TB), although improved in recent years with shorter, more tolerable regimens, remains largely standardized and based on limited drug susceptibility testing (DST). More individualized treatment with expanded DST access is likely to improve patient outcomes. To assess the potential of TB drug resistance prediction based on whole-genome sequencing (WGS) to provide more effective treatment regimens, we applied current South African treatment recommendations to a retrospective cohort of MDR/RR-TB patients from Khayelitsha, Cape Town. Routine DST and clinical data were used to retrospectively categorize patients into a recommended regimen, either a standardized short regimen or a longer individualized regimen. Potential regimen changes were then described with the addition of WGS-derived DST. WGS data were available for 1274 MDR/RR-TB patient treatment episodes across 2008 to 2017. Among 834 patients initially eligible for the shorter regimen, 385 (46%) may have benefited from reduced drug dosage or removing ineffective drugs when WGS data were considered. A further 187 (22%) patients may have benefited from more effective adjusted regimens. Among 440 patients initially eligible for a longer individualized regimen, 153 (35%) could have been switched to the short regimen. Overall, 305 (24%) patients had MDR/RR-TB with second-line TB drug resistance, where the availability of WGS-derived DST would have allowed more effective treatment individualization. These data suggest considerable benefits could accrue from routine access to WGS-derived resistance prediction. Advances in culture-free sequencing and expansion of the reference resistance mutation catalogue will increase the utility of WGS resistance prediction.

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid.

OBJECTIVES: To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. METHODS: The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. RESULTS: We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. CONCLUSIONS: In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST.

High mortality among patients hospitalized for drug-resistant tuberculosis with acquired second-line drug resistance and high HIV prevalence.

OBJECTIVES: We compared mortality between HIV-positive and HIV-negative South African adults with drug-resistant tuberculosis (DR-TB) and high incidence of acquired second-line drug resistance. METHODS: We performed a retrospective review of DR-TB patients with serial second-line TB drug susceptibility tests (2008-2015) who were hospitalized at a specialized TB hospital. We used Kaplan-Meier analysis and Cox models to examine associations with mortality. RESULTS: Of 245 patients, the median age was 33 years, 54% were male and 40% were HIV-positive, 96% of whom had ever received antiretroviral therapy (ART). At initial drug resistance detection, 99% of patients had resistance to at least rifampicin and isoniazid, and 18% had second-line drug resistance (fluoroquinolones and/or injectable drugs). At later testing, 88% of patients had acquired additional second-line drug resistance. Patient-initiated treatment interruptions (> 2 months) occurred in 47%. Mortality was 79%. Those with HIV had a shorter time to death (p = 0.02; log-rank): median survival time from DR-TB treatment initiation was 2.44 years [95% confidence interval (CI): 2.09-3.15] versus 3.99 years (95% CI: 3.12-4.75) for HIV-negative patients. HIV-positive patients who received ART within 6 months before DR-TB treatment had a higher mortality hazard than HIV-negative patients [adjusted hazard ratio (aHR) ratio = 1.82, 95% CI: 1.21-2.74]. By contrast, HIV-positive patients who did not receive ART within 6 months before DR-TB treatment did not have a significantly higher mortality hazard than HIV-negative patients (aHR = 1.09; 95% CI: 0.72-1.65), although those on ART had lower median CD4 counts than those not on ART (157 vs. 281 cells/µL, respectively; p = 0.02). CONCLUSIONS: A very high incidence of acquired second-line drug resistance and high overall mortality were observed, reinforcing the need to reduce the risk of acquired resistance and for more effective treatment.

Localization of EccA3 at the growing pole in Mycobacterium smegmatis.

BACKGROUND: Bacteria require specialized secretion systems for the export of molecules into the extracellular space to modify their environment and scavenge for nutrients. The ESX-3 secretion system is required by mycobacteria for iron homeostasis. The ESX-3 operon encodes for one cytoplasmic component (EccA3) and five membrane components (EccB3 - EccE3 and MycP3). In this study we sought to identify the sub-cellular location of EccA3 of the ESX-3 secretion system in mycobacteria. RESULTS: Fluorescently tagged EccA3 localized to a single pole in the majority of Mycobacterium smegmatis cells and time-lapse fluorescent microscopy identified this pole as the growing pole. Deletion of ESX-3 did not prevent polar localization of fluorescently tagged EccA3, suggesting that EccA3 unipolar localization is independent of other ESX-3 components. Affinity purification - mass spectrometry was used to identify EccA3 associated proteins which may contribute to the localization of EccA3 at the growing pole. EccA3 co-purified with fatty acid metabolism proteins (FAS, FadA3, KasA and KasB), mycolic acid synthesis proteins (UmaA, CmaA1), cell division proteins (FtsE and FtsZ), and cell shape and cell cycle proteins (MurS, CwsA and Wag31). Secretion system related proteins Ffh, SecA1, EccA1, and EspI were also identified. CONCLUSIONS: Time-lapse microscopy demonstrated that EccA3 is located at the growing pole in M. smegmatis. The co-purification of EccA3 with proteins known to be required for polar growth, mycolic acid synthesis, the Sec secretion system (SecA1), and the signal recognition particle pathway (Ffh) also suggests that EccA3 is located at the site of active cell growth.

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.

A treatment recommender clinical decision support system for personalized medicine: method development and proof-of-concept for drug resistant tuberculosis.

BACKGROUND: Personalized medicine tailors care based on the patient's or pathogen's genotypic and phenotypic characteristics. An automated Clinical Decision Support System (CDSS) could help translate the genotypic and phenotypic characteristics into optimal treatment and thus facilitate implementation of individualized treatment by less experienced physicians. METHODS: We developed a hybrid knowledge- and data-driven treatment recommender CDSS. Stakeholders and experts first define the knowledge base by identifying and quantifying drug and regimen features for the prototype model input. In an iterative manner, feedback from experts is harvested to generate model training datasets, machine learning methods are applied to identify complex relations and patterns in the data, and model performance is assessed by estimating the precision at one, mean reciprocal rank and mean average precision. Once the model performance no longer iteratively increases, a validation dataset is used to assess model overfitting. RESULTS: We applied the novel methodology to develop a treatment recommender CDSS for individualized treatment of drug resistant tuberculosis as a proof of concept. Using input from stakeholders and three rounds of expert feedback on a dataset of 355 patients with 129 unique drug resistance profiles, the model had a 95% precision at 1 indicating that the highest ranked treatment regimen was considered appropriate by the experts in 95% of cases. Use of a validation data set however suggested substantial model overfitting, with a reduction in precision at 1 to 78%. CONCLUSION: Our novel and flexible hybrid knowledge- and data-driven treatment recommender CDSS is a first step towards the automation of individualized treatment for personalized medicine. Further research should assess its value in fields other than drug resistant tuberculosis, develop solid statistical approaches to assess model performance, and evaluate their accuracy in real-life clinical settings.

A multi-phenotype genome-wide association study of clades causing tuberculosis in a Ghanaian- and South African cohort.

Despite decades of research and advancements in diagnostics and treatment, tuberculosis remains a major public health concern. New computational methods are needed to interrogate the intersection of host- and bacterial genomes. Paired host genotype datum and infecting bacterial isolate information were analysed for associations using a multinomial logistic regression framework implemented in SNPTest. A cohort of 853 admixed South African participants and a Ghanaian cohort of 1359 participants were included. Two directly genotyped variants, namely rs529920 and rs41472447, were identified in the Ghanaian cohort as being statistically significantly associated with risk for infection with strains of different members of the MTBC. Thus, a multinomial logistic regression using paired host-pathogen data may prove valuable for investigating the complex relationships driving infectious disease.

Xpert Mycobacterium tuberculosis/Rifampicin-Detected Rifampicin Resistance is a Suboptimal Surrogate for Multidrug-resistant Tuberculosis in Eastern Democratic Republic of the Congo: Diagnostic and Clinical Implications.

BACKGROUND: Rifampicin (RIF) resistance is highly correlated with isoniazid (INH) resistance and used as proxy for multidrug-resistant tuberculosis (MDR-TB). Using MTBDRplus as a comparator, we evaluated the predictive value of Xpert MTB/RIF (Xpert)-detected RIF resistance for MDR-TB in eastern Democratic Republic of the Congo (DRC). METHODS: We conducted a cross-sectional study involving data from new or retreatment pulmonary adult TB cases evaluated between July 2013 and December 2016. Separate, paired sputa for smear microscopy and MTBDRplus were collected. Xpert testing was performed subject to the availability of Xpert cartridges on sample remnants after microscopy. RESULTS: Among 353 patients, 193 (54.7%) were previously treated and 224 (63.5%) were MTBDRplus TB positive. Of the 224, 43 (19.2%) were RIF monoresistant, 11 (4.9%) were INH monoresistant, 53 (23.7%) had MDR-TB, and 117 (52.2%) were RIF and INH susceptible. Overall, among the 96 samples detected by MTBDRplus as RIF resistant, 53 (55.2%) had MDR-TB. Xpert testing was performed in 179 (50.7%) specimens; among these, 163 (91.1%) were TB positive and 73 (44.8%) RIF resistant. Only 45/73 (61.6%) Xpert-identified RIF-resistant isolates had concomitant MTBDRplus-detected INH resistance. Xpert had a sensitivity of 100.0% (95% CI, 92.1-100.0) for detecting RIF resistance but a positive-predictive value of only 61.6% (95% CI, 49.5-72.8) for MDR-TB. The most frequent mutations associated with RIF and INH resistance were S531L and S315T1, respectively. CONCLUSIONS: In this high-risk MDR-TB study population, Xpert had low positive-predictive value for the presence of MDR-TB. Comprehensive resistance testing for both INH and RIF should be performed in this setting.

Decentralized Care for Rifampin-Resistant Tuberculosis, Western Cape, South Africa.

In 2011, South Africa implemented a policy to decentralize treatment for rifampin-resistant tuberculosis (TB) to reduce durations of hospitalization and enable local treatment. We assessed policy implementation in Western Cape Province, where services expanded from 6 specialized TB hospitals to 406 facilities, by analyzing National Health Laboratory Service data on TB during 2012-2015. We calculated the percentage of patients who visited a TB hospital <1 year after rifampin-resistant TB diagnosis, the median duration of their hospitalizations, and the total distance between facilities visited. We assessed temporal changes with linear regression and stratified results by location. Of 2,878 patients, 65% were from Cape Town. In Cape Town, 29% visited a TB hospital; elsewhere, 68% visited a TB hospital. We found that hospitalizations and travel distances were shorter in Cape Town than in the surrounding areas.

Mycobacterium bovis Infection in Free-Ranging African Elephants.

Mycobacterium bovis infection in wildlife species occurs worldwide. However, few cases of M. bovis infection in captive elephants have been reported. We describe 2 incidental cases of bovine tuberculosis in free-ranging African elephants (Loxodonta africana) from a tuberculosis-endemic national park in South Africa and the epidemiologic implications of these infections.