Fast and low-cost decentralized surveillance of transmission of tuberculosis based on strain-specific PCRs tailored from whole genome sequencing data: a pilot study

Multicenter Preparedness Exercise Enables Rapid Development of Cluster-Specific PCR-Based Screening Assays from Bacterial Genomic Data

PCR-based screening assays targeting strain-specific genetic markers allow the timely detection and specific differentiation of bacterial strains. Especially in situations where an infection cluster occurs, fast assay development is crucial for supporting targeted control measures. However, the turnaround times (TATs) for assay setup may be high due to insufficient knowledge about screening assay methods, workflows, and software tools. Here, two blind-coded and quality-controlled ring trials were performed in which five German laboratories established PCR-based screening assays from genomic data that specifically target selected bacterial clusters within two bacterial monospecies sample panels. While the first ring trial was conducted without a time limit to train the participants and assess assay feasibility, in the second ring trial, a challenging time limit of 2 weeks was set to force fast assay development as soon as genomic data were available. During both ring trials, we detected high interlaboratory variability regarding the screening assay methods and targets, the TATs for assay setup, and the number of screening assays. The participants designed between one and four assays per cluster that targeted cluster-specific unique genetic sequences, genes, or single nucleotide variants using conventional PCRs, high-resolution melting assays, or TaqMan PCRs. Assays were established within the 2-week time limit, with TATs ranging from 4 to 13 days. TaqMan probe delivery times strongly influenced TATs. In summary, we demonstrate that a specific exercise improved the preparedness to develop functional cluster-specific PCR-based screening assays from bacterial genomic data. Furthermore, the parallel development of several assays enhances assay availability.

Rapid Identification of Relevant Microbial Strains by Identifying Multiple Marker Single Nucleotide Polymorphisms via Amplicon Sequencing: Epidemic Monkeypox Virus as a Proof of Concept

Despite the proven value of applying genomic data for epidemiological purposes, commonly used high-throughput sequencing formats are not adapted to the response times required to intervene and finally control outbreaks. In this study, we propose a fast alternative to whole-genome sequencing (WGS) to track relevant microbiological strains: nanopore sequencing of multiple amplicons including strain marker single nucleotide polymorphisms (SNPs). As a proof a concept, we evaluated the performance of our approach to offer a rapid response to the most recent public health global alarm, the monkeypox virus (MPXV) global outbreak. Through a multisequence alignment, a list of 42 SNPs were extracted as signature makers for this outbreak. Twenty primer pairs were designed to amplify in a multiplex PCR the regions including 22 of these SNPs. Amplicon pools were sequenced in a MinION device, and SNPs were called in real time by an in-house bioinformatic pipeline. A total of 120 specimens (95 MPXV-PCR positive, Ct values from 14 to 39) were selected. In 67.37% of the positive subset, all 22 SNPs were called. After excluding low viral load specimens, in 92% of samples ≥11 outbreak SNPs were called. No false positives were observed in any of the 25 negative specimens. The total turnaround time required for this strategy was 5 hours, and the cost per sample was 14 euros. Nanopore sequencing of multiple amplicons harboring signature SNPs escapes the targeting limitations of strain-specific PCRs and offers a powerful alternative to systematic WGS, paving the way to real-time genomic epidemiology and making immediate intervention possible to finally optimize transmission control. IMPORTANCE Nanopore sequencing of multiple amplicons harboring signature single nucleotide polymorphisms (SNPs) escapes the targeting limitations of strain-specific PCRs and offers a powerful alternative to systematic whole-genome analysis, paving the way to real-time genomic epidemiology and making immediate intervention possible to finally optimize transmission control.

Single Nucleotide Polymorphism-Based Real-Time PCR Screening Assay for Rapid Tracking of Bacterial Infection Clusters To Complement Whole-Genome Sequencing Efforts during Outbreak Investigations

