Microbiota profiling with long amplicons using Nanopore sequencing: full-length 16S rRNA gene and the 16S-ITS-23S of the rrn operon

GROND: a quality-checked and publicly available database of full-length 16S-ITS-23S rRNA operon sequences

Sequence comparison of 16S rRNA PCR amplicons is an established approach to taxonomically identify bacterial isolates and profile complex microbial communities. One potential application of recent advances in long-read sequencing technologies is to sequence entire rRNA operons and capture significantly more phylogenetic information compared to sequencing of the 16S rRNA (or regions thereof) alone, with the potential to increase the proportion of amplicons that can be reliably classified to lower taxonomic ranks. Here we describe GROND (Genome-derived Ribosomal Operon Database), a publicly available database of quality-checked 16S-ITS-23S rRNA operons, accompanied by multiple taxonomic classifications. GROND will aid researchers in analysis of their data and act as a standardised database to allow comparison of results between studies.

Benchmarking long-read genome sequence alignment tools for human genomics applications

Background

The utility of long-read genome sequencing platforms has been shown in many fields including whole genome assembly, metagenomics, and amplicon sequencing. Less clear is the applicability of long reads to reference-guided human genomics, which is the foundation of genomic medicine. Here, we benchmark available platform-agnostic alignment tools on datasets from nanopore and single-molecule real-time platforms to understand their suitability in producing a genome representation.

Results

For this study, we leveraged publicly-available data from sample NA12878 generated on Oxford Nanopore and sample NA24385 on Pacific Biosciences platforms. We employed state of the art sequence alignment tools including GraphMap2, long-read aligner (LRA), Minimap2, CoNvex Gap-cost alignMents for Long Reads (NGMLR), and Winnowmap2. Minimap2 and Winnowmap2 were computationally lightweight enough for use at scale, while GraphMap2 was not. NGMLR took a long time and required many resources, but produced alignments each time. LRA was fast, but only worked on Pacific Biosciences data. Each tool widely disagreed on which reads to leave unaligned, affecting the end genome coverage and the number of discoverable breakpoints. No alignment tool independently resolved all large structural variants (1,001-100,000 base pairs) present in the Database of Genome Variants (DGV) for sample NA12878 or the truthset for NA24385.

Conclusions

These results suggest a combined approach is needed for LRS alignments for human genomics. Specifically, leveraging alignments from three tools will be more effective in generating a complete picture of genomic variability. It should be best practice to use an analysis pipeline that generates alignments with both Minimap2 and Winnowmap2 as they are lightweight and yield different views of the genome. Depending on the question at hand, the data available, and the time constraints, NGMLR and LRA are good options for a third tool. If computational resources and time are not a factor for a given case or experiment, NGMLR will provide another view, and another chance to resolve a case. LRA, while fast, did not work on the nanopore data for our cluster, but PacBio results were promising in that those computations completed faster than Minimap2. Due to its significant burden on computational resources and slow run time, Graphmap2 is not an ideal tool for exploration of a whole human genome generated on a long-read sequencing platform.

Long-read MinION™ sequencing of 16S and 16S-ITS-23S rRNA genes provides species-level resolution of Lactobacillaceae in mixed communities

The Lactobacillaceae are lactic acid bacteria harnessed to deliver important outcomes across numerous industries, and their unambiguous, species-level identification from mixed community environments is an important endeavor. Amplicon-based metataxonomics using short-read sequencing of partial 16S rRNA gene regions is widely used to support this, however, the high genetic similarity among Lactobacillaceae species restricts our ability to confidently describe these communities even at genus level. Long-read sequencing (LRS) of the whole 16S rRNA gene or the near complete rRNA operon (16S-ITS-23S) has the potential to improve this. We explored species ambiguity amongst Lactobacillaceae using in-silico tool RibDif2, which identified allele overlap when various partial and complete 16S rRNA gene and 16S-ITS-23S rRNA regions were amplified. We subsequently implemented LRS by MinION™ to compare the capacity of V3-V4, 16S and 16S-ITS-23S rRNA amplicons to accurately describe the diversity of a 20-species Lactobacillaceae mock community in practice. In-silico analysis identified more instances of allele/species overlap with V3-V4 amplicons (n = 43) compared to the 16S rRNA gene (n = 11) and partial (n = up to 15) or complete (n = 0) 16S-ITS-23S rRNA amplicons. With subsequent LRS of a DNA mock community, 80% of target species were identified using V3-V4 amplicons whilst the 16S rRNA gene and 16S-ITS-23S rRNA region amplicons resulted in 95 and 100% of target species being identified. A considerable reduction in false-positive identifications was also seen with 16S rRNA gene (n = 3) and 16S-ITS-23S rRNA region (n = 9) amplicons compared with V3-V4 amplicons (n = 43). Whilst the target species affected by allele overlap in V3-V4 and 16S rRNA gene sequenced mock communities were predicted by RibDif2, unpredicted species ambiguity was observed in 16S-ITS-23S rRNA sequenced communities. Considering the average nucleotide identity (ANI) between ambiguous species (~97%) and the basecall accuracy of our MinION™ sequencing protocol (96.4%), the misassignment of reads between closely related taxa is to be expected. With basecall accuracy exceeding 99% for recent MinION™ releases, the increased species-level differentiating power promised by longer amplicons like the 16S-ITS-23S rRNA region, may soon be fully realized.

The newest Oxford Nanopore R10.4.1 full-length 16S rRNA sequencing enables the accurate resolution of species-level microbial community profiling

The long-read amplicon provides a species-level solution for the community. With the improvement of nanopore flowcells, the accuracy of Oxford Nanopore Technologies (ONT) R10.4.1 has been substantially enhanced, with an average of approximately 99%. To evaluate its effectiveness on amplicons, three types of microbiomes were analyzed by 16S ribosomal RNA (hereinafter referred to as "16S") amplicon sequencing using Novaseq, Pacbio sequel II, and Nanopore PromethION platforms (R9.4.1 and R10.4.1) in the current study. We showed the error rate, recall, precision, and bias index in the mock sample. The error rate of ONT R10.4.1 was greatly reduced, with a better recall in the case of the synthetic community. Meanwhile, in different types of environmental samples, ONT R10.4.1 analysis resulted in a composition similar to Pacbio data. We found that classification tools and databases influence ONT data. Based on these results, we conclude that the ONT R10.4.1 16S amplicon can also be used for application in environmental samples. IMPORTANCE The long-read amplicon supplies the community with a species-level solution. Due to the high error rate of nanopore sequencing early on, it has not been frequently used in 16S studies. Oxford Nanopore Technologies (ONT) introduced the R10.4.1 flowcell with Q20+ reagent to achieve more than 99% accuracy as sequencing technology advanced. However, there has been no published study on the performance of commercial PromethION sequencers with R10.4.1 flowcells on 16S sequencing or on the impact of accuracy improvement on taxonomy (R9.4.1 to R10.4.1) using 16S ONT data. In this study, three types of microbiomes were investigated by 16S ribosomal RNA (rRNA) amplicon sequencing using Novaseq, Pacbio sequel II, and Nanopore PromethION platforms (R9.4.1 and R10.4.1). In the mock sample, we displayed the error rate, recall, precision, and bias index. We observed that the error rate in ONT R10.4.1 is significantly lower, especially when deletions are involved. First and foremost, R10.4.1 and Pacific Bioscience platforms reveal a similar microbiome in environmental samples. This study shows that the R10.4.1 full-length 16S rRNA sequences allow for species identification of environmental microbiota.

Integrated Nematode Management in a World in Transition: Constraints, Policy, Processes, and Technologies for the Future

