Taxonomic and functional concordance between full-length ONT 16S and ONT shotgun metagenomics in the canine gut microbiome
- 1Szegedi Tudomanyegyetem Altalanos Orvostudomanyi Kar
Abstract
Background: Full-length Oxford Nanopore Technologies (ONT) 16S rRNA sequencing provides a scalable view of microbial community composition and can support phylogeny-based functional prediction, but it is not equivalent to shotgun metagenomics. We asked which biological conclusions are preserved when the same canine fecal specimens are profiled by full-length ONT 16S and ONT whole-genome shotgun (WGS) sequencing, and how their agreement depends on analytical scale, reference representation and classifier. Methods: Ninety-seven fecal specimens from 51 dogs were profiled with both assays from the same DNA extract. Functional profiles predicted from NanoASV/NanoPredict with PICRUSt2 were compared with WGS-supported KEGG Ortholog (KO) profiles generated by Kadath. Taxonomy was benchmarked in a source-genome-matched RefSeq universe and in a host-specific DogMAG universe using minitax and Kraken2. Agreement was evaluated at whole-profile, feature-abundance, detection, between-sample structure and biological-inference scales. Age-associated transfer was assessed with dog-aware continuous mixed models, grouped signed-score analyses and paired/dog-blocked PERMANOVA. Results: Functional whole-profile concordance was high: median within-sample CLR Spearman correlations ranged from 0.781 to 0.860 across developmental strata, while between-sample functional structure remained significant by Mantel (rho=0.543) and Procrustes (r=0.693; both p=0.001). Feature-wise transfer was substantially weaker (median KO-wise CLR Spearman=0.318). Continuous age-associated KO slopes showed substantial cross-assay concordance (Spearman=0.727; signed-score Spearman=0.753; direction agreement=77.9%), although 1,290/5,258 eligible KOs retained significant assay-by-age interactions. Taxonomically, exact genus/species abundance agreement was much lower than agreement in between-sample ecological structure. Host-specific DogMAG improved species-level median Spearman from 0.261 to 0.656 for minitax SpeciesEstimate and from 0.181 to 0.512 for Kraken2. The classifier effect was independent of reference choice: under both RefSeq and DogMAG, minitax yielded stronger 16S-WGS concordance than Kraken2, with all eight prespecified RefSeq paired genus/species endpoints and all 10 DogMAG primary paired endpoints significant after BH correction. The same ordering extended to developmental inference, with DogMAG genus/species age-slope concordance of 0.795/0.799 for SpeciesEstimate versus 0.693/0.702 for Kraken2. Taxonomic Aitchison PERMANOVA detected age-associated structure in every assay/reference/classifier/rank combination, whereas age-by-assay interactions were consistently significant but small (R2 approximately 1.1 to 2.2%). Stricter NanoASV identity thresholds removed substantial 16S abundance without improving species-level agreement. Conclusions: The extent of cross-assay agreement depends on the level of analysis. Full-length ONT 16S preserves broad functional organization, ecological structure and much of the direction of age-associated change, but exact fine-rank composition, individual-feature abundance and effect magnitude remain assay dependent. Host-specific reference representation substantially narrows the taxonomic gap, and classifier choice exerts an additional independent effect: within the same matched reference set, minitax consistently yields stronger 16S-WGS concordance than Kraken2 across abundance, detection, ecological-distance and developmental-inference endpoints. Full-length ONT 16S is therefore well suited to broad ecological screening and hypothesis generation, whereas WGS remains preferable when conclusions depend on quantitative fine-rank composition, directly supported gene content or precise feature-level effect estimates.
Related
- Genome-resolved metagenomics reveals the taxonomic and functional landscape of the rumen microbiome in Japanese Black × Holstein cattle and its associations with carcass traits
- Long-read shotgun metagenomics identifies distinct urogenital microbial community states associated with functional potential and urinary tract infection during pregnancy
- Integrated gene-catalog and genome-resolved metagenomics reveal taxonomic and metabolic differences in the gut microbiome in irritable bowel syndrome
- Metagenomics identifies distinct functional and taxonomic profiles between freshwaters and associated sediments in alpine habitats.
- Comparative Shotgun Metagenomics of Microbial Communities Associated with Two Geothermal Systems in Kazakhstan.
BibTeX
@misc{yaseen2026taxonomic,
title = {{Taxonomic and functional concordance between full-length ONT 16S and ONT shotgun metagenomics in the canine gut microbiome}},
author = {N. J. Yaseen and B. Kakuk and G. Gulyas and M. Asaduzzaman and Z. Boldogkoi and D. Tombacz},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.26.754683},
url = {https://www.biorxiv.org/content/10.64898/2026.09.26.754683v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Taxonomic and functional concordance between full-length ONT 16S and ONT shotgun metagenomics in the canine gut microbiome
AU - N. J. Yaseen
AU - B. Kakuk
AU - G. Gulyas
AU - M. Asaduzzaman
AU - Z. Boldogkoi
AU - D. Tombacz
PY - 2026
DA - 2026/09/26
PB - bioRxiv
DO - 10.64898/2026.09.26.754683
UR - https://www.biorxiv.org/content/10.64898/2026.09.26.754683v1
AB - Background: Full-length Oxford Nanopore Technologies (ONT) 16S rRNA sequencing provides a scalable view of microbial community composition and can support phylogeny-based functional prediction, but it is not equivalent to shotgun metagenomics. We asked which biological conclusions are preserved when the same canine fecal specimens are profiled by full-length ONT 16S and ONT whole-genome shotgun (WGS) sequencing, and how their agreement depends on analytical scale, reference representation and classifier. Methods: Ninety-seven fecal specimens from 51 dogs were profiled with both assays from the same DNA extract. Functional profiles predicted from NanoASV/NanoPredict with PICRUSt2 were compared with WGS-supported KEGG Ortholog (KO) profiles generated by Kadath. Taxonomy was benchmarked in a source-genome-matched RefSeq universe and in a host-specific DogMAG universe using minitax and Kraken2. Agreement was evaluated at whole-profile, feature-abundance, detection, between-sample structure and biological-inference scales. Age-associated transfer was assessed with dog-aware continuous mixed models, grouped signed-score analyses and paired/dog-blocked PERMANOVA. Results: Functional whole-profile concordance was high: median within-sample CLR Spearman correlations ranged from 0.781 to 0.860 across developmental strata, while between-sample functional structure remained significant by Mantel (rho=0.543) and Procrustes (r=0.693; both p=0.001). Feature-wise transfer was substantially weaker (median KO-wise CLR Spearman=0.318). Continuous age-associated KO slopes showed substantial cross-assay concordance (Spearman=0.727; signed-score Spearman=0.753; direction agreement=77.9%), although 1,290/5,258 eligible KOs retained significant assay-by-age interactions. Taxonomically, exact genus/species abundance agreement was much lower than agreement in between-sample ecological structure. Host-specific DogMAG improved species-level median Spearman from 0.261 to 0.656 for minitax SpeciesEstimate and from 0.181 to 0.512 for Kraken2. The classifier effect was independent of reference choice: under both RefSeq and DogMAG, minitax yielded stronger 16S-WGS concordance than Kraken2, with all eight prespecified RefSeq paired genus/species endpoints and all 10 DogMAG primary paired endpoints significant after BH correction. The same ordering extended to developmental inference, with DogMAG genus/species age-slope concordance of 0.795/0.799 for SpeciesEstimate versus 0.693/0.702 for Kraken2. Taxonomic Aitchison PERMANOVA detected age-associated structure in every assay/reference/classifier/rank combination, whereas age-by-assay interactions were consistently significant but small (R2 approximately 1.1 to 2.2%). Stricter NanoASV identity thresholds removed substantial 16S abundance without improving species-level agreement. Conclusions: The extent of cross-assay agreement depends on the level of analysis. Full-length ONT 16S preserves broad functional organization, ecological structure and much of the direction of age-associated change, but exact fine-rank composition, individual-feature abundance and effect magnitude remain assay dependent. Host-specific reference representation substantially narrows the taxonomic gap, and classifier choice exerts an additional independent effect: within the same matched reference set, minitax consistently yields stronger 16S-WGS concordance than Kraken2 across abundance, detection, ecological-distance and developmental-inference endpoints. Full-length ONT 16S is therefore well suited to broad ecological screening and hypothesis generation, whereas WGS remains preferable when conclusions depend on quantitative fine-rank composition, directly supported gene content or precise feature-level effect estimates.
ER -