Whole Genome Sequencing and Spatial Analysis Identifies Recent Tuberculosis Transmission Hotspots in Ghana

Genomic and Spatial Analysis on the Recent Transmission of Mycobacterium tuberculosis in Eastern China: A 10-Year Retrospective Population-Based Study

Purpose

Understanding the mode of Mycobacterium tuberculosis (M. tuberculosis) transmission is crucial for disease prevention and control. Compared to traditional genotyping methods, whole genome sequencing (WGS) provides higher resolution and comprehensive genetic information, enabling the tracing of infection sources and determining of transmission routes to resolve extensive tuberculosis (TB) outbreaks. We conducted a ten-year study on the transmission of M. tuberculosis in a population in eastern China.

Patients and Methods

We selected Lianyungang, an eastern city in China, as the study site. Patients diagnosed with active pulmonary TB from 2011 to 2020 were enrolled as the study subjects. We isolated and sequenced 2252 M. tuberculosis. Strains with pairwise genetic distances of less than 12 single nucleotide polymorphisms were defined as genomic clusters and which were considered recent transmissions. Kernel density estimation and K-function analysis were applied to explore the spatial distribution of recently transmitted strains.

Results

After excluding non-tuberculous mycobacteria and duplicated samples, 2114 strains were included in the final analysis. These strains comprised lineage 2 (1593, 75.35%) and 4 (521, 24.65%). There were 672 clustered strains, with a recent transmission rate of 31.79%. The logistic regression model showed that the risk of recent transmission was high in students [adjusted odds ratio (aOR): 2.68, 95% confidence interval (CI): 1.63-4.49, P<0.001] and people infected with L2.2.1 strains (aOR: 1.59, 95% CI: 1.20-2.12). Higher spatial aggregation of TB transmission has been concentrated in Haizhou, Donghai, and Guanyun for the past 10 years. Three outbreaks affecting 46 patients were spatially spaced, with 11 to 23 persons each. Different groups exhibited varying geographic distances between the initial and later cases.

Conclusion

There are areas with a high risk of transmission for M. tuberculosis in the research site, and the risk varies among different populations. Accurate prevention strategies targeted at specific regions and key populations can help curb the prevalence of TB.

Social-behavioral insights in understanding tuberculosis transmission pattern during the COVID-19 pandemic period in Kuala Lumpur, Malaysia: The MyTBNet study protocol

Many countries have reported increase of TB incidence during the COVID-19 pandemic period, which demands dire attention as it may threaten global effort to end TB transmission. Services, are among many others, were disrupted by the COVID-19 pandemic during the years 2020 and 2021; but its impact on the TB transmission is not well understood. This retrospective population-based molecular and epidemiological cohort study aims to determine the pattern of TB transmission in Kuala Lumpur (an area with high population density, moderate TB burden and high rates of COVID-19 cases) for the cohort of Pulmonary TB (PTB) cases notified from 2020 until 2021 and factors associated with clustering or clear epidemiologic linkage. This study will be carried out from 2022 until 2024. The study will utilise comparative phylogenetic analysis to determine the degree of relatedness between different isolates, based on the genomes similarities, and overlay this with epidemiological, clinical and social network data to enhance understanding of the social-behavioural dynamics of TB transmission. Mycobacterium tuberculosis complex (MTBC) cultures will be genotyped using Mycobacterial Interspersed Repetitive Unit Variable Number Tandem Repeats (MIRU-VNTR) and whole-genome sequence (WGS) for MTBC cluster isolates. Epidemiologic and genomic data will be overlaid on a social network constructed by means of interviews with patients, by using Social Network Analysis questionnaire, to determine the origins and transmission dynamics of the outbreak. The finding of this study would aid in the identification of TB transmission events, facilitating active case finding, TB screening, TB contact tracing, and the mapping of social contacts during critical period. This will contribute to building an effective preventive and preparedness strategy to interrupt TB transmission in Malaysia, tailored to the characteristics of the local population.

Define SNP thresholds for delineation of tuberculosis transmissions using whole-genome sequencing

