Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Baseline Predictors of Sputum Culture Conversion in Pulmonary Tuberculosis: Importance of Cavities, Smoking, Time to Detection and W-Beijing Genotype

  • Marianne E. Visser ,

    marianne.visser@gmail.com

    Affiliation School of Public Health, University of the Western Cape, Cape Town, South Africa

  • Michael C. Stead,

    Affiliations MRC Centre for Molecular and Cellular Biology, University of Stellenbosch, Cape Town, South Africa, Division of Medical Microbiology, Institute of Infectious Diseases and Molecular Medicine, University of Cape Town, Cape Town, South Africa

  • Gerhard Walzl,

    Affiliation MRC Centre for Molecular and Cellular Biology, University of Stellenbosch, Cape Town, South Africa

  • Rob Warren,

    Affiliation MRC Centre for Molecular and Cellular Biology, University of Stellenbosch, Cape Town, South Africa

  • Michael Schomaker,

    Affiliation School of Public Health and Family Medicine, Centre for Infectious Disease Epidemiology and Research, University of Cape Town, Cape Town, South Africa

  • Harleen M. S. Grewal,

    Affiliation Section of Microbiology and Immunology, The Gade Institute, University of Bergen, Haukeland University Hospital, Bergen, Norway

  • Elizabeth C. Swart,

    Affiliation Division of Dietetics, University of the Western Cape, Cape Town, South Africa

  • Gary Maartens

    Affiliation Division of Clinical Pharmacology, Department of Medicine, University of Cape Town, Cape Town, South Africa

Abstract

Background

Time to detection (TTD) on automated liquid mycobacterial cultures is an emerging biomarker of tuberculosis outcomes. The M. tuberculosis W-Beijing genotype is spreading globally, indicating a selective advantage. There is a paucity of data on the association between baseline TTD and W-Beijing genotype and tuberculosis outcomes.

Aim

To assess baseline predictors of failure of sputum culture conversion, within the first 2 months of antitubercular therapy, in participants with pulmonary tuberculosis.

Design

Between May 2005 and August 2008 we conducted a prospective cohort study of time to sputum culture conversion in ambulatory participants with first episodes of smear and culture positive pulmonary tuberculosis attending two primary care clinics in Cape Town, South Africa. Rifampicin resistance (diagnosed on phenotypic susceptibility testing) was an exclusion criterion. Sputum was collected weekly for 8 weeks for mycobacterial culture on liquid media (BACTEC MGIT 960). Due to missing data, multiple imputation was performed. Time to sputum culture conversion was analysed using a Cox-proportional hazards model. Bayesian model averaging determined the posterior effect probability for each variable.

Results

113 participants were enrolled (30.1% female, 10.5% HIV-infected, 44.2% W-Beijing genotype, and 89% cavities). On Kaplan Meier analysis 50.4% of participants underwent sputum culture conversion by 8 weeks. The following baseline factors were associated with slower sputum culture conversion: TTD (adjusted hazard ratio (aHR) = 1.11, 95% CI 1.02; 1.2), lung cavities (aHR = 0.13, 95% CI 0.02; 0.95), ever smoking (aHR = 0.32, 95% CI 0.1; 1.02) and the W-Beijing genotype (aHR = 0.51, 95% CI 0.25; 1.07). On Bayesian model averaging, posterior probability effects were strong for TTD, lung cavitation and smoking and moderate for W-Beijing genotype.

Conclusion

We found that baseline TTD, smoking, cavities and W-Beijing genotype were associated with delayed 2 month sputum culture. Larger studies are needed to confirm the relationship between the W-Beijing genotype and sputum culture conversion.

Introduction

