Predicting cyanobacterial oxic diazotrophy and prioritizing candidate protein families beyond known Nif proteins
Abstract
Diazotrophic cyanobacteria must reconcile oxygenic photosynthesis with the oxygen sensitivity of nitrogenase. Engineering studies have demonstrated nitrogenase activity in a non-diazotrophic cyanobacterium while highlighting the importance of accessory functions and oxygen management. Robust aerobic nitrogen fixation remains challenging, motivating systematic identification of candidate support functions beyond the canonical Nif machinery. We modeled nifHDK-defined diazotrophic potential from protein-family presence or absence in 426 cyanobacterial assemblies. Four classifiers evaluated with genus-blocked cross-validation achieved mean ROC AUCs of 0.853 to 0.893, with purity and genus-breadth filters fitted within training folds. Restricting the panel to 401 assemblies with current NCBI status gave AUCs of 0.806 to 0.881. Family-content logistic regression outperformed a genome-size-only baseline across ten genus partitions. Two inventories describe candidates across the full panel: 981 high-purity families and 473 families with positive residual associations (accompanied by three ubiquitous controls). Model importance orders the residual-association inventory. A separate heuristic score combines inventory membership with proteomics, condensate-driver ranks and genus-based morphotype proxies to prioritize 167 families for experiments. Proteomics annotations tracked detection opportunity; the condensate contrast between inventories attenuated after prevalence adjustment. These annotations help choose experiments but do not independently validate accessory functions. Protein-family content predicts nifHDK-defined diazotrophic potential, with performance sensitive to assembly inclusion, taxonomy, genome size and taxonomic weighting. The two candidate inventories distinguish high carrier purity from residual association, providing complementary resources for prioritizing previously uncharacterized functions associated with cyanobacterial oxic diazotrophy.
Related
- Predicting taxon-specific benthic cyanobacterial mat cover and anatoxin concentrations in northern California rivers
- RevPert: predicting candidate drivers of transcriptomic state transitions via gallery-native reverse perturbation
- ASOFormer: a Transformer-based model for predicting antisense oligonucleotide efficacy to support therapeutic candidate prioritization
- Plant–cyanobacterial symbioses: mutualism with emerging translational possibilities
- Title: The Cyanobacterial dnaX Intein Encodes an Out-of-Frame Homing Endonuclease and is Sporadically Distributed in the Phylum
BibTeX
@misc{young2026predicting,
title = {{Predicting cyanobacterial oxic diazotrophy and prioritizing candidate protein families beyond known Nif proteins}},
author = {J. Young and S. Aryal and D. Nelson and D. Raya and E. A. Gnimpieba Z. and R. Zhou},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.11.758377},
url = {https://www.biorxiv.org/content/10.64898/2026.10.11.758377v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Predicting cyanobacterial oxic diazotrophy and prioritizing candidate protein families beyond known Nif proteins
AU - J. Young
AU - S. Aryal
AU - D. Nelson
AU - D. Raya
AU - E. A. Gnimpieba Z.
AU - R. Zhou
PY - 2026
DA - 2026/10/11
PB - bioRxiv
DO - 10.64898/2026.10.11.758377
UR - https://www.biorxiv.org/content/10.64898/2026.10.11.758377v1
AB - Diazotrophic cyanobacteria must reconcile oxygenic photosynthesis with the oxygen sensitivity of nitrogenase. Engineering studies have demonstrated nitrogenase activity in a non-diazotrophic cyanobacterium while highlighting the importance of accessory functions and oxygen management. Robust aerobic nitrogen fixation remains challenging, motivating systematic identification of candidate support functions beyond the canonical Nif machinery. We modeled nifHDK-defined diazotrophic potential from protein-family presence or absence in 426 cyanobacterial assemblies. Four classifiers evaluated with genus-blocked cross-validation achieved mean ROC AUCs of 0.853 to 0.893, with purity and genus-breadth filters fitted within training folds. Restricting the panel to 401 assemblies with current NCBI status gave AUCs of 0.806 to 0.881. Family-content logistic regression outperformed a genome-size-only baseline across ten genus partitions. Two inventories describe candidates across the full panel: 981 high-purity families and 473 families with positive residual associations (accompanied by three ubiquitous controls). Model importance orders the residual-association inventory. A separate heuristic score combines inventory membership with proteomics, condensate-driver ranks and genus-based morphotype proxies to prioritize 167 families for experiments. Proteomics annotations tracked detection opportunity; the condensate contrast between inventories attenuated after prevalence adjustment. These annotations help choose experiments but do not independently validate accessory functions. Protein-family content predicts nifHDK-defined diazotrophic potential, with performance sensitive to assembly inclusion, taxonomy, genome size and taxonomic weighting. The two candidate inventories distinguish high carrier purity from residual association, providing complementary resources for prioritizing previously uncharacterized functions associated with cyanobacterial oxic diazotrophy.
ER -