Plant-parasitic nematodes are one of the most insidious pests limiting agricultural production, parasitizing mostly belowground and occasionally aboveground plant parts. They are an important and underestimated component of the estimated 30% yield loss inflicted on crops globally by biotic constraints. Nematode damage is intensified by interactions with biotic and abiotic factors constraints: soilborne pathogens, soil fertility degradation, reduced soil biodiversity, climate variability, and policies influencing the development of improved management options. This review focuses on the following topics: (a) biotic and abiotic constraints, (b) modification of production systems, (c) agricultural policies, (d) the microbiome, (e) genetic solutions, and (f) remote sensing. Improving integrated nematode management (INM) across all scales of agricultural production and along the Global North-Global South divide, where inequalities influence access to technology, is discussed. The importance of the integration of technological development in INM is critical to improving food security and human well-being in the future.

Boxwood phyllosphere fungal and bacterial communities and their differential responses to film-forming anti-desiccants

BACKGROUND

Anti-desiccant is a class of agrochemicals widely used to protect plants from water stresses, rapid temperature variations, heat and sunburn, frost and freeze damages, transplant shock, and pathogen and pest attack. Although anti-desiccants are generally considered non-toxic to organisms, it is unclear whether they may impact the phyllosphere microbial communities. In this study, three film-forming anti-desiccant products, TransFilm, Vapor Gard, and Wilt-Pruf were applied to the canopy of two boxwood cultivars 'Vardar Valley' and 'Justin Brouwers' on April 13 and August 26, 2021. Shoot samples were collected from boxwood plants treated with each of the three products, as well as nontreated control on June 16, August 26 (before the second treatment), and October 18. Microbial and plant genomic DNA was isolated together and 16S rRNA gene and the extended internal transcribed spacer regions were amplified with PCR and sequenced on a Nanopore MinION platform for bacterial and fungal identification.

RESULTS

Bacterial communities were more diverse than fungal communities. At the phylum level, the boxwood phyllosphere was dominated by Proteobacteria and Ascomycota; at the genus level, Methylobacterium and Shiraia were the most abundant bacteria and fungi, respectively. Among the three film-forming anti-desiccants, Vapor Gard and Wilt-Pruf had more impact than TransFilm on the microbial communities. Specifically, broader impacts were observed on fungal than bacterial community composition and structure, with most affected fungi being suppressed while bacteria promoted.

CONCLUSION

This study addressed several major knowledge gaps regarding boxwood phyllosphere microbiota and the impact of anti-desiccants on plant microbiome. We identified diverse microbial communities of boxwood, a major evergreen woody crop and an iconic landscape plant. We also found differential effects of three film-forming anti-desiccants on the composition and structure of bacterial and fungal communities. These findings advanced our understanding of the associated microbiome of this landmark plant, enabling growers to fully utilize the potentials of microbiome and three anti-desiccants in improving boxwood health and productivity.

Benchmarking assembly free nanopore read mappers to classify complex millipede gut microbiota via Oxford Nanopore Sequencing Technology

Millipedes are key players in recycling leaf litter into soil in tropical ecosystems. To elucidate their gut microbiota, we collected millipedes from different municipalities of Puerto Rico. Here we aim to benchmark which method is best for metagenomic skimming of this highly complex millipede microbiome. We sequenced the gut DNA with Oxford Nanopore Technologies' (ONT) MinION sequencer, then analyzed the data using MEGAN-LR, Kraken2 protein mode, Kraken2 nucleotide mode, GraphMap, and Minimap2 to classify these long ONT reads. From our two samples, we obtained a total of 87,110 and 99,749 ONT reads, respectively. Kraken2 nucleotide mode classified the most reads compared to all other methods at the phylum and class taxonomic level, classifying 75% of the reads in the two samples, the other methods failed to assign enough reads to either phylum or class to yield asymptotes in the taxa rarefaction curves indicating that they required more sequencing depth to fully classify this community. The community is hyper diverse with all methods classifying 20-50 phyla in the two samples. There was significant overlap in the reads used and phyla classified between the five methods benchmarked. Our results suggest that Kraken2 nucleotide mode is the most appropriate tool for the application of metagenomic skimming of this highly complex community.

RESCUE: a validated Nanopore pipeline to classify bacteria through long-read, 16S-ITS-23S rRNA sequencing

Despite the advent of third-generation sequencing technologies, modern bacterial ecology studies still use Illumina to sequence small (~400 bp) hypervariable regions of the 16S rRNA SSU for phylogenetic classification. By sequencing a larger region of the rRNA gene operons, the limitations and biases of sequencing small portions can be removed, allowing for more accurate classification with deeper taxonomic resolution. With Nanopore sequencing now providing raw simplex reads with quality scores above Q20 using the kit 12 chemistry, the ease, cost, and portability of Nanopore play a leading role in performing differential bacterial abundance analysis. Sequencing the near-entire rrn operon of bacteria and archaea enables the use of the universally conserved operon holding evolutionary polymorphisms for taxonomic resolution. Here, a reproducible and validated pipeline was developed, RRN-operon Enabled Species-level Classification Using EMU (RESCUE), to facilitate the sequencing of bacterial rrn operons and to support import into phyloseq. Benchmarking RESCUE showed that fully processed reads are now parallel or exceed the quality of Sanger, with median quality scores of approximately Q20+, using the R10.4 and Guppy SUP basecalling. The pipeline was validated through two complex mock samples, the use of multiple sample types, with actual Illumina data, and across four databases. RESCUE sequencing is shown to drastically improve classification to the species level for most taxa and resolves erroneous taxa caused by using short reads such as Illumina.

The presence of Enterococcus faecalis in saliva as a risk factor for endodontic infection

Aim

The aim of the present study was to investigate and correlate the prevalence of Enterococcus faecalis in saliva and in root canals with different pulpal and periapical conditions.

Methodology

Sixty-seven patients were divided into five groups based on pulpal and periapical tissue status: healthy vital teeth (HVT, n=7), healthy treated teeth without lesion (HTT, n=9), irreversible pulpitis (IP, n=13), necrosis (N, n=18), and post-treatment apical periodontitis (PTAP, n=20). Saliva, rubber dam, sterility control and pre-treatment root canal samples were collected and microbiologically processed by culture method. The phylogenetic relationship of E. faecalis isolates collected from root canals and saliva were investigated by whole genome sequencing. Fisher’s exact test was used to correlate the presence of E. faecalis in root canals or saliva with clinical and/or radiographic findings. Linear/logistic regression analyses were performed to establish the relationship between the presence of E. faecalis in root canals, saliva, and the status of periapical tissues.

Results

E. faecalis was found in 18 root canal and saliva samples. E. faecalis root canal isolates were recovered with the highest frequency from post-treatment apical periodontitis. The occurrence of E. faecalis in saliva was strongly associated with its detection in the root canals (P < 0.001). The pretreatment presence of E. faecalis in root canals was associated with significantly higher odds of having periapical lesions (OR=11.03; 95% CI, 1.27-95.70; p < 0.05). Saliva and root canal isolates from the same patient were highly correlated at the phylogenetic level (Jaccard index >0.95).

Conclusion

This pilot study confirms the role of E. faecalis in developing peri-radicular lesions in secondary endodontic infections and suggests that saliva could be the main source of infection. Further studies are needed to investigate the exact origin of this bacteria and its true role in the pathogenesis of secondary/persistent endodontic infections.

Genetic and Physiological Characterization of Soybean-Nodule-Derived Isolates from Bangladeshi Soils Revealed Diverse Array of Bacteria with Potential Bradyrhizobia for Biofertilizers

