Gene Essentiality Across Conditions and Predictors

Gene essentiality is dependence of growth on a gene’s function under a defined condition. This synthesis integrates transposon sequencing (TnSeq, which infers importance from mutant abundance), flux balance analysis (FBA, a constraint-based metabolic-modeling method), deletion phenotypes, annotation, pangenome conservation, orthology, pathway completeness, fitness modules, functional-dark-matter prioritization, metabolic capability versus dependency, genotype–condition prediction, metal fitness, respiratory wiring, and ecological context. Direct perturbation and measured growth provide the strongest evidence; conservation, annotation, modeling, and ecological associations are predictors or hypotheses rather than interchangeable definitions of essentiality. ₁

The Acinetobacter baylyi ADP1 Data Explorer compares tnseq, flux-balance-analysis, annotation, and pangenome status in Acinetobacter baylyi ADP1. ADP1 Deletion Collection Phenotype Analysis adds a 2,034-gene deletion analysis across 8 carbon sources, while Triple Essentiality Concordance Analysis for *Acinetobacter baylyi* ADP1 compares up to 5,852 genes across FBA, RB-TnSeq (random-barcode transposon sequencing), complete-gene knockout phenotypes, mutant growth, and proteomics. _{1 2 3} The wider corpus adds annotation-gap, metabolic, ecological, metal, module, antimicrobial-resistance, and dark-gene analyses. _{4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22} AMR Co-Fitness Support Networks and Fitness Cost of Antimicrobial Resistance Genes extend condition-specific interpretation across 28 and 25 organisms. _{23 24}

Condition-Specific Essentiality and Fitness

Essentiality changed with medium: 499 genes were essential on minimal media versus 346 on LB, and FBA flux classes changed between rich and minimal media for 177 of 866 genes, or 20%. ₁ The KBase Data Lakehouse atlas inventories 27,410,721 FitnessBrowser measurements and genes tested under an average of 121–187 conditions, although this inventory does not itself establish essentiality. ₆

The ADP1 deletion collection supports condition dependence: across 8 carbon sources, mean growth ratios ranged from 0.41 on urea to 1.36 on quinate; 625 genes, or 31% of the complete matrix, had condition-specificity scores ≥ 1.0. The urease complex was required on urea, aromatic degradation and pqq-biosynthesis on quinate, and the glyoxylate-shunt on acetate. Five principal components (PCA, a method summarizing correlated variation into orthogonal components) captured 82% of variance; the highest pairwise condition correlation was r = 0.58 and the median was r = 0.25. ₂

The triple-essentiality analysis supports the same pattern: across 8 carbon sources, 333 of 478 genes (70%) showed condition-specific defects, 10 (2%) were defective across all 8, and 135 (28%) showed no defect; mean pairwise defect correlation was 0.38, ranging from -0.03 to 1.0. The Q25 threshold classified 72% of genes as defective; defect rates were 73%, 74%, and 69% in the FBA-essential, FBA-variable, and FBA-blocked groups. The expected rate under independent conditions was 1 − 0.75⁸ = 90%, but aggregation depends on inter-condition correlation. ₃

Condition-specific effects are not predominantly accessory: strongly condition-specific genes were 77.3% core versus 70.3% for genes without such phenotypes (OR=1.78, p=1.8e-97), and 4,450 genes (2.7%) were neutral overall but critical in one condition. These ephemeral-niche genes were 3.0% of core genes, 1.7% of auxiliary genes, and 1.6% of singleton genes. This supports the hypothesis that core genes can be conditionally important but does not establish that all core genes are broadly adaptive. ₁₅

Among 757 experiments, 337 (44.5%) were field and 420 (55.5%) lab. Among 2,725 non-essential genes with fitness and pangenome links, genes important at fitness < -2 were 83.6% core for field-stress and 82.4% for field-core, versus 81.4% for lab-nutrient, 81.5% for lab-other, 73.4% for lab-antibiotic, and 71.2% for heavy-metals; field-stress, field-core, and lab-nutrient enrichment was significant after BH-FDR (Benjamini–Hochberg false-discovery-rate) correction. ₁₄

Capability–dependency analyses support condition-specific importance but refine capacity–dependence claims. One analysis identified 1,695 complete pathway-organism pairs: 267 (15.8%) latent, 547 (32.3%) intermediate, and 881 (51.9%) active; pathway category predicted class strongly (χ²=163.6, df=4, p=2.5×10⁻³⁴). Carbon-source pathways were latent in 217 of 892 cases (24.3%), versus 48 of 735 amino-acid-biosynthesis cases (6.5%), making carbon pathways 3.7× more likely to be latent. Latent fractions ranged from 0–31.6% across organisms, averaged 15.0%, and ranged from 4.7% to 21.1% across 16 threshold combinations (SD = 5.9 percentage points). ₁₉

A stricter four-way analysis of 161 organism–pathway pairs across 7 model bacteria and 23 GapMind pathways found 57 (35.4%) Active Dependency, 66 (41.0%) Latent Capability, 24 (14.9%) Incomplete but Important, and 14 (8.7%) Missing. Its composite importance score weighted essentiality at 40%, fitness breadth at 30%, and fitness magnitude at 30%. Mapping used Fitness Browser-native KEGG annotations through besthitkegg-to-keggmember-to-EC-to-KEGG-map-to-GapMind pathway links. All 66 Latent Capability pairs became fitness-important under at least one condition type, with nitrogen limitation, stress, and carbon limitation most frequent, although median-based subset thresholds can reclassify pairs by construction. ₂₅

The Web of Microbes bridge supports using metabolic context to test condition-specific fitness but qualifies production-to-dependency inference. For Pseudomonas FW300-N2E3, 19 produced metabolites matched Fitness Browser carbon- or nitrogen-source experiments; 5 were de novo products and 14 amplified metabolites. ₂₆ The archived snapshot recorded production without organism consumption: controls had D for detected in starting medium (742 observations) and N for not detected (1,023); organism observations had I for increased (1,338), E for emerged de novo (1,155), and N for no significant change (7,509). E and I were mutually exclusive across 10,744 observations, and no organism had a decrease or consumption action. ₂₆ Among ENIGMA isolates in R2A, E/(E+I) ranged from 15.2% to 32.4%, from 16 of 55 changes in Bacillus FW507-8R2A to 34 of 83 in Pseudomonas GW456-L13; this is a proposed metabolic-novelty phenotype, not an established essentiality predictor. ₂₆

Prediction Targets and Evaluation Boundaries

Essentiality analysis does not evaluate one universal task: lethal-versus-dispensable classification, quantitative growth impairment among viable mutants, continuous fitness ranking, and condition-specific defects have different targets and evaluation boundaries. The ADP1 triple-essentiality analysis demonstrates that changing the target changes the apparent performance of FBA, TnSeq, knockout assays, and proteomics. ₃ The report treats complete-gene knockout experiments as the best proxy for lethality, while fitness, expression, and growth rates reveal additional biology, so integration of multiple methods gives the most complete picture rather than any single assay standing in for essentiality. ₃