Tuberculosis cure after anti-tubercular therapy is best measured by bacteriological relapse within 2 years after completion of treatment. Because relapse is uncommon after short course therapy large sample sizes are required to detect differences in relapse rate, which, together with extended follow up periods, makes clinical trials using relapse as the primary endpoint very expensive. Sputum culture conversion after 2 months of treatment is a recognized surrogate biomarker of cure [1], [2]. Time to detection (TTD) on automated liquid mycobacterial cultures reflects mycobacterial load more accurately than sputum smear or culture grading, and has been shown to predict month 2 culture conversion, tuberculosis recurrence and relapse [3], [4].

Important factors associated with smear and culture non-conversion after 2 months include increasing age, chest radiographic features (cavitation and extent of disease), and higher sputum smear and culture grading at diagnosis [5][10].

There is limited evidence on the influence of M. tuberculosis genotype on the clinical outcome of tuberculosis. It has been suggested that the W-Beijing genotype has a selective advantage over other M. tuberculosis strains, since it is emerging worldwide and is linked to multidrug resistance in many areas [11], [12]. Data is conflicting regarding the role of the W-Beijing genotype on the risk of treatment failure or relapse [13][16].

We assessed baseline predictors, including W-Beijing genotype, of failure of sputum culture conversion within the first 2 months in participants with first episodes of smear positive pulmonary tuberculosis attending a primary care centre in Cape Town, South Africa.

Methods

Ethics statement

The trial protocol was approved by the University of Cape Town Ethics and Research Committee and written consent was obtained from all participants.

Study population and setting

Participants with first episodes of sputum smear positive pulmonary tuberculosis were enrolled in a randomized controlled trial of a micronutrient intervention (vitamin A and zinc), which has been reported elsewhere [17]. In brief, adults who attended two primary care tuberculosis clinics in Delft, Cape Town, South Africa between May 2005 and August 2008 were recruited. Rates of poverty and unemployment are high in Delft and the tuberculosis incidence rate increased from 496 to 737 cases per 100 000 during the study period. Study participants received standard anti-tubercular therapy for 5 days a week consisting of combination tablets contributing 600 mg rifampicin, 300 mg isoniazid, 1.6 g pyrazinamide and 1.1 g ethambutol for participants weighing 38 to 55 kg (Rifafour, Aventis Pharma Pty Ltd, Johannesburg). Doses were adjusted for participants weighing less than 38 kg or more than 55 kg. All participants were on directly observed therapy at the clinic or in the community. Research staff collected early morning sputum specimens weekly from participants for 8 weeks. The micronutrient intervention had no effect on . time to sputum smear or culture conversion up to 8 weeks on Kaplan Meier analysis (P = 0.15 and P = 0.38 respectively; log rank test) The current analysis is restricted to trial participants with a positive baseline sputum culture and with data on susceptibility to rifampicin and isoniazid. Participants with rifampicin resistance during the 8 week follow-up period, were excluded.

Study procedures and selection of isolates

In keeping with national guidelines at the time, none of our HIV-infected participants received anti-retroviral therapy, but all received co-trimoxazole prophylaxis. All participants were screened for misuse of alcohol with the Cut down, Annoyance, Guilt and Eye-opener (CAGE) questionnaire [18]. Smoking status was categorized as either never or ever smoked. Body weight and height of participants was determined. The extent and size of lung cavities at baseline was assessed independently by two pulmonologists experienced in the use of the Chest Radiograph Reading and Recording System (CRRS)- and disagreements on radiographic reading were resolved by consensus [19].

Prior to anti-tubercular therapy, trial participants supplied one spontaneous sputum specimen and thereafter one early morning unassisted sputum specimen every week, up to 8 weeks. Sputum specimens were examined by means of fluorescent microscopy (Auramine stain) and graded according to international standards [20]. All specimens were decontaminated and cultured on liquid media using the BACTEC MGIT 960 system (Becton Dickinson, Sparts, Maryland). Positive cultures were confirmed as M. tuberculosis complex using an in-house PCR assay [21]. The date of culture conversion was taken as the date of the first negative culture provided there were no subsequent positive cultures. Participants whose first negative culture occurred at week 8, were regarded as converters. Routine phenotypic drug susceptibility testing for isoniazid and rifampicin of isolates was carried out using the MGIT 960 system from January 2008; isolates cultured prior to that date were tested retrospectively, but some isolates had lost viability.

