Labs RSS GitHub • GPL-3.0
CRV Science Labs • bio-suite v0.2.0 (Rust / GPL-3.0)

Bioinformatics, Systematics & Emergent Genome Complexity

High-throughput multi-threaded Rust worker-agent pipelines for molecular phylogenetics, structural genomics, target discovery, and information-theoretic sequence complexity.

← Return to CRV Science Labs / Genomics & Complexity Theory
#3 ★ 1★ Baseline / Review Note Genomics & Complexity Theory

Horizontal Gene Transfer & Mobile Resistome Baseline Scan: Tetranucleotide KL-Divergence Pinpoints a 20.2 kb Acinetobacter-Derived blaNDM-1 Integrative Conjugative Element in Pseudomonas

Executing the CRV Science Labs Rust bio-suite hgt-scan pipeline across a 90 kb Pseudomonas chromosome benchmark resolves an unameliorated 20.22 kb tRNA-Gly integrative conjugative element (GI-ICE-PAE-01) carrying ISAba125-bracketed blaNDM-1, bleMBL, and sul1 resistance cassettes.

By Bryan P. White • 9 min read • High (Graduate / Technical Monograph)
🧪 Reproducible Worker-Agent Pipeline Provenance (GPL-3.0)
Repository: https://github.com/bpwhite/crvscience-labs.git  |  Pipeline ID: hgt-scan-v1
CLI Command: bio-suite hgt-scan --out runs/hgt-scan-report.json

1. Abstract & Information-Theoretic Compositional Anomaly Principle

Bacterial chromosomes maintain a species-specific genomic signature—an equilibrium k-mer frequency spectrum and synonymous codon usage bias enforced by host DNA replication, methylation, and mismatch-repair machinery[1]. When a clinical or environmental isolate acquires foreign DNA via a conjugative plasmid, prophage, or integrative conjugative element (ICE), the newly inserted block stands out as an information-theoretic anomaly before evolutionary amelioration gradually homogenizes its base composition toward the host background[2].

To benchmark Pipeline 3 (bio-suite hgt-scan), we executed the multi-threaded Rust worker-agent across a 90,000 bp Pseudomonas aeruginosa benchmark chromosome segment (NZ_CP089412.1_Pseudomonas_aeruginosa_ENV26, overall background GC = 60.27%, core host flank GC = 73.94%) using a 5,000 bp sliding window and 2,500 bp step size (35 windows total). For each window, the engine computes the 256-dimensional Laplace-smoothed tetranucleotide (4-mer) probability distribution and evaluates the Kullback-Leibler Divergence (DKL in bits) and 64-codon Relative Synonymous Codon Usage (RSCU) cosine distance against the whole-contig background, followed by 10 kb flanking mobility boundary interrogation (tRNAscan-SE 2.0[3], VirSorter2 / PhiSpy, ISEScan[4], and MOB-suite) and antimicrobial resistance cargo screening (AMRFinderPlus[5] and CARD[6]).

2. Sliding-Window Tetranucleotide KL-Divergence & RSCU Telemetry

Across the 35 overlapping 5 kb windows, windows HGT-W01 through HGT-W13 and HGT-W23 through HGT-W35 reflect the canonical GC-rich Pseudomonas host backbone (GC = 73.36%–74.54%, negative GC-skew = −0.0641 to −0.0797, baseline DKL = 0.3655–0.3836 bits). Between coordinates 35,000 bp and 55,000 bp (windows HGT-W15 through HGT-W21), the compositional profile undergoes a sharp phase transition: GC content drops to 12.16%–12.66%, GC-skew inverts to +0.4381 to +0.5635 (Z = +2.2240), and tetranucleotide Kullback-Leibler divergence surges more than fivefold to 1.9437–1.9759 bits (Z = +1.9640 to +2.0144) alongside an RSCU cosine distance of 0.4619–0.6902. Click Copy for Excel (TSV) to export the sliding-window compositional telemetry:

📊 Pipeline 3 (hgt-scan-v1) — 5 kb Sliding-Window Tetranucleotide KL-Divergence & RSCU Profile
Window ID Start (bp) End (bp) GC (%) GC Skew GC Skew Z 4-mer KL Div (bits) KL Div Z RSCU Cosine Dist Candidate Island
HGT-W010500073.94-0.06950.51940.3671-0.50250.0676false
HGT-W07150002000073.62-0.06490.49950.3711-0.49620.0676false
HGT-W13300003500073.88-0.07040.52330.3693-0.49910.0676false
HGT-W14325003750043.380.02440.11240.3276-0.56430.6377false
HGT-W15350004000012.660.51662.02071.94371.96400.6902true
HGT-W16375004250012.600.43811.68051.97592.01440.4619true
HGT-W17400004500012.460.54092.12601.96421.99610.5850true
HGT-W18425004750012.280.56352.22401.96501.99730.6811true
HGT-W19450005000012.160.49671.93451.97252.00910.4624true
HGT-W20475005250012.180.51071.99521.97272.00940.5835true
HGT-W21500005500012.520.49521.92801.97162.00760.6766true
HGT-W22525005750043.320.00920.17830.3320-0.55740.0200false
HGT-W23550006000074.04-0.06860.51550.3655-0.50500.5308false
HGT-W30725007750073.92-0.06820.51380.3673-0.50220.4670false
HGT-W35850009000073.60-0.06410.49600.3718-0.49520.5288false

