Indra Bergval

Combined Species Identification, Genotyping, and Drug Resistance Detection of Mycobacterium tuberculosis Cultures by MLPA on a Bead-Based Array

The population structure of Mycobacterium tuberculosis is typically clonal therefore genotypic lineages can be unequivocally identified by characteristic markers such as mutations or genomic deletions. In addition, drug resistance is mainly mediated by mutations. These issues make multiplexed detection of selected mutations potentially a very powerful tool to characterise Mycobacterium tuberculosis. We used Multiplex Ligation-dependent Probe Amplification (MLPA) to screen for dispersed mutations, which can be successfully applied to Mycobacterium tuberculosis as was previously shown. Here we selected 47 discriminative and informative markers and designed MLPA probes accordingly to allow analysis with a liquid bead array and robust reader (Luminex MAGPIX technology). To validate the bead-based MLPA, we screened a panel of 88 selected strains, previously characterised by other methods with the developed multiplex assay using automated positive and negative calling. In total 3059 characteristics were screened and 3034 (99.2%) were consistent with previous molecular characterizations, of which 2056 (67.2%) were directly supported by other molecular methods, and 978 (32.0%) were consistent with but not directly supported by previous molecular characterizations. Results directly conflicting or inconsistent with previous methods, were obtained for 25 (0.8%) of the characteristics tested. Here we report the validation of the bead-based MLPA and demonstrate its potential to simultaneously identify a range of drug resistance markers, discriminate the species within the Mycobacterium tuberculosis complex, determine the genetic lineage and detect and identify the clinically most relevant non-tuberculous mycobacterial species. The detection of multiple genetic markers in clinically derived Mycobacterium tuberculosis strains with a multiplex assay could reduce the number of TB-dedicated screening methods needed for full characterization. Additionally, as a proportion of the markers screened are specific to certain Mycobacterium tuberculosis lineages each profile can be checked for internal consistency. Strain characterization can allow selection of appropriate treatment and thereby improve treatment outcome and patient management.

Development of Multiplex Assay for Rapid Characterization of Mycobacterium tuberculosis

ABSTRACT We have developed a multiplex assay, based on multiplex ligation-dependent probe amplification (MLPA), that allows simultaneous detection of multiple drug resistance mutations and genotype-specific mutations at any location in the Mycobacterium tuberculosis genome. The assay was validated on a reference panel of well-characterized strains, and the results show that M. tuberculosis can be accurately characterized by our assay. Eighteen discriminatory markers identifying drug resistance (rpoB, katG, inhA, embB), members of the M. tuberculosis complex (16S rRNA, IS6110, TbD1), the principal genotypic group (katG, gyrA), and Haarlem and Beijing strains (ogt, mutT2, mutT4) were targeted. A sequence specificity of 100% was reached for 16 of the 18 selected genetic targets. In addition, a panel of 47 clinical M. tuberculosis isolates was tested by MLPA in order to determine the correlation between phenotypic drug resistance and MLPA and between spoligotyping and MLPA. Again, all mutations present in these isolates that were targeted by the 16 functional probes were identified. Resistance-associated mutations were detected by MLPA in 71% of the identified rifampin-resistant strains and in 80% of the phenotypically isoniazid-resistant strains. Furthermore, there was a perfect correlation between MLPA results and spoligotypes. When MLPA is used on confirmed M. tuberculosis clinical specimens, it can be a useful and informative instrument to aid in the detection of drug resistance, especially in laboratories where drug susceptibility testing is not common practice and where the rates of multidrug-resistant and extensively drug resistant tuberculosis are high. The flexibility and specificity of MLPA, along with the ability to simultaneously genotype and detect drug resistance mutations, make MLPA a promising tool for pathogen characterization.

Resistant mutants of Mycobacterium tuberculosis selected in vitro do not reflect the in vivo mechanism of isoniazid resistance.

Objectives The high prevalence of isoniazid-resistant Mycobacterium tuberculosis is often explained by a high mutation rate for this trait, although detailed information to support this theory is absent. We studied the development of isoniazid resistance in vitro, making use of a laboratory strain of M. tuberculosis. Methods Spontaneous isoniazid-resistant mutants were characterized by molecular methods allowing identification of the most commonly encountered resistance-conferring mutations. Additionally, we determined the in vitro mutation rates for isoniazid and rifampicin resistance, and characterized the genome of a triple-resistant strain. Results Results confirm that the in vitro mutation rate for isoniazid resistance (3.2 × 10−7 mutations/cell division) is much higher than the rate for rifampicin resistance (9.8 × 10−9 mutations/cell division). However, in the majority of the in vitro mutants katG was partially or completely deleted and neither of the two most common in vivo mutations, katG-S315T or inhA-C(-)15T, were found in 120 isogenic mutants. This implies that clinically prevalent resistance mutations were present in <0.8% of isoniazid-resistant strains selected in vitro (95% CI 0%–2.5%). The triple-resistant strain had acquired isoniazid resistance via a 49 kbp deletion, which included katG. Apart from previously identified resistance-conferring mutations, three additional point mutations were acquired during sequential selection steps. Conclusions These outcomes demonstrate that the in vivo mechanism of isoniazid resistance is not reflected by in vitro experiments. We therefore conclude that the high in vitro mutation rate for isoniazid resistance is not a satisfactory explanation for the fact that isoniazid monoresistance is significantly more widespread than monoresistance to rifampicin.