For facilitating tuberculosis (TB) control, we used a whole-genome sequencing (WGS)-based approach to delineate transmission networks in a country with an intermediate burden of TB. A cluster was defined as Mycobacterium tuberculosis isolates with identical genotypes, and an outbreak was defined as clustered cases with epidemiological links (epi-links). To refine a cluster predefined using space oligonucleotide typing and mycobacterial interspersed repetitive unit variable tandem repeat typing, we analyzed one pansusceptible TB (C1) and three multidrug-resistant (MDR)-TB (C2-C4) clusters from different scenarios. Pansusceptible TB cluster (C1) consisting of 28 cases had ≤5 single nucleotide polymorphisms (SNPs) difference between their isolates. C1 was a definite outbreak, with cases attending the same junior high school in 2012. Three MDR-TB clusters (C2-C4) with distinct genotypes were identified, each consisting of 12-22 cases. Some of the cases had either ≤5 or ≤15 SNPs difference with clear or probable epi-links. Of note, even though WGS could effectively assist TB contact tracing, we still observed missing epi-links in some cases within the same cluster. Our results showed that thresholds of ≤5 and ≤15 SNPs difference between isolates were used to categorize definite and probable TB transmission, respectively. Furthermore, a higher SNP threshold might be required to define an MDR-TB outbreak. WGS still needs to be combined with classical epidemiological methods for improving outbreak investigations. Importantly, different SNP thresholds have to be applied to define outbreaks.

IMPORTANCE

TB is a chronic disease. Depending on host factors and TB burden, clusters of cases may continue to increase for several years. Conventional genotyping methods overestimate TB transmission, hampering precise detection of outbreaks and comprehensive surveillance. WGS can be used to obtain SNP information of M. tuberculosis to improve discriminative limitations of conventional methods and to strengthen delineation of transmission networks. It is important to define the country-specific SNP thresholds for investigation of transmission. This study demonstrated the use of thresholds of ≤5 and ≤15 SNPs difference between isolates to categorize definite and probable transmission, respectively. Different SNP thresholds should be applied while a higher cutoff was required to define an MDR-TB outbreak. The utilization of SNP thresholds proves to be crucial for guiding public health interventions, eliminating the need for unnecessary public health actions, and potentially uncovering undisclosed TB transmissions.

Opinion review of drug resistant tuberculosis in West Africa: tackling the challenges for effective control

Drug-resistant (DR) tuberculosis (TB) is a major public health concern globally, complicating TB control and management efforts. West Africa has historically faced difficulty in combating DR-TB due to limited diagnostic skills, insufficient access to excellent healthcare, and ineffective healthcare systems. This has aided in the emergence and dissemination of DR Mycobacterium tuberculosis complex (MTBC) strains in the region. In the past, DR-TB patients faced insufficient resources, fragmented efforts, and suboptimal treatment outcomes. However, current efforts to combat DR-TB in the region are promising. These efforts include strengthening diagnostic capacities, improving access to quality healthcare services, and implementing evidence-based treatment regimens for DR-TB. Additionally, many West African National TB control programs are collaborating with international partners to scale up laboratory infrastructure, enhance surveillance systems, and promote infection control measures. Moreso, novel TB drugs and regimens, such as bedaquiline and delamanid, are being introduced to improve treatment outcomes for DR-TB cases. Despite these obstacles, there is optimism for the future of DR-TB control in West Africa. Investments are being made to improve healthcare systems, expand laboratory capacity, and support TB research and innovation. West African institutions are now supporting knowledge sharing, capacity building, and resource mobilization through collaborative initiatives such as the West African Network for TB, AIDS, and Malaria (WANETAM), the West African Health Organization (WAHO), and other regional or global partners. These efforts hold promise for improved diagnostics, optimized treatment regimens, and provide better patient outcomes in the future where drug-resistant TB in WA can be effectively controlled, reducing the burden of the disease, and improving the health outcomes of affected individuals.

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.

Botswana tuberculosis (TB) stakeholders broadly support scaling up next-generation whole genome sequencing: Ethical and practical considerations for Botswana and global health

Global health agencies are increasingly promoting the scale-up of next-generation whole genome sequencing (NG-WGS) of pathogens into infectious disease control programs, including for tuberculosis (TB). However, little is known about how stakeholders in low-to-middle income countries (LMICs) understand the ethics, benefits, and risks of these proposals. We conducted a qualitative study in Greater Gaborone, Botswana to learn how TB stakeholders there viewed a potential scale-up of NG-WGS into Botswana's TB program. We conducted 30 interviews and four deliberative dialogues with TB stakeholders based in Greater Gaborone, the country's largest city and capital. We created and showed participants an animated video series about a fictional family that experienced TB diagnosis, treatment, contact tracing, and data uses that were informed by NG-WGS. We analyzed transcripts using reflexive thematic analysis. We found broad support for the scale-up of TB NG-WGS in Botswana, owing to perceived benefits. Support was qualified with statements about ensuring adequate planning, resource-allocation, community and stakeholder engagement, capacity-building, and assessing ethical norms around publishing data. Our results suggest that scaling up NG-WGS for TB in Botswana would be supported by stakeholders there, contingent upon the government and other entities adequately investing in the initiative. These findings are relevant to other LMICs considering scale-ups of NG-WGS and related technologies for infectious diseases and suggest the need for sustained research into the acceptability of pathogen sequencing in other contexts.

