Geographic and environmental structure in microbial genomic patterns

This concept examines whether alph-aearth satellite-derived embeddings, environmental metadata, metagenomes, and pangenome data capture ecological variation at a resolution relevant to microbial genomes and gene content. It also addresses how harmonization, coverage, coordinate quality, phylogeny, sampling composition, and analytical choices reshape those inferences. Evidence comes from AlphaEarth Embeddings, Geography & Environment Explorer, Ecotype Reanalysis — Environmental vs Human-Associated Species, Ecotype Correlation Analysis, Functional Dark Matter — Experimentally Prioritized Novel Genetic Systems, Gene Function Ecological Agora, Genotype × Condition to Phenotype Prediction from ENIGMA Growth Curves, Harvard Forest Long-Term Warming — DNA vs RNA Functional Response, Lab Fitness Predicts Field Ecology at Oak Ridge, Lanthanide Methylotrophy Atlas — Summary, Global Biogeography of Environmental Bacterial Metal Resistance, Metal Resistance Ecology — Phylogenetic Conservation vs. Environmental Selection, Polyhydroxybutyrate Granule Formation Pathways — Distribution Across Clades and Environmental Selection, Plant Microbiome Ecotypes — Compartment-Specific Functional Guilds and Their Genetic Architecture, Prophage Gene Modules and Terminase-Defined Lineages Across Bacterial Phylogeny and Environmental Gradients, Carbon Source Utilization Predicts Ecology and Lifestyle in Pseudomonas, Soil Microbial Dark Matter and the Clay Shield Null Result, Soil Metal Concentrations Drive Functional Gene Shifts, pitfalls, SNIPE Defense System in the KBase Data Lakehouse Pangenome, Community Metabolic Ecology via NMDC × Pangenome Integration, and Metabolic Capability vs Dependency. _{1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21}

Scope and ecological validity

A central measurement problem is collection-site versus organismal microenvironment mismatch: geographic coordinates may describe where a sample was collected rather than the physical, chemical, or host-associated environment experienced by the organism. For host-associated bacteria, the relevant environment may be a host, tissue, body site, or other organism-proximal microenvironment. ₃ Thus, an available coordinate is not necessarily an ecologically valid predictor of gene content.

The AlphaEarth explorer supports this concern by showing that embeddings encode spatially autocorrelated context while metadata are unevenly representative of ecological conditions. AlphaEarth embeddings covered 83,287 genomes, 28.4% of 293,059 KBase pangenome genomes; 79,449 had valid (non-NaN) embeddings across all 64 dimensions and 3,838 had NaN in at least one dimension. ₁ The pitfalls audit supports treating alph-aearth embeddings as a partial sampling frame. In 3,838 of 83,287 genomes, or 4.6%, at least one of the 64 dimensions was NaN, which makes UMAP, cosine distance and other operations fail or silently return NaN. Embeddings also required valid latitude and longitude, which are often missing from NCBI metadata, especially for clinical isolates. ₁₈ The coverage section of the same digest states that only 83,227 of 293,059 genomes (28.4%) have environmental embeddings. That count differs from the 83,287 used in the digest’s NaN section and by the explorer, and neither source explains the discrepancy. _18,1 The cause of these NaNs is unknown. It may be missing satellite imagery at those coordinates, or coordinates that fall outside satellite coverage such as polar regions or ocean. Affected records must be filtered before any embedding-based analysis. _1,22 The pitfalls audit specifies dropping any genome with a missing embedding dimension before any embedding-based computation, which reduces the dataset from 83,287 to 79,449 genomes. ₁₈

Limited coverage is a recurring caveat wherever AlphaEarth is used. The explorer reports coverage of only 28.4% of all genomes, biased toward genomes with valid lat/lon metadata. The ecotype analysis likewise states that 28.4% genome coverage limited the environmental signal available for analysis. _1,3 The pangenome AMR atlas reports that only 28% of genomes have embeddings, biased toward genomes with geographic metadata. Functional-dark-matter tests describe embedding coverage of only 28% of genomes (83K/293K) and inconsistent NCBI isolation-source metadata, and both limit the power of biogeographic tests. _23,4 The pitfalls digest likewise gives ~28% (~83K/293K) and attributes the gap to sparse NCBI lat/lon. It advises a LEFT JOIN with an explicit coverage check, so that unmatched genomes are retained rather than silently dropped. ₁₈ The same ~28% (83K/293K) coverage was the reason a planned correlation of metabolic ecotypes with AlphaEarth environmental niche breadth was not run: coverage would significantly limit the sample size. This is a test that was never run, not a negative result. ₂₁ The ecological-agora audit covered an 18,989-species substrate. It found per-genome BioSample mapping at 100%, but the environmental fields were uneven: 76.8% isolation_source, 88.8% geo_loc_name (country/region text), 48.6% env_broad_scale and 47.3% lat_lon (numeric geographic coordinates). AlphaEarth 64-dim satellite-imagery env vectors covered 27.2% (5,157 species), and elsewhere the same report gives AlphaEarth as 27% species coverage. The reports do not reconcile their different coverage percentages, and the agora percentages describe its species substrate rather than all genomes. ₅ This supports cultivation-collection-bias-in-ecological-genomics: the embedded subset is a partial, metadata-selected sampling frame. _1,23

The embedded subset was strongly human-associated: Human clinical numbered 16,390 (20%), Human gut 13,466 (16%), and Human other 1,669 (2%), together 38%; Soil numbered 6,073 (7%), Marine 5,850 (7%), and Freshwater 5,840 (7%). Other contained 13,944 genomes (17%), including site-specific labels, generic water and bodily-fluid labels, and uncaptured clinical sites. ₁ The pitfalls audit found 52.7% unknown per-genome NCBI environment labels, or 94,957 of 180,025 genomes, although species-level majority-vote labels reached 91% coverage. ₁₈ The discoveries digest attributes the 38% human-associated share (clinical 20%, gut 16%, other 2%) to NCBI’s bias toward pathogen sequencing. Clinical isolates have good geographic metadata from epidemiological tracking, so they disproportionately have embeddings. ₂₂

This composition supports sampling-depth-and-downsampling-effects and refines ecological interpretation: environmental samples—Soil, Marine, Freshwater, Extreme, and Plant—showed mean cosine distance increasing from 0.27 below 100 km to 0.90 above 10,000 km, a 3.4× gradient, whereas human-associated samples increased from 0.37 to 0.75, a 2.0× gradient. ₁ A pooled analysis could therefore weaken an environmental signal when a large human-associated subset occupies relatively interchangeable urban environments, but this remains a confounding hypothesis rather than evidence that environmental similarity generally predicts gene content. ₁