Infection clusters of multidrug-resistant bacteria increase mortality and entail expensive infection control measures. Whereas whole-genome sequencing (WGS) is the current gold standard to confirm infection clusters, PCR-based assays targeting cluster-specific signatures, such as single nucleotide polymorphisms (SNPs) derived from WGS data, are more suitable to initially screen for cluster isolates within large sample sizes. Here, we evaluated four software tools (SeqSphere+, RUCS, Gegenees, and Find Differential Primers) regarding their efficiency to find SNPs within WGS data sets that were specific for two bacterial monospecies infection clusters but were absent from a WGS reference data set comprising several hundred diverse genotypes of the same bacterial species. Cluster-specific SNPs were subsequently used to establish a probe-based real-time PCR screening assay for in vitro differentiation between cluster and noncluster isolates. SeqSphere+ and RUCS found 2 and 24 SNPs for clusters 1 and 14 and 24 SNPs for cluster 2, respectively. However, some signatures detected by RUCS were not cluster specific. Interestingly, all SNPs identified by SeqSphere+ were also detected by RUCS. In contrast, analyses with the remaining tools either resulted in no SNPs (with Find Differential Primers) or failed (Gegenees). Design of six cluster-specific real-time PCR assays enabled reliable cluster screening in vitro. Our evaluation revealed that SeqSphere+ and RUCS identified cluster-specific SNPs that could be used for large-scale screening in surveillance samples via real-time PCR, thereby complementing WGS efforts. This faster and simplified approach for the surveillance of bacterial clusters will improve infection control measures and will enhance protection of patients and physicians. IMPORTANCE Infection clusters of multidrug-resistant bacteria threaten medical facilities worldwide and cause immense health care costs. In recent years, whole-genome sequencing (WGS) has been increasingly applied to detect and to further control bacterial clusters. However, as WGS is still expensive and time-consuming, its exclusive application for screening and confirmation of bacterial infection clusters contributes to high costs and enhanced turnaround times, which many hospitals cannot afford. Therefore, there is need for alternative methods that can enable further surveillance of bacterial clusters that are initially detected by WGS in a faster and more cost-efficient way. Here, we established a system based on real-time PCR that enables rapid large-scale sample screening for bacterial cluster isolates within 7 days after the initial detection of an infection cluster, thereby complementing WGS efforts. This faster and simplified surveillance of bacterial clusters will improve infection control measures and will enhance protection of patients and physicians.

Integrative transnational analysis to dissect tuberculosis transmission events along the migratory route from Africa to Europe

BACKGROUND

Growing international migration has increased the complexity of tuberculosis transmission patterns. Italy's decision to close its borders in 2018 made of Spain the new European porte entrée for migration from the Horn of Africa (HA). In one of the first rescues of migrants from this region at the end of 2018, tuberculosis was diagnosed in eight subjects, mainly unaccompanied minors.

METHODS

Mycobacterium tuberculosis isolates from these recently arrived migrants were analysed by Mycobacterial Interspersed Repetitive-Unit/Variable-Number of Tandem Repeat (MIRU-VNTR) and subsequent whole genome sequencing (WGS) analysis. Data were compared with those from collections from other European countries receiving migrants from the HA and a strain-specific PCR was applied for a fast searching of common strains. Infections in a cellular model were performed to assess strain virulence.

RESULTS

MIRU-VNTR analysis allowed identifying an epidemiological cluster involving three of the eight cases from Somalia (0 single-nucleotide polymorphisms between isolates, HA cluster). Following detailed interviews revealed that two of these cases had shared the same migratory route in most of the trip and had spent a long time at a detention camp in Libya. To confirm potential en route transmission for the three cases, we searched the same strain in collections from other European countries receiving migrants from the HA. MIRU-VNTR, WGS and a strain-specific PCR for the HA strain were applied. The same strain was identified in 12 cases from Eritrea diagnosed soon after their arrival in 2018 to the Netherlands, Belgium and Italy. Intracellular replication rate of the strain did not reveal abnormal virulence.

CONCLUSIONS

Our study suggests a potential en route transmission of a pan-susceptible strain, which caused at least 15 tuberculosis cases in Somalian and Eritrean migrants diagnosed in four different European countries.

Overlapping prison/community tuberculosis outbreaks in Costa Rica revealed by alternative analysis of suboptimal material

Costa Rica has a low incidence of tuberculosis. Thus, identifying transmission hotspots is key to implement interventions. A tuberculosis outbreak was suspected in a prison in Costa Rica. Given the suboptimal quality of the samples received in our laboratory in Madrid, we applied alternative schemes for their analysis. In the first scheme, we bypassed the standard approach of applying systematic mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) and used a strain-specific polymerase chain reaction (PCR) that allowed identifying a cluster involving six cases (C1). The second scheme followed the canonical MIRU-VNTR path coupled with a whole-genomic amplification step, by which a second unsuspected overlapping cluster (C2), was detected in the same prison. These findings justified the implementation of a surveillance programme adapted to local resources based on a tailored multiplex allele-specific oligonucleotide (ASO)-PCR targeting C1 and C2. Presence of the C2 strain at a different prison was determined. ASO-PCR was applied extensively and alerted to the active circulation of one of the strains within and beyond prisons. Our study shows that alternative methodological strategies may provide useful data in settings with lack of resources for performing systematic standard molecular epidemiology programmes and/or with suboptimal material for analysis.