Use of High-Resolution Geospatial and Genomic Data to Characterize Recent Tuberculosis Transmission, Botswana

Combining genomic and geospatial data can be useful for understanding Mycobacterium tuberculosis transmission in high-burden tuberculosis (TB) settings. We performed whole-genome sequencing on M. tuberculosis DNA extracted from sputum cultures from a population-based TB study conducted in Gaborone, Botswana, during 2012-2016. We determined spatial distribution of cases on the basis of shared genotypes among isolates. We considered clusters of isolates with ≤5 single-nucleotide polymorphisms identified by whole-genome sequencing to indicate recent transmission and clusters of ≥10 persons to be outbreaks. We obtained both molecular and geospatial data for 946/1,449 (65%) participants with culture-confirmed TB; 62 persons belonged to 5 outbreaks of 10-19 persons each. We detected geospatial clustering in just 2 of those 5 outbreaks, suggesting heterogeneous spatial patterns. Our findings indicate that targeted interventions applied in smaller geographic areas of high-burden TB identified using integrated genomic and geospatial data might help interrupt TB transmission during outbreaks.

An overview of tuberculosis outbreaks reported in the years 2011-2020

BACKGROUND

In many countries tuberculosis (TB) remains a highly prevalent disease and a major contributor to infectious disease mortality. The fight against TB requires surveillance of the population of strains circulating worldwide and the analysis of the prevalence of certain strains in populations. Nowadays, whole genome sequencing (WGS) allows for accurate tracking of TB transmission. Currently, there is a lack of a comprehensive summary of the characteristics of TB outbreaks.

METHODS

We systematically analyzed studies reporting TB outbreaks worldwide, monitored through WGS of Mycobacterium tuberculosis. We 1) mapped the reported outbreaks from 2011- 2020, 2) estimated the average size of the outbreaks, 3) indicated genetic lineages causing the outbreaks, and 4) determined drug-resistance patterns of M. tuberculosis strains involved in the outbreaks.

RESULTS

Most data originated from Europe, Asia, and North America. We found that TB outbreaks were reported throughout the globe, on all continents, and in countries with both high and low incidences. The detected outbreaks contained a median of five M. tuberculosis isolates. Most strains causing the outbreaks belonged to lineage four, more rarely to lineage two. Reported outbreak isolates were often drug resistant.

CONCLUSIONS

We conclude that more WGS surveillance of M. tuberculosis outbreaks is needed. Globally standardized procedures might improve the control of M. tuberculosis infections.

Spatial Distribution of Drug-Resistant Mycobacterium tuberculosis Infections in Rural Eastern Cape Province of South Africa

Tuberculosis (TB), an infectious airborne disease caused by Mycobacterium tuberculosis (Mtb), is a serious public health threat reported as the leading cause of morbidity and mortality worldwide. South Africa is a high-TB-burden country with TB being the highest infectious disease killer. This study investigated the distribution of Mtb mutations and spoligotypes in rural Eastern Cape Province. The Mtb isolates included were 1157 from DR-TB patients and analysed by LPA followed by spoligotyping of 441 isolates. The distribution of mutations and spoligotypes was done by spatial analysis. The rpoB gene had the highest number of mutations. The distribution of rpoB and katG mutations was more prevalent in four healthcare facilities, inhA mutations were more prevalent in three healthcare facilities, and heteroresistant isolates were more prevalent in five healthcare facilities. The Mtb was genetically diverse with Beijing more prevalent and largely distributed. Spatial analysis and mapping of gene mutations and spoligotypes revealed a better picture of distribution.

Quantifying Mycobacterium tuberculosis Transmission Dynamics Across Global Settings: A Systematic Analysis

The degree to which individual heterogeneity in the production of secondary cases ("superspreading") affects tuberculosis (TB) transmission has not been systematically studied. We searched for population-based or surveillance studies in which whole genome sequencing was used to estimate TB transmission and in which the size distributions of putative TB transmission clusters were enumerated. We fitted cluster-size-distribution data to a negative binomial branching process model to jointly infer the transmission parameters $R$ (the reproduction number) and the dispersion parameter, $k$, which quantifies the propensity of superspreading in a population (generally, lower values of $k$ ($<1.0$) suggest increased heterogeneity). Of 4,796 citations identified in our initial search, 9 studies from 8 global settings met the inclusion criteria (n = 5 studies of all TB; n = 4 studies of drug-resistant TB). Estimated $R$ values (range, 0.10-0.73) were below 1.0, consistent with declining epidemics in the included settings; estimated $k$ values were well below 1.0 (range, 0.02-0.48), indicating the presence of substantial individual-level heterogeneity in transmission across all settings. We estimated that a minority of cases (range, 2%-31%) drive the majority (80%) of ongoing TB transmission at the population level. Identifying sources of heterogeneity and accounting for them in TB control may have a considerable impact on mitigating TB transmission.

Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis

Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.