The explorer interprets the flatter human-associated gradient as an effect of satellite imagery: hospitals and clinics worldwide share a similar urban built environment. It cautions that the pooled “All samples” curve (2.0x ratio) blends the two signals and is dominated by the 38% human-associated fraction. On that reading, the embeddings are most informative for environmental microbiology samples and least informative for clinical isolates. ₁ The explorer proposed that this clinical bias may partially explain the weak ecotype signal (median partial correlation 0.0025), and that restricting the ecotype analysis to environmental samples could reveal a stronger signal. The discoveries digest recommended that any analysis using AlphaEarth embeddings as environment proxies should stratify by sample type, or at minimum exclude human-associated samples. _1,22 The pitfalls audit frames the rule slightly differently. It says to always stratify AlphaEarth analyses by environment category, at minimum computing results separately for environmental and human-associated samples. When embeddings serve as environment proxies, as in ecotype analysis, it suggests considering excluding human-associated samples entirely. ₁₈ Both are proposals rather than demonstrations, and the reanalysis described below contradicts the dilution explanation. ₂

The original ecotype analysis extracted 13,381 genomes across 224 species and produced correlations for 172 species. ₃ A partial correlation, defined here as association between two variables after accounting for a third, had a median value of 0.0025 for environment and 0.0143 for phylogeny. Phylogeny dominated in 60.5% of species versus 39.5% for environment; 12 species (7.0%) showed a significant positive environmental effect, 4 (2.3%) a significant negative effect, and 156 (90.7%) no significant environmental effect. ₃ These results support the interpretation that vertical inheritance generally dominates genome-wide gene-content similarity while environmental adaptation may be concentrated in selected gene subsets. ₃ The discoveries digest reports the same 172-species medians (0.0025 for environment, 0.0143 for phylogeny) and the same dominance shares (60.5% versus 39.5%). It also reports a null comparison: there was no significant difference between “environmental” bacteria, where lat/lon is meaningful, and host-associated bacteria (p=0.66). ₂₂ The reanalysis below supports this earlier null, because it also found no environmental advantage. ₂

The reanalysis contradicts the prediction that distinguishing or removing human-associated composition alone would reveal a stronger environmental gene-content signal. Among 224 species with at least 20 embedded genomes and at least 30% coverage, 106 (47%) were majority human-associated, 47 (21%) majority environmental, and 71 (32%) Mixed/Other. Across 183 species, environmental species had median partial correlation 0.051 versus 0.084 for human-associated species; Environmental > Human-associated gave U=1536 and p=0.83. Continuous fraction analyses found no relationship for environmental genomes (rho=-0.085, p=0.25) or human-associated genomes (rho=0.030, p=0.69). Environmental species had a higher NaN partial-correlation rate, 10/47 = 21% versus 7/100 = 7%; this filtering would, if anything, favor detection of a stronger environmental signal, which was not observed. ₂

The reanalysis used all embedded genomes, including up to 3,505 per species, whereas the original used diversity-maximizing downsampling with a maximum of 250 and different genome sets. Its overall median partial correlation was 0.081 across 183 species versus 0.003 in the original analysis; these absolute values are not comparable across methodologies. ₂ The weak signal therefore should not be interpreted as proof that host-associated organisms lack environmental structuring: collection-site variables may remain biologically misaligned with host microenvironments. ₃ The ecotype report also lists that AlphaEarth embeddings may not fully capture ecologically relevant environmental variation. ₃

Harmonization, coverage, and sampling frame

Coverage varied sharply by metadata field among the 83,287 AlphaEarth genomes: cleaned latitude and longitude were available for 83,286 genomes (100.0%), geo_loc_name for 83,270 (100.0%), isolation_source for 76,295 (91.6%), host for 53,091 (63.7%), env_broad_scale for 34,800 (41.8%), env_local_scale for 31,541 (37.9%), and env_medium for 31,483 (37.8%). ₁ The NCBI ncbi_env table contained 334 distinct harmonized attribute names across 4.1M rows; collection_date covered 273K genomes, geo_loc_name 272K, and isolation_source 245K. ₁

A keyword workflow mapped 5,774 unique isolation_source values into 12 broad categories and captured 71% of genomes with a label; 17% were Other and 12.5% Unknown. ₁ Structured env_broad_scale was cleaner but covered only 42% and included bare ENVO identifiers such as ENVO:00000428 and generic values such as not applicable and missing. ₁ Keyword matching offers broader coverage but can miss site-specific or non-English labels, whereas structured metadata offers greater semantic consistency at the cost of coverage. This supports ontology-and-category-schema-sensitivity and environmental-resistome, and refines interpretation of category-dependent ecological distributions. ₁ Unresolved habitat labels also limit community-level analyses. In an nmdc community metabolic study, 61 “Unknown” ecosystem samples have no resolved habitat identity. The report proposes, without establishing it, that they are a mix of soil subtypes or sediment environments. It names NMDC ENVO annotations or study-level metadata as ways to resolve them and improve power for its H2 test. ₂₀

The explorer supports treating cluster–environment correspondence as partly dependent on labeling decisions. The environmental-only reanalysis remains testable, but category harmonization and coordinate-quality decisions remain potential sources of ecological misclassification. ₁ The environmental tests also connect to pangenome-integration because they used pangenome gene-content similarity rather than a single phenotype. ₃

Geographic structure and coordinate quality

Across 50,000 genome pairs with good coordinates, mean embedding cosine distance increased from 0.41 for <100 km across 231 pairs, to 0.51 for 100–500 km across 1,058, 0.56 for 500–1K km across 2,016, 0.66 for 1K–2K km across 3,779, 0.78 for 2K–5K km across 5,824, 0.80 for 5K–10K km across 20,935, and 0.82 for 10K–20K km across 16,107. ₁ The relationship was strongest below 2,000 km and plateaued above 5,000 km, supporting environmental distance-decay rather than a simple raw-distance effect. ₁ The report interprets the embeddings as encoding geographic/environmental context derived from satellite imagery at each genome’s sampling location. It names likely features as climate, land use, vegetation type and urban/rural character, but these are inferred rather than measured. ₁

The explorer places this pattern against published work, which is literature context rather than project measurement. Pearman et al. (2024) found that macroalgal microbiome composition responds primarily to environmental variation rather than geographic separation. Liu et al. (2022) found that metal contamination explained more variation in mining-soil fungal community structure (4.16%) than geographic distance (1.21%). Zhang et al. (2022) found that rare and abundant bacterial subcommunities show similar distance-decay patterns but differ mechanistically. Wang et al. (2024) found that aridity shapes distinct soil-bacterial biogeographic patterns, with higher similarity-distance decay rates across environmental gradients. ₁ The report reads these studies as supporting its view that the embeddings capture environmental features more than raw distance. The support is indirect, because the cited studies concern community composition rather than satellite embeddings. ₁

Distance-decay also appears at metre scale in the subsurface. The 9 sso wells form a 3×3 grid spanning ~6 m. Across them, sediment-associated communities showed significant distance-decay of similarity (Mantel Spearman ρ = 0.323, p = 0.029, 9,999 permutations). Mean Bray-Curtis dissimilarity was 0.747, ranging from 0.558 (U3–M6) to 0.872 (U3–L9). NMDS stress was 0.067, and Procrustes correspondence with the physical grid was only marginal (m² = 0.379, p = 0.080). ₂₄ This refines the embedding result: it shows measurable community distance-decay within a grid far smaller than the <100 km nearest embedding bin. It links to subsurface-hydrogeological-zonation. _24,1

