Standard Genotyping Overestimates Transmission of Mycobacterium tuberculosis among Immigrants in a Low-Incidence Country

Transmission of drug-resistant Mycobacterium tuberculosis isolates between Finnish- and foreign-born cases, 2014-2021: A molecular epidemiological study

BACKGROUND

Data on the molecular epidemiology and transmission of drug-resistant Mycobacterium tuberculosis (MTB) in low-incidence settings with immigration from high-incidence settings is limited.

METHOD

We included 115 drug-resistant (DR) MTB isolates with whole-genome sequencing data isolated in Finland between 2014 and 2021. Potential transmission clusters were identified using a threshold of 12 single-nucleotide polymorphisms (SNPs). Highly related clusters were identified using a threshold of 5 SNPs.

RESULT

Of the 115 DR MTB isolates, 31 (27.0%) isolates were from Finnish-born cases and 84 (73.0%) were from foreign-born cases. The proportion of multidrug-resistant (MDR) MTB isolates (30/84, 35.7%) from foreign-born cases was higher than that of MDR MTB isolates from Finnish-born cases (8/31, 25.8%). Lineage 2 (40/115, 34.8%) and lineage 4 (40/115, 34.8%) were the most prevalent lineages. A total of 25 (21.7%) isolates were classified into eight potential transmission clusters (≤12 SNPs). Furthermore, five highly related clusters (≤5 SNPs) were identified, including three DR MTB isolates from Finnish-born cases and 14 DR isolates from foreign-born cases.

CONCLUSION

The risk of DR MTB transmission between Finnish- and foreign-born persons is not negligible. Further research on clustering analysis in drug-susceptible MTB is worth to inform tuberculosis management and control in low-incidence settings with increasing immigration.

Tuberculosis in Ukrainian War Refugees and Migrants in the Czech Republic and Slovakia: A Molecular Epidemiological Study

BACKGROUND

The war in Ukraine has led to significant migration to neighboring countries, raising public health concerns. Notable tuberculosis (TB) incidence rates in Ukraine emphasize the immediate requirement to prioritize approaches that interrupt the spread and prevent new infections.

METHODS

We conducted a prospective genomic surveillance study to assess migration's impact on TB epidemiology in the Czech Republic and Slovakia. Mycobacterium tuberculosis isolates from Ukrainian war refugees and migrants, collected from September 2021 to December 2022 were analyzed alongside 1574 isolates obtained from Ukraine, the Czech Republic, and Slovakia.

RESULTS

Our study revealed alarming results, with historically the highest number of Ukrainian tuberculosis patients detected in the host countries. The increasing number of cases of multidrug-resistant TB, significantly linked with Beijing lineage 2.2.1 (p < 0.0001), also presents substantial obstacles to control endeavors. The genomic analysis identified the three highly related genomic clusters, indicating the recent TB transmission among migrant populations. The largest clusters comprised war refugees diagnosed in the Czech Republic, TB patients from various regions of Ukraine, and incarcerated individuals diagnosed with pulmonary TB specialized facility in the Kharkiv region, Ukraine, pointing to a national transmission sequence that has persisted for over 14 years.

CONCLUSIONS

The data showed that most infections were likely the result of reactivation of latent disease or exposure to TB before migration rather than recent transmission occurring within the host country. However, close monitoring, appropriate treatment, careful surveillance, and social support are crucial in mitigating future risks, though there is currently no evidence of local transmission in EU countries.

SARS-CoV-2 lineage assignments using phylogenetic placement/UShER are superior to pangoLEARN machine-learning method

With the rapid spread and evolution of SARS-CoV-2, the ability to monitor its transmission and distinguish among viral lineages is critical for pandemic response efforts. The most commonly used software for the lineage assignment of newly isolated SARS-CoV-2 genomes is pangolin, which offers two methods of assignment, pangoLEARN and pUShER. PangoLEARN rapidly assigns lineages using a machine-learning algorithm, while pUShER performs a phylogenetic placement to identify the lineage corresponding to a newly sequenced genome. In a preliminary study, we observed that pangoLEARN (decision tree model), while substantially faster than pUShER, offered less consistency across different versions of pangolin v3. Here, we expand upon this analysis to include v3 and v4 of pangolin, which moved the default algorithm for lineage assignment from pangoLEARN in v3 to pUShER in v4, and perform a thorough analysis confirming that pUShER is not only more stable across versions but also more accurate. Our findings suggest that future lineage assignment algorithms for various pathogens should consider the value of phylogenetic placement.

Genomic transmission clusters and circulating lineages of Mycobacterium tuberculosis among refugees residing in refugee camps in Ethiopia

BACKGROUND

Understanding the transmission dynamics of Mycobacterium tuberculosis (Mtb) could benefit the design of tuberculosis (TB) prevention and control strategies for refugee populations. Whole Genome Sequencing (WGS) has not yet been used to document the Mtb transmission dynamics among refugees in Ethiopia. We applied WGS to accurately identify transmission clusters and Mtb lineages among TB cases in refugee camps in Ethiopia.

METHOD AND DESIGN

We conducted a cross-sectional study of 610 refugees in refugee camps in Ethiopia presenting with symptoms of TB. WGS data of 67 isolates was analyzed using the Maximum Accessible Genome for Mtb Analysis (MAGMA) pipeline; iTol and FigTree were used to visualize phylogenetic trees, lineages, and the presence of transmission clusters.

RESULTS

Mtb culture-positive refugees originated from South Sudan (52/67, 77.6%), Somalia (9/67, 13.4%). Eritrea (4/67, 6%), and Sudan (2/67, 3%). The majority (52, 77.6%) of the isolates belonged to Mtb lineage (L) 3, and one L9 was identified from a Somalian refugee. The vast majority (82%) of the isolates were pan-susceptible Mtb, and none were multi-drug-resistant (MDR)-TB. Based on the 5-single nucleotide polymorphisms cutoff, we identified eight potential transmission clusters containing 23.9% of the isolates. Contact investigation confirmed epidemiological links with either family or social interaction within the refugee camps or with neighboring refugee camps.

CONCLUSION

Four lineages (L1, L3, L4, and L9) were identified, with the majority of strains being L3, reflecting the Mtb L3 dominance in South Sudan, where the majority of refugees originated from. Recent transmission among refugees was relatively low (24%), likely due to the short study period. The improved understanding of the Mtb transmission dynamics using WGS in refugee camps could assist in designing effective TB control programs for refugees.

Whole genome sequencing improves the discrimination between Mycobacterium bovis strains on the southern border of Kruger National Park, South Africa

Background

Mycobacterium bovis forms part of the Mycobacterium tuberculosis complex and has an extensive host range and zoonotic potential. Various genotyping methods (e.g., spoligotyping) have been used to describe the molecular epidemiology of M. bovis. Advances in whole genome sequencing (WGS) have increased resolution to enable detection of genomic variants to the level of single nucleotide polymorphisms. This is especially relevant to One Health research on tuberculosis which benefits by being able to use WGS to identify epidemiologically linked cases, especially recent transmission. The use of WGS in molecular epidemiology has been extensively used in humans and cattle but is limited in wildlife. This approach appears to overcome the limitations of conventional genotyping methods due to lack of genetic diversity in M. bovis.

Methods

This pilot study investigated the spoligotype and WGS of M. bovis isolates (n = 7) from wildlife in Marloth Park (MP) and compared these with WGS data from other South African M. bovis isolates. In addition, the greater resolution of WGS was used to explore the phylogenetic relatedness of M. bovis isolates in neighbouring wildlife populations.

Results