In-Depth Study of a Nosocomial Outbreak Caused by Extensively Drug-Resistant Pseudomonas aeruginosa Using Whole Genome Sequencing Coupled With a Polymerase Chain Reaction Targeting Strain-Specific Single Nucleotide Polymorphisms

In 2013-2014, an outbreak involving 14 patients infected by an extensively drug-resistant strain of Pseudomonas aeruginosa was detected in a hospital in Madrid, Spain. Our objective was to evaluate an alternative strategy for investigating the outbreak in depth by means of molecular and genomic approaches. Pulsed-field gel electrophoresis (PFGE) was applied as a first-line approach, followed by a more refined whole genome sequencing analysis. Single nucleotide polymorphisms identified by whole genome sequencing were used to design a specific polymerase chain reaction (PCR) for screening unsuspected cases infected by the outbreak strain. Whole genome sequencing alerted us to the existence of greater genetic diversity than was initially assumed, splitting the PFGE-associated outbreak isolates into 4 groups, 2 of which represented coincidental transmission unrelated to the outbreak. A multiplex allele-specific PCR targeting outbreak-specific single nucleotide polymorphisms was applied to 290 isolates, which allowed us to identify 25 additional cases related to the outbreak during 2011-2017. Whole genome sequencing coupled with an outbreak-strain-specific PCR enabled us to markedly redefine the initial picture of the outbreak by 1) ruling out initially suspected cases, 2) defining likely independent coincidental transmission events, 3) predating the starting point of the outbreak, 4) capturing new unsuspected cases, and 5) revealing that the outbreak was still active.

Whole genome sequencing-based analysis of tuberculosis (TB) in migrants: rapid tools for cross-border surveillance and to distinguish between recent transmission in the host country and new importations

Background

The analysis of transmission of tuberculosis (TB) is challenging in areas with a large migrant population. Standard genotyping may fail to differentiate transmission within the host country from new importations, which is key from an epidemiological perspective.

Aim

To propose a new strategy to simplify and optimise cross-border surveillance of tuberculosis and to distinguish between recent transmission in the host country and new importations

Methods

We selected 10 clusters, defined by 24-locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR), from a population in Spain rich in migrants from eastern Europe, north Africa and west Africa and reanalysed 66 isolates by whole-genome sequencing (WGS). A multiplex-allele-specific PCR was designed to target strain-specific marker single nucleotide polymorphisms (SNPs), identified from WGS data, to optimise the surveillance of the most complex cluster.

Results

In five of 10 clusters not all isolates showed the short genetic distances expected for recent transmission and revealed a higher number of SNPs, thus suggesting independent importations of prevalent strains in the country of origin. In the most complex cluster, rich in Moroccan cases, a multiplex allele-specific oligonucleotide-PCR (ASO-PCR) targeting the marker SNPs for the transmission subcluster enabled us to prospectively identify new secondary cases. The ASO-PCR-based strategy was transferred and applied in Morocco, demonstrating that the strain was prevalent in the country.

Conclusion

We provide a new model for optimising the analysis of cross-border surveillance of TB transmission in the scenario of global migration.

Trends of Two Epidemic Multidrug-Resistant Strains of Mycobacterium tuberculosis in Argentina Disclosed by Tailored Molecular Strategy

Two Mycobacterium tuberculosis strains-M (sublineage 4.1) and Ra (sublineage 4.3)-have long prevailed in Argentina among patients with multidrug-resistant tuberculosis (MDR-TB). Recently, budget constraints have hampered the surveillance of MDR-TB transmission. Based on whole-genome sequence analysis, we used M- and Ra-specific single nucleotide polymorphisms to tailor two multiplex allele-specific polymerase chain reactions (PCRs), which we applied to 252 stored isolates (95% of all newly diagnosed MDR-TB cases countrywide, 2015-2017). Compared with the latest data available (2007-2009), the M strain has receded (80/324 to 20/252, P < 0.0001), particularly among cross-border migrants (12/58 to 0/53, P = 0.0003) and HIV-infected people (30/97 to 7/74, P = 0.0007), but it still accounts for 4/12 new cases of extensively drug-resistant TB. Differently, the Ra strain remained stable in frequency (39/324 to 33/252) and contributed marginally to the extensive drug-resistance load (1/12). Our novel strategy disclosed recent trends of the two major MDR-TB strains, providing meaningful data to allocate control interventions more efficiently.