For lethal-versus-dispensable classification, refined FBA-versus-complete-knockout comparisons covered 724 genes in rich media, with recall = 60.8%, precision = 64.0%, specificity = 79.7%, F1 = 0.624, and Cohen’s kappa = 0.486; the minimal-media comparison covered 833 genes, with recall = 65.6%, precision = 69.2%, specificity = 78.9%, F1 = 0.673, and Cohen’s kappa = 0.493. Confusion-matrix counts were True Positives = 164, False Positives = 92, False Negatives = 106, and True Negatives = 362 in rich media, and True Positives = 229, False Positives = 102, False Negatives = 120, and True Negatives = 382 in minimal media. ₃ FBA performed better in minimal media (F1=0.67) than rich media (F1=0.62), which the report attributes tentatively to tighter metabolic constraints in minimal media rather than to a demonstrated mechanism. ₃ These results support moderate FBA concordance at the lethal boundary but do not generalize to quantitative growth defects among TnSeq-dispensable genes, where FBA class was not associated with defect status (p=0.63). ₃ On that basis the report recommends FBA only as a first-pass screen, citing ~60-70% accuracy (κ≈0.49), better performance in minimal than rich media, and experimental validation of all predictions. ₃

The quantitative-growth analysis was restricted to 478 genes that were TnSeq-dispensable and had FBA and mutant-growth data. FBA class was not associated with growth-defect status: χ² = 0.93, p = 0.63, 2 df; defect rates were 73.1% (57/78) for FBA-essential, 73.5% (150/204) for FBA-variable, and 69.4% (136/196) for FBA-blocked genes; the Kruskal–Wallis test on mean growth rates was H = 1.67, p = 0.43. The null association persisted at Q10, Q15, Q20, Q25, Q30, and Q35, whose chi-squared p-values were 0.985, 0.887, 0.405, 0.629, 0.313, and 0.058; only at Q40 did a marginal association appear (p = 0.048), driven by FBA-essential versus FBA-blocked defect rates of 90% vs 81%, while the threshold-independent Kruskal–Wallis test on continuous mean growth rates remained non-significant (H = 1.67, p = 0.43). ₃ The reported threshold-by-threshold values for percent with defect, FBA-essential rate, FBA-variable rate, and FBA-blocked rate were: Q10, 41%, 41%, 42%, 41%; Q15, 56%, 55%, 57%, 55%; Q20, 65%, 65%, 68%, 62%; Q25, 72%, 73%, 74%, 69%; Q30, 77%, 79%, 79%, 73%; Q35, 82%, 85%, 85%, 77%; and Q40, 85%, 90%, 88%, 81%. ₃ Within the FBA-essential set, 57/78 (73%) showed growth defects — matching the 72% background rate among all 478 genes — and 21/78 (27%) grew normally, which the report reads as FBA over-prediction; among FBA-blocked genes with zero predicted flux, 136/196 (69%) showed defects, read as FBA under-prediction, and 60/196 (31%) grew normally, consistent with FBA. ₃

Continuous fitness retained ranking information better than binary essentiality fraction for knockout-lethality prediction: inverted fitness produced AUC (Area Under the ROC Curve, a threshold-independent classification measure) = 0.700 in rich media and 0.725 in minimal media, versus essentiality-fraction AUC = 0.344 across 1,933 rich-media genes and 0.403 across 2,011 minimal-media genes, both worse than random; minimal-media log2 proteomics reached AUC = 0.743. Among dispensable genes, growth rates vary from 10% to 100% of wild-type, so a binary label discards most of that range. This refines interpretation of condition-specific-fitness because reducing continuous outcomes to essential-versus-dispensable labels loses information. ₃

RB-TnSeq and complete-gene knockout assays disagreed strongly when both were treated as binary essentiality tests in rich media. Across RB-TnSeq thresholds of 0.01, 0.025, 0.05, 0.10, and 0.20, Cohen’s kappa was -0.139, -0.122, -0.081, -0.024, and -0.014. At the 0.05 threshold, 1,933 genes were compared: recall = 7.9%, precision = 5.8%, specificity = 82.8%, F1 = 0.067, and kappa = -0.081; both-essential numbered 18 genes (0.9%), both-dispensable 1,411 (73.0%), knockout-essential/TnSeq-dispensable 211 (10.9%), and knockout-dispensable/TnSeq-essential 293 (15.2%). ₃ The report’s own summary of this result is systematic TnSeq–knockout discordance (κ<0) with fitness a better predictor (AUC=0.70-0.73) than essentiality_fraction (AUC=0.34-0.40); these are the stated ranges rather than recomputed values. ₃

The report interprets this discordance as consistent with transposon insertion measuring fitness cost while potentially preserving partial function, whereas complete deletion measures lethality. Partial or truncated protein production, read-through transcription, retained functional domains, and condition aggregation are proposed mechanisms, not resolved demonstrations. The knockout ground truth was the adp1-single-gene-deletion-collection of de Berardinis et al. (2008), which supplied essentiality calls for rich and minimal media and mutant growth rates on 8 carbon sources. ₃ Rich-media (LB) knockout data covered 2,953 genes, including 346 essential and 2,607 dispensable genes; the merged minimal-media set covered 3,092 genes, including 499 essential genes, with minimal-media calls prioritized and rich-media calls used as fallback. ₃

Proteomics supplied an independent continuous evidence axis rather than a condition-matched causal test. Across 2,288 genes measured over 7 Acinetobacter strains, essential genes (n = 464) had mean log2 expression = 28.43 ± 2.94 versus 25.73 ± 2.96 for dispensable genes (n = 1,824), a difference of 2.70 log2 units or 6.5-fold; Mann–Whitney U p = 9.91×10⁻⁵⁹. Expression correlated with knockout essentiality at Pearson r = 0.345 (p = 5.32×10⁻⁶⁵) and Spearman ρ = 0.338 (p = 3.28×10⁻⁶²), with ROC AUC = 0.743, the best overall continuous predictor in this comparison and comparable to fitness. In minimal media the report pairs proteomics (AUC=0.74) with FBA (F1=0.67) as complementary measures, expression validating essentiality while FBA predicts metabolic requirement. This proteomics analysis was not part of the original study and provides independent validation of essentiality through expression levels; averaging across 7 strains rather than the exact knockout condition makes it supporting evidence, not condition-matched causation, and complements multi-omics-integration. ₃

Fitness Importance, Conservation, and Accessory Essentiality

The conservation result is a weak but reproducible gradient. Across approximately 194,000 genes from 43 bacteria, essential genes were 82% core (n=27,693), genes often sick in >10% of experiments 78% core (n=15,989), mixed genes 70% core (n=20,739), sometimes-sick genes 72% core (n=25,201), always-neutral genes 66% core (n=94,889), and sometimes-beneficial genes 70% core (n=9,705). By strongest measured effect, essential genes were 82.2% core, genes with min_fit < -3 were 77.7% core, and genes with min_fit from -1 to 0 were 66.4% core. Fitness breadth had a positive but weak association with conservation (Spearman rho=0.086, p=8.1e-230): essential genes were 82% core, genes affecting 20+ experiments 79%, genes affecting 6–20 experiments 73%, genes affecting 1–5 experiments 71%, and genes with 0 experiments 66%. ₁₅ The central digest records the same cohort as 194,216 protein-coding genes in 43 bacteria and reads the 16-percentage-point gradient from essential genes at 82% core to always-neutral genes at 66% core — reproduced when binning by strongest negative fitness effect (min_fit < -3 → 78% core vs min_fit near 0 → 66%) and by fitness breadth (important in 20+ conditions → 79% core vs 0 conditions → 66%) — as establishing that the essentiality–conservation link is quantitative rather than binary. ₂₇