Mycobacterium tuberculosis population structure determines the outcome of genetics-based second-line drug resistance testing

ABSTRACT The global emergence of multidrug-resistant tuberculosis has highlighted the need for the development of rapid tests to identify resistance to second-line antituberculosis drugs. Resistance to fluoroquinolones and aminoglycosides develops through nonsynonymous single nucleotide polymorphisms in the gyrA and gyrB genes and the rrs gene, respectively. Using DNA sequencing as the gold standard for the detection of mutations conferring resistance, in conjunction with spoligotyping, we demonstrated heteroresistance in 25% and 16.3% of Mycobacterium tuberculosis isolates resistant to ofloxacin and amikacin, respectively. Characterization of follow-up isolates from the same patients showed that the population structure of clones may change during treatment, suggesting different phases in the emergence of resistance. The presence of underlying mutant clones was identified in isolates which failed to show a correlation between phenotypic resistance and mutation in the gyrA or rrs gene. These clones harbored previously described mutations in either the gyrA or rrs gene, suggesting that rare mutations conferring resistance to ofloxacin or amikacin may not be as important as was previously thought. We concluded that the absence of a correlation between genotypic and phenotypic resistance implies an early phase in the emergence of resistance within the patient. Thus, the diagnostic utility of genetics-based drug susceptibility tests will depend on the proportion of patients whose bacilli are in the process of acquiring resistance in the study setting. These data have implications for the interpretation of molecular and microbiological diagnostic tests for patients with drug-susceptible and drug-resistant tuberculosis who fail to respond to treatment and for those with discordant results.

Pre-existing isoniazid resistance, but not the genotype of Mycobacterium tuberculosis drives rifampicin resistance codon preference in vitro.

Both the probability of a mutation occurring and the ability of the mutant to persist will influence the distribution of mutants that arise in a population. We studied the interaction of these factors for the in vitro selection of rifampicin (RIF)-resistant mutants of Mycobacterium tuberculosis. We characterised two series of spontaneous RIF-resistant in vitro mutants from isoniazid (INH)-sensitive and -resistant laboratory strains and clinical isolates, representing various M. tuberculosis genotypes. The first series were selected from multiple parallel 1 ml cultures and the second from single 10 ml cultures. RIF-resistant mutants were screened by Multiplex Ligation-dependent Probe Amplification (MLPA) or by sequencing the rpoB gene. For all strains the mutation rate for RIF resistance was determined with a fluctuation assay. The most striking observation was a shift towards rpoB-S531L (TCG→TTG) mutations in a panel of laboratory-generated INH-resistant mutants selected from the 10-ml cultures (p<0.001). All tested strains showed similar mutation rates (1.33×10−8 to 2.49×10−7) except one of the laboratory-generated INH mutants with a mutation rate measured at 5.71×10−7, more than 10 times higher than that of the INH susceptible parental strain (5.46–7.44×10−8). No significant, systematic difference in the spectrum of rpoB-mutations between strains of different genotypes was observed. The dramatic shift towards rpoB-S531L in our INH-resistant laboratory mutants suggests that the relative fitness of resistant mutants can dramatically impact the distribution of (subsequent) mutations that accumulate in a M. tuberculosis population, at least in vitro. We conclude that, against specific genetic backgrounds, certain resistance mutations are particularly likely to spread. Molecular screening for these (combinations of) mutations in clinical isolates could rapidly identify these particular pathogenic strains. We therefore recommend that isolates are screened for the distribution of resistance mutations, especially in regions that are highly endemic for (multi)drug resistant tuberculosis.

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates

BackgroundMultiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains.ResultsTo validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal.ConclusionBased on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Evaluation of SNP-based genotyping to monitor tuberculosis control in a high MDR-TB setting

