A Predominant Variable-Number Tandem-Repeat Cluster of Mycobacterium tuberculosis Isolates among Asylum Seekers in the Netherlands and Denmark, Deciphered by Whole-Genome Sequencing

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.

Pathways of TB Transmission in Children-A Systematic Review of Molecular Epidemiological Studies

The widespread paradigm that younger children usually do not transmit M. tuberculosis complex (Mtbc) to their contacts has not yet been proven by genotypically confirmed transmissions. Therefore, we undertook a systematic review of molecular-epidemiological studies to investigate documented source and secondary TB (tuberculosis) cases among children. We searched the literature published before August 2022 using PubMed, Cochrane, and Google Scholar databases. PRISMA statement was used for systematic review. Of 312 records retrieved, 39 studies including children aged below 15 years offered epidemiological links between cluster members. In the 39 studies from 16 countries, 225 children were reported as cluster members of whom the overwhelming majority were infected by adults. Only 3 children-of those were 2 children aged below 10-were reported to be the definite source cases of 11 other children and 1 adult with genotypically matched Mtbc isolates. To date, molecular-epidemiological studies involving children with verified transmission links are scarce. As far as the heterogeneity of the studies we identified allows, we could conclude that the results confirm the paradigm that children aged below 10 hardly ever transmit Mtbc to others. The true extent of TB transmission through children may, however, be underestimated by those selected studies.

Optimization and Validation of Candida auris Short Tandem Repeat Analysis

Candida auris is an easily transmissible yeast with resistance to different antifungal compounds. Outbreaks of C. auris are mostly observed in intensive care units. To take adequate measures during an outbreak, it is essential to understand the transmission route, which requires isolate genotyping. In 2019, a short tandem repeat (STR) genotyping analysis was developed for C. auris. To determine the discriminatory power of this method, we performed STR analysis of 171 isolates with known whole-genome sequencing (WGS) data using Illumina reads, and we compared their resolutions. We found that STR analysis separated the 171 isolates into four clades (clades I to IV), as was also seen with WGS analysis. Then, to improve the separation of isolates in clade IV, the STR assay was optimized by the addition of 2 STR markers. With this improved STR assay, a total of 32 different genotypes were identified, while all isolates with differences of >50 single-nucleotide polymorphisms (SNPs) were separated by at least 1 STR marker. Altogether, we optimized and validated the C. auris STR panel for clades I to IV and established its discriminatory power, compared to WGS SNP analysis using Illumina reads. IMPORTANCE The emerging fungal pathogen Candida auris poses a threat to public health, mainly causing outbreaks in intensive care units. Genotyping is essential for investigating potential outbreaks and preventing further spread. Previously, we developed a STR genotyping scheme for rapid and high-resolution genotyping, and WGS SNP outcomes for some isolates were compared to STR data. Here, we compared WGS SNP and STR outcomes for a larger sample cohort. Also, we optimized the resolution of this typing scheme with the addition of 2 STR markers. Altogether, we validated and optimized this rapid, reliable, and high-resolution typing scheme for C. auris.

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.

An implementation science approach to evaluating pathogen whole genome sequencing in public health

BACKGROUND

Pathogen whole genome sequencing (WGS) is being incorporated into public health surveillance and disease control systems worldwide and has the potential to make significant contributions to infectious disease surveillance, outbreak investigation and infection prevention and control. However, to date, there are limited data regarding (i) the optimal models for integration of genomic data into epidemiological investigations and (ii) how to quantify and evaluate public health impacts resulting from genomic epidemiological investigations.

METHODS

We developed the Pathogen Genomics in Public HeAlth Surveillance Evaluation (PG-PHASE) Framework to guide examination of the use of WGS in public health surveillance and disease control. We illustrate the use of this framework with three pathogens as case studies: Listeria monocytogenes, Mycobacterium tuberculosis and SARS-CoV-2.

RESULTS

The framework utilises an adaptable whole-of-system approach towards understanding how interconnected elements in the public health application of pathogen genomics contribute to public health processes and outcomes. The three phases of the PG-PHASE Framework are designed to support understanding of WGS laboratory processes, analysis, reporting and data sharing, and how genomic data are utilised in public health practice across all stages, from the decision to send an isolate or sample for sequencing to the use of sequence data in public health surveillance, investigation and decision-making. Importantly, the phases can be used separately or in conjunction, depending on the need of the evaluator. Subsequent to conducting evaluation underpinned by the framework, avenues may be developed for strategic investment or interventions to improve utilisation of whole genome sequencing.