These values support a relationship between functional importance, breadth of measured effects, and conservation, but the 16-percentage-point difference between 82% core among essential genes and 66% among always-neutral genes limits conservation as a standalone predictor. ₁₅

A separate linkage analysis supports the direction while using different filters: it identified 27,693 putative essential genes, 18.6% of 148,826 protein-coding genes across 33 organisms, with a per-organism range of 12.9-28.9%. Essential genes were 86.1% core versus 81.2% for non-essential genes, with a median odds ratio of 1.56 for enrichment in the core genome — an odds ratio, which is not a probability ratio and should not be read as a likelihood multiplier; 18 of 33 organisms showed significant enrichment by Fisher’s exact test at BH-FDR q < 0.05. Essentiality was inferred from RB-TnSeq under the library-construction and growth conditions represented in the Fitness Browser. ₂₈ The central digest reports the same 86.1% versus 81.2% contrast with median OR=1.56 and 18 of 33 organisms significant after BH-FDR, and reads the enrichment as real but modest because most genes in well-characterized bacteria are core regardless of essentiality, with the signal strongest in organisms with larger clades, where more genomes make core classification more reliable. ₂₇ Of 44 of 48 Fitness Browser organisms mapped to pangenome species clades, the analysis produced 177,863 gene-to-cluster links with 100.0% median protein identity and 94.2% median gene coverage. Thirty-four organisms had at least 90% coverage, and 33 entered downstream analysis because Dyella79 was excluded for a locus-tag mismatch. Four organisms—Cola, Kang, Magneto, and SB2B—were unmatched because their species had too few genomes in GTDB for pangenome construction. ₂₈

Of the linked genes, 145,821 (82.0%) were core and 32,042 (18.0%) auxiliary; 7,574 singleton genes were a subset of the auxiliary category. Strongest enrichment signals were Methanococcus maripaludis S2 (OR=5.21), Ralstonia syzygii PSI07 (OR=3.41), and Marinobacter adhaerens (OR=3.08). ₂₈

Essential genes can therefore be strain-specific, auxiliary, or highly divergent. The essential-auxiliary category contained 3,683 genes: 13.4% were enzymes and 38.2% hypothetical, with leading subsystems including ribosomes, DNA replication, type 4 secretion, and plasmid replication. These profiles support strain-specific variants of conserved machinery or mobile-element-associated functions, but do not demonstrate compensation for missing core functions. ₂₈ As secondary literature context rather than new measurement, the report aligns its essential-core, essential-auxiliary, and essential-unmapped categories with Rosconi et al. (2022), who studied essentiality across 36 streptococcus-pneumoniae strains and concluded that the pan-genome makes gene essentiality strain-dependent, distinguishing universal essential, core strain-specific essential, and accessory essential genes. ₂₈

The broader inventory contains 1,011,650,903 KBase pangenome gene rows spanning 293,059 genomes, 27,690 GTDB species clades or pangenomes, and 132.5M gene clusters. A separate 48-organism analysis found 41,059 of 221,005 genes essential (18.6%), with organism-level rates from 12.2% in Pedo557 to 29.7% in Magneto. Across 17,222 ortholog families, 859 (5.0%) were universally essential, 4,799 (27.9%) variably essential, and 11,564 (67.1%) never essential. _{6 27 12} Antimicrobial-resistance genes fall rarely into the putatively essential class: only 4.6% of AMR genes were absent from the fitness matrices, substantially lower than the ~14% background essential rate estimated in a prior analysis on a different organism set, so the contrast is across datasets rather than within one. ₂₄

Universally essential genes were 91.7% core, variably essential genes 88.9%, never-essential genes 81.7%, and orphan essential genes 49.5%; counts were 6,963, 90,844, 47,140, and 4,683. Of 17,222 ortholog families, 839 of 859 universally essential families were strict single-copy families with copy ratio ≤1.5 and no non-essential paralogs. Essentiality penetrance and core fraction had weak positive correlation (rho=0.123, p=1.6e-17). ₁₂

Fifteen named ortholog families were essential in every one of the 48 organisms with no exceptions: ribosomal proteins rpsC, rplW, rplK, rplB, rplA, rplF, rps11, rpsJ, rpsI, and rpsM; chaperonin groEL; CTP synthase pyrG; translation elongation factor G fusA; valyl-tRNA synthetase valS; and geranyltranstransferase SelGGPS — described as the irreducible essential core of bacterial life. _{12 27} The central digest labels that cohort as 48 bacteria, whereas the source report states the panel spans Proteobacteria, Bacteroidetes, Firmicutes, and Archaea and names Methanococcus among its organisms, so the bacteria-only phrasing understates the taxonomic span of the same analysis. _{27 12}

How “universal” should be read is fixed by the classification’s own table definition: a family counts as universally essential when it is essential in every organism where the family has members, not when it is present in all 48. ₁₂ The same report’s family classification of 17,222 ortholog families across 48 bacteria gives 859 (5.0%) universally essential, 4,799 (27.9%) variably essential, and 11,564 (67.1%) never essential, with 839 of the universally essential families strict single-copy (copy ratio ≤1.5, no non-essential paralogs) and 20 containing paralogs. ₁₂ Variably essential families had a median essentiality penetrance of 33%, with 813 families more than 50% essential and 704 families less than 10% essential, so most essential gene families are essential in only a minority of organisms. ₁₂

The report reads variable essentiality as the dominant mode rather than the exception, and extends Rosconi et al. (2022) — who demonstrated across 36 streptococcus-pneumoniae strains that the pan-genome makes essentiality “strain-dependent and evolvable” — from within-species to cross-species variation across 48 phylogenetically diverse species spanning Proteobacteria, Bacteroidetes, Firmicutes, and Archaea, offering as an example that the same gene can be essential in Methanococcus but dispensable in Pseudomonas depending on metabolic context. ₁₂ The mechanisms it proposes for variable essentiality — presence of paralogs, alternative pathways, or compensatory functions in the accessory genome — are explanations rather than tested causes, and its inference that only the 859 universally essential families are reliable broad-spectrum antibiotic targets is an interpretation, not a drug-efficacy result. ₁₂

Universally essential genes were 91.7% core versus 80.7% for non-essential genes, 71% of universally essential families were 100% core across all genomes in their species, and orphan essentials were only 49.5% core, marking them as strain-specific essential functions. _{12 27} The 80.7% non-essential figure and the 81.7% never-essential figure above come from different gene groupings and must not be merged. ₁₂ For the orphan essentials, recent acquisition, rapid evolution, and lineage-specific innovation are proposed rather than demonstrated explanations; their low core rate suggests many are not conserved even within their own species, and the report compares them with the 149 of 473 genes (31%) of unknown function in the minimal Mycoplasma genome of Hutchison et al. (2016), concluding that this dark matter of essential genes is larger across diverse bacteria. ₁₂