The phylogenetic analyses showed the closest relatives to the seven isolates from MP were isolates from wildlife in Kruger National Park (KNP), which shares a border with MP. However, WGS data indicated that the KNP and MP isolates formed two distinct clades, even though they had similar spoligotypes and identical in silico genetic regions of difference profiles.

Conclusions

Mycobacterium bovis isolates from MP were hypothesized to be directly linked to KNP wildlife, based on spoligotyping. However, WGS indicated more complex epidemiology. The presence of two distinct clades which were genetically distinct (SNP distance of 19-47) and suggested multiple transmission events. Therefore, WGS provided new insight into the molecular epidemiology of the M. bovis isolates from MP and their relationship to isolates from KNP. This approach will facilitate greater understanding of M. bovis transmission at wildlife-livestock-human interfaces and advances One Health research on tuberculosis, especially across different host species.

Implications for TB control among migrants in large cities in China: A prospective population-based genomic epidemiology study in Shenzhen

Internal migrants are a challenge for TB control in large Chinese cities and understanding this epidemiology is crucial for designing effective control and prevention strategies. We conducted a prospective genomic epidemiological study of culture-positive TB patients diagnosed between June 1, 2018 and May 31, 2021 in the Longhua District of Shenzhen. Treatment status was obtained from local and national TB registries and all isolates were sequenced. Genomic clusters were defined as strains differing by ≤12 SNPs. Risk factors for clustering were identified with multivariable analysis and then Bayesian models and TransPhylo were used to infer the timing of transmission within clusters. Of the 2277 culture-positive patients, 70.1% (1596/2277) were migrants: 72.1% (1043/1446) of the migrants patients developed TB within two years of arriving in Longhua; 38.8% within 6 months of arriving; and 12.3% (104/843) had TB symptoms when they arrived. Only 15.4% of Longhua strains were in genomic clusters. More than one third (33.6%) of patients were not treated in Shenzhen but were involved in nearly one third of the recent transmission events. Clustering was associated with migrants not treated in Shenzhen, males, and teachers/trainers. TB in Longhua is prinicipally due to reactivation of infections in migrants, but a proportion may have had clinical or incipient TB upon arrival in the district. Patients diagnosed but not treated in Longhua were involved in recent local TB transmission. Controlling TB in Shenzhen will require strategies to comprehensively diagnose and treat active TB in the internal migrant population.

Implementation of national whole-genome sequencing of Mycobacterium tuberculosis, National Public Health Laboratory, Singapore, 2019-2022

The National Tuberculosis Programme (NTBP) monitors the occurrence and spread of tuberculosis (TB) and multidrug-resistant TB (MDR-TB) in Singapore. Since 2020, whole-genome sequencing (WGS) of Mycobacterium tuberculosis isolates has been performed at the National Public Health Laboratory (NPHL) for genomic surveillance, replacing spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeats analysis (MIRU-VNTR). Four thousand three hundred and seven samples were sequenced from 2014 to January 2023, initially as research projects and later developed into a comprehensive public health surveillance programme. Currently, all newly diagnosed culture-positive cases of TB in Singapore are prospectively sent for WGS, which is used to perform lineage classification, predict drug resistance profiles and infer genetic relationships between TB isolates. This paper describes NPHL's operational and technical experiences with implementing the WGS service in an urban TB-endemic setting, focusing on cluster detection and genomic drug susceptibility testing (DST). Cluster detection: WGS has been used to guide contact tracing by detecting clusters and discovering unknown transmission networks. Examples have been clusters in a daycare centre, housing apartment blocks and a horse-racing betting centre. Genomic DST: genomic DST prediction (gDST) identifies mutations in core genes known to be associated with TB drug resistance catalogued in the TBProfiler drug resistance mutation database. Mutations are reported with confidence scores according to a standardized approach referencing NPHL's internal gDST confidence database, which is adapted from the World Health Organization (WHO) TB drug mutation catalogue. Phenotypic-genomic concordance was observed for the first-line drugs ranging from 2959/2998 (98.7 %) (ethambutol) to 2983/2996 (99.6 %) (rifampicin). Aspects of internal database management, reporting standards and caveats in results interpretation are discussed.

Mycobacterium tuberculosis epidemiology in Oman: whole-genome sequencing uncovers transmission pathways

Tuberculosis (TB) originating from expatriates that hail from high TB-burden countries is hypothesized to play a role in continued TB transmission in Oman. Here, we used whole-genome sequencing (WGS) to assess national TB transmission dynamics. The annual incidence per 100,000 population per year was calculated for nationals and expatriates. A convenience sample of Mycobacterium tuberculosis (MTB) isolates from 2018 to 2019 was sequenced and analyzed with publicly available TB sequences from Bangladesh, Tanzania, the Philippines, India, and Pakistan. Relatedness was assessed by generating core-genome single nucleotide polymorphism (SNP) distances. The incidence of TB was five cases per 100,000 persons in 2018 and seven cases per 100,000 persons in 2020 (R2 = 0.34, P = 0.60). Incidence among nationals was 3.9 per 100,000 persons in 2018 and 3.5 per 100,000 persons in 2020 (R2 = 0.20, P = 0.70), and incidence among expatriates was 7.2 per 100,000 persons in 2018 and 12.7 per 100,000 persons in 2020 (R2 = 0.74, P = 0.34). Sixty-eight local MTB isolates were sequenced and analyzed with 393 global isolates. Isolates belonged to nine distinct spoligotypes. Two isolates, originating from an expatriate and an Omani national, were grouped into a WGS-based cluster (SNP distance < 12), which was corroborated by an epidemiological investigation. Relatedness of local and global isolates (SNP distance < 100) was also seen. The relatedness between MTB strains in Oman and those in expatriate countries of origin can aid inform TB control policy. Our results provide evidence that WGS can complement epidemiological analysis to achieve the End TB strategy goal in Oman. IMPORTANCE Tuberculosis (TB) incidence in Oman remains above national program control targets. TB transmission originating from expatriates from high TB-burden countries has been hypothesized to play a role. We used whole-genome sequencing (WGS) to assess TB transmission dynamics between expatriates and Omani nationals to inform TB control efforts. Available Mycobacterium tuberculosis isolates from 2018 to 2019 underwent WGS and analysis with publicly available TB sequences from Bangladesh, the Philippines, India, and Pakistan to assess for genetic relatedness. Our analysis revealed evidence of previously unrecognized transmission between an expatriate and an Omani national, which was corroborated by epidemiological investigation. Analysis of local and global isolates revealed evidence of distant relatedness between local and global isolates. Our results provide evidence that WGS can complement classic public health surveillance to inform targeted interventions to achieve the End TB strategy goal in Oman.

Model-based Analysis of Tuberculosis Genotype Clusters in the United States Reveals High Degree of Heterogeneity in Transmission and State-level Differences Across California, Florida, New York, and Texas

BACKGROUND

Reductions in tuberculosis (TB) transmission have been instrumental in lowering TB incidence in the United States. Sustaining and augmenting these reductions are key public health priorities.

METHODS

We fit mechanistic transmission models to distributions of genotype clusters of TB cases reported to the Centers for Disease Control and Prevention during 2012-2016 in the United States and separately in California, Florida, New York, and Texas. We estimated the mean number of secondary cases generated per infectious case (R0) and individual-level heterogeneity in R0 at state and national levels and assessed how different definitions of clustering affected these estimates.

RESULTS

