Chromosomal and integrative routes of horizontal gene transfer (HGT) describe gene dissemination associated with chromosomal neighborhoods, integrative mobilizable elements (IMEs), integrative conjugative elements, and type IV secretion system (T4SS) or conjugative machinery rather than plasmid carriage. 1
Evidence from Environmental MAGs
Among 30,497 high-quality environmental metagenome-assembled genomes (MAGs), 6,652 (21.8%) carried T4SS or conjugative machinery under a multi-marker definition including VirB4/6/8/9/10/11, VirD4, TraI/D, TrwB, and TraG. 1 T4SS-positive genomes had 10× higher mobile genetic element (MGE) density than other genomes, with p<0.001. 1
CAZy genes were not detected on plasmids by ICEfinder, while 12 IMEs occurred among the top 100 accumulators. 1 This combination supports, but does not establish, the hypothesis that chromosomal or integrative mechanisms contribute to dissemination of carbohydrate-active enzyme diversity when plasmid mobilization is not detected. 1
T4SS–CAZy Neighborhoods
The report identified 92 CAZy families with elevated co-occurrence within ≤10 kb of T4SS loci, although threshold validation remained pending. 1 GT2 glycosyltransferases were the top hit, occurring in 767 genomes; the accompanying figure is given only in abbreviated form as an average of 5,041 bp for the T4SS co-occurrence distance, and does not support a gene-length reading. 1
In a parsed set of 376 genomes, T4SS occurred in 503 GT2-neighborhood entries and GT2 occurred in 495 entries, confirming co-localization at the contig level. 1 GH23, a murein lytic transglycosylase family, occurred 106 times in GT2 neighborhoods, suggesting that cell-wall-remodeling genes cluster with GT2–T4SS syntenic loci. 1
The associations were enriched in marine sediment (OR=5.5, q<10⁻⁹⁸), barley rhizosphere (OR=10.4), and maize rhizosphere (OR=4.1). 1 A 70/30 discovery/validation split reproduced the enrichment patterns across both sets, supporting the robustness of the observed associations while not resolving their mechanism. 1
Cross-Phylum HGT Evidence
The GT2 gene tree contained 77 detected HGT events, including 32 normalized high-confidence cross-phylum events. 1 Node_4915 was the strongest event, spanning 8 phyla at a maximum divergence of 4.843. 1 It contained 35 genes, was 82.9% syntenic, and included WOR-3, Desulfobacterota, Patescibacteria, Bacteroidota, Firmicutes_A, Methanobacteriota, Bdellovibrionota, and Acidobacteriota. 1
Among the 77 events, 65 spanned 2 phyla and 12 spanned ≥3 phyla (15.6%). 1 Divergence and synteny were negatively correlated (Spearman ρ = −0.615, p<0.001), with more phylogenetically distant events having lower syntenic percentage. 1 These gene-tree and synteny results provide positive evidence for cross-phylum transfer, but the report treats the proposed T4SS-mediated mechanism as an observational hypothesis requiring experimental validation. 1
The cross-species COG analysis supports the broader HGT interpretation: across 32 species spanning 9 phyla and 357,623 genes, novel or singleton genes were enriched in COG L, the mobile-element category, by +10.88% with 100% consistency. 2 This functional pattern is consistent with mobile elements being a major source of genomic novelty, but it refines rather than proves the T4SS claim: it supports mobile-element-associated innovation at the comparative-genomic level without identifying T4SS or an integrative route as the causal mechanism. 2
The costly-plus-dispensable analysis supports this broader mobile-element interpretation: among 5,526 genes across 142,190 genes from 43 bacteria, costly, non-conserved genes were 7.45× more likely than costly-conserved genes to contain mobile-element keywords (OR=7.45, p=4.6e-71), and had poorer annotation and narrower ortholog breadth. 3 Their enrichment in transposase, integrase, phage, insertion-sequence, recombinase, and prophage signatures is consistent with recent HGT-associated genome expansion, but it refines the present claim rather than identifying a chromosomal, IME, or T4SS route: the analysis did not establish which mobile mechanism transferred these genes. 3
The gene-function ecological atlas supports a chromosomal or integrative reading while making the route’s provenance explicit: its synthesis reports 3 of 3 confirmed hypotheses anchored in expected biomes and phenotypes, with pre-registered KOs (KEGG Orthology gene families) found not to be phage-borne, and with integrative conjugative elements (ICEs) and chromosomal recombination implicated by literature rather than by the project’s own data. 4
Chromosomal Integration versus Mobile Carriage of Resistance Genes
Counter-intuitively, non-core (putatively acquired) antimicrobial resistance (AMR) genes showed stronger phylogenetic signal (median r = 0.222) than core (intrinsic) genes (median r = 0.117), a highly significant difference (paired t-test: t = -8.35, p = 7.0e-16, n = 489). 5 The authors interpret this as suggesting that so-called “acquired” AMR genes are more often inherited clonally within lineages than randomly acquired via horizontal gene transfer. 5 The same report records an alternative, statistical explanation for part of the gap: core AMR genes, being nearly universal (>=95% prevalence by definition), have near-zero Jaccard distances with little variance, which inherently suppresses distance-based correlation metrics like the Mantel r (a Mantel test assesses correlation between two distance matrices by permutation) — an artifact that partially explains the lower core signal independent of biology. 5
The clonal reading rests on a proposed model rather than a direct measurement: the report advances an alternative in which a lineage acquires resistance genes (likely via mobile genetic elements) that then become stably integrated into the genome and are vertically inherited, creating clonal AMR lineages — a model not directly established by its Mantel-test associations. 5 Its literature context separately attributes to Maier et al. (2025) and others the finding that AMR gene transfer is largely confined within closely related lineages, and to Maier et al. (2025) specifically a negative influence of genetic compatibility on cross-species AMR transfer; these are secondary-literature findings rather than this project’s own measurements. 5
As reported for E. coli by Vanacker, Lenuzza & Rasigade (2023), plasmid-borne horizontally transferred resistance genes are less costly than chromosomal mutations; the fitness-cost project invokes that external result to explain why its own transposon-based cost estimate is modest, reasoning that many Fitness Browser AMR genes are likely acquired determinants such as beta-lactamases and acetyltransferases rather than costly target modifications — an interpretation, not a measured comparison of carriage routes. 6 Low-cost plasmid carriage and chromosome-like phylogenetic signal therefore constrain different stages of the same process, acquisition and retention, and neither alone identifies the route. 6,5
How fast resistance genes move in and out of genomes is not settled by this corpus: the AMR pangenome atlas poses temporal analysis — using phylogenetic tree distances to estimate AMR gene gain/loss rates and ask whether acquired resistance genes are being gained faster than lost — as a future question and reports no rate estimate. 7
Association with Metal-Resistance Functions
GT2-neighborhood MAGs (n=376) had a mean of 0.045 metal-resistance types, compared with 0.004 for non-GT2 MAGs (n=260,276; Mann–Whitney p=8.6e-27). 1 Genomes with GT2 in T4SS-proximal neighborhoods carried 11× more metal-resistance genes. 1 This finding supports a link between putative chromosomal or integrative HGT hubs and resistance-function breadth, extending the evidence summarized in environmental-resistome and metal-cross-resistance. 1
The metal cross-resistance analysis refines this interpretation by showing that metal-important genes can share a broadly positive fitness architecture: 98.1% of 317 organism–metal-pair observations were positive across 28 organisms, while general-stress genes, metal-shared genes, and metal-specific genes differed only modestly in mean pangenome core fraction (92.0%, 91.0%, and 89.8%, respectively). 8 Thus, the broader resistance phenotype near GT2–T4SS neighborhoods may reflect shared stress functions as well as transferred or co-localized resistance genes; the study did not include non-metal stress controls that would separate these explanations. 8
Tensions and Limits
The absence of CAZy genes on plasmids in the ICEfinder analysis does not demonstrate that plasmids never carry or mobilize these genes, because the result is a database- and detection-dependent negative observation. 1 Likewise, T4SS–CAZy proximity, MGE density, phylogenetic incongruence, and IME occurrence do not independently prove that T4SS machinery caused the observed transfers. 1 The report therefore supports a chromosomal or integrative transfer model without establishing a direct mechanistic route. 1
The costly-plus-dispensable result refines this limit: mobile-element enrichment, short gene length, poor annotation, and narrow distribution support HGT-associated novelty, but the evidence is compatible with multiple mobile routes and does not distinguish chromosomal or integrative transfer from other mechanisms. 3 Its burden classification is also sensitive to noise because “burden” is defined as max_fit > 1 in any experiment. 3
The metal cross-resistance study provides a related limitation rather than a direct contradiction: multi-metal tolerance scores did not correlate with BacDive isolation from metal environments at the Fitness Browser species scale (Spearman rho approximately -0.02, p > 0.8), but matching and strain collapsing left 20 independent species, making the test underpowered. 8 This refines the environmental-resistome interpretation: local GT2-neighborhood enrichment remains a measured association, but it should not be treated as evidence that cross-resistance signatures generally predict metal-associated isolation environments. 1,8
The COG analysis is limited by approximately 70% annotation coverage, so unassigned genes may skew functional distributions; its 32-species comparison may also miss phylum-specific patterns. 2 These limitations further constrain using COG L enrichment to identify the specific transfer route. 2
Two projects read resistance-gene phylogenetic patterns in different directions, and at different taxonomic scales. The metal-ecology analysis states that ecological range is deeply canalized by lineage identity while metal AMR is more labile, likely because AMR genes reside on mobile genetic elements and are frequently exchanged among distantly related taxa — an attribution that is an interpretation, not a direct measurement of horizontal gene transfer in that analysis. 9 This sits in tension with the strain-variation reading that non-core AMR genes are more often inherited clonally within lineages than randomly acquired via horizontal gene transfer, though the two are not strictly exclusive: exchange among distant taxa and subsequent clonal retention within a lineage describe different stages and scales and could coexist. 5 Both positions are interpretations layered on correlational evidence, so the tension should not be averaged away: cross-taxon lability and within-lineage retention are separable predictions that would need gain/loss-rate estimates on a shared tree to distinguish. 9,5
Donor identification at deep taxonomic ranks was deferred in the gene-function ecological atlas because per-CDS sequence is not present in the queryable schemas of the KBase Data Lakehouse, and composition-based donor inference would require external FASTA plus codon profiling; the project therefore reports at a joint Innovator-Exchange label that leaves broker and open donors undistinguished. 4
Its gene-neighborhood mobile-genetic-element (MGE) cargo scan also cannot detect ICE-mediated transfer of Bacteroidota polysaccharide utilization loci (PUL), because bakta product-keyword matching detects MGE machinery genes but does not detect cargo genes carried on MGEs without synteny analysis. 4 These two limits refine the central caveat of this page: a route that goes undetected may simply be invisible to the available schemas and annotation keywords. 4
The report identified four pending validation tasks: a synteny-threshold permutation test using unfiltered Spark data, BLAST validation of Node_4915 against NCBI nr, a housekeeping-gene null baseline, and biome-enrichment factorization using θ = OR(T4SS-CAZy) / [OR(T4SS) × OR(CAZy)]. 1
These findings extend pangenome-integration by connecting cross-phylum GT2 transfer evidence with chromosomal and integrative genomic context. 1 The COG result further supports this interpretation by associating genomic novelty with mobile-element functions, while leaving the T4SS-mediated route unresolved. 2 The costly-plus-dispensable analysis supports the mobile-novelty component but leaves the route unresolved, while the metal cross-resistance result further refines the resistance association by identifying shared-stress architecture as an alternative explanation that must be controlled in transfer-focused analyses. 3,8 The source reports are summarized at T4SS–CAZy Environmental HGT, COG Functional Category Analysis, Gene-Resolution Metal Cross-Resistance Across Diverse Bacteria, The 5,526 Costly + Dispensable Genes, Within-Species AMR Strain Variation, Fitness Cost of Antimicrobial Resistance Genes, Pan-Bacterial AMR Gene Landscape, Gene Function Ecological Agora, and Metal Resistance Ecology — Phylogenetic Conservation vs. Environmental Selection.
Open Directions
- Use the unfiltered Spark data for a synteny-threshold permutation test to determine whether the ≤10 kb T4SS–CAZy co-occurrence exceeds a distance-matched null expectation. 1
- BLAST Node_4915 against NCBI nr to test whether its 8-phylum distribution and maximum divergence of 4.843 are consistent with homologous transfer rather than annotation or tree-reconstruction artifacts. 1
- Compare GT2–T4SS neighborhoods with housekeeping-gene neighborhoods as a null baseline to test whether the observed cross-phylum signal is specific to GT2-associated loci. 1
- Apply the proposed biome-enrichment factorization, θ = OR(T4SS-CAZy) / [OR(T4SS) × OR(CAZy)], to test whether co-occurrence reflects an interaction beyond the independent distributions of T4SS and CAZy genes. 1
- Combine ICEfinder calls, IME annotations, contig context, and long-read or closed-genome validation to determine whether the observed GT2 neighborhoods are chromosomal, integrative, or plasmid-associated. 1
- Partition the COG L-enriched novel-gene signal by mobile-element subtype and compare it with T4SS/IME calls across the 32-species dataset to test whether mobile novelty is specifically coupled to chromosomal or integrative machinery. 2
- Link costly-plus-dispensable genes to contig context, ICEfinder/IME calls, and T4SS neighborhoods to test whether their mobile-element enrichment is specifically associated with chromosomal or integrative routes rather than mobile elements generally. 3
- Add non-metal stress controls and normalize metal experiments by concentration relative to MIC to test whether GT2-neighborhood resistance breadth reflects metal-specific transfer or a shared-stress fitness program. 8
- Carry out the deferred temporal analysis — estimating AMR gene gain/loss rates from phylogenetic tree distances — to test whether acquired resistance genes are gained faster than they are lost, the quantity neither the atlas nor the strain-variation analysis reports. 7
- Pair the non-core AMR phylogenetic-signal set (median r = 0.222; paired t-test n = 489) with contig-level ICEfinder and IME calls to test whether high-signal resistance genes sit in chromosomal or integrative contexts rather than on plasmids. 5,1
- Add synteny analysis to bakta-based MGE-machinery scans so that cargo genes carried on ICEs, including Bacteroidota PUL, become callable, closing the gap that currently prevents ICE-mediated transfer from being detected. 4
- Make per-CDS sequence available alongside the queryable schemas so composition-based donor inference can split the joint Innovator-Exchange label into broker and open donors at deep taxonomic ranks. 4