Whole genome sequencing for drug resistance profile prediction in Mycobacterium tuberculosis

Whole genome sequencing allows rapid detection of drug-resistant isolates. However, the availability of high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data has thus far been limited.We determined drug resistance profiles of 176 genetically diverse clinical isolates from Democratic Republic of the Congo, Ivory Coast, Peru, Thailand and Switzerland by quantitative phenotypic DST for 11 antituberculous drugs using the BD BACTEC MGIT 960 system and 7H10 agar dilution to generate a cross-validated phenotypic DST readout. We compared DST results with predicted drug resistance profiles inferred by whole genome sequencing.Classification of strains by the two phenotypic DST methods into resistotype/wild type populations was concordant in 73-99 % of cases, depending on the drug. Our data suggests that the established critical concentration (5 mg/L) for ethambutol resistance (MGIT 960 system) is too high and may misclassify strains as susceptible, compared to 7H10 agar dilution. Increased minimal inhibitory concentrations were explained by mutations identified by whole genome sequencing. Using whole genome sequences, we were able to predict quantitative drug resistance levels for the majority of drug resistance mutations. Predicting quantitative levels of drug resistance by whole genome sequencing was partially limited due to incompletely understood drug resistance mechanisms. The overall sensitivity and specificity of whole genome-based DST were 86.8 % and 94.5 %, respectively.Despite some limitations, whole genome sequencing has the potential to infer resistance profiles without the need for time-consuming phenotypic methods.