ngsLCA — A toolkit for fast and flexible lowest common ancestor inference and taxonomic profiling of metagenomic data

Research output: Contribution to journalJournal articleResearchpeer-review

Metagenomic data generated from environmental samples is increasingly common in the analysis of modern and ancient biological communities. To obtain taxonomic profiles from this type of data, DNA sequences are aligned against large genomic reference databases and the lowest common ancestor (LCA) needs to be inferred for each sequence with multiple alignments. To date, efforts have mainly focused on improving the speed, sensitivity and specificity of alignment tools, and little effort has been applied to the LCA algorithm that generates the taxonomic profiles from alignments. We present ngsLCA, a command-line toolkit with two separate modules: the main program (in C/C++) performing LCA inference, and an R package for generating tables and visualisations of the taxonomic profiles. ngsLCA processed large datasets in BAM/SAM alignment format 4–11 times faster and used less memory compared to other available programs. It is compatible with the NCBI taxonomy and has flexible parameter settings. Furthermore, the toolkit offers functions for filtering, contamination removal, taxonomic clustering, and multiple ways of visualising the generated taxonomic profiles. ngsLCA bridges a gap in current metagenomic analyses by supplying a computationally light, easy-to-use, accurate, fast and flexible LCA algorithm with R functions for processing and illustrating the taxonomic profiles.

Original languageEnglish
JournalMethods in Ecology and Evolution
Number of pages10
ISSN2041-210X
DOIs
Publication statusE-pub ahead of print - 13 Oct 2022

Bibliographical note

Publisher Copyright:
© 2022 The Authors. Methods in Ecology and Evolution published by John Wiley & Sons Ltd on behalf of British Ecological Society.

    Research areas

  • environmental DNA (eDNA), lowest common ancestor (LCA), metagenomics, next-generation sequencing, sedimentary ancient DNA (sedaDNA), shotgun sequencing, taxonomic profiling, toolkit

ID: 323855167