Genetic and physiological characterization of bacteria derived from nodules of leguminous plants in the exploration of biofertilizer is of paramount importance from agricultural and environmental perspectives. Phylogenetic analysis of the 16S rRNA gene of 84 isolates derived from Bangladeshi soils revealed an unpredictably diverse array of nodule-forming and endosymbiotic bacteria-mostly belonging to the genus Bradyrhizobium. A sequence analysis of the symbiotic genes (nifH and nodD1) revealed similarities with the 16S rRNA gene tree, with few discrepancies. A phylogenetic analysis of the partial rrn operon (16S-ITS-23S) and multi-locus sequence analysis of atpD, glnII, and gyrB identified that the Bradyrhizobium isolates belonged to Bradyrhizobium diazoefficiens, Bradyrhizobium elkanii, Bradyrhizobium liaoningense and Bradyrhizobium yuanmingense species. In the pot experiment, several isolates showed better activity than B. diazoefficiens USDA110, and the Bho-P2-B2-S1-51 isolate of B. liaoningense showed significantly higher acetylene reduction activity in both Glycine max cv. Enrei and Binasoybean-3 varieties and biomass production increased by 9% in the Binasoybean-3 variety. Tha-P2-B1-S1-68 isolate of B. diazoefficiens significantly enhanced shoot length and induced 10% biomass production in Binasoybean-3. These isolates grew at 1-4% NaCl concentration and pH 4.5-10 and survived at 45 °C, making the isolates potential candidates for eco-friendly soybean biofertilizers in salty and tropical regions.

Metagenomics Next Generation Sequencing (mNGS): An Exciting Tool for Early and Accurate Diagnostic of Fungal Pathogens in Plants

Crop output is directly impacted by infections, with fungi as the major plant pathogens, making accurate diagnosis of these threats crucial. Developing technology and multidisciplinary approaches are turning to genomic analyses in addition to traditional culture methods in diagnostics of fungal plant pathogens. The metagenomic next-generation sequencing (mNGS) method is preferred for genotyping identification of organisms, identification at the species level, illumination of metabolic pathways, and determination of microbiota. Moreover, the data obtained so far show that this new approach is promising as an emerging new trend in fungal disease detection. Another approach covered by mNGS technologies, known as metabarcoding, enables use of specific markers specific to a genetic region and allows for genotypic identification by facilitating the sequencing of certain regions. Although the core concept of mNGS remains constant across applications, the specific sequencing methods and bioinformatics tools used to analyze the data differ. In this review, we focus on how mNGS technology, including metabarcoding, is applied for detecting fungal pathogens and its promising developments for the future.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

Deep Metagenomic Sequencing for Endophthalmitis Pathogen Detection Using a Nanopore Platform

PURPOSE

To evaluate the utility of nanopore sequencing for identifying potential causative pathogens in endophthalmitis, comparing culture results against full-length 16S rRNA nanopore sequencing (16S Nanopore), whole genome nanopore sequencing (Nanopore WGS), and Illumina (Illumina WGS).

DESIGN

Cross-sectional diagnostic comparison.

METHODS

Patients with clinically suspected endophthalmitis underwent intraocular vitreous biopsy as per standard care. Clinical samples were cultured by conventional methods, together with full-length 16S rRNA and WGS using nanopore and Illumina sequencing platforms.

RESULTS

Of 23 patients (median age 68.5 years [range 47-88]; 14 males [61%]), 18 cases were culture-positive. Nanopore sequencing identified the same cultured organism in all of the culture-positive cases and identified potential pathogens in two culture-negative cases (40%). Nanopore WGS was able to additionally detect the presence of bacteriophages in three samples. The agreements at genus level between culture and 16S Nanopore, Nanopore WGS, and Illumina WGS were 75%, 100%, and 78%, respectively.

CONCLUSIONS

Whole genome sequencing has higher sensitivity and provides a viable alternative to culture and 16S sequencing for detecting potential pathogens in endophthalmitis. Moreover, WGS has the ability to detect other potential pathogens in culture-negative cases. Whilst Nanopore and Illumina WGS provide comparable data, nanopore sequencing provides potential for cost-effective point-of-care diagnostics.

16S-ITGDB: An Integrated Database for Improving Species Classification of Prokaryotic 16S Ribosomal RNA Sequences

Analyzing 16S ribosomal RNA (rRNA) sequences allows researchers to elucidate the prokaryotic composition of an environment. In recent years, third-generation sequencing technology has provided opportunities for researchers to perform full-length sequence analysis of bacterial 16S rRNA. RDP, SILVA, and Greengenes are the most widely used 16S rRNA databases. Many 16S rRNA classifiers have used these databases as a reference for taxonomic assignment tasks. However, some of the prokaryotic taxonomies only exist in one of the three databases. Furthermore, Greengenes and SILVA include a considerable number of taxonomies that do not have the resolution to the species level, which has limited the classifiers’ performance. In order to improve the accuracy of taxonomic assignment at the species level for full-length 16S rRNA sequences, we manually curated the three databases and removed the sequences that did not have a species name. We then established a taxonomy-based integrated database by considering both taxonomies and sequences from all three 16S rRNA databases and validated it by a mock community. Results showed that our taxonomy-based integrated database had improved taxonomic resolution to the species level. The integrated database and the related datasets are available at https://github.com/yphsieh/ItgDB.

Investigation of memory-enhancing effects of Streptococcus thermophilus EG007 in mice and elucidating molecular and metagenomic characteristics using nanopore sequencing

Over the past decades, accumulating evidences have highlighted the gut microbiota as a key player in the brain functioning via microbiota–gut–brain axis, and accordingly, the beneficial role of several probiotic strains in cognitive ability also have been actively investigated. However, the majority of the research have demonstrated the effects against age-related cognitive decline or neurological disease. To this end, we aimed to investigate lactic acid bacteria strains having beneficial effects on the cognitive function of healthy young mice and elucidate underlying characteristics by carrying out nanopore sequencing-based genomics and metagenomics analysis. 8-week consumption of Streptococcus thermophilus EG007 demonstrated marked enhancements in behavior tests assessing short-term spatial and non-spatial learning and memory. It was revealed that EG007 possessed genes encoding various metabolites beneficial for a health condition in many aspects, including gamma-aminobutyric acid producing system, a neurotransmitter associated with mood and stress response. Also, by utilizing 16S–23S rRNA operon as a taxonomic marker, we identified more accurate species-level compositional changes in gut microbiota, which was increase of certain species, previously reported to have associations with mental health or down-regulation of inflammation or infection-related species. Moreover, correlation analysis revealed that the EG007-mediated altered microbiota had a significant correlation with the memory traits.

Microbial Identification Using rRNA Operon Region: Database and Tool for Metataxonomics with Long-Read Sequence

