Contributions
Metabolic Capability vs Dependency crossed GapMind pathway-completeness predictions with random barcode transposon sequencing (RB-TnSeq; genome-wide pooled mutant fitness profiling) data from the Fitness Browser and classified 161 organism-pathway pairs across 7 bacteria and 23 GapMind pathways 1. It found 57 (35.4%) Active Dependency pairs, 66 (41.0%) Latent Capability, 24 (14.9%) Incomplete but Important, and 14 (8.7%) Missing; the report suggests that the Incomplete but Important pairs point to annotation gaps or salvage routes 1. The project scored importance with a composite of 40% essentiality, 30% fitness breadth, and 30% fitness magnitude, and counted a pathway as complete when GapMind predicted it complete in >50% of conspecific genomes 1. It found mean core gene completeness of 0.986 for Active Dependencies versus 0.975 for Latent Capabilities, and attributed the small gap partly to its use of well-studied model organisms with near-complete, well-annotated core genomes 1. Only 7 of the 48 Fitness Browser organisms had matching GapMind genome data, which limited the analysis to 161 pairs, and the pilot step noted that Escherichia coli was absent from GapMind 1.
The project stratified fitness effects by condition type and found that all 66 Latent Capability pairs became fitness-important under at least one condition type, most often nitrogen limitation, stress, or carbon limitation 1. The report cautions that this result is partly a consequence of its median-based importance threshold, which splits pathways approximately 50/50 by construction, so that some pathways will inevitably cross it when many condition subsets are tested 1. It proposes calibrating the threshold against known essentials from essential_metabolome, and states that until then not every reclassification should be treated as biological confirmation 1. From this result the report draws a Black Queen Hypothesis implication that remains a hypothesis: pathway gene loss may be constrained by the breadth of conditions a species encounters 1.
The project extended the analysis to 2,810 GTDB (Genome Taxonomy Database) species that had at least 10 genomes in the pangenome database, covering 80 GapMind pathways with data extracted from the KBase Data Lakehouse spanning 293,000 genomes 1. It found that variable pathway count (pathways present in 10-90% of genomes) correlated with pangenome openness (fraction of accessory gene clusters) at raw Spearman rho=0.327, p=7.2e-71, and at partial rho=0.530, p=2.83e-203 after controlling for genome count 1. The report’s reading that pathway gain and loss is a major driver of pangenome dynamics remains an inference from correlation, not a tested mechanism 1. Within genera, the signal was positive in 13 of 18 genera tested but significant at p<0.05 in only 5 (Clostridium, Eubacterium, Mesorhizobium, Pseudomonas, Streptomyces), with the strongest correlations being rho=0.534 for Clostridium and rho=0.506 for Eubacterium, and smaller per-genus sample sizes limited statistical power for many groups 1. The project did not perform phylum-level stratification, which would require joining the GTDB taxonomy table, and used genus-level grouping from GTDB species clade names as a conservative phylogenetic control 1.
The project compared core-only with all-gene GapMind completeness and found the largest accessory-dependent gaps in amino acid biosynthesis: 0.146 for leucine and valine, 0.141 for arginine, and 0.140 for lysine and threonine 1. The report summarizes this as approximately 14% of species-level completeness for these pathways coming from accessory genes 1. It interprets this pattern as consistent with Black Queen dynamics and possible community-level sharing, while noting that the analysis does not demonstrate metabolite exchange 1. Among 225 species with at least 50 genomes and at least 3 variable pathways, Jaccard-based hierarchical clustering found a median of 4 metabolic ecotypes per species and a maximum of 8, seen in Alistipes onderdonkii and Barnesiella intestinihominis 1. Ecotype count correlated with pangenome openness at raw Spearman rho=0.262, p=6.8e-05, and at partial Spearman rho=0.322, p=8.0e-07 after controlling for genome count 1.