The report also describes its universally essential set as smaller than computationally conserved sets because it is restricted to genes where transposon disruption is lethal across all 48 tested organisms, which it calls more stringent than computational conservation. ₁₂ That description does not match the same report’s table definition of universally essential as essential in every organism where the family has members, so the stringency criterion is an unresolved inconsistency inside the source rather than an established threshold; and because essentiality here is called from absent transposon mutants, such calls supply only an upper bound on true essentiality. _{12 29} On the stringency claim the report positions itself as complementary to Duffield et al. (2010), who identified 52 conserved essential proteins across 14 bacterial genomes as drug targets by computational prediction, by supplying experimental RB-TnSeq essentiality across 48 organisms. ₁₂ Cross-referencing the Database of Essential Genes to identify families essential in non-Fitness Browser organisms is a proposed comparison, not an analysis this report performed. ₁₂

Conservation, Annotation, and Functional Novelty

Essential genes were more annotation-rich: 33% had COG assignments versus 5% of dispensable genes, and 92% had KEGG KO assignments versus 53%; approximately 8% lacked KO assignments and were candidate novel essential functions. In ADP1, essential genes were also more likely to be in the core pangenome, consistent with the pattern that conserved genes tend to be essential. ₁ AlphaFold analysis found median MSA depth 15,308 for core clusters versus 5,299 for auxiliary+singleton and 5,527 for auxiliary non-singleton clusters; across 38,051,842 bridged pairs, MSA depth and domain-hit count had Spearman ρ=0.7563. Yet 415,603 core clusters had MSA depth <10 and 68.9% were hypothetical, and only 29.3% of 132,531,501 clusters bridged to AlphaFold. ₃₀

Direct linkage refines the annotation gradient: essential-core contained 22,751 genes and was the most enzyme-rich and best-annotated class at 41.9% enzymes, 87% with known function, and 13.0% hypothetical; essential-auxiliary contained 3,683 genes and was 13.4% enzymes and 38.2% hypothetical; essential-unmapped contained 1,259 strain-specific essentials with no pangenome cluster match and was 18.2% enzymes and 44.7% hypothetical, the least characterized class. Essential-core was enriched in Protein Metabolism by +13.7 percentage points, Cofactors/Vitamins by +6.2%, Cell Wall by +3.9%, and Fatty Acid biosynthesis by +3.1%, and depleted in Carbohydrates by -7.9%, Amino Acids by -5.6%, and Membrane Transport by -4.0%, relative to non-essential genes; the report reads the depleted categories as functions that tend to be conditionally important rather than universally essential, which is its inference rather than a measured essentiality call. Essential-auxiliary genes included ribosomes, DNA replication, type 4 secretion, and plasmid-replication functions; essential-unmapped genes included divergent ribosomal proteins L34, L36, S11, and S12, translation factors, transposases, and DNA-binding proteins, which the report proposes are likely recently acquired or highly divergent variants of core functions — an interpretation, not a demonstrated origin. ₂₈ The central digest reports the same three-class contrast — essential-core at 41.9% enzymes and only 13.0% hypothetical, the well-characterized metabolic backbone of ribosomes, DNA replication, cell wall, and cofactor biosynthesis; essential-auxiliary at 13.4% enzymes and 38.2% hypothetical; and essential-unmapped at 44.7% hypothetical, which it names prime targets for functional discovery. ₂₇ The digest states the functional-category differences explicitly in percentage points — Protein Metabolism +13.7 pp, Cofactors/Vitamins +6.2 pp, and Cell Wall +3.9 pp enrichment, and Carbohydrates -7.9 pp, Amino Acids -5.6 pp, and Membrane Transport -4.0 pp depletion — and likewise reads the depleted categories as conditionally important rather than universally essential. ₂₇ The project report writes those same differences with a percent sign (for example Cofactors/Vitamins (+6.2%), Carbohydrates (-7.9%)), so the two notations refer to the one set of percentage-point differences between class frequencies and must not be read as interchangeable with percentages. _{28 27}

A poorly conserved essential class remains: 7,084 essential genes had no detectable ortholog in any other Fitness Browser organism, and 58.7% were hypothetical versus 8.2% of universally essential genes, which both the project report and the central digest read as the least characterized yet most functionally important genes in each organism and prime candidates for novel biology. _{12 27} Hypothetical essential genes numbered 8,297, or 20.2% of all essentials — a set that includes these hypothetical orphans rather than standing apart from them. Of these, 3,912 had orthologs and were predictable by module transfer, whereas 4,385 orphan genes were not. A further analysis transferred 1,382 function predictions to hypothetical essential genes from fitness-module context in non-essential orthologs. ₁₂

The functional-dark-matter analysis found 57,011 dark genes (24.9% of 228,709 genes), including 9,557 essential dark genes; 39,532 (69.3%) had pangenome links, 12,686 were accessory, and 511 were both accessory and strongly fitness-active. ₁₇ The truly-dark analysis sharpens this result: of 57,011 dark genes, 33,105 were reclassified by Bakta, 6,427 remained hypothetical in both pipelines, and 17,479 lacked pangenome links. Among strongly phenotyped genes, stress conditions were reported at 28.7% versus 43.2% for annotation-lag genes (OR = 0.53, p < 0.001); the results section reported stress at 43.3% versus 54.7%, carbon-source conditions at 13.7% versus 21.7%, and motility at 1.2% versus 3.2%. ₃₁

Among 39,532 pangenome-linked dark genes, persistent hypothetical genes had an essential fraction of 18.0% versus 13.4% for annotation-lag genes, core fraction of 43.1% versus 72.7%, ortholog presence of 29.3% versus 63.7%, and median ortholog breadth of 1 versus 4 organisms. Median length was 121 versus 194 amino acids and mean GC content was 0.542 versus 0.584; length d = −0.432, core OR = 0.284, essential OR = 1.420, GC d = −0.395, ortholog-presence OR = 0.236, and ortholog-breadth d = −1.072, with corresponding p-values < 1e-100, < 1e-100, < 1e-10, 2e-115, < 1e-100, and < 1e-100. Thus the aggregate 49.5% core fraction for orphan essentials is refined by a shorter, narrower-breadth persistent-dark subset, although short genes are also harder to annotate and measure. _{31 12}

Most persistent-dark genes still had sequence-level clues: 79.4% had UniRef50 links, 84.7% database cross-references, 4.0% Pfam hits, and 4.6% KEGG KOs; eggNOG-mapper supplied partial signal for 43.5%, while 17.0% had COG categories and 55.4% of those assignments were category S, function unknown. Only 246 genes (3.8%) had no clues; 3,867 (60.2%) had minimal sequence identifiers, 711 (11.1%) had partial Pfam, COG, or eggNOG function, and 1,603 (24.9%) had phenotype or module evidence without functional annotation. ₃₁

Essential dark genes are invisible to standard RB-TnSeq: the 9,557 essential dark genes have no viable transposon mutants and therefore zero rows in the genefitness table. ₁₇ Of 17,344 scored dark genes, 9,557 (55%) were essential, making them the majority of phenotype-bearing dark matter even though they score poorly under the fitness-centric framework and were absent from its top 100; they were nonetheless ranked separately, so their low fitness-based scores reflect the scoring framework rather than lesser importance. ₁₇

That separate prioritization used five evidence dimensions that do not require fitness magnitudes — gene neighbor context (0.25), cross-organism conservation (0.20), phylogenetic breadth (0.20), domain annotations (0.15), and crispri (CRISPR interference, transcriptional knockdown) tractability (0.20) — and ranked all 9,557 essential dark genes, producing 50 ranked candidates amenable to CRISPRi knockdown experiments. ₁₇