Recent development of long-read sequencing platforms has enabled researchers to explore bacterial community structure through analysis of full-length 16S rRNA gene (∼1,500 bp) or 16S-ITS-23S rRNA operon region (∼4,300 bp), resulting in higher taxonomic resolution than short-read sequencing platforms. Despite the potential of long-read sequencing in metagenomics, resources and protocols for this technology are scarce. Here, we describe MIrROR, the database and analysis tool for metataxonomics using the bacterial 16S-ITS-23S rRNA operon region. We collected 16S-ITS-23S rRNA operon sequences extracted from bacterial genomes from NCBI GenBank and performed curation. A total of 97,781 16S-ITS-23S rRNA operon sequences covering 9,485 species from 43,653 genomes were obtained. For user convenience, we provide an analysis tool based on a mapping strategy that can be used for taxonomic profiling with MIrROR database. To benchmark MIrROR, we compared performance against publicly available databases and tool with mock communities and simulated data sets. Our platform showed promising results in terms of the number of species covered and the accuracy of classification. To encourage active 16S-ITS-23S rRNA operon analysis in the field, BLAST function and taxonomic profiling results with 16S-ITS-23S rRNA operon studies, which have been reported as BioProject on NCBI are provided. MIrROR (http://mirror.egnome.co.kr/) will be a useful platform for researchers who want to perform high-resolution metagenome analysis with a cost-effective sequencer such as MinION from Oxford Nanopore Technologies. IMPORTANCE Metabarcoding is a powerful tool to investigate community diversity in an economic and efficient way by amplifying a specific gene marker region. With the advancement of long-read sequencing technologies, the field of metabarcoding has entered a new phase. The technologies have brought a need for development in several areas, including new markers that long-read can cover, database for the markers, tools that reflect long-read characteristics, and compatibility with downstream analysis tools. By constructing MIrROR, we met the need for a database and tools for the 16S-ITS-23S rRNA operon region, which has recently been shown to have sufficient resolution at the species level. Bacterial community analysis using the 16S-ITS-23S rRNA operon region with MIrROR will provide new insights from various research fields.

Rock-inhabiting fungi: terminology, diversity, evolution and adaptation mechanisms

Rock-inhabiting fungi (RIF) constitute an ecological group associated with terrestrial rocks. This association is generally restricted to the persistent colonisation of rocks and peculiar morphological features based on melanisation and slow growth, which endow RIF with significance in eukaryotic biology, special status in ecology, and exotic potential in biotechnology. There is a need to achieve a better understanding of the hidden biodiversity, antistress biology, origin and convergent evolution of RIF, which will facilitate cultural relic preservation, exploitation of the biogeochemical cycle of rock elements and biotechnology applications. This review focuses on summarising the current knowledge of rock-inhabiting fungi, with particular reference to terminology, biodiversity and geographic distribution, origin and evolution, and stress adaptation mechanisms. We especially teased out the definition through summing up the terms related to rock-inhabting fungi, and also provided a checklist of rock-inhabiting fungal taxa recorded following updated classification schemes.

A ribosomal operon database and MegaBLAST settings for strain-level resolution of microbiomes

Current methods to characterize microbial communities generally employ sequencing of the 16S rRNA gene (<500 bp) with high accuracy (∼99%) but limited phylogenetic resolution. However, long-read sequencing now allows for the profiling of near-full-length ribosomal operons (16S-ITS-23S rRNA genes) on platforms such as the Oxford Nanopore MinION. Here, we describe an rRNA operon database with >300 ,000 entries, representing >10 ,000 prokaryotic species and ∼ 150, 000 strains. Additionally, BLAST parameters were identified for strain-level resolution using in silico mutated, mock rRNA operon sequences (70-95% identity) from four bacterial phyla and two members of the Euryarchaeota, mimicking MinION reads. MegaBLAST settings were determined that required <3 s per read on a Mac Mini with strain-level resolution for sequences with >84% identity. These settings were tested on rRNA operon libraries from the human respiratory tract, farm/forest soils and marine sponges ( n = 1, 322, 818 reads for all sample sets). Most rRNA operon reads in this data set yielded best BLAST hits (95 ± 8%). However, only 38-82% of library reads were compatible with strain-level resolution, reflecting the dominance of human/biomedical-associated prokaryotic entries in the database. Since the MinION and the Mac Mini are both portable, this study demonstrates the possibility of rapid strain-level microbiome analysis in the field.

A Round Trip to the Desert: In situ Nanopore Sequencing Informs Targeted Bioprospecting

Bioprospecting expeditions are often performed in remote locations, in order to access previously unexplored samples. Nevertheless, the actual potential of those samples is only assessed once scientists are back in the laboratory, where a time-consuming screening must take place. This work evaluates the suitability of using Nanopore sequencing during a journey to the Tabernas Desert (Spain) for forecasting the potential of specific samples in terms of bacterial diversity and prevalence of radiation- and desiccation-resistant taxa, which were the target of the bioprospecting activities. Samples collected during the first day were analyzed through 16S rRNA gene sequencing using a mobile laboratory. Results enabled the identification of locations showing the greatest and the least potential, and a second, informed sampling was performed focusing on those sites. After finishing the expedition, a culture collection of 166 strains belonging to 50 different genera was established. Overall, Nanopore and culturing data correlated well, since samples holding a greater potential at the microbiome level also yielded a more interesting set of microbial isolates, whereas samples showing less biodiversity resulted in a reduced (and redundant) set of culturable bacteria. Thus, we anticipate that portable sequencers hold potential as key, easy-to-use tools for in situ-informed bioprospecting strategies.

First record of Lobothorax typus from Odisha coast, India, with molecular characterization and observation of its biphasic moulting

Two distinct morphologies of Cymothoidae isopod, Lobothorax typus were collected from the marine water of Bay of Bengal, Goapalpur-on-Sea as the first record of this parasite from coastal water of Odisha, India. All specimens were found attached to the buccal region of different individuals of the same host fish Trichiurus lepturus. With the aid of COI gene sequencing and morphological analysis, the individuals were found to be conspecific. The most prominent variation among the two morphologies includes the size of 5th pereonite and pleon length to total body length ratio. These variations are as a result of the biphasic moulting process. Maximum Likelihood tree analysis based on COI gene sequences concluded the monophyletic taxonomy of different buccal attaching genera under the family Cymothoidae which is in congruence with their morphological divergence.

Long-read sequencing to interrogate strain-level variation among adherent-invasive Escherichia coli isolated from human intestinal tissue

Adherent-invasive Escherichia coli (AIEC) is a pathovar linked to inflammatory bowel diseases (IBD), especially Crohn’s disease, and colorectal cancer. AIEC are genetically diverse, and in the absence of a universal molecular signature, are defined by in vitro functional attributes. The relative ability of difference AIEC strains to colonize, persist, and induce inflammation in an IBD-susceptible host is unresolved. To evaluate strain-level variation among tissue-associated E. coli in the intestines, we develop a long-read sequencing approach to identify AIEC by strain that excludes host DNA. We use this approach to distinguish genetically similar strains and assess their fitness in colonizing the intestine. Here we have assembled complete genomes using long-read nanopore sequencing for a model AIEC strain, NC101, and seven strains isolated from the intestinal mucosa of Crohn’s disease and non-Crohn’s tissues. We show these strains can colonize the intestine of IBD susceptible mice and induce inflammatory cytokines from cultured macrophages. We demonstrate that these strains can be quantified and distinguished in the presence of 99.5% mammalian DNA and from within a fecal population. Analysis of global genomic structure and specific sequence variation within the ribosomal RNA operon provides a framework for efficiently tracking strain-level variation of closely-related E. coli and likely other commensal/pathogenic bacteria impacting intestinal inflammation in experimental settings and IBD patients.

A Journey on the Skin Microbiome: Pitfalls and Opportunities

The human skin microbiota is essential for maintaining homeostasis and ensuring barrier functions. Over the years, the characterization of its composition and taxonomic diversity has reached outstanding goals, with more than 10 million bacterial genes collected and cataloged. Nevertheless, the study of the skin microbiota presents specific challenges that need to be addressed in study design. Benchmarking procedures and reproducible and robust analysis workflows for increasing comparability among studies are required. For various reasons and because of specific technical problems, these issues have been investigated in gut microbiota studies, but they have been largely overlooked for skin microbiota. After a short description of the skin microbiota, the review tackles methodological aspects and their pitfalls, covering NGS approaches and high throughput culture-based techniques. Recent insights into the "core" and "transient" types of skin microbiota and how the manipulation of these communities can prevent or combat skin diseases are also covered. Finally, this review includes an overview of the main dermatological diseases, the changes in the microbiota composition associated with them, and the recommended skin sampling procedures. The last section focuses on topical and oral probiotics to improve and maintain skin health, considering their possible applications for skin diseases.

Comprehensive Wet-Bench and Bioinformatics Workflow for Complex Microbiota Using Oxford Nanopore Technologies

The advent of high-throughput sequencing techniques has recently provided an astonishing insight into the composition and function of the human microbiome. Next-generation sequencing (NGS) has become the gold standard for advanced microbiome analysis; however, 3rd generation real-time sequencing, such as Oxford Nanopore Technologies (ONT), enables rapid sequencing from several kilobases to >2 Mb with high resolution. Despite the wide availability and the enormous potential for clinical and translational applications, ONT is poorly standardized in terms of sampling and storage conditions, DNA extraction, library creation, and bioinformatic classification. Here, we present a comprehensive analysis pipeline with sampling, storage, DNA extraction, library preparation, and bioinformatic evaluation for complex microbiomes sequenced with ONT. Our findings from buccal and rectal swabs and DNA extraction experiments indicate that methods that were approved for NGS microbiome analysis cannot be simply adapted to ONT. We recommend using swabs and DNA extractions protocols with extended washing steps. Both 16S rRNA and metagenomic sequencing achieved reliable and reproducible results. Our benchmarking experiments reveal thresholds for analysis parameters that achieved excellent precision, recall, and area under the precision recall values and is superior to existing classifiers (Kraken2, Kaiju, and MetaMaps). Hence, our workflow provides an experimental and bioinformatic pipeline to perform a highly accurate analysis of complex microbial structures from buccal and rectal swabs.

IMPORTANCE Advanced microbiome analysis relies on sequencing of short DNA fragments from microorganisms like bacteria, fungi, and viruses. More recently, long fragment DNA sequencing of 3rd generation sequencing has gained increasing importance and can be rapidly conducted within a few hours due to its potential real-time sequencing. However, the analysis and correct identification of the microbiome relies on a multitude of factors, such as the method of sampling, DNA extraction, sequencing, and bioinformatic analysis. Scientists have used different protocols in the past that do not allow us to compare results across different studies and research fields. Here, we provide a comprehensive workflow from DNA extraction, sequencing, and bioinformatic workflow that allows rapid and accurate analysis of human buccal and rectal swabs with reproducible protocols. This workflow can be readily applied by many scientists from various research fields that aim to use long-fragment microbiome sequencing.

Nanopore Sequencing Provides Rapid and Reliable Insight Into Microbial Profiles of Intensive Care Units

Fast and accurate identification of pathogens is an essential task in healthcare settings. Second-generation sequencing platforms such as Illumina have greatly expanded the capacity with which different organisms can be detected in hospital samples, and third-generation nanopore-driven sequencing devices such as Oxford Nanopore's minION have recently emerged as ideal sequencing platforms for routine healthcare surveillance due to their long-read capacity and high portability. Despite its great potential, protocols and analysis pipelines for nanopore sequencing are still being extensively validated. In this work, we assess the ability of nanopore sequencing to provide reliable community profiles based on 16S rRNA sequencing in comparison to traditional Illumina platforms using samples collected from Intensive Care Units of a hospital in Brazil. While our results demonstrate that lower throughputs may be a shortcoming of the method in more complex samples, we show that the use of single-use Flongle flowcells in nanopore sequencing runs can provide insightful information on the community composition in healthcare settings.

Lactobacillus acidophilus LA14 Alleviates Liver Injury

Although the probiotic Lactobacillus acidophilus LA14 is used worldwide, its effect on liver diseases remains unelucidated. Here, 32 rats were divided into four groups, gavaged with L. acidophilus LA14 (3 × 10⁹ CFU) or phosphate-buffered saline for 7 days, and then intraperitoneally injected with d-galactosamine or saline. After 24 h, blood, liver, ileum, and feces samples were collected for liver injury, inflammation, intestinal barrier, gut microbiota, metabolome, and transcriptome analyses. Pretreatment with L. acidophilus LA14 alleviated the d-galactosamine-induced elevation of serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), and bile acids; mitigated the histological injury to the liver and gut; and suppressed the inflammatory cytokines macrophage inflammatory protein 1α (MIP-1α), MIP-3α, and MCP-1. L. acidophilus LA14 also ameliorated the d-galactosamine-induced dysbiosis of the gut microbiota and metabolism, such as the enrichment of Bacteroides sp. strain dnLKV3 and the depletion of Streptococcus, butanoic acid, and N-acetyl-d-glucosamine. The underlying mechanism of L. acidophilus LA14 included prevention of not only the d-galactosamine-induced upregulation of infection- and tumor-related pathways but also the d-galactosamine-induced downregulation of antioxidation-related pathways during this process, as reflected by the liver transcriptome and proteome analyses. Furthermore, the administration of L. acidophilus LA14 to healthy rats did not alter the tested liver indicators but significantly enriched the beneficial Lactobacillus and Bifidobacterium species, promoted metabolism and regulated pathways to improve immunity. The ability of L. acidophilus LA14 to alleviate liver injury was further confirmed with an acetaminophen-induced mouse model. These results might provide a reference for future studies on the application of L. acidophilus LA14 for the prevention of liver injury.