Coordinates are metadata, not ground truth. Of 83,286 genomes with coordinates, 50,109 (60.2%) were Good, 30,469 (36.6%) Suspicious cluster, and 2,708 (3.3%) Low precision. Suspicious clusters included potentially institutional addresses and legitimate field sites: Rifle, Colorado (39.54, -107.78; 1,883 genomes), Saanich Inlet (48.36, -123.30; 1,529), Siberian soda lakes (52.11, 79.17; 812), Pittsburgh (40.44, -79.97; 630), and Lima (-12.0, -77.0; 1,641). ₁ This supports treating coordinate quality as part of ecological validity: removing clustered coordinates could remove well-sampled environmental campaigns, while retaining them could preserve clinical or institutional concentration. ₁ The Suspicious cluster flag is a heuristic for institutional addresses: shared coordinates with >50 genomes and >10 species. Legitimate field sites it flagged include Rifle (DOE IFRC), Saanich Inlet (an oceanographic time series) and Siberian soda lakes (extremophile campaigns). The discoveries digest recommends checking isolation_source homogeneity at each location to distinguish field sites from institutional addresses. _1,22 The ecotype analysis adds that geographic coordinates from NCBI metadata are often missing or imprecise, which reduces the quality of environmental associations. ₃ The pitfalls audit lists the failure modes: ncbi_lat_lon in gtdb_metadata is most often NULL. Otherwise it is often a malformed string such as “missing” or “not collected”, or a low-precision value such as “37.0 N 122.0 W”. ₁₈

The reanalysis refines the geographic result: environmental embeddings had a 3.4× geographic-signal ratio versus 2.0× for human-associated embeddings, yet the stronger environmental geographic signal did not produce stronger species-level gene-content correlations. The reanalysis offers two untested explanations. First, embedding similarity is not ecological relevance: the climate, vegetation and land-use variation that environmental embeddings capture may not strongly predict which genes a bacterium has. Second, human-associated species such as klebsiella or Enterococcus have global epidemiological patterns in which different lineages dominate different regions. Embedding–gene-content associations in these species could therefore be real but epidemiological rather than ecological. Both remain hypotheses. ₂ Within one constant-protocol soil study, local vegetation and geography remained associated with eukaryotic signal, whereas cross-study environment prediction had out-of-study R²=-0.30 and within-study five-fold R²=+0.17 ± 0.06. ₂₅ The pitfalls audit found a pooled near-versus-far ratio of 2.0×, versus 3.4× for environmental samples and 2.0× for human-associated samples. It warns that the pooled curve blends two populations and therefore underestimates the true geographic signal in the embeddings. ₁₈

The within-study contrast came from neon soil. Eukaryotic fraction differed across 11 local-vegetation levels (Kruskal H=119.1, p=7.6×10⁻²¹). Medians were 23% in sedge/forb herbaceous soil, 14% in emergent wetland, 13% in dwarf scrub, 2% in evergreen forest, and ≈0 in deciduous forest, cropland and pasture. Eukaryotic fraction also differed across 47 sites (Kruskal H=310.4, p=2.4×10⁻⁴⁶). It was highest at Arctic tundra sites (Utqiaġvik 30%, Caribou-Poker Creeks 22%, Toolik 17%) versus temperate forests (2–8%). ₂₅ As literature context, the report cites Holm et al. (2025, PMID 40708004), who link vegetation, land use and geography to soil microbial variation. ₂₅ This supports environmental and geographic structure within a constant protocol. The negative cross-study R² refines that structure as study-confounded; see study-batch-confounding-of-environmental-associations. ₂₅

The global metal-resistance analysis supports spatial structure but refines it as sampling-sensitive: among 289 5° grid cells with at least 5 MAGs, 11 hotspots had OR>2 and q<0.05 and 3 coldspots had OR<0.5 and q<0.05. The strongest hotspot was the Atacama/Andean region of Chile at lat=-25°, lon=-70°, with 21.8% prevalence, OR=9.83, q=7.6e-12, and n=101 MAGs. ₁₀ Of 260,652 environmental MAGs, 30,497 remained after host-associated-biome filtering and 22,356 had usable coordinates; 2.8% carried at least one metal-resistance type. Soil had 5.8% prevalence, OR=5.05, q=1.8e-82; marine had 1.2%, OR=0.20, q=9.2e-80; rhizosphere had 1.9%, OR=0.66, q=0.30 (NS). ₁₀

Embedding clusters, ecological compartments, and gene content

UMAP, a nonlinear dimensional-reduction method, projected the 64-dimensional embeddings into two dimensions, and DBSCAN, a density-based clustering method, identified 320 clusters. Many were dominated by one harmonized environment category; Air, Extreme, and Plant concentrated in a few clusters, while Human gut and Human clinical genomes were distributed across many clusters, potentially reflecting geographic substructure. Phylum structure was partly confounded with environment; Campylobacterota were predominantly gut-associated. ₁ The ecological-agora analysis refines this result by grouping the 4,633 P1B species that had AlphaEarth 64-dim embeddings into k=10 environmental clusters with kmeans-clustering (k-means). The report does not reconcile this count with the 5,157 species it lists as covered. Of 90 cyanobacteriia species, 35.6% fell in cluster 0 (marine+sponge) and 17.8% in clusters 5+8, a marine-leaning profile. Of 50 mycobacteriaceae species, 34.0% fell in cluster 1 (gut+sludge) and 20.0% in cluster 3 (gut), a host-clustered profile. The 729 bacteroidota species spread across clusters 0/1/3 at ~18% each, consistent with a wide gut and marine niche range. ₅ Its figure figures/p4d1_alphaearth_env_cluster_panel.png shows the k=10 panel, with Cyanobacteriia C0 marine at 35.6%. All three pre-registered atlas findings grounded cleanly in expected biomes at p < 10⁻¹¹: mycolic-acid Innovator-Isolated, PSII Innovator-Exchange and Bacteroidota PUL. BacDive phenotype anchors also confirmed two of the three, and AlphaEarth env-cluster concentration confirmed all three. This is concentration consistent with expectation, not a causal test. ₅ The biome assignments behind these anchors used mgnify as the primary source, ENVO as the secondary and isolation_source as the tertiary for the 18,989 P1B species representatives. Clade × expected-biome enrichment was tested with one-sided Fisher’s exact tests against the atlas distribution. ₅ The same clusters differed in per-species mean recent-gain counts across regulatory, metabolic and mixed KO categories: 3,178, 10,891 and 6,418 in cluster 5 (mouse-gut + marine + acid mine drainage); 3,091, 10,063 and 6,008 in cluster 8 (marine + human-gut + soil); 2,652, 8,081 and 5,780 in cluster 1 (gut + sludge); 2,627, 8,173 and 5,446 in cluster 3 (gut); 1,303, 5,314 and 3,087 in cluster 9 (soil + grassland); and 490, 1,475 and 1,100 in cluster 6 (hypersaline + Microcystis). The report reads mixed-biome clusters 5 and 8 as carrying the highest recent-gain density and hypersaline/extremophile cluster 6 the lowest. This links environmental clustering to gene-function-acquisition-depth, but it is a descriptive cluster-level contrast, not a causal test. ₅