The three highest-ranked essential candidates each carry candidate-specific supporting evidence: escherichia-coli K-12 Keio locus 14796, score 0.875, a YbeY domain, conserved in 5+ phyla with a synteny-confirmed operon and predicted ion transport; shewanella-oneidensis MR-1 locus 200382, score 0.874, RimP_N/DUF150_C domains, a predicted ribosome assembly factor in a conserved neighborhood; and klebsiella-oxytoca BWI76_RS08540, score 0.865, OmpA/TIGR02802 domains, predicted cell division and conserved across Enterobacteriaceae. All three functions are predictions, not established assignments. ₁₇

All top-50 candidates carry high-confidence functional hypotheses derived from operon context, each paired with a specific CRISPRi experiment design naming target organism, sgRNA target, growth condition, expected phenotype, and validation strategy — designs, not validated functional assignments. ₁₇ The ranking’s inputs are clues rather than confirmation: candidates score highest where detectable domains providing structural clues, membership in widely conserved ortholog groups confirming importance beyond one organism, conserved gene neighborhoods with annotated partners providing functional context, and residence in CRISPRi-amenable organisms converge, making them experimentally testable despite being essential. ₁₇

The conservation evidence behind these rankings is a lower bound: synteny was assessed only in a five-gene window within the 48 Fitness Browser genomes, whereas tools such as STRING v12 and EFI-GNT analyze thousands of genomes with gene-fusion and shared-phylogenetic-profile scoring, so true conservation rates are likely higher under broader taxonomic sampling. ₁₇ The co-fitness signal is unavailable for essential genes altogether, because no viable mutants means no fitness profiles and hence no co-fitness data. ₁₇

Experimental recommendations are therefore stratified by gene class: CRISPRi knockdown for essential genes in tractable organisms, condition-specific TnSeq for fitness-active genes, and broad phenotypic screens for true knowledge gaps, with top conditions derived from fitness data or keyword inference. ₁₇ The 8,900 essential genes in the covering set cannot be studied by transposon knockout at all, which is why the recommended approach for them is CRISPRi knockdown — Mobile-CRISPRi for non-model organisms, following Peters et al. 2019 — with growth curves under standard and stress conditions; the covering set assigns those genes to specific organisms, the highest-priority targets being Keio:14796, MR-1:200382, and Koxy:BWI76_RS08540. ₁₇

Two figures document this prioritization: fig19_essential_scores.png shows essential dark gene score distributions across 5 dimensions, and fig20_essential_top20.png shows the top 20 essential dark gene candidates with score breakdown. ₁₇

The module-conservation result supports core enrichment at the module level: across 1,116 ICA modules in 32 organisms, module genes were 86.0% core versus 81.5% for all genes (OR=1.46, p=1.6e-87); among 974 modules with at least 3 mapped genes, 577 (59%) were >90% core, 349 (36%) mixed, and 48 (5%) <50% core. Module-family breadth did not predict conservation (Spearman rho=-0.01, p=0.914). Thirty-eight accessory module families had <50% core genes; horizontal transfer or niche-specific operons remain hypotheses. ₂₂

Predictor Agreement, Metabolic Capacity, and Ecological Context

Among 866 genes with FBA and TnSeq calls, 639 (73.8%) were concordant and 227 discordant; the report reads the discordant 26% as candidates for metabolic-model refinement or as reflecting regulatory effects not captured by FBA, and proposes testing whether those 227 genes are enriched for specific pathways or regulatory functions. ₁ Among 478 TnSeq-dispensable genes, FBA class was not associated with quantitative growth-defect status (χ²=0.93, p=0.63, 2 df; Kruskal–Wallis H=1.67, p=0.43); refined FBA-versus-knockout comparisons gave rich-media F1=0.624 and kappa=0.486 across 724 genes, and minimal-media F1=0.673 and kappa=0.493 across 833. _{1 3}

RB-TnSeq thresholds contradicted knockout essentiality in rich media: at 0.05, recall was 7.9%, precision 5.8%, specificity 82.8%, F1 0.067, and kappa -0.081 across 1,933 genes; kappa ranged from -0.139 to -0.014 across thresholds 0.01–0.20. FBA missed complex-level essentiality: 10/13 Complex I subunits had quinate-specific defects despite 1.76× higher predicted Complex I flux on aromatic substrates (0.55 vs 0.31) and 0% predicted essentiality, while 30/51 network genes lacked FBA mappings. _{3 5}

The respiratory study refines this model–measurement discordance. Quinate required Complex I, acetate required Complex I plus cytochrome bo3, ACIAD3522, and additional components, lactate specifically required cytochrome bo3, glucose had no specifically required respiratory component, and urea was generally demanding. Complex I growth ratios were 0.37 on quinate, 1.44 on glucose, and 0.49 on acetate; ACIAD3522 was dispensable on quinate and glucose but had a growth ratio of 0.013 on acetate, where the report describes it as lethal. ₃₂

ADP1 has three parallel NADH dehydrogenases: 13-subunit proton-pumping Complex I, single-subunit non-proton-pumping ndh 2, and single-subunit ACIAD3522, an NADH-FMN oxidoreductase. NDH-2 is ACIAD_RS16420 with KO K03885; it is TnSeq-dispensable, absent from the deletion collection, standalone rather than operonic, and core. FBA predicts zero NDH-2 flux on standard carbon sources because growth optimization routes NADH through Complex I. Complex I, NDH-2, and ACIAD3522 had mean protein levels of 27.6 at the 66th percentile, 27.0 at the 59th percentile, and 26.2 at the 48th percentile, respectively, against a genome median of 26.4; the spread was 1.4 units. After retaining organisms with 1–2 NDH-2 hits and excluding those with more than 2, 5 of 14 organisms had validated NDH-2, and their mean Complex I aromatic deficit was −0.297 versus −0.156 without validated NDH-2 (p=0.52). ₃₂

Fitness modules were stable at 17–52 modules per organism across 32 bacteria; 94.2% had significantly elevated within-module cofitness, mean |r|=0.34 versus background 0.12, and mean genomic-adjacency enrichment was 22.7×. Ortholog transfer achieved 95.8% precision and 91.2% coverage for gene-level KO prediction, whereas Module-ICA had <1% KO-level precision but captured process-level co-regulation. With 20% of KEGG-annotated genes withheld, ortholog transfer achieved 95.8% strict precision, 91.2% coverage, and 0.934 F1; domain prediction achieved 29.1%, 66.6%, and 0.401; Module-ICA achieved <1% and 23.3%; and cofitness voting achieved <1% and 73.0%. _{27 16}

GapMind completeness and Fitness Browser measurements were compared through proxy mappings, so Active Dependency, Incomplete but Important, Latent Capability, and Missing are pathway-level evidence rather than direct gene-level essentiality calls. The metal atlas likewise found that a 1,286-KEGG-KO presence signature did not predict metal fitness. _{25 20}

GapMind found 17 of 18 amino-acid biosynthesis pathways complete in all 7 successfully mapped organisms; serine biosynthesis was complete in 6 of 7, with DvH the sole exception. The seven-organism sample contained 6 organisms with 18/18 pathways and DvH with 17/18 (94.4%), but only 7 of 45 organisms mapped successfully. ₁₃

