Back

Metatranscriptomic Insights into Microbial Dynamics Prior to Disease Onset in Atlantic Salmon Aquaculture

Harding, A.; Diblasi, C.; Manousi, D.; Soung Kwak, J.; Hazlerigg, D. G.; Morkore, T.; Snippen, L. G.; de Leon, A. V.-P.; Saitou, M.

2025-04-15 microbiology
10.1101/2025.04.15.648526 bioRxiv
Show abstract

Poly(A)-selected RNA-seq datasets are routinely generated in aquaculture research, yet the microbial information contained in unmapped reads is seldom explored due to the low abundance of nonhost transcripts and concerns about contamination. In this study, we repurposed Atlantic salmon gill RNA-seq data to assess whether meaningful microbial signals can be recovered using a contamination-aware and compositionally appropriate framework. Unmapped reads were analyzed with a custom Kraken2 database composed exclusively of complete, circularized salmon-associated bacterial genomes together with all available Atlantic salmon assemblies and the human genome. Although microbial sequences represented only a small fraction of total reads, 21 genera were detectable across samples. Genus-level profiles, Jaccard-based ordination, and ANCOM-BC analyses consistently revealed clear differences between tanks, whereas no associations were observed for sex or survival status. Three species exhibited significant tank-specific effects, indicating that environmental factors contributed the strongest detectable structure in the data. The limited microbial diversity recovered here reflects the expected constraints of poly(A)-enriched libraries, yet the results demonstrate that unmapped reads from host-derived RNA-seq can still provide informative environmental signatures when analyzed with curated reference databases and compositional statistical approaches. This strategy offers a practical means to extract exploratory microbiome information from existing transcriptomic datasets.

Matching journals

The top 8 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.