Study participants and samples
Subjects were recruited as part of the Gambian national TB prevalence study, a multistage cluster survey to estimate the burden of smear positive and bacteriologically confirmed PTB, conducted between December 2011 to January 2013. During the survey, one or more sputum samples were obtained from participants with respiratory symptoms and/or abnormal chest X-ray suspected of PTB. Smear microscopy, decontamination, liquid cultures and isolation of MTBC were performed using recommended standards and techniques as previously described15. The final radiological diagnosis was determined through consensus by a pulmonologist and radiologist.
From the parent survey, 903 suspected cases of NTM were identified from AFB positive cultures that tested negative for MTBC with the rapid identification assay BD MGIT TM TBc Identification (Becton, Dickinson and Company, USA). This assay tests for the MTBC-specific Mycobacterium Protein Target 64 (MPT64) antigen. All 903 were selected for inclusion in our study (Figure 1). Following retrieval of stored sputa, only samples for 575 (63.7%) of 903 participants yielded AFB positive cultures, all of which were selected for analysis.
The study population was categorized into four categories: 1 - no respiratory symptoms with abnormal chest X-ray (CXR); 2 - no respiratory symptoms with normal CXR; 3 - respiratory symptoms with abnormal CXR, and 4 - respiratory symptoms with normal CXR. Participants with NTM positive samples were assessed for clinically relevant disease by application of the ATS/IDSA diagnostic criteria11 and classified as confirmed, probable, suspected and non-PNTM disease.
Laboratory procedures for NTM identification
All laboratory procedures were performed at the TB Reference Laboratory, Medical Research Council Unit The Gambia at the London School of Hygiene & Tropical Medicine (MRCG at LSHTM). Sputa corresponding to AFB positive cultures and MPT64 rapid test negative in the original prevalence survey, were retrieved. Briefly, 0.5 mL of stored, decontaminated sputum samples from each suspected NTM case were cultured in Mycobacteria Growth Indicator Tubes (MGIT) supplemented with 0.8mL of BACTECTM MGITTM growth supplement and BBLTM MGITTM PANTATM antibiotic16. MGIT tubes were incubated at 37oC in the automated BACTEC MGIT 960TM system (Becton Dickinson Diagnostic Instrument Systems, Franklin Lakes, United States of America) until flagged as positive. Samples that failed to show any growth after 42 days of incubation in the machine were removed and classified as negative based on the manufacturer’s protocol. All cultures confirmed as acid-fast bacilli (AFB) with Ziehl-Neelsen (ZN) staining were classified as suspected NTM while non-AFB cultures were considered contaminants and excluded from the study.
16S rRNA gene PCRs and sequencing
Mycobacterial deoxyribonucleic acid (DNA) was extracted from MGIT cultures for 16S rRNA gene Polymerase Chain Reactions (PCR) by the boiled lysate method described by Aldous et al17. Briefly, 0.5 mL of culture was centrifuged at 10,000 rpm for 15 minutes. The resulting pellets were resuspended in 0.1mL sterile Tris ethylene diamine tetra acetic acid (TE) buffer® (Sigma Aldrich, St. Louis, Missouri, United States) and heated to 99oC in a heat block for 20 minutes followed by sonication for 15 minutes. Tubes were spun at 14,000 rpm for 5 minutes following which the supernatant was used for mycobacteria species-specific 16S rRNA PCRs and 2.5 μl of extracted DNA was added to 22.5 μl of master mix. We used primer pairs P1 (TGCTTAACACATGCAAGTCG) and P2 (TCTCTAGACGCGTCCTGTGC) to amplify regions A and B of the mycobacteria 16S rRNA gene. The PCR comprised a denaturation step of 5 minutes at 95°C, then 45 cycles of 95°C for 45 seconds, 56°C for 45 seconds, 72°C for 45 seconds and a final extension step at 72°C for 10 minutes. Amplified DNA products were viewed on a 1% (w/v) agarose gel stained with 500 ng/μl ethidium bromide.
16S rRNA gene sequencing was performed on purified PCR products by the Sanger chain-termination DNA sequencing method18. The obtained sequences were edited and analysed in BioEdit version 7.2.5 Ibis Biosciences http://www.mbio.ncsu.edu/BioEdit/bioedit.html. Background and ambiguous non-standard nucleotides were removed and the correct DNA nucleotides inserted. Consensus sequences were computed by matching forward and reverse traces of each sample using SeqTrace 0.9.0 software19 http://seqtrace.googlecode.com. We filtered low quality base calls and performed end trimming. High quality finished unknown sequences were exported in Fasta format for downstream analyses. NTM and MTBC 16S rRNA reference sequences derived from American Type Culture Collection (ATCC) isolates were downloaded from Gen Bank to illustrate the phylogenetic relationship of unknown sequences and identify them. These sequences were aligned together with unknown sequences in BioEdit software using ClustalW2 Multiple Sequence Alignment tool. A Phylogenetic tree was constructed from the aligned sequences based on regions A and B of the Mycobacteria gene using maximum likelihood algorithm. Unknown sequences were further identified based on their relatedness to NTM and MTBC reference strains on the phylogenetic tree. All sequences were made publicly available in GenBank with accession numbers KX607141.1 to KX607408.1.
Data management and statistical analysis
Data collection and management for the parent study have been described previously15. All data analyses for this study were carried out with Stata version 12 (Stata Corp, College Station, TX). The main outcome of the study was the prevalence of NTM amongst AFB positives cultures/MPT64 negative strains in pulmonary samples of a nationally representative population investigated for PTB in The Gambia. NTM prevalence and other categorical variables were summarised using frequency counts and proportions.
Following retrieval and re-culture, samples for 328 of the 903 originally AFB culture positives in the parent study failed to grow. To account for this disparity, NTM prevalence was estimated using two models: Model 1 in a complete case analysis used only results from the 575 participants with positive cultures after retrieval of stored sputum samples. Model 2 assumes the negative results for 328 of 903 participants, are truly negative.
Associations between NTM prevalence and categorical risk factors were examined using logistic regression models. Variables reaching statistical significance (p<0.05) in univariate analyses were included in multivariable logistic regression models to obtain odds ratios and their 95% confidence intervals. All models adjusted for clustering as per the parent study design.