Self-Sufficiency, Anaerobic Toolkit, and Cultivation Bias in Clay-Confined Cultured Bacterial Genomes

Recent literature on the deep terrestrial subsurface (Beaver & Neufeld 2024 review; Bagnoud 2016 Mont Terri Opalinus; Mitzscherling 2023 Opalinus rock-attached communities; Engel 2019 Grimsel bentonite) generates testable predictions for what cultured bacterial genomes from clay-confined sites should encode: greater biosynthetic completeness ("self-sufficiency"), a recurring anaerobic-respiration toolkit, and — because cultivable porewater organisms differ qualitatively from rock-attached Geobacter/Geothrix lineages — a cultivation-driven porewater bias in any genome-resolved cohort. We test these predictions on the ~25–40 BERDL cultured genomes traceable to clay-confined origins (Mont Terri Opalinus boreholes, bentonite formations, kaolin lenses, BacDive deep-clay strains) versus a phylogenetically matched soil/sediment baseline drawn from the 5,151 species linked to soil biosamples in kbase_ke_pangenome.ncbi_env.

This project is the genome-resolved complement to enigma_sso_asv_ecology, which characterized in-situ subsurface community structure via 16S ASVs at Oak Ridge.

Cohort

The cohort was assembled from kbase_ke_pangenome.ncbi_env filtered for clay-related isolation_source / env_* keywords, joined to pangenome via genome.ncbi_biosample_id. Final cohort sizes after CheckM ≥80% / ≤5% filtering:

