T4SS–CAZy Environmental HGT

Raw report: t4ss_cazy_environmental_hgt Feeds into: chromosomal-and-integrative-gene-transfer, environmental-resistome, gene-cooccurrence-ecological-guilds, gene-function-acquisition-depth, horizontal-gene-transfer-driven-innovation, metal-cross-resistance, pangenome-integration, phb-granule-ecology, scale-dependent-mobile-element-associations

Overview

This preliminary report analyzes type IV secretion system (T4SS) and conjugative machinery in high-quality environmental metagenome-assembled genomes (MAGs), carbohydrate-active enzyme (CAZy) gene neighborhoods, and horizontal gene transfer (HGT) in GT2 glycosyltransferases. Core analyses NB01–NB04 are complete; NB05 threshold validation and NB06 manuscript figures remain pending. The report labels the work preliminary on this basis, yet it also presents NB05 analysis results, summarized below. ₁

Key Findings

• T4SS prevalence: 6,652 of 30,497 high-quality environmental MAGs (21.8%) carry T4SS or conjugative machinery, using a multi-marker definition comprising VirB4/6/8/9/10/11, VirD4, TraI/D, TrwB, and TraG. ₁
• CAZy co-occurrence: 92 CAZy families show elevated co-occurrence with T4SS loci at distances of ≤10 kb, although threshold validation is pending. GT2 glycosyltransferases are the top hit, occurring in 767 genomes with a reported average of 5,041 bp. ₁
• Biome enrichment: T4SS–CAZy associations are enriched in marine sediment (OR=5.5, q<10⁻⁹⁸), barley rhizosphere (OR=10.4), and maize rhizosphere (OR=4.1). ₁
• GT2 HGT: The GT2 gene tree contains 77 detected HGT events, including 32 normalized high-confidence cross-phylum events. The strongest event, Node_4915, spans 8 phyla at a maximum divergence of 4.843. ₁
• Transfer context: CAZy genes were not detected on plasmids by ICEfinder; 12 integrative mobilizable elements (IMEs) occurred among the top 100 accumulators. T4SS-positive genomes have 10× higher mobile genetic element (MGE) density than other genomes (p<0.001), consistent with chromosomal or integrative transfer rather than plasmid mobilization. ₁
• Robustness: A 70/30 discovery/validation split reproduced the enrichment patterns across both sets. ₁

NB05 HGT Characterization

Among 77 events in Detected_HGT_Events.csv, 65 span 2 phyla and 12 span ≥3 phyla (15.6%). Node_4915 contains 35 genes, 82.9% syntenic, spans 8 phyla—WOR-3, Desulfobacterota, Patescibacteria, Bacteroidota, Firmicutes_A, Methanobacteriota, Bdellovibrionota, and Acidobacteriota—and has Max_Divergence = 4.843. Divergence and synteny are negatively correlated (Spearman ρ = −0.615, p<0.001), with more phylogenetically distant events having lower syntenic percentage. The report interprets this as consistent with sequence divergence after transfer and calls it the expected biological pattern; this is the report’s interpretation, not an established mechanism. ₁

The most-involved phyla across all detected events are Firmicutes_A (27 events), Pseudomonadota (22), Bacillota_A (19), and Actinomycetota (10). The report includes figures/fig_nb05_hgt_scatter.png. ₁

GT2 Neighborhood and Metal Resistance Results

In 376 genomes, the GT2 neighborhood was parsed as a list-format region: T4SS occurred in 503 neighborhood entries and GT2 in 495, confirming syntenic co-localization at the contig level. GH23, a murein lytic transglycosylase family, was the second most common CAZy family in GT2 neighborhoods, with 106 occurrences, suggesting that cell-wall-remodeling genes cluster with GT2–T4SS syntenic loci; this clustering is suggested rather than mechanistically established. The report includes figures/fig_nb05_neighbourhood_functions.png for the neighborhood 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). The report separately states that genomes with GT2 in T4SS-proximal neighborhoods carry 11× more metal-resistance genes. It presents this as an independent link between CAZy–T4SS synteny and a metal-resistance niche-breadth hypothesis; describing these genomes as hubs for GT2 horizontal transfer is interpretive, and the association is observational. ₁

Interpretation and Caveats

The observed T4SS co-localization with GT2 glycosyltransferase cassettes and phylogenetic incongruence in the GT2 gene tree provide positive evidence for 32 cross-phylum HGT events. The absence of plasmid-borne CAZy genes in the ICEfinder analysis, together with 10× higher MGE density in T4SS-positive genomes, is consistent with chromosomal or integrative mechanisms involving IMEs or integrative conjugative elements rather than plasmid mobilization. These results support the hypothesis that T4SS machinery mediates environmental dissemination of carbohydrate-active enzyme diversity across phylogenetically distant bacteria, but all associations are observational and require experimental validation for mechanistic confirmation. ₁

The report identifies four pending validation tasks: a synteny-threshold permutation test requiring 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)]. ₁

Slots Into

• pangenome-integration — cross-phylum GT2 HGT events and the proposed chromosomal or integrative transfer route extend the corpus’s treatment of gene mobility across genomic backgrounds. ₁
• environmental-resistome — T4SS-proximal GT2 neighborhoods are associated with 11× more metal-resistance genes, linking environmental HGT niches to resistance functions. ₁
• metal-cross-resistance — the 0.045 versus 0.004 metal-resistance-type comparison provides a new association between CAZy–T4SS neighborhoods and metal-resistance breadth. ₁
• chromosomal-and-integrative-gene-transfer — T4SS prevalence in 21.8% of high-quality environmental MAGs, the absence of plasmid-borne CAZy genes in ICEfinder output, 12 IMEs among the top 100 accumulators and 10× higher MGE density in T4SS-positive genomes are read as consistent with chromosomal or integrative transfer, pending validation. ₁
• horizontal-gene-transfer-driven-innovation — 77 GT2 HGT events (32 high-confidence cross-phylum), Node_4915 spanning 8 phyla, and the divergence–synteny correlation (Spearman ρ = −0.615) give gene-tree evidence for cross-phylum CAZy transfer. ₁
• scale-dependent-mobile-element-associations — CAZy–T4SS co-occurrence is defined at ≤10 kb with threshold validation still pending, and contig-level neighborhood co-localization is a distinct scale from genome-level metal-resistance comparisons. ₁
• gene-cooccurrence-ecological-guilds — GT2-neighborhood MAGs carry more metal-resistance types (0.045 versus 0.004), framed by the report as a niche-breadth hypothesis rather than a demonstrated guild. ₁