In clusters of genotypically linked TB cases that occurred within a state over a 5-year period (reference scenario), the estimated R0 was 0.29 (95% confidence interval [CI], .28-.31) in the United States. Transmission was highly heterogeneous; 0.24% of simulated cases with individual R0 >10 generated 19% of all recent secondary transmissions. R0 estimate was 0.16 (95% CI, .15-.17) when a cluster was defined as cases occurring within the same county over a 3-year period. Transmission varied across states: estimated R0s were 0.34 (95% CI, .3-.4) in California, 0.28 (95% CI, .24-.36) in Florida, 0.19 (95% CI, .15-.27) in New York, and 0.38 (95% CI, .33-.46) in Texas.

CONCLUSIONS

TB transmission in the United States is characterized by pronounced heterogeneity at the individual and state levels. Improving detection of transmission clusters through incorporation of whole-genome sequencing and identifying the drivers of this heterogeneity will be essential to reducing TB transmission.

Developing customized stepwise MIRU-VNTR typing for tuberculosis surveillance in Georgia

Introduction

Mycobacterial Interspersed Repetitive Units–Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia.

Methods

We used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method.

Results

For the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each.

Conclusion

Our study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

Insights into the Complexity of a Dormant Mycobacterium tuberculosis Cluster Once Transmission Is Resumed

Genotyping tools help identify the complexity in Mycobacterium tuberculosis transmission clusters. We carried out a thorough analysis of the epidemiological and bacteriological complexity of a cluster in Almería, Spain. The cluster, initially associated with Moroccan migrants and with no secondary cases identified in 4 years, then reappeared in Spanish-born individuals. In one case, two Mycobacterium tuberculosis clonal variants were identified. We reanalyzed the cluster, supported by the characterization of multiple cultured isolates and respiratory specimens, whole-genome sequencing, and epidemiological case interviews. Our findings showed that the cluster, which was initially thought to have restarted activity with just a single case harboring a small degree of within-host diversity, was in fact currently growing due to coincidental reactivation of past exposures, with clonal diversity transmitted throughout the cluster. In one case, within-host diversity was amplified, probably due to prolonged diagnostic delay. IMPORTANCE The precise study of the dynamics of tuberculosis transmission in socio-epidemiologically complex scenarios may require more thorough analysis than the standard molecular epidemiology strategies. Our study illustrates the epidemiological and bacteriological complexity present in a transmission cluster in a challenging epidemiological setting with a high proportion of migrant cases. The combination of whole-genome sequencing, refined and refocused epidemiological interviews, and in-depth analysis of the bacterial composition of sputa and cultured isolates was crucial in order to correctly reinterpret the true nature of this cluster. Our global approach allowed us to reinterpret correctly the unnoticed epidemiological and bacteriological complexity involved in the Mycobacterium tuberculosis transmission event under study, which had been overlooked by the usual molecular epidemiology approaches.

What Can Genetic Relatedness Tell Us About Risk Factors for Tuberculosis Transmission?

BACKGROUND

To stop tuberculosis (TB), the leading infectious cause of death globally, we need to better understand transmission risk factors. Although many studies have identified associations between individual-level covariates and pathogen genetic relatedness, few have identified characteristics of transmission pairs or explored how closely covariates associated with genetic relatedness mirror those associated with transmission.

METHODS

We simulated a TB-like outbreak with pathogen genetic data and estimated odds ratios (ORs) to correlate each covariate and genetic relatedness. We used a naive Bayes approach to modify the genetic links and nonlinks to resemble the true links and nonlinks more closely and estimated modified ORs with this approach. We compared these two sets of ORs with the true ORs for transmission. Finally, we applied this method to TB data in Hamburg, Germany, and Massachusetts, USA, to find pair-level covariates associated with transmission.

RESULTS

Using simulations, we found that associations between covariates and genetic relatedness had the same relative magnitudes and directions as the true associations with transmission, but biased absolute magnitudes. Modifying the genetic links and nonlinks reduced the bias and increased the confidence interval widths, more accurately capturing error. In Hamburg and Massachusetts, pairs were more likely to be probable transmission links if they lived in closer proximity, had a shorter time between observations, or had shared ethnicity, social risk factors, drug resistance, or genotypes.

CONCLUSIONS

We developed a method to improve the use of genetic relatedness as a proxy for transmission, and aid in understanding TB transmission dynamics in low-burden settings.

Using Whole-genome Sequencing to Determine the Timing of Secondary Tuberculosis in British Columbia, Canada

Combined with epidemiological data, whole-genome sequencing (WGS) can help better resolve individual tuberculosis (TB) transmission events to a degree not possible with traditional genotyping. We combine WGS data with patient-level data to calculate the timing of secondary TB among contacts of people diagnosed with active TB in British Columbia, Canada.

Use of next generation sequencing technologies for the diagnosis and epidemiology of infectious diseases

For the first time, next generation sequencing technologies provide access to genomic information at a price and scale that allow their implementation in routine clinical practice and epidemiology. While there are still many obstacles to their implementation, there are also multiple examples of their major advantages compared with previous methods. Their main advantage is that a single determination allows epidemiological information on the causative microorganism to be obtained simultaneously, as well as its resistance profile, although these advantages vary according to the pathogen under study. This review discusses several examples of the clinical and epidemiological use of next generation sequencing applied to complete genomes and microbiomes and reflects on its future in clinical practice.

Molecular epidemiology and drug-resistance of tuberculosis in Luodian revealed by whole genome sequencing

In this study, we aimed to investigate the molecular epidemiology and drug-resistance profiles of tuberculosis (TB) in Luodian, an area with highest TB incidence and limited healthcare resources in Guizhou, China. The passive case finding strategy was used to identify suspected pulmonary TB with symptoms, and individuals with positive Mycobacterium tuberculosis (MTB) culture were enrolled from May 22, 2018 to April 21, 2019. All the 107 cases except three came from nine towns, including 55.1% from Longping and Bianyang. The phylogeny tree showed that 53.3% of strains were Lineage 2 (Beijing genotype), while 46.7% were Lineage 4 (Euro-American genotype). Among Lineage 2 strains, 66.7% were of "modern" Beijing type. Seven clusters with genomic distance within 12 SNPs were identified. The clusters included 14 strains, accounting for a clustering rate of 13.1%. The distance separating the clustered cases was between 2.1 and 71.0 km (Km), with an average paired distance of 21.8 Km (interquartile range, 2.8-38.0 Km). Based on the gene mutations associated with drug-resistance, we predicted that 4.8% of strains were resistant to isoniazid, 3.7% to rifampicin, and 3.7% to streptomycin; only one strain (0.9%) had multidrug resistance (MDR). This study found low drug-resistance rates in Luodian, and the sub-lineage of the "modern" Beijing branch has recent expansion in Luodian. This work may also serve as a genomic baseline to assess the evolution and spread of MTB in Guizhou.

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.

Assessment of closely related Mycobacterium tuberculosis variants with different transmission success and in vitro infection dynamics

Whole genome sequencing (WGS) is able to differentiate closely related Mycobacterium tuberculosis variants within the same transmission cluster. Our aim was to evaluate if this higher discriminatory power may help identify and characterize more actively transmitted variants and understand the factors behind their success. We selected a robust MIRU-VNTR-defined cluster from Almería, Spain (22 cases throughout 2003–2019). WGS allowed discriminating, within the same epidemiological setting, between a successfully transmitted variant and seven closely related variants that did not lead to secondary cases, or were involved in self-limiting transmission (one single secondary case). Intramacrophagic growth of representative variants was evaluated in an in vitro infection model using U937 cells. Intramacrophage multiplication ratios (CFUs at Day 4/CFUs at Day 0) were higher for the actively transmitted variant (range 5.3–10.7) than for the unsuccessfully transmitted closely related variants (1.5–3.95). Two SNPs, mapping at the DNA binding domain of DnaA and at kdpD, were found to be specific of the successful variant.