The essential-metabolome analysis takes the universally essential gene families as input but records them as 859 families across 45 organisms, so its organism denominator differs from the 48-organism panel and the two must not be merged. _{13 12} Its own caveat is that the RB-TnSeq experiments were conducted in rich media, where nutrient supplementation can make biosynthetic genes appear non-essential, so the data cannot distinguish “essential for biosynthesis” from “essential for viability”. ₁₃ Whether complete pathways are actually essential for viability, which pathways are essential only under specific growth conditions, and whether pathway gaps cluster by phylogeny or ecology all remain unresolved. ₁₃ The larger analysis found 267 of 1,695 complete pathway-organism pairs latent, especially carbon pathways, while amino-acid biosynthesis was usually active. DvH’s apparent serine auxotrophy remains a hypothesis requiring testing. ₁₉

Across approximately 293,000 genomes, the newer pathway analysis assessed 80 GapMind pathways—18 amino-acid biosynthesis and 62 carbon-source utilization pathways. Core-only versus all-gene completeness differed by 0.146 for leucine and valine, 0.141 for arginine, 0.140 for lysine, and 0.140 for threonine biosynthesis; corresponding all-gene/core-only values were 0.614/0.468, 0.614/0.468, 0.613/0.472, 0.804/0.664, and 0.803/0.663. These results are consistent with accessory contributions and possible community sharing, but do not demonstrate metabolite exchange or Black Queen dynamics. ₂₅

The capability–conservation relationship is not simple. Among 755 active dependencies, mean pathway conservation was 0.829; among 508 intermediate pathways it was 0.907; and among 248 latent capabilities it was 0.869. The active-versus-latent test gave Mann–Whitney U p=0.94 and rank-biserial r=0.052. Across 2,810 GTDB species with at least 10 genomes, variable pathway count correlated with pangenome openness at raw Spearman rho=0.327, p=7.2e-71, and partial Spearman rho=0.530, p=2.83e-203 after controlling for genome count. The signal was positive in 13 of 18 genera and significant at p<0.05 in 5 of 18: Clostridium, Eubacterium, Mesorhizobium, Pseudomonas, and Streptomyces; strongest within-genus correlations were rho=0.534 for Clostridium and rho=0.506 for Eubacterium. _{19 25}

All 10 target species formed metabolic clusters with silhouette scores >0.2, ranging from 0.35 for PALSA-747 sp. to 0.89 for Salmonella enterica. Cluster membership was associated with isolation environment in Salmonella enterica (χ²=1570.2, df=25, p<0.0001) and Phenylobacterium sp. (χ²=12.2, df=1, p=0.0005), but not in Stutzerimonas, Alteromonas, Pelagibacter, or Prochlorococcus. Strongest pathway heterogeneity occurred in valine/leucine biosynthesis (0.60), tryptophan biosynthesis (0.51), and lysine/threonine biosynthesis (0.49). These associations lack explicit phylogenetic correction and may reflect metadata or annotation limitations. ₁₉

Among 225 species with at least 50 genomes and at least 3 variable pathways, binary completeness profiles across 80 pathways produced a median of 4 metabolic ecotypes per species and a maximum of 8 in Alistipes onderdonkii and Barnesiella intestinihominis. Ecotype count correlated with pangenome openness at raw Spearman rho=0.262, p=6.8e-05, and partial rho=0.322, p=8.0e-07 after controlling for genome count. Genome fluidity is a plausible explanation, but sampling bias remains. ₂₅

The ENIGMA genotype–phenotype analysis supports interaction-aware models but contradicts expectations of general genotype-to-strain-fitness predictability. Across 46,389 genome × condition pairs from 727 genomes and 363 conditions, binary-growth AUC was 0.620 under genus-blocked holdout; amino acids, nucleosides, and carbon sources performed better (AUC 0.775, 0.780, and 0.695) than metals, antibiotics, and nitrogen (0.605, 0.619, and 0.435). KO × condition interactions increased mean AUC from 0.620 to 0.653 and improved 80 of 106 held-out genera; 95 of 343 individually testable conditions achieved AUC > 0.75. Continuous µmax, lag, and max_A had negative held-out R². ₁₈

GapMind achieved AUC 0.646, accuracy 78.8%, and coverage 24.3% on 118 testable pairs; matched-condition Carbon Source Phenotypes transfer achieved AUC 0.800 and accuracy 76.8% at 23% coverage, with internal five-fold AUC 0.858. Approximately 76% of conditions lacked either coverage source, and performance fell to approximately AUC 0.63. Gene content therefore predicts some capabilities but does not establish within-strain essentiality. ₁₈

Tensions

The two conservation–essentiality integrations report different essential-gene core fractions: 82% (and 82.2% under strongest-effect grouping) across approximately 194,000 genes from 43 bacteria, versus 86.1% among 148,826 genes from 33 organisms. _{15 28} This methodological tension involves datasets, organism filters, pangenome mappings, and denominators; the estimates must not be averaged.

The pitfalls audit names the denominator explicitly: computing percent core over all genes, including unmapped ones, gives essential = 82% core, whereas computing it over mapped genes only gives 86% core, and the same inconsistency once produced figures showing 86%/78% while the prose said 82%/66%. This identifies a likely source of the gap without settling which cohort estimate is correct, so the two analyses must be recomputed under a single documented denominator rather than averaged. ₂₉

ADP1 analyses associate essentiality with core status and richer annotation, whereas 415,603 core clusters had MSA depth <10 and 68.9% were hypothetical. Core enrichment also coexists with 66% core status among always-neutral genes, 24.4% positive-fitness core genes versus 19.9% accessory genes, and function-specific reversals. _{1 30 8 9 15}

The direct linkage adds a sampling tension: essential genes were 86.1% core versus 81.2% for non-essential genes, but clades with only 2 genomes can have trivially high core fractions and pangenome coverage varies. The main Escherichia coli clade was absent because it contained too many genomes; Keio, an E. coli BW25113 strain, mapped to the small s__Escherichia_coli_E clade at only 26.1% coverage. ₂₈

Persistent hypothetical genes were less core and had narrower ortholog breadth than annotation-lag genes, but had a higher essential fraction, 18.0% versus 13.4%. This does not overturn the aggregate conservation–essentiality gradient; it shows that essentiality can persist in lineage-restricted genes and that short-gene measurement and annotation biases must be separated from biological novelty. ₃₁

Module genes were 86.0% core versus 81.5% for all genes, but zero essential genes appeared in modules because essential genes are invisible to insertion-based ICA. This is a measurement tension rather than a biological contradiction. ₂₂

Pathway-level conservation conflicts with the gene-level gradient: latent complete pathways had mean conservation 0.869 versus 0.829 for active dependencies, with p=0.94. The newer analysis found Active Dependencies at mean core completeness 0.986 versus 0.975 for Latent Capabilities in only 7 matched model organisms, while its larger analysis links pathway variability—not pathway dependency—to openness. _{19 25}

Broad metal-important genes were 87.4% core versus 76.9% baseline, whereas DvH heavy-metal genes were 71.2% core and not significantly enriched after correction. Metal-specific genes were 84.8% core pooled, less enriched than general sick genes at 90.2%, while metal+stress genes were 94.3% core. These results should not be averaged because definitions, thresholds, organism sets, and stress histories differ. _{20 14 21}

