NanoAmpli-Seq: a workflow for amplicon sequencing for mixed microbial communities on the nanopore sequencing platform

A species-level identification pipeline for human gut microbiota based on the V3-V4 regions of 16S rRNA

16S rRNA gene sequencing is pivotal for identifying bacterial species in microbiome studies, especially using the V3-V4 hypervariable regions. A fixed 98.5% similarity threshold is often applied for species-level identification, but this approach can cause misclassification due to varying thresholds among species. To address this, our study integrated data from SILVA, NCBI, and LPSN databases, extracting V3-V4 region sequences and supplementing them with 16S rRNA sequences from 1,082 human gut samples. This resulted in a non-redundant amplicon sequence variants (ASVs) database specific to the V3-V4 regions (positions 341-806). Utilizing this database, we identified flexible classification thresholds for 674 families, 3,661 genera, and 15,735 species, finding clear thresholds for 87.09% of families and 98.38% of genera. For the 896 most common human gut species, we established precise taxonomic thresholds. To leverage these findings, we developed the asvtax pipeline, which applies flexible thresholds for more accurate taxonomic classification, notably improving the identification of new ASVs. The asvtax pipeline not only enhances the precision of species-level classification but also provides a robust framework for analyzing complex microbial communities, facilitating more reliable ecological and functional interpretations in microbiome research.

Breaking free from references: a consensus-based approach for community profiling with long amplicon nanopore data

