SeCAT (Sequence Consensus Amplicon Trimming): a simulation-informed Nextflow pipeline for cross-study harmonisation of 16S rRNA metabarcoding datasets
- 1Nature-based Solutions Research Centre, Aquatic Research Facility, University of Derby
Abstract
Global syntheses of 16S rRNA metabarcoding data rely on combining datasets, yet amplicon sequence variants generated with different primer pairs can rarely be compared directly. Existing approaches correct batch effects statistically or aggregate taxonomically; none first establishes whether the sequences are positionally comparable. We present SeCAT (Sequence Consensus Amplicon Trimming), a Nextflow pipeline that maps every studys variants onto a 16S rRNA full-length reference alignment, determines the largest window shared across studies, and tests by simulation whether trimming to that window degrades community structure before any data are merged. Applied to thirteen marine studies sequenced with three primer sets, twelve converged on a consensus window with a median trimmed read length of 249 bp (9,582 alignment columns). Degradation tests across twelve studies and six taxonomic ranks (72 study-rank combinations) returned 31 passes, 38 cautions and three confirmed failures. Eleven studies (682 samples after filtering) were harmonised (trimmed to the window and merged), collapsing 20,891 primer-specific variants into 12,667 MetaASVs, trimmed sequences deduplicated across studies. Harmonisation preserved community structure at family level (Mantel r = 1.00; Procrustes m{superscript 2} < 0.001) without adding between-study variance (|{Delta}R{superscript 2}| < 0.001). Alpha diversity, including observed richness, was recovered across all levels (intraclass correlations 0.93-1.00), whereas sequence-level congruence was moderate (Procrustes m{superscript 2} = 0.52), as expected where between-study similarity at sequence level exists only after harmonisation. The harmonised output supports synthesis from aggregated ranks down to exact sequences. By establishing positional comparability before merging, SeCAT opens the archive of primer-fragmented 16S rRNA surveys to sequence-level synthesis.
Related
- SedaDNA metabarcoding reveals hidden histories of within-species sequence variation
- Local ETM+ to OLI EVI Harmonisation in Tropical Forest Environments
- Diagnose Harmonisation (DHARM): An open-source tool for the harmonisation and analysis of batch effects in MRI-derived measures.
- TaxHarmonizeR: A data-driven taxonomic harmonisation tool for ecological datasets
- Waterborne eDNA metabarcoding of the upper Llano River watershed (Texas, USA) detects regionally unreported green plant taxa
BibTeX
@misc{terrey2026secat,
title = {{SeCAT (Sequence Consensus Amplicon Trimming): a simulation-informed Nextflow pipeline for cross-study harmonisation of 16S rRNA metabarcoding datasets}},
author = {D. Terrey and P. Huber and E. Pereira and M. C. Pelayo and R. M. W. Ferguson and A. J. Dumbrell and M. Sweet and A. Bani and S. Metz},
howpublished = {bioRxiv},
year = {2026},
month = nov,
doi = {10.64898/2026.11.03.763908},
url = {https://www.biorxiv.org/content/10.64898/2026.11.03.763908v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - SeCAT (Sequence Consensus Amplicon Trimming): a simulation-informed Nextflow pipeline for cross-study harmonisation of 16S rRNA metabarcoding datasets
AU - D. Terrey
AU - P. Huber
AU - E. Pereira
AU - M. C. Pelayo
AU - R. M. W. Ferguson
AU - A. J. Dumbrell
AU - M. Sweet
AU - A. Bani
AU - S. Metz
PY - 2026
DA - 2026/11/03
PB - bioRxiv
DO - 10.64898/2026.11.03.763908
UR - https://www.biorxiv.org/content/10.64898/2026.11.03.763908v1
AB - Global syntheses of 16S rRNA metabarcoding data rely on combining datasets, yet amplicon sequence variants generated with different primer pairs can rarely be compared directly. Existing approaches correct batch effects statistically or aggregate taxonomically; none first establishes whether the sequences are positionally comparable. We present SeCAT (Sequence Consensus Amplicon Trimming), a Nextflow pipeline that maps every studys variants onto a 16S rRNA full-length reference alignment, determines the largest window shared across studies, and tests by simulation whether trimming to that window degrades community structure before any data are merged. Applied to thirteen marine studies sequenced with three primer sets, twelve converged on a consensus window with a median trimmed read length of 249 bp (9,582 alignment columns). Degradation tests across twelve studies and six taxonomic ranks (72 study-rank combinations) returned 31 passes, 38 cautions and three confirmed failures. Eleven studies (682 samples after filtering) were harmonised (trimmed to the window and merged), collapsing 20,891 primer-specific variants into 12,667 MetaASVs, trimmed sequences deduplicated across studies. Harmonisation preserved community structure at family level (Mantel r = 1.00; Procrustes m{superscript 2} < 0.001) without adding between-study variance (|{Delta}R{superscript 2}| < 0.001). Alpha diversity, including observed richness, was recovered across all levels (intraclass correlations 0.93-1.00), whereas sequence-level congruence was moderate (Procrustes m{superscript 2} = 0.52), as expected where between-study similarity at sequence level exists only after harmonisation. The harmonised output supports synthesis from aggregated ranks down to exact sequences. By establishing positional comparability before merging, SeCAT opens the archive of primer-fragmented 16S rRNA surveys to sequence-level synthesis.
ER -