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.

Phylodynamic analysis reveals disparate transmission dynamics of Mycobacterium tuberculosis-complex lineages in Botswana

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

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Differential rates of Mycobacterium tuberculosis transmission associate with host–pathogen sympatry

Transmission dynamics and phylogeography of Mycobacterium tuberculosis in China based on whole-genome phylogenetic analysis

Genetic Composition and Evolution of the Prevalent Mycobacterium Tuberculosis Lineages 2 and 4 in the Chinese and Zhejiang Province Populations

Local Adaptation in Populations of Mycobacterium Tuberculosis Endemic to the Indian Ocean Rim

Mycobacterium Tuberculosis Lineage 4 Comprises Globally Distributed and Geographically Restricted Sublineages

New Mycobacterium tuberculosis Beijing clonal complexes in China revealed by phylogenetic and Bayesian population structure analyses of 24-loci MIRU-VNTRs

Bayesian Population Structure Analysis Reveals Presence of Phylogeographically Specific Sublineages Within Previously Ill-Defined T Group of Mycobacterium Tuberculosis.

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

Back-to-Africa introductions of Mycobacterium tuberculosis as the main cause of tuberculosis in Dar es Salaam, Tanzania

Evolution and epidemic success of Mycobacterium tuberculosis in eastern China: evidence from a prospective study

The Bioinformatics Analysis Of Comparative Genomics Of Mycobacterium Tuberculosis Complex (Mtbc) Provides Insight Into Dissimilarities Between Intraspecific Groups Differingin Host Association, Virulence, And Epitope Diversity

Distribution and transmission of M. tuberculosis in a high-HIV prevalence city in Malawi: a genomic and spatial analysis

Local adaptive evolution of two distinct clades of Beijing and T families of Mycobacterium tuberculosis in Chongqing: a Bayesian population structure and phylogenetic study

Mycobacterium tuberculosis complex drug-resistance, phylogenetics, and evolution in Nigeria: Comparison with Ghana and Cameroon

Onset of infectiousness explains differences in transmissibility across lineages

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread

A novel Ancestral Beijing sublineage of Mycobacterium tuberculosis suggests the transition site to Modern Beijing sublineages

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