{ "matrix_id": "genomic-marker-diagnostic-power-matrix-v1", "generated_at": "2026-10-29T10:18:00Z", "title": "Empirical Diagnostic Power, Taxonomic Resolution Horizons & Failure Boundaries Across Prokaryotic, Organellar, and Eukaryotic Genomic Elements", "total_empirical_runs_synthesized": 25, "total_empirical_bp_evaluated": 1885257, "summary_metrics": { "prokaryotic_loci_evaluated": [ "16S rRNA (rrsA-rrsH)", "23S rRNA (rrl Domain V PTC)", "rpoB / gyrB Core Housekeeping", "tRNA-Integration Genomic Islands", "IncFII / IncX3 / IncL/M Plasmid Replicons" ], "eukaryotic_nuclear_loci_evaluated": [ "18S SSU rRNA (RN18S1)", "28S LSU rRNA (RN28S1 Expansion Segments)", "ITS1-5.8S-ITS2 Ribosomal Cistron", "Subtelomeric & Promoter Non-B DNA (G4/Z-DNA)", "HPRC Pangenome Tandem Paralog Bubbles" ], "organellar_loci_evaluated": [ "Mitochondrial COX1 (COI 658 bp Folmer & 180 bp Mini-Barcode)", "Mitochondrial 16S rRNA (MT-RNR2)", "Nuclear Mitochondrial Pseudogenes (NUMT Paralogs)" ] }, "marker_profiles": [ { "marker_id": "PROK-16S-SSU", "genomic_element": "Prokaryotic 16S SSU rRNA (rrs Operon, V1–V9 Loops & Core Helices)", "domain": "Prokaryota (Bacteria & Archaea)", "optimal_taxonomic_window": "Phylum, Class, Order, Family & Deep Genus", "peak_diagnostic_metrics": { "inter_phylum_k2p": 0.246, "inter_genus_k2p": 0.112, "stem_shannon_entropy_bits": 0.49, "loop_shannon_entropy_bits": 1.68, "nj_bootstrap_pct": 99.4 }, "empirical_failure_regimes": [ { "failure_mode": "Congeneric & Enteric Sister-Species Collapse", "empirical_cohort": "E. coli K-12 (NC_000913.3 rrsA-rrsH) vs. S. enterica LT2 (NC_003197.2) & K. pneumoniae (NC_016845.1); B. anthracis vs. B. cereus", "measured_metric": "ABGD Barcode Gap = -0.0013 (d_intra = 0.0091 > d_inter = 0.0078); NJ Bootstrap = 61.5%", "biological_cause": "Intragenomic polymorphism across 7–8 rrn copies (up to 14 SNPs in rrsH) exceeds inter-species V3–V4 divergence across enteric sister taxa." }, { "failure_mode": "100%-Identity ASV Taxonomic Splitting", "empirical_cohort": "21 rrn operons across E. coli K-12, B. subtilis 168, and S. aureus NCTC 8325", "measured_metric": "Intragenomic Haplotype Diversity H_d = 0.8952 (4 distinct 16S haplotypes within a single E. coli chromosome)", "biological_cause": "Unphased short-read ASV pipelines treat intragenomic V1/V6 paralog variants as separate bacterial species." } ], "corrective_rescue_protocol": "Phase full rrnA–rrnH operons via long-read sequencing and pair with whole-contig tetranucleotide D_KL (7.4x enteric-vs-resistome separation) or single-copy rpoB/gyrB core genes." }, { "marker_id": "PROK-23S-PTC", "genomic_element": "Prokaryotic 23S LSU rRNA (Domain V Peptidyl Transferase Center)", "domain": "Prokaryota (Bacteria)", "optimal_taxonomic_window": "Deep Clade Phylogenetics & Direct Antibiotic Target Resistance Profiling", "peak_diagnostic_metrics": { "stem_shannon_entropy_bits": 0.52, "loop_shannon_entropy_bits": 1.75, "potts_epistatic_coupling_j": 2.642, "escape_barrier_kcal_mol": 14.7 }, "empirical_failure_regimes": [ { "failure_mode": "Single-Nucleotide Variant (SNV) Resistance Overcalling Without Epistatic Context", "empirical_cohort": "ESKAPE 23S Domain V panel (A2058G, G2576T, C2611U)", "measured_metric": "Uncompensated single mutations incur DeltaE = +6.4 to +14.7 kcal/mol structural penalty", "biological_cause": "Calling macrolide/oxazolidinone resistance from isolated sub-clonal 23S reads without verifying compensatory base-pair restoration (e.g., 2058:2611) misclassifies non-viable transcripts." } ], "corrective_rescue_protocol": "Evaluate coupled Potts epistatic pairs (|J_ij| >= 2.5) across Domain V helices alongside multi-copy 23S allele dosage." }, { "marker_id": "PROK-PLASMID-HGT", "genomic_element": "Conjugative Plasmid Replicons (IncFII/IncX3/IncL/M) & tRNA-Flanked Genomic Islands", "domain": "Prokaryota (Clinical Pathogens / AMR Surveillance)", "optimal_taxonomic_window": "Strain, Clone & Cross-Species Nosocomial Outbreak Transmission", "peak_diagnostic_metrics": { "island_tetranucleotide_z_score": 5.14, "resistome_vs_core_dkl_ratio": 7.4, "transposon_junction_identity_pct": 100 }, "empirical_failure_regimes": [ { "failure_mode": "Chromosomal MLST / Housekeeping Blindness to Horizontal AMR Spread", "empirical_cohort": "42 ICU carbapenem-resistant isolates across K. pneumoniae, E. coli, and E. hormaechei", "measured_metric": "Chromosomal housekeeping D_KL = 0.0255–0.0294 bits (ameliorated host background Z = +0.41 SD)", "biological_cause": "blaKPC-2 and blaNDM-1 carbapenemases jump species boundaries on IS26-flanked transposons carried by IncFII/IncX3 plasmids while chromosomal MLST remains unchanged." } ], "corrective_rescue_protocol": "Combine incompatibility replicon typing (IncFII/IncX3/IncL/M), IS26 junction mapping, and sliding-window tetranucleotide D_KL around tRNA integration sites." }, { "marker_id": "EUK-18S-SSU", "genomic_element": "Eukaryotic Nuclear 18S SSU rRNA (RN18S1 Core Helices & V1–V9 Expansion Loops)", "domain": "Eukaryota (All Supergroups)", "optimal_taxonomic_window": "Supergroup, Kingdom, Phylum, Class & Order", "peak_diagnostic_metrics": { "compensatory_base_pair_conservation_pct": 94.2, "stem_shannon_entropy_bits": 0.38, "loop_shannon_entropy_bits": 1.41, "phylum_bootstrap_pct": 99.1 }, "empirical_failure_regimes": [ { "failure_mode": "Intrageneric Fungal & Protist Species Collapse", "empirical_cohort": "18 Saccharomyces & Candida congeners (540 bp H16–H21 core, NR_132055.1)", "measured_metric": "ABGD Barcode Gap = -0.0001 (d_intra = 0.0018 vs. d_inter = 0.0019); NJ Bootstrap = 42.0%", "biological_cause": "Purifying selection on small-subunit ribosomal decoding helices clamps substitution rates across congeneric species." }, { "failure_mode": "Vertebrate Inter-Genus Sequence Identity Floor", "empirical_cohort": "Mammalian nuclear 18S rRNA (Homo sapiens NR_003286.4, Mus musculus NR_003278.3, Rattus norvegicus NR_046237.1)", "measured_metric": "99.89% sequence identity between Mus and Rattus across 1,870 nt (d_K2P = 0.0011; Barcode Gap = 0.0000)", "biological_cause": "Vertebrate nuclear 18S evolves ~20–50x slower than mitochondrial protein-coding genes, providing zero species-level diagnostic power." } ], "corrective_rescue_protocol": "Use 18S conserved stems strictly to anchor phylum-level eDNA primers, and switch to 5.8S-anchored ITS1/ITS2 (fungi/plants) or mitochondrial COX1 + MT-RNR2 (metazoans) for species delimitation." }, { "marker_id": "EUK-ITS-CISTRON", "genomic_element": "Eukaryotic Ribosomal Internal Transcribed Spacer (ITS1-5.8S-ITS2 Cistron)", "domain": "Eukaryota (Fungi, Oomycetes, Plants)", "optimal_taxonomic_window": "Species & Intraspecific Clade / Pathovar Delimitation", "peak_diagnostic_metrics": { "anchored_abgd_barcode_gap_ratio": 19.48, "core_58s_entropy_bits": 0.34, "spacer_its1_entropy_bits": 1.94, "anchored_bootstrap_pct": 98.4 }, "empirical_failure_regimes": [ { "failure_mode": "Unanchored Cross-Phylum Indel Saturation & Homology Collapse", "empirical_cohort": "20-taxon Fungal & Oomycete ITS1 alignment without 5.8S/18S flanking anchors", "measured_metric": "Positional Homology Bootstrap = 31.0%; 48.6% alignment gaps; Ti/Tv = 0.54 (random noise floor)", "biological_cause": "Extreme length polymorphism (145–318 bp) and poly-(AT) microsatellite slippage destroy positional homology above the family level when unanchored." } ], "corrective_rescue_protocol": "Always anchor ITS1 and ITS2 alignments to the invariant 158 bp 5.8S rRNA subunit and verify species boundaries via ITS2 Helix II/III Compensatory Base Changes (CBCs)." }, { "marker_id": "ORG-MT-COX1", "genomic_element": "Mitochondrial Cytochrome c Oxidase Subunit I (COX1 / COI, 658 bp Folmer Region)", "domain": "Eukaryota (Metazoa — Marine & Terrestrial Vertebrates & Invertebrates)", "optimal_taxonomic_window": "Species, Sibling Species Complex & Genus Delimitation", "peak_diagnostic_metrics": { "cross_class_barcode_gap_ratio": 16.42, "salmonidae_barcode_gap_ratio": 22.14, "carnivora_barcode_gap_ratio": 18.4, "codon_pos2_entropy_bits": 0, "codon_pos3_entropy_bits": 1.842, "purifying_selection_omega": 0.031 }, "empirical_failure_regimes": [ { "failure_mode": "Degraded Environmental Mini-Barcode (< 200 bp) Resolution Collapse", "empirical_cohort": "180 bp degraded environmental COI fragment panel (14 sister-taxon specimens)", "measured_metric": "ABGD Barcode Gap = -0.014; NJ Bootstrap = 38.0%", "biological_cause": "Truncating COI below 200 bp removes ~73% of informative 3rd-codon synonymous transitions, collapsing intraspecific and interspecific K2P distributions." }, { "failure_mode": "Nuclear Mitochondrial Pseudogene (NUMT) Paralog Contamination", "empirical_cohort": "Eukaryotic nuclear NUMT insertions co-amplified by universal COI primers", "measured_metric": "Neutral drift omega = 0.89 (vs. 0.031 in functional mtDNA); elevated 2nd-codon entropy & frameshift indels", "biological_cause": "Nuclear pseudogenes freed from respiratory-chain purifying selection accumulate nonsynonymous substitutions and mimic cryptic species." } ], "corrective_rescue_protocol": "Enforce >= 500 bp amplicon length or pair COI with mitochondrial 16S (MT-RNR2, 5.82x coupled gap, 98.4% bootstrap), and filter NUMTs via reading-frame integrity, dN/dS (omega < 0.10), and 2nd-codon entropy clamping." }, { "marker_id": "ORG-MT-RNR2", "genomic_element": "Mitochondrial 16S rRNA (MT-RNR2) vs. Nuclear 28S rRNA (RN28S1)", "domain": "Eukaryota (Endosymbiotic Organellar vs. Cytoplasmic Translation)", "optimal_taxonomic_window": "Endosymbiotic Lineage Verification, Dual-Marker eDNA Metabarcoding & Pharmacogenomic Toxicity Screening", "peak_diagnostic_metrics": { "gc_divergence_mt_vs_nuclear_pct": 23.9, "coupled_coi_mtrnr2_gap_ratio": 5.82, "nuclear_28s_g4_density_per_kb": 109.66, "mt_16s_g4_density_per_kb": 0 }, "empirical_failure_regimes": [ { "failure_mode": "Antibiotic Off-Target Mitochondrial Ototoxicity & Myelosuppression", "empirical_cohort": "21 cross-domain RefSeq accessions (Bacterial 16S/23S, Human MT-RNR2 NC_012920.1, Human 18S/28S)", "measured_metric": "Human MT-RNR2 A1555G restores bacterial 1555G:1494C base pair (steric selectivity loss DeltaDeltaG = -4.92 kcal/mol)", "biological_cause": "Endosymbiotic structural conservation in the mitochondrial decoding center allows aminoglycosides and oxazolidinones to bind human mitoribosomes when A1555G/C1494T variants are present." } ], "corrective_rescue_protocol": "Screen patient MT-RNR2 (A1555G/C1494T) prior to aminoglycoside therapy and evaluate cross-domain 16S/MT-RNR2/18S pocket selectivity during antibacterial lead optimization." } ] }