During 2009, spoligotyping was performed on the majority of serial isolates collected during the trial period [22]. Strain-specific PCR assays were performed to detect M. tuberculosis infections by W-Beijing and non-W-Beijing isolates [23]. In order to optimize the baseline data, we captured the strain type for all the available baseline isolates, as well as those within the first 2 weeks of treatment, in order to categorize participants as either W-Beijing or non-W-Beijing strain types. Ten participants who had two strain types, one of which was W-Beijing, were categorized as having the W-Beijing strain for the analysis.

Statistical analysis

Data was captured in duplicate in Microsoft Office Excel 2003 and validated with SAS 9.2 software. To account for missing data, multiple imputation was conducted with the Amelia II software package [24],a reliable imputation procedure implemented in the statistical software package R,and five imputed data sets were generated. This number of data sets is generally considered adequate [25].

Cox proportional hazard regression modeling with the imputed datasets was used to estimate the hazard ratio of culture conversion for the following pre-determined baseline variables: age, gender, smoking status, alcohol abuse, sputum smear grading, sputum TTD, presence and extent of cavitation, M. tuberculosis genotype, HIV-infection, body mass index, haemoglobin, total lymphocyte count, serum C-Reactive Protein (CRP) (at baseline and change in concentrations within the first 2 weeks) and albumin. Hazard ratios are reported together with 95% confidence intervals based on Rubin's rules [25]. The proportional hazards assumption was verified by testing interaction effects of analysis time with baseline variables (α = 0.05) and graphically via log-log-plots [26].

To confirm the stability of our results we conducted a sensitivity analysis using Bayesian model averaging to each imputed data set (M = 20), and the results were combined by Rubins' rules [25], [27]. Model averaging is an alternative to model selection and combines estimates from different models, which contain different sets of variables, with the aim of achieving stable estimates that incorporate modelling uncertainty in addition to sampling uncertainty: Models that have a higher posterior probability to be correct receive a higher weight [26]. The combined weighted estimates can then be interpreted as a stabilized hazard ratio of culture conversion, which, together with adjusted standard errors, reflects the uncertainty in the modeling process. The posterior effect probabilities presented in our analysis state the probability that the hazard ratio in the Cox regression model for a variable is not one, thus providing a summary measure of how likely it is that a baseline variable has an effect and ranks the importance of the baseline variables.

Results

Baseline characteristics

A total of 154 participants with new smear-positive pulmonary tuberculosis were enrolled in the micronutrient trial, of whom 113 fulfilled the inclusion criteria for the current study (Fig. 1). Four of the 113 sputum isolates were mono-resistant to isoniazid. Complete data on all the variables investigated in the current analysis was available in 72 patients. The baseline characteristics of these participants before and after multiple imputation, is shown in Table 1. There were 11 HIV-infected participants (5 males, 6 females) and HIV-status was not determined in 9 who were lost to follow-up before the test could be done. All participants underwent chest radiograph examination at baseline, but 22 could not be located at the clinic site either due to mis-filing or loss (e.g. chest radiographs sent with participants who were referred to hospitals). Strain data was unavailable for seven participants whose sputum cultures were discarded before spoligotyping could be performed.

thumbnail
Figure 1. Participant selection from the parent micronutrient trial for the sputum culture conversion analysis.

https://doi.org/10.1371/journal.pone.0029588.g001

thumbnail
Table 1. Baseline characteristics of participants with drug sensitive smear-positive pulmonary tuberculosis.#

https://doi.org/10.1371/journal.pone.0029588.t001

Participants had a high baseline bacillary burden, as reflected by the large proportion with a high sputum smear grading and a short median duration of TTD of 7 days (Table 1). Lung cavities were common and there was a high prevalence of smoking and alcohol abuse among participants. Spoligotyping showed that almost half of our participants had the W-Beijing strain type.

