Polyhydroxybutyrate (PHB) biosynthesis is unevenly distributed across bacterial and archaeal diversity, enriched in environmentally variable habitats, associated with larger genomes, and often accessory in phylogenetically unusual lineages. These results support environmental selection and widespread horizontal gene transfer (HGT), while showing that the apparent relationship between PHB and environmental niche breadth is largely explained by genome size. 1
This synthesis is based on the Polyhydroxybutyrate Granule Formation Pathways — Distribution Across Clades and Environmental Selection analysis of 27,690 gtdb species, pangenome annotations, environmental metadata, AlphaEarth embeddings, and nmdc metagenomes. 1 The Prophage Gene Modules and Terminase-Defined Lineages Across Bacterial Phylogeny and Environmental Gradients provides a related cross-project comparison: across 27,702 bacterial species, genome size was the dominant predictor of prophage-module composition, but environment remained significant after genome-size stratification and phylogenetic controls. This supports and refines the PHB interpretation by showing that genome size can dominate broad ecological associations without eliminating pathway- or module-specific environmental structure. 2 The Plant Microbiome Ecotypes — Compartment-Specific Functional Guilds and Their Genetic Architecture further supports treating genome size as a confounder: marker count correlated with genome size at r = 0.44, p < 1e-300, and plant-associated dual-nature classifications varied from 54% in the lowest genome-size quartile to 87% in the top three quartiles. 3
Pathway distribution and completeness
phaC, the PHA synthase that catalyzes the committed step of PHB biosynthesis and serves as the pathway’s presence marker, occurred in 21.9% of 27,690 GTDB species, while 21.7% carried a complete PHB pathway defined as phaC plus phaA and phaB. 1 The analysis identified 118,513 PHB-related gene clusters across 19,496 species using eggNOG annotations. 1 The near-identical prevalence of phaC and complete pathways indicates that nearly all phaC-carrying species also possess the upstream biosynthetic enzymes. 1
PHB prevalence was highest in Pseudomonadota at 60.9% (4,544/7,456 species), followed by Myxococcota at 52.9% (83/157), Halobacteriota at 34.4% (89/259), Thermoproteota at 27.9% (65/233), Desulfobacterota at 21.9% (70/319), and Actinomycetota at 18.5% (587/3,172). 1 Campylobacterota (0/271), Gemmatimonadota (0/102), Nanoarchaeota (0/135), and Marinisomatota (0/90) had no detected PHB, while Patescibacteria had 0.4% prevalence (4/981). 1
The 6,067 phaC-carrying species spanned 95 of 142 GTDB phyla. 1 Pseudomonadota accounted for 74.9% of phaC-carrying species (4,544/6,067), and the five most PHB-enriched phyla accounted for 85.5% of phaC-positive species while comprising 30.4% of total species diversity. 1 The 85.5% figure does not agree with the report’s own per-phylum carrier counts for those five phyla, which do not appear to add up to that share of the 6,067 phaC-positive species. The report leaves this discrepancy unresolved, so the 85.5% value should be treated with caution. 1 At the order level, 23 orders with more than 20 species exceeded 50% phaC prevalence; Azospirillales and Rhodospirillales each reached 100%, Caulobacterales reached 99.3%, and Sphingomonadales reached 95.9%, whereas 117 orders were below 10%. 1
The gene-level census was keyed to KEGG orthologs (KOs). It found 86,318 phaA (K00626, beta-ketothiolase, pleiotropic) clusters across 17,969 species, 11,792 phaC (K03821, PHA synthase, committed step) clusters across 6,067 species, 9,617 phaB (K00023, acetoacetyl-CoA reductase, pleiotropic) clusters across 6,977 species, 6,130 phaP (K14205, phasin) clusters across 4,571 species, 4,656 phaZ (K05973, PHB depolymerase) clusters across 3,151 species, and 0 phaR (K18080, transcriptional regulator) clusters across 0 species. 1 Species-level categories were complete (6,005), synthase-only (62), precursors-only (12,869), accessory-only (555), and absent (8,199). 1 The 46.5% precursors-only category likely overestimates partial PHB capability because phaA and phaB also participate in general fatty-acid and thiolase metabolism. 1 The absence of phaR may reflect incomplete annotation or misassignment of K18080 rather than biological absence. 1
Environmental selection and the feast/famine hypothesis
The classical feast/famine hypothesis (Dawes & Senior 1973) proposes that PHB is most advantageous under temporal carbon fluctuation. The report frames its 27,690-species survey as the first pan-bacterial genomic test of this hypothesis and says its data strongly support it. That framing is the report’s own: the tests below are correlational and do not control for phylogeny. 1
PHB prevalence followed an environmental gradient: plant-associated species had 44.0% phaC prevalence (625 total), soil species had 43.6% (1,484 total), wastewater/engineered species had 34.5% (1,124 total), freshwater species had 25.5% (3,263 total), sediment species had 20.1% (1,020 total), marine species had 18.7% (3,010 total), human-associated species had 11.1% (1,237 total), human clinical species had 7.4% (2,472 total), and animal-associated species had 3.3% (3,711 total). 1
A chi-squared test of PHB presence against environmental variability category yielded chi2 = 1,656.36, p ~ 0, dof = 2. 1 This strongly supports the hypothesis that PHB genes are enriched in temporally variable feast/famine environments and depleted in stable or host-associated environments. 1 PHB also has functions beyond carbon storage, including stress resistance, redox balance, and cryoprotection, so environmental variability is supported as a selective force but is not necessarily the sole driver of PHB distribution. 1
The report sets these results beside cited prior work, which is not new measurement from this corpus. One cited estimate puts intracellular carbon storage (including PHB) in soil as a significant carbon pool, at 0.19-0.46x extractable microbial biomass. A cited survey of PHA depolymerase (degradation) genes across 3,078 metagenomes found uneven distribution, with the highest frequency in wastewater and the lowest in marine environments. That pattern broadly parallels the synthesis-gene gradient reported here, but it concerns degradation genes rather than phaC. PHB has been reported to confer multi-stress resistance (UV, osmotic, oxidative, temperature, freezing) beyond simple carbon storage. A cited methanotroph study found lineage-specific PHB capacity: all Type II methanotrophs were positive and all Type I were negative. The report reads its own family-level result (41 enriched, 62 depleted families within phyla) as showing that PHB is phylogenetically structured even at fine taxonomic scales. 1
The environmental enrichment supports, rather than replaces, the broader environment-embedding-geography interpretation of environmental specialization. 1 The prophage analysis supports this distinction: environment significantly affected prophage-module composition beyond host-family phylogeny and genome-size quartile, with PERMANOVA, a permutation-based multivariate analysis of variance, values of p=0.01 for each factor; genome size was dominant (F=212.99), followed by environment (F=30.04) and phylogeny (F=6.17). 2 Thus, environmental structure can persist after accounting for genome size, but its strength and biological interpretation remain feature-specific rather than automatically generalizable from prophage modules to PHB. 2
After stratification by genome size, high-variability environments still had higher phaC prevalence in every quartile: Q1 (0.4–1.8 Mbp) had 10.7% versus 2.5%, a 4.4x enrichment with p = 1.18 x 10^-11; Q2 (1.8–2.5 Mbp) had 18.5% versus 4.0%, a 4.6x enrichment with p = 3.25 x 10^-47; Q3 (2.5–3.7 Mbp) had 34.3% versus 10.9%, a 3.1x enrichment with p = 4.17 x 10^-70; and Q4 (3.7–14.0 Mbp) had 50.9% versus 37.4%, a 1.4x enrichment with p = 3.21 x 10^-14. 1 Even among the smallest genomes (Q1, 0.4–1.8 Mbp), species from high-variability environments were 4.4x more likely to carry phaC than those from low-variability environments. The report takes this as an environmental enrichment signal that is independent of genome size. The quartile tests do not control for shared ancestry, so the signal is shown to be independent of genome size but not of phylogeny. 1 The report’s interpretation summarizes this result as environmental enrichment within all four genome size quartiles (1.4–4.6x, all p < 1e-11). That summary conflicts with the quartile table’s Q1 p-value of 1.18 x 10^-11, which is not below 1e-11. Both values are kept here and are not reconciled. 1
Family-level tests found substantial heterogeneity in phaC prevalence. Among 248 families with >=10 species, Fisher’s exact tests (exact tests of association in 2×2 contingency tables), Bonferroni-corrected at alpha = 0.05, identified 41 enriched, 62 depleted, and 145 nonsignificant families. 1 Enriched families were concentrated in Pseudomonadota (Burkholderiaceae, Burkholderiaceae_B, Rhodocyclaceae, Caulobacteraceae, Sphingomonadaceae, Xanthobacteraceae, Rhodanobacteraceae, Legionellaceae). They also appeared in Actinomycetota (Mycobacteriaceae, Nocardioidaceae), Halobacteriota (Haloarculaceae, Natrialbaceae), Thermoproteota (Nitrosopumilaceae), Cyanobacteriota (Microcystaceae), and Bacillota (Bacillaceae_G). 1
The environments of these families support the species-level gradient: enriched families skewed toward freshwater (5/41) and wastewater (3/41), whereas depleted families skewed toward marine (18/62) and host-associated (12/62) environments. 1 This family-level support is weak, however, because the “other_unknown” category covered 30/41 enriched and 26/62 depleted families, reflecting sparse environmental metadata for many species. 1
Plant-microbiome results refine the interpretation of environmental structure by showing that compartment effects can be statistically real but weak: refined root, rhizosphere, and phyllosphere functional profiles had PERMANOVA R² = 0.071, while location-only db-RDA, distance-based redundancy analysis, R² = 0.060. 3 The same analysis attributed much of an earlier R² = 0.527 signal to taxonomic sampling and dominant clades, since excluding three genome-rich species per compartment reduced R² to 0.072. 3 This does not contradict the PHB quartile result, but it cautions that environmental labels and taxonomic composition can create large apparent effects unless genome size, phylogeny, and sampling structure are controlled. 3
Genome size and niche breadth
Among 2,008 species with sufficient AlphaEarth environmental-embedding representation (a minimum of 5 genomes per species with embeddings, with within-species embedding variance used as a proxy for niche breadth), PHB-positive species had median embedding variance of 0.3295 (n = 531), compared with 0.2472 for PHB-negative species (n = 1,477); the Mann–Whitney U statistic was 446,546 with p = 1.88 x 10^-6. 1 This difference held after restricting the comparison to phyla containing both PHB+ and PHB− representatives. 1 The raw PHB–niche-breadth association was rho = 0.106 with p = 1.77 x 10^-06. 1
PHB-positive species had median genome size of 4.34 Mbp versus 2.44 Mbp for PHB-negative species, with rank-biserial r = -0.592 and p ~ 0. 1 Elsewhere the report says PHB+ species have 1.9 Mbp larger genomes on average. The values it reports are medians, so that statement does not establish a mean difference. 1 Genome size correlated with niche breadth at rho = 0.302 with p = 1.5e-43. 1 After controlling for genome size using partial Spearman correlation, the PHB–niche-breadth association became partial rho = -0.047 with p = 0.037, representing a 56.3% reduction in effect size and a sign reversal. 1
Thus, the analysis refines the raw association: genome size, rather than PHB specifically, is the primary correlate of the apparent niche-breadth relationship. 1 The report concludes that the niche-breadth hypothesis is substantially confounded by genome size and cannot serve as independent evidence for the feast/famine hypothesis. Its reasoning is that larger genomes encode more metabolic pathways in general, and PHB is one of many pathways enriched in larger-genome species. 1 The environmental enrichment of PHB across genome-size quartiles nevertheless supports a genome-size-independent contribution from environmental variability. 1 The prophage result supports the genome-size caution: prophage cluster count correlated with genome size at rho=0.717, and genome size was the strongest factor in module composition, while within every genome-size quartile environmental effects on prophage module count remained significant with all p < 6.5e-78. 2 This cross-project agreement strengthens the case for treating genome size as a major confounder, not as evidence that all pathway-level environmental associations are spurious. 2
The report explains this confounding with prior comparative genomics. Genome size scales linearly with metabolic enzyme count (r > 0.9; Vieira-Silva & Rocha 2010) and predicts habitat breadth in soil bacteria (Barberán et al. 2014). Fluctuating environments are argued to select for larger genomes (Bentkowski et al. 2015). Stable oligotrophic environments are argued to select for streamlined genomes that shed “optional” pathways, including PHB, glycogen, and polyphosphate storage (Giovannoni et al. 2005, 2014). From this the report proposes a causal chain: environment variability → genome size → metabolic pathway count → PHB presence. In that chain, PHB presence is partly a downstream consequence of genome size rather than an independent adaptation. As support, the report cites the collapse of the PHB–niche-breadth association after controlling for genome size: partial rho = -0.047 versus raw rho = 0.106, where the raw Spearman correlation had p = 1.77 x 10^-06. Its summary also says PHB+ species appear to occupy broader environmental niches (raw p = 1.88e-06). That p-value comes from the Mann–Whitney comparison of embedding variance, not from the Spearman correlation, and the report says this association is largely explained by genome size. The chain is a proposed hypothesis, not a directly tested mediation model. The environmental enrichment within quartiles also shows that genome size does not account for all of the environmental signal. 1
The plant-microbiome study supports and extends this caution: in its marker analysis, approximately 33% of dual-nature species shifted to neutral after genome-size normalization, so broad marker-presence classifications can partly reflect genome capacity rather than ecological function. 3 Phylogenetic logistic regression or phylogenetic independent contrasts remain necessary because PHB presence is correlated with phylogeny and genome size. 1
The AlphaEarth result should be interpreted cautiously because 83K/293K genomes, or 28%, had embeddings, and the 2,008 analyzed species represented 7.2% of total species diversity. 1 The analysis may therefore be biased toward better-sampled lineages. 1 The Ecotype Correlation Analysis supports and refines this caution: across 172 species, the median partial correlation of environmental similarity with gene-content similarity was 0.0025 versus 0.0143 for phylogeny, with phylogeny dominant in 60.5% of species; environmental effects were significant and positive in 12 species (7.0%), significant and negative in 4 (2.3%), and non-significant in 156 (90.7%). 4 AlphaEarth embeddings covered only 28.4% of genomes in that analysis, and host-associated coordinates could represent collection sites rather than biological microenvironments. 4 These whole-gene-content results do not contradict the PHB-specific environmental enrichment: they refine it by suggesting that environmental selection may be concentrated in particular pathways rather than expressed as a strong genome-wide signal. 4
Accessory status and horizontal transfer
Phylogenetically discordant phaC distribution identified 311 potential HGT acquisition events in phaC-carrying species from families with less than 20% phaC prevalence and 278 potential HGT loss events in phaC-lacking species from families with more than 80% prevalence. 1 Among the 311 discordant phaC-positive species, 60.1% carried phaC as accessory genome, compared with 32.3% of all phaC-carrying species. 1 This elevated accessory fraction supports ongoing horizontal acquisition, but core/accessory status is a proxy rather than a directly reconstructed gene tree. 1 The analysis classified a gene cluster as core if it was present in >90% of genomes within a species and as accessory otherwise. The 60.1% versus 32.3% contrast is therefore indirect, pangenome-scale evidence for horizontal transfer of phaC, consistent with earlier single-species case studies cited by the report. 1 The cited prior work on PHA-pathway horizontal transfer includes Halomonas. The report presents its inferred transfer signals across diverse families (Lachnospiraceae, Chitinophagaceae, Pelagibacteraceae, Enterobacteriaceae) as extending that work. Its description of this as extensive HGT rests on inference from accessory status. The 311 potential acquisition events, with a 60.1% accessory rate (vs 32.3% overall), remain putative. 1
The principal putative recipient families were Lachnospiraceae with 38 acquisitions and 3.1% family phaC prevalence, Chitinophagaceae with 22 and 12.4%, Pelagibacteraceae with 13 and 5.3%, Enterobacteriaceae with 11 and 2.3%, and Planococcaceae with 10 and 12.8%. 1 The report also highlights SAR324, Bacillota_A, and Eremiobacterota as lineages with a high accessory phaC fraction and suggests they are active phaC recipients via HGT. This is a hypothesis inferred from accessory status, not a demonstrated transfer. 1 Overall, 5,371 species carried phaC as core and 1,959 as accessory, with some species carrying both core and accessory copies. 1
The HGT interpretation complements pangenome-integration and gene-function-acquisition-depth by connecting accessory pangenome status to a pathway-level ecological pattern. 1 It remains a hypothesis about transfer history until a directly reconstructed phaC gene tree identifies incongruent branches. 1
The T4SS–CAZy Environmental HGT supports and refines this indirect-HGT interpretation with an independent environmental gene-mobility analysis. Among 30,497 high-quality environmental MAGs, 6,652 (21.8%) carried type IV secretion system (T4SS) or conjugative machinery, and 92 carbohydrate-active enzyme (CAZy) families showed elevated co-occurrence with T4SS loci at distances of ≤10 kb; threshold validation remains pending. 5 GT2 glycosyltransferases showed 77 detected HGT events, including 32 normalized high-confidence cross-phylum events, with the strongest event spanning 8 phyla at a maximum divergence of 4.843. 5 These observations make chromosomal or integrative transfer a plausible route for environmental dissemination of accessory carbohydrate functions, but they do not directly establish transfer of phaC or validate the PHB-specific acquisition counts. 5
The same analysis found no CAZy genes on plasmids by ICEfinder; 12 integrative mobilizable elements occurred among the top 100 accumulators, and T4SS-positive genomes had 10× higher mobile genetic element density than other genomes (p<0.001). 5 This refines the PHB HGT hypothesis by identifying a possible integrative or chromosomal mechanism rather than simple plasmid mobilization, while the observational design and pending validation prevent mechanistic confirmation. 5
The prophage analysis supports the plausibility of a modular-transfer interpretation but does not directly validate PHB transfer: packaging, lysis, and lysogenic-regulation modules were consistently co-localized, whereas integration and anti-defense functions were more dispersed, and module-level environmental effects were decoupled from individual TerL-defined lineages. 2 This refines rather than strengthens the specific phaC HGT claim, because the prophage result demonstrates modular exchange as an ecological pattern while leaving the PHB transfer history unresolved. 2
Plant-microbiome evidence further refines the scale of the HGT inference. Species carrying plant-interaction marker singletons were 16 times more likely to carry transposase or integrase singletons (Fisher OR = 15.95, p = 8.8e-20), but marker singletons had an overall enrichment ratio of 0.78 relative to the genomic average, with Wilcoxon p < 1e-300. 3 At genus scale, plant-associated genera had a median 3.7 mobile elements per genome versus 2.8 in non-plant genera (Mann–Whitney p = 1.49×10⁻⁵). 3 The contrast between singleton co-occurrence, overall lower marker mobility, and genus-level mobilome enrichment supports the existing caution that accessory status is only an indirect transfer proxy and that transfer signals can reverse with aggregation scale. 3
NMDC cross-validation
A two-tier taxonomy mapping connected 3,014/3,492 NMDC taxon columns (86.3%) to GTDB genera with known PHB status: 2,336 columns through the gtdb_metadata NCBI-taxid-to-GTDB-genus bridge and 678 through direct genus-name matching. 1 PHB inference scores were calculated for 6,365 NMDC metagenomic samples as abundance-weighted sums of genus-level phaC prevalence, with a median 87.2% of taxonomic abundance matched to pangenome genera. 1 These scores are taxonomy-based inferences of community PHB capacity, not direct measurements of PHB genes or granules. 1
PHB inference scores correlated with depth at rho = -0.119 and p = 1.14 x 10^-21, temperature at rho = +0.088 and p = 1.86 x 10^-12, maximum depth at rho = +0.076 and p = 1.15 x 10^-9, minimum depth at rho = -0.055 and p = 1.05 x 10^-5, pH at rho = +0.049 and p = 7.95 x 10^-5, and ammonium nitrogen at rho = +0.044 and p = 4.49 x 10^-4. 1 The report’s summary rounds these to rho = -0.12, p = 1.1e-21 for sample depth and rho = 0.09, p = 1.9e-12 for temperature. 1 These effect sizes were modest, with |rho| < 0.12, and the measured abiotic values were point-in-time observations rather than direct measures of temporal variability. 1 The abiotic-correlation output contains 21 Spearman correlations of PHB score with abiotic variables, 6 of them significant at p<0.001. 1 The report reads the negative depth correlation as consistent with greater carbon-input variability at the surface and the positive temperature correlation as consistent with more dynamic carbon cycling, but both readings are inferential. 1
In genus-level cross-validation, 693 genera matched between pangenome and NMDC metagenomes. 1 Genera with at least 50% phaC prevalence had significantly higher NMDC abundance than genera with lower prevalence by Mann–Whitney testing (p = 8.41 x 10^-22). 1 The most abundant PHB-high genera included Mycobacterium, Pseudomonas, Cupriavidus, Burkholderia, and Methylobacterium. 1
The NMDC results support ecological relevance of the pangenome pattern, but they do not directly validate the feast/famine mechanism because NMDC studies were biased toward terrestrial and soil environments and their abiotic measurements did not quantify temporal variability. 1 The prophage study supports this cross-collection validation strategy: its taxonomy-based inference across 6,365 NMDC samples achieved 87.2% median matching coverage and found 57 significant module–abiotic correlations at FDR < 0.05, while agreeing with pangenome enrichment for head, tail, and anti-defense modules. 2 This refines the evidentiary status of the PHB NMDC analysis: abundance-weighted taxonomy bridges can reproduce collection-level functional patterns, but they remain indirect and do not establish pathway mechanism. 2
Tensions and analytical limits
The results support both a genome-size explanation for the apparent PHB–niche-breadth association and an environmental-selection explanation for PHB prevalence. 1 These claims address different contrasts: genome size largely explains the cross-species embedding association, whereas environmental enrichment persists within all four genome-size quartiles. 1
The ecotype analysis refines the scope of the environmental-selection claim rather than resolving it: environmental similarity was generally weak as a predictor of whole-genome gene-content similarity, while the PHB analysis found pathway prevalence enrichment across environmental-variability categories and genome-size quartiles. 4,1 The two analyses use different outcomes and environmental representations, so the genome-wide result cannot test whether PHB is an environmentally selected exception or whether the observed PHB pattern reflects residual phylogenetic and sampling structure. 4,1
The prophage evidence supports the first side of this tension by independently showing genome-size dominance alongside residual environmental structure, but it does not resolve whether PHB’s environmental enrichment is causal: prophage environmental effects were concentrated in structurally variable and anti-defense modules, whereas near-universal modules masked enrichment, and prophage annotation lacked dedicated detection. 2 Therefore, module-level environmental sensitivity cannot be assumed to imply the same mechanism for PHB. 2
The T4SS–CAZy analysis refines the transfer side of this tension but does not resolve it: it found syntenic and phylogenetic evidence for environmental transfer of GT2 carbohydrate functions, yet the relevant genes were not detected on plasmids and the threshold, BLAST, housekeeping-gene, and biome-factorization validations remain pending. 5 These results support a testable integrative-transfer hypothesis for environmental carbohydrate-gene dissemination, not a directly demonstrated mechanism for PHB or phaC. 5
Environmental interpretation is limited by sparse metadata: 34.9% of species had “other_unknown” as their primary environment. 1 The pangenome also included complete genomes and MAGs with variable genome quality and gene-detection rates. 1 The attempted PHA synthase class analysis classified all 11,792 phaC clusters as “other_pfam” because the eggNOG PFAMs column contained domain names such as Abhydrolase_1 and PhaC_N rather than the expected Pfam accession IDs PF00561 and PF07167. 1 Partial-correlation assumptions and imprecise or missing geographic coordinates additionally limited the ecotype analysis, particularly for host-associated organisms. 4
Plant-compartment labels were inferred primarily from NCBI isolation_source metadata; only 7,995 of 293,059 genomes, or 2.7%, had plant-associated annotations, and the plant study found only 11.7% Jaccard overlap between pangenome-derived plant genera and MGnify rhizosphere genera. 3 This supports the existing warning that collection metadata and metagenomic detection do not necessarily represent the same ecological phenomenon. 3
Open Directions
- Reconstruct a directly aligned phaC gene tree and test gene-tree/species-tree incongruence to distinguish horizontal acquisition from lineage sorting or annotation error. 1
- Apply phylogenetic logistic regression or phylogenetic independent contrasts to the 27,690-species pangenome table to test whether environmental variability predicts phaC presence after shared ancestry and genome size are controlled. 1
- Map PF00561 to Abhydrolase_1 and PF07167 to PhaC_N, then classify the 11,792 phaC clusters by PHA synthase class to test whether ecological distributions differ among synthase classes. 1
- Combine repeated environmental measurements with NMDC taxonomic abundances and PHB inference scores to test temporal variability directly rather than using point-in-time abiotic measurements. 1
- Reanalyze the AlphaEarth association after expanding embedding coverage beyond 83K/293K genomes and testing lineage-balanced sampling, asking whether the partial rho = -0.047 association is reproducible across underrepresented clades. 1
- Use genome-quality filtering and independent pathway evidence to test whether the 46.5% precursors-only category reflects genuine partial PHB capability or pleiotropic phaA/phaB metabolism. 1
- Test whether environmental associations are stronger for PHB and other specific functional modules than for whole-genome gene content by combining lineage-balanced pangenomes, COG functional categories, direct environmental metadata, and nonlinear or phylogenetically informed models. 4
- Apply the prophage study’s constrained, genome-size-stratified module-enrichment design to PHB and other accessory pathways, asking whether pathway-specific environmental associations persist after matched host phylogeny and genome-size controls. 2
- Use plant-compartment genomes to test whether PHB prevalence and phaC core/accessory status differ after matching genome size, phylogeny, and metadata quality, rather than treating plant association as a homogeneous environmental category. 3
- Compare direct mobile-element annotations with phaC gene-tree incongruence across matched genome- and genus-level analyses to determine whether transfer signals change with aggregation scale. 3
- Apply the pending T4SS–CAZy synteny permutation test, Node_4915 BLAST validation, housekeeping-gene null baseline, and biome-enrichment factorization to determine whether the observed integrative-transfer signal is robust; then test the resulting transfer signatures against phaC neighborhoods without treating GT2 as a direct proxy for PHB. 5
- Use the full 64-dimensional AlphaEarth embeddings, rather than variance alone, with dimensionality reduction to map PHB+ versus PHB- species in environmental space. The aim is to find the specific environmental axes that separate them and to test whether any environmental signal survives beyond the variance metric, which genome size confounds. The report proposes this; it is not a reported result. 1
- Use the NMDC
study_tableecosystem columns (ecosystem, ecosystem_category, ecosystem_type, ecosystem_subtype, specific_ecosystem) and theenv_triads_flattenedENVO terms to classify NMDC samples in more detail. This would test whether PHB inference scores track ecosystem type beyond point-in-time abiotic variables. The report proposes this; it has not been done. 1