ABSTRACT
Whole Genome Sequencing (WGS) is a promising tool in the global fight against tuberculosis (TB). The aim of this study was to evaluate the use of WGS in routine conditions for detection of drug resistance markers and transmission clusters in a multidrug-resistant TB hot-spot area in Peru. For this, 140 drug-resistant Mycobacterium tuberculosis strains from Lima and Callao were prospectively selected and processed through routine (GenoType MTBDRsl and BACTEC MGIT) and WGS workflows, simultaneously. Resistance was determined in accordance with the World Health Organization mutation catalogue. Agreements between WGS and BACTEC results were calculated for rifampicin, isoniazid, pyrazinamide, moxifloxacin, levofloxacin, amikacin and capreomycin. Transmission clusters were determined using different cut-off values of Single Nucleotide Polymorphism differences. 100% (140/140) of strains had valid WGS results for 13 anti-TB drugs. However, the availability of final, definitive phenotypic BACTEC MGIT results varied by drug with 10-17% of invalid results for the seven compared drugs. The median time to obtain results of WGS for the complete set of drugs was 11.5 days, compared to 28.6-52.6 days for the routine workflow. Overall categorical agreement by WGS and BACTEC MGIT for the compared drugs was 96.5%. Kappa index was good (0.65≤k≤1.00), except for moxifloxacin, but the sensitivity and specificity values were high for all cases. 97.9% (137/140) of strains were characterized with only one sublineage (134 belonging to "lineage 4" and 3 to "lineage 2"), and 2.1% (3/140) were mixed strains presenting two different sublineages. Clustering rates of 3.6% (5/140), 17.9% (25/140) and 22.1% (31/140) were obtained for 5, 10 and 12 SNP cut-off values, respectively. In conclusion, routine WGS has a high diagnostic accuracy to detect resistance against key current anti-TB drugs, allowing results to be obtained through a single analysis and helping to cut quickly the chain of transmission of drug-resistant TB in Peru.
Subject(s)
Antitubercular Agents , Mycobacterium tuberculosis , Tuberculosis, Multidrug-Resistant , Whole Genome Sequencing , Mycobacterium tuberculosis/genetics , Mycobacterium tuberculosis/drug effects , Mycobacterium tuberculosis/isolation & purification , Peru/epidemiology , Tuberculosis, Multidrug-Resistant/microbiology , Tuberculosis, Multidrug-Resistant/epidemiology , Tuberculosis, Multidrug-Resistant/drug therapy , Whole Genome Sequencing/methods , Humans , Antitubercular Agents/pharmacology , Antitubercular Agents/therapeutic use , Polymorphism, Single Nucleotide , Drug Resistance, Multiple, Bacterial/genetics , Microbial Sensitivity Tests , Genome, Bacterial , Male , FemaleABSTRACT
BACKGROUND: Tuberculosis (TB) is a communicable, preventable and curable disease caused by the bacterium Mycobacterium tuberculosis (MTB). Peru is amongst the 30 countries with the highest burden of multidrug-resistant tuberculosis (MDR-TB) worldwide. In the fight against drug-resistant tuberculosis, the UKMYC6 microdilution plate was developed and validated by the CRyPTIC project. The objective of the study was to evaluate the use of the broth microdilution (BMD) plate methodology for susceptibility testing of drug-resistant MTB strains in Peru. METHODS: MTB strains isolated between 2015 and 2018 in Peru were used. 496 nationally-representative strains determined as drug-resistant by the routine 7H10 Agar Proportion Method (APM) were included in the present study. The Minimum Inhibitory Concentration (MIC) of 13 antituberculosis drugs were determined for each strain using the UKMYC6 microdilution plates. Diagnostic agreement between APM and BMD plate methodology was determined for rifampicin, isoniazid, ethambutol, ethionamide, kanamycin and levofloxacin. Phenotypes were set using binary (or ternary) classification based on Epidemiological cut-off values (ECOFF/ECV) proposed by the CRyPTIC project. Whole Genome Sequencing (WGS) was performed on strains with discrepant results between both methods. RESULTS: MIC distributions were determined for 13 first- and second-line anti-TB drugs, including new (bedaquiline, delamanid) and repurposed (clofazimine, linezolid) agents. MIC results were available for 80% (397/496) of the strains at 14 days and the remainder at 21 days. The comparative analysis determined a good agreement (0.64 ≤ k ≤ 0.79) for the drugs rifampicin, ethambutol, ethionamide and kanamycin, and the best agreement (k > 0.8) for isoniazid and levofloxacin. Overall, 12% of MIC values were above the UKMYC6 plate dilution ranges, most notably for the drugs rifampicin and rifabutin. No strain presented MICs higher than the ECOFF/ECV values for the new or repurposed drugs. Discrepant analysis using genotypic susceptibility testing by WGS supported half of the results obtained by APM (52%, 93/179) and half of those obtained by BMD plate methodology (48%, 86/179). CONCLUSIONS: The BMD methodology using the UKMYC6 plate allows the complete susceptibility characterization, through the determination of MICs, of drug-resistant MTB strains in Peru. This methodology shows good diagnostic performances for rifampicin, isoniazid, ethambutol, ethionamide, kanamycin and levofloxacin. It also allows for the characterization of MICs for other drugs used in previous years against tuberculosis, as well as for new and repurposed drugs recently introduced worldwide.
Subject(s)
Mycobacterium tuberculosis , Tuberculosis, Lymph Node , Tuberculosis, Multidrug-Resistant , Antitubercular Agents/pharmacology , Ethambutol/pharmacology , Ethionamide , Humans , Isoniazid , Kanamycin , Levofloxacin , Microbial Sensitivity Tests , Peru , Rifampin/pharmacology , Tuberculosis, Multidrug-Resistant/microbiologyABSTRACT
Peru has the highest burden of multidrug-resistant tuberculosis in the Americas region. Since 1999, the annual number of extensively drug-resistant tuberculosis (XDR-TB) Peruvian cases has been increasing, becoming a public health challenge. The objective of this study was to perform genomic characterization of Mycobacterium tuberculosis strains obtained from Peruvian patients with XDR-TB diagnosed from 2011 to 2015 in Peru. Whole genome sequencing (WGS) was performed on 68 XDR-TB strains from different regions of Peru. 58 (85.3%) strains came from the most populated districts of Lima and Callao. Concerning the lineages, 62 (91.2%) strains belonged to the Euro-American Lineage, while the remaining 6 (8.8%) strains belonged to the East-Asian Lineage. Most strains (90%) had high-confidence resistance mutations according to pre-established WHO-confident grading system. Discordant results between microbiological and molecular methodologies were caused by mutations outside the hotspot regions analysed by commercial molecular assays (rpoB I491F and inhA S94A). Cluster analysis using a cut-off ≤ 10 SNPs revealed that only 23 (34%) strains evidenced recent transmission links. This study highlights the relevance and utility of WGS as a high-resolution approach to predict drug resistance, analyse transmission of strains between groups, and determine evolutionary patterns of circulating XDR-TB strains in the country.
Subject(s)
Drug Resistance, Multiple, Bacterial/genetics , Extensively Drug-Resistant Tuberculosis/microbiology , Genome, Bacterial/genetics , Mycobacterium tuberculosis/genetics , Adolescent , Adult , Antitubercular Agents/therapeutic use , Drug Resistance, Multiple, Bacterial/drug effects , Extensively Drug-Resistant Tuberculosis/drug therapy , Female , Humans , Male , Microbial Sensitivity Tests/methods , Middle Aged , Mutation/genetics , Mycobacterium tuberculosis/drug effects , Peru , Tuberculosis, Multidrug-Resistant/drug therapy , Tuberculosis, Multidrug-Resistant/microbiology , Whole Genome Sequencing/methods , Young AdultABSTRACT
OBJECTIVE.: To determine the genetic structure of drug-resistant strains of Mycobacterium tuberculosis that circulated throughout Peru during the years 2011-2015, by using haplotypes obtained from a line probe assay. MATERIALS AND METHODS.: A total of 6589 samples that were admitted to the Instituto Nacional de Salud for routine diagnosis using the GenoType® MTBDRplus v2 assay were analyzed during the study period. Resistant haplotypes were created by concatenating 21 polymorphic sites of the evaluated genes using the line probe assay; and the association analysis was carried out with phenotypes obtained by the 7H10 agar ratio method. RESULTS.: The most frequent mutations were: rpoB S531L (55.4%) and rpoB D516V (18.5%) for rifampicin resistance, and katG S315T (59.5%) and inhA c-15t (25.7%) for isoniazid resistance. We obtained 13 representative haplotypes (87.8% of analyzed samples), 6 corresponded to the multidrug-resistant genotype, 4 to the isoniazid mono-resistant genotype and 3 to the rifampicin mono-resistant genotype. Eighteen regions and the province of Callao showed high haplotype diversity; four showed moderate diversity and two showed low diversity. CONCLUSIONS.: Most regions showed high haplotype diversity; in addition, most drug-resistant strains of Mycobacterium tuberculosis were concentrated in the cities of Lima and Callao. Likewise, drug-resistant Mycobacterium tuberculosis strains circulating in Peru mainly contain the genetic markers with the highest prevalence worldwide, which are associated with resistance to rifampicin and isoniazid.
OBJETIVO.: Determinar la estructura genética de las cepas drogorresistentes de Mycobacterium tuberculosis que circularon en todo el Perú durante los años 2011-2015 a través de haplotipos obtenidos de un ensayo con sondas en línea. MATERIALES Y MÉTODOS.: Se analizaron 6589 muestras que ingresaron al Instituto Nacional de Salud para el diagnóstico rutinario mediante el ensayo GenoType® MTBDRplus v2, durante el periodo de estudio. Se crearon haplotipos resistentes mediante la concatenación de 21 sitios polimórficos de los genes evaluados por el ensayo con sondas en línea, y se realizó el análisis de asociación con fenotipos obtenidos por el método de proporciones agar 7H10. RESULTADOS.: Las mutaciones de mayores frecuencias fueron: rpoB S531L (55,4%) y rpoB D516V (18,5%) para la resistencia a rifampicina, y katG S315T (59,5%) e inhA c-15t (25,7%) para la resistencia a isoniacida. Se obtuvieron 13 haplotipos representativos (87,8% de muestras analizadas) de los cuales seis correspondieron al genotipo multidrogorresistente, cuatro al genotipo monorresistente a isoniacida y tres al genotipo monorresistente a rifampicina. Dieciocho departamentos, y la provincia del Callao, presentaron una alta diversidad haplotípica; cuatro presentaron moderada diversidad y dos presentaron baja diversidad. CONCLUSIONES.: Existe una alta diversidad haplotípica en la mayoría de los departamentos, además de una concentración de las cepas de Mycobacterium tuberculosis drogorresistentes en las ciudades de Lima y Callao. Asimismo, las cepas de Mycobacterium tuberculosis con perfil drogorresistente que circulan en el Perú contienen principalmente los marcadores genéticos de mayor prevalencia a nivel mundial asociados con la resistencia frente a rifampicina e isoniacida.