CONCLUSIONS

Comprehensive evaluation is critical to support health departments, public health laboratories and other stakeholders to successfully incorporate microbial genomics into public health practice. The PG-PHASE Framework aims to assist public health laboratories, health departments and authorities who are either considering transitioning to whole genome sequencing or intending to assess the integration of WGS in public health practice, including the capacity to detect and respond to outbreaks and associated costs, challenges and facilitators in the utilisation of microbial genomics and public health impacts.

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.

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.

Towards elimination of childhood and adolescent tuberculosis in the Netherlands: an epidemiological time-series analysis of national surveillance data

BACKGROUND

Tuberculosis (TB) in children and adolescents is a sentinel event for ongoing transmission. In the Netherlands, epidemiological characteristics of childhood and adolescent TB have not been fully evaluated. Therefore, we aimed to assess TB epidemiology within this population to provide guidance for TB elimination.

METHODS

A retrospective time-series analysis using national surveillance data from 1993-2018 was performed in children (aged <15 years) and adolescents (aged 15-19 years) with TB. Poisson regression models offset with log-population size were used to estimate notification rates and rate ratios. Trends in notification rates were estimated using average annual percentage changes (AAPC) based on the segmented linear regression analysis.

RESULTS

Among 3899 children and adolescents with TB notified during 1993-2018, 2418 (62%) were foreign-born (725 (41.3%) out of 1755 children and 1693 (78.9%) out of 2144 adolescents). The overall notification rate in children was 2.3 per 100 000 person-years, declining steadily during the study period (AAPC -10.9%, 95% CI -12.6--9.1). In adolescents, the overall notification rate was 8.4 per 100 000 person-years, strongly increasing during 1993-2001 and 2012-2018. Compared to Dutch-born children and adolescents, substantially higher notification rates were observed among African-born children and adolescents (116.8 and 316.6 per 100 000 person-years, respectively). Additionally, an increasing trend was observed in African-born adolescents (AAPC 18.5%, 95% CI 11.9-25.5). Among the foreign-born population, those from countries in the horn of Africa contributed most to the TB caseload.

CONCLUSION

TB notification rate among children was low and constantly declining across different demographic groups. However, heterogeneities were shown in adolescents, with an increasing trend in the foreign-born, particularly those from Africa.

Climate change impact on migration, travel, travel destinations and the tourism industry

Background: Climate change is not only increasing ambient temperature but also accelerating the frequency, duration and intensity of extreme weather and climate events, such as heavy precipitation and droughts, and causing sea level rise, which can lead to population displacement. Climate change-related reductions in land productivity and habitability and in food and water security can also interact with demographic, economic and social factors to increase migration. In addition to migration, climate change has also implications for travel and the risk of disease. This article discusses the impact of climate change on migration and travel with implications for public health practice. Methods: Literature review. Results: Migrants may be at increased risk of communicable and non-communicable diseases, due to factors in their country of origin and their country of destination or conditions that they experience during migration. Although migration has not been a significant driver of communicable disease outbreaks to date, public health authorities need to ensure that effective screening and vaccination programmes for priority communicable diseases are in place.Population growth coupled with socio-economic development is increasing travel and tourism, and advances in technology have increased global connectivity and reduced the time required to cover long distances. At the same time, as a result of climate change, many temperate regions, including high-income countries, are now suitable for vector-borne disease transmission. This is providing opportunities for importation of vectors and pathogens from endemic areas that can lead to cases or outbreaks of communicable diseases with which health professionals may be unfamiliar. Conclusion: Health systems need to be prepared for the potential population health consequences of migration, travel and tourism and the impact of climate change on these. Integrated surveillance, early detection of cases and other public health interventions are critical to protect population health and prevent and control communicabledisease outbreaks.

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.

Prevalence and Characterization of Heterogeneous Variable-Number Tandem-Repeat Clusters Comprising Drug-Susceptible and/or Variable Resistant Mycobacterium tuberculosis Complex Isolates in the Netherlands from 2004 to 2016

