Estimating Gene Prevalence with Incomplete Genomes: Application to the Prochlorococcus Pangenome
- 1University of Washington
Abstract
Gene prevalence is an essential statistic in population genomics, defined as the proportion of a population carrying at least one ortholog of a given gene family. Common approaches to estimating gene prevalence either assume complete genomes or exclude genomes below a completeness threshold, resulting in biased and high variance estimates that fail to take advantage of all available data. Here we developed a maximum-likelihood framework that explicitly models genome completeness in order to estimate gene prevalence and associated confidence intervals from collections of complete and incomplete genomes. We applied the methodology to marine Prochlorococcus, a globally abundant photosynthetic bacterium for which the majority of available genomes are incomplete. Our methodology produced gene prevalence estimates with associated confidence intervals for a total of 52,500 gene families, more than 4 times the number quantifiable using complete genomes alone. We find a core of 714 genes are likely universal to the Prochlorococcus genus and a substantial majority of 49,828 gene families occur at low prevalence. Our results reveal a Prochlorococcus pangenome dominated by rare flexible genes and provide a general framework for incorporating incomplete genomes into quantitative analyses of microbial genomic diversity.
Related
- Estimating the Prevalence of Malignant Hyperthermia Susceptibility Based on Causative RYR1 and CACNA1S Variants
- A meta-analysis-guided sampling plan for estimating the prevalence of 10 zoonotic and 4 transboundary animal diseases across India
- Estimating marginal association between paired outcomes in clustered data with an application to analyzing caries and periodontal disease data from the National Health and Nutrition Examination Survey
- Estimating Three‐Dimensional Posture to Correct Drone‐Based Measurements of Free‐Swimming Marine Mammals: An Application to Dugongs in Shallow Waters
- Defining traits of low-light adapted Prochlorococcus inhabiting surface waters of the Equatorial Pacific Ocean
BibTeX
@misc{blaskowski2026estimating,
title = {{Estimating Gene Prevalence with Incomplete Genomes: Application to the Prochlorococcus Pangenome}},
author = {S. Blaskowski and L. G. Christie and Y. Raut and J. Fuhrman and J. Bien and E. V. Armbrust},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.19.760200},
url = {https://www.biorxiv.org/content/10.64898/2026.10.19.760200v1},
note = {CC BY-NC},
}
RIS
TY - UNPB
TI - Estimating Gene Prevalence with Incomplete Genomes: Application to the Prochlorococcus Pangenome
AU - S. Blaskowski
AU - L. G. Christie
AU - Y. Raut
AU - J. Fuhrman
AU - J. Bien
AU - E. V. Armbrust
PY - 2026
DA - 2026/10/19
PB - bioRxiv
DO - 10.64898/2026.10.19.760200
UR - https://www.biorxiv.org/content/10.64898/2026.10.19.760200v1
AB - Gene prevalence is an essential statistic in population genomics, defined as the proportion of a population carrying at least one ortholog of a given gene family. Common approaches to estimating gene prevalence either assume complete genomes or exclude genomes below a completeness threshold, resulting in biased and high variance estimates that fail to take advantage of all available data. Here we developed a maximum-likelihood framework that explicitly models genome completeness in order to estimate gene prevalence and associated confidence intervals from collections of complete and incomplete genomes. We applied the methodology to marine Prochlorococcus, a globally abundant photosynthetic bacterium for which the majority of available genomes are incomplete. Our methodology produced gene prevalence estimates with associated confidence intervals for a total of 52,500 gene families, more than 4 times the number quantifiable using complete genomes alone. We find a core of 714 genes are likely universal to the Prochlorococcus genus and a substantial majority of 49,828 gene families occur at low prevalence. Our results reveal a Prochlorococcus pangenome dominated by rare flexible genes and provide a general framework for incorporating incomplete genomes into quantitative analyses of microbial genomic diversity.
ER -