IMPORTANCE The probiotic Lactobacillus acidophilus LA14 is widely used, but its effect on liver diseases has not been elucidated. We explored the protective effect of L. acidophilus LA14 on the liver using rats with d-galactosamine-induced liver injury. Pretreatment with L. acidophilus LA14 alleviated the d-galactosamine-induced elevation of serum ALT, AST, ALP, and bile acids, mitigated the histological injury to the liver and gut, and suppressed the inflammatory cytokines MIP-1α, MIP-3α, and MCP-1. These effects were correlated with the modulations of the gut microbiome, metabolome, and hepatic gene expression induced by L. acidophilus LA14. Moreover, the ability of L. acidophilus LA14 to alleviate liver injury was further confirmed with an acetaminophen-induced mouse model. These results might provide a reference for future studies on the application of L. acidophilus LA14 for the prevention of liver injury.

Establishment and assessment of an amplicon sequencing method targeting the 16S-ITS-23S rRNA operon for analysis of the equine gut microbiome

Microbial communities are commonly studied by using amplicon sequencing of part of the 16S rRNA gene. Sequencing of the full-length 16S rRNA gene can provide higher taxonomic resolution and accuracy. To obtain even higher taxonomic resolution, with as few false-positives as possible, we assessed a method using long amplicon sequencing targeting the rRNA operon combined with a CCMetagen pipeline. Taxonomic assignment had > 90% accuracy at the species level in a mock sample and at the family level in equine fecal samples, generating similar taxonomic composition as shotgun sequencing. The rRNA operon amplicon sequencing of equine fecal samples underestimated compositional percentages of bacterial strains containing unlinked rRNA genes by a fourth to a third, but unlinked rRNA genes had a limited effect on the overall results. The rRNA operon amplicon sequencing with the A519F + U2428R primer set was able to detect some kind of archaeal genomes such as Methanobacteriales and Methanomicrobiales, whereas full-length 16S rRNA with 27F + 1492R could not. Therefore, we conclude that amplicon sequencing targeting the rRNA operon captures more detailed variations of equine microbiota.

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.

Effect of chlorpyrifos on the earthworm Eudrilus euginae and their gut microbiome by toxicological and metagenomic analysis

The earthworms are important soil invertebrates and play a crucial role in pedogenesis. The application of pesticides and prolonged exposure to pesticides causes mortality of earthworms apart from profoundly affecting the resident gut microbiome. The microbiome plays a significant effect on the metabolic processes associated with earthworms. The pesticide Chlorpyrifos (CPF) was studied for its toxicity on Eudrilus euginae by toxicity studies. The LC₅₀ value of filter paper contact test and acute toxicity test was 3.8 mg/mL and 180 mg/kg. The prolonged exposure of earthworms to pesticide on reproductive toxicity resulted in the mortality of earthworms and absence of cocoon formation. Further, the effects of CPF on the whole gut microbiome of E. euginae was analyzed using a long amplicon Nanopore sequencing. Results indicated no fluctuations with Firmicutes and Bacteroidetes, that were found to be dominant at bacterial phyla level while at the genus level, remarkable differences were noticed. Clostridium dominated the earthworm gut prior to CPF exposure while Bacillus dominated after exposure. Similarly, the fungal members such as Ascomycota and Basidiomycota were observed to dominate the gut of earthworm at the phyla level before and after exposure to CPF. In contrast, Clavispora (65%) was the dominant genus before CPF exposure and Taloromyces (42%) dominated after the CPF exposure. Our study demonstrates the effect of CPF on the mortality of E. euginae while the amplicon sequencing established the unique microbiome of the gut in response to the CPF exposure.

