Unsupervised decomposition of massively parallel TraDIS-Xpress fitness data from Escherichia coli identifies distinct functional and expression modules of antimicrobial resistance determinants
- 1Quadram Institute Bioscience
Abstract
Linking genotype to phenotype is a central goal in biology, one way to do this is using transposon insertion sequencing (TIS) where the fitness impacts of mutations across a genome are assessed in different conditions. Traditional analyses use comparison of each gene response independently in control versus test conditions, so a compendium of hundreds of test conditions returns hundreds of separate ranked gene lists. To unlock the potential information hidden within large data sets, we applied independent component (ICA) analysis to 754 Escherichia coli BW25113 TraDIS-Xpress samples. These samples included a series of antibiotic and other stress condition exposures. ICA is designed to identify hidden structure within complex data and to extract meaningfully linked components. In this analysis each annotated gene was resolved into a gene body and 198 bp 5' and 3' features so that expression effects driven by the transposon's outwards facing promoter enter the matrix separately from gene disruption effects. The decomposition returned 99 components, which we term fxModules, which were defined as containing elements with similar responses across the whole compendium rather than per any individual experiment. Of the 754 stress samples, the antibiotic arm included 110 samples derived from exposure to 14 separate agents at multiple concentrations. Many of these have been highly studied allowing us to explore the accuracy of predicted responsive fxModules. This analysis showed excellent predictive power of the approach and the fxModules responsive to individual drugs typically reproduced known mechanisms at a resolution superior to analysis of single gene sets. For example, exposure to fosfomycin identified responsive fxModules containing genes which are known to be involved in uptake of the drug, cAMP-CRP control and exposure to gentamicin identified oxidative phosphorylation split across three fxModules responding with different magnitudes and directions at different concentrations. This recapitulates recent work studying the landscape of E. coli genes sensitive to this drug with high precision without any prior knowledge of functional annotation. Similarly precise predictions were made for other drugs. Responsive modules also identified novel groups of genes and linkages beyond the boundaries of defined pathways as being responsive to various stresses demonstrating the ability to predict new genotype-phenotype links. Inclusion of the regulatory regions also resolved fxModules which were entirely contributed by promoter regions, for example for chloramphenicol genes responsible for efflux were identified by signals in intergenic control regions, and unidirectional insertions in gene bodies, a signal that a gene level matrix would not be able to predict with accuracy. These data show that ICA decomposition provides an efficient and accurate annotation independent approach to identifying sets of genes with related functional impacts in any given condition in the compendium. The functional modular association of genes also help in predicting gene functions of unknown genes based on the group function. This approach allows rapid identification of novel genotype-phenotype links and unlocks the power of multiple comparisons of transposon mutagenesis data.
Related
- Whole-genome sequencing-based characterization of antimicrobial resistance determinants and phenotype–genotype concordance in uropathogenic Escherichia coli isolates from ICU patients in Saudi Arabia
- Genomic insight into the antibiotic-resistant determinants of Escherichia coli from small ruminants and manure in North West Province, South Africa
- Decoding candidemia in critically ill patients: unsupervised clustering identifies three unique phenotypes
- Pathogenome and Plasmid-Borne Antimicrobial Resistance Phenotypes in a Multidrug-Resistant O111:H8 Shiga Toxin-Producing Escherichia coli Strain
- Antimicrobial resistance in environment and wastewaters of a swine slaughterhouse – tracking ESKAPE bacteria and ESBL-/carbapenemase-producing Escherichia coli
BibTeX
@misc{yasir2026unsupervised,
title = {{Unsupervised decomposition of massively parallel TraDIS-Xpress fitness data from Escherichia coli identifies distinct functional and expression modules of antimicrobial resistance determinants}},
author = {M. Yasir and A. K. Turner and A. Telatin and S. Bastkowski and E. R. Holden and I. G. Charles and M. A. Webber},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.11.758300},
url = {https://www.biorxiv.org/content/10.64898/2026.10.11.758300v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Unsupervised decomposition of massively parallel TraDIS-Xpress fitness data from Escherichia coli identifies distinct functional and expression modules of antimicrobial resistance determinants
AU - M. Yasir
AU - A. K. Turner
AU - A. Telatin
AU - S. Bastkowski
AU - E. R. Holden
AU - I. G. Charles
AU - M. A. Webber
PY - 2026
DA - 2026/10/11
PB - bioRxiv
DO - 10.64898/2026.10.11.758300
UR - https://www.biorxiv.org/content/10.64898/2026.10.11.758300v1
AB - Linking genotype to phenotype is a central goal in biology, one way to do this is using transposon insertion sequencing (TIS) where the fitness impacts of mutations across a genome are assessed in different conditions. Traditional analyses use comparison of each gene response independently in control versus test conditions, so a compendium of hundreds of test conditions returns hundreds of separate ranked gene lists. To unlock the potential information hidden within large data sets, we applied independent component (ICA) analysis to 754 Escherichia coli BW25113 TraDIS-Xpress samples. These samples included a series of antibiotic and other stress condition exposures. ICA is designed to identify hidden structure within complex data and to extract meaningfully linked components. In this analysis each annotated gene was resolved into a gene body and 198 bp 5' and 3' features so that expression effects driven by the transposon's outwards facing promoter enter the matrix separately from gene disruption effects. The decomposition returned 99 components, which we term fxModules, which were defined as containing elements with similar responses across the whole compendium rather than per any individual experiment. Of the 754 stress samples, the antibiotic arm included 110 samples derived from exposure to 14 separate agents at multiple concentrations. Many of these have been highly studied allowing us to explore the accuracy of predicted responsive fxModules. This analysis showed excellent predictive power of the approach and the fxModules responsive to individual drugs typically reproduced known mechanisms at a resolution superior to analysis of single gene sets. For example, exposure to fosfomycin identified responsive fxModules containing genes which are known to be involved in uptake of the drug, cAMP-CRP control and exposure to gentamicin identified oxidative phosphorylation split across three fxModules responding with different magnitudes and directions at different concentrations. This recapitulates recent work studying the landscape of E. coli genes sensitive to this drug with high precision without any prior knowledge of functional annotation. Similarly precise predictions were made for other drugs. Responsive modules also identified novel groups of genes and linkages beyond the boundaries of defined pathways as being responsive to various stresses demonstrating the ability to predict new genotype-phenotype links. Inclusion of the regulatory regions also resolved fxModules which were entirely contributed by promoter regions, for example for chloramphenicol genes responsible for efflux were identified by signals in intergenic control regions, and unidirectional insertions in gene bodies, a signal that a gene level matrix would not be able to predict with accuracy. These data show that ICA decomposition provides an efficient and accurate annotation independent approach to identifying sets of genes with related functional impacts in any given condition in the compendium. The functional modular association of genes also help in predicting gene functions of unknown genes based on the group function. This approach allows rapid identification of novel genotype-phenotype links and unlocks the power of multiple comparisons of transposon mutagenesis data.
ER -