UMAP and DBSCAN are exploratory rather than direct demonstrations of ecological units. UMAP depends on n_neighbors and min_dist; DBSCAN produced 320 clusters with eps=0.5, which may be too fine-grained, and coarser clustering could better match environment categories. ₁ This refines ecotype-clustering-validity: clusters require validation in the original 64-dimensional space, across parameter settings, against independently harmonized environmental variables, and through reproducible genomic or functional differences. ₁

ENIGMA identified 8 metabolic guilds across 20 taxonomic orders from 7,167 KO presence/absence features across 123 strains; guilds aligned with, but were not identical to, taxonomy. ₆ Across 46,389 genome × condition pairs from 727 genomes and 363 conditions, binary growth AUC was 0.620 under genus-blocked validation and KO × condition interactions increased mean AUC to 0.653 in 80 of 106 held-out genera. Amino acids and nucleosides had AUC 0.775 and 0.780, metals 0.605, antibiotics 0.619, and nitrogen 0.435; continuous µmax, lag, and yield-related max_A had negative R² under genus-blocked holdout. ₆ ENIGMA also places local co-occurrence on a global environmental axis, using 464K global 16S samples. The acidic rhodanobacter–ralstonia–Dyella co-occurrence cluster (Cluster B) had mean pH 5.43, versus 6.78 for the neutral-pH Brevundimonas–Caulobacter–sphingomonas cluster (Cluster A), a delta of -1.35. Mean temperature was 22.6 C versus 15.7 C, reported as +6.9 C. The report’s summary gives the pH means as 5.4 and 6.8 while keeping the 1.35-unit difference. The acidic cluster occurred in 3,720 exclusive global samples versus 179,023 for the neutral-pH cluster, reported as 48x rarer. Its dominant environments were aquatic 39%, soil 23% and plant 18%, versus aquatic 33%, plant 23% and soil 22% for the neutral-pH cluster. ₆ The report interprets this as a global pH-driven niche partition that mirrors the oak-ridge-field-research-center contamination gradient, where nitric acid leachate lowers pH in plume wells. It calls this a fundamental microbial niche axis rather than site-specific adaptation. That reading is observational and extrapolated, so it suggests the hypothesis of a pH niche partition rather than demonstrating one. The report also expects Cluster B strains to show different pH-dependent growth profiles and cluster-distinguishing genomic features to predict pH tolerance; both expectations are proposed, not validated. ₆

Pseudomonas supports environmental effects on pathway content but refines them toward deep-clade structure. Among 12,732 genomes across 433 species clades, Pseudomonas s.s. contained 19 species and 6,905 genomes, Pseudomonas_E 398 species and 5,687 genomes, and P. aeruginosa 6,760 genomes (53%). Isolation sources were classifiable for 8,171/12,732 (64.2%); among 54 free-living and plant-associated species, carbon profiles associated with environment under 999 permutations (p=0.006), but Random Forest balanced accuracy was 0.408 +/- 0.169 against a 0.250 chance baseline. ₁₅ Among 7 Pseudomonas s.s. and 189 Pseudomonas_E species with at least 5 genomes, 43 of 62 pathways differed after Mann–Whitney testing with Benjamini–Hochberg false-discovery-rate (FDR) correction (q<0.05). The first five PCs captured 74.9% of variance; PC1 explained 31.2% and PC2 17.9%; between-environment distance was 2.054 versus 1.890 within, p=0.006. ₁₅

The plant analysis supports compartment structure but contradicts strong species-level separation. In a refined 17-marker panel across 607 species, PERMANOVA gave pseudo-F=23.2, p=0.001, R²=0.071; PERMDISP gave F=15.6, p=0.001; db-RDA, distance-based redundancy analysis, gave location-only R²=0.060, p=0.001. Earlier R²=0.527 fell to R²=0.072 after excluding three genome-rich species per compartment, an 86% loss. ₁₃ Root-associated species were dominated by Rhizobium, Mesorhizobium, and Bradyrhizobium; phyllosphere species by Sphingomonas and Methylobacterium; Pseudomonas_E was the only genus significantly represented across all compartments. ACC deaminase, T3SS, nitrogen fixation, and quorum sensing had OR=69.3, 65.6, 14.5, and 24.1, respectively, but marker presence does not establish activity or phenotype. ₁₃

The PHB analysis supports environmental context while refining it with genome-size control. Across 27,690 GTDB species, 21.9% carried phaC and 21.7% had a complete PHB pathway; prevalence was 44.0% plant-associated, 43.6% soil, 34.5% wastewater/engineered, 25.5% freshwater, 20.1% sediment, 18.7% marine, 11.1% human-associated, 7.4% human clinical, and 3.3% animal-associated. High-variability environments exceeded low-variability environments by 4.4×, 4.6×, 3.1×, and 1.4× across genome-size quartiles Q1–Q4, with p=1.18×10^-11, 3.25×10^-47, 4.17×10^-70, and 3.21×10^-14. Yet a niche-breadth test used AlphaEarth 64-dimensional embedding variance as a proxy, in 2,008 species with a minimum of 5 genomes per species with embeddings. PHB-positive species had median embedding variance 0.3295 (n = 531) versus 0.2472 for PHB-negative species (n = 1,477) (Mann–Whitney U=446,546, p=1.88×10^-6). The difference held after restricting to phyla containing both PHB+ and PHB− representatives. Genome size correlated with niche breadth at rho=0.302, p=1.5e-43, but partial rho after genome-size control was -0.047, p=0.037, a 56.3% reduction and sign reversal. ₁₂ The 2,008 species represent 7.2% of species diversity. With only 28% of genomes (83K/293K) embedded, the niche-breadth analysis may be skewed toward better-sampled lineages. The report proposes, but did not perform, mapping PHB+ and PHB− species across all 64 embedding dimensions with dimensionality reduction to identify the specific environmental axes that separate them. ₁₂

Prophage modules show the same layered structure. In 1,773 species, PERMANOVA on Bray-Curtis module composition found significant effects of genome-size quartile, environment, and family phylogeny, each p=0.01; genome size dominated (F=212.99), followed by environment (F=30.04) and phylogeny (F=6.17), and genome size correlated with cluster count at rho=0.717. Environmental effects remained within every genome-size quartile, all Kruskal–Wallis p<6.5e-78. ₁₄ AlphaEarth embeddings add a continuous axis: among 2,008 species with ≥5 embedded genomes, environmental niche breadth correlated with prophage module count after controlling for genome size (partial Spearman rho=0.468, p=8.41e-110). ₁₄ The report cautions that only 28% of genomes have embeddings, biasing the subsample toward clinically and environmentally well-sampled lineages. ₁₄ Across 27,702 bacterial species, 4,005,537 eggNOG-derived prophage-associated clusters were found; packaging and lysogenic regulation occurred in 100.0% of species, lysis 99.9%, integration 99.1%, anti-defense 64.3%, head morphogenesis 56.1%, and tail 55.6%; 34.9% carried all 7 modules. ₁₄

