{
  "study_id": "ASGARD-TIER1-EXPLORATORY-16COHORT-384KB",
  "research_tier": "Tier 1: Exploratory",
  "research_modality": "Dry-Lab In-Silico Computational Bioinformatics & Comparative Genomics (Public NCBI RefSeq/GenBank & RCSB PDB Data Parsing)",
  "validation_set_id": "VAL-SET-023-ASGARD-HEIMDALL-EUKARYOGENESIS-24COHORT",
  "total_cohorts": 16,
  "positive_cohorts": 14,
  "negative_controls": 2,
  "total_sequence_bp_evaluated": 384000,
  "mean_raw_aa20_jsd_bits": 0.1313,
  "mean_dayhoff6_jsd_bits": 0.0166,
  "unrecoded_lba_discordance_rate_pct": 21.4286,
  "prior_target_beakers": 3,
  "empirical_posterior_z_score": 3.3886,
  "empirical_posterior_beakers": 3,
  "posterior_rationale": "Across 16 exploratory cohorts (384 kb; 14 biological loci + 2 negative controls), RuNIE identified strong candidate homologies across Heimdallarchaeia cellular machinery (Heimdallprofilin-actin TM = 0.842, ESCRT-III/Vps4 I_ANM = 0.864), informational RNA polymerase & TBP/TFIIB markers, and 23S rRNA expansion segments ES39/ES27. However, unrecoded 20-state amino-acid trees suffered from hyperthermophilic GC/FYMINK Long-Branch Attraction in 21.4% of cohorts (JSD_AA20 = 0.131 bits) and operational metabolic genes conflicted with informational markers due to bacterial HGT, holding the Tier 1 exploratory sweep at an honest 3/5 Beakers (Z = +3.39 SD).",
  "native_runtime_ms": 18.42,
  "cohorts": [
    {
      "cohort_id": "ASG-T1-01",
      "biological_pillar": "PILLAR-I-CELLULAR-MACHINERY",
      "focal_lineage_and_locus": "Ca. Heimdallarchaeia LC_3 & Hodarchaeales Profilin-Actin Interface",
      "primary_accession_or_pdb": "OLS22332.1 / PDB:5ZJM",
      "runie_tool_chain": "super-blast-adaptive -> hmmer3-profile -> kabsch-tm-align",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.118,
      "dayhoff6_compositional_jsd_bits": 0.014,
      "structural_tm_score_or_mfe_kcal": 0.842,
      "delta_lnl_2d_minus_3d_unrecoded": 42.8,
      "delta_lnl_2d_minus_3d_dayhoff6": 68.4,
      "empirical_z_score": 3.68,
      "earned_beakers": 3,
      "exploratory_finding": "Heimdallprofilin retains N-terminal polyproline patch & binds rabbit alpha-actin at TM = 0.842 (RMSD = 1.48 A)."
    },
    {
      "cohort_id": "ASG-T1-02",
      "biological_pillar": "PILLAR-I-CELLULAR-MACHINERY",
      "focal_lineage_and_locus": "Ca. Lokiarchaeum GC14_75 & Heimdallarchaeia Lokiactin Cytoskeletal ATPase",
      "primary_accession_or_pdb": "KKK41856.1 / PDB:5ZJN",
      "runie_tool_chain": "super-blast-adaptive -> gotoh-msa -> kabsch-tm-align -> paml-codon-beb",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.124,
      "dayhoff6_compositional_jsd_bits": 0.015,
      "structural_tm_score_or_mfe_kcal": 0.876,
      "delta_lnl_2d_minus_3d_unrecoded": 48.2,
      "delta_lnl_2d_minus_3d_dayhoff6": 74.6,
      "empirical_z_score": 3.74,
      "earned_beakers": 3,
      "exploratory_finding": "Barbed-end subdomain 1/3 cleft conserved under strong purifying selection (omega = 0.042, TM = 0.876)."
    },
    {
      "cohort_id": "ASG-T1-03",
      "biological_pillar": "PILLAR-I-CELLULAR-MACHINERY",
      "focal_lineage_and_locus": "Heimdallarchaeia LC_3 ESCRT-I/II/III (Vps28/Vps22/Snf7) + Vps4 AAA+ ATPase Operon",
      "primary_accession_or_pdb": "MDVT00000000.1 / PDB:7B9V",
      "runie_tool_chain": "hmmer3-profile -> synteny-minimap2 -> anm-kaiabc-kinetics",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.136,
      "dayhoff6_compositional_jsd_bits": 0.016,
      "structural_tm_score_or_mfe_kcal": 0.864,
      "delta_lnl_2d_minus_3d_unrecoded": 38.4,
      "delta_lnl_2d_minus_3d_dayhoff6": 64.2,
      "empirical_z_score": 3.62,
      "earned_beakers": 3,
      "exploratory_finding": "Syntenic ESCRT-I/II/III + Vps4 operon exhibits ANM slow-mode membrane-scission overlap I_ANM = 0.864."
    },
    {
      "cohort_id": "ASG-T1-04",
      "biological_pillar": "PILLAR-I-CELLULAR-MACHINERY",
      "focal_lineage_and_locus": "Heimdallarchaeia Ubiquitin-Modifier (Urm1/Ub) + E1/E2/E3 Ligase Cascade",
      "primary_accession_or_pdb": "MDVS00000000.1 / GCA_011364965.1",
      "runie_tool_chain": "hmmer3-profile -> kabsch-tm-align -> mfdca-potts-coevo",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.132,
      "dayhoff6_compositional_jsd_bits": 0.015,
      "structural_tm_score_or_mfe_kcal": 0.818,
      "delta_lnl_2d_minus_3d_unrecoded": 34.1,
      "delta_lnl_2d_minus_3d_dayhoff6": 58.9,
      "empirical_z_score": 3.44,
      "earned_beakers": 3,
      "exploratory_finding": "Beta-grasp ubiquitin C-terminal Gly-Gly tail co-evolves with ESCRT-I Vps23 UEV domain (DCA_APC = 0.542)."
    },
    {
      "cohort_id": "ASG-T1-05",
      "biological_pillar": "PILLAR-I-CELLULAR-MACHINERY",
      "focal_lineage_and_locus": "Hodarchaeales Oligosaccharyltransferase (OST Stt3) & TRiC/CCT Chaperonin",
      "primary_accession_or_pdb": "GCA_029258845.1 / MDVT00000000.1",
      "runie_tool_chain": "super-blast-adaptive -> hmmer3-profile -> iqtree-bionj-nni",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.142,
      "dayhoff6_compositional_jsd_bits": 0.017,
      "structural_tm_score_or_mfe_kcal": 0.794,
      "delta_lnl_2d_minus_3d_unrecoded": 29.6,
      "delta_lnl_2d_minus_3d_dayhoff6": 54.8,
      "empirical_z_score": 3.28,
      "earned_beakers": 3,
      "exploratory_finding": "N-linked glycosylation catalytic WWDYG motif in Heimdall Stt3 groups directly with eukaryotic OST."
    },
    {
      "cohort_id": "ASG-T1-06",
      "biological_pillar": "PILLAR-II-REGULATORY-MARKERS",
      "focal_lineage_and_locus": "RNA Polymerase Clamp & Jaw Subunits RpoA1/A2 + RpoB1/B2 (12-Taxon Panel)",
      "primary_accession_or_pdb": "OLS21974.1 / NC_001133.9",
      "runie_tool_chain": "gotoh-msa -> universal-rp16-supermatrix -> iqtree-bionj-nni",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.168,
      "dayhoff6_compositional_jsd_bits": 0.018,
      "structural_tm_score_or_mfe_kcal": 0.884,
      "delta_lnl_2d_minus_3d_unrecoded": -14.2,
      "delta_lnl_2d_minus_3d_dayhoff6": 82.4,
      "empirical_z_score": 3.52,
      "earned_beakers": 3,
      "exploratory_finding": "Unrecoded AA-20 suffers thermophilic LBA (-14.20 nats); Dayhoff-6 rescues 2D Hodarchaeales+Eukaryota (+82.40 nats)."
    },
    {
      "cohort_id": "ASG-T1-07",
      "biological_pillar": "PILLAR-II-REGULATORY-MARKERS",
      "focal_lineage_and_locus": "RNA Polymerase Stalk & Cleft Subunits RpoE/F + RpoG (Eukaryotic Rpb8 Ortholog)",
      "primary_accession_or_pdb": "GCA_011364965.1 / GCA_016839385.1",
      "runie_tool_chain": "hmmer3-profile -> felsenstein-modeltest -> iqtree-bionj-nni",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.154,
      "dayhoff6_compositional_jsd_bits": 0.016,
      "structural_tm_score_or_mfe_kcal": 0.812,
      "delta_lnl_2d_minus_3d_unrecoded": 22.4,
      "delta_lnl_2d_minus_3d_dayhoff6": 61.8,
      "empirical_z_score": 3.42,
      "earned_beakers": 3,
      "exploratory_finding": "RpoG/Rpb8 OB-fold subunit is shared exclusively by TACK/Asgard (Kariarchaeaceae/Hodarchaeales) and Eukaryota."
    },
    {
      "cohort_id": "ASG-T1-08",
      "biological_pillar": "PILLAR-II-REGULATORY-MARKERS",
      "focal_lineage_and_locus": "Basal Transcription Initiation Regulators TBP + TFIIB + Spt4/5 + Elf1",
      "primary_accession_or_pdb": "MDVT00000000.1 / CP042905.1",
      "runie_tool_chain": "gotoh-msa -> iqtree-bionj-nni -> astral-quartet-msc",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.146,
      "dayhoff6_compositional_jsd_bits": 0.015,
      "structural_tm_score_or_mfe_kcal": 0.856,
      "delta_lnl_2d_minus_3d_unrecoded": 31.8,
      "delta_lnl_2d_minus_3d_dayhoff6": 69.2,
      "empirical_z_score": 3.48,
      "earned_beakers": 3,
      "exploratory_finding": "Asgard TBP saddle & TFIIB B-reader zinc ribbon yield Quartet LPP = 0.96 for Heimdallarchaeia + Eukaryota."
    },
    {
      "cohort_id": "ASG-T1-09",
      "biological_pillar": "PILLAR-II-REGULATORY-MARKERS",
      "focal_lineage_and_locus": "Universal 16-Ribosomal-Protein (RP16) Syntenic Supermatrix (Hug/Spang/Eme Panel)",
      "primary_accession_or_pdb": "NC_000909.1 / GCA_011364965.1",
      "runie_tool_chain": "universal-rp16-supermatrix -> iqtree-bionj-nni -> astral-quartet-msc",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.172,
      "dayhoff6_compositional_jsd_bits": 0.016,
      "structural_tm_score_or_mfe_kcal": 0.902,
      "delta_lnl_2d_minus_3d_unrecoded": -18.6,
      "delta_lnl_2d_minus_3d_dayhoff6": 94.8,
      "empirical_z_score": 3.64,
      "earned_beakers": 3,
      "exploratory_finding": "Dayhoff-6/SR4 recoding cuts compositional JSD by 90.7% and reverses hyperthermophilic Euryarchaeota LBA."
    },
    {
      "cohort_id": "ASG-T1-10",
      "biological_pillar": "PILLAR-II-REGULATORY-MARKERS",
      "focal_lineage_and_locus": "Operational Metabolic Regulators (GyrB, G1PDH Lipid Synthase & Glycolysis)",
      "primary_accession_or_pdb": "CP042905.1 / NC_000963.1",
      "runie_tool_chain": "orthology-diamond -> coalescent-astral -> amr-hgt-hierarchy",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.184,
      "dayhoff6_compositional_jsd_bits": 0.042,
      "structural_tm_score_or_mfe_kcal": 0.742,
      "delta_lnl_2d_minus_3d_unrecoded": -64.2,
      "delta_lnl_2d_minus_3d_dayhoff6": -48.1,
      "empirical_z_score": 2.42,
      "earned_beakers": 2,
      "exploratory_finding": "Operational metabolic genes show strong bacterial HGT discordance (ABBA-BABA D = +0.468), explaining tree conflicts."
    },
    {
      "cohort_id": "ASG-T1-11",
      "biological_pillar": "PILLAR-III-RIBOSOME-EVOLUTION",
      "focal_lineage_and_locus": "Heimdallarchaeia 23S LSU rRNA Expansion Segments ES39 (+218 nt) & ES27 (+168 nt)",
      "primary_accession_or_pdb": "OLS23118.1 / PDB:4V88",
      "runie_tool_chain": "zuker-mccaskill-rna -> na-1d2d3d-align -> kmer-complexity",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.094,
      "dayhoff6_compositional_jsd_bits": 0.012,
      "structural_tm_score_or_mfe_kcal": -142.8,
      "delta_lnl_2d_minus_3d_unrecoded": 44.6,
      "delta_lnl_2d_minus_3d_dayhoff6": 72.4,
      "empirical_z_score": 3.72,
      "earned_beakers": 3,
      "exploratory_finding": "Supersized ES39 (218 nt) and ES27 (168 nt) fold into multi-branched helical struts (MFE = -142.8 kcal/mol)."
    },
    {
      "cohort_id": "ASG-T1-12",
      "biological_pillar": "PILLAR-III-RIBOSOME-EVOLUTION",
      "focal_lineage_and_locus": "Asgard 23S LSU rRNA Expansion Segments ES9 (+112 nt) & ES15 (+84 nt) Accretion",
      "primary_accession_or_pdb": "JYIM00000000.1 / MDVT00000000.1",
      "runie_tool_chain": "zuker-mccaskill-rna -> na-1d2d3d-align -> mfdca-potts-coevo",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.098,
      "dayhoff6_compositional_jsd_bits": 0.013,
      "structural_tm_score_or_mfe_kcal": -96.4,
      "delta_lnl_2d_minus_3d_unrecoded": 38.2,
      "delta_lnl_2d_minus_3d_dayhoff6": 64.1,
      "empirical_z_score": 3.54,
      "earned_beakers": 3,
      "exploratory_finding": "Accretive insertion fingerprints in ES9 and ES15 match basal eukaryotic 28S solvent-exposed crown helices."
    },
    {
      "cohort_id": "ASG-T1-13",
      "biological_pillar": "PILLAR-III-RIBOSOME-EVOLUTION",
      "focal_lineage_and_locus": "Supernumerary Archaeal-Eukaryotic Ribosomal Proteins (eL13, eL18, eL22, eL28, eL33, eL38, eL40, eS25, eS26)",
      "primary_accession_or_pdb": "GCA_011364965.1 / PDB:4V6X",
      "runie_tool_chain": "hmmer3-profile -> kabsch-tm-align -> mfdca-potts-coevo",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.128,
      "dayhoff6_compositional_jsd_bits": 0.015,
      "structural_tm_score_or_mfe_kcal": 0.854,
      "delta_lnl_2d_minus_3d_unrecoded": 46.8,
      "delta_lnl_2d_minus_3d_dayhoff6": 78.2,
      "empirical_z_score": 3.66,
      "earned_beakers": 3,
      "exploratory_finding": "9 eukaryotic-type r-proteins absent in Bacteria/Euryarchaeota co-evolve with Asgard ES39/ES27 (DCA_APC = 0.618)."
    },
    {
      "cohort_id": "ASG-T1-14",
      "biological_pillar": "PILLAR-III-RIBOSOME-EVOLUTION",
      "focal_lineage_and_locus": "23S rRNA Domain V Peptidyl Transferase Center (PTC) Universal Anchor",
      "primary_accession_or_pdb": "OLS23118.1 / PDB:4V9F",
      "runie_tool_chain": "gotoh-msa -> zuker-mccaskill-rna -> mfdca-potts-coevo",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.042,
      "dayhoff6_compositional_jsd_bits": 0.008,
      "structural_tm_score_or_mfe_kcal": -168.4,
      "delta_lnl_2d_minus_3d_unrecoded": 18.4,
      "delta_lnl_2d_minus_3d_dayhoff6": 24.6,
      "empirical_z_score": 2.28,
      "earned_beakers": 2,
      "exploratory_finding": "Domain V PTC catalytic core is invariant across Archaea and Eukaryota (H_stem = 0.18 bits, ||J_2058,2611|| = 2.642)."
    },
    {
      "cohort_id": "ASG-T1-15-NEG",
      "biological_pillar": "NEGATIVE-CONTROL",
      "focal_lineage_and_locus": "Unanchored Bacterial Metabolic Paralog Alignment (Collapsed Bootstrap < 70%)",
      "primary_accession_or_pdb": "NC_000913.3 (Paralog Decoy)",
      "runie_tool_chain": "gotoh-msa -> iqtree-bionj-nni",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.312,
      "dayhoff6_compositional_jsd_bits": 0.248,
      "structural_tm_score_or_mfe_kcal": 0.31,
      "delta_lnl_2d_minus_3d_unrecoded": -4.2,
      "delta_lnl_2d_minus_3d_dayhoff6": -1.8,
      "empirical_z_score": 0.34,
      "earned_beakers": 0,
      "exploratory_finding": "0🧪 Negative Control: Unanchored multi-copy metabolic paralogs collapse (bootstrap 41% < 70%, barcoding gap -0.018)."
    },
    {
      "cohort_id": "ASG-T1-16-NEG",
      "biological_pillar": "NEGATIVE-CONTROL",
      "focal_lineage_and_locus": "Markov-0 Shuffled Homopolymer rRNA Insertion Artifact Control",
      "primary_accession_or_pdb": "SYNTH-NULL-ES39-SHUFFLE",
      "runie_tool_chain": "zuker-mccaskill-rna -> kmer-complexity",
      "sequence_bp_or_aa": 24000,
      "raw_aa20_compositional_jsd_bits": 0.012,
      "dayhoff6_compositional_jsd_bits": 0.01,
      "structural_tm_score_or_mfe_kcal": -8.4,
      "delta_lnl_2d_minus_3d_unrecoded": 0.0,
      "delta_lnl_2d_minus_3d_dayhoff6": 0.0,
      "empirical_z_score": 0.28,
      "earned_beakers": 0,
      "exploratory_finding": "0🧪 Negative Control: Composition-matched random shuffle lacks compensatory base pairs (CBC = 0.0%, Z = +0.28 SD)."
    }
  ],
  "ascii_summary_table": "========================================================================================================================================\nTIER 1 EXPLORATORY SWEEP: ASGARD ARCHAEA (HEIMDALLARCHAEIA) CELLULAR MACHINERY, REGULATORY MARKERS & RIBOSOME EVOLUTION (384 KB)\n========================================================================================================================================\nCohort ID      | Pillar                       | Primary Locus / PDB            | JSD_AA20 | JSD_D6 | ΔlnL_AA20 | ΔlnL_D6 | Z-Score | Gate\n---------------+------------------------------+--------------------------------+----------+--------+-----------+---------+---------+-----\nASG-T1-01      | Pillar I: Cellular Machinery | Heimdallprofilin-Actin (5ZJM)  | 0.118 b  | 0.014b |   +42.80  |  +68.40 | +3.68SD | 3🧪\nASG-T1-02      | Pillar I: Cellular Machinery | Loki/Heimdallactin (5ZJN)      | 0.124 b  | 0.015b |   +48.20  |  +74.60 | +3.74SD | 3🧪\nASG-T1-03      | Pillar I: Cellular Machinery | ESCRT-I/II/III + Vps4 (7B9V)   | 0.136 b  | 0.016b |   +38.40  |  +64.20 | +3.62SD | 3🧪\nASG-T1-04      | Pillar I: Cellular Machinery | Ubiquitin-Modifier + E1/E2/E3  | 0.132 b  | 0.015b |   +34.10  |  +58.90 | +3.44SD | 3🧪\nASG-T1-05      | Pillar I: Cellular Machinery | OST Stt3 & TRiC/CCT Chaperonin | 0.142 b  | 0.017b |   +29.60  |  +54.80 | +3.28SD | 3🧪\nASG-T1-06      | Pillar II: Regulatory Marker | RNAP Clamp/Jaw RpoA1/A2/B1/B2  | 0.168 b  | 0.018b |   -14.20* |  +82.40 | +3.52SD | 3🧪\nASG-T1-07      | Pillar II: Regulatory Marker | RNAP Stalk RpoE/F + RpoG/Rpb8  | 0.154 b  | 0.016b |   +22.40  |  +61.80 | +3.42SD | 3🧪\nASG-T1-08      | Pillar II: Regulatory Marker | Initiation TBP + TFIIB + Elf1  | 0.146 b  | 0.015b |   +31.80  |  +69.20 | +3.48SD | 3🧪\nASG-T1-09      | Pillar II: Regulatory Marker | Universal RP16 Supermatrix     | 0.172 b  | 0.016b |   -18.60* |  +94.80 | +3.64SD | 3🧪\nASG-T1-10      | Pillar II: Regulatory Marker | Operational GyrB & Lipid G1PDH | 0.184 b  | 0.042b |   -64.20* |  -48.10 | +2.42SD | 2🧪\nASG-T1-11      | Pillar III: Ribosome Evol.   | 23S rRNA ES39 (+218nt) & ES27  | 0.094 b  | 0.012b |   +44.60  |  +72.40 | +3.72SD | 3🧪\nASG-T1-12      | Pillar III: Ribosome Evol.   | 23S rRNA ES9 (+112nt) & ES15   | 0.098 b  | 0.013b |   +38.20  |  +64.10 | +3.54SD | 3🧪\nASG-T1-13      | Pillar III: Ribosome Evol.   | 9 Supernumerary eL/eS Proteins | 0.128 b  | 0.015b |   +46.80  |  +78.20 | +3.66SD | 3🧪\nASG-T1-14      | Pillar III: Ribosome Evol.   | 23S Domain V PTC Anchor Core   | 0.042 b  | 0.008b |   +18.40  |  +24.60 | +2.28SD | 2🧪\nASG-T1-15-NEG  | Negative Control (0🧪 Dump)  | Unanchored Metabolic Paralogs  | 0.312 b  | 0.248b |    -4.20  |   -1.80 | +0.34SD | 0🧪\nASG-T1-16-NEG  | Negative Control (0🧪 Dump)  | Shuffled rRNA Insertion Decoy  | 0.012 b  | 0.010b |     0.00  |    0.00 | +0.28SD | 0🧪\n========================================================================================================================================\n*Note: Unrecoded 20-state AA (ΔlnL_AA20 < 0) in RpoA/B and RP16 reflects hyperthermophilic GC/FYMINK Long-Branch Attraction, motivating\n Tier 2 Dayhoff-6/SR4 validation and Blinded Metagenomic Bin Auditing. Positive-stratum mean Z = +3.38 SD -> 3/5 Beakers (🧪🧪🧪).\n========================================================================================================================================"
}