Primer, Pipelines, Parameters: Issues in 16S rRNA Gene Sequencing

Short-amplicon 16S rRNA gene sequencing is currently the method of choice for studies investigating microbiomes. However, comparative studies on differences in procedures are scarce. We sequenced human stool samples and mock communities with increasing complexity using a variety of commonly used protocols. Short amplicons targeting different variable regions (V-regions) or ranges thereof (V1-V2, V1-V3, V3-V4, V4, V4-V5, V6-V8, and V7-V9) were investigated for differences in the composition outcome due to primer choices. Next, the influence of clustering (operational taxonomic units [OTUs], zero-radius OTUs [zOTUs], and amplicon sequence variants [ASVs]), different databases (GreenGenes, the Ribosomal Database Project, Silva, the genomic-based 16S rRNA Database, and The All-Species Living Tree), and bioinformatic settings on taxonomic assignment were also investigated. We present a systematic comparison across all typically used V-regions using well-established primers. While it is known that the primer choice has a significant influence on the resulting microbial composition, we show that microbial profiles generated using different primer pairs need independent validation of performance. Further, comparing data sets across V-regions using different databases might be misleading due to differences in nomenclature (e.g., Enterorhabdus versus Adlercreutzia) and varying precisions in classification down to genus level. Overall, specific but important taxa are not picked up by certain primer pairs (e.g., Bacteroidetes is missed using primers 515F-944R) or due to the database used (e.g., Acetatifactor in GreenGenes and the genomic-based 16S rRNA Database). We found that appropriate truncation of amplicons is essential and different truncated-length combinations should be tested for each study. Finally, specific mock communities of sufficient and adequate complexity are highly recommended.

IMPORTANCE In 16S rRNA gene sequencing, certain bacterial genera were found to be underrepresented or even missing in taxonomic profiles when using unsuitable primer combinations, outdated reference databases, or inadequate pipeline settings. Concerning the last, quality thresholds as well as bioinformatic settings (i.e., clustering approach, analysis pipeline, and specific adjustments such as truncation) are responsible for a number of observed differences between studies. Conclusions drawn by comparing one data set to another (e.g., between publications) appear to be problematic and require independent cross-validation using matching V-regions and uniform data processing. Therefore, we highlight the importance of a thought-out study design including sufficiently complex mock standards and appropriate V-region choice for the sample of interest. The use of processing pipelines and parameters must be tested beforehand.

Evaluation of full-length nanopore 16S sequencing for detection of pathogens in microbial keratitis

BACKGROUND

Microbial keratitis is a leading cause of preventable blindness worldwide. Conventional sampling and culture techniques are time-consuming, with over 40% of cases being culture-negative. Nanopore sequencing technology is portable and capable of generating long sequencing reads in real-time. The aim of this study is to evaluate the potential of nanopore sequencing directly from clinical samples for the diagnosis of bacterial microbial keratitis.

METHODS

Using full-length 16S rRNA amplicon sequences from a defined mock microbial community, we evaluated and benchmarked our bioinformatics analysis pipeline for taxonomic assignment on three different 16S rRNA databases (NCBI 16S RefSeq, RDP and SILVA) with clustering at 97%, 99% and 100% similarities. Next, we optimised the sample collection using an ex vivo porcine model of microbial keratitis to compare DNA recovery rates of 12 different collection methods: 21-gauge needle, PTFE membrane (4 mm and 6 mm), Isohelix™ SK-2S, Sugi® Eyespear, Cotton, Rayon, Dryswab™, Hydraflock®, Albumin-coated, Purflock®, Purfoam and Polyester swabs. As a proof-of-concept study, we then used the sampling technique that provided the highest DNA recovery, along with the optimised bioinformatics pipeline, to prospectively collected samples from patients with suspected microbial keratitis. The resulting nanopore sequencing results were then compared to standard microbiology culture methods.

RESULTS

We found that applying alignment filtering to nanopore sequencing reads and aligning to the NCBI 16S RefSeq database at 100% similarity provided the most accurate bacterial taxa assignment. DNA concentration recovery rates differed significantly between the collection methods (p < 0.001), with the Sugi® Eyespear swab providing the highest mean rank of DNA concentration. Then, applying the optimised collection method and bioinformatics pipeline directly to samples from two patients with suspected microbial keratitis, sequencing results from Patient A were in agreement with culture results, whilst Patient B, with negative culture results and previous antibiotic use, showed agreement between nanopore and Illumina Miseq sequencing results.

CONCLUSION

We have optimised collection methods and demonstrated a novel workflow for identification of bacterial microbial keratitis using full-length 16S nanopore sequencing.

Rapid and real-time identification of fungi up to species level with long amplicon nanopore sequencing from clinical samples

The availability of long-read technologies, like Oxford Nanopore Technologies, provides the opportunity to sequence longer fragments of the fungal ribosomal operon, up to 6 Kb (18S-ITS1-5.8S-ITS2-28S) and to improve the taxonomy assignment of the communities up to species level and in real-time. We assess the applicability for taxonomic assignment of amplicons targeting a 3.5 Kb region (V3 18S-ITS1-5.8S-ITS2-28S D2) and a 6 Kb region (V1 18S-ITS1-5.8S-ITS2-28S D12) with the What's in my pot (WIMP) classifier. We used the ZymoBIOMICSTM mock community and different microbiological fungal cultures as positive controls. Long amplicon sequencing correctly identified Saccharomyces cerevisiae and Cryptococcus neoformans from the mock community and Malassezia pachydermatis, Microsporum canis and Aspergillus fumigatus from the microbiological cultures. Besides, we identified Rhodotorula graminis in a culture mislabelled as Candida spp. We applied the same approach to external otitis in dogs. Malassezia was the dominant fungal genus in dogs' ear skin, whereas Ma. pachydermatis was the main species in the healthy sample. Conversely, we identified a higher representation of Ma. globosa and Ma. sympodialis in otitis affected samples. We demonstrate the suitability of long ribosomal amplicons to characterize the fungal community of complex samples, either healthy or with clinical signs of infection.

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.

Real-Time Culture-Independent Microbial Profiling Onboard the International Space Station Using Nanopore Sequencing

For the past two decades, microbial monitoring of the International Space Station (ISS) has relied on culture-dependent methods that require return to Earth for analysis. This has a number of limitations, with the most significant being bias towards the detection of culturable organisms and the inherent delay between sample collection and ground-based analysis. In recent years, portable and easy-to-use molecular-based tools, such as Oxford Nanopore Technologies’ MinION™ sequencer and miniPCR bio’s miniPCR™ thermal cycler, have been validated onboard the ISS. Here, we report on the development, validation, and implementation of a swab-to-sequencer method that provides a culture-independent solution to real-time microbial profiling onboard the ISS. Method development focused on analysis of swabs collected in a low-biomass environment with limited facility resources and stringent controls on allowed processes and reagents. ISS-optimized procedures included enzymatic DNA extraction from a swab tip, bead-based purifications, altered buffers, and the use of miniPCR and the MinION. Validation was conducted through extensive ground-based assessments comparing current standard culture-dependent and newly developed culture-independent methods. Similar microbial distributions were observed between the two methods; however, as expected, the culture-independent data revealed microbial profiles with greater diversity. Protocol optimization and verification was established during NASA Extreme Environment Mission Operations (NEEMO) analog missions 21 and 22, respectively. Unique microbial profiles obtained from analog testing validated the swab-to-sequencer method in an extreme environment. Finally, four independent swab-to-sequencer experiments were conducted onboard the ISS by two crewmembers. Microorganisms identified from ISS swabs were consistent with historical culture-based data, and primarily consisted of commonly observed human-associated microbes. This simplified method has been streamlined for high ease-of-use for a non-trained crew to complete in an extreme environment, thereby enabling environmental and human health diagnostics in real-time as future missions take us beyond low-Earth orbit.

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.