Mycobacterium tuberculosis (Mtb) lineage identification and typing of clinical isolates in general is performed only retrospectively. The results are rarely linked to drug susceptibility testing (DST) or patient data. Consequently, the association between Mtb lineage, (multi)drug resistance and treatment history is not fully explored at the local level. Here we evaluated a new SNP based typing assay. We furthermore assessed the added value of genotyping of Mtb isolates for epidemiological purposes and guidance of tuberculosis (TB) control. Mtb lineage, DST profile and treatment history were determined for 399 samples at the National TB Reference Laboratory (NRL) in Tbilisi, Georgia by local staff. Data was shared electronically and analysis was performed remotely. Out of 399 isolates, 74 (74/399, 18.5%) were at least multidrug resistant (MDR)-TB, of which 63 (63/74, 85.1%) were members of three different Mtb Beijing lineages. Previous treatment was reported in 38/74 (51.4%) MDR(+) patients. The availability of this data allows associations with lineages. Notably, multidrug resistant TB was more strongly associated with the Beijing lineage than treatment history. Of all MDR-TB Beijing strains 56.7% (42/74) were members of a genetic cluster. This is most easily explained by (ongoing) MDR-TB transmission rather than drug resistance amplification. This knowledge is useful when designing intervention strategies for MDR-TB. Our study provides an example that on-site integrated Mtb genotyping is realistic and could support TB control activities.

Beijing Lineage of MDR Mycobacterium tuberculosis in Bulgaria, 2007–2011

To assess the spread of the Mycobacterium tuberculosis Beijing genotype among patients with multidrug-resistant and extensively resistant tuberculosis in Bulgaria, we genotyped 188 (72%) of 261 microbiologically confirmed resistant isolates obtained during 2007–2011. The estimated prevalence of the Beijing genotype among these patients was 3.2%.

Methylation in Mycobacterium tuberculosis is lineage specific with associated mutations present globally.

DNA methylation is an epigenetic modification of the genome involved in regulating crucial cellular processes, including transcription and chromosome stability. Advances in PacBio sequencing technologies can be used to robustly reveal methylation sites. The methylome of the Mycobacterium tuberculosis complex is poorly understood but may be involved in virulence, hypoxic survival and the emergence of drug resistance. In the most extensive study to date, we characterise the methylome across the 4 major lineages of M. tuberculosis and 2 lineages of M. africanum, the leading causes of tuberculosis disease in humans. We reveal lineage-specific methylated motifs and strain-specific mutations that are abundant globally and likely to explain loss of function in the respective methyltransferases. Our work provides a set of sixteen new complete reference genomes for the Mycobacterium tuberculosis complex, including complete lineage 5 genomes. Insights into lineage-specific methylomes will further elucidate underlying biological mechanisms and other important phenotypes of the epi-genome.

Pyrazinamide resistance-conferring mutations in pncA and the transmission of multidrug resistant TB in Georgia

BackgroundThe ongoing epidemic of multidrug-resistant tuberculosis (MDR-TB) in Georgia highlights the need for more effective control strategies. A new regimen to treat MDR-TB that includes pyrazinamide (PZA) is currently being evaluated and PZA resistance status will largely influence the success of current and future treatment strategies. PZA susceptibility testing was not routinely performed at the National Reference Laboratory (NRL) in Tbilisi between 2010 and September 2015. We here provide a first insight into the prevalence of PZA resistant TB in this region.MethodsPhenotypic susceptibility to PZA was determined in a convenience collection of well-characterised TB patient isolates collected at the NRL in Tbilisi between 2012 and 2013. In addition, the pncA gene was sequenced and whole genome sequencing was performed on two isolates.ResultsOut of 57 isolates tested 33 (57.9%) showed phenotypic drug resistance to PZA and had a single pncA mutation. All of these 33 isolates were MDR-TB strains. pncA mutations were absent in all but one of the 24 PZA susceptible isolate. In total we found 18 polymorphisms in the pncA gene. From the two major MDR-TB clusters represented (94–32 and 100–32), 10 of 15, 67.0% and 13 of 14, 93.0% strains, respectively were PZA resistant. We also identified a member of the potentially highly transmissive clade A strain carrying the characteristic I6L substitution in PncA. Another strain with the same MLVA type as the clade A strain acquired a different mutation in pncA and was genetically more distantly related suggesting that different branches of this particular lineage have been introduced into this region.ConclusionIn this high MDR-TB setting more than half of the tested MDR-TB isolates were resistant to PZA. As PZA is part of current and planned MDR-TB treatment regimens this is alarming and deserves the attention of health authorities. Based on our typing and sequence analysis results we conclude that PZA resistance is the result of primary transmission as well as acquisition within the patient and recommend prospective genotyping and PZA resistance testing in high MDR-TB settings. This is of utmost importance in order to preserve bacterial susceptibility to PZA to help protect (new) second line drugs in PZA containing regimens.