By week 6, nine percent of participants were lost-to follow-up and by week 8, this increased to 15%. Two cultures from participants who had not undergone culture conversion up to week 7, had contaminated cultures at week 8 and were categorized as not having converted. On Kaplan Meier analysis 50.4% of participants had undergone culture conversion by 8 weeks (data not shown).

Predictors of delayed sputum culture conversion

Cox-proportional hazard regression analysis showed that TTD, the presence of lung cavities at diagnosis and ever smokers, were all associated with delayed culture conversion in univariate analysis. After adjustment for co-variates, TTD and lung cavities remained significant (Table 2). Although the W-Beijing genotype was not significant in univariate analysis, the confidence interval was suggestive of a delayed effect on conversion, even after adjustment for co-variates. INH mono-resistance was not associated with conversion in univariate analysis (HR = 0.59; 95%CI, 0.082; 4.28) and was not explored further due to the limited number of participants.

thumbnail
Table 2. Cox proportional hazards regression analysis of baseline variables associated with sputum culture conversion after multiple imputation*.

https://doi.org/10.1371/journal.pone.0029588.t002

The estimated hazard ratios of the variables remained similar after model averaging (estimates not shown) confirming the stability of our main results. After Bayesian model averaging ever smokers, the presence of lung cavities and TTD had a posterior effect probability of greater than 80%, thus having a rather strong effect on sputum conversion, even after taking model uncertainty into account. A moderate effect of W-Beijing genotype was also present (Table 2).

Discussion

We found that a shorter baseline TTD, tobacco smoking, the presence of lung cavities and W-Beijing genotype were associated with delayed sputum culture conversion within the first 2 months in adults with smear-positive pulmonary tuberculosis. Bayesian model averaging showed that TTD, lung cavitation and smoking had strong effects and the W-Beijing genotype a moderate effect. The association we found between W-Beijing and sputum culture conversion has only been reported previously among a small number of participants, who participated in a multicentre trial [28]. Our finding that a shorter TTD is associated with delayed culture conversion adds to the very limited data on this potential surrogate marker [3], [4]. An increase of 3 days in TTD at baseline was associated with a 40% increase in the likelihood of conversion. Cavitation is a well- known factor associated with delayed sputum culture conversion [3], [5] and relapse [29].

Our data suggests that W-Beijing strains may be associated with delayed sputum culture conversion within the first 2 months, after adjusting for baseline bacillary load, the presence of lung cavities and the effect of smoking. Since the association between W-Beijing and culture conversion was largely unaffected by the adjustment for bacillary load at diagnosis and cavitatory disease, it suggests that other factors such as the host response to anti-tubercular treatment, may be important. Almost half of our participants had the W-Beijing genotype at baseline or within the first 2 weeks of treatment. The proportion of W-Beijing strains reported among adults and children with tuberculosis in Cape Town has ranged from 17 to 33% [30][32] and data from a paediatric referral hospital shows a marked increase in the prevalence of W-Beijing strains among children with tuberculosis between 2000 and 2003 [32]. One of the properties accounting for the rapid spread of the W-Beijing strains in our area may be operating by prolonging transmission in treated patients.

There is conflicting data on the role of smoking in sputum culture conversion [33]. Our unadjusted analyses showed that ever smokers had a significant longer time to culture conversion, compared with never-smokers during the first 2 months of treatment (HR = 0.45 (95% CI: 0.25–0.82)). This finding is similar to that of two other clinical trials [3], [34]. Retrospective data show that smokers are more likely to have clinical symptoms of tuberculosis such as cough and dyspnoea, as well as cavitary lesions, compared with non-smokers [35], [36]. When we adjusted for the presence of cavities and other confounders, the statistical significance of ever smoking was somewhat reduced, but smoking still had a strong effect in the final model, confirmed by model averaging estimates.

