Application of machine learning techniques to tuberculosis drug resistance analysis.
Kouchaki S., Yang Y., Walker TM., Walker AS., Wilson DJ., Peto TEA., Crook DW., Clifton DA., CRyPTIC consortium None.
Motivation:Timely identification of Mycobacterium tuberculosis (MTB) resistance to existing drugs is vital to decrease mortality and prevent the amplification of existing antibiotic resistance. Machine learning methods have been widely applied for timely predicting resistance of MTB given a specific drug and identifying resistance markers. However, they have been not validated on a large cohort of MTB samples from multi-centers across the world in terms of resistance prediction and resistance marker identification. Summary:Several machine learning classifiers and linear dimension reduction techniques were developed and compared for a cohort of 13402 isolates collected from 16 countries across six continents and tested 11 drugs. Results:Compared to conventional molecular diagnostic test, area under curve (AUC) of the best machine learning classifier increased for all drugs especially by 23.11%, 15.22%, and 10.14% for pyrazinamide (PZA), ciprofloxacin (CIP), and ofloxacin (OFX) respectively (p < 0.01). Logistic regression (LR) and gradient tree boosting (GBT) found to perform better than other techniques. Moreover, LR/GBT with a sparse principal component analysis/non-negative matrix factorisation step compared with the classifier alone enhanced the best performance in terms of F1-score by 12.54%, 4.61%, 7.45%, and 9.58% for amikacin (AK), moxifloxacin (MOX), OFX, and capreomycin (CAP) respectively, as well increasing AUC for AK and CAP. Results provided a comprehensive comparison of various techniques and confirmed the application of machine learning for better prediction of the large diverse TB data. Furthermore, mutation ranking showed the possibility of finding new resistance/susceptible markers. Availability:The source code can be found at http://www.robots.ox.ac.uk/~davidc/code.php. Supplementary information:Supplementary data are available at Bioinformatics online.