Effects of migration on tuberculosis epidemiological indicators in low and medium tuberculosis incidence countries: A systematic review

BACKGROUND

Tuberculosis (TB) remains one of the top ten causes of death each year globally. While the risk of migrant TB is linked to the TB incidence in their country of origin, the migration process can increase the TB risk.

OBJECTIVE

We aimed to synthesis the evidence on key differences in the epidemiological profile of TB between migrants from high TB incidence birth countries and non-migrants resident in low to medium incidence TB countries.

METHODS

We conducted a systematic review where the population was all active TB cases in countries with low to medium TB incidence (<40/100,000 population), the exposure was migration to a low or medium TB incidence country and the comparator was non-migrant TB cases in low or medium incidence countries. Overall proportions were compared between migrants and non-migrants, using Fisher's exact test. Meta-analysis of proportions was carried out for the primary outcome (active TB) while meta-analyses of odds ratios (ORs) were performed using a random effects model for secondary outcomes; sputum-smear positivity, any first line drug resistance, multi-drug resistance (MDR), clustered cases, HIV coinfections and successful treatment. Heterogeneity was evaluated and sources were investigated using subgroup and sensitivity analysis.

RESULTS

Significant differences were found in the overall proportions of high TB incidence migrants and non-migrants for MDR cases, clustered cases, HIV coinfections and successful treatment, as well as a significant difference in the OR among MDR cases (3.91).

CONCLUSION

This review has demonstrated significant differences in key epidemiological indicators between high TB incidence migrants and non-migrants, indicating policy implications.

Clinical and public health utility of Mycobacterium tuberculosis whole genome sequencing

The World Health Organization (WHO) estimates that around 10 million people develop tuberculosis (TB) every year, with 1.5 million deaths attributed to TB in 2019 (World Health Organization, 2020). The majority of the disease burden occurs in low-income countries, where access to diagnostics and tailored treatment remains problematic. The current COVID-19 pandemic further threatens to impact global TB control by diverting resources, reducing notifications and hence significantly increasing deaths attributable to TB (World Health Organization, 2020). Whole genome sequencing (WGS) is becoming increasingly accessible, and has particular value in the diagnosis and management of TB disease (Cabibbe et al., 2018; Meehan et al., 2019). Not only does it have the potential to give more rapid and complete information on drug-resistance, but the high discriminatory power it offers allows detection of clusters and transmission pathways, as well as likely contamination events, mixed infections and to differentiate between re-infection and relapse with much greater confidence than previous typing methods.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Capture and visualization of live Mycobacterium tuberculosis bacilli from tuberculosis patient bioaerosols

Interrupting transmission is an attractive anti-tuberculosis (TB) strategy but it remains underexplored owing to our poor understanding of the events surrounding transfer of Mycobacterium tuberculosis (Mtb) between hosts. Determining when live, infectious Mtb bacilli are released and by whom has proven especially challenging. Consequently, transmission chains are inferred only retrospectively, when new cases are diagnosed. This process, which relies on molecular analyses of Mtb isolates for epidemiological fingerprinting, is confounded by the prolonged infectious period of TB and the potential for transmission from transient exposures. We developed a Respiratory Aerosol Sampling Chamber (RASC) equipped with high-efficiency filtration and sampling technologies for liquid-capture of all particulate matter (including Mtb) released during respiration and non-induced cough. Combining the mycobacterial cell wall probe, DMN-trehalose, with fluorescence microscopy of RASC-captured bioaerosols, we detected and quantified putative live Mtb bacilli in bioaerosol samples arrayed in nanowell devices. The RASC enabled non-invasive capture and isolation of viable Mtb from bioaerosol within 24 hours of collection. A median 14 live Mtb bacilli (range 0-36) were isolated in single-cell format from 90% of confirmed TB patients following 60 minutes bioaerosol sampling. This represented a significant increase over previous estimates of transmission potential, implying that many more organisms might be released daily than commonly assumed. Moreover, variations in DMN-trehalose incorporation profiles suggested metabolic heterogeneity in aerosolized Mtb. Finally, preliminary analyses indicated the capacity for serial image capture and analysis of nanowell-arrayed bacilli for periods extending into weeks. These observations support the application of this technology to longstanding questions in TB transmission including the propensity for asymptomatic transmission, the impact of TB treatment on Mtb bioaerosol release, and the physiological state of aerosolized bacilli.

Whole genomic sequencing based genotyping reveals a specific X3 sublineage restricted to Mexico and related with multidrug resistance

Whole genome sequencing (WGS) has been shown to be superior to traditional procedures of genotyping in tuberculosis (TB), nevertheless, reports of its use in drug resistant TB (DR-TB) isolates circulating in Mexico, are practically unknown. Considering the above the main of this work was to identify and characterize the lineages and genomic transmission clusters present in 67 DR-TB isolates circulating in southeastern Mexico. The results show the presence of three major lineages: L1 (3%), L2 (3%) and L4 (94%), the last one included 16 sublineages. Sublineage 4.1.1.3 (X3) was predominant in 18 (27%) of the isolates, including one genomic cluster, formed by eleven multidrug resistant isolates and sharing the SIT 3278, which seems to be restricted to Mexico. By the use of WGS, it was possible to identify the high prevalence of L4 and a high number of sublineages circulating in the region, also was recognized the presence of a novel X3 sublineage, formed exclusively by multidrug resistant isolates and with restrictive circulation in Mexico for at least the past 17 years.

The impact of HIV infection on tuberculosis transmission in a country with low tuberculosis incidence: a national retrospective study using molecular epidemiology

BACKGROUND

HIV is known to increase the likelihood of reactivation of latent tuberculosis to active TB disease; however, its impact on tuberculosis infectiousness and consequent transmission is unclear, particularly in low-incidence settings.

METHODS

National surveillance data from England, Wales and Northern Ireland on tuberculosis cases in adults from 2010 to 2014, strain typed using 24-locus mycobacterial-interspersed-repetitive-units-variable-number-tandem-repeats was used retrospectively to identify clusters of tuberculosis cases, subdivided into 'first' and 'subsequent' cases. Firstly, we used zero-inflated Poisson regression models to examine the association between HIV status and the number of subsequent clustered cases (a surrogate for tuberculosis infectiousness) in a strain type cluster. Secondly, we used logistic regression to examine the association between HIV status and the likelihood of being a subsequent case in a cluster (a surrogate for recent acquisition of tuberculosis infection) compared to the first case or a non-clustered case (a surrogate for reactivation of latent infection).

RESULTS

