Evolutionary dynamics of the insertion sequence IS6110 in the Mycobacterium tuberculosis complex
- 1Swiss TPH
Abstract
Insertion sequences (IS) are the most common type of transposable element in prokaryotes and shape the structure of genomes through transposition and by providing a substrate for recombination. Despite the mutational impact of IS, the evolutionary dynamics of most elements in host species remain unknown. Here we study the dynamics of IS6110 in 10,000 strains of the Mycobacterium tuberculosis complex (MTBC). We developed a tool that allows the detection and comparison of IS insertions from short reads without using a reference genome. Using ancestral state reconstruction (ASR) on presence-absence patterns of IS6110, we describe the distribution of copy numbers (CNs) in the MTBC, infer birth rates of the element, and identify genomic regions with large numbers of parallel IS6110 insertions. Copy numbers in the MTBC range from 1 in some clades to more than 30 in strains of La3 (M. orygis). IS6110 birth rates scale approximately linearly with copy number and are elevated on terminal branches, consistent with the delayed action of purifying selection. A key characteristic of IS6110 is its occurrence in hotspots: the 5% most frequently targeted regions account for half of all independent insertion events. The motif 5'-TCTCAAAW-3' is enriched around target sites and in hotspots, suggesting that the accumulation of insertions in these regions results through a combination of non-random insertion and purifying selection in other regions. To conclude the study, we propose a niche constraints model according to which the distribution of IS6110 in the MTBC is governed by the rarity of regions that have both suitable DNA properties and little functional value for the host.
Related
- Variation in multiple classes of simple sequence repeats can alter drug susceptibility in Mycobacterium tuberculosis
- Two potent, functionally redundant MazF toxins inactivate the 16S rRNA anti-Shine-Dalgarno sequence and are differentially expressed in Mycobacterium tuberculosis
- Mycobacterium tuberculosis gene expression dynamics in host environments with single-cell resolution
- Schistosoma mansoni infection enhances insertion sequence activity in the host gut microbiome
- Genome-wide evolutionary shifts shape the emergence of dominant circulating clones in Mycobacterium abscessus
BibTeX
@misc{stritt2026evolutionary,
title = {{Evolutionary dynamics of the insertion sequence IS6110 in the Mycobacterium tuberculosis complex}},
author = {C. Stritt and C. Loiseau and S. Kalkan and S. Borrell and D. Brites and S. Gagneux},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.24.753865},
url = {https://www.biorxiv.org/content/10.64898/2026.09.24.753865v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Evolutionary dynamics of the insertion sequence IS6110 in the Mycobacterium tuberculosis complex
AU - C. Stritt
AU - C. Loiseau
AU - S. Kalkan
AU - S. Borrell
AU - D. Brites
AU - S. Gagneux
PY - 2026
DA - 2026/09/24
PB - bioRxiv
DO - 10.64898/2026.09.24.753865
UR - https://www.biorxiv.org/content/10.64898/2026.09.24.753865v1
AB - Insertion sequences (IS) are the most common type of transposable element in prokaryotes and shape the structure of genomes through transposition and by providing a substrate for recombination. Despite the mutational impact of IS, the evolutionary dynamics of most elements in host species remain unknown. Here we study the dynamics of IS6110 in 10,000 strains of the Mycobacterium tuberculosis complex (MTBC). We developed a tool that allows the detection and comparison of IS insertions from short reads without using a reference genome. Using ancestral state reconstruction (ASR) on presence-absence patterns of IS6110, we describe the distribution of copy numbers (CNs) in the MTBC, infer birth rates of the element, and identify genomic regions with large numbers of parallel IS6110 insertions. Copy numbers in the MTBC range from 1 in some clades to more than 30 in strains of La3 (M. orygis). IS6110 birth rates scale approximately linearly with copy number and are elevated on terminal branches, consistent with the delayed action of purifying selection. A key characteristic of IS6110 is its occurrence in hotspots: the 5% most frequently targeted regions account for half of all independent insertion events. The motif 5'-TCTCAAAW-3' is enriched around target sites and in hotspots, suggesting that the accumulation of insertions in these regions results through a combination of non-random insertion and purifying selection in other regions. To conclude the study, we propose a niche constraints model according to which the distribution of IS6110 in the MTBC is governed by the rarity of regions that have both suitable DNA properties and little functional value for the host.
ER -