The prophage study refines complete-lineage interpretations toward modules: MMseqs2 at 70% AAI grouped 38,085 TerL sequences from 11,789 species into 10,991 lineages; 6,921 (63%) were singletons. Thresholds of 50% and 80% AAI gave 4,001 and 16,283 lineages. Of 824 lineages with at least 5 species, 325 were specialists and 499 generalists, but 0/500 tests showed significant environment-specific enrichment after FDR correction. ₁₄

The SNIPE result supports testing accessory defense systems but refines marker choice. DUF4041/PF13250 occurred in 4,572 clusters across 1,696 species and 33 bacterial and archaeal phyla; 13.3% were core, 30.7% accessory, and 56.1% singleton, so 86.7% were accessory or singleton. Among species with AlphaEarth embeddings, SNIPE-bearing n=1,069 and non-SNIPE n=13,977; Welch tests found 22 of 64 AlphaEarth dimensions significantly different at Bonferroni-corrected p < 0.05, largest at A19 (Cohen’s d=0.26, p=5.0e-14), with |d|=0.11–0.26. The report reads this as a consistent but modest environmental shift rather than strict habitat segregation, and PCA likewise showed SNIPE species broadly distributed with a detectable environmental bias. ₁₉ On this basis the report rejects random distribution (H0) for both taxonomy and environmental-niche association. It did not map any dimension, including A19, to a named environmental variable. It describes AlphaEarth dimensions as derived from satellite remote-sensing imagery at sample lat/lon (Goormaghtigh et al. 2025), capturing signals such as vegetation indices, temperature, precipitation and land cover. It lists AlphaEarth as 64-dimensional vectors for 83K genomes at 28.4% coverage, biased toward environmental isolates with lat/lon data. ₁₉ The corrected architecture is PF13250/DUF4041 plus PF13455/Mug113, not canonical GIY-YIG PF01541: zero surveyed clusters contained both DUF4041/PF13250 and PF01541, while 54/4,572 carried Mug113. ₁₉

Environmental functional associations

Across 51,748 soil samples and nine metals—copper, cobalt, chromium, nickel, zinc, lead, arsenic, cadmium, and mercury—the soil-metal analysis identified 2,355 significant COG–metal associations at FDR<0.05; COG means Cluster of Orthologous Groups functional category and FDR means false discovery rate. Chromium and lead produced the strongest signals, with ABC and RND transporters and biosynthesis genes prominent. In 7,566 samples matched to nearby KBase genomes within 10 km, copper associations were strongest positive for cell-division and nucleotide-transport categories BQ, FQ, and FK and strongest negative for energy-production categories DI and CE. These observational results refine environmental-metal inference toward possible energetic trade-offs, not demonstrated mechanisms. ₁₇ Many genomes within 10 km may not be co-located with the copper measurements. The report therefore requires sensitivity analyses at 5 km and 20 km before COG–copper associations can be attributed to the measured sites. ₁₇

A conditional db-RDA gave R²=0.799 and p=0.005 using 999 permutations after conditioning on batch and project effects; the unconditional metal-only R² was not reported. Biome-stratified PGLS found distinct metal–COG relationships in soil, marine, and wastewater, supporting environment-specific rather than universal responses. ₁₇ The 2,355 discoveries among 3,915 implied tests—nine metals × 435 COGs—represent a 60% discovery rate; co-contamination makes tests non-independent, so Benjamini–Hochberg FDR may be anti-conservative and true FDR higher. Effect sizes were not systematically reported. ₁₇

The global metal analysis supports environmental composition while refining the signal from broad prevalence to sampled biome and geographic differences. ₁₀ MicrobeAtlas found Levins’ B_std niche breadth phylogenetically conserved across 1,264 bacterial genera (Pagel’s λ=0.787, p=7.9×10^-102), and habitat range even more conserved (λ=0.909, p=1.4×10^-157). In 606 genera with AMR annotations, metal type diversity associated with B_std in PGLS (β=+0.021, SE=0.0056, p=1.5×10^-4), while total AMR cluster count and core AMR fraction were not significant. The association remained after pangenome coverage (β=+0.0204, SE=0.0057, p=3.4×10^-4) and genome size (β=+0.0218, SE=0.0061, p=3.6×10^-4), but strict 5% within-environment prevalence reduced the sample to 379 genera and gave β=+0.0166, p=0.092. ₁₁

AMR analyses extend embedding associations to resistance genes. Among 2,659 species with both AMR data and AlphaEarth embeddings, 52 of 64 embedding dimensions significantly correlated with AMR diversity (FDR, false discovery rate, < 0.05; top: A34, rho = +0.24). A Mantel test, which correlates two distance matrices, found that environmental distance (cosine distance of embeddings) predicted AMR profile distance, calculated as Bray-Curtis dissimilarity of mechanism fractions (r = 0.098, p = 0.001). The coupling was strongest for clinical species (r = 0.177), moderate for soil (r = 0.129) and weakest for aquatic species (r = 0.061). ₂₆ The report states that this continuous analysis confirms its discrete NB02 findings but adds limited interpretability, because embedding dimensions are opaque and genome coverage is 28%. ₂₆ Separately, among 2,684 species with ≥3 genomes and embeddings, environmental diversity correlated with AMR count (Spearman rho=0.466, p=1.6e-144). Those AMR genes were less core (rho=-0.173, p=1.8e-19). ₂₃ The strongest coupling in the clinical stratum sits in tension with the explorer’s expectation that embeddings are least informative for clinical isolates. It could instead reflect the epidemiological rather than ecological structure proposed by the reanalysis, which remains untested. These results also feed environmental-resistome. _26,1,2

Functional-dark-matter tests covered 151 accessory dark-gene clusters across 31 species; 10 of 137 environmental-category tests and 1 of 67 AlphaEarth tests were significant at FDR<0.05. Nitrogen-source carriers correlated with total nitrogen (ρ=+0.109, n=1,231, FDR=2.3e-4) and ammonium nitrogen (ρ=+0.231, n=1,230, FDR=8.0e-16), pH carriers with pH (ρ=+0.157, n=4,366, FDR=7.4e-25), and anaerobic carriers negatively with dissolved oxygen (ρ=-0.298, n=272, FDR=1.5e-6). ₄ The report places the strongest within-species signals in Pseudomonas and P. syringae, with NMDC community-level correlations across 6,365 metagenomic samples as corroboration. Its figure fig10_embedding_umap.png shows an AlphaEarth embedding UMAP of carrier versus non-carrier genomes. ₄

The discoveries digest reports a stronger embedding-derived association with pathway completeness. AlphaEarth niche breadth correlated with pathway completeness (r=0.392, p=7.1e-70). Embedding variance correlated more strongly (r=0.412, p=1.8e-77), and geographic distance alone more weakly (r=0.360, p=1.8e-58). ₂₂ The digest interprets this as species with broader ecological niches requiring more complete metabolic toolkits. That interpretation remains a hypothesis. The PHB analysis found that genome size correlated with niche breadth (rho=0.302). There, the partial rho for the PHB–niche-breadth association after controlling for genome size was -0.047. Genome size could therefore confound the completeness association as well. _22,12