We included 18,864 strain-typed cases, 2238 were the first cases of clusters and 8471 were subsequent cases. Seven hundred and fifty-nine (4%) were HIV-positive. Outcome 1: HIV-positive pulmonary tuberculosis cases who were the first in a cluster had fewer subsequent cases associated with them (mean 0.6, multivariable incidence rate ratio [IRR] 0.75 [0.65-0.86]) than those HIV-negative (mean 1.1). Extra-pulmonary tuberculosis (EPTB) cases with HIV were less likely to be the first case in a cluster compared to HIV-negative EPTB cases. EPTB cases who were the first case had a higher mean number of subsequent cases (mean 2.5, IRR (3.62 [3.12-4.19]) than those HIV-negative (mean 0.6). Outcome 2: tuberculosis cases with HIV co-infection were less likely to be a subsequent case in a cluster (odds ratio 0.82 [0.69-0.98]), compared to being the first or a non-clustered case.

CONCLUSIONS

Outcome 1: pulmonary tuberculosis-HIV patients were less infectious than those without HIV. EPTB patients with HIV who were the first case in a cluster had a higher number of subsequent cases and thus may be markers of other undetected cases, discoverable by contact investigations. Outcome 2: tuberculosis in HIV-positive individuals was more likely due to reactivation than recent infection, compared to those who were HIV-negative.

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

Extensive Homoplasy but No Evidence of Convergent Evolution of Repeat Numbers at MIRU Loci in Modern Mycobacterium tuberculosis Lineages

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis. MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis. In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes.

Genomic epidemiology of Mycobacterium tuberculosis in Santa Catarina, Southern Brazil

Mycobacterium tuberculosis (M.tb), the pathogen responsible for tuberculosis (TB) poses as the major cause of death among infectious diseases. The knowledge about the molecular diversity of M.tb enables the implementation of more effective surveillance and control measures and, nowadays, Whole Genome Sequencing (WGS) holds the potential to produce high-resolution epidemiological data in a high-throughput manner. Florianópolis, the state capital of Santa Catarina (SC) in south Brazil, shows a high TB incidence (46.0/100,000). Here we carried out a WGS-based evaluation of the M.tb strain diversity, drug-resistance and ongoing transmission in the capital metropolitan region. Resistance to isoniazid, rifampicin, streptomycin was identified respectively in 4.0% (n = 6), 2.0% (n = 3) and 1.3% (n = 2) of the 151 studied strains by WGS. Besides, resistance to pyrazinamide and ethambutol was detected in 0.7% (n = 1) and reistance to ethionamide and fluoroquinolone (FQ) in 1.3% (n = 2), while a single (0.7%) multidrug-resistant (MDR) strain was identified. SNP-based typing classified all isolates into M.tb Lineage 4, with high proportion of sublineages LAM (60.3%), T (16.4%) and Haarlem (7.9%). The average core-genome distance between isolates was 420.3 SNPs, with 43.7% of all isolates grouped across 22 genomic clusters thereby showing the presence of important ongoing TB transmission events. Most clusters were geographically distributed across the study setting which highlights the need for an urgent interruption of these large transmission chains. The data conveyed by this study shows the presence of important and uncontrolled TB transmission in the metropolitan area and provides precise data to support TB control measures in this region.

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.

In vitro activity of apramycin against multidrug-, carbapenem- and aminoglycoside-resistant Enterobacteriaceae and Acinetobacter baumannii

Objectives

Widespread antimicrobial resistance often limits the availability of therapeutic options to only a few last-resort drugs that are themselves challenged by emerging resistance and adverse side effects. Apramycin, an aminoglycoside antibiotic, has a unique chemical structure that evades almost all resistance mechanisms including the RNA methyltransferases frequently encountered in carbapenemase-producing clinical isolates. This study evaluates the in vitro activity of apramycin against multidrug-, carbapenem- and aminoglycoside-resistant Enterobacteriaceae and Acinetobacter baumannii, and provides a rationale for its superior antibacterial activity in the presence of aminoglycoside resistance determinants.

Methods

A thorough antibacterial assessment of apramycin with 1232 clinical isolates from Europe, Asia, Africa and South America was performed by standard CLSI broth microdilution testing. WGS and susceptibility testing with an engineered panel of aminoglycoside resistance-conferring determinants were used to provide a mechanistic rationale for the breadth of apramycin activity.

Results

MIC distributions and MIC₉₀ values demonstrated broad antibacterial activity of apramycin against Escherichia coli, Klebsiella pneumoniae, Enterobacter spp., Morganella morganii, Citrobacter freundii, Providencia spp., Proteus mirabilis, Serratia marcescens and A. baumannii. Genotypic analysis revealed the variety of aminoglycoside-modifying enzymes and rRNA methyltransferases that rendered a remarkable proportion of clinical isolates resistant to standard-of-care aminoglycosides, but not to apramycin. Screening a panel of engineered strains each with a single well-defined resistance mechanism further demonstrated a lack of cross-resistance to gentamicin, amikacin, tobramycin and plazomicin.

Conclusions

Its superior breadth of activity renders apramycin a promising drug candidate for the treatment of systemic Gram-negative infections that are resistant to treatment with other aminoglycoside antibiotics.

Tuberculosis in elephants: Origins and evidence of interspecies transmission

Tuberculosis (TB) is a devastating disease in elephants caused by either Mycobacterium tuberculosis or M. bovis. It is an ancient disease, and TB in elephants was first reported over two millennia ago in Sri Lanka. Outbreaks of TB worldwide, in captive and free-ranging elephant populations, have been recorded. Interspecies transmission of TB among elephants and humans has been confirmed in several geographic localities using spoligotyping, MIRU-VNTR analysis, and/or comparative genomics. Active surveillance of TB in wild and captive elephants and their handlers is necessary to prevent TB transmission at the elephant-human interface and to aid in the conservation of Asian and African elephants. In this review, we present an overview of diagnosis, reports of TB outbreaks in the past 25 years, TB in wild elephants, its transmission, and possible prevention and control strategies that can be applied at the elephant-human interface.

Genome sequence analysis of multidrug-resistant Mycobacterium tuberculosis from Malaysia

Mycobacterium tuberculosis (MTB) is commonly used as a model to study pathogenicity and multiple drug resistance in bacteria. These MTB characteristics are highly dependent on the evolution and phylogeography of the bacterium. In this paper, we describe 15 new genomes of multidrug-resistant MTB (MDRTB) from Malaysia. The assessments and annotations on the genome assemblies suggest that strain differences are due to lineages and horizontal gene transfer during the course of evolution. The genomes show mutations listed in current drug resistance databases and global MTB collections. This genome data will augment existing information available for comparative genomic studies to understand MTB drug resistance mechanisms and evolution.

Measurement(s)	DNA • genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Factor Type(s)	strain
Sample Characteristic - Organism	Mycobacterium tuberculosis
Sample Characteristic - Location	Malaysia

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12130008

A model of tuberculosis clustering in low incidence countries reveals more transmission in the United Kingdom than the Netherlands between 2010 and 2015

Tuberculosis (TB) remains a public health threat in low TB incidence countries, through a combination of reactivated disease and onward transmission. Using surveillance data from the United Kingdom (UK) and the Netherlands (NL), we demonstrate a simple and predictable relationship between the probability of observing a cluster and its size (the number of cases with a single genotype). We demonstrate that the full range of observed cluster sizes can be described using a modified branching process model with the individual reproduction number following a Poisson lognormal distribution. We estimate that, on average, between 2010 and 2015, a TB case generated 0.41 (95% CrI 0.30,0.60) secondary cases in the UK, and 0.24 (0.14,0.48) secondary cases in the NL. A majority of cases did not generate any secondary cases. Recent transmission accounted for 39% (26%,60%) of UK cases and 23%(13%,37%) of NL cases. We predict that reducing UK transmission rates to those observed in the NL would result in 538(266,818) fewer cases annually in the UK. In conclusion, while TB in low incidence countries is strongly associated with reactivated infections, we demonstrate that recent transmission remains sufficient to warrant policies aimed at limiting local TB spread.

Characterization of Polymorphisms Associated with Multidrug-Resistant Tuberculosis by Whole Genomic Sequencing: A Preliminary Report from Mexico

Whole genome sequencing (WGS) has been proposed as a tool for the diagnosis of drug resistance in tuberculosis (TB); however, there have been few studies on its effectiveness in countries with significantly high drug resistance rates. This study therefore aimed to evaluate the effectiveness of WGS to identify mutations related to drug resistance in TB isolates from an endemic region of Mexico. The results showed that, of 35 multidrug-resistant isolates analyzed, the values of congruence found between the phenotypic drug susceptibility testing and polymorphisms were 94% for isoniazid, 97% for rifampicin, 90% for ethambutol, and 82% for pyrazinamide. It was also possible to identify eight isolates as potential pre-extensive drug resistant (XDR) and one as XDR. Twenty nine isolates were classified within L4 and two transmission clusters were identified. The results show the potential utility of WGS for predicting resistance against first- and second-line drugs, as well as providing a phylogenetic characterization of TB drug-resistant isolates circulating in Mexico.

Computational pan-genome mapping and pairwise SNP-distance improve detection of Mycobacterium tuberculosis transmission clusters

Next-generation sequencing based base-by-base distance measures have become an integral complement to epidemiological investigation of infectious disease outbreaks. This study introduces PANPASCO, a computational pan-genome mapping based, pairwise distance method that is highly sensitive to differences between cases, even when located in regions of lineage specific reference genomes. We show that our approach is superior to previously published methods in several datasets and across different Mycobacterium tuberculosis lineages, as its characteristics allow the comparison of a high number of diverse samples in one analysis—a scenario that becomes more and more likely with the increased usage of whole-genome sequencing in transmission surveillance.

Tuberculosis still is a threat to global health. It is essential to detect and interrupt transmissions to stop the spread of this infectious disease. With the rising use of next-generation sequencing methods, its application in the surveillance of Mycobacterium tuberculosis has become increasingly important in the last years. The main goal of molecular surveillance is the identification of patient-patient transmission and cluster detection. The mutation rate of M. tuberculosis is very low and stable. Therefore, many existing methods for comparative analysis of isolates provide inadequate results since their resolution is too limited. There is a need for a method that takes every detectable difference into account. We developed PANPASCO, a novel approach for comparing pairs of isolates using all genomic information available for each pair. We combine improved SNP-distance calculation with the use of a pan-genome incorporating more than 100 M. tuberculosis reference genomes representing lineages 1-4 for read mapping prior to variant detection. We thereby enable the collective analysis and comparison of similar and diverse isolates associated with different M. tuberculosis strains.

Diversity of Mycobacterium tuberculosis in the Middle Fly District of Western Province, Papua New Guinea: microbead-based spoligotyping using DNA from Ziehl-Neelsen-stained microscopy preparations

Tuberculosis remains the world's leading cause of death from an infectious agent, and is a serious health problem in Papua New Guinea (PNG) with an estimated 36,000 new cases each year. This study describes the genetic diversity of Mycobacterium tuberculosis among tuberculosis patients in the Balimo/Bamu region in the Middle Fly District of Western Province in PNG, and investigates rifampicin resistance-associated mutations. Archived Ziehl-Neelsen-stained sputum smears were used to conduct microbead-based spoligotyping and assess genotypic resistance. Among the 162 samples included, 80 (49.4%) generated spoligotyping patterns (n = 23), belonging predominantly to the L2 Lineage (44%) and the L4 Lineage (30%). This is consistent with what has been found in other PNG regions geographically distant from Middle Fly District of Western Province, but is different from neighbouring South-East Asian countries. Rifampicin resistance was identified in 7.8% of the successfully sequenced samples, with all resistant samples belonging to the L2/Beijing Lineage. A high prevalence of mixed L2/L4 profiles was suggestive of polyclonal infection in the region, although this would need to be confirmed. The method described here could be a game-changer in resource-limited countries where large numbers of archived smear slides could be used for retrospective (and prospective) studies of M. tuberculosis genetic epidemiology.

Molecular Epidemiology of Tuberculosis in British Columbia, Canada: A 10-Year Retrospective Study

Background

Understanding regional molecular epidemiology allows for the development of more efficient tuberculosis prevention strategies in low-incidence settings.

Methods

We analyzed 24-locus mycobacterial interspersed repetitive-unit–variable-number tandem repeat (MIRU-VNTR) genotyping for 2290 Mycobacterium tuberculosis clinical isolates collected in the province of British Columbia (BC), Canada, in 2005–2014. Laboratory data for each isolate were linked to case-level clinical and demographic data. These data were used to describe the molecular epidemiology of tuberculosis across the province.

Results

We detected >1500 distinct genotypes across the 4 major M. tuberculosis lineages, reflecting BC’s diverse population. Disease site and clustering rates varied across lineages, and MIRU-VNTR was used to group the 2290 isolates into 189 clusters (2–70 isolates per cluster), with an overall clustering rate of 42.4% and an estimated local transmission rate of 34.1%. Risk factors for clustering varied between Canadian-born and foreign-born individuals; the former had increased odds (odds ratio, 7.8; 95% confidence interval [CI], 6.2–9.6) of belonging to a genotypic cluster, although nearly one-quarter of clusters included both Canadian- and foreign-born persons. Large clusters (≥10 cases) occurred more frequently within the M. tuberculosis Euro-American lineage, and individual-level risk factors associated with belonging to a large cluster included being Canadian born (adjusted odds ratio, 3.3; 95% CI, 2.3–4.8), residing in a rural area (2.3; 1.2–4.5), and illicit drug use (2.0; 1.2–3.4).

Conclusions

Although tuberculosis in BC largely arises through reactivation of latent tuberculosis in foreign-born persons, locally transmitted infections occur in discrete populations with distinct disease and risk factor profiles, representing groups for targeted interventions.

Molecular epidemiology of M. tuberculosis in Ethiopia: A systematic review and meta-analysis

The molecular epidemiology of Mycobacterium tuberculosis (M. tuberculosis, Mtb) is poorly documented in Ethiopia. The data that exists has not yet been collected in an overview metadata form. Thus, this review summarizes available literature on the genomic diversity, geospatial distribution and transmission patterns of Mtb lineages (L) and sublineages in Ethiopia. Spoligotyping and Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR) based articles were identified from MEDLINE via PubMed and Scopus. The last date of article search was done on 12th February 2019. Articles were selected following the PRISMA flow diagram. The proportion of (sub)lineages was summarized at national level and further disaggregated by region. Clustering and recent transmission index (RTI) were determined using metan command and random effect meta-analysis model. The meta-analysis was computed using Stata 14 (Stata Corp. College Station, TX, USA). Among 4371 clinical isolates, 99.5% were Mtb and 0.5% were M. bovis. Proportionally, L4, L3, L1 and L7 made up 62.3%, 21.7%, 7.9% and 3.4% of the total isolates, respectively. Among sublineages, L4.2. ETH/SIT149, L4.10/SIT53, L3. ETH1/SIT25 and L4.6/SIT37 were the leading clustered isolates accounting for 14.4%, 9.7%, 7.2% and 5.5%, respectively. Based on MIRU-VNTR, the rate of clustering was 41% and the secondary case rate from a single source case was estimated at 29%. Clustering and recent transmission index was higher in eastern and southwestern Ethiopia compared with the northwestern part of the country. High level of genetic diversity with a high rate of clustering was noted which collectively mirrored the phenomena of micro-epidemics and super-spreading. The largest set of clustered strains deserves special attention and further characterization using whole genome sequencing (WGS) to better understand the evolution, genomic diversity and transmission dynamics of Mtb.