We were not able to demonstrate any significant effect of either baseline CRP concentrations, or change in concentrations after 2 weeks on treatment, on the risk of sputum culture conversion. CRP is a non-specific acute phase reactant and a marker of macrophage activation in active tuberculosis that is increased at diagnosis in proportion to the extent of disease and decreases during and after treatment [37][39]. Pre-treatment serum CRP concentrations have a high sensitivity for the diagnosis of active tuberculosis among smear-negative tuberculosis suspects [40] and was recently shown to be a risk factor for sputum smear conversion among participants in a micronutrient trial in Nigeria [9].

Our study has several limitations. First, there were missing observations in our data. We dealt with this by imputing missing data, which has been shown to be superior to complete case analysis (in which only subjects with no missing values are analyzed). If data is missing at random, it means that the missingness of data depends on observed co-variables or the outcome [41]. We find this assumption reasonable in our study, since the missing data related mainly to data capturing issues. Therefore, multiple imputation yields unbiased estimates that correctly reflect the uncertainty related to the missing data.

Second, a large proportion of our study participants did not convert during the 2 month period (49.6%), resulting in less statistical power. Similar high rates of culture positivity have been observed at 2 months at African sites using liquid culture media [4], [9], [42]. One of the strengths of our analysis is that we measured the time to sputum culture conversion, instead of a binary endpoint of culture conversion. Third, we did not assess the presence of a BCG scar among our participants. The absence of a BCG scar was recently shown to be an important factor in culture non-conversion among a Tanzanian cohort [8]. On the contrary, it has been suggested that the widespread use of the BCG vaccine may aid in the spread of W-Beijing strains, since it offers poor protection against W-Beijing strains in experimental models [11]. The BCG coverage among our local infants is high (99%) [43]. Fourth, rates of smoking and alcohol misuse were high in our study, but similar to those reported in another study from our region [4]. Variables associated with sputum conversion may be different in populations with lower rates of smoking and alcohol abuse.

Delayed sputum culture conversion increases the risk of relapse and prolongs the period of infectiousness. Therefore our findings about factors delaying culture conversion have several implications for public health. Smoking is associated with many health hazards, including increasing the risk of developing active tuberculosis [33]. Global tobacco use is increasing [44] and remains high in most countries with a high burden of tuberculosis [45], [46]. Effective tobacco control measures are needed, particularly in low-income countries, where smoking rates are highest among poorer and less educated communities [47]. Earlier detection of tuberculosis should reduce the number of patients presenting with cavitary disease and those with a higher sputum bacillary load and thus contribute to earlier sputum conversion. Finally, as mentioned above, BCG vaccination offers particularly poor coverage of W-Beijing strains. Newer tuberculosis vaccines will hopefully have greater activity against all strains, including W-Beijing.

In conclusion, a shorter baseline TTD, tobacco smoking, the presence of lung cavities, and W-Beijing genotype were associated with delayed sputum culture conversion within 2 months in adults with smear-positive pulmonary tuberculosis. Larger studies in high-burden settings are needed to confirm the relationship between the W-Beijing genotype and sputum culture conversion.

Acknowledgments

We would like to thank our field staff as well as the clinic staff at the Delft and Delft-South Community Health Centres, for their valuable contributions towards the micronutrient trial. M.C.S. was responsible for strain typing under supervision of Dr. H. Segal at the Institute for Infectious Disease and Molecular Medicine, University of Cape Town and R.W. We would also like to thank the following people and institutions: Dr. A. Whitelaw, V. January and K. Mentoor from the NHLS Groote Schuur Microbiology Laboratory; M. Marais from the Nutritional Intervention Unit, MRC for the biochemical analyses; A. Menezes, Dr. N. Chegou and M.de Kock from the MRC Centre of Molecular and Cellular Biology at the University of Stellenbosch for the processing and storing of all blood samples; Dr. R. Dawson and Dr. K. Narunsky from the UCT Lung Institute for the evaluation of the chest radiographs.

Author Contributions