Microbiota - The Unseen Players in Adult Asthmatic Airways

Modulation of human lung airway physiology by commensal microbiota has become one of the key mechanisms involved in the pathogenesis of adult asthma. Recent evidence suggests that the composition of respiratory microbiota plays a significant role in the manifestation of adult asthma; however, scientific evidence about the relationship between airway microbial diversity and phenotypes of adult asthma is limited. Further research is needed to understand the interactions between the airway microbiota and host immune response to develop microbiota-based strategies in management of adult asthma. This study reviews the advances in culture-independent methods for detection of airway microbiome, the current data about airway microbiota in healthy individuals and in adult patients with asthma with a focus on bacterial communities, and the future research directions in airway microbiome.

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.

High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic Nanopore sequencing

The rise of antimicrobial-resistant Neisseria gonorrhoeae is a significant public health concern. Against this background, rapid culture-independent diagnostics may allow targeted treatment and prevent onward transmission. We have previously shown metagenomic sequencing of urine samples from men with urethral gonorrhea can recover near-complete N. gonorrhoeae genomes. However, disentangling the N. gonorrhoeae genome from metagenomic samples and robustly identifying antimicrobial resistance determinants from error-prone Nanopore sequencing is a substantial bioinformatics challenge. Here, we show an N. gonorrhoeae diagnostic workflow for analysis of metagenomic sequencing data obtained from clinical samples using R9.4.1 Nanopore sequencing. We compared results from simulated and clinical infections with data from known reference strains and Illumina sequencing of isolates cultured from the same patients. We evaluated three Nanopore variant callers and developed a random forest classifier to filter called SNPs. Clair was the most suitable variant caller after SNP filtering. A minimum depth of 20× reads was required to confidently identify resistant determinants over the entire genome. Our findings show that metagenomic Nanopore sequencing can provide reliable diagnostic information in N. gonorrhoeae infection.

A lab in the field: applications of real-time, in situ metagenomic sequencing

High-throughput metagenomic sequencing is considered one of the main technologies fostering the development of microbial ecology. Widely used second-generation sequencers have enabled the analysis of extremely diverse microbial communities, the discovery of novel gene functions, and the comprehension of the metabolic interconnections established among microbial consortia. However, the high cost of the sequencers and the complexity of library preparation and sequencing protocols still hamper the application of metagenomic sequencing in a vast range of real-life applications. In this context, the emergence of portable, third-generation sequencers is becoming a popular alternative for the rapid analysis of microbial communities in particular scenarios, due to their low cost, simplicity of operation, and rapid yield of results. This review discusses the main applications of real-time, in situ metagenomic sequencing developed to date, highlighting the relevance of this technology in current challenges (such as the management of global pathogen outbreaks) and in the next future of industry and clinical diagnosis.

Bacterial Vaginosis: Current Diagnostic Avenues and Future Opportunities

A healthy female genital tract harbors a microbiome dominated by lactic acid and hydrogen peroxide producing bacteria, which provide protection against infections by maintaining a low pH. Changes in the bacterial compositions of the vaginal microbiome can lead to bacterial vaginosis (BV), which is often associated with vaginal inflammation. Bacterial vaginosis increases the risk of acquiring sexually transmitted infections (STIs) like human immunodeficiency virus (HIV) and affects women's reproductive health negatively. In pregnant women, BV can lead to chorioamnionitis and adverse pregnancy outcomes, including preterm premature rupture of the membranes and preterm birth. In order to manage BV effectively, good diagnostic procedures are required. Traditionally clinical and microscopic methods have been used to diagnose BV; however, these methods require skilled staff and time and suffer from reduced sensitivity and specificity. New diagnostics, including highly sensitive and specific point-of-care (POC) tests, treatment modalities and vaccines can be developed based on the identification of biomarkers from the growing pool of vaginal microbiome and vaginal metabolome data. In this review the current and future diagnostic avenues will be discussed.

Untangling Species-Level Composition of Complex Bacterial Communities through a Novel Metagenomic Approach

16S small-subunit (SSU) rRNA gene-based bacterial profiling is the gold standard for cost-effective taxonomic reconstruction of complex bacterial populations down to the genus level. However, it has been proven ineffective in clinical and research settings requiring higher taxonomic resolution. We therefore developed a bacterial profiling method based on the internal transcribed spacer (ITS) region employing optimized primers and a comprehensive ITS database for accurate cataloguing of bacterial communities at (sub)species resolution. Performance of the microbial ITS profiling pipeline was tested through analysis of host-associated, food, and environmental matrices, while its efficacy in clinical settings was assessed through analysis of mucosal biopsy specimens of colorectal cancer, leading to the identification of putative novel biomarkers. The data collected indicate that the proposed pipeline represents a major step forward in cost-effective identification and screening of microbial biomarkers at (sub)species level, with relevant impact in research, industrial, and clinical settings.IMPORTANCE We developed a novel method for accurate cataloguing of bacterial communities at (sub)species level involving amplification of the internal transcribed spacer (ITS) region through optimized primers, followed by next-generation sequencing and taxonomic classification of amplicons by means of a comprehensive database of bacterial ITS sequences. Host-associated, food, and environmental matrices were employed to test the performance of the microbial ITS profiling pipeline. Moreover, mucosal biopsy samples from colorectal cancer patients were analyzed to demonstrate the scientific relevance of this profiling approach in a clinical setting through identification of putative novel biomarkers. The results indicate that the ITS-based profiling pipeline proposed here represents a key metagenomic tool with major relevance for research, industrial, and clinical settings.

Does Intraspecific Variation in rDNA Copy Number Affect Analysis of Microbial Communities?

Amplicon sequencing of partial regions of the ribosomal RNA loci (rDNA) is widely used to profile microbial communities. However, the rDNA is dynamic and can exhibit substantial interspecific and intraspecific variation in copy number in prokaryotes and, especially, in microbial eukaryotes. As change in rDNA copy number is a common response to environmental change, rDNA copy number is not necessarily a property of a species. Variation in rDNA copy number, especially the capacity for large intraspecific changes driven by external cues, complicates analyses of rDNA amplicon sequence data. We highlight the need to (i) interpret amplicon sequence data in light of possible interspecific and intraspecific variation, and (ii) examine the potential plasticity in rDNA copy number as an important ecological factor to better understand how microbial communities are structured in heterogeneous environments.

Host genotype and exercise exhibit species-level selection for members of the gut bacterial communities in the mouse digestive system

The mammalian gut microbiome can potentially impact host health and disease state. It is known that the mouse-genome, eating-behavior, and exercise-status promotes higher taxonomic rank-level alterations (e.g. family to phyla-level) of the gut microbiota. Here, host genotype or activity status was investigated to determine if selection of individual bacterial species or strains could be discerned within the murine digestive system. For this study, the fecal bacterial community of adenylyl cyclase 5 knock-out (AC5KO, n = 7) mice or their wild-type (WT, n = 10) littermates under exercise or sedentary conditions were profiled by sequencing rRNA operons. AC5KO mice were chosen since this genotype displays enhanced longevity/exercise capacity and protects against cardiovascular/metabolic disease. Profiling of rRNA operons using the Oxford MinION yielded 65,706 2-D sequences (after size selection of 3.7-5.7 kb) which were screened against an NCBI 16S rRNA gene database. These sequences were binned into 1,566 different best BLAST hits (BBHs) and counted for each mouse sample. Non-metric multidimensional scaling (NMDS) of the gut microbial community demonstrated clustering by physical activity (p = 0.001) but not by host genotype. Additionally, sequence similarity and phylogenetic analysis demonstrated that different bacterial species (closely related to Muribaculum intestinale and Parasutterella excrementihominis) inhabit AC5KO or WT mice depending on activity status. Other bacterial species of the gut microbiota did not follow such patterning (e.g. Turicibacter sanguinis and Turicimonas muris). Our results support the need of improved taxonomic resolution for better characterization of bacterial communities to deepen our understanding of the role of the gut microbiome on host health.