Whole-genome sequencing for rapid, reliable and routine investigation of Mycobacterium tuberculosis transmission in local communities

Contact investigations following the diagnosis of active tuberculosis (TB) are paramount for the control of the disease. Epidemiological data are very powerful for contact tracing but might be delayed and/or difficult to integrate, especially in the setting of multiple contact-tracing investigations. The aim of this study was to address the added-value of whole-genome sequencing (WGS) to routine local TB surveillance systems. From November 2016 to July 2017, the local TB surveillance system identified three clusters that could constitute a unique larger outbreak. Epidemiological and clinical information were integrated with WGS genotyping data of Mycobacterium tuberculosis strains obtained using a simple DNA extraction method coupled with sequencing using an Illumina MiSeq platform and an in-house bioinformatics pipeline for single nucleotide polymorphism (SNP) analysis. Epidemiological investigations identified three putative TB clusters potentially interrelated including eight patients with active TB. Seven M. tuberculosis isolates were available and analysed by WGS. Using a 5-SNP threshold to define recent transmission, WGS-based genotyping supported the occurrence of the three clusters as well as a link between clusters 1 and 2 (SNP ≤1), constituting a larger outbreak. This outbreak was clearly delineated by refuting a potential link with the third cluster (SNP >500). Genotyping data did not support the belonging of patient 7 to any studied cluster. This study illustrates the usefulness of WGS genotyping for routine TB surveillance systems in local communities to rapidly confirm or disprove epidemiological hypotheses and delineate TB clusters, especially in the context of multiple contact-tracing investigations.

