Comparative characterization of the mucosal eukaryotic virome across independent whole-metagenome shotgun sequencing cohorts
Abstract
The human mucosal eukaryotic virome remains incompletely characterized because methodological heterogeneity limits direct comparisons across independent whole-metagenome shotgun sequencing (WMS) datasets. To address this limitation, we reanalyzed publicly available DNA-based WMS datasets representing healthy body-site samples together with nasopharyngeal carcinoma (NPC), human papillomavirus (HPV)-associated disease, and oral squamous cell carcinoma (OSCC) using the same ViromeScan-based analytical approach. Viral count profiles were converted to within-sample relative abundances for Bray–Curtis analyses, host-unmapped read depth was included as a covariate in PERMANOVA, and alpha-diversity findings were evaluated using repeated rarefaction. Anatomical site remained associated with eukaryotic virome diversity and community structure in the HMP cohort after accounting for sequencing depth, and site-associated differences were also observed in the independent NPC cohort analyzed using the same approach. HPV infection status was also associated with viral community structure after depth adjustment, whereas DOI-defined OSCC groups showed no significant global differences in alpha diversity, normalized viral recovery, or community structure and yielded only exploratory taxon-level signals. Across datasets, viral recovery varied with host-unmapped sequencing depth, with substantial between-sample heterogeneity. Together, these findings show that consistent application of publicly available analytical tools can support comparative characterization of eukaryotic virome profiles from existing WMS datasets, while also highlighting the importance of sequencing depth and cautious interpretation of taxon-level associations.
Related
- NanoPAMDA enables rapid PAM characterization of metagenome-derived Cas9 effectors
- Integrating multiplex digital PCR and long-read sequencing toward standardized genome integrity and purity characterization of rAAV9.
- An ONT sequencing framework for hantavirus genome assembly, variant analysis, and defective viral genome characterization from whole-segment amplicons
- RNA sequencing of antennae and comparative analysis of key gene families in Apis mellifera workers exhibiting drastic differences in propolis yield
- Transcriptomic Characterization of Terminal Complement Complex-bound Cells in Human Choroid Using Single-Cell RNA Sequencing
BibTeX
@article{na2026comparative,
title = {{Comparative characterization of the mucosal eukaryotic virome across independent whole-metagenome shotgun sequencing cohorts}},
author = {Hee Sam Na and Jung Hwa Park and Tae Sung Kim and Na Rae Choi},
journal = {Scientific Reports},
year = {2026},
month = oct,
doi = {10.1038/s41598-026-73604-6},
url = {https://www.nature.com/articles/s41598-026-73604-6},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Comparative characterization of the mucosal eukaryotic virome across independent whole-metagenome shotgun sequencing cohorts
AU - Hee Sam Na
AU - Jung Hwa Park
AU - Tae Sung Kim
AU - Na Rae Choi
PY - 2026
DA - 2026/10/01
JO - Scientific Reports
DO - 10.1038/s41598-026-73604-6
UR - https://www.nature.com/articles/s41598-026-73604-6
AB - The human mucosal eukaryotic virome remains incompletely characterized because methodological heterogeneity limits direct comparisons across independent whole-metagenome shotgun sequencing (WMS) datasets. To address this limitation, we reanalyzed publicly available DNA-based WMS datasets representing healthy body-site samples together with nasopharyngeal carcinoma (NPC), human papillomavirus (HPV)-associated disease, and oral squamous cell carcinoma (OSCC) using the same ViromeScan-based analytical approach. Viral count profiles were converted to within-sample relative abundances for Bray–Curtis analyses, host-unmapped read depth was included as a covariate in PERMANOVA, and alpha-diversity findings were evaluated using repeated rarefaction. Anatomical site remained associated with eukaryotic virome diversity and community structure in the HMP cohort after accounting for sequencing depth, and site-associated differences were also observed in the independent NPC cohort analyzed using the same approach. HPV infection status was also associated with viral community structure after depth adjustment, whereas DOI-defined OSCC groups showed no significant global differences in alpha diversity, normalized viral recovery, or community structure and yielded only exploratory taxon-level signals. Across datasets, viral recovery varied with host-unmapped sequencing depth, with substantial between-sample heterogeneity. Together, these findings show that consistent application of publicly available analytical tools can support comparative characterization of eukaryotic virome profiles from existing WMS datasets, while also highlighting the importance of sequencing depth and cautious interpretation of taxon-level associations.
SN - 2045-2322
ER -