The respiratory study adds a model–measurement tension: Complex I was required on quinate despite theoretical quinate production of 0.57 NADH per carbon, while it was dispensable on glucose at 1.50 NADH per carbon; FBA predicted zero NDH-2 flux because it optimizes ATP yield, and the cross-species NDH-2 comparison found no significant compensatory pattern (p=0.52). _{32 27}

Web of Microbes production matches cannot establish uptake, pathway completeness, or gene essentiality: the snapshot has no organism consumption actions, and 107 formula-only matches expanded to 900 candidate molecules. Likewise, only 859 of 17,222 ortholog families were universally essential, while 4,799 were variably essential and 11,564 never essential; 17 of 18 amino-acid pathways were complete in all 7 mapped organisms, yet DvH lacked predicted serine biosynthesis. _{26 27 13 19}

ADP1 showed moderate FBA–knockout concordance, whereas annotation-gap analysis found 42.5% baseline FBA accuracy across 574 organism–carbon-source combinations; FBA lacked mappings for 30/51 aromatic-network genes and predicted 0% Complex I essentiality despite 1.76× higher aromatic flux. Module-ICA had <1% strict KO precision despite strong cofitness, versus 95.8% for ortholog transfer; these are different prediction targets rather than a biological contradiction. _{3 4 27 16}

The counter-ion analysis found 39.8% overlap between metal-important and NaCl-stress genes, while removing shared-stress genes preserved core enrichment for 12 of 14 metals and the original 87.4% core, OR=2.08 conclusion. A separate analysis found 14.7% of metal-important genes sick under osmotic stress; the 2.7× discrepancy reflects stricter thresholds and different organism sets. _{11 21}

Interpretation and Limits

Direct TnSeq calls, measured condition differences, and quantified FBA–TnSeq comparisons provide the strongest evidence. No gene had measurements across all six database modalities, and incomplete overlap limits separation of biological essentiality from coverage, annotation, model, and measurement effects. Deletion measurements were single-timepoint growth ratios across only 8 carbon sources; 499 essential genes and 316 incomplete-data genes were excluded from the complete matrix. _{1 2}

Fitness measurements were biased toward rich media and standard stresses, ecological niches were underrepresented, single-gene knockouts did not capture epistasis, and the 43-bacterium cohort was primarily Pseudomonadota. Singleton or novel genes may lack fitness data because of poor transposon coverage rather than true neutrality. ₁₅ Only 48 organisms were fitness-tested among 293,000 genomes with pathway predictions; latent-capability rates were threshold-sensitive, and pathway-level conservation cannot detect partial erosion or gene redundancy. The newer study was limited to 7 matched organisms for Tier 1 comparison, 80 pathways, KEGG-based mappings, and genome-count sampling bias. _{19 25}

The conservation linkage extended taxonomic breadth but retained sampling limits: 44 of 48 organisms were linked, 34 had at least 90% coverage, and only 33 were used downstream; four organisms were unmatched because their GTDB species clades had too few genomes. Ten organisms were excluded from Phase 2 because they had <90% DIAMOND coverage, and Dyella79 was excluded because the Fitness Browser table used N515DRAFT_* locus tags while protein sequences used ABZR86_RS*, producing a 0% join rate. ₂₈

The essential-gene set is an upper bound because missing transposon insertions can reflect short genes, low-complexity regions, or scaffold-edge locations rather than essentiality; essential genes were slightly shorter on average, consistent with insertion bias. RB-TnSeq essentiality was typically measured in rich-media library-construction conditions, so genes essential only under stress would be missed. _{28 12} The reason is partly schema-level: across all 45 tables of the Fitness Browser dataset there is no explicit essentiality column, insertion count, or has-fitness-data flag, and gene.type encodes feature type (1=CDS, 5=pseudo, 7=rRNA, 2=tRNA) rather than essentiality, so putative essentials must be identified by comparing type=‘1’ genes against genefitness, which yields an upper bound on true essentiality. ₂₉ The converse error is equally operational: querying only genefitness silently omits approximately 14.3% of protein-coding genes — those with no transposon insertions and hence no fitness entries, a group reported as 82% core — so any fitness-versus-conservation analysis must add them back from the gene table or state explicitly that the most extreme fitness category is missing. ₂₉

Pangenome coverage limits conservation comparisons. Clades containing only 2 genomes can have trivially high core fractions, because presence in both equals 100% core. BBH orthology is conservative and can miss paralogs, gene fusions, and distant homologs, so some apparent orphan essentials may have undetected divergent orthologs. _{28 12}

The Web of Microbes source was a 2018 archived snapshot containing 37 organisms across 5 ENIGMA-funded projects, 589 metabolites, and 10,744 action observations; 332 metabolites, or 56.4%, were unidentified with an Unk_ prefix. It came from a single laboratory, included only 20 experimental organisms, may not represent the current GNPS2-hosted resource, and lacked consumption data. ₂₆

The truly-dark study found that 17,479 of 57,011 dark genes lacked pangenome links, leaving a 31% coverage gap; at the observed 16.3% truly-dark rate among linked genes, approximately 2,841 additional truly dark genes may exist among the unlinked set, including 2,208 with strong fitness phenotypes. Bakta may produce false-negative functional calls, BBH orthology covered only 32 of 48 Fitness Browser organisms, GC deviation is an imperfect HGT proxy, short genes are harder both to annotate and measure, and polar effects may contribute to strong phenotypes. ₃₁

The respiratory study is limited because NDH-2 has no growth data; its stoichiometric explanation uses theoretical rather than measured flux; only 4 organisms lacked NDH-2 in the cross-species comparison; annotation variation may miss orthologs; ACIAD3522 may have another function; text matching produced false-positive NDH-2 calls; and the planned pangenome KO co-occurrence analysis was not performed. ₃₂

Other limits include failed organism-to-clade linkage caused by boolean-string taxids, Fitness Browser taxonomic bias, missing or inconsistent environment metadata, AlphaEarth coverage of 83K/293K genomes (28%), and limited NMDC validation. The ENIGMA CORAL survey contained 6,705 genomes, 15,015 genes, 4,346 field samples across 596 Oak Ridge locations, and 213,044 ASVs, but no DvH gene-level fitness data. _{19 17 14}

The ENIGMA dataset contained 27,632 curves from 123 strains, 195 molecules, 303 plates, and 7.57M timepoints; 15,227 (55.1%) showed no growth, 9,861 (35.7%) were fit-ok, and approximately 9% failed fitting. Only 275 of 486 strain × condition Fitness Browser anchors showed measurable growth, a strain-name collision produced 1,751 spurious clinical genomes, and verified GTDB linkages fell from 32 to 20. ₁₈

The pitfalls audit adds operational limits: KBase Data Lakehouse is migrating from Delta to Iceberg; tenant and database names require explicit resolution; REST requests can return 504, 524, 503, or empty responses; direct Spark SQL is preferred; numeric fields and Spark DECIMAL values require explicit casting; billion-row pangenome tables require key filters and server-side aggregation; a 1 GB Spark Connect result cap makes large .toPandas() collections unsafe; and disabling automatic broadcast joins caused an NB10 job joining 13.7M rows with 18,989 species-taxonomy rows to hang for 17+ minutes. Short ENIGMA strain names are not globally unique, MetaPhlAn3 cohort names require synonymy reconciliation, and erroneous joins can create extreme apparent abundance differences or incorrect genus matches. ₂₉