Anchor_deep genera include exactly the lineages Bagnoud (2016) identified as recurrent across 7 Mont Terri boreholes: Desulfosporosinus (×2), BRH-c8a (Peptococcaceae c8a in Bagnoud's nomenclature; ×2), BRH-c4a (Desulfotomaculales), Lutibacter + BRH-c54 (Bacteroidota), Roseovarius (Rhodobacterales), and Stenotrophomonas (the bentonite isolate).

H1 — Self-sufficiency (NB04)

Mann–Whitney U two-sided. Effect size = Cohen's d.

H2 — Anaerobic toolkit (NB05)

Per-marker Fisher (anchor_deep vs soil_baseline, BH-FDR adjusted):

Within-phylum control (Bacillota_B):

H3 — Porewater-vs-rock-attached signature (NB06)

Marker class breakdown (CheckM-filtered):

Tests against Mitzscherling (2023) rock-attached null (SR ~0.2%, IR ~7%):

Pairwise cohort Fisher (SR_complete only):

Literature Context

• The H3 result directly replicates the Bagnoud (2016) Mont Terri porewater paradigm at the population-genome level. Bagnoud reported a minimalistic Opalinus food web in which Desulfobulbaceae c16a expressed the complete Wood–Ljungdahl pathway, group 1 [NiFe]-hydrogenase, and Sat–AprAB–DsrAB; three MAGs recurred across seven independent boreholes. Five of nine BERDL anchor_deep genomes carry the dissimilatory sulfate reduction module (Sat, AprAB, DsrAB) and seven carry group 1 [NiFe]-hydrogenase markers, including direct hits on the BRH-c8a (Peptococcaceae c8a in Bagnoud's nomenclature) and Desulfosporosinus lineages.

• The H3 result also directly diverges from the Mitzscherling et al. (2023) rock-attached community profile (SRB <0.2%, IRB 4.3–10.2% dominated by Geobacter and Geothrix). Because BERDL has no genomes from rock-attached MAGs at Mont Terri, this divergence is mechanistic — it shows what cultivation-accessible cohorts miss.

• The H2 within-phylum decomposition matches Beaver & Neufeld's (2024) prediction that "Bacillota dominate deeper isolated fluids ... in part because they favor the reductive acetyl-CoA (Wood–Ljungdahl) pathway and form spores." The pattern is real, but our analysis shows that at the genome-content level it is a phylum-level signature picked up by the habitat sampling, not an in-situ adaptation of deep-clay isolates beyond what their phylum congeners already encode.

• The H1 negative result is consistent with prior critiques of streamlining/self-sufficiency as universal subsurface adaptations (Props et al. 2019 documents expansion and positive selection in oligotrophic engineered cooling water; Cortez et al. 2022 shows streamlining is acidophile-specific). Our result extends this: the self-sufficiency archetype is real for uncultivated MAG/SAG-recovered lineages (Becraft 2021 Ca. Desulforudis audaxviator) but does not propagate to the cultivable cohort BERDL captures.

Novel Contribution

1. First quantitative test of the porewater-vs-rock-attached signature dichotomy in cultured pangenome data. Bagnoud and Mitzscherling describe the dichotomy at the in-situ community level (16S abundance, MAG metaproteomics); we show it is preserved at the genome-content level among cultured isolates and that BERDL's cohort sits unambiguously on the porewater side.

2. Direct measurable cultivation bias diagnostic. Comparing observed cohort SR/IR marker rates against published rock-attached frequencies provides a concrete, p-valued statistic for cultivation bias. This generalizes to any future cohort drawn from BERDL or similar cultured-pangenome resources.

3. Phylogenetic decomposition of the anaerobic toolkit signal. Cohort-level toolkit enrichment (1.89 vs 0.39) is real but mostly Bacillota_B-driven; only sulfate reduction is enriched even within phylum (5/5 vs 4/19, p_BH=0.04). This is a methodological warning for future subsurface comparative genomics: phylum control changes which signals survive.

4. Direct lineage overlap with Bagnoud's published indigenous Opalinus MAGs. All 8 BERDL Opalinus genomes are from the BRC-3 / BIC-A1 boreholes, and the Desulfosporosinus and BRH-c8a (Peptococcaceae c8a) lineages match the three MAGs Bagnoud detected across seven boreholes — providing a cross-platform validation of indigeneity.

Limitations

• Small anchor cohort. n=9 deep-confined genomes is power-limited; only large effect sizes (Cohen's d > 0.7) are reliably detectable in unfiltered comparisons. Reported p-values for marginal effects (e.g., H1 within-Bacillota_B p=0.07) should be treated as descriptive.

• Cultivation bias is the headline finding, not a confounder. Because the cohort is cultured-only, all conclusions apply to "cultivable porewater-cultured deep-clay isolates," not to the full Mont Terri / bentonite microbial community. CPR / DPANN episymbionts (Bell 2022) and rock-attached Geobacter / Geothrix lineages are essentially absent. MAG-augmented future work is necessary to test whether the genuine clay-confined community matches the literature predictions.

• Compartment annotation is text-based. Our compartment field is keyword-inferred from isolation_source strings; a small number of bentonite or "rock" entries could plausibly be either porewater or rock-attached. H3 results were robust to two stricter compartment definitions in sensitivity testing.

• GapMind ceiling effect. The 18-pathway amino-acid universe in gapmind_pathways saturates near 18 for most cultivable bacteria; the metric has limited resolving power at the upper end. A finer-grained self-sufficiency metric (e.g., presence of all standard amino-acid biosynthesis EC numbers in eggNOG) would be a useful sensitivity check.

• eggNOG cluster-level annotations propagate within ≥90% AAI clusters. Strain-level marker variants (e.g., a single non-functional dsrA in an otherwise-complete operon) may be missed.

1. MAG-augmented expansion. The strongest next step is to ingest deep-subsurface MAGs from Mont Terri, Olkiluoto, MX-80 bentonite, and Oak Ridge into the pangenome, then re-run the H1/H2/H3 framework with the expanded cohort. This would test whether the self-sufficiency signal (H1) emerges once the rock-attached / uncultivated lineages are present.

2. Direct sub-cohort comparison: BRC-3 porewater vs BIC-A1 borehole. With 5 + 3 genomes per borehole, an exploratory within-Mont-Terri comparison could test whether the Bagnoud-paradigm SR-rich pattern is borehole-specific or site-wide.

3. Apply the porewater-bias diagnostic to other subsurface cohorts. Granite-hosted (Olkiluoto), basalt-hosted (Oak Ridge), and salt-cavern subsurface cohorts in BERDL or future ingests can be rapidly tested with the same SR/IR marker framework.

4. Genus-level analysis within Bacillota_B. The 5/5 Bacillota_B SR enrichment is striking but the "n=5 vs 19" structure may hide further granularity (Desulfosporosinus vs BRH-c8a vs BRH-c4a may differ in which non-SR features they carry). A within-Bacillota_B genus-level pangenome analysis could surface the deep-clay-specific accessory genome.

5. Cross-link to Bagnoud's metaproteomics evidence. Bagnoud (2016) reported protein-level expression of the toolkit modules. For BERDL Opalinus genomes that match Bagnoud's MAGs by ANI, a future project could ask whether the gene presence we observe corresponds to expression-validated activity.

Sources

Generated Data

Full bibliography in references.md. Primary citations supporting findings:

• Bagnoud A et al. (2016). Reconstructing a hydrogen-driven microbial metabolic network in Opalinus Clay rock. Nat Commun 7:12770. PMID: 27739431.
• Beaver RC, Neufeld JD (2024). Microbial ecology of the deep terrestrial subsurface. ISME J 18(1):wrae091. PMID: 38780093.
• Becraft ED et al. (2021). Evolutionary stasis of a deep subsurface microbial lineage. ISME J 15(10):2830–2842. PMID: 33824425.
• Bell E et al. (2022). Active anaerobic methane oxidation and sulfur disproportionation in the deep terrestrial subsurface. ISME J 16(5):1583–1593. PMID: 35173296.
• Beller HR et al. (2012). Genomic and physiological characterization of the chromate-reducing, aquifer-derived Firmicute Pelosinus sp. strain HCF1. Appl Environ Microbiol 78(24):8791–8800. PMID: 23064329.
• Cortez D et al. (2022). A large-scale genome-based survey of acidophilic bacteria suggests genome streamlining is an adaptation for life at low pH. Front Microbiol 13:803241. PMID: 35387071.
• Engel K et al. (2019). Stability of microbial community profiles associated with compacted bentonite from the Grimsel Underground Research Laboratory. mSphere 4(6):e00601-19. PMID: 31852805.
• Mitzscherling J et al. (2023). Clay-associated microbial communities and their relevance for a nuclear waste repository in the Opalinus Clay rock formation. MicrobiologyOpen 12(4):e1370. PMID: 37642485.
• Props R et al. (2019). Gene expansion and positive selection as bacterial adaptations to oligotrophic conditions. mSphere 4(1):e00011-19. PMID: 30728279.