3. Mobility Boundary Mapping & Resistome Cargo Interrogation (GI-ICE-PAE-01)

In Step 2 and Step 3, bio-suite hgt-scan scans the 10 kb boundaries flanking the compositional anomaly (34,890–55,110 bp; 20.22 kb total) and screens internal open reading frames against AMRFinderPlus, CARD, and VFDB. The island is precisely bracketed by an 18 bp attL/attR target-site duplication at the 3′ end of tRNA-Gly(GCC) alongside a XerC/IntI1 tyrosine site-specific recombinase (PF00589) and a MOBP1 conjugative relaxase, classifying GI-ICE-PAE-01 as an active integrative conjugative element:

📋 GI-ICE-PAE-01 Flanking Mobility Markers & Mobilized Resistome / Virulence Cargo
Feature / ORF ID Start (bp) End (bp) Strand Symbol / Element Category / Drug Class Database / Source Identity (%) E-Value
BOUNDARY-013489034965+tRNA-Gly(GCC) attLMobility Integration SitetRNAscan-SE 2.0100.001.4e-19
BOUNDARY-023542036410+intI1 / XerCTyrosine Recombinase (PF00589)VirSorter2 / PhiSpy98.403.2e-64
BOUNDARY-034120042285+ISAba125DDE Transposase (Tn125 composite)ISEScan99.808.5e-52
ENV26_ORF_00414245043263+blaNDM-1Carbapenems / Beta-LactamsAMRFinderPlus / CARD ARO:300058999.630.0
ENV26_ORF_00424331043675+bleMBLGlycopeptide / BleomycinCARD ARO:3002819100.001.1e-84
ENV26_ORF_00464680047610+sul1Sulfonamides (Dihydropteroate synthase)CARD ARO:300041099.280.0
ENV26_ORF_00494920050640+feOAB-tonB-virSiderophore Iron AcquisitionVFDB VFG00184296.404.5e-140
BOUNDARY-045180053120-mobP1 / traGConjugative Relaxase & T4SS CouplingMOB-suite / PlasmidFinder97.102.1e-48
BOUNDARY-055511055128+tRNA-Gly(GCC) attR18 bp Direct Repeat RestorationtRNAscan-SE 2.0100.001.8e-18

4. Mechanistic Synthesis & 1★ Signal Rating Assessment

The co-occurrence of an unameliorated AT-rich tetranucleotide signature (DKL = 1.9759 bits) with an upstream ISAba125 insertion sequence and downstream bleMBL gene reproduces the canonical Tn125 transposon architecture originally mobilized from Acinetobacter baumannii into Enterobacterales and Pseudomonas spp.[5, 6]. Because the island has not yet undergone mutational amelioration toward the high-GC Pseudomonas background, bio-suite hgt-scan flags the transfer event directly from primary sequence composition prior to phenotypic susceptibility testing.

Signal Rating Assessment (1★ Baseline / Review Note): This report is rated 1★ because it serves as a deterministic baseline validation of the bio-suite hgt-scan pipeline on a single 90 kb reference benchmark contig (n = 35 windows). Achieving a 3★–5★ discovery rating requires executing the scanner across newly sequenced, uncharacterized clinical or environmental metagenome-assembled genomes (MAGs) and confirming novel mobile resistance elements via long-read junction spanning and conjugation assays.

References & Primary Sources

  1. [1] Karlin, S., Mrázek, J., & Campbell, A. M. (1997). Compositional biases of bacterial genomes and evolutionary implications. Journal of Bacteriology, 179(12), 3899–3913. https://doi.org/10.1128/jb.179.12.3899-3913.1997
  2. [2] Lawrence, J. G., & Ochman, H. (1997). Amelioration of bacterial genomes: rates of change and exchange. Journal of Molecular Evolution, 44(4), 383–397. https://doi.org/10.1007/PL00006158
  3. [3] Chan, P. P., Lin, B. Y., Mak, A. J., & Lowe, T. M. (2021). tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes. Nucleic Acids Research, 49(16), 9077–9096. https://doi.org/10.1093/nar/gkab688
  4. [4] Xie, Z., & Tang, H. (2017). ISEScan: automated identification of insertion sequence elements in prokaryotic genomes. Bioinformatics, 33(21), 3340–3347. https://doi.org/10.1093/bioinformatics/btx433
  5. [5] Feldgarden, M., Brover, V., Gonzalez-Escalona, N., et al. (2021). AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Scientific Reports, 11, 12728. https://doi.org/10.1038/s41598-021-91456-0
  6. [6] Alcock, B. P., Huynh, W., Chalil, R., et al. (2023). CARD 2023: expanded curation, support for machine learning, and resistome prediction at the Comprehensive Antibiotic Resistance Database. Nucleic Acids Research, 51(D1), D690–D699. https://doi.org/10.1093/nar/gkac920