Whole Genome Sequencing for Surveillance of Diphtheria in Low Incidence Settings

Corynebacterium diphtheriae (C. diphtheriae) is a relatively rare pathogen in most Western countries. While toxin producing strains can cause pharyngeal diphtheria with potentially fatal outcomes, the more common presentation is wound infections. The diphtheria toxin is encoded on a prophage and can also be carried by Corynebacterium ulcerans and Corynebacterium pseudotuberculosis. Currently, across Europe, infections are mainly diagnosed in travelers and refugees from regions where diphtheria is more endemic, patients from urban areas with poor hygiene, and intravenous drug users. About half of the cases are non-toxin producing isolates. Rapid identification of the bacterial pathogen and toxin production is a critical element of patient and outbreak management. Beside the immediate clinical management of the patient, public health agencies should be informed of toxigenic C. diphtheriae diagnoses as soon as possible. The collection of case-related epidemiological data from the patient is often challenging due to language barriers and social circumstances. However, information on patient contacts, vaccine status and travel/refugee route, where appropriate, is critical, and should be documented. In addition, isolates should be characterized using high resolution typing, in order to identify transmissions and outbreaks. In recent years, whole genome sequencing (WGS) has become the gold standard of high-resolution typing methods, allowing detailed investigations of pathogen transmissions. De-centralized sequencing strategies with redundancy in sequencing capacities, followed by data exchange may be a valuable future option, especially since WGS becomes more available and portable. In this context, the sharing of sequence data, using public available platforms, is essential. A close interaction between microbiology laboratories, treating physicians, refugee centers, social workers, and public health officials is a key element in successful management of suspected outbreaks. Analyzing bacterial isolates at reference centers may further help to provide more specialized microbiological techniques and to standardize information, but this is also more time consuming during an outbreak. Centralized communication strategies between public health agencies and laboratories helps considerably in establishing and coordinating effective surveillance and infection control. We review the current literature on high-resolution typing of C. diphtheriae and share our own experience with the coordination of a Swiss-German outbreak.

Status and potential of bacterial genomics for public health practice: a scoping review

BACKGROUND

Next-generation sequencing (NGS) is increasingly being translated into routine public health practice, affecting the surveillance and control of many pathogens. The purpose of this scoping review is to identify and characterize the recent literature concerning the application of bacterial pathogen genomics for public health practice and to assess the added value, challenges, and needs related to its implementation from an epidemiologist's perspective.

METHODS

In this scoping review, a systematic PubMed search with forward and backward snowballing was performed to identify manuscripts in English published between January 2015 and September 2018. Included studies had to describe the application of NGS on bacterial isolates within a public health setting. The studied pathogen, year of publication, country, number of isolates, sampling fraction, setting, public health application, study aim, level of implementation, time orientation of the NGS analyses, and key findings were extracted from each study. Due to a large heterogeneity of settings, applications, pathogens, and study measurements, a descriptive narrative synthesis of the eligible studies was performed.

RESULTS

Out of the 275 included articles, 164 were outbreak investigations, 70 focused on strategy-oriented surveillance, and 41 on control-oriented surveillance. Main applications included the use of whole-genome sequencing (WGS) data for (1) source tracing, (2) early outbreak detection, (3) unraveling transmission dynamics, (4) monitoring drug resistance, (5) detecting cross-border transmission events, (6) identifying the emergence of strains with enhanced virulence or zoonotic potential, and (7) assessing the impact of prevention and control programs. The superior resolution over conventional typing methods to infer transmission routes was reported as an added value, as well as the ability to simultaneously characterize the resistome and virulome of the studied pathogen. However, the full potential of pathogen genomics can only be reached through its integration with high-quality contextual data.

CONCLUSIONS

For several pathogens, it is time for a shift from proof-of-concept studies to routine use of WGS during outbreak investigations and surveillance activities. However, some implementation challenges from the epidemiologist's perspective remain, such as data integration, quality of contextual data, sampling strategies, and meaningful interpretations. Interdisciplinary, inter-sectoral, and international collaborations are key for an appropriate genomics-informed surveillance.

A nationwide study of multidrug-resistant tuberculosis in Portugal 2014-2017 using epidemiological and molecular clustering analyses

Background

Increasing multidrug-resistant tuberculosis (MDR-TB) incidence is a major threat against TB eradication worldwide. We aim to conduct a detailed MDR-TB study in Portugal, an European country with endemic TB, combining genetic analysis and epidemiological data, in order to assess the efficiency of public health containment of MRD-TB in the country.

Methods

We used published MIRU-VNTR data, that we reanalysed using a phylogenetic analysis to better describe MDR-TB cases transmission occurring in Portugal from 2014 to 2017, further enriched with epidemiological data of these cases.

Results

We show an MDR-TB transmission scenario, where MDR strains likely arose and are transmitted within local chains. 63% of strains were clustered, suggesting high primary transmission (estimated as 50% using MIRU-VNTR data and 15% considering epidemiological links). These values are higher than those observed across Europe and even for sensitive strains in Portugal using similar methodologies. MDR-TB cases are associated with individuals born in Portugal and evolutionary analysis suggests a local evolution of strains. Consistently the sublineage LAM, the most common in sensitive strains in Europe, is the more frequent in Portugal in contrast with the remaining European MDR-TB picture where immigrant-associated Beijing strains are more common.

Conclusions

Despite efforts to track and contain MDR-TB strains in Portugal, their transmission patterns are still as uncontrolled as that of sensitive strains, stressing the need to reinforce surveillance and containment strategies.

Electronic supplementary material

The online version of this article (10.1186/s12879-019-4189-7) contains supplementary material, which is available to authorized users.

Whole genomic sequencing as a tool for diagnosis of drug and multidrug-resistance tuberculosis in an endemic region in Mexico

Background

Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico.

Methods

WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized.

Findings

Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition.

Conclusions

The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis.

The problem with defining foreign birth as a risk factor in tuberculosis epidemiology studies

OBJECTIVE

To examine how stratifying persons born outside Canada according to tuberculosis (TB) incidence in their birth country and other demographic factors refines our understanding of TB epidemiology and local TB transmission.

BACKGROUND

Population-level TB surveillance programs and research studies in low incidence settings often report all persons born outside the country in which the study is conducted as "foreign-born"-a single label for a highly diverse population with variable TB risks. This may mask important TB epidemiologic trends and not accurately reflect local transmission patterns.