Diagnosis of tuberculosis among COVID-19 suspected cases in Ghana

Background

Tuberculosis (TB) and COVID-19 pandemics are both diseases of public health threat globally. Both diseases are caused by pathogens that infect mainly the respiratory system, and are involved in airborne transmission; they also share some clinical signs and symptoms. We, therefore, took advantage of collected sputum samples at the early stage of COVID-19 outbreak in Ghana to conduct differential diagnoses of long-standing endemic respiratory illness, particularly tuberculosis.

Methodology

Sputum samples collected through the enhanced national surveys from suspected COVID-19 patients and contact tracing cases were analyzed for TB. The sputum samples were processed using Cepheid’s GeneXpert MTB/RIF assay in pools of 4 samples to determine the presence of Mycobacterium tuberculosis complex. Positive pools were then decoupled and analyzed individually. Details of positive TB samples were forwarded to the NTP for appropriate case management.

Results

Seven-hundred and seventy-four sputum samples were analyzed for Mycobacterium tuberculosis in both suspected COVID-19 cases (679/774, 87.7%) and their contacts (95/774, 12.3%). A total of 111 (14.3%) were diagnosed with SARS CoV-2 infection and six (0.8%) out of the 774 individuals tested positive for pulmonary tuberculosis: five (83.3%) males and one female (16.7%). Drug susceptibility analysis identified 1 (16.7%) rifampicin-resistant tuberculosis case. Out of the six TB positive cases, 2 (33.3%) tested positive for COVID-19 indicating a coinfection. Stratifying by demography, three out of the six (50%) were from the Ayawaso West District. All positive cases received appropriate treatment at the respective sub-district according to the national guidelines.

Conclusion

Our findings highlight the need for differential diagnosis among COVID-19 suspected cases and regular active TB surveillance in TB endemic settings.

Whole-Genome Sequencing Reveals Recent Transmission of Multidrug-Resistant Mycobacterium tuberculosis CAS1-Kili Strains in Lusaka, Zambia

Globally, tuberculosis (TB) is a major cause of death due to antimicrobial resistance. Mycobacterium tuberculosis CAS1-Kili strains that belong to lineage 3 (Central Asian Strain, CAS) were previously implicated in the spread of multidrug-resistant (MDR)-TB in Lusaka, Zambia. Thus, we investigated recent transmission of those strains by whole-genome sequencing (WGS) with Illumina MiSeq platform. Twelve MDR CAS1-Kili isolates clustered by traditional methods (MIRU-VNTR and spoligotyping) were used. A total of 92% (11/12) of isolates belonged to a cluster (≤12 SNPs) while 50% (6/12) were involved in recent transmission events, as they differed by ≤5 SNPs. All the isolates had KatG Ser315Thr (isoniazid resistance), EmbB Met306 substitutions (ethambutol resistance) and several kinds of rpoB mutations (rifampicin resistance). WGS also revealed compensatory mutations including a novel deletion in embA regulatory region (−35A > del). Several strains shared the same combinations of drug-resistance-associated mutations indicating transmission of MDR strains. Zambian strains belonged to the same clade as Tanzanian, Malawian and European strains, although most of those were pan-drug-susceptible. Hence, complimentary use of WGS to traditional epidemiological methods provides an in-depth insight on transmission and drug resistance patterns which can guide targeted control measures to stop the spread of MDR-TB.

A Saudi Arabian Public Health Perspective of Tuberculosis

Tuberculosis is a global health challenge due to its spreading potential. The Kingdom of Saudi Arabia (KSA) faces a challenge in the spread of tuberculosis from migrant workers, but the foremost threat is the huge number of pilgrims who travel to visit sacred sites of the Islamic world located in the holy cities of Makkah and Al Madina. Pilgrims visit throughout the year but especially in the months of Ramadan and Zul-Hijah. The rise of resistance in Mycobacterium tuberculosis is an established global phenomenon that makes such large congregations likely hotspots in the dissemination and spread of disease at a global level. Although very stringent and effective measures exist, the threat remains due to the ever-changing dynamics of this highly pathogenic disease. This overview primarily highlights the current public health challenges posed by this disease to the Saudi health system, which needs to be highlighted not only to the concerned authorities of KSA, but also to the concerned global quarters since the pilgrims and migrants come from all parts of the world with a majority coming from high tuberculosis-burdened countries.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.