Open Directions

• Perform condition-matched RB-TnSeq and complete-gene knockout experiments, then compare binary lethality, continuous fitness, and quantitative growth within the same media and genes to determine which discordances arise from assay modality. ₃
• Analyze insertion position, retained domains, and transcript structure for the 211 knockout-essential/TnSeq-dispensable genes to test whether partial gene products explain the discordance. ₃
• Reanalyze the 499 minimal-media-essential and 346 LB-essential genes with condition-matched FBA and TnSeq calls, focusing on the 177 genes whose flux classes changed; test the 625 condition-specific genes with replicated growth measurements and an expanded carbon-source panel. _{1 2}
• Reanalyze the eight-carbon-source data with explicitly condition-specific labels and dependence-aware aggregation to determine how inter-condition correlation changes essentiality estimates. ₃
• Fit combined FBA-plus-fitness-plus-proteomics predictors separately for lethal-versus-dispensable classification and quantitative growth ranking, then compare metrics appropriate to each target. ₃
• Add measured trace aromatic compounds to the FBA media definitions and rerun condition-specific flux simulations to test whether environmental assumptions explain the aromatic-degradation errors. ₃
• Compare universal, variable, orphan, and never-essential ortholog families across matched condition panels, and characterize the 7,084 orphan essentials for lineage-specific innovation versus missed orthology. _{27 12}
• Extend conservation analysis from binary essentiality to condition-specific fitness, especially genes with fitness < -2 under stress; use the 177,863 Fitness Browser–pangenome links to correlate mean fitness effects with core fraction and analyze the 3,683 essential-auxiliary genes for compensation of missing core functions. ₂₈
• Reanalyze fitness effects with condition-stratified distributions and interaction terms; test whether core genes’ heavier positive and negative tails persist within matched assays, whether epistasis explains additional bidirectional effects, and whether the 82% versus 66% contrast transfers outside rich media and standard stresses. ₁₅
• Combine transposon coverage metrics with singleton and auxiliary-gene estimates to distinguish true neutrality from callability-limited missing data, and compare conservation–fitness relationships across bacterial lineages beyond primarily Proteobacteria representation. ₁₅
• Harmonize the two essentiality–core analyses with identical organism filters, pangenome clade definitions, coverage thresholds, and essentiality denominators to determine why the core fractions are 82% versus 86.1%. _{15 28}
• Expand GapMind mapping from 7 to the remaining 38 organisms, including Escherichia coli K-12; compare pathway-level conservation with gene-level essentiality after direct per-step mapping and recompute latent rates across the 16 threshold combinations. _{13 19}
• Expand the four-way capability–dependency analysis beyond the 7 matched organisms, quantify KEGG-mapping failures using direct per-step gene evidence, and test whether all 66 Latent Capability reclassifications persist under preregistered thresholds. ₂₅
• Compare variable pathway count and ecotype number with pangenome openness using phylogenetic independent contrasts and genome-count-matched resampling across the 2,810 species; compare metabolic ecotypes with environment metadata using phylogenetic correction. _{25 19}
• Test Complex I, PQQ, iron-acquisition, and NADH-flux hypotheses with matched ADP1 perturbations; construct an ADP1 NDH-2 deletion mutant, measure NADH/NAD⁺ ratios on quinate, glucose, acetate, lactate, and urea, and test whether respiratory requirements follow flux rate rather than total NADH yield. _{5 32}
• Expand K03885 and K00330–K00343 searches across 27K species using KO-based orthology, test NDH-2/Complex I co-occurrence, characterize ACIAD3522 biochemically, and reanalyze quinate-versus-succinate proteomics. ₃₂
• Analyze the 3,683 essential-auxiliary genes with functional annotation and genomic-context methods to test whether they compensate for missing core functions or instead represent mobile-element-associated dependencies. ₂₈
• Combine BBH, profile-based homology, genomic-context analysis, and module-informed CRISPRi (CRISPR interference, transcriptional knockdown) tests to distinguish genuinely lineage-specific essential genes from divergent or missed homologs among the 7,084 orphan essentials. ₁₂
• Compare ICA modules with ortholog-transfer predictions while varying experiment count, component caps, membership thresholds, and Pfam/KEGG enrichment rules; test whether the 38 accessory module families are enriched for mobile elements, operon structure, or niche-associated environments. _{16 22}
• Test the 9,557 essential dark genes with CRISPRi, validate the 998 dark-gene operon pairs supported by conserved synteny and strong cofitness, and fold the 6,427 truly dark genes into structure-, Foldseek-, and condition-matched growth validation beginning with the 100 ranked candidates. _{17 31}
• Extend pangenome linkage to the 17,479 unlinked dark genes and test whether the estimated approximately 2,841 additional truly dark genes reproduce linked-set essentiality and phenotype distributions; test the Methanococcus-focused cofitness signal reported as r > 0.97 and characterize contiguous dark islands. ₃₁
• Refit genotype–condition models after ChEBI-based canonicalization and compare high-AUC substrates with condition-matched TnSeq and knockout essentiality. ₁₈
• Reanalyze the metal atlas with organism-specific tolerance normalization, phylogenetic independent contrasts, and removal of genes important for non-metal stresses; resolve locusId mismatches and repeat metal-specificity estimates across sick-rate thresholds from 1–20%. _{20 21}
• Obtain current Web of Microbes or GNPS2/Northen data with consumption actions; curate exact identities for the 107 formula-only ModelSEED matches, build a pathway-to-metabolite lookup table for GapMind, and test E/(E+I) against accessory genes and pangenome openness after strain-resolved matching. ₂₆
• Rebuild organism and strain joins with assembly accessions, genus checks, NCBI-taxid-backed synonymy, and explicit GTDB-version tracking; repeat environment-linked essentiality analyses after AlphaEarth NaN filtering and stratification by environmental versus human-associated samples, quantifying the effect of 52.7% unknown per-genome environment labels. ₂₉
• Execute large-scale joins with live catalog discovery, explicit casts, key-filtered Spark SQL, and checkpointed outputs, comparing results against legacy namespace and REST workflows. ₂₉
• Test whether any gene families are universally essential across all 43 organisms with fitness data, the cross-organism essential-family question the conservation–fitness synthesis explicitly leaves open. ₈
• Probe tol pal system function in Caulobacter using conditional-depletion strains, since Tol-Pal essentiality there complicates straightforward depletion experiments while remaining testable in such strains. ₇
• Cross-reference the universally essential families against the Database of Essential Genes to identify families essential in organisms outside the Fitness Browser panel, closing the gap between experimentally measured and computationally conserved essentiality sets. ₁₂
• Test whether complete pathways are actually essential for viability, identify pathways essential only under specific growth conditions, and determine whether pathway gaps cluster by phylogeny or ecology, using condition-varied growth data rather than rich-media RB-TnSeq alone. ₁₃
• Execute the CRISPRi designs attached to the top 50 essential dark-gene candidates, starting with Keio:14796, MR-1:200382, and Koxy:BWI76_RS08540, to convert operon-context hypotheses into measured phenotypes for the covering set’s 8,900 essential genes that transposon knockout cannot reach. ₁₇