Microbiota profiling with long amplicons using Nanopore sequencing: full-length 16S rRNA gene and whole rrn operon [version 1; referees: 2 approved, 1 approved with reservations]

Abstract

Background: Profiling the microbiome of low-biomass samples is challenging for metagenomics since these samples often contain DNA from other sources, such as the host or the environment. The usual approach is sequencing specific hypervariable regions of the 16S rRNA gene, which fails to assign taxonomy to genus and species level. Here, we aim to assess long-amplicon PCR-based approaches for assigning taxonomy at the genus and species level. We use Nanopore sequencing with two different markers: full-length 16S rRNA (~1,500 bp) and the whole rrn operon (16S rRNA–ITS–23S rRNA; 4,500 bp). Methods: We sequenced a clinical isolate of Staphylococcus pseudintermedius, two mock communities (HM-783D, Bei Resources; D6306, ZymoBIOMICS™) and two pools of low-biomass samples (dog skin from either the chin or dorsal back), using the MinION™ sequencer 1D PCR barcoding kit. Sequences were pre-processed, and data were analyzed using the WIMP workflow on EPI2ME or Minimap2 software with rrn database. Results: The full-length 16S rRNA and the rrn operon were used to retrieve the microbiota composition at the genus and species level from the bacterial isolate, mock communities and complex skin samples. For the Staphylococcus pseudintermedius isolate, when using EPI2ME, the amplicons were assigned to the correct bacterial species in ~98% of the cases with the rrn operon marker, and in ~68% of the cases with the 16S rRNA gene. In both skin microbiota samples, we detected many species with an environmental origin. In chin, we found different Pseudomonas species in high abundance, whereas in dorsal skin there were more taxa with lower abundances. Conclusions: Both full-length 16S rRNA and the rrn operon retrieved the microbiota composition of simple and complex microbial communities, even from the low-biomass samples such as dog skin. For an increased resolution at the species level, using the rrn operon would be the best choice.