Third-generation sequencing platforms, such as Oxford Nanopore Technology (ONT), have made it possible to characterize communities through the sequencing of long amplicons. While this theoretically allows for an increased taxonomic resolution compared to short-read sequencing platforms such as Illumina, the high error rate remains problematic for accurately identifying the community members present within a sample. Here, we present and validate CONCOMPRA, a tool that allows the detection of closely related strains within a community by drafting and mapping to consensus sequences. We show that CONCOMPRA outperforms several other tools for profiling bacterial communities using full-length 16S rRNA gene sequencing. Since CONCOMPRA does not rely on a sequence database for profiling communities, it is applicable to systems and amplicons for which little to no reference data exists. Our validation test shows that the amplification of long PCR products is likely to produce chimeric byproducts that inflate alpha diversity and skew community structure, stressing the importance of chimera detection. CONCOMPRA is available on GitHub (https://github.com/willem-stock/CONCOMPRA).

High accuracy meets high throughput for near full-length 16S ribosomal RNA amplicon sequencing on the Nanopore platform

Small subunit (SSU) ribosomal RNA (rRNA) gene amplicon sequencing is a foundational method in microbial ecology. Currently, short-read platforms are commonly employed for high-throughput applications of SSU rRNA amplicon sequencing, but at the cost of poor taxonomic classification due to limited fragment lengths. The Oxford Nanopore Technologies (ONT) platform can sequence full-length SSU rRNA genes, but its lower raw-read accuracy has so-far limited accurate taxonomic classification and de novo feature generation. Here, we present a sequencing workflow, termed ssUMI, that combines unique molecular identifier (UMI)-based error correction with newer (R10.4+) ONT chemistry and sample barcoding to enable high throughput near full-length SSU rRNA (e.g. 16S rRNA) amplicon sequencing. The ssUMI workflow generated near full-length 16S rRNA consensus sequences with 99.99% mean accuracy using a minimum subread coverage of 3×, surpassing the accuracy of Illumina short reads. The consensus sequences generated with ssUMI were used to produce error-free de novo sequence features with no false positives with two microbial community standards. In contrast, Nanopore raw reads produced erroneous de novo sequence features, indicating that UMI-based error correction is currently necessary for high-accuracy microbial profiling with R10.4+ ONT sequencing chemistries. We showcase the cost-competitive scalability of the ssUMI workflow by sequencing 87 time-series wastewater samples and 27 human gut samples, obtaining quantitative ecological insights that were missed by short-read amplicon sequencing. ssUMI, therefore, enables accurate and low-cost full-length 16S rRNA amplicon sequencing on Nanopore, improving accessibility to high-resolution microbiome science.

Primed and ready: nanopore metabarcoding can now recover highly accurate consensus barcodes that are generally indel-free

BACKGROUND

DNA metabarcoding applies high-throughput sequencing approaches to generate numerous DNA barcodes from mixed sample pools for mass species identification and community characterisation. To date, however, most metabarcoding studies employ second-generation sequencing platforms like Illumina, which are limited by short read lengths and longer turnaround times. While third-generation platforms such as the MinION (Oxford Nanopore Technologies) can sequence longer reads and even in real-time, application of these platforms for metabarcoding has remained limited possibly due to the relatively high read error rates as well as the paucity of specialised software for processing such reads.

RESULTS

We show that this is no longer the case by performing nanopore-based, cytochrome c oxidase subunit I (COI) metabarcoding on 34 zooplankton bulk samples, and benchmarking the results against conventional Illumina MiSeq sequencing. Nanopore R10.3 sequencing chemistry and super accurate (SUP) basecalling model reduced raw read error rates to ~ 4%, and consensus calling with amplicon_sorter (without further error correction) generated metabarcodes that were ≤ 1% erroneous. Although Illumina recovered a higher number of molecular operational taxonomic units (MOTUs) than nanopore sequencing (589 vs. 471), we found no significant differences in the zooplankton communities inferred between the sequencing platforms. Importantly, 406 of 444 (91.4%) shared MOTUs between Illumina and nanopore were also found to be free of indel errors, and 85% of the zooplankton richness could be recovered after just 12-15 h of sequencing.

CONCLUSION

Our results demonstrate that nanopore sequencing can generate metabarcodes with Illumina-like accuracy, and we are the first study to show that nanopore metabarcodes are almost always indel-free. We also show that nanopore metabarcoding is viable for characterising species-rich communities rapidly, and that the same ecological conclusions can be obtained regardless of the sequencing platform used. Collectively, our study inspires confidence in nanopore sequencing and paves the way for greater utilisation of nanopore technology in various metabarcoding applications.

R2C2 + UMI: Combining concatemeric and unique molecular identifier-based consensus sequencing enables ultra-accurate sequencing of amplicons on Oxford Nanopore Technologies sequencers

The sequencing of PCR amplicons is a core application of high-throughput sequencing technology. Using unique molecular identifiers (UMIs), individual amplified molecules can be sequenced to very high accuracy on an Illumina sequencer. However, Illumina sequencers have limited read length and are therefore restricted to sequencing amplicons shorter than 600 bp unless using inefficient synthetic long-read approaches. Native long-read sequencers from Pacific Biosciences and Oxford Nanopore Technologies can, using consensus read approaches, match or exceed Illumina quality while achieving much longer read lengths. Using a circularization-based concatemeric consensus sequencing approach (R2C2) paired with UMIs (R2C2 + UMI), we show that we can sequence an ∼550-nt antibody heavy chain (Immunoglobulin heavy chain - IGH) and an ∼1,500-nt 16S amplicons at accuracies up to and exceeding Q50 (<1 error in 100,000 sequenced bases), which exceeds accuracies of UMI-supported Illumina-paired sequencing as well as synthetic long-read approaches.

Environmental DNA metabarcoding of freshwater fish in Malaysian tropical rivers using short-read nanopore sequencing as a potential biomonitoring tool

The approach of combining cost-effective nanopore sequencing and emerging environmental DNA (eDNA) metabarcoding could prove to be a promising tool for biodiversity documentation, especially in Malaysia. Given the substantial funding constraints in recent years, especially in relation to the country's biodiversity, many researchers have been limited to conduct restricted research without extended monitoring periods, potentially hindering comprehensive surveys and could compromise the conservation efforts. Therefore, the present study aimed to evaluate the application of eDNA metabarcoding on freshwater fish using short reads generated through nanopore sequencing. This assessment focused on species detection in three selected rivers within the Endau Rompin Landscape in Malaysia. Additionally, the study compared levels of species detection between eDNA metabarcoding and conventional sampling methods, examined the effectiveness of primer choice, and applied both metabarcoding and shotgun sequencing to the eDNA approach. We successfully identified a total of 22 and 71 species with an identification threshold of >97% and >90%, respectively, through the MinION platform. The eDNA metabarcoding approach detected over 13% more freshwater fish species than when the conventional method was used. Notably, the distinction in freshwater fish detection between eDNA primers for 12S rRNA and cytochrome oxidase I was insignificant. The cost for eDNA metabarcoding proved to be more effective compared to conventional sampling with cost reduction at 33.4%. With favourable cost-effectiveness and increased species detection, eDNA metabarcoding could complement existing methods, enhance holistic diversity documentation for targeted habitats and facilitate effective conservation planning.

Exploiting long read sequencing to detect azole fungicide resistance mutations in Pyrenophora teres using unique molecular identifiers

Resistance to fungicides is a global challenge as target proteins under selection can evolve rapidly, reducing fungicide efficacy. To manage resistance, detection technologies must be fast and flexible enough to cope with a rapidly increasing number of mutations. The most important agricultural fungicides are azoles that target the ergosterol biosynthetic enzyme sterol 14α-demethylase (CYP51). Mutations associated with azole resistance in the Cyp51 promoter and coding sequence can co-occur in the same allele at different positions and codons, increasing the complexity of resistance detection. Resistance mutations arise rapidly and cannot be detected using traditional amplification-based methods if they are not known. To capture the complexity of azole resistance in two net blotch pathogens of barley we used the Oxford Nanopore MinION to sequence the promoter and coding sequence of Cyp51A. This approach detected all currently known mutations from biologically complex samples increasing the simplicity of resistance detection as multiple alleles can be profiled in a single assay. With the mobility and decreasing cost of long read sequencing, we demonstrate this approach is broadly applicable for characterizing resistance within known agrochemical target sites.

Application of Biological Nanopore Sequencing Technology in the Detection of Microorganisms†

Environmental pollution and the spread of pathogenic microorganisms pose a significant threat to the health of humans and the planet. Thus, understanding and detecting microorganisms is crucial for maintaining a healthy living environment. Nanopore sequencing is a single‐molecule detection method developed in the 1990s that has revolutionized various research fields. It offers several advantages over traditional sequencing methods, including low cost, label‐free, time‐saving detection speed, long sequencing reading, real‐time monitoring, convenient carrying, and other significant advantages. In this review, we summarize the technical principles and characteristics of nanopore sequencing and discuss its applications in amplicon sequencing, metagenome sequencing, and whole‐genome sequencing of environmental microorganisms, as well as its in situ application under some special circumstances. We also analyze the advantages and challenges of nanopore sequencing in microbiology research. Overall, nanopore sequencing has the potential to greatly enhance the detection and understanding of microorganisms in environmental research, but further developments are needed to overcome the current challenges.

An integrative protocol for one-step PCR amplicon library construction and accurate demultiplexing of pooled sequencing data

High-throughput sequencing of amplicons has been widely used to precisely and efficiently identify species compositions and analyze community structures, greatly promoting biological studies involving large amounts of complex samples, especially those involving environmental and pathogen-monitoring ones. Commercial library preparation kits for amplicon sequencing, which generally require multiple steps, including adapter ligation and indexing, are expensive and time-consuming, especially for applications at a large scale. To overcome these limitations, a "one-step PCR approach" has been previously proposed for constructions of amplicon libraries using long fusion primers. However, efficient amplifications of target genes and accurate demultiplexing of pooled sequencing data remain to be addressed. To tackle these, we present an integrative protocol for one-step PCR amplicon library construction (OSPALC). High-quality reads have been generated by this approach to reliably identify species compositions of mock bacterial communities and environmental samples. With this protocol, the amplicon library is constructed through one regular PCR with long primers, and the total cost per DNA/cDNA sample decreases to just 7% of the typical cost via the multi-step PCR approach. Empirically tested primers and optimized PCR conditions to construct OSPALC libraries for 16S rDNA V4 regions are demonstrated as a case study. Tools to design primers targeting at any genomic regions are also presented. In principle, OSPALC can be readily applied to construct amplicon libraries of any target genes using DNA or RNA samples, and will facilitate research in numerous fields.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00182-1.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Single-worm long-read sequencing reveals genome diversity in free-living nematodes

Obtaining sufficient genetic material from a limited biological source is currently the primary operational bottleneck in studies investigating biodiversity and genome evolution. In this study, we employed multiple displacement amplification (MDA) and Smartseq2 to amplify nanograms of genomic DNA and mRNA, respectively, from individual Caenorhabditis elegans. Although reduced genome coverage was observed in repetitive regions, we produced assemblies covering 98% of the reference genome using long-read sequences generated with Oxford Nanopore Technologies (ONT). Annotation with the sequenced transcriptome coupled with the available assembly revealed that gene predictions were more accurate, complete and contained far fewer false positives than de novo transcriptome assembly approaches. We sampled and sequenced the genomes and transcriptomes of 13 nematodes from early-branching species in Chromadoria, Dorylaimia and Enoplia. The basal Chromadoria and Enoplia species had larger genome sizes, ranging from 136.6 to 738.8 Mb, compared with those in the other clades. Nine mitogenomes were fully assembled, and displayed a complete lack of synteny to other species. Phylogenomic analyses based on the new annotations revealed strong support for Enoplia as sister to the rest of Nematoda. Our result demonstrates the robustness of MDA in combination with ONT, paving the way for the study of genome diversity in the phylum Nematoda and beyond.

R2C2+UMI: Combining concatemeric consensus sequencing with unique molecular identifiers enables ultra-accurate sequencing of amplicons on Oxford Nanopore Technologies sequencers

The sequencing of PCR amplicons is a core application of high-throughput sequencing technology. Using unique molecular identifiers (UMIs), individual amplified molecules can be sequenced to very high accuracy on an Illumina sequencer. However, Illumina sequencers have limited read length and are therefore restricted to sequencing amplicons shorter than 600bp unless using inefficient synthetic long-read approaches. Native long-read sequencers from Pacific Biosciences and Oxford Nanopore Technologies can, using consensus read approaches, match or exceed Illumina quality while achieving much longer read lengths. Using a circularization-based concatemeric consensus sequencing approach (R2C2) paired with UMIs (R2C2+UMI) we show that we can sequence ~550nt antibody heavy-chain (IGH) and ~1500nt 16S amplicons at accuracies up to and exceeding Q50 (<1 error in 100,0000 sequenced bases), which exceeds accuracies of UMI-supported Illumina paired sequencing as well as synthetic long-read approaches.

Analysis of Lung Microbiome in COVID-19 Patients during Time of Hospitalization

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is the pathogenic agent of the rapidly spreading pneumonia called coronavirus disease 2019 (COVID-19), primarily infects the respiratory and digestive tract. Several studies have indicated the alterations of the bacterial microbiome in the lower respiratory tract during viral infection. However, both bacterial and fungal microbiota in the lung of COVID-19 patients remained to be explored.

METHODS

In this study, we conducted nanopore sequencing analyses of the lower respiratory tract samples from 38 COVID-19 patients and 26 non-COVID-19 pneumonia controls. Both bacterial and fungal microbiome diversities and microbiota abundances in the lung were compared.

RESULTS

Our results revealed significant differences in lung microbiome between COVID-19 patients and non-COVID-19 controls, which were strongly associated with SARS-CoV-2 infection and clinical status. COVID-19 patients exhibited a notably higher abundance of opportunistic pathogens, particularly Acinetobacter baumannii and Candida spp. Furthermore, the potential pathogens enriched in COVID-19 patients were positively correlated with inflammation indicators.

CONCLUSIONS

Our study highlights the differences in lung microbiome diversity and composition between COVID-19 patients and non-COVID-19 patients. This may contribute to predicting co-pathogens and selecting optimal treatments for respiratory infections caused by SARS-CoV-2.

Nanopore Is Preferable over Illumina for 16S Amplicon Sequencing of the Gut Microbiota When Species-Level Taxonomic Classification, Accurate Estimation of Richness, or Focus on Rare Taxa Is Required

Nanopore sequencing is a promising technology used for 16S rRNA gene amplicon sequencing as it can provide full-length 16S reads and has a low up-front cost that allows research groups to set up their own sequencing workflows. To assess whether Nanopore with the improved error rate of the Kit 12 chemistry should be adopted as the preferred sequencing technology instead of Illumina for 16S amplicon sequencing of the gut microbiota, we used a mock community and human faecal samples to compare diversity, richness, and species-level community structure, as well as the replicability of the results. Nanopore had less noise, better accuracy with the mock community, a higher proportion of reads from the faecal samples classified to species, and better replicability. The difference between the Nanopore and Illumina results of the faecal bacterial community structure was significant but small compared to the variation between samples. The results show that Nanopore is a better choice for 16S rRNA gene amplicon sequencing when the focus is on species-level taxonomic resolution, the investigation of rare taxa, or an accurate estimation of richness. Illumina 16S sequencing should be reserved for communities with many unknown species, and for studies that require the resolution of amplicon sequence variants.

Species-Level Characterization of the Microbiome in Breast Tissues with Different Malignancy and Hormone-Receptor Statuses Using Nanopore Sequencing

Unambiguous evidence indicates that microbes are closely linked to various human diseases, including cancer. Most prior work investigating the microbiome of breast tissue describes an association between compositional differences of microbial species in benign and malignant tissues, but few studies have examined the relative abundance of microbial communities within human breast tissue at the species level. In this work, a total of 44 breast tissue samples including benign and malignant tissues with adjacent normal breast tissue pairs were collected, and Oxford Nanopore long-read sequencing was employed to assess breast tissue microbial signatures. Nearly 900 bacterial species were detected from the four dominant phyla: Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes. The bacteria with the highest abundance in all breast tissues was Ralstonia pickettii, and its relative abundance increased with decreasing malignancy. We further examined the breast-tissue microbiome composition with different hormone-receptor statuses, and the relative abundance of the genus Pseudomonas increased most significantly in breast tissues. Our study provides a rationale for exploring microbiomes associated with breast carcinogenesis and cancer development. Further large-cohort investigation of the breast microbiome is necessary to characterize a microbial risk signature and develop potential microbial-based prevention therapies.

High-resolution phylogenetic and population genetic analysis of microbial communities with RoC-ITS

Microbial communities are inter-connected systems of incredible complexity and dynamism that play crucial roles in health, energy, and the environment. To better understand microbial communities and how they respond to change, it is important to know which microbes are present and their relative abundances at the greatest taxonomic resolution possible. Here, we describe a novel protocol (RoC-ITS) that uses the single-molecule Nanopore sequencing platform to assay the composition of microbial communities at the subspecies designation. Using rolling-circle amplification, this methodology produces long-read sequences from a circular construct containing the complete 16S ribosomal gene and the neighboring internally transcribed spacer (ITS). These long reads can be used to generate a high-fidelity circular consensus sequence. Generally, the ribosomal 16S gene provides phylogenetic information down to the species-level, while the much less conserved ITS region contains strain-level information. When linked together, this combination of markers allows for the identification of individual ribosomal units within a specific organism and the assessment of their relative stoichiometry, as well as the ability to monitor subtle shifts in microbial community composition with a single generic assay. We applied RoC-ITS to an artificial microbial community that was also sequenced using the Illumina platform, to assess its accuracy in quantifying the relative abundance and identity of each species.

Optimization of Oxford Nanopore Technology Sequencing Workflow for Detection of Amplicons in Real Time Using ONT-DART Tool

An optimized, well-tested and validated targeted genomic sequencing-based high-throughput assay is currently not available ready for routine biodefense and biosurveillance applications. Earlier, we addressed this gap by developing and establishing baseline comparisons of a multiplex end-point Polymerase Chain Reaction (PCR) assay followed by Oxford Nanopore Technology (ONT) based amplicon sequencing to real time PCR and customized data processing. Here, we expand upon this effort by identifying the optimal ONT library preparation method for integration into a novel software platform ONT-DART (ONT-Detection of Amplicons in Real-Time). ONT-DART is a dockerized, real-time, amplicon-sequence analysis workflow that is used to reproducibly process and filter read data to support actionable amplicon detection calls based on alignment metrics, within sample statistics, and no-template control data. This analysis pipeline was used to compare four ONT library preparation protocols using R9 and Flongle (FL) flow cells. The two 4-Primer methods tested required the shortest preparation times (5.5 and 6.5 h) for 48 libraries but provided lower fidelity data. The Native Barcoding and Ligation methods required longer preparation times of 8 and 12 h, respectively, and resulted in higher overall data quality. On average, data derived from R9 flow cells produced true positive calls for target organisms more than twice as fast as the lower throughput FL flow cells. These results suggest that utilizing the R9 flowcell with an ONT Native Barcoding amplicon library method in combination with ONT-DART platform analytics provides the best sequencing-based alternative to current PCR-based biodetection methods.

Long-Amplicon Single-Molecule Sequencing Reveals Novel, Trait-Associated Variants of VERNALIZATION1 Homoeologs in Hexaploid Wheat

The gene VERNALIZATION1 (VRN1) is a key controller of vernalization requirement in wheat. The genome of hexaploid wheat (Triticum aestivum) harbors three homoeologous VRN1 loci on chromosomes 5A, 5B, and 5D. Structural sequence variants including small and large deletions and insertions and single nucleotide polymorphisms (SNPs) in the three homoeologous VRN1 genes not only play an important role in the control of vernalization requirement, but also have been reported to be associated with other yield related traits of wheat. Here we used single-molecule sequencing of barcoded long-amplicons to assay the full-length sequences (∼13 kbp plus 700 bp from the promoter sequence) of the three homoeologous VRN1 genes in a panel of 192 predominantly European winter wheat cultivars. Long read sequences revealed previously undetected duplications, insertions and single-nucleotide polymorphisms in the three homoeologous VRN1 genes. All the polymorphisms were confirmed by Sanger sequencing. Sequence analysis showed the predominance of the winter alleles vrn-A1, vrn-B1, and vrn-D1 across the investigated cultivars. Associations of SNPs and structural variations within the three VRN1 genes with 20 economically relevant traits including yield, nodal root-angle index and quality related traits were evaluated at the levels of alleles, haplotypes, and copy number variants. Cultivars carrying structural variants within VRN1 genes showed lower grain yield, protein yield and biomass compared to those with intact genes. Cultivars carrying a single vrn-A1 copy and a unique haplotype with a high number of SNPs were found to have elevated grain yield, kernels per spike and kernels per m2 along with lower grain sedimentation values. In addition, we detected a novel SNP polymorphism within the G-quadruplex region of the promoter of vrn-A1 that was associated with deeper roots in winter wheat. Our findings show that multiplex, single-molecule long-amplicon sequencing is a useful tool for detecting variants in target genes within large plant populations, and can be used to simultaneously assay sequence variants among target multiple gene homoeologs in polyploid crops. Numerous novel VRN1 haplotypes and alleles were identified that showed significantly associations to economically important traits. These polymorphisms were converted into PCR or KASP assays for use in marker-assisted breeding.

Seedborne Cercospora beticola Can Initiate Cercospora Leaf Spot from Sugar Beet (Beta vulgaris) Fruit Tissue

Cercospora leaf spot (CLS) is a globally important disease of sugar beet (Beta vulgaris) caused by the fungus Cercospora beticola. Long-distance movement of C. beticola has been indirectly evidenced in recent population genetic studies, suggesting potential dispersal via seed. Commercial sugar beet "seed" consists of the reproductive fruit (true seed surrounded by maternal pericarp tissue) coated in artificial pellet material. In this study, we confirmed the presence of viable C. beticola in sugar beet fruit for 10 of 37 tested seed lots. All isolates harbored the G143A mutation associated with quinone outside inhibitor resistance, and 32 of 38 isolates had reduced demethylation inhibitor sensitivity (EC₅₀ > 1 µg/ml). Planting of commercial sugar beet seed demonstrated the ability of seedborne inoculum to initiate CLS in sugar beet. C. beticola DNA was detected in DNA isolated from xylem sap, suggesting the vascular system is used to systemically colonize the host. We established nuclear ribosomal internal transcribed spacer region amplicon sequencing using the MinION platform to detect fungi in sugar beet fruit. Fungal sequences from 19 different genera were identified from 11 different sugar beet seed lots, but Fusarium, Alternaria, and Cercospora were consistently the three most dominant taxa, comprising an average of 93% relative read abundance over 11 seed lots. We also present evidence that C. beticola resides in the pericarp of sugar beet fruit rather than the true seed. The presence of seedborne inoculum should be considered when implementing integrated disease management strategies for CLS of sugar beet in the future.

Clinical Performance of Nanopore Targeted Sequencing for Diagnosing Infectious Diseases

The gold standard for confirming bacterial infections is culture-positive, which has a long sample-to-result turnaround time and poor sensitivity for unculturable and fastidious pathogens; therefore, it is hard to guide early, targeted antimicrobial therapy and reduce overuse of broad-spectrum antibiotics. Nanopore targeted sequencing (NTS) is reported to be advantageous in detection speed and range over culture in prior published reports. However, investigation of the clinical performance of NTS is deficient at present. Thus, we assessed the feasibility of NTS for the first time with cohort and systematic comparisons with traditional culture assays and PCR followed by Sanger sequencing. This retrospective study was performed on 472 samples, including 6 specimen types from 436 patients, to evaluate the clinical performance of NTS designed for identifying the microbial composition of various infections. Of these samples, 86.7% were found to be NTS positive, which was significantly higher than culture-positive (26.7%). A total of 425 significant human opportunistic bacteria and fungi detected by NTS were selected to go through validation with PCR followed by Sanger sequencing. The average accuracy rate was 85.2% (maximum 100% created by Cryptococcus neoformans, the last one 66.7% provided by both Staphylococcus haemolyticus and Moraxella osloensis, minimum 0% produced by Burkholderia cepacia). The accuracy rate also varied with sample type; the highest accuracy rate was found in pleural and ascites fluid (95.8%) followed by bronchoalveolar lavage fluid (88.7%), urine (86.8%), and wound secretions (85.0%), while the lowest was present in cerebrospinal fluid (58.8%). NTS had a diagnostic sensitivity of 94.5% and specificity of 31.8%. The positive and negative predictive values of NTS were 79.9% and 66.7%, respectively. For diagnosis of infectious diseases, the sensitivity was greatly increased by 56.7% in NTS compared with culture (94.5% vs 37.8%). Therefore, NTS can accurately detect the causative pathogens in infectious samples, particularly in pleural and ascites fluid, bronchoalveolar lavage fluid, urine, and wound secretions, with a short turnaround time of 8-14 h, and might innovatively contribute to personalizing antibiotic treatments for individuals with standardized protocols in clinical practices. IMPORTANCE Nanopore targeted sequencing (NTS) is reported to be advantageous in detection speed and range over culture in prior published reports. Investigation of the clinical performance of NTS is deficient at present. In our study, cohort and systematic comparisons among three assays (culture, NTS, and Sanger sequencing) were analyzed retrospectively for the first time. We found that NTS undoubtedly has incomparable advantages in accurately detecting the causative pathogens in infectious samples, particularly in pleural and ascites fluid, bronchoalveolar lavage fluid, urine, and wound secretions, with a short turnaround time of 8-14 h. For sterile specimens like blood and cerebrospinal fluid (CSF), the NTS outcomes should be validated using other nucleic acid based detection technology. Overall, NTS might innovatively contribute to guiding early, targeted antimicrobial therapy with lower cost and reduce overuse of broad-spectrum antibiotics.

Amplicon_sorter: A tool for reference‐free amplicon sorting based on sequence similarity and for building consensus sequences

Oxford Nanopore Technologies (ONT) is a third‐generation sequencing technology that is gaining popularity in ecological research for its portable and low‐cost sequencing possibilities. Although the technology excels at long‐read sequencing, it can also be applied to sequence amplicons. The downside of ONT is the low quality of the raw reads. Hence, generating a high‐quality consensus sequence is still a challenge. We present Amplicon_sorter, a tool for reference‐free sorting of ONT sequenced amplicons based on their similarity in sequence and length and for building solid consensus sequences.

Nanopore sequencing technology, bioinformatics and applications

Rapid advances in nanopore technologies for sequencing single long DNA and RNA molecules have led to substantial improvements in accuracy, read length and throughput. These breakthroughs have required extensive development of experimental and bioinformatics methods to fully exploit nanopore long reads for investigations of genomes, transcriptomes, epigenomes and epitranscriptomes. Nanopore sequencing is being applied in genome assembly, full-length transcript detection and base modification detection and in more specialized areas, such as rapid clinical diagnoses and outbreak surveillance. Many opportunities remain for improving data quality and analytical approaches through the development of new nanopores, base-calling methods and experimental protocols tailored to particular applications.

Nanopore-Based Surveillance of Zoonotic Bacterial Pathogens in Farm-Dwelling Peridomestic Rodents

The effective control of rodent populations on farms is crucial for food safety, as rodents are reservoirs and vectors for several zoonotic pathogens. Clear links have been identified between rodents and farm-level outbreaks of pathogens throughout Europe and Asia; however, comparatively little research has been devoted to studying the rodent–agricultural interface in the USA. Here, we address this knowledge gap by metabarcoding bacterial communities of rodent pests collected from Minnesota and Wisconsin food animal farms. We leveraged the Oxford Nanopore MinION sequencer to provide a rapid real-time survey of putative zoonotic foodborne pathogens, among others. Rodents were live trapped (n = 90) from three dairy and mixed animal farms. DNA extraction was performed on 63 rodent colons along with 2 shrew colons included as outgroups in the study. Full-length 16S amplicon sequencing was performed. Our farm-level rodent-metabarcoding data indicate the presence of multiple foodborne pathogens, including Salmonella spp., Campylobacter spp., Staphylococcus aureus, and Clostridium spp., along with many mastitis pathogens circulating within five rodent species (Microtus pennsylvanicus, Mus musculus, Peromyscus leucopus, Peromyscus maniculatus, and Rattus norvegicus) and a shrew (Blarina brevicauda). Interestingly, we observed a higher abundance of enteric pathogens (e.g., Salmonella) in shrew feces compared to the rodents analyzed in our study. Knowledge gained from our research efforts will directly inform and improve farm-level biosecurity efforts and public health interventions to reduce future outbreaks of foodborne and zoonotic disease.

Triangulation of microbial fingerprinting in anaerobic digestion reveals consistent fingerprinting profiles

The anaerobic digestion microbiome has been puzzling us since the dawn of molecular methods for mixed microbial community analysis. Monitoring of the anaerobic digestion microbiome can either take place via a non-targeted holistic evaluation of the microbial community through fingerprinting or by targeted monitoring of selected taxa. Here, we compared four different microbial community fingerprinting methods, i.e., amplicon sequencing, metaproteomics, metabolomics and cytomics, in their ability to characterise the full-scale anaerobic digestion microbiome. Cytometric fingerprinting through cytomics reflects a, for anaerobic digestion, novel, single cell-based approach of direct microbial community fingerprinting by flow cytometry. Three different digester types, i.e., sludge digesters, digesters treating agro-industrial waste and dry anaerobic digesters, each reflected different operational parameters. The α-diversity analysis yielded inconsistent results, especially for richness, across the different methods. In contrast, β-diversity analysis resulted in comparable profiles, even when translated into phyla or functions, with clear separation of the three digester types. In-depth analysis of each method's features i.e., operational taxonomic units, metaproteins, metabolites, and cytometric traits, yielded certain similar features, yet, also some clear differences between the different methods, which was related to the complexity of the anaerobic digestion process. In conclusion, cytometric fingerprinting through flow cytometry is a reliable, fast method for holistic monitoring of the anaerobic digestion microbiome, and the complementary identification of key features through other methods could give rise to a direct interpretation of anaerobic digestion process performance.

Evaluation of DNA extraction protocols from liquid-based cytology specimens for studying cervical microbiota

Cervical microbiota (CM) are considered an important factor affecting the progression of cervical intraepithelial neoplasia (CIN) and are implicated in the persistence of human papillomavirus (HPV). Collection of liquid-based cytology (LBC) samples is routine for cervical cancer screening and HPV genotyping and can be used for long-term cytological biobanking. We sought to determine whether it is possible to access microbial DNA from LBC specimens, and compared the performance of four different extraction protocols: (ZymoBIOMICS DNA Miniprep Kit; QIAamp PowerFecal Pro DNA Kit; QIAamp DNA Mini Kit; and IndiSpin Pathogen Kit) and their ability to capture the diversity of CM from LBC specimens. LBC specimens from 20 patients (stored for 716 ± 105 days) with CIN values of 2 or 3 were each aliquoted for each of the four kits. Loss of microbial diversity due to long-term LBC storage could not be assessed due to lack of fresh LBC samples. Comparisons with other types of cervical sampling were not performed. We observed that all DNA extraction kits provided equivalent accessibility to the cervical microbial DNA within stored LBC samples. Approximately 80% microbial genera were shared among all DNA extraction protocols. Potential kit contaminants were observed as well. Variation between individuals was a significantly greater influence on the observed microbial composition than was the method of DNA extraction. We also observed that HPV16 was significantly associated with community types that were not dominated by Lactobacillus iners.

Skin microbiome alterations in seborrheic dermatitis and dandruff: A systematic review

Seborrheic dermatitis (SD) and dandruff (DF) are common chronic inflammatory skin diseases characterized by recurrent greasy scales, sometimes with erythema and itchiness. Although the exact pathophysiology of the disease is still unclear, current theories highlight the role of microbes on the skin surface in the pathogenesis of SD. Here, we conducted a systematic review to investigate the skin microbiome alterations in patients with SD/DF. We searched Medline/PubMed, Embase and Web of Science for research studies published in English between 1 January 2000 and 31 December 2020. A total of 12 studies with 706 SD/DF samples and 379 healthy samples were included in this study. The scalp and face were predominated by the fungi of Ascomycota and Basidiomycota and the bacteria of Actinobacteria and Firmicutes. In general, the included studies demonstrated an increased Malassezia restricta/Malassezia globosa ratio and a reduction in the Cutibaterium/Staphylococcus ratio in the setting of SD/DF. Staphylococcus was associated with epidermal barrier damage, including elevated levels of trans-epidermal water loss and pH, while Cutibacterium had a positive correlation with water content. Malassezia was also found to be related to an increased itching score and disease severity. Further studies focusing on the interactions between various microbes and the host and microbes can help us to better understand the pathogenesis of SD/DF.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome

Anelloviruses are a ubiquitous component of healthy human viromes and remain highly prevalent after being acquired early in life. The full extent of "anellome" diversity and its evolutionary dynamics remain unexplored. We employed in-depth sequencing of blood-transfusion donor(s)-recipient pairs coupled with public genomic resources for a large-scale assembly of anellovirus genomes and used the data to characterize global and personal anellovirus diversity through time. The breadth of the anellome is much greater than previously appreciated, and individuals harbor unique anellomes and transmit lineages that can persist for several months within a diverse milieu of endemic host lineages. Anellovirus sequence diversity is shaped by extensive recombination at all levels of divergence, hindering traditional phylogenetic analyses. Our findings illuminate the transmission dynamics and vast diversity of anelloviruses and set the foundation for future studies to characterize their biology.

Investigation of Fats, Oils, and Grease Co-digestion With Food Waste in Anaerobic Membrane Bioreactors and the Associated Microbial Community Using MinION Sequencing

Co-digestion of fats, oils, and grease (FOG) with food waste (FW) can improve the energy recovery in anaerobic membrane bioreactors (AnMBRs). Here, we investigated the effect of co-digestion of FW and FOG in AnMBRs at fat mass loading of 0.5, 0.75, and 1.0 kg m^–3 day^–1 with a constant organic loading rate of 5.0 gCOD L^–1 day^–1 in both a single-phase (SP) and two-phase (TP) configuration. A separate mono-digestion of FW at an identical organic loading rate was used as the benchmark. During co-digestion, higher daily biogas production, ranging from 4.0 to 12.0%, was observed in the two-phase methane phase (TP-MP) reactor compared to the SP reactor, but the difference was statistically insignificant (p > 0.05) due to the high variability in daily biogas production. However, the co-digestion of FW with FOG at 1.0 kg m^–3 day^–1 fat loading rate significantly (p < 0.05) improved daily biogas production in both the SP (11.0%) and TP (13.0%) reactors compared to the mono-digestion of FW. Microbial community analyses using cDNA-based MinION sequencing of weekly biomass samples from the AnMBRs revealed the prevalence of Lactobacillus (92.2–95.7% relative activity) and Anaerolineaceae (13.3–57.5% relative activity), which are known as fermenters and fatty acid degraders. Syntrophic fatty acid oxidizers were mostly present in the SP and TP-MP reactors, possibly because of the low pH and short solid retention time (SRT) in the acid phase digesters. A greater abundance of the mcrA gene copies (and methanogens) was observed in the SP and MP reactors compared to the acid-phase (AP) reactors. This study demonstrates that FW and FOG can be effectively co-digested in AnMBRs and is expected to inform full-scale decisions on the optimum fat loading rate.

Is Oxford Nanopore sequencing ready for analyzing complex microbiomes?

This minireview will discuss the improvements in Oxford Nanopore (Oxford; sequencing technology that make the MinION a viable platform for microbial ecology studies. Specific issues being addressed are the increase in sequence accuracy from 65 to 96.5% during the last 5 years, the ability to obtain a quantifiable/predictive signal from the MinION with respect to target molecule abundance, simple-to-use GUI-based pathways for data analysis and the modest additional equipment needs for sequencing in the field. Coupling these recent improvements with the low capital costs for equipment and the reasonable per sample cost makes MinION sequencing an attractive option for virtually any laboratory.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

A pilot study of red complex and three genera subgingival microbiome in periodontitis subjects with and without diabetes, evaluated by MinION platform

Background: Subgingival niche is one biofilm habitat containing rich microbiota, which plays an active role in maintaining the health of periodontal tissue and determining host response. As such, a study of changing subgingival biofilms is important for understanding the effect of a systemic condition. In this study, we compared the occurrence of six bacteria cohabiting in the subgingival area of periodontitis subjects, with (DP, n = 8) and without (NDP, n = 4) diabetes. Methods: The six genus and species of targeted bacteria were confirmed by 16S rRNA amplicon sequencing on MinION nanopore platform. Descriptive statistic was used to describe the obtained data. Results: We found that the six genus and species of targeted bacteria were detected but in different quantities in either group's periodontal pocket. Our data showed that Tannerella forsythia was the most abundant species in subgingival biofilms of the DP group of the red complex bacteria. In contrast, Aggregatibacter sp., which belongs to the phylum of proteobacteria, was present at a relatively lower level. In contrast, Fusobacterium sp., which belongs to orange complex bacteria, showed relative similarities in subgingival biofilms of both groups tested, while Veillonella sp., were abundant in the DP groups.  Conclusions: Our data show that the diversity of classic periodontopathogens increased in the subgingival niche of periodontitis subjects with diabetes. It is the first study in Indonesia to apply MinION-based, full-length 16S rRNA amplicon sequencing in periodontitis patients with and without diabetes.

A pilot study of red complex and three genera subgingival microbiome in periodontitis subjects with and without diabetes, evaluated by MinION platform

Background: Subgingival niche is one biofilm habitat containing rich microbiota, which plays an active role in maintaining the health of periodontal tissue and determining host response. As such, a study of changing subgingival biofilms is important for understanding the effect of a systemic condition. In this study, we compared the occurrence of six bacteria cohabiting in the subgingival area of periodontitis subjects, with (DP, n = 8) and without (NDP, n = 4) diabetes.

Methods: The six genus and species of targeted bacteria were confirmed by 16S rRNA amplicon sequencing on MinION nanopore platform. Descriptive statistic was used to describe the obtained data.

Results: We found that the six genus and species of targeted bacteria were detected but in different quantities in either group's periodontal pocket. Our data showed that Tannerella forsythia was the most abundant species in subgingival biofilms of the DP group of the red complex bacteria. In contrast, Aggregatibacter sp., which belongs to the phylum of proteobacteria, was present at a relatively lower level. In contrast, Fusobacterium sp., which belongs to orange complex bacteria, showed relative similarities in subgingival biofilms of both groups tested, while Veillonella sp., were abundant in the DP groups. 

Conclusions: Our data show that the diversity of classic periodontopathogens increased in the subgingival niche of periodontitis subjects with diabetes. It is the first study in Indonesia to apply MinION-based, full-length 16S rRNA amplicon sequencing in periodontitis patients with and without diabetes.

Freshwater monitoring by nanopore sequencing

While traditional microbiological freshwater tests focus on the detection of specific bacterial indicator species, including pathogens, direct tracing of all aquatic DNA through metagenomics poses a profound alternative. Yet, in situ metagenomic water surveys face substantial challenges in cost and logistics. Here, we present a simple, fast, cost-effective and remotely accessible freshwater diagnostics workflow centred around the portable nanopore sequencing technology. Using defined compositions and spatiotemporal microbiota from surface water of an example river in Cambridge (UK), we provide optimised experimental and bioinformatics guidelines, including a benchmark with twelve taxonomic classification tools for nanopore sequences. We find that nanopore metagenomics can depict the hydrological core microbiome and fine temporal gradients in line with complementary physicochemical measurements. In a public health context, these data feature relevant sewage signals and pathogen maps at species level resolution. We anticipate that this framework will gather momentum for new environmental monitoring initiatives using portable devices.

High-accuracy long-read amplicon sequences using unique molecular identifiers with Nanopore or PacBio sequencing

High-throughput amplicon sequencing of large genomic regions remains challenging for short-read technologies. Here, we report a high-throughput amplicon sequencing approach combining unique molecular identifiers (UMIs) with Oxford Nanopore Technologies (ONT) or Pacific Biosciences circular consensus sequencing, yielding high-accuracy single-molecule consensus sequences of large genomic regions. We applied our approach to sequence ribosomal RNA operon amplicons (~4,500 bp) and genomic sequences (>10,000 bp) of reference microbial communities in which we observed a chimera rate <0.02%. To reach a mean UMI consensus error rate <0.01%, a UMI read coverage of 15× (ONT R10.3), 25× (ONT R9.4.1) and 3× (Pacific Biosciences circular consensus sequencing) is needed, which provides a mean error rate of 0.0042%, 0.0041% and 0.0007%, respectively.

Quantitative PCR based detection system for cyanobacterial geosmin/2-methylisoborneol (2-MIB) events in drinking water sources: Current status and challenges

Taste and odor (T&O) are an important issue in drinking water, aquaculture, recreation and a few other associated industries, and cyanobacteria-relevant geosmin and 2-methylisoborneol (2-MIB) are the two most commonly detected T&O compounds worldwide. A rise in the cyanobacterial blooms and associated geosmin/2-MIB episodes due to anthropogenic activities as well as climate change has led to global concerns for drinking water quality. The increasing awareness for the safe drinking, aquaculture or recreational water systems has boost the demand for rapid, robust, on-site early detection and monitoring system for cyanobacterial geosmin/2-MIB events. In past years, research has indicated quantitative PCR (qPCR) as one of the promising tools for detection of geosmin/2-MIB episodes. It offers advantages of detecting the source organism even at very low concentrations, distinction of odor-producing cyanobacterial strains from non-producers and evaluation of odor producing potential of the cyanobacteria at much faster rates compared to conventional techniques.The present review aims at examining the current status of developed qPCR primers and probes in identifying and detecting the cyanobacterial blooms along with geosmin/2-MIB events. Among the more than 100 articles about cyanobacteria associated geosmin/2-MIB in drinking water systems published after 1990, limited reports (approx. 10 each for geosmin and 2-MIB) focused on qPCR detection and its application in the field. Based on the review of literature, a comprehensive open access global cyanobacterial geosmin/2-MIB events database (CyanoGM Explorer) is curated. It acts as a single platform to access updated information related to origin and geographical distribution of geosmin/2-MIB events, cyanobacterial producers, frequency, and techniques associated with the monitoring of the events. Although a total of 132 cyanobacterial strains from 21 genera and 72 cyanobacterial strains from 13 genera have been reported for geosmin and 2-MIB production, respectively, only 58 geosmin and 28 2-MIB synthesis regions have been assembled in the NCBI database. Based on the identity, geosmin sequences were found to be more diverse in the geosmin synthase conserved/primer design region, compared to 2-MIB synthesis region, hindering the design of universal primers/probes. Emerging technologies such as the bioelectronic nose, Surface Enhanced Raman Scattering (SERS), and nanopore sequencing are discussed for future applications in early on-site detection of geosmin/2-MIB and producers. In the end, the paper also highlights various challenges in applying qPCR as a universal system of monitoring and development of response system for geosmin/2-MIB episodes.

Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology

BACKGROUND

Amplicon sequencing of phylogenetic marker genes, e.g., 16S, 18S, or ITS ribosomal RNA sequences, is still the most commonly used method to determine the composition of microbial communities. Microbial ecologists often have expert knowledge on their biological question and data analysis in general, and most research institutes have computational infrastructures to use the bioinformatics command line tools and workflows for amplicon sequencing analysis, but requirements of bioinformatics skills often limit the efficient and up-to-date use of computational resources.

RESULTS

We present dadasnake, a user-friendly, 1-command Snakemake pipeline that wraps the preprocessing of sequencing reads and the delineation of exact sequence variants by using the favorably benchmarked and widely used DADA2 algorithm with a taxonomic classification and the post-processing of the resultant tables, including hand-off in standard formats. The suitability of the provided default configurations is demonstrated using mock community data from bacteria and archaea, as well as fungi.

CONCLUSIONS

By use of Snakemake, dadasnake makes efficient use of high-performance computing infrastructures. Easy user configuration guarantees flexibility of all steps, including the processing of data from multiple sequencing platforms. It is easy to install dadasnake via conda environments. dadasnake is available at https://github.com/a-h-b/dadasnake.

UMI-linked consensus sequencing enables phylogenetic analysis of directed evolution

The success of protein evolution campaigns is strongly dependent on the sequence context in which mutations are introduced, stemming from pervasive non-additive interactions between a protein's amino acids ('intra-gene epistasis'). Our limited understanding of such epistasis hinders the correct prediction of the functional contributions and adaptive potential of mutations. Here we present a straightforward unique molecular identifier (UMI)-linked consensus sequencing workflow (UMIC-seq) that simplifies mapping of evolutionary trajectories based on full-length sequences. Attaching UMIs to gene variants allows accurate consensus generation for closely related genes with nanopore sequencing. We exemplify the utility of this approach by reconstructing the artificial phylogeny emerging in three rounds of directed evolution of an amine dehydrogenase biocatalyst via ultrahigh throughput droplet screening. Uniquely, we are able to identify lineages and their founding variant, as well as non-additive interactions between mutations within a full gene showing sign epistasis. Access to deep and accurate long reads will facilitate prediction of key beneficial mutations and adaptive potential based on in silico analysis of large sequence datasets.

Interpretations of Environmental Microbial Community Studies Are Biased by the Selected 16S rRNA (Gene) Amplicon Sequencing Pipeline

One of the major methods to identify microbial community composition, to unravel microbial population dynamics, and to explore microbial diversity in environmental samples is high-throughput DNA- or RNA-based 16S rRNA (gene) amplicon sequencing in combination with bioinformatics analyses. However, focusing on environmental samples from contrasting habitats, it was not systematically evaluated (i) which analysis methods provide results that reflect reality most accurately, (ii) how the interpretations of microbial community studies are biased by different analysis methods and (iii) if the most optimal analysis workflow can be implemented in an easy-to-use pipeline. Here, we compared the performance of 16S rRNA (gene) amplicon sequencing analysis tools (i.e., Mothur, QIIME1, QIIME2, and MEGAN) using three mock datasets with known microbial community composition that differed in sequencing quality, species number and abundance distribution (i.e., even or uneven), and phylogenetic diversity (i.e., closely related or well-separated amplicon sequences). Our results showed that QIIME2 outcompeted all other investigated tools in sequence recovery (>10 times fewer false positives), taxonomic assignments (>22% better F-score) and diversity estimates (>5% better assessment), suggesting that this approach is able to reflect the in situ microbial community most accurately. Further analysis of 24 environmental datasets obtained from four contrasting terrestrial and freshwater sites revealed dramatic differences in the resulting microbial community composition for all pipelines at genus level. For instance, at the investigated river water sites Sphaerotilus was only reported when using QIIME1 (8% abundance) and Agitococcus with QIIME1 or QIIME2 (2 or 3% abundance, respectively), but both genera remained undetected when analyzed with Mothur or MEGAN. Since these abundant taxa probably have implications for important biogeochemical cycles (e.g., nitrate and sulfate reduction) at these sites, their detection and semi-quantitative enumeration is crucial for valid interpretations. A high-performance computing conformant workflow was constructed to allow FAIR (Findable, Accessible, Interoperable, and Re-usable) 16S rRNA (gene) amplicon sequence analysis starting from raw sequence files, using the most optimal methods identified in our study. Our presented workflow should be considered for future studies, thereby facilitating the analysis of high-throughput 16S rRNA (gene) sequencing data substantially, while maximizing reliability and confidence in microbial community data analysis.

A Sample-to-Report Solution for Taxonomic Identification of Cultured Bacteria in the Clinical Setting Based on Nanopore Sequencing

Amplicon sequencing of the 16S rRNA gene is commonly used for the identification of bacterial isolates in diagnostic laboratories and mostly relies on the Sanger sequencing method. The latter, however, suffers from a number of limitations, with the most significant being the inability to resolve mixed amplicons when closely related species are coamplified from a mixed culture. This often leads to either increased turnaround time or absence of usable sequence data. Short-read next-generation sequencing (NGS) technologies could solve the mixed amplicon issue but would lack both cost efficiency at low throughput and fast turnaround times. Nanopore sequencing developed by Oxford Nanopore Technologies (ONT) could solve those issues by enabling a flexible number of samples per run and an adjustable sequencing time. Here, we report on the development of a standardized laboratory workflow combined with a fully automated analysis pipeline LORCAN (long read consensus analysis), which together provide a sample-to-report solution for amplicon sequencing and taxonomic identification of the resulting consensus sequences. Validation of the approach was conducted on a panel of reference strains and on clinical samples consisting of single or mixed rRNA amplicons associated with various bacterial genera by direct comparison to the corresponding Sanger sequences. Additionally, simulated read and amplicon mixtures were used to assess LORCAN's behavior when dealing with samples with known cross-contamination levels. We demonstrate that by combining ONT amplicon sequencing results with LORCAN, the accuracy of Sanger sequencing can be closely matched (>99.6% sequence identity) and that mixed samples can be resolved at the single-base resolution level. The presented approach has the potential to significantly improve the flexibility, reliability, and availability of amplicon sequencing in diagnostic settings.

Profiling of Oral Bacterial Communities

The profiling of bacterial communities by the sequencing of housekeeping genes such as that encoding the small subunit ribosomal RNA has revealed the extensive diversity of bacterial life on earth. Standard protocols have been developed and are widely used for this application, but individual habitats may require modification of methods. This review discusses the sequencing and analysis methods most appropriate for the study of the bacterial component of the human oral microbiota. If possible, DNA should be extracted from samples soon after collection. If samples have to be stored for practical reasons, precautions to avoid DNA degradation on freezing should be taken. A critical aspect of profiling oral bacterial communities is the choice of region of the 16S rRNA gene for sequencing. The V1-V2 region provides the best discrimination between species of the genus Streptococcus, the most common genus in the mouth and important in health and disease. The MiSeq platform is most commonly used for sequencing, but long-read technologies are now becoming available that should improve the resolution of analyses. There are a variety of well-established data analysis pipelines available, including mothur and QIIME, which identify sequence reads as phylotypes by comparing them to reference data sets or grouping them into operational taxonomic units. DADA2 has improved sequence error correction capabilities and resolves reads to unique variants. Two curated oral 16S rRNA databases are available: HOMD and CORE. Expert interpretation of community profiles is required, both to detect the presence of contaminating DNA, which is commonly present in the reagents used in analysis, and to differentiate oral and nonoral bacteria and determine the significance of findings. Despite advances in shotgun whole-genome metagenomic methods, oral bacterial community profiling via 16S rRNA sequence analysis remains a valuable technique for the characterization of oral bacterial populations.

Bacterial and Eukaryotic Small-Subunit Amplicon Data Do Not Provide a Quantitative Picture of Microbial Communities, but They Are Reliable in the Context of Ecological Interpretations

High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

High-throughput sequencing (HTS) of gene amplicons is a preferred method of assessing microbial community composition, because it rapidly provides information from a large number of samples at high taxonomic resolution and low costs. However, mock community studies show that HTS data poorly reflect the actual relative abundances of individual phylotypes, casting doubt on the reliability of subsequent statistical analysis and data interpretation. We investigated how accurately HTS data reflect the variability of bacterial and eukaryotic community composition and their relationship with environmental factors in natural samples. For this, we compared results of HTS from three independent aquatic time series (n = 883) with those from an established, quantitative microscopic method (catalyzed reporter deposition-fluorescence in situ hybridization [CARD-FISH]). Relative abundances obtained by CARD-FISH and HTS disagreed for most bacterial and eukaryotic phylotypes. Nevertheless, the two methods identified the same environmental drivers to shape bacterial and eukaryotic communities. Our results show that amplicon data do provide reliable information for their ecological interpretations. Yet, when studying specific phylogenetic groups, it is advisable to combine HTS with quantification using microscopy and/or the addition of internal standards.

IMPORTANCE High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

High-throughput sequencing for community analysis: the promise of DNA barcoding to uncover diversity, relatedness, abundances and interactions in spider communities

Large-scale studies on community ecology are highly desirable but often difficult to accomplish due to the considerable investment of time, labor and, money required to characterize richness, abundance, relatedness, and interactions. Nonetheless, such large-scale perspectives are necessary for understanding the composition, dynamics, and resilience of biological communities. Small invertebrates play a central role in ecosystems, occupying critical positions in the food web and performing a broad variety of ecological functions. However, it has been particularly difficult to adequately characterize communities of these animals because of their exceptionally high diversity and abundance. Spiders in particular fulfill key roles as both predator and prey in terrestrial food webs and are hence an important focus of ecological studies. In recent years, large-scale community analyses have benefitted tremendously from advances in DNA barcoding technology. High-throughput sequencing (HTS), particularly DNA metabarcoding, enables community-wide analyses of diversity and interactions at unprecedented scales and at a fraction of the cost that was previously possible. Here, we review the current state of the application of these technologies to the analysis of spider communities. We discuss amplicon-based DNA barcoding and metabarcoding for the analysis of community diversity and molecular gut content analysis for assessing predator-prey relationships. We also highlight applications of the third generation sequencing technology for long read and portable DNA barcoding. We then address the development of theoretical frameworks for community-level studies, and finally highlight critical gaps and future directions for DNA analysis of spider communities.

Nanopore sequencing of microbial communities reveals the potential role of sea lice as a reservoir for fish pathogens

Caligus rogercresseyi is a copepod ectoparasite with a high prevalence in salmon farms in Chile, causing severe welfare and economic concerns to the sector. Information on the parasite's underpinning mechanisms to support its life strategy is recently being investigated. Due to the critical role of microbiota, this study aimed to characterize the microbiota community associated with C. rogercresseyi from different regions with salmon aquaculture in Chile. Using third-generation sequencing with Nanopore technology (MinION) the full 16S rRNA gene from sea lice obtained from 8 areas distributed over the three main aquaculture regions were sequenced. Microbiota of the parasite is mainly comprised of members of phyla Proteobacteria and Bacteroidetes, and a core microbiota community with 147 taxonomical features was identified, and it was present in sea lice from the three regions. This community accounted for 19% of total identified taxa but more than 70% of the total taxonomical abundance, indicating a strong presence in the parasite. Several taxa with bioactive compound secretory capacity were identified, such as members of genus Pseudoalteromonas and Dokdonia, suggesting a possible role of the lice microbiota during the host infestation processes. Furthermore, the microbiota community was differentially associated with the salmon production, where several potential pathogens such as Vibrio, Tenacibaculum, and Aeromonas in Los Lagos, Aysén, and Magallanes region were identified. Notably, the Chilean salmon industry was initially established in the Los Lagos region but it's currently moving to the south, where different oceanographic conditions coexist with lice populations. The results originated by this study will serve as foundation to investigate putative role of sea lice as vectors for fish pathogens and also as reservoirs for antibiotic-resistant genes.

Computational methods for 16S metabarcoding studies using Nanopore sequencing data

Assessment of bacterial diversity through sequencing of 16S ribosomal RNA (16S rRNA) genes has been an approach widely used in environmental microbiology, particularly since the advent of high-throughput sequencing technologies. An additional innovation introduced by these technologies was the need of developing new strategies to manage and investigate the massive amount of sequencing data generated. This situation stimulated the rapid expansion of the field of bioinformatics with the release of new tools to be applied to the downstream analysis and interpretation of sequencing data mainly generated using Illumina technology. In recent years, a third generation of sequencing technologies has been developed and have been applied in parallel and complementarily to the former sequencing strategies. In particular, Oxford Nanopore Technologies (ONT) introduced nanopore sequencing which has become very popular among molecular ecologists. Nanopore technology offers a low price, portability and fast sequencing throughput. This powerful technology has been recently tested for 16S rRNA analyses showing promising results. However, compared with previous technologies, there is a scarcity of bioinformatic tools and protocols designed specifically for the analysis of Nanopore 16S sequences. Due its notable characteristics, researchers have recently started performing assessments regarding the suitability MinION on 16S rRNA sequencing studies, and have obtained remarkable results. Here we present a review of the state-of-the-art of MinION technology applied to microbiome studies, the current possible application and main challenges for its use on 16S rRNA metabarcoding.

Reintroducing mothur: 10 Years Later

More than 10 years ago, we published the paper describing the mothur software package in Applied and Environmental Microbiology Our goal was to create a comprehensive package that allowed users to analyze amplicon sequence data using the most robust methods available. mothur has helped lead the community through the ongoing sequencing revolution and continues to provide this service to the microbial ecology community. Beyond its success and impact on the field, mothur's development exposed a series of observations that are generally translatable across science. Perhaps the observation that stands out the most is that all science is done in the context of prevailing ideas and available technologies. Although it is easy to criticize choices that were made 10 years ago through a modern lens, if we were to wait for all of the possible limitations to be solved before proceeding, science would stall. Even preceding the development of mothur, it was necessary to address the most important problems and work backwards to other problems that limited access to robust sequence analysis tools. At the same time, we strive to expand mothur's capabilities in a data-driven manner to incorporate new ideas and accommodate changes in data and desires of the research community. It has been edifying to see the benefit that a simple set of tools can bring to so many other researchers.

Portable sequencer in the fight against infectious disease

Infectious disease is still a major threat in the world today. Five decades ago, it was considered soon to be eradicated, but the adaptation of pathogens to environmental pressure, such as antimicrobials, encouraged the emergence and reemergence of infectious disease. The fight with infectious disease starts with prevention, diagnosis, and treatment. Diagnosis can be upheld by observing the cause of disease under the microscope or detecting the presence of nucleic acid and proteins of the pathogens. The molecular techniques span from classical polymerase chain reaction (PCR) to sequencing the nucleic acid composition. Here, we are reviewing the works have been undertaken to utilize a portable sequencer, MinION, in various aspects of infectious disease management.

A rapid approach to profiling diverse fungal communities using the MinION™ nanopore sequencer

The Oxford Nanopore Technologies MinION™ sequencer holds the capability to generate long amplicon reads; however, only a small amount of information is available regarding methodological approaches and the ability to identify a broad diversity of fungal taxa. To assess capabilities, three fungal mock communities were sequenced, each of which had varying ratios of 16 taxa. The data were processed through our selected pipeline. The MinION recovered all mock community members, when mixed at equal ratios. When a taxon was represented at a lower ratio, it was not recovered or decreased in relative abundance. Despite high error rates, highly accurate consensus sequences can be derived. This methodological approach identified all mock community taxa, demonstrating the MinION can be used as a practical alternative to profile fungal communities.

Simultaneous detection and comprehensive analysis of HPV and microbiome status of a cervical liquid-based cytology sample using Nanopore MinION sequencing

Human papillomavirus (HPV) is a major pathogen that causes cervical cancer and many other related diseases. HPV infection related cervical microbiome could be an induce factor of cervical cancer. However, it is uncommon to find a single test on the market that can simultaneously provide information on both HPV and the microbiome. Herein, a novel method was developed in this study to simultaneously detect HPV infection and microbiota composition promptly and accurately. It provides a new and simple way to detect vaginal pathogen situation and also provide valuable information for clinical diagnose. This approach combined multiplex PCR, which targeted both HPV16 E6E7 and full-length 16S rRNA, and Nanopore sequencing to generate enough information to understand the vagina condition of patients. One HPV positive liquid-based cytology (LBC) sample was sequenced and analyzed. After comparing with Illumina sequencing, the results from Nanopore showed a similar microbiome composition. An instant sequencing evaluation showed that 15 min sequencing is enough to identify the top 10 most abundant bacteria. Moreover, two HPV integration sites were identified and verified by Sanger sequencing. This approach has many potential applications in pathogen detection and can potentially aid in providing a more rapid clinical diagnosis.