Screening of inmates transferred to Spain reveals a Peruvian prison as a reservoir of persistent Mycobacterium tuberculosis MDR strains and mixed infections

It is relevant to evaluate MDR-tuberculosis in prisons and its impact on the global epidemiology of this disease. However, systematic molecular epidemiology programs in prisons are lacking. A health-screening program performed on arrival for inmates transferred from Peruvian prisons to Spain led to the diagnosis of five MDR-TB cases from one of the biggest prisons in Latin America. They grouped into two MIRU-VNTR-clusters (Callao-1 and Callao-2), suggesting a reservoir of two prevalent MDR strains. A high-rate of overexposure was deduced because one of the five cases was coinfected by a pansusceptible strain. Callao-1 strain was also identified in 2018 in a community case in Spain who had been in the same Peruvian prison in 2002-5. A strain-specific-PCR tailored from WGS data was implemented in Peru, allowing the confirmation that these strains were currently responsible for the majority of the MDR cases in that prison, including a new mixed infection.

A Mycobacterium tuberculosis Beijing strain persists at high rates and extends its geographic boundaries 20 years after importation

Transmission of Beijing Mycobacterium tuberculosis can be investigated based on genotypic analysis of clinical isolates. A Beijing strain began to spread on Gran Canaria Island, Spain, at the end of the last century. In 1996, only 3 years after its importation to the island, its frequency had increased to 27.1% of all the isolates. The strain was tracked during the following years, and the most recent data obtained corresponded to 2007-8, when its presence continued to be alarming (21%). In the current study, we updated data on the distribution of this strain 20 years (2013–2014) after it was first detected on the island and extended the analysis for the first time to all the mycobacteriology laboratories covering the population of the Canary Island archipelago. Rapid updating was enabled by means of 2 different strain-specific PCRs: one targeting a peculiar feature of the strain, which was identified based on an IS6110 copy mapping in the Rv2180c gene, and a newly defined strain-specific single nucleotide polymorphism, which was identified by whole-genome sequencing. The results showed that the strain has remained highly prevalent (20.90% of all isolates), has spread throughout the neighbouring islands, and has also reached high representativeness in them (11–32%).

The Evolution of Genotyping Strategies To Detect, Analyze, and Control Transmission of Tuberculosis

The introduction of genotypic tools to analyze Mycobacterium tuberculosis isolates has transformed our knowledge of the transmission dynamics of this pathogen. We discuss the development of the laboratory methods that have been applied in recent years to study the epidemiology of M. tuberculosis. This review integrates two approaches: on the one hand, it considers how genotyping techniques have evolved over the years; and on the other, it looks at how the way we think these techniques should be applied has changed. We begin by examining the application of fingerprinting tools to suspected outbreaks only, before moving on to universal genotyping schemes, and finally we describe the latest real-time strategies used in molecular epidemiology. We also analyze refined approaches to obtaining epidemiological data from patients and to increasing the discriminatory power of genotyping by techniques based on genomic characterization. Finally, we review the development of integrative solutions to reconcile the speed of PCR-based methods with the high discriminatory power of whole-genome sequencing in easily implemented formats adapted to low-resource settings. Our analysis of future considerations highlights the need to bring together the three key elements of high-quality surveillance of transmission in tuberculosis, namely, speed, precision, and ease of implementation.

Molecular epidemiology of tuberculosis

The application of genotyping tools allowed us to discriminate between the Mycobacterium tuberculosis isolates obtained in the laboratory. The differentiation between single strains opened the door to molecular epidemiology studies, which had helped us to progress in our knowledge of how this pathogen is transmitted in the progressively more complex socio-epidemiological scenario. The genetic stability of this microorganism led to develop specific methodologies, which are thoroughly revised in this chapter. In addition to their application in epidemiology, we review, how they can offer a response to different diagnostic and clinical challenges. Finally, we focus on describing the novel genomic revolution we are experiencing in the analysis of tuberculosis, the methodology in which it is based and the novel possibilities it offers, including new routes of integrating both the molecular and genomic languages in innovative post-genomic proposals, better suited to our real-life context.