The variable-number tandem-repeat (VNTR) typing method is used to study tuberculosis (TB) transmission. Clustering of Mycobacterium tuberculosis isolates with identical VNTR patterns is assumed to reflect recent transmission. Hence, clusters are thought to be homogeneous regarding antibiotic resistance. In practice, however, heterogeneous clusters are also identified. This study investigates the prevalence and characteristics of heterogeneous VNTR clusters and assesses whether isolates in these clusters remain clustered when subjected to whole-genome sequencing (WGS). In the period from 2004 to 2016, 9,072 isolates were included. Demographic and epidemiological linkage data were obtained from the Netherlands Tuberculosis Register. VNTR clusters were defined as homogeneous when isolates shared identical resistance profiles or as heterogeneous if both susceptible and (variable) resistant isolates were found. Multivariate logistic regression analysis was performed to identify factors associated with heterogeneous clustering. Isolates from 2016 were subjected to WGS, and a genetic distance of 12 single nucleotide polymorphisms (SNPs) was used as the cutoff for WGS clustering. In total, 4,661/9,072 (51%) isolates were clustered into 985 different VNTR clusters, of which 217 (22%) were heterogeneous. Patient characteristics associated with heterogeneous clustering were non-Dutch ethnicity (odds ratio [OR], 1.46 [95% confidence interval {CI}, 1.22 to 1.75]), asylum seeker (OR, 1.51 [95% CI, 1.24 to 1.85]), extrapulmonary TB (OR, 1.26 [95% CI, 1.09 to 1.46]), previous TB diagnosis (OR, 1.38 [95% CI, 1.04 to 1.82]), and not being a contact of a TB patient (OR, 1.35 [95% CI, 1.08 to 1.69]). With WGS, 34% of heterogeneous and 78% of homogeneous isolates from 2016 remained clustered. Heterogeneous VNTR clusters are common but seem to be explained by a substantial degree of false clustering by VNTR typing compared to WGS.

Genome-based transmission modelling separates imported tuberculosis from recent transmission within an immigrant population

In many countries the incidence of tuberculosis (TB) is low and is largely shaped by immigrant populations from high-burden countries. This is the case in Norway, where more than 80 % of TB cases are found among immigrants from high-incidence countries. A variable latent period, low rates of evolution and structured social networks make separating import from within-border transmission a major conundrum to TB control efforts in many low-incidence countries. Clinical Mycobacterium tuberculosis isolates belonging to an unusually large genotype cluster associated with people born in the Horn of Africa have been identified in Norway over the last two decades. We modelled transmission based on whole-genome sequence data to estimate infection times for individual patients. By contrasting these estimates with time of arrival in Norway, we estimate on a case-by-case basis whether patients were likely to have been infected before or after arrival. Independent import was responsible for the majority of cases, but we estimate that about one-quarter of the patients had contracted TB in Norway. This study illuminates the transmission dynamics within an immigrant community. Our approach is broadly applicable to many settings where TB control programmes can benefit from understanding when and where patients acquired TB.

Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: A population-based study

BACKGROUND

Patients with Mycobacterium tuberculosis isolates sharing identical DNA fingerprint patterns can be epidemiologically linked. However, municipal health services in the Netherlands are able to confirm an epidemiological link in only around 23% of the patients with isolates clustered by the conventional variable number of tandem repeat (VNTR) genotyping. This research aims to investigate whether whole genome sequencing (WGS) is a more reliable predictor of epidemiological links between tuberculosis patients than VNTR genotyping.

METHODS

VNTR genotyping and WGS were performed in parallel on all Mycobacterium tuberculosis complex isolates received at the Netherlands National Institute for Public Health and the Environment in 2016. Isolates were clustered by VNTR when they shared identical 24-loci VNTR patterns; isolates were assigned to a WGS cluster when the pair-wise genetic distance was ≤ 12 single nucleotide polymorphisms (SNPs). Cluster investigation was performed by municipal health services on all isolates clustered by VNTR in 2016. The proportion of epidemiological links identified among patients clustered by either method was calculated.

RESULTS

In total, 535 isolates were genotyped, of which 25% (134/535) were clustered by VNTR and 14% (76/535) by WGS; the concordance between both typing methods was 86%. The proportion of epidemiological links among WGS clustered cases (57%) was twice as common than among VNTR clustered cases (31%).

CONCLUSION

When WGS was applied, the number of clustered isolates was halved, while all epidemiologically linked cases remained clustered. WGS is therefore a more reliable tool to predict epidemiological links between tuberculosis cases than VNTR genotyping and will allow more efficient transmission tracing, as epidemiological investigations based on false clustering can be avoided.