Species-level evaluation of the human respiratory microbiome

BACKGROUND

Changes to human respiratory tract microbiome may contribute significantly to the progression of respiratory diseases. However, there are few studies examining the relative abundance of microbial communities at the species level along the human respiratory tract.

FINDINGS

Bronchoalveolar lavage, throat swab, mouth rinse, and nasal swab samples were collected from 5 participants. Bacterial ribosomal operons were sequenced using the Oxford Nanopore MinION to determine the relative abundance of bacterial species in 4 compartments along the respiratory tract. More than 1.8 million raw operon reads were obtained from the participants with ∼600,000 rRNA reads passing quality assurance/quality control (70-95% identify; >1,200 bp alignment) by Discontiguous MegaBLAST against the EZ BioCloud 16S rRNA gene database. Nearly 3,600 bacterial species were detected overall (>750 bacterial species within the 5 dominant phyla: Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, and Fusobacteria. The relative abundance of bacterial species along the respiratory tract indicated that most microbes (95%) were being passively transported from outside into the lung. However, a small percentage (<5%) of bacterial species were at higher abundance within the lavage samples. The most abundant lung-enriched bacterial species were Veillonella dispar and Veillonella atypica while the most abundant mouth-associated bacterial species were Streptococcus infantis and Streptococcus mitis.

CONCLUSIONS

Most bacteria detected in lower respiratory samples do not seem to colonize the lung. However, >100 bacterial species were found to be enriched in bronchoalveolar lavage samples (compared to mouth/nose) and may play a substantial role in lung health.

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.

Body lice and bed bug co-infestation in an emergency department patient, Ohio, USA

Body lice and bed bugs are hematophagous insects that parasitize humans. Body lice are established vectors of several bacterial pathogens (e.g. Bartonella quintana, Borrelia recurrentis). Bed bugs are biologically competent vectors of some of the same agents, but their vectorial capacity for these in nature is unclear. In particular, a lack of exposure to louse-borne pathogens in bed bugs in the field could be a factor that limits their contribution to transmission. Here, we describe a case of a patient seen in an urban emergency department who was suffering from infestation with both body lice and bed bugs. Insects were collected from the patient and tested for the presence of louse-borne bacterial pathogens using 16S rRNA gene amplicon sequencing. Although no Bartonella, Borrelia, or Rickettsia were detected, this case provides evidence of ecological overlap between body lice and bed bugs and highlights several potential risk factors for co-infestation. The ecological relationships between bed bugs, body lice, and louse-borne bacteria should be further investigated in the field to determine the frequency of co-infestations and identify possible instances of pathogen infection in bed bugs.

Unlinked rRNA genes are widespread among bacteria and archaea

Ribosomes are essential to cellular life and the genes for their RNA components are the most conserved and transcribed genes in bacteria and archaea. Ribosomal RNA genes are typically organized into a single operon, an arrangement thought to facilitate gene regulation. In reality, some bacteria and archaea do not share this canonical rRNA arrangement-their 16S and 23S rRNA genes are separated across the genome and referred to as "unlinked". This rearrangement has previously been treated as an anomaly or a byproduct of genome degradation in intracellular bacteria. Here, we leverage complete genome and long-read metagenomic data to show that unlinked 16S and 23S rRNA genes are more common than previously thought. Unlinked rRNA genes occur in many phyla, most significantly within Deinococcus-Thermus, Chloroflexi, and Planctomycetes, and occur in differential frequencies across natural environments. We found that up to 41% of rRNA genes in soil were unlinked, in contrast to the human gut, where all sequenced rRNA genes were linked. The frequency of unlinked rRNA genes may reflect meaningful life history traits, as they tend to be associated with a mix of slow-growing free-living species and intracellular species. We speculate that unlinked rRNA genes may confer selective advantages in some environments, though the specific nature of these advantages remains undetermined and worthy of further investigation. More generally, the prevalence of unlinked rRNA genes in poorly-studied taxa serves as a reminder that paradigms derived from model organisms do not necessarily extend to the broader diversity of bacteria and archaea.

Relative Performance of MinION (Oxford Nanopore Technologies) versus Sequel (Pacific Biosciences) Third-Generation Sequencing Instruments in Identification of Agricultural and Forest Fungal Pathogens

Culture-based molecular identification methods have revolutionized detection of pathogens, yet these methods are slow and may yield inconclusive results from environmental materials. The second-generation sequencing tools have much-improved precision and sensitivity of detection, but these analyses are costly and may take several days to months. Of the third-generation sequencing techniques, the portable MinION device (Oxford Nanopore Technologies) has received much attention because of its small size and possibility of rapid analysis at reasonable cost. Here, we compare the relative performances of two third-generation sequencing instruments, MinION and Sequel (Pacific Biosciences), in identification and diagnostics of fungal and oomycete pathogens from conifer (Pinaceae) needles and potato (Solanum tuberosum) leaves and tubers. We demonstrate that the Sequel instrument is efficient for metabarcoding of complex samples, whereas MinION is not suited for this purpose due to a high error rate and multiple biases. However, we find that MinION can be utilized for rapid and accurate identification of dominant pathogenic organisms and other associated organisms from plant tissues following both amplicon-based and PCR-free metagenomics approaches. Using the metagenomics approach with shortened DNA extraction and incubation times, we performed the entire MinION workflow, from sample preparation through DNA extraction, sequencing, bioinformatics, and interpretation, in 2.5 h. We advocate the use of MinION for rapid diagnostics of pathogens and potentially other organisms, but care needs to be taken to control or account for multiple potential technical biases.IMPORTANCE Microbial pathogens cause enormous losses to agriculture and forestry, but current combined culturing- and molecular identification-based detection methods are too slow for rapid identification and application of countermeasures. Here, we develop new and rapid protocols for Oxford Nanopore MinION-based third-generation diagnostics of plant pathogens that greatly improve the speed of diagnostics. However, due to high error rate and technical biases in MinION, the Pacific BioSciences Sequel platform is more useful for in-depth amplicon-based biodiversity monitoring (metabarcoding) from complex environmental samples.

A Comparison of Three Different Bioinformatics Analyses of the 16S-23S rRNA Encoding Region for Bacterial Identification

Rapid and reliable identification of bacterial pathogens directly from patient samples is required for optimizing antimicrobial therapy. Although Sanger sequencing of the 16S ribosomal RNA (rRNA) gene is used as a molecular method, species identification and discrimination is not always achievable for bacteria as their 16S rRNA genes have sometimes high sequence homology. Recently, next generation sequencing (NGS) of the 16S–23S rRNA encoding region has been proposed for reliable identification of pathogens directly from patient samples. However, data analysis is laborious and time-consuming and a database for the complete 16S–23S rRNA encoding region is not available. Therefore, a better, faster, and stronger approach is needed for NGS data analysis of the 16S–23S rRNA encoding region. We compared speed and diagnostic accuracy of different data analysis approaches: de novo assembly followed by Basic Local Alignment Search Tool (BLAST), operational taxonomic unit (OTU) clustering, or mapping using an in-house developed 16S–23S rRNA encoding region database for the identification of bacterial species. De novo assembly followed by BLAST using the in-house database was superior to the other methods, resulting in the shortest turnaround time (2 h and 5 min), approximately 2 h less than OTU clustering and 4.5 h less than mapping, and a sensitivity of 80%. Mapping was the slowest and most laborious data analysis approach with a sensitivity of 60%, whereas OTU clustering was the least laborious approach with 70% sensitivity. Although the in-house database requires more sequence entries to improve the sensitivity, the combination of de novo assembly and BLAST currently appears to be the optimal approach for data analysis.