The plant study found 5,341 of 5,671 eggNOG ortholog groups significant in raw plant-versus-non-plant tests, but only 50 enriched families retained significance after phylum control. COG3569 had Fisher OR=8.92 and phylogenetically controlled OR=6.01, q=7.7e-242, occurring in 54% of plant-associated versus 12% of non-plant species. Beneficial plant-growth-promotion clusters had a 64.6% core fraction versus 45.2% pathogenic and 46.8% genome-wide; U=83,567,419, p=3.38e-125, bootstrap 95% CI [0.089, 0.106]. However, ten Pseudomonas strains showed biocontrol traits largely accessory, so the result is taxonomically scoped. ₁₃

The proposed whole-genome test remains whether environmental samples with similar AlphaEarth embeddings share more accessory genes after controlling for phylogeny. The original median partial correlation was 0.0025, while the reanalysis found no stronger environmental than human-associated association; these are not directly comparable because the reanalysis used all embedded genomes, including up to 3,505 per species, whereas the original used diversity-maximizing downsampling with a maximum of 250 and different genome sets. _{1 2} The pitfalls audit requires genome_id rather than gtdb_taxonomy_id, shared representations such as COG, KEGG orthologs, or Pfam, genus checks, and assembly accessions such as GCF_*. ₁₈

Selected modules may retain environmental signal where whole-genome similarity does not. Across 15 phylogenetically stratified species, 44/95 module-species tests (46.3%) showed higher within-module co-occurrence than expected, with mean contig co-localization 0.769; packaging, lysis, and lysogenic regulation were more co-localized, integration and anti-defense more dispersed. ₁₄ Plant complementarity contradicted a strong co-occurrence model: corrected pair-level complementarity was slightly below random, Cohen’s d≈-0.4 and permutation p<0.001; prevalence-weighted d=-0.39. BacDive and GapMind agreed at 83.1%, but GapMind core completeness was 0% across compartments and the NMDC C-score analysis found 0 PGP-dominant, 3 pathogen-dominant, and 66 dual-or-mixed genera among 69. ₁₃

Additional contrasts and provenance limitations

Across 293,059 GTDB r214 genomes, only 7,995 (2.7%) had plant-associated annotations covering 1,136 species; endophytes comprised 29 species, and BacDive identified 2,482 plant-related strains but added zero genome upgrades beyond ncbi_env. ₁₃ The functional-dark-matter collection contained 37/48 organisms from Pseudomonadota, and its top 500 candidates included 417 Gammaproteobacteria but no Archaea, Actinobacteria, or Epsilonproteobacteria. ₄ In 2,759 NMDC prokaryote-targeted metagenome runs, environmental metadata predicted eukaryotic read fraction across studies, but the association was largely confounded by study or batch. ₂₅

Harvard Forest refines distance-decay with a controlled single-site contrast. In 42 biosamples from a 25-year +5°C warming experiment, direct-sample taxonomy reproduced Actinobacteria-up/Acidobacteria-down in organic soil, while mineral soil showed no phylum-level changes after FDR correction. PERMANOVA attributed 7.6% of variance to warming (p=0.069), 30.6% to soil horizon (p=0.0002), and 41% to treatment × horizon (p=0.0002). Carbon-cycling KOs comprised 5 of 57 heated-up DNA hits versus 428 of 12,806 other hits, Fisher OR=2.78, one-sided p=0.042; approximately 39% of DNA KOs were organic-only, mineral-only, or sign-flipping. ₇

Oak Ridge supports environmental context but refines laboratory-to-field transfer: across 108 enigma-coral groundwater sites, 16s-amplicon-sequencing detected 14 of 26 Fitness Browser genera; Sphingomonas, Pseudomonas, and Caulobacter occurred at 93%, 91%, and 82%. Five of 11 sufficiently prevalent genera had uranium associations after BH-FDR correction, in both directions, while aggregate laboratory metal-tolerance score had rho=0.503, p=0.095, n=12. ₈

The lanthanide atlas supports environmental association while refining its scope: xoxF occurred in 3,690 genomes versus 195 with mxaF, ratio 18.92:1, xoxF fraction 0.9498 with 95% CI [0.9425, 0.9558], and one-sided p=7.6×10^-22. Soil_sediment had 6.84% xoxF, OR=1.92, p_BH=6.1×10^-39; marine 4.76%, OR=1.31, p_BH=7.8×10^-7; REE-impacted samples 10.81%, OR=3.51, p_BH=0.082; REE-acid-mine-drainage contained only 37 MAGs. ₉ AlphaEarth covered 39.5% (1,457 / 3,690) of xoxF genomes. That is better than the 28% pangenome baseline but leaves 60.5% of xoxF carriers without environmental coordinates. Embedding-based niche analysis was scoped out of the report, although a coverage-restricted supplementary analysis is feasible. ₉

The soil-frontier analysis supports coverage bias while refining its direction. Across 5,441 soil samples with clay, mine proximity, nighttime lights, uranium, and functional gene counts, the proposed clay shield was not supported globally: Soil & Climate, Geochemical, and Industrial models had out-of-sample R²=-0.205 ± 0.197, -0.331 ± 0.071, and -0.221 ± 0.042. Low- and high-clay cross-validation R² were -0.268 and -0.292, difference 0.024, 95% CI [-0.423, 0.161]; clay importance was approximately 0.14 but did not improve high-clay prediction. ₁₆

The Genomic Discovery Index (GDI), OTU Richness / (Mean Genome Completeness + 1), gave forest 902.36, cropland 890.82, grassland 503.42, and wetland 525.13 at 1° spatial bins. Frontier areas with GDI>1000 had mean pH=6.74 versus 5.94 in mapped areas, a +0.8 pH-unit gap. GDI is a novel index without published precedent and can equal 902 even with zero genomes because richness and the denominator can be separated from actual completeness; forest and cropland differ by only 1.3% without uncertainty estimates. ₁₆ This supports genomic-under-representation but does not distinguish biological novelty from sequencing, assembly, annotation, spatial-shift, or outlier effects. ₁₆

Provenance limitations are pervasive. ncbi_env is entity–attribute–value data. Obtaining genome metadata requires three steps: join genome.ncbi_biosample_id to ncbi_env.accession, filter attribute_name to isolation_source, collection_date, geo_loc_name and host, and pivot from long to wide format (see cross-tenant-data-bridging). Missing classifications cannot be interpreted as biological absence. ₁₈ The global metal analysis had a 30.8% per-sample coordinate gap, although all 532 MAG grid cells represented in its coordinate dataset had ENA coverage. ₁₀ MicrobeAtlas found pangenome coverage correlated with metal type diversity at r=0.35, p=2.6×10^-19; single-genome genera had lower metal type counts (Mann–Whitney p=5.1×10^-10), only 16.8% of reads joined genus-level AMR data, and rarefaction retained a positive direction in all 200 iterations but only 57.5% were Bonferroni-significant at p<0.0083. ₁₁