Urgent Implementation in a Hospital Setting of a Strategy To Rule Out Secondary Cases Caused by Imported Extensively Drug-Resistant Mycobacterium tuberculosis Strains at Diagnosis

Current migratory movements require new strategies for rapidly tracking the transmission of high-risk imported Mycobacterium tuberculosis strains. Whole-genome sequencing (WGS) enables us to identify single-nucleotide polymorphisms (SNPs) and therefore design PCRs to track specific relevant strains. However, fast implementation of these strategies in the hospital setting is difficult because professionals working in diagnostics, molecular epidemiology, and genomics are generally at separate institutions. In this study, we describe the urgent implementation of a system that integrates genomics and molecular tools in a genuine high-risk epidemiological alert involving 2 independent importations of extensively drug resistant (XDR) and pre-XDR Beijing M. tuberculosis strains from Russia into Spain. Both cases involved commercial sex workers with long-standing tuberculosis (TB). The system was based on strain-specific PCRs tailored from WGS data that were transferred to the local node that was managing the epidemiological alert. The optimized tests were available for prospective implementation in the local node 33 working days after receiving the primary cultures of the XDR strains and were applied to all 42 new incident cases. An interpretable result was obtained in each case (directly from sputum for 27 stain-positive cases) and corresponded to the amplification profiles for strains other than the targeted pre-XDR and XDR strains, which made it possible to prospectively rule out transmission of these high-risk strains at diagnosis.

[Tuberculosis and immigration]

The incidence of tuberculosis worldwide is declining. However, in Western countries this decline is slower due to the impact of immigration. Tuberculosis in the immigrant population is related to health status in the country of origin and with overcrowding and poverty conditions in the host country. Immigrants with tuberculosis are younger, have a higher prevalence of extrapulmonary forms, greater proportion of drug resistance and higher treatment default rates than those of natives. New molecular techniques not only reduce diagnostic delay time but also allow the rapid identification of resistances and improve knowledge of transmission patterns. It is necessary to implement measures to improve treatment compliance in this population group like facilitating access to health card, the use of fixed-dose combination drugs, the participation of cultural mediators and community health workers and gratuity of drugs.

Ultrafast Assessment of the Presence of a High-Risk Mycobacterium tuberculosis Strain in a Population

A persistent 8-year infection by a Beijing Mycobacterium tuberculosis strain from a previous outbreak after importation from West Africa obliged us to investigate secondary cases. We developed a multiplex PCR method based on whole-genome sequencing to target strain-specific single nucleotide polymorphisms (SNPs). In 1 week, we analyzed 868 isolates stored over 6 years. Only 2 cases (immigrants from Guinea Conakry) harbored the strain, which ruled out transmission-despite opportunities-and challenged some of the advantages associated with Beijing strains.

Data for action: collection and use of local data to end tuberculosis

Accelerating progress in the fight against tuberculosis will require a drastic shift from a strategy focused on control to one focused on elimination. Successful disease elimination campaigns are characterised by locally tailored responses that are informed by appropriate data. To develop such a response to tuberculosis, we suggest a three-step process that includes improved collection and use of existing programmatic data, collection of additional data (eg, geographic information, drug resistance, and risk factors) to inform tailored responses, and targeted collection of novel data (eg, sequencing data, targeted surveys, and contact investigations) to improve understanding of tuberculosis transmission dynamics. Development of a locally targeted response for tuberculosis will require substantial investment to reconfigure existing systems, coupled with additional empirical data to evaluate the effectiveness of specific approaches. Without adoption of an elimination strategy that uses local data to target hotspots of transmission, ambitious targets to end tuberculosis will almost certainly remain unmet.

Clinical value of whole-genome sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) is now common as a result of new technologies that can rapidly sequence a complete bacterial genome for US$500 or less. Many studies have addressed questions about tuberculosis with WGS, and knowing the sequence of the entire genome, rather than only a few fragments, has greatly increased the precision of molecular epidemiology and contact tracing. Additionally, topics such as the mutation rate, drug resistance, the target of new drugs, and the phylogeny and evolution of the Mycobacterium tuberculosis complex bacteria have been elucidated by WGS. Nonetheless, WGS has not explained differences in transmissibility between strains, or why some strains are more virulent than others or more prone to development of multidrug resistance. With advances in technology, WGS of clinical specimens could become routine in high-income countries; however, its relevance will probably depend on easy to use software to efficiently process the sequences produced and accessible genomic databases that can be mined in future studies.