Ribovore: ribosomal RNA sequence analysis for GenBank submissions and database curation

Alejandro A Schäffer; Richard McVeigh; Barbara Robbertse; Conrad L Schoch; Anjanette Johnston; Beverly A Underwood; Ilene Karsch-Mizrachi; Eric P Nawrocki

doi:10.1186/s12859-021-04316-z

Ribovore: ribosomal RNA sequence analysis for GenBank submissions and database curation

BMC Bioinformatics. 2021 Aug 12;22(1):400. doi: 10.1186/s12859-021-04316-z.

Authors

Alejandro A Schäffer^{1

2}, Richard McVeigh², Barbara Robbertse², Conrad L Schoch², Anjanette Johnston², Beverly A Underwood², Ilene Karsch-Mizrachi², Eric P Nawrocki³

Affiliations

¹ Cancer Data Science Laboratory, National Cancer Institute, National Institutes of Health, Bethesda, MD, 20892, USA.
² National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, 20894, USA.
³ National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, 20894, USA. nawrocke@ncbi.nlm.nih.gov.

Abstract

Background: The DNA sequences encoding ribosomal RNA genes (rRNAs) are commonly used as markers to identify species, including in metagenomics samples that may combine many organismal communities. The 16S small subunit ribosomal RNA (SSU rRNA) gene is typically used to identify bacterial and archaeal species. The nuclear 18S SSU rRNA gene, and 28S large subunit (LSU) rRNA gene have been used as DNA barcodes and for phylogenetic studies in different eukaryote taxonomic groups. Because of their popularity, the National Center for Biotechnology Information (NCBI) receives a disproportionate number of rRNA sequence submissions and BLAST queries. These sequences vary in quality, length, origin (nuclear, mitochondria, plastid), and organism source and can represent any region of the ribosomal cistron.

Results: To improve the timely verification of quality, origin and loci boundaries, we developed Ribovore, a software package for sequence analysis of rRNA sequences. The ribotyper and ribosensor programs are used to validate incoming sequences of bacterial and archaeal SSU rRNA. The ribodbmaker program is used to create high-quality datasets of rRNAs from different taxonomic groups. Key algorithmic steps include comparing candidate sequences against rRNA sequence profile hidden Markov models (HMMs) and covariance models of rRNA sequence and secondary-structure conservation, as well as other tests. Nine freely available blastn rRNA databases created and maintained with Ribovore are used for checking incoming GenBank submissions and used by the blastn browser interface at NCBI. Since 2018, Ribovore has been used to analyze more than 50 million prokaryotic SSU rRNA sequences submitted to GenBank, and to select at least 10,435 fungal rRNA RefSeq records from type material of 8350 taxa.

Conclusion: Ribovore combines single-sequence and profile-based methods to improve GenBank processing and analysis of rRNA sequences. It is a standalone, portable, and extensible software package for the alignment, classification and validation of rRNA sequences. Researchers planning on submitting SSU rRNA sequences to GenBank are encouraged to download and use Ribovore to analyze their sequences prior to submission to determine which sequences are likely to be automatically accepted into GenBank.

Keywords: Alignment; Annotation; Ribosomal RNA; ncRNA.

MeSH terms

DNA, Ribosomal
Databases, Nucleic Acid*
Phylogeny
RNA, Ribosomal*
RNA, Ribosomal, 16S / genetics
RNA, Ribosomal, 18S / genetics
Sequence Analysis, RNA

Substances

DNA, Ribosomal
RNA, Ribosomal
RNA, Ribosomal, 16S
RNA, Ribosomal, 18S