The ecological-agora analysis found 12 of 32 genus-level mismatches caused by ENIGMA strain-name collisions, including one introducing 1,751 spurious clinical genomes; GTDB-Tk and genus checks reduced verified linkages from 32 to 20. Harvard Forest’s abiotic_features field was all zeros because of an NMDC parsing artifact, leaving +5°C as the only environmental contrast. Oak Ridge was limited by genus-level 16S resolution, temporal mismatch, aggregation, and uncontrolled pH, dissolved oxygen, carbon, and other covariates. _5,6,7,8

Other annotation limitations include 12 of 22 plant marker Pfams absent from bakta_pfam_domains despite abundance in interproscan_domains; the ten non-zero Bakta Pfams returned only 10–35% of InterProScan counts. Versioned identifiers recovered 19,364 hits across 7,962 species for ten marker Pfams. ₁₃ Keyword-based Pseudomonas classification from free-text isolation sources had approximately 6.7% unknown and approximately 29.1% other, and the report notes that such misclassification could attenuate true ecological signal; majority voting can obscure generalists, GapMind omits toluene, naphthalene, and benzoate degradation, and phylogenetic non-independence was not explicitly controlled. ₁₅

SNIPE annotations may miss divergent homologues; only 54/4,572 DUF4041 clusters had Mug113 co-annotation and DUF4041 may occur in non-SNIPE proteins. ₁₉ The planned COG V defense-gene-density join could not be completed because it required Spark Connect for a 93M × 132M row join; NMDC returned only 2 samples for the queried ecosystem format, and PhageFoundry had no Klebsiella fitness/TnSeq tables. ₁₉ The prophage bridge assumed genus-level conservation, had 87.2% median matching coverage, and only 28% of genomes had AlphaEarth embeddings. ₁₄ In the NMDC workflow, 1,067 of 2,759 runs were pooled from multiple biosamples, so biosample joins could create pseudo-replication; GOTTCHA2 supplied usable eukaryotic signal, whereas Kraken2 and Centrifuge were prokaryote-restricted in that deployment. ₂₅

Tensions

Embeddings show clear geographic distance-decay, yet environment similarity had median partial correlation 0.0025 with gene content. _1,3 The environmental-only reanalysis found median 0.051 versus 0.084 for human-associated species, U=1536, p=0.83, while the ecotype analysis found environmental median 0.0025 and phylogenetic median 0.0143. _2,3 These estimates cannot be reconciled by averaging: the studies used different genome sets, species filters, and downsampling. _2,1

The original and reanalysis partial-correlation medians also differ: 0.0025 for the original environment measure versus 0.081 across all 183 species in the reanalysis. _3,2 This contradicts any direct comparison of their absolute effect sizes, but the discrepancy is attributed to different genome sets, full-genome extraction, and downsampling procedures. ₂

The structured env_broad_scale field is cleaner but covers only 42%; keyword isolation_source harmonization covers 71% of genomes with a label but leaves 17% Other and 12.5% Unknown. ₁ Coordinates add a parallel tension because distance-decay may reflect environmental context, epidemiological structure, institutional clustering, or approximate coordinates. _1,2 The sources also disagree on which isolates the embedded subset favours. The pitfalls digest says valid coordinates are often missing, especially for clinical isolates. The discoveries digest says clinical isolates have good geographic metadata and disproportionately have embeddings. The SNIPE report describes coverage as biased toward environmental isolates with lat/lon data, and the prophage report describes bias toward clinically and environmentally well-sampled lineages. _18,22,19,14 The pitfalls digest itself gives 83,227 embedded genomes in one section and 83,287 in another. ₁₈

Pseudomonas carbon profiles were environmentally associated (p=0.006), but balanced accuracy was 0.408 +/- 0.169 and PCA primarily separated Pseudomonas s.s. from Pseudomonas_E. ₁₅ PHB showed raw breadth association rho=0.106, p=1.77×10^-06, but genome-size-controlled partial rho=-0.047, p=0.037. ₁₂

The soil-metal study creates a tension between strong conditional multivariate association and uncertain attribution: conditional db-RDA R²=0.799, p=0.005, but unconditional metal-only R² was not reported; chromium, copper, lead, and zinc co-vary, the 60% discovery rate may be inflated by correlated tests, and proximity within 10 km does not guarantee co-location. ₁₇ MicrobeAtlas produced β=+0.021, p=1.5×10^-4, but strict prevalence gave p=0.092 and groundwater-specific fold enrichment was null (rho=+0.042, p=0.242). ₁₁

Prophage module effects remained significant within genome-size quartiles, but genome size dominated PERMANOVA (F=212.99), calls were annotation-based, only 28% of genomes had embeddings, and NMDC burden was genus-inferred. ₁₄ Human-associated module enrichments contrasted with 0/500 significant TerL-lineage enrichments, supporting modular exchange but not independent adaptation of complete phage lineages. ₁₄ Genome-size control also behaves differently across traits. After genome-size control, prophage module count kept a partial Spearman rho=0.468 with AlphaEarth niche breadth, whereas the PHB–niche-breadth partial rho was -0.047. Genome size therefore cannot be treated as a uniform explanation of embedding-derived niche-breadth associations. _14,12

Plant compartment PERMANOVA had refined R²=0.071, while earlier R²=0.527 collapsed to 0.072 after removing genome-rich species. Plant-associated mobilome burden was higher (3.7 versus 2.8), but marker singletons were depleted (enrichment ratio 0.78), and corrected complementarity was negative (d≈-0.4). ₁₃

SNIPE was statistically detectable in 22/64 AlphaEarth dimensions, but its largest effect was d=0.26, AlphaEarth covered only 28.4% of genomes, and no dimension had a biological interpretation. ₁₉ Its 86.7% accessory-plus-singleton distribution supports mobile defense-island turnover, but only 54/4,572 DUF4041 clusters carried Mug113 and DUF4041 can occur outside SNIPE. ₁₉

The soil-frontier models all had negative out-of-sample R² and low- versus high-clay difference 0.024 with 95% CI [-0.423, 0.161]. ₁₆ This supports the conclusion that the clay shield was not demonstrated globally, but contradicts treating negative prediction as evidence of genuine biological unpredictability because spatial distributional shift, outlier leverage, and model failure were not separated. ₁₆ Forest GDI=902.36 and cropland GDI=890.82 were jointly highest, but the 1.3% difference lacks bootstrap confidence intervals and GDI can be high when completeness is zero. ₁₆

The ecological-agora, ENIGMA, Harvard Forest, Oak Ridge, and lanthanide results retain their respective tensions: cluster alignment is not causal; global Pseudomonas representation differs from ENIGMA’s environmental clade; treatment alone explained 7.6% (p=0.069) versus treatment × horizon 41%; laboratory metal tolerance was non-significant in field abundance; and REE-impacted samples did not pass FDR while REE-AMD contained only 37 MAGs. _5,6,7,8,9

The pitfalls audit found K=4 LDA ecotypes with a 33-species Tier-A list, but held-out-species Jaccard values were 0.230 for E1 and 0.064 for E3; independent within-substudy evidence reduced the list to 3 candidates. ₁₈ This contradicts treating the unvalidated list as stable and supports held-out-feature clustering and independent validation.

Open Directions

