Abstract:Whole genome sequencing (WGS) has become the main tool for studying the transmission of Mycobacterium tuberculosis complex (MTBC) strains; however, the clonal expansion of one strain often limits its application in local MTBC outbreaks. The use of an alternative reference genome and the inclusion of repetitive regions in the analysis could potentially increase the resolution, but the added value has not yet been defined. Here, we leveraged short and long WGS read data of a previously reported MTBC outbreak in the Colombian Amazon Region to analyze possible transmission chains among 74 patients in the indigenous setting of Puerto Nariño (March to October 2016). In total, 90.5% (67/74) of the patients were infected with one distinct MTBC strain belonging to lineage 4.3.3. Employing a reference genome from an outbreak strain and highly confident single nucleotide polymorphisms (SNPs) in repetitive genomic regions, e.g., the proline-glutamic acid/proline-proline-glutamic-acid (PE/PPE) gene family, increased the phylogenetic resolution compared to a classical H37Rv reference mapping approach. Specifically, the number of differentiating SNPs increased from 890 to 1,094, which resulted in a more granular transmission network as judged by an increasing number of individual nodes in a maximum parsimony tree, i.e., 5 versus 9 nodes. We also found in 29.9% (20/67) of the outbreak isolates, heterogenous alleles at phylogenetically informative sites, suggesting that these patients are infected with more than one clone. In conclusion, customized SNP calling thresholds and employment of a local reference genome for a mapping approach can improve the phylogenetic resolution in highly clonal MTBC populations and help elucidate within-host MTBC diversity. IMPORTANCE The Colombian Amazon around Puerto Nariño has a high tuberculosis burden with a prevalence of 1,267/100,000 people in 2016. Recently, an outbreak of Mycobacterium tuberculosis complex (MTBC) bacteria among the indigenous populations was identified with classical MTBC genotyping methods. Here, we employed a whole-genome sequencing-based outbreak investigation in order to improve the phylogenetic resolution and gain new insights into the transmission dynamics in this remote Colombian Amazon Region. The inclusion of well-supported single nucleotide polymorphisms in repetitive regions and a de novo-assembled local reference genome provided a more granular picture of the circulating outbreak strain and revealed new transmission chains. Multiple patients from different settlements were possibly infected with at least two different clones in this high-incidence setting. Thus, our results have the potential to improve molecular surveillance studies in other high-burden settings, especially regions with few clonal multidrug-resistant (MDR) MTBC lineages/clades.

Population structure of Mycobacterium tuberculosis in El Oro: a first insight into Ecuador-Peru tuberculosis transmission.

A first insight into tuberculosis transmission at the border of Ecuador and Colombia: a retrospective study of the population structure of Mycobacterium tuberculosis in Esmeraldas province

Molecular epidemiology of Mycobacterium tuberculosis in Ecuador: Recent advances and future challenges

Longitudinal Molecular Epidemiological Study of Multidrug Resistant Tb and Transition to Extensively Drug Resistance in Eastern Rual China

Prevalence, drug resistance, and genotypic diversity of the RDRio subfamily of Mycobacterium tuberculosis in Ecuador: a retrospective analysis for years 2012–2016

Who Transmits Tuberculosis to Whom: A Cross-Sectional Analysis of a Cohort Study in Lima, Peru

Whole genome analysis of extensively drug resistant Mycobacterium tuberculosis strains in Peru

Population-based sequencing of Mycobacterium tuberculosis reveals how current population dynamics are shaped by past epidemics

Transmission Dynamics of a Mycobacterium tuberculosis Complex Outbreak in an Indigenous Population in the Colombian Amazon Region

Genetic diversity and drug susceptibility of Mycobacterium tuberculosis in a city with a high prevalence of drug resistant tuberculosis from Southeast of Mexico

Pan-genome association study of Mycobacterium tuberculosis lineage-4 revealed specific genes related to the high and low prevalence of the disease in patients from the North-Eastern area of Medellín, Colombia

Comparative whole-genome sequence analysis of Mycobacterium tuberculosis isolated from pulmonary tuberculosis and tuberculous lymphadenitis patients in Northwest Ethiopia

A Genome-Focused Investigation Reveals the Emergence of a Mycobacterium tuberculosis Strain Related to Multidrug-Resistant Tuberculosis in the Amazon Region of Brazil

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread

Whole-genome sequencing-based genetic diversity, transmission dynamics, and drug-resistant mutations in Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

Whole Genome Sequencing of Mycobacterium tuberculosis under routine conditions in a high-burden area of multidrug-resistant tuberculosis in Peru

Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay

Whole genome sequencing identifies bacterial factors affecting transmission of multidrug-resistant tuberculosis in a high-prevalence setting

Evolutionary Trajectories and Transmission Dynamics of Multidrug-Resistant <i>mycobacterium Tuberculosis</i> in Tibet, China

Assessing a transmission network of Mycobacterium tuberculosis in an African city using single nucleotide polymorphism threshold analysis

Transmission Ofmycobacterium Tuberculosisin China: A Population-Based Molecular Epidemiologic Study