METHODS

We used population-level data from two large cohorts in British Columbia (BC), Canada: an immigration cohort (n = 337,492 permanent residents to BC) and a genotyping cohort (n = 2290 culture-confirmed active TB cases). We stratified active TB case counts, incidence rates, and genotypic clustering (an indicator of TB transmission) in BC by birth country TB incidence, age at immigration, and years since arrival.

RESULTS

Persons from high-incidence countries had a 12-fold higher TB incidence than those emigrating from low-incidence settings. Estimates of local transmission, as captured by genotyping, versus reactivation of latent TB infection acquired outside Canada varied when data were stratified by birthplace TB incidence, as did patient-level characteristics of individuals in each group, such as age and years between immigration and diagnosis.

CONCLUSION

Categorizing persons beyond simply "foreign-born", particularly in the context of TB epidemiologic and molecular data, is needed for a more accurate understanding of TB rates and patterns of transmission.

Whole-Genome Sequencing as Tool for Investigating International Tuberculosis Outbreaks: A Systematic Review

Background: Whole-genome sequencing (WGS) can support the investigation of tuberculosis (TB) outbreaks. The technique has been applied to estimate the timing and directionality of transmission and to exclude cases from an investigation. This review assesses how WGS was applied in international outbreak investigations and discusses the advantages and challenges of the application of WGS.

Methods: Databases were searched for reports on international TB outbreak investigations. Information was extracted on: Why was WGS applied?; How was WGS applied?; Organizational issues; WGS methodology; What was learned/what were the implications of the WGS investigation?; and challenges and lessons learned.

Results: Three studies reporting on international outbreak investigations were identified. Retrospective WGS sequencing was performed in all studies and prospective typing in two to study TB transmission. In one study, WGS data were produced centrally (i.e., in one laboratory) and analysis was done centrally. In two studies, WGS data production was done in a decentralized manner, and analysis was centralized in one laboratory. Three groups of professionals were involved in the international outbreak investigation: public health authorities, laboratory experts, and clinicians. The reported WGS methodology applied differed between the studies in some aspects, e.g., sequencing platform; quality measures, percentage of the reference genome covered, and the mean genomic coverage; analysis, use of a reference genome or de novo assembly; and software used for alignment and analysis. In all three studies, in-house scripts were used for variance calling, and the single nucleotide polymorphism (SNP) approach was used for analysis. All outbreak investigation reports stated that WGS refuted suspected transmission events and provided supporting evidence for epidemiological data. Several challenges were reported of which most were not related to WGS. The only challenge related to WGS was the timeframe of getting WGS data if WGS is not routinely performed.

Conclusions: WGS was considered a useful addition in international TB outbreak investigations. Further standardization of the WGS methodology and good structures for international collaboration and coordination are needed to take full advantage of this new technology. Whether the use of WGS results in earlier detection of cases and thus limits transmission still needs to be determined.

Role and value of whole genome sequencing in studying tuberculosis transmission

BACKGROUND

Tuberculosis (TB) remains a serious public health threat worldwide. Theoretically ultimate resolution of whole genome sequencing (WGS) for Mycobacterium tuberculosis complex (MTBC) strain classification makes this technology very attractive for epidemiological investigations.

OBJECTIVES

To summarize the evidence available in peer-reviewed publications on the role and place of WGS in detection of TB transmission.

SOURCES

A total of 69 peer-reviewed publications identified in Pubmed database.

CONTENT

Evidence from >30 publications suggests that a cut-off value of fewer than six single nucleotide polymorphisms between strains efficiently excludes cases that are not the result of recent transmission and could be used for the identification of drug-sensitive isolates involved in direct human-to-human TB transmission. Sensitivity of WGS to identify epidemiologically linked isolates is high, reaching 100% in eight studies with specificity (17%-95%) highly dependent on the settings. Drug resistance and specific phylogenetic lineages may be associated with accelerated mutation rates affecting genetic distances. WGS can be potentially used to distinguish between true relapses and re-infections but in high-incidence low-diversity settings this would require consideration of epidemiological links and minority alleles. Data from four studies looking into within-host diversity highlight a need for developing criteria for acceptance or rejection of WGS relatedness results depending on the proportion of minority alleles.

IMPLICATIONS

WGS will potentially allow for more targeted public health actions preventing unnecessary investigations of false clusters. Consensus on standardization of raw data quality control processing criteria, analytical pipelines and reporting language is yet to be reached.

Universal genotyping reveals province-level differences in the molecular epidemiology of tuberculosis

OBJECTIVES

Compare the molecular epidemiology of tuberculosis (TB) between two large Canadian provinces-Ontario and British Columbia (BC)-to identify genotypic clusters within and across both provinces, allowing for an improved understanding of genotype data and providing context to more accurately identify clusters representing local transmission.

DESIGN

We compared 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping for 3,314 Ontario and 1,602 BC clinical Mycobacterium tuberculosis isolates collected from 2008 through 2014. Laboratory data for each isolate was linked to case-level records to obtain clinical and demographic data.

RESULTS

The demographic characteristics of persons with TB varied between provinces, most notably in the proportion of persons born outside Canada, which was reflected in the large number of unique genotypes (n = 3,461). The proportion of clustered isolates was significantly higher in BC. Substantial clustering amongst non-Lineage 4 TB strains was observed within and across the provinces. Only two large clusters (≥10 cases/cluster) representing within province transmission had interprovincial genotype matches.

CONCLUSION

We recommend expanding analysis of shared genotypes to include neighbouring jurisdictions, and implementing whole genome sequencing to improve identification of TB transmission, recognize outbreaks, and monitor changing trends in TB epidemiology.

Whole-Genome Sequencing for Drug Resistance Profile Prediction in Mycobacterium tuberculosis

Whole-genome sequencing allows rapid detection of drug-resistant Mycobacterium tuberculosis isolates. However, the availability of high-quality data linking quantitative phenotypic drug susceptibility testing (DST) and genomic data have thus far been limited. We determined drug resistance profiles of 176 genetically diverse clinical M. tuberculosis isolates from the 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 to 99% of cases, depending on the drug. Our data suggest that the established critical concentration (5 mg/liter) for ethambutol resistance (MGIT 960 system) is too high and misclassifies strains as susceptible, unlike 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.

Transmission dynamics study of tuberculosis isolates with whole genome sequencing in southern Sweden

Epidemiological contact tracing complemented with genotyping of clinical Mycobacterium tuberculosis isolates is important for understanding disease transmission. In Sweden, tuberculosis (TB) is mostly reported in migrant and homeless where epidemiologic contact tracing could pose a problem. This study compared epidemiologic linking with genotyping in a low burden country. Mycobacterium tuberculosis isolates (n = 93) collected at Scania University Hospital in Southern Sweden were analysed with the standard genotyping method mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) and the results were compared with whole genome sequencing (WGS). Using a maximum of twelve single nucleotide polymorphisms (SNPs) as the upper threshold of genomic relatedness noted among hosts, we identified 18 clusters with WGS comprising 52 patients with overall pairwise genetic maximum distances ranging from zero to nine SNPs. MIRU-VNTR and WGS clustered the same isolates, although the distribution differed depending on MIRU-VNTR limitations. Both genotyping techniques identified clusters where epidemiologic linking was insufficient, although WGS had higher correlation with epidemiologic data. To summarize, WGS provided better resolution of transmission than MIRU-VNTR in a setting with low TB incidence. WGS predicted epidemiologic links better which could consolidate and correct the epidemiologically linked cases, avoiding thus false clustering.

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%).