Conceived and designed the experiments: MEV HMSG GW RW ECS GM. Performed the experiments: MEV MCS. Analyzed the data: MEV MCS MS GM. Contributed reagents/materials/analysis tools: HMSG GW RW MS. Wrote the paper: MEV MCS GW RW HMSG MS ECS GM.

References

  1. 1. Mitchison DA (1993) Assessment of new sterilizing drugs for treating pulmonary tuberculosis by culture at 2 months. Am Rev Respir Dis 147: 1062–1063.
  2. 2. Wallis RS, Doherty TM, Onyebujoh P, Vahedi M, Laang H, et al. (2009) Biomarkers for tuberculosis disease activity, cure, and relapse. Lancet Infect Dis 9: 162–172.
  3. 3. Dorman SE, Johnson JL, Goldberg S, Muzanye G, Padayatchi N, et al. (2009) Substitution of moxifloxacin for isoniazid during intensive phase treatment of pulmonary tuberculosis. Am J Resp Crit Care Med 180: 273–280.
  4. 4. Hesseling AC, Walzl G, Enarson EA, Carroll NM, Duncan K, et al. (2010) Baseline sputum time to detection predicts month 2 culture conversion and recurrent tuberculosis in HIV-uninfected patients. Int J Tuberc Lung Dis 14(5): 560–570.
  5. 5. Telzak EE, Barkat AF, Pollard CL, Turett GS, Justman JE, et al. (1997) Factors influencing time to sputum conversion among patients with smear-positive pulmonary tuberculosis. Clin Infec Dis 25: 666–670.
  6. 6. Güler M, Ünsal E, Dursun B, Aydin Ö, Capan N (2007) Factors influencing sputum smear and culture conversion time among patients with new case pulmonary tuberculosis. Int J Clin Prac 61(2): 231–235.
  7. 7. Banu Rekha VV, Balasubramanian R, Swaminathan S, Ramachandran R, Rahman F, et al. (2007) Sputum conversion at the end of intensive phase of Category-1 regimen in the treatment of pulmonary tuberculosis patients with diabetes mellitus or HIV infection: An analysis of risk factors. Indian J Med Res 126: 452–458.
  8. 8. Wang J, Lee L, Yu C, Chien Y, Yang P, et al. (2009) Factors influencing time to smear conversion in patients with smear-positive pulmonary tuberculosis. Respirology 14: 1012–1019.
  9. 9. Jeremiah K, PrayGod G, Faurholt-Jepsen , Range N, Andersen AB, et al. (2010) BCG vaccination status may predict sputum conversion in patients with pulmonary tuberculosis: a new consideration for an old vaccine? Thorax 65(12): 1072–1076.
  10. 10. Lawson L, Thacher TD, Yassin MA, Onuoha NA, Usman A, et al. (2010) Randomized controlled trial of zinc and vitamin A as co-adjuvants for the treatment of pulmonary tuberculosis. Trop Med Int Health 15(12): 1481–1490.
  11. 11. Nicol MP, Wilkinson RJ (2008) The clinical consequences of strain diversity in Mycobacterium tuberculosis. Trans R Soc Trop Med Hyg 102(10): 955–965.
  12. 12. Parwati I, van Crevel R, van Soolingen D (2010) Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains. Lancet Infect Dis 10(2): 103–111.
  13. 13. Lan NTN, Lien HTK, Tung LB, Borgdorff MW, Kremer K, et al. (2006) Mycobacterium tuberculosis Beijing Genotype and risk for treatment failure and relapse, Vietnam. Emerg Infect Dis 12: 736–743.
  14. 14. Burman WJ, Bliven EE, Cowan L, Bozeman L, Nahid P, et al. (2009) Relapse associated with active disease caused by Beijing strain of Mycobacterium tuberculosis. Emerg Infect Dis 15(7): 1061–1067.
  15. 15. Parwati I, AlisjahbanaB , Apriani L, Soetikno RD, Ottenhoff TH, et al. (2010) Mycobacterium tuberculosis Beijing genotype is an independent risk factor for tuberculosis treatment failure in Indonesia. J Infect Dis 201: 553–557.
  16. 16. Buu TN, Huyen MNT, van Soolingen D, Lan NTN, Quy HT, et al. (2010) The Mycobacterium tuberculosis Beijing genotype does not affect tuberculosis treatment failure in Vietnam. Clin Infect Dis 51(8): 879–886.
  17. 17. Visser ME, Grewal HMS, Swart EC, Dhansay MA, Walzl G, et al. (2011) The effect of vitamin A and zinc supplementation on treatment outcomes in pulmonary tuberculosis: a randomised controlled trial. Am J Clin Nutr 93(1): 93–100.
  18. 18. Schoeman JH, Parry CDH, Lombard CJ, Klopper HJ (1994) Assessment of alcohol-screening instruments in tuberculosis patients. Tuber Lung Dis 75: 371–376.
  19. 19. Den Boon S, Bateman ED, Enarson DA, Borgdorff MW, Verver CJ, et al. (2005) Development and evaluation of a new chest radiograph reading and recording system for epidemiological surveys of tuberculosis and lung disease. Int J Tuberc Lung Dis 9(10): 1088–1096.
  20. 20. World Health Organization, Regional Office for South -East Asia (2006) Chapter 17 Tuberculosis.
  21. 21. De Wit D, Steyn L, Shoemaker S, Sogin M (1990) Direct detection of Mycobacterium tuberculosis in clinical specimens by DNA amplification. J Clin Microbiol 28: 2437–2441.
  22. 22. Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, et al. (1997) Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol 35(4): 907–914.
  23. 23. Mallard K, McNerney R, Crampin AC, Houben R, Ndlovu R, et al. (2010) Molecular Detection of Mixed Infections of Mycobacterium tuberculosis Strains in Sputum Samples from Patients in Karonga District, Malawi. J Clin Microbiol 48(12): 4512–4518.
  24. 24. Honaker J, King G, Blackwell M (2011) Amelia II: A program for missing data. R Package version 1.5–5. Available: http://www.gking.harvard.edu/amelia. Accessed 2011 Nov 29.
  25. 25. Little RJA, Rubin DB (2002) Statistical analysis with missing data. Second edition. New York: Wiley. 408 p.
  26. 26. Volinsky CT, Madigan D, Raftery AE, Kronmal RA (1997) Bayesian model averaging in proportional hazard models. Assessing the risk of a stroke. Applied Statistics- Journal of the Royal Statistical Society Series C 46(4): 433–448.
  27. 27. Schomaker M, Wan ATK, Heumann C (2010) Frequentist model averaging with missing observations. Computational Statistics & Data analysis 54: 3336–3347.
  28. 28. Nahid P, Bliven EE, Kim EY, Mac Kenzie WR, Stout JE, et al. (2010) Influence of M. tuberculosis lineage variability within a clinical trial for pulmonary tuberculosis. PLoS ONE 5(5): e10753.
  29. 29. Benator D, Bhattacharaya M, Bozeman L, Burman W, Cantazaro A, et al. (2002) Rifapentine and isoniazid once a week versus rifampicin and isoniazid twice a week for treatment of drug-susceptible pulmonary tuberculosis in HIV-negative patients: a randomized trial. Lancet 360(9332): 528–534.
  30. 30. Richardson M, van Lill SWP, van der Spuy GD, Munch Z, Booysen CN, et al. (2002) Historic and recent events contribute to the disease dynamics of Beijing-like Mycobacterium tuberculosis isolates in a high incidence region. Int J Tuberc Lung Dis 6(11): 1001–1011.
  31. 31. Marais BJ, Victor RC, Hesseling AC, Barnard M, Jordaan A, et al. (2006) Beijing and Haarlem genotypes are overrepresented among children with drug-resistant tuberculosis in the Western Cape Province of South Africa. J Clin Microbiol 44(10): 3539–3543.
  32. 32. Cowley D, Govender D, February B, Wolfe M, Steyn L, et al. (2008) Recent and rapid emergence of W-Beijing strains of Mycobacterium tuberculosis in Cape Town, South Africa. Clin Infect Dis 47: 1252–1259.
  33. 33. Slama K, Chian C-Y, Enarson DA, Hassmiller K, Fanning A, et al. (2007) Tobacco and tuberculosis: a qualitative systematic review and meta-analysis. Int J Tuberc Lung Dis 11(10): 1049–1061.
  34. 34. Durban Immunotherapy group (1999) Immunotherapy with Micobacterium vaccae in patients with newly diagnosed pulmonary tuberculosis: a randomized controlled trial. Lancet 354: 116–119.
  35. 35. Leung CC, Yew WW, Chan CK, Tam CM, Lam CW, et al. (2003) Smoking and tuberculosis in Hong Kong. Int J Tuberc Lung Dis 7(10): 980–986.
  36. 36. Altet-Gomez MN, Alcaide J, Godoy P, Romero MA, del Rey H (2005) Clinical and epidemiological aspects of smoking and tuberculosis: a study of 13,038 cases. Int J Tuberc Lung Dis 9(4): 430–436.
  37. 37. Scott GM, Murphy PG, Gemidjioglu ME (1990) Predicting deterioration of treated tuberculosis by corticosteroid reserve and C-reactive protein. J Infect 21(1): 61–69.
  38. 38. Plit ML, Theron AJ, Fickl H, van Rensburg CE, Pendel S, et al. (1998) Influence of antimicrobial chemotherapy and smoking status on the plasma concentrations of vitamin C, Vitamin E, beta-carotene, acute phase reactants, iron and lipid peroxides in patients with pulmonary tuberculosis. Int J Tuberc Lung Dis 2(7): 590–596.
  39. 39. Djoba Siawaya JF, Bapela NB, Ronacher K, Veenstra H, Kidd M, et al. (2008) Immune parameters as markers of tuberculosis extent of disease and early prediction of anti-tuberculosis chemotherapy response. J Infect 56: 340–347.
  40. 40. Wilson D, Badri M, Maartens G (2011) Performance of serum C-Reactive protein as a screening test for smear-negative tuberculosis in an ambulatory high HIV prevalence population. PLoS ONE 6(1): e15248.
  41. 41. Van der Heijden GJMG, Donders ART, Stijnen T, Moons KGM (2006) Imputation of missing values is superior to complete case analysis and the missing-indicator method in multivariable diagnostic research: A clinical example. J Clin Epidemiol 59: 1102–1109.
  42. 42. Rustomjee R, Lienhardt C, Kanyok T, Davies GR, Levin J, et al. (2008) A Phase II study of the sterilising activities of ofloxacin, gatifloxacin and moxifloxacin in pulmonary tuberculosis. Int J Tuberc Lung Dis 12(2): 128–138.
  43. 43. Corrigal J (2006) Vaccination coverage of the Western Cape Province, 2005. Provincial Government of the Western Cape, Cape Town, South Africa.
  44. 44. Food and Agriculture Organization of the United Nations (2003) Projections of tobacco production, consumption and trade to the year 2010. FAO, Rome, Italy.
  45. 45. World Health Organization (2011) WHO report on the global tobacco epidemic. World Health Organization, Geneva. Available: http://www.who.int/tobacco. Accessed 2011 Oct 25.
  46. 46. World Health Organization (2011) Global tuberculosis control 2011. World Health Organization, Geneva. Available: http://www.who.int/tb/publications/global_report/en/www.who. Accessed 2011 Oct 25.
  47. 47. Hosseinpoor AR, Parker LA, d'Espaignet , Chatterji S (2011) Social determinants of smoking in low- and middle-income countries: results from the World Health survey. PLoS ONE 6(5): e20331.