• Compare keyword isolation_source with env_broad_scale for the 34,800 and 76,295 covered genomes using confusion matrices, category-specific embedding distances, and alternative ENVO mappings; determine which schema preserves ecological separation with fewer ambiguous assignments. ₁
• Reclassify the 30,469 Suspicious cluster and 2,708 Low precision coordinates using isolation_source homogeneity, species composition, coordinate precision, and field-site evidence; test the 0.41-to-0.82 distance pattern after reclassification. ₁
• Reweight or downsample Human clinical, Human gut, Human other, Soil, Marine, Freshwater, Extreme, and Plant categories, then repeat geographic-distance and gene-content analyses. ₁
• Directly compare full-genome extraction with diversity-maximized downsampling, including genome count as a covariate, to resolve the 27× partial-correlation discrepancy; repeat environmental-only analyses and test whether median association exceeds 0.0025. _2,1
• Use the 13,381-genome metadata set with direct environmental metadata and alternative environmental embedding distances to test whether non-AlphaEarth representations improve gene-content associations after accounting for phylogenetic distance. ₃
• Combine direct environmental metadata with embedding distances and partial-correlation or nonlinear distance-based methods to determine whether coarse descriptors or linear modeling hide environmental effects. ₃
• Stratify the 172-species result by host-associated versus environmental lifestyle, replace collection coordinates with host, tissue, body-site, or other organism-proximal metadata where available, and test whether the environmental effect changes. ₃
• Partition gene clusters into COG functional categories, including V-Defense and L-Mobile, and test whether category-specific associations with environment exceed the whole-genome signal. ₃
• For species with identified ecotype clusters, compare gene content between clusters using organism-proximal environmental metadata to determine whether ecotype- or locus-level associations are missed by collection-site coordinates. ₃
• Refit ecotypes with held-out-feature clustering, leave-one-species-out refitting, pathway/EC-number clustering, and independent validation; test whether E1 and E3 Jaccard values remain below 0.3 and whether the independent list remains at 3 candidates. ₁₈
• Re-run UMAP and DBSCAN across alternative n_neighbors, min_dist, and eps values and test whether environment-dominated clusters persist in the original 64-dimensional space. ₁
• Correlate A00–A63 with latitude, temperature, precipitation, NDVI, and land cover; analyze the 3,838 NaN-containing records. ₁
• Re-test environmental associations with GroupKFold by study, within-study contrasts, pooled versus single-biosample NMDC runs, and carrier-versus-non-carrier controls. _25,4
• Reconcile ENIGMA and global profiles using GTDB-validated identifiers, ChEBI-ID condition canonicalization, and study/site-blocked validation. _6,18
• Reanalyze Harvard Forest with horizon-matched DNA, RNA, and metabolite contrasts, restored temperature, pH, and nitrogen covariates, and multi-site comparisons. ₇
• Match Oak Ridge assemblies to Fitness Browser organisms at species or strain level and test uranium-specific fitness while controlling for pH, redox, carbon, and date. ₈
• Reanalyze xoxF/mxaF with phylogeny-aware models, ORF-integrity checks, Bakta validation, and expanded REE-AMD sampling. ₉
• Complete the metal-resistance map, normalize hotspots by log(n_MAGs), check sample-accession prefixes, and stratify the 30.8% coordinate gap. ₁₀
• Reanalyze MicrobeAtlas with finer environments, phylogenetic mixed models, sampling-effort covariates, covered-read fractions, contamination metadata, and accession-prefix controls. ₁₁
• For PHB, map PF00561/PF07167, reconstruct a phaC tree, and apply phylogenetic logistic regression to distinguish transfer, lineage, genome-size, and environmental effects. ₁₂
• Extend plant analyses with phylogenetic GLMMs, broader tree coverage, genome-size matching, independent metagenomes/transcriptomes, and reaction-level complementarity after recovering the 62/322 lost genera. ₁₃
• Test prophage effects with geNomad or VIBRANT, intact-versus-domesticated classification, lineage and genome-size controls, and assembly-level validation of the NMDC bridge; resolve module-versus-lineage effects with TerL trees, host phylogeny, and defense-island analyses. ₁₄
• Extend Pseudomonas profiling with aromatic-degradation modules and PGLS or phylogenetic logistic regression; test within-species ecotypes in P. fluorescens and P. putida against kescience-fitnessbrowser. ₁₅
• Test SNIPE using PF13250/DUF4041 plus PF13455 rather than PF01541; quantify DUF4041-only false positives, search divergent homologues, and compare coverage- and phylogeny-matched species. ₁₉
• Test whether SNIPE predicts A19 and other AlphaEarth dimensions, pH, temperature, precipitation, land cover, or NMDC ecosystem metadata after controlling for lineage, genome size, sampling effort, and coordinate quality; determine whether the effect generalizes beyond n=1,069 and n=13,977. ₁₉
• Test the proposed SNIPE–ManYZ trade-off experimentally in Klebsiella using strain-specific mutant libraries or curated Tn-Seq data; current co-occurrence evidence cannot establish mechanism. ₁₉
• Decompose soil negative R² using spatial blocking, train/test distribution-shift diagnostics, fold-level MSE leverage checks, and outlier sensitivity analyses. ₁₆
• Recompute GDI after uniform 16S rarefaction, report richness and completeness separately, bootstrap biome-level 95% CIs, and control the +0.8 pH gap for sample number per pH bin. ₁₆
• Rerun soil-metal models as COG ~ Cr | Cu + Zn + Pb, report rho distributions and unconditional db-RDA R², assess dependence-aware FDR and effect sizes, and test copper proximity at 5 km, 10 km, and 20 km with Moran’s I and SEVM. ₁₇
• Classify significant metal-associated COGs into resistance, stress, membrane, energy, and unknown categories and compare soil, marine, and wastewater models to distinguish metal-specific from generic multi-metal contamination. ₁₇
• Build a provenance audit that discovers live KBase Data Lakehouse namespaces and schemas, validates tenant and identifier mappings, uses genome_id and sample_id bridge keys, and records missingness before environmental joins. ₁₈
• Run the coverage-restricted embedding niche analysis on the 1,457 AlphaEarth-covered xoxF genomes against matched non-carriers to test whether lanthanide methylotrophy tracks environmental context beyond categorical biomes. ₉
• Test the ENIGMA pH partition directly by measuring pH-dependent growth of Cluster A and Cluster B strains and asking whether cluster-distinguishing genomic features predict pH tolerance under genus-blocked validation. ₆
• Execute the unrun correlation of metabolic ecotypes with AlphaEarth niche breadth on the embedded subset, reporting sample size and coverage bias explicitly. ₂₁
• Resolve the 61 “Unknown” NMDC ecosystem samples with NMDC ENVO annotations or study-level metadata and re-run the habitat comparison to test whether finer soil-subtype structure emerges. ₂₀
• Map PHB+ and PHB− species across all 64 AlphaEarth dimensions, and repeat the prophage module-count and PHB niche-breadth partial correlations with identical genome-size and phylogenetic controls to explain their divergent post-control behaviour. _12,14
• Re-run the agora recent-gain cluster contrast with phylogenetic and genome-size controls to test whether mixed-biome clusters 5 and 8 genuinely carry higher recent-gain density than cluster 6. ₅